The Genome in a Bottle Consortium, hosted by the National Institute of Standards and Technology (NIST) is creating reference materials and data for human genome sequencing, as well as methods for genome comparison and benchmarking. Here, we describe a large, diverse set of sequencing data for seven human genomes; five are current or candidate NIST Reference Materials. The pilot genome, NA12878, has been released as NIST RM 8398. We also describe data from two Personal Genome Project trios, one of Ashkenazim Jewish ancestry and one of Chinese ancestry. The data come from 12 technologies: BioNano Genomics, Complete Genomics paired-end and LFR, Ion Proton exome, Oxford Nanopore, Pacific Biosciences, SOLiD, 10X Genomics GemCode WGS, and Illumina exome and WGS paired-end, mate-pair, and synthetic long reads. Cell lines, DNA, and data from these individuals are publicly available. Therefore, we expect these data to be useful for revealing novel information about the human genome and improving sequencing technologies, SNP, indel, and structural variant calling, and de novo assembly.

RNAi screens have implicated hundreds of host proteins as HIV-1 dependency factors (HDFs). While informative, these early studies overlap poorly due to false positives and false negatives. To ameliorate these issues, we combined information from the existing HDF screens together with new screens performed with multiple orthologous RNAi reagents (MORR). In addition to being traditionally validated, the MORR screens and the historical HDF screens were quantitatively integrated by the adaptation of an established analysis program, RIGER, for the collective interpretation of each gene's phenotypic significance. False positives were addressed by the removal of poorly expressed candidates through gene expression filtering, as well as with GESS, which identifies off-target effects. This workflow produced a quantitatively integrated network of genes that modulate HIV-1 replication. We further investigated the roles of GOLGI49, SEC13, and COG in HIV-1 replication. Collectively, the MORR-RIGER method minimized the caveats of RNAi screening and improved our understanding of HIV-1-host cell interactions.

N(6)-methyladenosine (m6A) is the most abundant modified base in eukaryotic mRNA and has been linked to diverse effects on mRNA fate. Current mapping approaches localize m6A residues to transcript regions 100-200 nt long but cannot identify precise m6A positions on a transcriptome-wide level. Here we developed m6A individual-nucleotide-resolution cross-linking and immunoprecipitation (miCLIP) and used it to demonstrate that antibodies to m6A can induce specific mutational signatures at m6A residues after ultraviolet light-induced antibody-RNA cross-linking and reverse transcription. We found that these antibodies similarly induced mutational signatures at N(6),2'-O-dimethyladenosine (m6Am), a modification found at the first nucleotide of certain mRNAs. Using these signatures, we mapped m6A and m6Am at single-nucleotide resolution in human and mouse mRNA and identified small nucleolar RNAs (snoRNAs) as a new class of m6A-containing non-coding RNAs (ncRNAs).

DNA methylation is a chemical modification of cytosine bases that is pivotal for gene regulation, cellular specification and cancer development. Here, we describe an R package, methylKit, that rapidly analyzes genome-wide cytosine epigenetic profiles from high-throughput methylation and hydroxymethylation sequencing experiments. methylKit includes functions for clustering, sample quality visualization, differential methylation analysis and annotation features, thus automating and simplifying many of the steps for discerning statistically significant bases or regions of DNA methylation. Finally, we demonstrate methylKit on breast cancer data, in which we find statistically significant regions of differential methylation and stratify tumor subtypes. methylKit is available at http://code.google.com/p/methylkit.

The IFITMs inhibit influenza A virus (IAV) replication in vitro and in vivo. Here, we establish that the antimycotic heptaen, amphotericin B (AmphoB), prevents IFITM3-mediated restriction of IAV, thereby increasing viral replication. Consistent with its neutralization of IFITM3, a clinical preparation of AmphoB, AmBisome, reduces the majority of interferon's protective effect against IAV in vitro. Mechanistic studies reveal that IFITM1 decreases host-membrane fluidity, suggesting both a possible mechanism for IFITM-mediated restriction and its negation by AmphoB. Notably, we reveal that mice treated with AmBisome succumbed to a normally mild IAV infection, similar to animals deficient in Ifitm3. Therefore, patients receiving antifungal therapy with clinical preparations of AmphoB may be functionally immunocompromised and thus more vulnerable to influenza, as well as other IFITM3-restricted viral infections.

