Molecular and genetic dissection of recursive splicing

Intronic ratchet points (RPs) are abundant within long introns in the Drosophila genome and consist of juxtaposed splice acceptor and splice donor (SD) sites. Although they appear to encompass zero-nucleotide exons, we recently clarified that intronic recursive splicing (RS) requires a cryptic exon at the RP (an RS-exon), which is subsequently always skipped and thus absent from mRNA. In addition, Drosophila encodes a smaller set of expressed exons bearing features of RS. Here, we investigate mechanisms that regulate the choice between RP and RS-exon SDs. First, analysis of Drosophila RP SD mutants demonstrates that SD competition suppresses inclusion of cryptic exons in endogenous contexts. Second, characterization of RS-exon reporters implicates exonic sequences as influencing choice of RS-exon usage. Using RS-exon swap and mutagenesis assays, we show exonic sequences can determine RS-exon inclusion. Finally, we provide evidence that splicing can suppress utilization of RP SDs to enable RS-exon expression. Overall, multiple factors can influence splicing of Drosophila RS-exons, which usually result in their complete suppression as zero-nucleotide RPs, but occasionally yield translated RS-exons.

Maternal methionine supplementation during gestation alters alternative splicing and DNA methylation in bovine skeletal muscle

Background The evaluation of alternative splicing, including differential isoform expression and differential exon usage, can provide some insights on the transcriptional changes that occur in response to environmental perturbations. Maternal nutrition is considered a major intrauterine regulator of fetal developmental programming. The objective of this study was to assess potential changes in splicing events in the longissimus dorsi muscle of beef calves gestated under control or methionine-rich diets. RNA sequencing and whole-genome bisulfite sequencing were used to evaluate muscle transcriptome and methylome, respectively. Results Alternative splicing patterns were significantly altered by maternal methionine supplementation. Most of the altered genes were directly implicated in muscle development, muscle physiology, ATP activities, RNA splicing and DNA methylation, among other functions. Interestingly, there was a significant association between DNA methylation and differential exon usage. Indeed, among the set of genes that showed differential exon usage, significant differences in methylation level were detected between significant and non-significant exons, and between contiguous and non-contiguous introns to significant exons. Conclusions Overall, our findings provide evidence that a prenatal diet rich in methyl donors can significantly alter the offspring transcriptome, including changes in isoform expression and exon usage, and some of these changes are mediated by changes in DNA methylation.

Characterization of testis-specific LINC01016 gene reveals isoform-specific roles in controlling biological processes

Long noncoding RNAs (lncRNAs) have emerged as critical regulators of biological processes. However, the aberrant expression of an isoform from the same lncRNA gene could lead to RNA with altered functions due to changes in their conformations, leading to diseases. Here, we describe a detailed characterization of the gene which encodes long intergenic non-protein coding RNA 01016 (LINC01016, a.k.a., LncRNA1195) with a focus on its structure, exon usage, and expression in human and macaque tissues. In this study, we show that it is among the highly expressed lncRNAs in the testis, exclusively conserved among non-human primates, suggesting its recent evolution and is expressed and processed into 12 distinct RNAs in testis, cervix, and uterus tissues. Further, we integrate de novo annotation of expressed LINC01016 transcripts and isoform-dependent gene expression analyses to show that human LINC01016 is a multi-exon gene, processed through differential exon usage with isoform-specific roles. Furthermore, in cervical, testicular, and uterine cancers, LINC01016 isoforms are differentially expressed, and their expression is predictive of survival in these cancers. The study has revealed an essential aspect of lncRNA biology, which is rarely associated with coding RNAs that lncRNA genes are precisely processed to generate isoforms with distinct biological roles in specific tissues.

Nanopore sequencing of brain-derived full-length circRNAs reveals circRNA-specific exon usage, intron retention and microexons

Circular RNA (circRNA) is a class of covalently joined non-coding RNAs with functional roles in a wide variety of cellular processes. Their composition shows extensive overlap with exons found in linear mRNAs making it difficult to delineate their composition using short-read RNA sequencing, particularly for long and multi-exonic circRNAs. Here, we use long-read nanopore sequencing of nicked circRNAs (circNick-LRS) and characterize a total of 18,266 and 39,623 circRNAs in human and mouse brain, respectively. We further develop an approach for targeted long-read sequencing of a panel of circRNAs (circPanel-LRS), eliminating the need for prior circRNA enrichment and find >30 circRNA isoforms on average per targeted locus. Our data show that circRNAs exhibit a large number of splicing events such as novel exons, intron retention and microexons that preferentially occur in circRNAs. We propose that altered exon usage in circRNAs may reflect resistance to nonsense-mediated decay in the absence of translation.

The atypical RNA-binding protein Taf15 regulates dorsoanterior neural development through diverse mechanisms in Xenopus tropicalis

