ERG activates a stem-like proliferation-differentiation program in prostate epithelial cells with mixed basal-luminal identity

To gain insight into how ERG translocations cause prostate cancer, we performed single cell transcriptional profiling of an autochthonous mouse model at an early stage of disease initiation. Despite broad expression of ERG in all prostate epithelial cells, proliferation was enriched in a small, stem-like population with mixed-luminal basal identity (called intermediate cells). Through a series of lineage tracing and primary prostate tissue transplantation experiments, we find that tumor initiating activity resides in a subpopulation of basal cells that co-express the luminal genes Tmprss2 and Nkx3.1 (called BasalLum) but not in the larger population of classical Krt8+ luminal cells. Upon ERG activation, BasalLum cells give rise to the highly proliferative intermediate state, which subsequently transitions to the larger population of Krt8+ luminal cells characteristic of ERG-positive human cancers. Furthermore, this proliferative population is characterized by an ERG-specific chromatin state enriched for NFkB, AP-1, STAT and NFAT binding, with implications for TF cooperativity. The fact that the proliferative potential of ERG is enriched in a small stem-like population implicates the chromatin context of these cells as a critical variable for unmasking its oncogenic activity.

Reprogramming tumor-associated macrophages to outcompete endovascular endothelial progenitor cells and suppress tumor neoangiogenesis

Tumors develop by invoking a supportive environment characterized by aberrant angiogenesis and infiltration of tumor-associated macrophages (TAMs). In a transgenic model of breast cancer, we found that TAMs localized to the tumor parenchyma and were smaller than mammary tissue macrophages. TAMs had low activity of the metabolic regulator mammalian/mechanistic target of rapamycin complex 1 (mTORC1), and depletion of negative regulator of mTORC1 signaling, tuberous sclerosis complex 1 (TSC1), in TAMs inhibited tumor growth in a manner independent of adaptive lymphocytes. Whereas wild-type TAMs exhibited inflammatory and angiogenic gene expression profiles, TSC1-deficient TAMs had a pro-resolving phenotype. TSC1-deficient TAMs relocated to a perivascular niche, depleted protein C receptor (PROCR)-expressing endovascular endothelial progenitor cells, and rectified the hyperpermeable blood vasculature, causing tumor tissue hypoxia and cancer cell death. TSC1-deficient TAMs were metabolically active and effectively eliminated PROCR-expressing endothelial cells in cell competition experiments. Thus, TAMs exhibit a TSC1-dependent mTORC1-low state, and increasing mTORC1 signaling promotes a pro-resolving state that suppresses tumor growth, defining an innate immune tumor suppression pathway that may be exploited for cancer immunotherapy.

Abstract 2048: Mutations in FOXA1 alter chromatin remodeling and cell fate in prostate organoids

Genomic analysis of targeted patient tumor sequencing identified frequent mutations, 41% in prostate cancer (Li, et al., 2020) in the gene FOXA1, a developmentally important pioneer transcription factor (TF) in mammary and prostate tissues. Previous work by our group and others has shown that these FOXA1 mutations alter global chromatin accessibility and promote growth in prostate cells (Adams, 2019), but the underlying molecular details, including the identity of partner TFs, remain unclear. To address this topic, we generated mouse prostate organoids expressing Foxa1 alleles harboring three distinct classes of mutations: (i) overexpression of WT Foxa1 (reflecting focal gene amplification seen in tumors), (ii) a series of mutants within the Wing2region of the forkhead binding domain (FHBD) and (iii) a mutant bearing a stop codon after the FHBD to represent a series of C-terminal truncation mutants. We performed single nucleus multiome sequencing to obtain gene expression (snRNA-seq) and chromatin accessibility (snATAC-seq) readouts from the same individual nuclei. Whereas each Foxa1 mutant has distinct, often mutant-specific features, several themes emerge. These include alterations in the relative proportion of stem-like (L2) luminal cells vs secretory (L1) luminal cells as well as changes in luminal or basal gene signatures, increased androgen receptor signaling output, and enrichment for motifs of distinct classes of partner TFs. For example, cells expressing the truncation mutant show gain in the accessibility of Gata3 and Pou2f1 TF binding motifs, as well as enhanced numbers of L1-like luminal cells. Functional studies demonstrate that Pou2f1 is specifically required for the pro-luminal phenotype in cells expressing the truncation mutant whereas Gata3 plays a more general pro-luminal role. Correlations in motif accessibility and transcription factor expression across single cells further revealed a composite androgen receptor (AR)-FOXA1 motif enriched in the pro-luminal truncation mutant, while the canonical AR motif was enriched in pro-basal cell mutants. Finally, Foxa1 mutants cooperative with Trp53 and Pten loss in orthotopic prostate tumorigenicity assays, most strikingly manifest by reversion of the basal-like features characteristic of Trp53/Pten loss tumors to Ck8+ luminal adenocarcinoma histology, mirroring that seen in Foxa1-mutant human prostate cancer tumors in mice. Thus, mutant Foxa1 alleles cooperate with canonical prostate cancer tumor suppressors and alter the histologic phenotype of prostate cancers in mice through the activation of basal or luminal lineage differentiation programs.

