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.

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 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.

Chemosensing is a primary sense in nature, however little is known about its mechanism in Chelicerata. As a model organism we used the mite Varroa destructor, a key parasite of honeybees. Here we describe a transcriptomic analysis of two physiological stages for the Varroa foreleg, the site of primary olfactory organ. The transcriptomic analysis revealed transcripts of chemosensory related genes belonging to several groups. These include Niemann-Pick disease protein, type C2 (NPC2), gustatory receptors (GRs), ionotropic receptors (IRs), sensory neuron membrane proteins (SNMPs) and odorant binding proteins (OBP). However, no insect odorant receptors (ORs) and odorant co-receptors (ORcos) were found. In addition, we identified a homolog of the most ancient IR co-receptor, IR25a, in Varroa as well as in other members of Acari. High expression of this transcript in the mite's forelegs, while not detectable in the other pairs of legs, suggests a function for this IR25a-like in Varroa chemosensing.

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.

Advancing our understanding of the connections among groundwater, food, and climate is critical to meet global food demands while optimizing water resources usage. However, our understanding of the linkages among groundwater, food, and climate is still limited. Here, we offer a Bayesian framework to simulate crop yield at a regional scale and quantify its relationships and associated uncertainty with climate, groundwater, agricultural, and energy-related variables. We implemented the framework in the rice-producing regions of Louisiana from 1960-2015. To build a parsimonious model, we used a probability-based variable selection approach to detect the key drivers of rice yield. Rice yield increased, groundwater declined, and area planted declined or did not change over 56yrs. The number of irrigation wells, groundwater level, air temperature, and area planted were found to be the key drivers of rice yield. The regression coefficients showed that rice yield was positively related to groundwater level, and negatively related to area planted and the number of irrigation wells. The limited influence of N fertilizer was noted on rice yield for the period when fertilizer data were available. The inverse relationship between rice yield and area planted pointed to the adaption of efficient crop management practices that maintained or increased yield, despite the decline in area planted. The farmers' ability to install irrigation wells during droughts sustained the yields over long-term but not short-term. This decline in rice yield in response to drought over the short-term might explain the negative relation between yield and irrigation wells. Overall, this work highlighted the uncertainty in relationships between rice yield and key drivers and quantified the intimate connection between food and groundwater. This work may have implications for managing two highly competing commodities (i.e., groundwater and food) in agricultural regions.

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 an ongoing microbial tracking investigation of the International Space Station (ISS), several Sphingomonas strains were isolated. Based on the 16S rRNA gene sequence, phylogenetic analysis identified the ISS strains as Sphingomonas sanguinis (n = 2) and one strain isolated from the Kennedy Space Center cleanroom (used to assemble various Mars mission spacecraft components) as Sphingomonas paucimobilis. Metagenomic sequence analyses of different ISS locations identified 23 Sphingomonas species. An abundance of shotgun metagenomic reads were detected for S. sanguinis in the location from where the ISS strains were isolated. A complete metagenome-assembled genome was generated from the shotgun reads metagenome, and its comparison with the whole-genome sequences (WGS) of the ISS S. sanguinis isolates revealed that they were highly similar. In addition to the phylogeny, the WGS of these Sphingomonas strains were compared with the WGS of the type strains to elucidate genes that can potentially aid in plant growth promotion. Furthermore, the WGS comparison of these strains with the well-characterized Sphingomonas sp. LK11, an arid desert strain, identified several genes responsible for the production of phytohormones and for stress tolerance. Production of one of the phytohormones, indole-3-acetic acid, was further confirmed in the ISS strains using liquid chromatography-mass spectrometry. Pathways associated with phosphate uptake, metabolism, and solubilization in soil were conserved across all the S. sanguinis and S. paucimobilis strains tested. Furthermore, genes thought to promote plant resistance to abiotic stress, including heat/cold shock response, heavy metal resistance, and oxidative and osmotic stress resistance, appear to be present in these space-related S. sanguinis and S. paucimobilis strains. Characterizing these biotechnologically important microorganisms found on the ISS and harnessing their key features will aid in the development of self-sustainable long-term space missions in the future. IMPORTANCESphingomonas is ubiquitous in nature, including the anthropogenically contaminated extreme environments. Members of the Sphingomonas genus have been identified as potential candidates for space biomining beyond earth. This study describes the isolation and identification of Sphingomonas members from the ISS, which are capable of producing the phytohormone indole-3-acetic acid. Microbial production of phytohormones will help future in situ studies, grow plants beyond low earth orbit, and establish self-sustainable life support systems. Beyond phytohormone production, stable genomic elements of abiotic stress resistance, heavy metal resistance, and oxidative and osmotic stress resistance were identified, rendering the ISS Sphingomonas isolate a strong candidate for biotechnology-related applications.

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.