5-Formylcytosine is an activating epigenetic mark for RNA Pol III during zygotic reprogramming.

5-Methylcytosine (5mC) is an established epigenetic mark in vertebrate genomic DNA, but whether its oxidation intermediates formed during TET-mediated DNA demethylation possess an instructive role of their own that is also physiologically relevant remains unresolved. Here, we reveal a 5-formylcytosine (5fC) nuclear chromocenter, which transiently forms during zygotic genome activation (ZGA) in Xenopus and mouse embryos. We identify this chromocenter as the perinucleolar compartment, a structure associated with RNA Pol III transcription. In Xenopus embryos, 5fC is highly enriched on Pol III target genes activated at ZGA, notably at oocyte-type tandem arrayed tRNA genes. By manipulating Tet and Tdg enzymes, we show that 5fC is required as a regulatory mark to promote Pol III recruitment as well as tRNA expression. Concordantly, 5fC modification of a tRNA transgene enhances its expression in vivo. The results establish 5fC as an activating epigenetic mark during zygotic reprogramming of Pol III gene expression.

Genome-wide Nucleotide-Resolution Mapping of DNA Replication Patterns, Single-Strand Breaks, and Lesions by GLOE-Seq.

DNA single-strand breaks (SSBs) are among the most common lesions in the genome, arising spontaneously and as intermediates of many DNA transactions. Nevertheless, in contrast to double-strand breaks (DSBs), their distribution in the genome has hardly been addressed in a meaningful way. We now present a technique based on genome-wide ligation of 3'-OH ends followed by sequencing (GLOE-Seq) and an associated computational pipeline designed for capturing SSBs but versatile enough to be applied to any lesion convertible into a free 3'-OH terminus. We demonstrate its applicability to mapping of Okazaki fragments without prior size selection and provide insight into the relative contributions of DNA ligase 1 and ligase 3 to Okazaki fragment maturation in human cells. In addition, our analysis reveals biases and asymmetries in the distribution of spontaneous SSBs in yeast and human chromatin, distinct from the patterns of DSBs.

Extensive nuclear gyration and pervasive non-genic transcription during primordial germ cell development in zebrafish.

Primordial germ cells (PGCs) are the precursors of germ cells, which migrate to the genital ridge during early development. Relatively little is known about PGCs after their migration. We studied this post-migratory stage using microscopy and sequencing techniques, and found that many PGC-specific genes, including genes known to induce PGC fate in the mouse, are only activated several days after migration. At this same time point, PGC nuclei become extremely gyrated, displaying general broad opening of chromatin and high levels of intergenic transcription. This is accompanied by changes in nuage morphology, expression of large loci (PGC-expressed non-coding RNA loci, PERLs) that are enriched for retro-transposons and piRNAs, and a rise in piRNA biogenesis signatures. Interestingly, no nuclear Piwi protein could be detected at any time point, indicating that the zebrafish piRNA pathway is fully cytoplasmic. Our data show that the post-migratory stage of zebrafish PGCs holds many cues to both germ cell fate establishment and piRNA pathway activation.

Signalling codes for the maintenance and lineage commitment of embryonic gastric epithelial progenitors.

The identity of embryonic gastric epithelial progenitors is unknown. We used single-cell RNA-sequencing, genetic lineage tracing and organoid assays to assess whether Axin2- and Lgr5-expressing cells are gastric progenitors in the developing mouse stomach. We show that Axin2+ cells represent a transient population of embryonic epithelial cells in the forestomach. Lgr5+ cells generate both glandular corpus and squamous forestomach organoids ex vivo Only Lgr5+ progenitors give rise to zymogenic cells in culture. Modulating the activity of the WNT, BMP and Notch pathways in vivo and ex vivo, we found that WNTs are essential for the maintenance of Lgr5+ epithelial cells. Notch prevents differentiation of the embryonic epithelial cells along all secretory lineages and hence ensures their maintenance. Whereas WNTs promote differentiation of the embryonic progenitors along the zymogenic cell lineage, BMPs enhance their differentiation along the parietal lineage. In contrast, WNTs and BMPs are required to suppress differentiation of embryonic gastric epithelium along the pit cell lineage. Thus, coordinated action of the WNT, BMP and Notch pathways controls cell fate determination in the embryonic gastric epithelium.