Citation Format: Erik Ladewig, Abbas Nazir, Christina Leslie, Charles Sawyers. Mutations in FOXA1 alter chromatin remodeling and cell fate in prostate organoids [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 2048.

The oncogenic PI3K-induced transcriptomic landscape reveals key functions in splicing and gene expression regulation

The PI3K pathway regulates proliferation, survival, and metabolism and is frequently activated across human cancers. A comprehensive elucidation of how this signaling pathway controls transcriptional and co-transcriptional processes could provide new insights into the key functions of PI3K signaling in cancer. Here, we undertook a transcriptomic approach to investigate genome-wide gene expression and transcription factor (TF) activity changes, as well as splicing and isoform usage dynamics, downstream of PI3K. These analyses uncovered widespread alternatively spliced (AS) isoforms linked to proliferation, metabolism, and splicing in PIK3CA mutant cells, which were reversed by inhibition of PI3Kα. Analysis of paired tumor biopsies from PIK3CA-mutated breast cancer patients undergoing treatment with PI3Kα inhibitors identified widespread splicing alterations that affect specific isoforms in common with the preclinical models, and these alterations, namely PTK2/FRNK and AFMID isoforms, were validated as functional drivers of cancer cell growth or migration. Mechanistically, isoform-specific splicing factors mediated PI3K-dependent RNA splicing. Treatment with splicing inhibitors rendered breast cancer cells more sensitive to the PI3Kα inhibitor alpelisib, resulting in greater growth inhibition than alpelisib alone. This study provides the first comprehensive analysis of widespread splicing alterations driven by oncogenic PI3K in breast cancer. The atlas of PI3K-mediated splicing programs establishes a key role for the PI3K pathway in regulating splicing, opening new avenues for exploiting PI3K signaling as a therapeutic vulnerability in breast cancer.

Genetic mechanisms of HLA-I loss and immune escape in diffuse large B cell lymphoma

Fifty percent of diffuse large B cell lymphoma (DLBCL) cases lack cell-surface expression of the class I major histocompatibility complex (MHC-I), thus escaping recognition by cytotoxic T cells. Here we show that, across B cell lymphomas, loss of MHC-I, but not MHC-II, is preferentially restricted to DLBCL. To identify the involved mechanisms, we performed whole exome and targeted HLA deep-sequencing in 74 DLBCL samples, and found somatic inactivation of B2M and the HLA-I loci in 80% (34 of 42) of MHC-INEG tumors. Furthermore, 70% (22 of 32) of MHC-IPOS DLBCLs harbored monoallelic HLA-I genetic alterations (MHC-IPOS/mono), indicating allele-specific inactivation. MHC-INEG and MHC-IPOS/mono cases harbored significantly higher mutational burden and inferred neoantigen load, suggesting potential coselection of HLA-I loss and sustained neoantigen production. Notably, the analysis of >500,000 individuals across different cancer types revealed common germline HLA-I homozygosity, preferentially in DLBCL. In mice, germinal-center B cells lacking HLA-I expression did not progress to lymphoma and were counterselected in the context of oncogene-driven lymphomagenesis, suggesting that additional events are needed to license immune evasion. These results suggest a multistep process of HLA-I loss in DLBCL development including both germline and somatic events, and have direct implications for the pathogenesis and immunotherapeutic targeting of this disease.

Genomic Alterations in PIK3CA-Mutated Breast Cancer Result in mTORC1 Activation and Limit the Sensitivity to PI3Kα Inhibitors

PI3Kα inhibitors have shown clinical activity in PIK3CA-mutated estrogen receptor-positive (ER⁺) patients with breast cancer. Using whole genome CRISPR/Cas9 sgRNA knockout screens, we identified and validated several negative regulators of mTORC1 whose loss confers resistance to PI3Kα inhibition. Among the top candidates were TSC1, TSC2, TBC1D7, AKT1S1, STK11, MARK2, PDE7A, DEPDC5, NPRL2, NPRL3, C12orf66, SZT2, and ITFG2. Loss of these genes invariably results in sustained mTOR signaling under pharmacologic inhibition of the PI3K-AKT pathway. Moreover, resistance could be prevented or overcome by mTOR inhibition, confirming the causative role of sustained mTOR activity in limiting the sensitivity to PI3Kα inhibition. Cumulatively, genomic alterations affecting these genes are identified in about 15% of PIK3CA-mutated breast tumors and appear to be mutually exclusive. This study improves our understanding of the role of mTOR signaling restoration in leading to resistance to PI3Kα inhibition and proposes therapeutic strategies to prevent or revert this resistance. SIGNIFICANCE: These findings show that genetic lesions of multiple negative regulators of mTORC1 could limit the efficacy of PI3Kα inhibitors in breast cancer, which may guide patient selection strategies for future clinical trials.

FOXA1 Mutations Reveal Distinct Chromatin Profiles and Influence Therapeutic Response in Breast Cancer

Mutations in the pioneer transcription factor FOXA1 are a hallmark of estrogen receptor-positive (ER⁺) breast cancers. Examining FOXA1 in ∼5,000 breast cancer patients identifies several hotspot mutations in the Wing2 region and a breast cancer-specific mutation SY242CS, located in the third β strand. Using a clinico-genomically curated cohort, together with breast cancer models, we find that FOXA1 mutations associate with a lower response to aromatase inhibitors. Mechanistically, Wing2 mutations display increased chromatin binding at ER loci upon estrogen stimulation, and an enhanced ER-mediated transcription without changes in chromatin accessibility. In contrast, SY242CS shows neomorphic properties that include the ability to open distinct chromatin regions and activate an alternative cistrome and transcriptome. Structural modeling predicts that SY242CS confers a conformational change that mediates stable binding to a non-canonical DNA motif. Taken together, our results provide insights into how FOXA1 mutations perturb its function to dictate cancer progression and therapeutic response.

Distinct genomic profile and specific targeted drug responses in adult cerebellar glioblastoma

Background

Despite extensive efforts on the genomic characterization of gliomas, very few studies have reported the genetic alterations of cerebellar glioblastoma (C-GBM), a rare and lethal disease. Here, we provide a systematic study of C-GBM to better understand its specific genomic features.

Methods

We collected a cohort of C-GBM patients and compared patient demographics and tumor pathologies with supratentorial glioblastoma (S-GBM). To uncover the molecular characteristics, we performed DNA and mRNA sequencing and DNA methylation arrays on 19, 6, and 4 C-GBM cases, respectively. Moreover, chemical drug screening was conducted to identify potential therapeutic options for C-GBMs.

Results

Despite differing anatomical origins of C-GBM and S-GBM, neither histological, cytological, nor patient demographics appeared significantly different between the 2 types. However, we observed striking differences in mutational patterns, including frequent alterations of ATRX, PDGFRA, NF1, and RAS and absence of EGFR alterations in C-GBM. These results show a distinct evolutionary path in C-GBM, suggesting specific therapeutic targeted options. Targeted-drug screening revealed that C-GBMs were more responsive to mitogen-activated protein kinase kinase (MEK) inhibitor and resistant to epidermal growth factor receptor inhibitors than S-GBMs. Also, differential expression analysis indicated that C-GBMs may have originated from oligodendrocyte progenitor cells, suggesting that different types of cells can undergo malignant transformation according to their location in brain. Master regulator analysis with differentially expressed genes between C-GBM and proneural S-GBM revealed NR4A1 as a potential therapeutic target.

Conclusions

Our results imply that unique gliomagenesis mechanisms occur in adult cerebellum and new treatment strategies are needed to provide greater therapeutic benefits for C-GBM patients.

Key Points

1. Distinct genomic profiles of 19 adult cerebellar GBMs were characterized. 2. MEK inhibitor was highly sensitive to cerebellar GBM compared with supratentorial GBM. 3. Master regulator analysis revealed NR4A1 as a potential therapeutic target in cerebellar GBM.

Double PIK3CA mutations in cis increase oncogenicity and sensitivity to PI3Kα inhibitors

Activating mutations in PIK3CA are frequent in human breast cancer, and phosphoinositide 3-kinase alpha (PI3Kα) inhibitors have been approved for therapy. To characterize determinants of sensitivity to these agents, we analyzed PIK3CA-mutant cancer genomes and observed the presence of multiple PIK3CA mutations in 12 to 15% of breast cancers and other tumor types, most of which (95%) are double mutations. Double PIK3CA mutations are in cis on the same allele and result in increased PI3K activity, enhanced downstream signaling, increased cell proliferation, and tumor growth. The biochemical mechanisms of dual mutations include increased disruption of p110α binding to the inhibitory subunit p85α, which relieves its catalytic inhibition, and increased p110α membrane lipid binding. Double PIK3CA mutations predict increased sensitivity to PI3Kα inhibitors compared with single-hotspot mutations.

Abstract P3-03-01: Double PIK3CA mutations in cis drive oncogene addiction and enhance sensitivity to PI3K alpha inhibitors in breast cancer

Activating mutations in PIK3CA, the gene coding for the catalytic subunit (p110α) of phosphoinositide-3-kinase (PI3K), are the most frequent oncogenic alterations in estrogen receptor-positive (ER+) breast cancer and are also prevalent in other tumor types. PI3Kα inhibitors including alpelisib have recently been shown to be clinically active in ER+ PIK3CA mutant breast cancer. To characterize determinants of sensitivity to these agents, we undertook a comprehensive analysis of PIK3CA mutant cancer genomes and observed the presence of double PIK3CA mutations in 12-15% of breast cancer and other tumor types. These double PIK3CA mutations are clonal, located in cis on the same allele, and are composed of a single hotspot mutation combined with a recurrent second-site mutation. Double PIK3CA mutations in cis result in increased PI3K activity and downstream signaling together with enhanced cell proliferation and tumor growth compared to single hotspot mutations. The biochemical mechanisms underlying this increased oncogenicity include increased disruption of p110α binding to the inhibitory subunit p85α, which relieves its catalytic inhibition, and increased membrane lipid binding. Finally, these double PIK3CA mutations predict for increased sensitivity to PI3Kα inhibitors compared to single hotspot mutations (e.g. E545K or H1047R) in experimental models and in patients with ER+ PIK3CA mutant metastatic breast cancer from the SANDPIPER randomized phase III clinical trial. These findings implicate double PIK3CA mutations in cis as a novel mechanism of oncogene addiction relative to single hotspot mutations, providing a rationale to develop PI3Kα inhibitors for the therapy of double PIK3CA mutant cancers.

Citation Format: Neil Vasan, Pedram Razavi, Jared L Johnson, Hong Shao, Timothy Wilson, Frauke Schimmoller, Hardik Shah, Alesia Antoine, Erik Ladewig, Alexander Gorelick, Ting-Yu Lin, Eneda Toska, Guotai Xu, Abiha Kazmi, Matthew T Chang, Barry S. Taylor, Maura N. Dickler, Komal Jhaveri, Sarat Chandarlapaty, Raul Rabadan, Ed Reznik, Melissa L Smith, Robert Sebra, Lori Friedman, Lewis C Cantley, Maurizio Scaltriti, José Baselga. Double PIK3CA mutations in cis drive oncogene addiction and enhance sensitivity to PI3K alpha inhibitors in breast cancer [abstract]. In: Proceedings of the 2019 San Antonio Breast Cancer Symposium; 2019 Dec 10-14; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2020;80(4 Suppl):Abstract nr P3-03-01.

Abstract 3917: Double PIK3CA mutations in cis enhance oncogene activation and sensitivity to PI3K alpha inhibitors in breast cancer

Activating mutations in PIK3CA, the gene coding for the catalytic subunit (p110α) of phosphoinositide 3-kinase (PI3K), are the most frequent oncogenic alterations in all cancers, including estrogen receptor-positive (ER+) breast cancer. There are many distinct oncogenic PIK3CA mutations including major hotspot mutations (e.g. E542K, E545K, H1047R) which are predictive of response to PI3Kα inhibitors. While responses to PI3Kα inhibitors can be variable, some patients with single hotspot mutant tumors derive a deep and prolonged clinical benefit. In a comprehensive analysis of the genomic data available in breast cancer, we observed double PIK3CA mutations and investigated their potential biological relevance and potential correlation with sensitivity to PI3Kα inhibitors. We detected double PIK3CA mutations in 10-15% of all PIK3CA-mutant cancers including breast cancer. Double PIK3CA mutations in breast cancer are defined by a pattern of co-expression of a major hotspot mutation associated with a recurrent second-site PIK3CA mutation (‘minor mutation’). We have found that these double PIK3CA mutations are compound, that is in cis on the same allele, using single molecule real time sequencing (SMRT-seq) of fresh breast tumor samples. Double compound PIK3CA mutations result in increased PI3K activity and downstream signaling compared to single hotspot mutants in nontransformed cells and in ER+ breast cancer cells. Double compound mutations also increase cell proliferation and xenograft growth compared to single hotspot mutants. Biochemical experiments using recombinant PI3K compound mutant protein complexes reveal a combination mechanism of action, with simultaneous PI3K destabilization and increased membrane binding compared to single hotspot mutants, leading to increased kinase activity. Importantly, these compound mutations predict for increased sensitivity to PI3Kα inhibitors compared to single hotspot mutants in both preclinical models and also in patients with breast cancer. Together, our data support that double compound PIK3CA mutations enhance the oncogenicity of single hotspot PIK3CA mutations and this increased dependency results in increased sensitivity to PI3Kα inhibitors compared to single hotspot mutations. We propose that double compound activating mutations are an additional mechanism of oncogenic transformation in addition to other known mechanisms such as gene amplification, point mutation, or gene translocation. We speculate that double compound mutant PIK3CA may function as a clinical biomarker of increased sensitivity to PI3K-directed targeted therapies and may improve the therapeutic window of PI3K inhibitors in ER+ breast cancer and other PIK3CA-mutant tumor histologies.

Citation Format: Neil Vasan, Pedram Razavi, Jared Johnson, Hong Shao, Hardik Shah, Alesia Antoine, Erik Ladewig, Alexander Gorelick, Eneda Toska, Guotai Xu, Abiha Kazmi, Matthew T. Chang, Barry S. Taylor, Maura N. Dickler, Komal Jhaveri, Raul Rabadan, Ed Reznik, Melissa L. Smith, Robert Sebra, Lewis C. Cantley, Maurizio Scaltriti, Jose Baselga. Double PIK3CA mutations in cis enhance oncogene activation and sensitivity to PI3K alpha inhibitors in breast cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 3917.

Abstract 4346: PI3K pathway mediated splicing defects in ER+ breast cancers

Activating mutations in PIK3CA, the gene encoding for the catalytic subunit (p100a), are the most common oncogenic alterations in estrogen receptor-positive (ER+). The majority of PIK3CA mutations occur within two hot spots; exons 9 and exon 20 which encode the helical (E545K) and kinase domains (H1047R), respectively. These mutations result in hyperactivation of the PI3K/AKT/mTOR pathway and provide the rationale for the development of inhibitors targeting the PI3K pathway. To this end, PI3K α-specific inhibitors are showing antitumor activity in patients with PIK3CA-mutant, ER-positive breast cancer. Evidence of alternative mRNA regulation and splicing in various cancers has been described in the literature. Although, the majority of events appear to have unknown clinical significance, there is evidence alternative splicing can lead to drug resistance. The role of the PI3K pathway on transcriptional regulation and mRNA processing is not well studied. Through transcriptome analysis and cell assays in mouse MEF and human MCF10A cells we demonstrate splicing defects accrue as a result of the PI3K pathway activation by the PIK3CA H1047R mutation. PI3Ka pathway inhibitors were able to restore wildtype exon inclusion levels in these mutants, suggesting a potential clinical benefit. We propose that PIK3CA H1047R imposes mRNA differential isoform regulation by acting through the most commonly mutated PI3K pathway in ER+ breast cancers and that such mutations are targetable with PI3K pathway inhibitors.

Citation Format: Erik Ladewig, Lauren Fairchild, Maurizio Scaltriti, Christina Leslie, Eneda Toska, Jose Baselga. PI3K pathway mediated splicing defects in ER+ breast cancers [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 4346.

Nuclear Proximity of Mtr4 to RNA Exosome Restricts DNA Mutational Asymmetry

The distribution of sense and antisense strand DNA mutations on transcribed duplex DNA contributes to the development of immune and neural systems along with the progression of cancer. Because developmentally matured B cells undergo biologically programmed strand-specific DNA mutagenesis at focal DNA/RNA hybrid structures, they make a convenient system to investigate strand-specific mutagenesis mechanisms. We demonstrate that the sense and antisense strand DNA mutagenesis at the immunoglobulin heavy chain locus and some other regions of the B cell genome depends upon localized RNA processing protein complex formation in the nucleus. Both the physical proximity and coupled activities of RNA helicase Mtr4 (and senataxin) with the noncoding RNA processing function of RNA exosome determine the strand-specific distribution of DNA mutations. Our study suggests that strand-specific DNA mutagenesis-associated mechanisms will play major roles in other undiscovered aspects of organismic development.

Clonal evolution of glioblastoma under therapy

Glioblastoma (GBM) is the most common and aggressive primary brain tumor. To better understand how GBM evolves, we analyzed longitudinal genomic and transcriptomic data from 114 patients. The analysis shows a highly branched evolutionary pattern in which 63% of patients experience expression-based subtype changes. The branching pattern, together with estimates of evolutionary rate, suggests that relapse-associated clones typically existed years before diagnosis. Fifteen percent of tumors present hypermutation at relapse in highly expressed genes, with a clear mutational signature. We find that 11% of recurrence tumors harbor mutations in LTBP4, which encodes a protein binding to TGF-β. Silencing LTBP4 in GBM cells leads to suppression of TGF-β activity and decreased cell proliferation. In recurrent GBM with wild-type IDH1, high LTBP4 expression is associated with worse prognosis, highlighting the TGF-β pathway as a potential therapeutic target in GBM.

Diversity of miRNAs, siRNAs, and piRNAs across 25 Drosophila cell lines

We expanded the knowledge base for Drosophila cell line transcriptomes by deeply sequencing their small RNAs. In total, we analyzed more than 1 billion raw reads from 53 libraries across 25 cell lines. We verify reproducibility of biological replicate data sets, determine common and distinct aspects of miRNA expression across cell lines, and infer the global impact of miRNAs on cell line transcriptomes. We next characterize their commonalities and differences in endo-siRNA populations. Interestingly, most cell lines exhibit enhanced TE-siRNA production relative to tissues, suggesting this as a common aspect of cell immortalization. We also broadly extend annotations of cis-NAT-siRNA loci, identifying ones with common expression across diverse cells and tissues, as well as cell-restricted loci. Finally, we characterize small RNAs in a set of ovary-derived cell lines, including somatic cells (OSS and OSC) and a mixed germline/somatic cell population (fGS/OSS) that exhibits ping-pong piRNA signatures. Collectively, the ovary data reveal new genic piRNA loci, including unusual configurations of piRNA-generating regions. Together with the companion analysis of mRNAs described in a previous study, these small RNA data provide comprehensive information on the transcriptional landscape of diverse Drosophila cell lines. These data should encourage broader usage of fly cell lines, beyond the few that are presently in common usage.

Comparative analysis of the transcriptome across distant species

The transcriptome is the readout of the genome. Identifying common features in it across distant species can reveal fundamental principles. To this end, the ENCODE and modENCODE consortia have generated large amounts of matched RNA-sequencing data for human, worm and fly. Uniform processing and comprehensive annotation of these data allow comparison across metazoan phyla, extending beyond earlier within-phylum transcriptome comparisons and revealing ancient, conserved features. Specifically, we discover co-expression modules shared across animals, many of which are enriched in developmental genes. Moreover, we use expression patterns to align the stages in worm and fly development and find a novel pairing between worm embryo and fly pupae, in addition to the embryo-to-embryo and larvae-to-larvae pairings. Furthermore, we find that the extent of non-canonical, non-coding transcription is similar in each organism, per base pair. Finally, we find in all three organisms that the gene-expression levels, both coding and non-coding, can be quantitatively predicted from chromatin features at the promoter using a 'universal model' based on a single set of organism-independent parameters.

Functional small RNAs are generated from select miRNA hairpin loops in flies and mammals

In the canonical animal microRNA (miRNA) pathway, Drosha generates ∼60- to 70-nucleotide pre-miRNA hairpins that are cleaved by Dicer into small RNA duplexes that load into Argonaute proteins, which retain a single mature strand in the active complex. The terminal loops of some miRNA hairpins regulate processing efficiency, but once liberated by Dicer, they are generally considered nonfunctional by-products. Here, we show that specific miRNA loops accumulate in effector Argonaute complexes in Drosophila and mediate miRNA-type repression. This was unexpected, since endogenous loading of Argonaute proteins was believed to occur exclusively via small RNA duplexes. Using in vitro assays, which recapitulate Argonaute-specific loop loading from synthetic pre-miRNAs and even single-stranded oligoribonucleotides corresponding to miRNA loops, we reveal that the loop-loading mechanism is distinct from duplex loading. We also show that miRNA loops loaded into the miRNA effector AGO1 are subject to 3' resection, and structure-function analyses indicate selectivity of loop loading. Finally, we demonstrate that select miRNA loops in mammals are similarly loaded into Argonaute complexes and direct target repression. Altogether, we reveal a conserved mechanism that yields functional RNAs from miRNA loop regions, broadening the repertoire of Argonaute-dependent regulatory RNAs and providing evidence for functionality of endogenous ssRNA species.

Discovery of hundreds of mirtrons in mouse and human small RNA data

Atypical miRNA substrates do not fit criteria often used to annotate canonical miRNAs, and can escape the notice of miRNA genefinders. Recent analyses expanded the catalogs of invertebrate splicing-derived miRNAs (“mirtrons”), but only a few tens of mammalian mirtrons have been recognized to date. We performed meta-analysis of 737 mouse and human small RNA data sets comprising 2.83 billion raw reads. Using strict and conservative criteria, we provide confident annotation for 237 mouse and 240 human splicing-derived miRNAs, the vast majority of which are novel genes. These comprise three classes of splicing-derived miRNAs in mammals: conventional mirtrons, 5′-tailed mirtrons, and 3′-tailed mirtrons. In addition, we segregated several hundred additional human and mouse loci with candidate (and often compelling) evidence. Most of these loci arose relatively recently in their respective lineages. Nevertheless, some members in each of the three mirtron classes are conserved, indicating their incorporation into beneficial regulatory networks. We also provide the first Northern validation for mammalian mirtrons, and demonstrate Dicer-dependent association of mature miRNAs from all three classes of mirtrons with Ago2. The recognition of hundreds of mammalian mirtrons provides a new foundation for understanding the scope and evolutionary dynamics of Dicer substrates in mammals.

Common and distinct patterns of terminal modifications to mirtrons and canonical microRNAs

Nucleotide modifications to microRNAs or their precursors can influence their processing and/or activity. A challenge to their analysis is the lack of independent references for the termini generated by primary processing; typically, these are empirically assigned as the most abundant mapped reads. Mirtrons offer such an independent measure since these microRNA hairpins are defined by splicing. Consequently, mirtron-derived reads that deviate from splice sites reflect modification following primary processing. Analysis in Drosophila revealed multiple modification patterns, including select alterations of 5' termini, many 3' resection events, and unexpectedly abundant 3' untemplated monouridylation. Resections occur on mature AGO1-loaded species, whereas uridylation occurs on pre-miRNAs but is compatible with dicing and AGO1 loading. Strikingly, we found many mirtrons whose modified reads are more abundant than those produced by primary processing. In some cases, these abundant modified reads matched the genome owing to fortuitous uridines in downstream flanking exons, thus highlighting the value of an independent reference for the primary-processed sequence. We could further extend the principle of abundant and preferred uridylation of mirtrons, relative to canonical pre-miRNAs, to Caenorhabditis elegans, mouse, and human. Finally, we found that 3' resection occurs broadly across AGO1-loaded canonical miRNA and star species. Altogether, these findings substantially broaden the complexity of terminal modification pathways acting upon small regulatory RNAs.