Splicing profile by capture RNA-seq identifies pathogenic germline variants in tumor suppressor genes.

Germline variants in tumor suppressor genes (TSGs) can result in RNA mis-splicing and predisposition to cancer. However, identification of variants that impact splicing remains a challenge, contributing to a substantial proportion of patients with suspected hereditary cancer syndromes remaining without a molecular diagnosis. To address this, we used capture RNA-sequencing (RNA-seq) to generate a splicing profile of 18 TSGs (APC, ATM, BRCA1, BRCA2, BRIP1, CDH1, CHEK2, MLH1, MSH2, MSH6, MUTYH, NF1, PALB2, PMS2, PTEN, RAD51C, RAD51D, and TP53) in 345 whole-blood samples from healthy donors. We subsequently demonstrated that this approach can detect mis-splicing by comparing splicing profiles from the control dataset to profiles generated from whole blood of individuals previously identified with pathogenic germline splicing variants in these genes. To assess the utility of our TSG splicing profile to prospectively identify pathogenic splicing variants, we performed concurrent capture DNA and RNA-seq in a cohort of 1000 patients with suspected hereditary cancer syndromes. This approach improved the diagnostic yield in this cohort, resulting in a 9.1% relative increase in the detection of pathogenic variants, demonstrating the utility of performing simultaneous DNA and RNA genetic testing in a clinical context.

Transposable element expression in tumors is associated with immune infiltration and increased antigenicity.

Profound global loss of DNA methylation is a hallmark of many cancers. One potential consequence of this is the reactivation of transposable elements (TEs) which could stimulate the immune system via cell-intrinsic antiviral responses. Here, we develop REdiscoverTE, a computational method for quantifying genome-wide TE expression in RNA sequencing data. Using The Cancer Genome Atlas database, we observe increased expression of over 400 TE subfamilies, of which 262 appear to result from a proximal loss of DNA methylation. The most recurrent TEs are among the evolutionarily youngest in the genome, predominantly expressed from intergenic loci, and associated with antiviral or DNA damage responses. Treatment of glioblastoma cells with a demethylation agent results in both increased TE expression and de novo presentation of TE-derived peptides on MHC class I molecules. Therapeutic reactivation of tumor-specific TEs may synergize with immunotherapy by inducing inflammation and the display of potentially immunogenic neoantigens.

mountainClimber Identifies Alternative Transcription Start and Polyadenylation Sites in RNA-Seq.

Alternative transcription start (ATS) and alternative polyadenylation (APA) create alternative RNA isoforms and modulate many aspects of RNA expression and protein production. However, ATS and APA remain difficult to detect in RNA sequencing (RNA-seq). Here, we developed mountainClimber, a de novo cumulative-sum-based approach to identify ATS and APA as change points. Unlike many existing methods, mountainClimber runs on a single sample and identifies multiple ATS or APA sites anywhere in the transcript. We analyzed 2,342 GTEx samples (36 tissues, 215 individuals) and found that tissue type is the predominant driver of transcript end variations. 75% and 65% of genes exhibited differential APA and ATS across tissues, respectively. In particular, testis displayed longer 5' untranslated regions (UTRs) and shorter 3' UTRs, often in genes related to testis-specific biology. Overall, we report the largest study of transcript ends across human tissues to our knowledge. mountainClimber is available at github.com/gxiaolab/mountainClimber.

Wnt11 regulates cardiac chamber development and disease during perinatal maturation.

Ventricular chamber growth and development during perinatal circulatory transition is critical for functional adaptation of the heart. However, the chamber-specific programs of neonatal heart growth are poorly understood. We used integrated systems genomic and functional biology analyses of the perinatal chamber specific transcriptome and we identified Wnt11 as a prominent regulator of chamber-specific proliferation. Importantly, downregulation of Wnt11 expression was associated with cyanotic congenital heart defect (CHD) phenotypes and correlated with O2 saturation levels in hypoxemic infants with Tetralogy of Fallot (TOF). Perinatal hypoxia treatment in mice suppressed Wnt11 expression and induced myocyte proliferation more robustly in the right ventricle, modulating Rb1 protein activity. Wnt11 inactivation was sufficient to induce myocyte proliferation in perinatal mouse hearts and reduced Rb1 protein and phosphorylation in neonatal cardiomyocytes. Finally, downregulated Wnt11 in hypoxemic TOF infantile hearts was associated with Rb1 suppression and induction of proliferation markers. This study revealed a previously uncharacterized function of Wnt11-mediated signaling as an important player in programming the chamber-specific growth of the neonatal heart. This function influences the chamber-specific development and pathogenesis in response to hypoxia and cyanotic CHDs. Defining the underlying regulatory mechanism may yield chamber-specific therapies for infants born with CHDs.

