Summary By analyzing gene expression data in gliobastoma in combination with matched microRNA profiles, we have uncovered a post-transcriptional regulation layer of surprising magnitude, comprising over 248,000 microRNA (miR)-mediated interactions. These include ~7,000 genes whose transcripts act as miR ‘sponges’ and 148 genes that act through alternative, non-sponge interactions. Biochemical analyses in cell lines confirmed that this network regulates established drivers of tumor initiation and subtype, including PTEN, PDGFRA, RB1, VEGFA, STAT3, and RUNX1, suggesting that these interactions mediate crosstalk between canonical oncogenic pathways. RNA silencing of 13 microRNA-mediated PTEN regulators, whose locus deletions are predictive of PTEN expression variability, was sufficient to downregulate PTEN in a 3′ UTR-dependent manner and to increase tumor-cell growth rates. Thus, this miR-mediated network provides a mechanistic, experimentally validated rationale for the loss of PTEN expression in a large number of glioma samples with an intact PTEN locus.

In mammals, extensive chromatin reorganization is essential for reprogramming terminally committed gametes to a totipotent state during preimplantation development. However, the global chromatin landscape and its dynamics in this period remain unexplored. Here we report a genome-wide map of accessible chromatin in mouse preimplantation embryos using an improved assay for transposase-accessible chromatin with high throughput sequencing (ATAC-seq) approach with CRISPR/Cas9-assisted mitochondrial DNA depletion. We show that despite extensive parental asymmetry in DNA methylomes, the chromatin accessibility between the parental genomes is globally comparable after major zygotic genome activation (ZGA). Accessible chromatin in early embryos is widely shaped by transposable elements and overlaps extensively with putative cis-regulatory sequences. Unexpectedly, accessible chromatin is also found near the transcription end sites of active genes. By integrating the maps of cis-regulatory elements and single-cell transcriptomes, we construct the regulatory network of early development, which helps to identify the key modulators for lineage specification. Finally, we find that the activities of cis-regulatory elements and their associated open chromatin diminished before major ZGA. Surprisingly, we observed many loci showing non-canonical, large open chromatin domains over the entire transcribed units in minor ZGA, supporting the presence of an unusually permissive chromatin state. Together, these data reveal a unique spatiotemporal chromatin configuration that accompanies early mammalian development.

Spermatogenesis is a differentiation process during which diploid spermatogonial stem cells (SSCs) produce haploid spermatozoa. This highly specialized process is precisely controlled at the transcriptional, posttranscriptional, and translational levels. Here we report that N-methyladenosine (mA), an epitranscriptomic mark regulating gene expression, plays essential roles during spermatogenesis. We present comprehensive mA mRNA methylomes of mouse spermatogenic cells from five developmental stages: undifferentiated spermatogonia, type A spermatogonia, preleptotene spermatocytes, pachytene/diplotene spermatocytes, and round spermatids. Germ cell-specific inactivation of the mA RNA methyltransferase Mettl3 or Mettl14 with Vasa-Cre causes loss of mA and depletion of SSCs. mA depletion dysregulates translation of transcripts that are required for SSC proliferation/differentiation. Combined deletion of Mettl3 and Mettl14 in advanced germ cells with Stra8-GFPCre disrupts spermiogenesis, whereas mice with single deletion of either Mettl3 or Mettl14 in advanced germ cells show normal spermatogenesis. The spermatids from double-mutant mice exhibit impaired translation of haploid-specific genes that are essential for spermiogenesis. This study highlights crucial roles of mRNA mA modification in germline development, potentially ensuring coordinated translation at different stages of spermatogenesis.

The intestine plays an important role in nutrient digestion and absorption, microbe defense, and hormone secretion. Although major cell types have been identified in the mouse intestinal epithelium, cell type–specific markers and functional assignments are largely unavailable for human intestine. Here, our single-cell RNA-seq analyses of 14,537 epithelial cells from human ileum, colon, and rectum reveal different nutrient absorption preferences in the small and large intestine, suggest the existence of Paneth-like cells in the large intestine, and identify potential new marker genes for human transient-amplifying cells and goblet cells. We have validated some of these insights by quantitative PCR, immunofluorescence, and functional analyses. Furthermore, we show both common and differential features of the cellular landscapes between the human and mouse ilea. Therefore, our data provide the basis for detailed characterization of human intestine cell constitution and functions, which would be helpful for a better understanding of human intestine disorders, such as inflammatory bowel disease and intestinal tumorigenesis.