The RNA modification N6-methyladenosine (m6A) post-transcriptionally regulates RNA function. The cellular machinery that controls m6A includes methyltransferases and demethylases that add or remove this modification, as well as m6A-binding YTHDF proteins that promote the translation or degradation of m6A-modified mRNA. We demonstrate that m6A modulates infection by hepatitis C virus (HCV). Depletion of m6A methyltransferases or an m6A demethylase, respectively, increases or decreases infectious HCV particle production. During HCV infection, YTHDF proteins relocalize to lipid droplets, sites of viral assembly, and their depletion increases infectious viral particles. We further mapped m6A sites across the HCV genome and determined that inactivating m6A in one viral genomic region increases viral titer without affecting RNA replication. Additional mapping of m6A on the RNA genomes of other Flaviviridae, including dengue, Zika, yellow fever, and West Nile virus, identifies conserved regions modified by m6A. Altogether, this work identifies m6A as a conserved regulatory mark across Flaviviridae genomes.

JAK1 is a critical effector of pro-inflammatory cytokine signaling and plays important roles in immune function, while abnormal JAK1 activity has been linked to immunological and neoplastic diseases. Specific functions of JAK1 in the context of hematopoiesis, and specifically within hematopoietic stem cells (HSCs), have not clearly been delineated. Here, we show that conditional Jak1 loss in HSCs reduces their self-renewal and markedly alters lymphoid/myeloid differentiation in vivo. Jak1-deficient HSCs exhibit decreased competitiveness in vivo and are unable to rescue hematopoiesis in the setting of myelosuppression. They exhibit increased quiescence, an inability to enter the cell cycle in response to hematopoietic stress, and a marked reduction in cytokine sensing, including in response to type I interferons and IL-3. Moreover, Jak1 loss is not fully rescued by expression of a constitutively active Jak2 allele. Together, these data highlight an essential role for Jak1 in HSC homeostasis and stress responses.

Methylation of the N(6) position of adenosine (m(6)A) is a posttranscriptional modification of RNA with poorly understood prevalence and physiological relevance. The recent discovery that FTO, an obesity risk gene, encodes an m(6)A demethylase implicates m(6)A as an important regulator of physiological processes. Here, we present a method for transcriptome-wide m(6)A localization, which combines m(6)A-specific methylated RNA immunoprecipitation with next-generation sequencing (MeRIP-Seq). We use this method to identify mRNAs of 7,676 mammalian genes that contain m(6)A, indicating that m(6)A is a common base modification of mRNA. The m(6)A modification exhibits tissue-specific regulation and is markedly increased throughout brain development. We find that m(6)A sites are enriched near stop codons and in 3' UTRs, and we uncover an association between m(6)A residues and microRNA-binding sites within 3' UTRs. These findings provide a resource for identifying transcripts that are substrates for adenosine methylation and reveal insights into the epigenetic regulation of the mammalian transcriptome.

Mechanisms regulating mammalian meiotic progression are poorly understood. Here we identify mouse YTHDC2 as a critical component. A screen yielded a sterile mutant, 'ketu', caused by a Ythdc2 missense mutation. Mutant germ cells enter meiosis but proceed prematurely to aberrant metaphase and apoptosis, and display defects in transitioning from spermatogonial to meiotic gene expression programs. ketu phenocopies mutants lacking MEIOC, a YTHDC2 partner. Consistent with roles in post-transcriptional regulation, YTHDC2 is cytoplasmic, has 3'→5' RNA helicase activity in vitro, and has similarity within its YTH domain to an N6-methyladenosine recognition pocket. Orthologs are present throughout metazoans, but are diverged in nematodes and, more dramatically, Drosophilidae, where Bgcn is descended from a Ythdc2 gene duplication. We also uncover similarity between MEIOC and Bam, a Bgcn partner unique to schizophoran flies. We propose that regulation of gene expression by YTHDC2-MEIOC is an evolutionarily ancient strategy for controlling the germline transition into meiosis.

In contrast to infections with human immunodeficiency virus (HIV) in humans and simian immunodeficiency virus (SIV) in macaques, SIV infection of a natural host, sooty mangabeys (Cercocebus atys), is non-pathogenic despite high viraemia. Here we sequenced and assembled the genome of a captive sooty mangabey. We conducted genome-wide comparative analyses of transcript assemblies from C. atys and AIDS-susceptible species, such as humans and macaques, to identify candidates for host genetic factors that influence susceptibility. We identified several immune-related genes in the genome of C. atys that show substantial sequence divergence from macaques or humans. One of these sequence divergences, a C-terminal frameshift in the toll-like receptor-4 (TLR4) gene of C. atys, is associated with a blunted in vitro response to TLR-4 ligands. In addition, we found a major structural change in exons 3-4 of the immune-regulatory protein intercellular adhesion molecule 2 (ICAM-2); expression of this variant leads to reduced cell surface expression of ICAM-2. These data provide a resource for comparative genomic studies of HIV and/or SIV pathogenesis and may help to elucidate the mechanisms by which SIV-infected sooty mangabeys avoid AIDS.

