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.

Follicular helper T cells (Tfh cells) are required for T cell help to B cells, and BCL6 is the defining transcription factor of Tfh cells. However, the functions of BCL6 in Tfh cells have largely remained unclear. Here we defined the BCL6 cistrome in primary human germinal center Tfh cells to assess mechanisms of BCL6 regulation of CD4 T cells, comparing and contrasting BCL6 function in T and B cells. BCL6 primarily acts as a repressor in Tfh cells, and BCL6 binding was associated with control of Tfh cell migration and repression of alternative cell fates. Interestingly, although some BCL6-bound genes possessed BCL6 DNA-binding motifs, many BCL6-bound loci were instead characterized by the presence of DNA motifs for AP1 or STAT. AP1 complexes are key positive downstream mediators of TCR signaling and external stimuli. We show that BCL6 can directly bind AP1, and BCL6 depends on AP1 for recruitment to BCL6-binding sites with AP1 motifs, suggesting that BCL6 subverts AP1 activity. These findings reveal that BCL6 has broad and multifaceted effects on Tfh biology and provide insight into how this master regulator mediates distinct cell context-dependent phenotypes.

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

Germinal center (GC) B cells feature repression of many gene enhancers to establish their characteristic transcriptome. Here we show that conditional deletion of Lsd1 in GCs significantly impaired GC formation, associated with failure to repress immune synapse genes linked to GC exit, which are also direct targets of the transcriptional repressor BCL6. We found that BCL6 directly binds LSD1 and recruits it primarily to intergenic and intronic enhancers. Conditional deletion of Lsd1 suppressed GC hyperplasia caused by constitutive expression of BCL6 and significantly delayed BCL6-driven lymphomagenesis. Administration of catalytic inhibitors of LSD1 had little effect on GC formation or GC-derived lymphoma cells. Using a CRISPR-Cas9 domain screen, we found instead that the LSD1 Tower domain was critical for dependence on LSD1 in GC-derived B cells. These results indicate an essential role for LSD1 in the humoral immune response, where it modulates enhancer function by forming repression complexes with BCL6.

Hyperthermophilic methanogens are often H2 limited in hot subseafloor environments, and their survival may be due in part to physiological adaptations to low H2 conditions and interspecies H2 transfer. The hyperthermophilic methanogen Methanocaldococcus jannaschii was grown in monoculture at high (80 to 83 μM) and low (15 to 27 μM) aqueous H2 concentrations and in coculture with the hyperthermophilic H2 producer Thermococcus paralvinellae The purpose was to measure changes in growth and CH4 production kinetics, CH4 fractionation, and gene expression in M. jannaschii with changes in H2 flux. Growth and cell-specific CH4 production rates of M. jannaschii decreased with decreasing H2 availability and decreased further in coculture. However, cell yield (cells produced per mole of CH4 produced) increased 6-fold when M. jannaschii was grown in coculture rather than monoculture. Relative to high H2 concentrations, isotopic fractionation of CO2 to CH4 (εCO2-CH4) was 16‰ larger for cultures grown at low H2 concentrations and 45‰ and 56‰ larger for M. jannaschii growth in coculture on maltose and formate, respectively. Gene expression analyses showed H2-dependent methylene-tetrahydromethanopterin (H4MPT) dehydrogenase expression decreased and coenzyme F420-dependent methylene-H4MPT dehydrogenase expression increased with decreasing H2 availability and in coculture growth. In coculture, gene expression decreased for membrane-bound ATP synthase and hydrogenase. The results suggest that H2 availability significantly affects the CH4 and biomass production and CH4 fractionation by hyperthermophilic methanogens in their native habitats.IMPORTANCE Hyperthermophilic methanogens and H2-producing heterotrophs are collocated in high-temperature subseafloor environments, such as petroleum reservoirs, mid-ocean ridge flanks, and hydrothermal vents. Abiotic flux of H2 can be very low in these environments, and there is a gap in our knowledge about the origin of CH4 in these habitats. In the hyperthermophile Methanocaldococcus jannaschii, growth yields increased as H2 flux, growth rates, and CH4 production rates decreased. The same trend was observed increasingly with interspecies H2 transfer between M. jannaschii and the hyperthermophilic H2 producer Thermococcus paralvinellae With decreasing H2 availability, isotopic fractionation of carbon during methanogenesis increased, resulting in isotopically more negative CH4 with a concomitant decrease in H2-dependent methylene-tetrahydromethanopterin dehydrogenase gene expression and increase in F420-dependent methylene-tetrahydromethanopterin dehydrogenase gene expression. The significance of our research is in understanding the nature of hyperthermophilic interspecies H2 transfer and identifying biogeochemical and molecular markers for assessing the physiological state of methanogens and possible source of CH4 in natural environments.

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.

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.

ALK positive diffuse large B-cell lymphomas (DLBCL) are a distinct lymphoma subtype associated with a poor outcome. Most of them feature a t(2;17) encoding a clathrin (CLTC)-ALK fusion protein. The contribution of deregulated ALK-activity in the pathogenesis and maintenance of these DLBCLs is not yet known. We established and characterized the first CLTC-ALK positive DLBCL cell line (LM1). LM1 formed tumors in NOD-SCID mice. The selective ALK inhibitor NVP-TAE684 inhibited growth of LM1 cells in vitro at nanomolar concentrations. NVP-TAE684 repressed ALK-activated signalling pathways and induced apoptosis of LM1 DLBCL cells. Inhibition of ALK-activity resulted in sustained tumor regression in the xenotransplant tumor model. These data indicate a role of CLTC-ALK in the maintenance of the malignant phenotype thereby providing a rationale therapeutic target for these otherwise refractory tumors.

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.