ABSTRACT The FET family of atypical RNA-binding proteins includes Fused in sarcoma (FUS), Ewing's sarcoma (EWS) and the TATA-binding protein-associate factor 15 (TAF15). FET proteins are highly conserved, suggesting specialized requirements for each protein. Fus regulates splicing of transcripts required for mesoderm differentiation and cell adhesion in Xenopus, but the roles of Ews and Taf15 remain unknown. Here, we analyze the roles of maternally deposited and zygotically transcribed Taf15, which is essential for the correct development of dorsoanterior neural tissues. By measuring changes in exon usage and transcript abundance from Taf15-depleted embryos, we found that Taf15 may regulate dorsoanterior neural development through fgfr4 and ventx2.1. Taf15 uses distinct mechanisms to downregulate Fgfr4 expression, namely retention of a single intron within fgfr4 when maternal and zygotic Taf15 is depleted, and reduction in the total fgfr4 transcript when zygotic Taf15 alone is depleted. The two mechanisms of gene regulation (post-transcriptional versus transcriptional) suggest that Taf15-mediated gene regulation is target and co-factor dependent, contingent on the milieu of factors that are present at different stages of development.

Altered splicing associated with the pathology of inflammatory bowel disease

Background Aberrant splicing of individual genes is a well-known mechanism promoting pathology for a wide range of conditions, but disease is less commonly attributed to global disruption of exon usage. To explore the possible association of aberrant splicing with inflammatory bowel disease, we developed a pipeline for quantifying transcript abundance and exon inclusion transcriptome-wide and applied it to a dataset of ileal and rectal biopsies, both obtained in duplicate from 34 pediatric or young adult cases of ulcerative colitis and Crohn’s disease. Results Expression and splicing covary to some extent, and eight individuals exhibited aberrant profiles that can be explained by altered ratios of epithelial to stromal and immune cells. Ancestry-related biases in alternative splicing accounting for 5% of the variance were also observed, in part also related to cell-type proportions. In addition, two individuals were identified who had 284 exons with significantly divergent percent spliced in exons, including in the established IBD risk gene CEACAM1, which caused their ileal samples to resemble the rectum. Conclusions These results imply that quantitative differences in splice usage contribute to the pathology of inflammatory bowel disease in a previously unrecognized manner.

Streamlining differential exon and 3′ UTR usage with diffUTR

Background Despite the importance of alternative poly-adenylation and 3′ UTR length for a variety of biological phenomena, there are limited means of detecting UTR changes from standard transcriptomic data. Results We present the diffUTR Bioconductor package which streamlines and improves upon differential exon usage (DEU) analyses, and leverages existing DEU tools and alternative poly-adenylation site databases to enable differential 3′ UTR usage analysis. We demonstrate the diffUTR features and show that it is more flexible and more accurate than state-of-the-art alternatives, both in simulations and in real data. Conclusions diffUTR enables differential 3′ UTR analysis and more generally facilitates DEU and the exploration of their results.

Gene structure, differential exon usage, and expression of the testis long intergenic non-protein coding RNA 1016 in humans reveals isoform-specific roles in controlling biological processes

Long noncoding RNAs (lncRNAs) have emerged as critical regulators of biological processes. The constant expansion of newly-identified lncRNA genes requires that each one be comprehensively annotated to understand its molecular functions. Here, we describe a detailed characterization of the gene which encodes long intergenic non-protein coding RNA 01016 (LINC01016, a.k.a., LncRNA1195) with a focus on its structure, exon usage, and expression in human and macaque tissues. In this study, we show that it is exclusively conserved among non-human primates, suggesting its recent evolution and is expressed and processed into 12 distinct RNAs in testis, cervix, and uterus tissues. Further, we integrate de novo annotation of expressed LINC01016 transcripts and isoform-dependent gene expression analyses to show that human LINC01016 is a multi-exon gene, processed through differential exon usage with isoform-specific functions. Furthermore, in gynecological cancers, such as cervical squamous cell carcinoma and uterine corpus endometrial carcinoma, LINC01016 is downregulated; however, its higher expression is predictive of relapse-free survival in these cancers. Collectively, these analyses reveal that, unlike coding RNAs, lncRNA isoforms are differentially regulated and precisely processed in specific tissues to perform distinct biological roles.

Selective Activation of CNS and Reference PPARGC1A Promoters Is Associated with Distinct Gene Programs Relevant for Neurodegenerative Diseases

The transcriptional regulator peroxisome proliferator activated receptor gamma coactivator 1A (PGC-1α), encoded by PPARGC1A, has been linked to neurodegenerative diseases. Recently discovered CNS-specific PPARGC1A transcripts are initiated far upstream of the reference promoter, spliced to exon 2 of the reference gene, and are more abundant than reference gene transcripts in post-mortem human brain samples. The proteins translated from the CNS and reference transcripts differ only at their N-terminal regions. To dissect functional differences between CNS-specific isoforms and reference proteins, we used clustered regularly interspaced short palindromic repeats transcriptional activation (CRISPRa) for selective endogenous activation of the CNS or the reference promoters in SH-SY5Y cells. Expression and/or exon usage of the targets was ascertained by RNA sequencing. Compared to controls, more differentially expressed genes were observed after activation of the CNS than the reference gene promoter, while the magnitude of alternative exon usage was comparable between activation of the two promoters. Promoter-selective associations were observed with canonical signaling pathways, mitochondrial and nervous system functions and neurological diseases. The distinct N-terminal as well as the shared downstream regions of PGC-1α isoforms affect the exon usage of numerous genes. Furthermore, associations of risk genes of amyotrophic lateral sclerosis and Parkinson’s disease were noted with differentially expressed genes resulting from the activation of the CNS and reference gene promoter, respectively. Thus, CNS-specific isoforms markedly amplify the biological functions of PPARGC1A and CNS-specific isoforms and reference proteins have common, complementary and selective functions relevant for neurodegenerative diseases.