The Release 6 reference sequence of the Drosophila melanogaster genome.

Drosophila melanogaster plays an important role in molecular, genetic, and genomic studies of heredity, development, metabolism, behavior, and human disease. The initial reference genome sequence reported more than a decade ago had a profound impact on progress in Drosophila research, and improving the accuracy and completeness of this sequence continues to be important to further progress. We previously described improvement of the 117-Mb sequence in the euchromatic portion of the genome and 21 Mb in the heterochromatic portion, using a whole-genome shotgun assembly, BAC physical mapping, and clone-based finishing. Here, we report an improved reference sequence of the single-copy and middle-repetitive regions of the genome, produced using cytogenetic mapping to mitotic and polytene chromosomes, clone-based finishing and BAC fingerprint verification, ordering of scaffolds by alignment to cDNA sequences, incorporation of other map and sequence data, and validation by whole-genome optical restriction mapping. These data substantially improve the accuracy and completeness of the reference sequence and the order and orientation of sequence scaffolds into chromosome arm assemblies. Representation of the Y chromosome and other heterochromatic regions is particularly improved. The new 143.9-Mb reference sequence, designated Release 6, effectively exhausts clone-based technologies for mapping and sequencing. Highly repeat-rich regions, including large satellite blocks and functional elements such as the ribosomal RNA genes and the centromeres, are largely inaccessible to current sequencing and assembly methods and remain poorly represented. Further significant improvements will require sequencing technologies that do not depend on molecular cloning and that produce very long reads.

Frequent mutation of histone-modifying genes in non-Hodgkin lymphoma.

Follicular lymphoma (FL) and diffuse large B-cell lymphoma (DLBCL) are the two most common non-Hodgkin lymphomas (NHLs). Here we sequenced tumour and matched normal DNA from 13 DLBCL cases and one FL case to identify genes with mutations in B-cell NHL. We analysed RNA-seq data from these and another 113 NHLs to identify genes with candidate mutations, and then re-sequenced tumour and matched normal DNA from these cases to confirm 109 genes with multiple somatic mutations. Genes with roles in histone modification were frequent targets of somatic mutation. For example, 32% of DLBCL and 89% of FL cases had somatic mutations in MLL2, which encodes a histone methyltransferase, and 11.4% and 13.4% of DLBCL and FL cases, respectively, had mutations in MEF2B, a calcium-regulated gene that cooperates with CREBBP and EP300 in acetylating histones. Our analysis suggests a previously unappreciated disruption of chromatin biology in lymphomagenesis.

Analysis of FOXO1 mutations in diffuse large B-cell lymphoma.

Diffuse large B-cell lymphoma (DLBCL) accounts for 30% to 40% of newly diagnosed lymphomas and has an overall cure rate of approximately 60%. Previously, we observed FOXO1 mutations in non-Hodgkin lymphoma patient samples. To explore the effects of FOXO1 mutations, we assessed FOXO1 status in 279 DLBCL patient samples and 22 DLBCL-derived cell lines. FOXO1 mutations were found in 8.6% (24/279) of DLBCL cases: 92.3% (24/26) of mutations were in the first exon, 46.2% (12/26) were recurrent mutations affecting the N-terminal region, and another 38.5% (10/26) affected the Forkhead DNA binding domain. Recurrent mutations in the N-terminal region resulted in diminished T24 phosphorylation, loss of interaction with 14-3-3, and nuclear retention. FOXO1 mutation was associated with decreased overall survival in patients treated with rituximab, cyclophosphamide, doxorubicin, vincristine, and prednisone (P = .037), independent of cell of origin (COO) and the revised International Prognostic Index (R-IPI). This association was particularly evident (P = .003) in patients in the low-risk R-IPI categories. The independent relationship of mutations in FOXO1 to survival, transcending the prognostic influence of the R-IPI and COO, indicates that FOXO1 mutation is a novel prognostic factor that plays an important role in DLBCL pathogenesis.