Transcription factors are frequently altered in leukaemia through chromosomal translocation, mutation or aberrant expression. AML1-ETO, a fusion protein generated by the t(8;21) translocation in acute myeloid leukaemia, is a transcription factor implicated in both gene repression and activation. AML1-ETO oligomerization, mediated by the NHR2 domain, is critical for leukaemogenesis, making it important to identify co-regulatory factors that 'read' the NHR2 oligomerization and contribute to leukaemogenesis. Here we show that, in human leukaemic cells, AML1-ETO resides in and functions through a stable AML1-ETO-containing transcription factor complex (AETFC) that contains several haematopoietic transcription (co)factors. These AETFC components stabilize the complex through multivalent interactions, provide multiple DNA-binding domains for diverse target genes, co-localize genome wide, cooperatively regulate gene expression, and contribute to leukaemogenesis. Within the AETFC complex, AML1-ETO oligomerization is required for a specific interaction between the oligomerized NHR2 domain and a novel NHR2-binding (N2B) motif in E proteins. Crystallographic analysis of the NHR2-N2B complex reveals a unique interaction pattern in which an N2B peptide makes direct contact with side chains of two NHR2 domains as a dimer, providing a novel model of how dimeric/oligomeric transcription factors create a new protein-binding interface through dimerization/oligomerization. Intriguingly, disruption of this interaction by point mutations abrogates AML1-ETO-induced haematopoietic stem/progenitor cell self-renewal and leukaemogenesis. These results reveal new mechanisms of action of AML1-ETO, and provide a potential therapeutic target in t(8;21)-positive acute myeloid leukaemia.

The Mixed Lineage Leukemia (MLL) protein is an important epigenetic regulator required for the maintenance of gene activation during development. MLL chromosomal translocations produce novel fusion proteins that cause aggressive leukemias in humans. Individual MLL fusion proteins have distinct leukemic phenotypes even when expressed in the same cell type, but how this distinction is delineated on a molecular level is poorly understood. Here, we highlight a unique molecular mechanism whereby the RUNX1 gene is directly activated by MLL-AF4 and the RUNX1 protein interacts with the product of the reciprocal AF4-MLL translocation. These results support a mechanism of transformation whereby two oncogenic fusion proteins cooperate by activating a target gene and then modulating the function of its downstream product.

Mutations in epigenetic modifiers and signaling factors often co-occur in myeloid malignancies, including TET2 and NRAS mutations. Concurrent Tet2 loss and NrasG12D expression in hematopoietic cells induced myeloid transformation, with a fully penetrant, lethal chronic myelomonocytic leukemia (CMML), which was serially transplantable. Tet2 loss and Nras mutation cooperatively led to decrease in negative regulators of mitogen-activated protein kinase (MAPK) activation, including Spry2, thereby causing synergistic activation of MAPK signaling by epigenetic silencing. Tet2/Nras double-mutant leukemia showed preferential sensitivity to MAPK kinase (MEK) inhibition in both mouse model and patient samples. These data provide insights into how epigenetic and signaling mutations cooperate in myeloid transformation and provide a rationale for mechanism-based therapy in CMML patients with these high-risk genetic lesions.

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.

The TET2 DNA hydroxymethyltransferase is frequently disrupted by somatic mutations in diffuse large B cell lymphomas (DLBCLs), a tumor that originates from germinal center (GC) B cells. Here, we show that TET2 deficiency leads to DNA hypermethylation of regulatory elements in GC B cells, associated with silencing of the respective genes. This hypermethylation affects the binding of transcription factors including those involved in exit from the GC reaction and involves pathways such as B cell receptor, antigen presentation, CD40, and others. Normal GC B cells manifest a typical hypomethylation signature, which is caused by AID, the enzyme that mediates somatic hypermutation. However, AID-induced demethylation is markedly impaired in TET2-deficient GC B cells, suggesting that AID epigenetic effects are partially dependent on TET2. Last, we find that TET2 mutant DLBCLs also manifest the aberrant TET2-deficient GC DNA methylation signature, suggesting that this epigenetic pattern is maintained during and contributes to lymphomagenesis.

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.

BCOR is a component of a variant Polycomb repressive complex 1 (PRC1.1). PRC1 and PRC2 complexes together constitute a major gene regulatory system critical for appropriate cellular differentiation. The gene is upregulated in germinal center (GC) B cells and mutated in a number of hematologic malignancies. We report BCOR inactivating alterations in 4/7 classic Hodgkin lymphoma (cHL) cell lines, subclonal somatic mutations in Hodgkin and Reed-Sternberg (HRS) cells of 4/10 cHL cases, and deletions in HRS cells of 7/17 primary cHL cases. In mice, conditional loss of Bcor driven by AID-Cre in GC B cells resulted in gene expression changes of 46 genes (>2-fold) including upregulated Lef1 that encodes a transcription factor responsible for establishing T-cell identity and Il9r (interleukin-9 receptor), an important member of the cytokine network in cHL. Our findings suggest a role for BCOR loss in cHL pathogenesis and GC-B cell homeostasis.