Decoding the Long Noncoding RNA During Cardiac Maturation: A Roadmap for Functional Discovery.

BACKGROUND: Cardiac maturation during perinatal transition of heart is critical for functional adaptation to hemodynamic load and nutrient environment. Perturbation in this process has major implications in congenital heart defects. Transcriptome programming during perinatal stages is an important information but incomplete in current literature, particularly, the expression profiles of the long noncoding RNAs (lncRNAs) are not fully elucidated. METHODS AND RESULTS: From comprehensive analysis of transcriptomes derived from neonatal mouse heart left and right ventricles, a total of 45 167 unique transcripts were identified, including 21 916 known and 2033 novel lncRNAs. Among these lncRNAs, 196 exhibited significant dynamic regulation along maturation process. By implementing parallel weighted gene co-expression network analysis of mRNA and lncRNA data sets, several lncRNA modules coordinately expressed in a developmental manner similar to protein coding genes, while few lncRNAs revealed chamber-specific patterns. Out of 2262 lncRNAs located within 50 kb of protein coding genes, 5% significantly correlate with the expression of their neighboring genes. The impact of Ppp1r1b-lncRNA on the corresponding partner gene Tcap was validated in cultured myoblasts. This concordant regulation was also conserved in human infantile hearts. Furthermore, the Ppp1r1b-lncRNA/Tcap expression ratio was identified as a molecular signature that differentiated congenital heart defect phenotypes. CONCLUSIONS: The study provides the first high-resolution landscape on neonatal cardiac lncRNAs and reveals their potential interaction with mRNA transcriptome during cardiac maturation. Ppp1r1b-lncRNA was identified as a regulator of Tcap expression, with dynamic interaction in postnatal cardiac development and congenital heart defects.

Global analyses of endonucleolytic cleavage in mammals reveal expanded repertoires of cleavage-inducing small RNAs and their targets.

In mammals, small RNAs are important players in post-transcriptional gene regulation. While their roles in mRNA destabilization and translational repression are well appreciated, their involvement in endonucleolytic cleavage of target RNAs is poorly understood. Very few microRNAs are known to guide RNA cleavage. Endogenous small interfering RNAs are expected to induce target cleavage, but their target genes remain largely unknown. We report a systematic study of small RNA-mediated endonucleolytic cleavage in mouse through integrative analysis of small RNA and degradome sequencing data without imposing any bias toward known small RNAs. Hundreds of small cleavage-inducing RNAs and their cognate target genes were identified, significantly expanding the repertoire of known small RNA-guided cleavage events. Strikingly, both small RNAs and their target sites demonstrated significant overlap with retrotransposons, providing evidence for the long-standing speculation that retrotransposable elements in mRNAs are leveraged as signals for gene targeting. Furthermore, our analysis showed that the RNA cleavage pathway is also present in human cells but affecting a different repertoire of retrotransposons. These results show that small RNA-guided cleavage is more widespread than previously appreciated. Their impact on retrotransposons in non-coding regions shed light on important aspects of mammalian gene regulation.

The Sm protein methyltransferase PRMT5 is not required for primordial germ cell specification in mice.

PRMT5 is a type II protein arginine methyltransferase with roles in stem cell biology, reprograming, cancer and neurogenesis. During embryogenesis in the mouse, it was hypothesized that PRMT5 functions with the master germline determinant BLIMP1 to promote primordial germ cell (PGC) specification. Using a Blimp1-Cre germline conditional knockout, we discovered that Prmt5 has no major role in murine germline specification, or the first global epigenetic reprograming event involving depletion of cytosine methylation from DNA and histone H3 lysine 9 dimethylation from chromatin. Instead, we discovered that PRMT5 functions at the conclusion of PGC reprograming I to promote proliferation, survival and expression of the gonadal germline program as marked by MVH. We show that PRMT5 regulates gene expression by promoting methylation of the Sm spliceosomal proteins and significantly altering the spliced repertoire of RNAs in mammalian embryonic cells and primordial cells.