Mutational and structural analysis of diffuse large B-cell lymphoma using whole-genome sequencing.

Diffuse large B-cell lymphoma (DLBCL) is a genetically heterogeneous cancer composed of at least 2 molecular subtypes that differ in gene expression and distribution of mutations. Recently, application of genome/exome sequencing and RNA-seq to DLBCL has revealed numerous genes that are recurrent targets of somatic point mutation in this disease. Here we provide a whole-genome-sequencing-based perspective of DLBCL mutational complexity by characterizing 40 de novo DLBCL cases and 13 DLBCL cell lines and combining these data with DNA copy number analysis and RNA-seq from an extended cohort of 96 cases. Our analysis identified widespread genomic rearrangements including evidence for chromothripsis as well as the presence of known and novel fusion transcripts. We uncovered new gene targets of recurrent somatic point mutations and genes that are targeted by focal somatic deletions in this disease. We highlight the recurrence of germinal center B-cell-restricted mutations affecting genes that encode the S1P receptor and 2 small GTPases (GNA13 and GNAI2) that together converge on regulation of B-cell homing. We further analyzed our data to approximate the relative temporal order in which some recurrent mutations were acquired and demonstrate that ongoing acquisition of mutations and intratumoral clonal heterogeneity are common features of DLBCL. This study further improves our understanding of the processes and pathways involved in lymphomagenesis, and some of the pathways mutated here may indicate new avenues for therapeutic intervention.

Aberrant patterns of H3K4 and H3K27 histone lysine methylation occur across subgroups in medulloblastoma.

Recent sequencing efforts have described the mutational landscape of the pediatric brain tumor medulloblastoma. Although MLL2 is among the most frequent somatic single nucleotide variants (SNV), the clinical and biological significance of these mutations remains uncharacterized. Through targeted re-sequencing, we identified mutations of MLL2 in 8 % (14/175) of MBs, the majority of which were loss of function. Notably, we also report mutations affecting the MLL2-binding partner KDM6A, in 4 % (7/175) of tumors. While MLL2 mutations were independent of age, gender, histological subtype, M-stage or molecular subgroup, KDM6A mutations were most commonly identified in Group 4 MBs, and were mutually exclusive with MLL2 mutations. Immunohistochemical staining for H3K4me3 and H3K27me3, the chromatin effectors of MLL2 and KDM6A activity, respectively, demonstrated alterations of the histone code in 24 % (53/220) of MBs across all subgroups. Correlating these MLL2- and KDM6A-driven histone marks with prognosis, we identified populations of MB with improved (K4+/K27-) and dismal (K4-/K27-) outcomes, observed primarily within Group 3 and 4 MBs. Group 3 and 4 MBs demonstrate somatic copy number aberrations, and transcriptional profiles that converge on modifiers of H3K27-methylation (EZH2, KDM6A, KDM6B), leading to silencing of PRC2-target genes. As PRC2-mediated aberrant methylation of H3K27 has recently been targeted for therapy in other diseases, it represents an actionable target for a substantial percentage of medulloblastoma patients with aggressive forms of the disease.

An unbiased AAV-STARR-seq screen revealing the enhancer activity map of genomic regions in the mouse brain in vivo.

Enhancers are important cis-regulatory elements controlling cell-type specific expression patterns of genes. Furthermore, combinations of enhancers and minimal promoters are utilized to construct small, artificial promoters for gene delivery vectors. Large-scale functional screening methodology to construct genomic maps of enhancer activities has been successfully established in cultured cell lines, however, not yet applied to terminally differentiated cells and tissues in a living animal. Here, we transposed the Self-Transcribing Active Regulatory Region Sequencing (STARR-seq) technique to the mouse brain using adeno-associated-viruses (AAV) for the delivery of a highly complex screening library tiling entire genomic regions and covering in total 3 Mb of the mouse genome. We identified 483 sequences with enhancer activity, including sequences that were not predicted by DNA accessibility or histone marks. Characterizing the expression patterns of fluorescent reporters controlled by nine candidate sequences, we observed differential expression patterns also in sparse cell types. Together, our study provides an entry point for the unbiased study of enhancer activities in organisms during health and disease.