The High Andean Paramo ecosystem is a unique neotropical mountain biome considered a diversity and evolutionary hotspot. Lichens, which are complex symbiotic structures that contain diverse commensal microbial communities, are prevalent in Paramos. There they play vital roles in soil formation and mineral fixation. In this study we analyzed the microbiomes of seven lichen genera in Colombian Paramos using 16S rRNA gene amplicon sequencing and provide the first description of the bacterial communities associated with Cora and Hypotrachyna lichens. Paramo lichen microbiomes varied in diversity indexes and number of OTUs, but were composed predominantly by the phyla Acidobacteria, Actinobacteria, Bacteroidetes, Cyanobacteria, Proteobacteria, and Verrucomicrobia. In the case of Cora and Cladonia, the microbiomes were distinguished based on the identity of the lichen host. While the majority of the lichen-associated microorganisms were not present in all lichens sampled, sixteen taxa shared among this diverse group of lichens suggest a core lichen microbiome that broadens our concept of these symbiotic structures. Additionally, we identified strains producing compounds active against clinically relevant microbial strains. These results indicate that lichen microbiomes from the Paramo ecosystem are diverse and host-specific but share a taxonomic core and can be a source of new bacterial taxa and antimicrobials.

Tumors are composed of many different cell types including cancer cells, fibroblasts, and immune cells. Dissecting functional metabolic differences between cell types within a mixed population can be challenging due to the rapid turnover of metabolites relative to the time needed to isolate cells. To overcome this challenge, we traced isotope-labeled nutrients into macromolecules that turn over more slowly than metabolites. This approach was used to assess differences between cancer cell and fibroblast metabolism in murine pancreatic cancer organoid-fibroblast co-cultures and tumors. Pancreatic cancer cells exhibited increased pyruvate carboxylation relative to fibroblasts, and this flux depended on both pyruvate carboxylase and malic enzyme 1 activity. Consequently, expression of both enzymes in cancer cells was necessary for organoid and tumor growth, demonstrating that dissecting the metabolism of specific cell populations within heterogeneous systems can identify dependencies that may not be evident from studying isolated cells in culture or bulk tissue.

The discovery of the 2,000-year-old Dead Sea Scrolls had an incomparable impact on the historical understanding of Judaism and Christianity. "Piecing together" scroll fragments is like solving jigsaw puzzles with an unknown number of missing parts. We used the fact that most scrolls are made from animal skins to "fingerprint" pieces based on DNA sequences. Genetic sorting of the scrolls illuminates their textual relationship and historical significance. Disambiguating the contested relationship between Jeremiah fragments supplies evidence that some scrolls were brought to the Qumran caves from elsewhere; significantly, they demonstrate that divergent versions of Jeremiah circulated in parallel throughout Israel (ancient Judea). Similarly, patterns discovered in non-biblical scrolls, particularly the Songs of the Sabbath Sacrifice, suggest that the Qumran scrolls represent the broader cultural milieu of the period. Finally, genetic analysis divorces debated fragments from the Qumran scrolls. Our study demonstrates that interdisciplinary approaches enrich the scholar's toolkit.

In addition to mutations, epigenetic alterations are important contributors to malignant transformation and tumor progression. The aim of this work was to identify epigenetic events in which promoter or gene body DNA methylation induces gene expression changes that drive melanocyte malignant transformation and metastasis. We previously developed a linear mouse model of melanoma progression consisting of spontaneously immortalized melanocytes, premalignant melanocytes, a nonmetastatic tumorigenic, and a metastatic cell line. Here, through the integrative analysis of methylome and transcriptome data, we identified the relationship between promoter and/or gene body DNA methylation alterations and gene expression in early, intermediate, and late stages of melanoma progression. We identified adenylate cyclase type 3 (Adcy3) and inositol polyphosphate 4-phosphatase type II (Inpp4b), which affect tumor growth and metastatic potential, respectively. Importantly, the gene expression and DNA methylation profiles found in this murine model of melanoma progression were correlated with available clinical data from large population-based primary melanoma cohorts, revealing potential prognostic markers.

With the rapid advancement of sequencing technologies, next generation sequencing (NGS) analysis has been widely applied in cancer genomics research. More recently, NGS has been adopted in clinical oncology to advance personalized medicine. Clinical applications of precision oncology require accurate tests that can distinguish tumor-specific mutations from artifacts introduced during NGS processes or data analysis. Therefore, there is an urgent need to develop best practices in cancer mutation detection using NGS and the need for standard reference data sets for systematically measuring accuracy and reproducibility across platforms and methods. Within the SEQC2 consortium context, we established paired tumor-normal reference samples and generated whole-genome (WGS) and whole-exome sequencing (WES) data using sixteen library protocols, seven sequencing platforms at six different centers. We systematically interrogated somatic mutations in the reference samples to identify factors affecting detection reproducibility and accuracy in cancer genomes. These large cross-platform/site WGS and WES datasets using well-characterized reference samples will represent a powerful resource for benchmarking NGS technologies, bioinformatics pipelines, and for the cancer genomics studies.