R-loops are functional structures in chromatin comprising one single-stranded DNA and a DNA:RNA hybrid. Here, we report ssDRIP-seq, a single-strand DNA ligation-based library preparation technique for genome-wide identification of R-loops. When applied in Arabidopsis, ssDRIP-seq exhibits high efficiency, low bias and strand specificity. We found that Arabidopsis R-loops are enriched by both AT and GC skews, and are formed in the sense and antisense orientations. R-loops are strongly enriched in gene promoters and gene bodies, and are highly associated with noncoding RNA and repetitive genomic regions. Furthermore, R-loops are negatively correlated with CG DNA hypermethylation, and are prevalent in regions with multiple chromatin modifications, showing strong correlations with the activated and repressed gene loci. Our analyses indicate that R-loops are common features in the Arabidopsis genome and suggest that the R-loops play diverse roles in genome organization and gene regulation, thereby providing insights into plant nuclear genome formation and function.

The closely regulated process of mRNA translation is crucial for precise control of protein abundance and quality. Ribosome profiling, a combination of ribosome foot-printing and RNA deep sequencing, has been used in a large variety of studies to quantify genome-wide mRNA translation. Here, we developed Xtail, an analysis pipeline tailored for ribosome profiling data that comprehensively and accurately identifies differentially translated genes in pairwise comparisons. Applied on simulated and real datasets, Xtail exhibits high sensitivity with minimal false-positive rates, outperforming existing methods in the accuracy of quantifying differential translations. With published ribosome profiling datasets, Xtail does not only reveal differentially translated genes that make biological sense, but also uncovers new events of differential translation in human cancer cells on mTOR signalling perturbation and in human primary macrophages on interferon gamma (IFN-γ) treatment. This demonstrates the value of Xtail in providing novel insights into the molecular mechanisms that involve translational dysregulations.

Elevated expression of RNA binding protein HNRNPC has been reported in cancer cells, while the essentialness and functions of HNRNPC in tumors were not clear. We showed that repression of HNRNPC in the breast cancer cells MCF7 and T47D inhibited cell proliferation and tumor growth. Our computational inference of the key pathways and extensive experimental investigations revealed that the cascade of interferon responses mediated by RIG‐I was responsible for such tumor‐inhibitory effect. Interestingly, repression of HNRNPC resulted in accumulation of endogenous double‐stranded RNA (dsRNA), the binding ligand of RIG‐I. These up‐regulated dsRNA species were highly enriched by Alu sequences and mostly originated from pre‐mRNA introns that harbor the known HNRNPC binding sites. Such source of dsRNA is different than the recently well‐characterized endogenous retroviruses that encode dsRNA. In summary, essentialness of HNRNPC in the breast cancer cells was attributed to its function in controlling the endogenous dsRNA and the down‐stream interferon response. This is a novel extension from the previous understandings about HNRNPC in binding with introns and regulating RNA splicing.

Germinal centers (GCs) support high-affinity, long-lived humoral immunity. How memory B cells develop in GCs is not clear. Through the use of a cell-cycle-reporting system, we identified GC-derived memory precursor cells (GC-MP cells) that had quit cycling and reached G0 phase while in the GC, exhibited memory-associated phenotypes with signs of affinity maturation and localized toward the GC border. After being transferred into adoptive hosts, GC-MP cells reconstituted a secondary response like genuine memory B cells. GC-MP cells expressed the interleukin 9 (IL-9) receptor and responded to IL-9. Acute treatment with IL-9 or antibody to IL-9 accelerated or retarded the positioning of GC-MP cells toward the GC edge and exit from the GC, and enhanced or inhibited the development of memory B cells, which required B cell-intrinsic responsiveness to IL-9. Follicular helper T cells (T cells) produced IL-9, and deletion of IL-9 from T cells or, more specifically, from GC T cells led to impaired memory formation of B cells. Therefore, the GC development of memory B cells is promoted by T cell-derived IL-9.