Lake Hillier is a hypersaline lake known for its distinctive bright pink color. The cause of this phenomenon in other hypersaline sites has been attributed to halophiles, Dunaliella, and Salinibacter, however, a systematic analysis of the microbial communities, their functional features, and the prevalence of pigment-producing-metabolisms has not been previously studied. Through metagenomic sequencing and culture-based approaches, our results evidence that Lake Hillier is composed of a diverse set of microorganisms including archaea, bacteria, algae, and viruses. Our data indicate that the microbiome in Lake Hillier is composed of multiple pigment-producer microbes, including Dunaliella, Salinibacter, Halobacillus, Psychroflexus, Halorubrum, many of which are cataloged as polyextremophiles. Additionally, we estimated the diversity of metabolic pathways in the lake and determined that many of these are related to pigment production. We reconstructed complete or partial genomes for 21 discrete bacteria (N = 14) and archaea (N = 7), only 2 of which could be taxonomically annotated to previously observed species. Our findings provide the first metagenomic study to decipher the source of the pink color of Australia's Lake Hillier. The study of this pink hypersaline environment is evidence of a microbial consortium of pigment producers, a repertoire of polyextremophiles, a core microbiome and potentially novel species.

The cancer genome is commonly altered with thousands of structural rearrangements including insertions, deletions, translocation, inversions, duplications, and copy number variations. Thus, structural variant (SV) characterization plays a paramount role in cancer target identification, oncology diagnostics, and personalized medicine. As part of the SEQC2 Consortium effort, the present study established and evaluated a consensus SV call set using a breast cancer reference cell line and matched normal control derived from the same donor, which were used in our companion benchmarking studies as reference samples.

Here, we performed shotgun metagenome sequencing of swab samples collected on floors at a train station in Narita City, Chiba, Japan. The taxonomic analysis revealed that Actinobacteria and Proteobacteria were the dominant phyla. The data will contribute to insight into the microbiome community on the surfaces of urban built environments.

Erythropoiesis involves complex interrelated molecular signals influencing cell survival, differentiation, and enucleation. Diseases associated with ineffective erythropoiesis, such as β-thalassemias, exhibit erythroid expansion and defective enucleation. Clear mechanistic determinants of what make erythropoiesis effective are lacking. We previously demonstrated that exogenous transferrin ameliorates ineffective erythropoiesis in β-thalassemic mice. In the current work, we utilize transferrin treatment to elucidate a molecular signature of ineffective erythropoiesis in β-thalassemia. We hypothesize that compensatory mechanisms are required in β-thalassemic erythropoiesis to prevent apoptosis and enhance enucleation. We identify pleckstrin-2-a STAT5-dependent lipid binding protein downstream of erythropoietin-as an important regulatory node. We demonstrate that partial loss of pleckstrin-2 leads to worsening ineffective erythropoiesis and pleckstrin-2 knockout leads to embryonic lethality in β-thalassemic mice. In addition, the membrane-associated active form of pleckstrin-2 occurs at an earlier stage during β-thalassemic erythropoiesis. Furthermore, membrane-associated activated pleckstrin-2 decreases cofilin mitochondrial localization in β-thalassemic erythroblasts and pleckstrin-2 knockdown in vitro induces cofilin-mediated apoptosis in β-thalassemic erythroblasts. Lastly, pleckstrin-2 enhances enucleation by interacting with and activating RacGTPases in β-thalassemic erythroblasts. This data elucidates the important compensatory role of pleckstrin-2 in β-thalassemia and provides support for the development of targeted therapeutics in diseases of ineffective erythropoiesis.

The emergence and spread of SARS-CoV-2 lineage B.1.1.7, first detected in the United Kingdom, has become a global public health concern because of its increased transmissibility. Over 2500 COVID-19 cases associated with this variant have been detected in the US since December 2020, but the extent of establishment is relatively unknown. Using travel, genomic, and diagnostic data, we highlight the primary ports of entry for B.1.1.7 in the US and locations of possible underreporting of B.1.1.7 cases. Furthermore, we found evidence for many independent B.1.1.7 establishments starting in early December 2020, followed by interstate spread by the end of the month. Finally, we project that B.1.1.7 will be the dominant lineage in many states by mid to late March. Thus, genomic surveillance for B.1.1.7 and other variants urgently needs to be enhanced to better inform the public health response.