The functional consequences of intron retention: alternative splicing coupled to NMD as a regulator of gene expression

Genome-Wide Analysis of Alternative Splicing during Host-Virus Interactions in Chicken

The chicken is a model animal for the study of evolution, immunity and development. In addition to their use as a model organism, chickens also represent an important agricultural product. Pathogen invasion has already been shown to modulate the expression of hundreds of genes, but the role of alternative splicing in avian virus infection remains unclear. We used RNA-seq data to analyze virus-induced changes in the alternative splicing of Gallus gallus, and found that a large number of alternative splicing events were induced by virus infection both in vivo and in vitro. Virus-responsive alternative splicing events preferentially occurred in genes involved in metabolism and transport. Many of the alternatively spliced transcripts were also expressed from genes with a function relating to splicing or immune response, suggesting a potential impact of virus infection on pre-mRNA splicing and immune gene regulation. Moreover, exon skipping was the most frequent AS event in chickens during virus infection. This is the first report describing a genome-wide analysis of alternative splicing in chicken and contributes to the genomic resources available for studying host-virus interaction in this species. Our analysis fills an important knowledge gap in understanding the extent of genome-wide alternative splicing dynamics occurring during avian virus infection and provides the impetus for the further exploration of AS in chicken defense signaling and homeostasis.

Full-length transcript sequencing of human and mouse cerebral cortex identifies widespread isoform diversity and alternative splicing

Alternative splicing is a post-transcriptional regulatory mechanism producing distinct mRNA molecules from a single pre-mRNA with a prominent role in the development and function of the central nervous system. We used long-read isoform sequencing to generate full-length transcript sequences in the human and mouse cortex. We identify novel transcripts not present in existing genome annotations, including transcripts mapping to putative novel (unannotated) genes and fusion transcripts incorporating exons from multiple genes. Global patterns of transcript diversity are similar between human and mouse cortex, although certain genes are characterized by striking differences between species. We also identify developmental changes in alternative splicing, with differential transcript usage between human fetal and adult cortex. Our data confirm the importance of alternative splicing in the cortex, dramatically increasing transcriptional diversity and representing an important mechanism underpinning gene regulation in the brain. We provide transcript-level data for human and mouse cortex as a resource to the scientific community.

•

There is widespread transcript diversity in the cortex and many novel transcripts

•

Some genes display big differences in isoform number between human and mouse cortex

•

There is evidence of differential transcript usage between human fetal and adult cortex

•

There are many novel isoforms of genes associated with human brain disease

CYP11B1 variants influence skeletal maturation via alternative splicing

We performed genome-wide association study meta-analysis to identify genetic determinants of skeletal age (SA) deviating in multiple growth disorders. The joint meta-analysis (N = 4557) in two multiethnic cohorts of school-aged children identified one locus, CYP11B1 (expression confined to the adrenal gland), robustly associated with SA (rs6471570-A; β = 0.14; P = 6.2 × 10-12). rs6410 (a synonymous variant in the first exon of CYP11B1 in high LD with rs6471570), was prioritized for functional follow-up being second most significant and the one closest to the first intron-exon boundary. In 208 adrenal RNA-seq samples from GTEx, C-allele of rs6410 was associated with intron 3 retention (P = 8.11 × 10-40), exon 4 inclusion (P = 4.29 × 10-34), and decreased exon 3 and 5 splicing (P = 7.85 × 10-43), replicated using RT-PCR in 15 adrenal samples. As CYP11B1 encodes 11-β-hydroxylase, involved in adrenal glucocorticoid and mineralocorticoid biosynthesis, our findings highlight the role of adrenal steroidogenesis in SA in healthy children, suggesting alternative splicing as a likely underlying mechanism.

Genome-Wide Analysis of Light-Regulated Alternative Splicing in Artemisia annua L

Artemisinin is currently the most effective ingredient in the treatment of malaria, which is thus of great significance to study the genetic regulation of Artemisia annua. Alternative splicing (AS) is a regulatory process that increases the complexity of transcriptome and proteome. The most common mechanism of alternative splicing (AS) in plant is intron retention (IR). However, little is known about whether the IR isoforms produced by light play roles in regulating biosynthetic pathways. In this work we would explore how the level of AS in A. annua responds to light regulation. We obtained a new dataset of AS by analyzing full-length transcripts using both Illumina- and single molecule real-time (SMRT)-based RNA-seq as well as analyzing AS on various tissues. A total of 5,854 IR isoforms were identified, with IR accounting for the highest proportion (48.48%), affirming that IR is the most common mechanism of AS. We found that the number of up-regulated IR isoforms (1534/1378, blue and red light, respectively) was more than twice that of down-regulated (636/682) after treatment of blue or red light. In the artemisinin biosynthetic pathway, 10 genes produced 16 differentially expressed IR isoforms. This work demonstrated that the differential expression of IR isoforms induced by light has the potential to regulate sesquiterpenoid biosynthesis. This study also provides high accuracy full-length transcripts, which can be a valuable genetic resource for further research of A. annua, including areas of development, breeding, and biosynthesis of active compounds.

The RNA binding protein Quaking represses splicing of the Fibronectin EDA exon and downregulates the interferon response

Quaking (QKI) controls RNA metabolism in many biological processes including innate immunity, where its roles remain incompletely understood. To illuminate these roles, we performed genome scale transcriptome profiling in QKI knockout cells with or without poly(I:C) transfection, a double-stranded RNA analog that mimics viral infection. Analysis of RNA-sequencing data shows that QKI knockout upregulates genes induced by interferons, suggesting that QKI is an immune suppressor. Furthermore, differential splicing analysis shows that QKI primarily controls cassette exons, and among these events, we noted that QKI silences splicing of the extra domain A (EDA) exon in fibronectin (FN1) transcripts. QKI knockout results in elevated production and secretion of FN1-EDA protein, which is a known activator of interferons. Consistent with an upregulation of the interferon response in QKI knockout cells, our results show reduced production of dengue virus-2 and Japanese encephalitis virus in these cells. In conclusion, we demonstrate that QKI downregulates the interferon system and attenuates the antiviral state.

Recent research progress on circular RNAs: Biogenesis, properties, functions, and therapeutic potential

Circular RNAs (circRNAs), an emerging family member of RNAs, have gained importance in research due to their new functional roles in cellular physiology and disease progression. circRNAs are usually available in a wide range of cells and have shown tissue-specific expression as well as developmental specific expression. circRNAs are characterized by structural stability, conservation, and high abundance in the cell. In this review, we discuss the different models of biogenesis. The properties of circRNAs such as localization, structure and conserved pattern, stability, and expression specificity are also been illustrated. Furthermore, we discuss the biological functions of circRNAs such as microRNA (miRNA) sponging, cell cycle regulation, cell-to-cell communication, transcription regulation, translational regulation, disease diagnosis, and therapeutic potential. Finally, we discuss the recent research progress and future perspective of circRNAs. This review provides an understanding of potential diagnostic markers and the therapeutic potential of circRNAs, which are emerging daily.

Extensive Changes in Transcription Dynamics Reflected on Alternative Splicing Events in Systemic Lupus Erythematosus Patients

In addition to increasing the complexity of the transcriptional output, alternative RNA splicing can lead to the reduction of mRNA translation or the production of non-functional or malfunctional proteins, thus representing a vital component of the gene regulation process. Herein, we set out to detect and characterize alternative splicing events that occur in whole-blood samples of patients with Systemic Lupus Erythematosus (SLE) as compared to healthy counterparts. Through the implementation of a computational pipeline on published RNA-sequencing data, we identified extensive changes in the transcription dynamics affecting a large number of genes. We found a predominance of intron retention events, with the majority introducing premature stop codons, suggestive of gene repression, in both inactive and active SLE patient samples. Alternative splicing affected a distinct set of genes from the ones detected as differentially expressed in the same comparisons, while alternatively spliced genes tended to reside in genome areas associated with increased gene co-expression. Functional analysis of genes affected by alternative splicing pointed towards particular functions related to metabolism and histone acetylation as of potential interest. Together, our findings underline the importance of incorporating alternative splicing analyses in the context of molecular characterization of complex diseases such as SLE.

Genome-Wide Identification and Characterization of Amino Acid Polyamine Organocation Transporter Family Genes Reveal Their Role in Fecundity Regulation in a Brown Planthopper Species (Nilaparvata lugens)

The brown planthopper (BPH), Nilaparvata lugens Stål (Hemiptera:Delphacidae), is one of the most destructive pests of rice worldwide. As a sap-feeding insect, the BPH is incapable of synthesizing several amino acids which are essential for normal growth and development. Therefore, the insects have to acquire these amino acids from dietary sources or their endosymbionts, in which amino acid transporters (AATs) play a crucial role by enabling the movement of amino acids into and out of insect cells. In this study, a common amino acid transporter gene family of amino acid/polyamine/organocation (APC) was identified in BPHs and analyzed. Based on a homology search and conserved functional domain recognition, 20 putative APC transporters were identified in the BPH genome. Molecular trait analysis showed that the verified BPH APC family members were highly variable in protein features, conserved motif distribution patterns, and exon/intron organization. Phylogenetic analysis of five hemipteran species revealed an evolutionary pattern of interfamily conservation and lineage-specific expansion of this gene family. Moreover, stage- and tissue-specific expression analysis revealed diverse expression patterns in the 20 BPH APC transporter genes. Lastly, a potential BPH fecundity regulatory gene of NlAPC09 was identified and shown to participate in the fecundity regulation through the use of quantitative polymerase chain reaction (qPCR) and RNA inference experiments. Our results provide a basis for further functional investigations of APC transporters in BPH.

The RNA-binding protein LUC7L2 mediates MITA/STING intron retention to negatively regulate innate antiviral response

MITA (also known as STING) is an ER-located adaptor protein, which mediates DNA-triggered innate immune response and is critically involved in autoimmune diseases and tumorigenesis. MITA is regulated by post-translational modifications, but how post-transcriptional mechanisms are involved in the regulation of MITA is still largely unknown. Here, we identified the RNA-binding protein LUC7L2 as a negative regulator of DNA virus-triggered innate immune response. LUC7L2-deficient mice exhibited resistance to lethal herpes simplex virus 1 (HSV-1) infection and reduced HSV-1 loads in the brain. Mechanistically, LUC7L2 directly bound to intron 3 of MITA precursor messenger RNA, inhibited its splicing and promoted its nonsense-mediated decay, leading to its downregulation at protein level. LUC7L2-deficient cells had markedly increased MITA level, leading to heightened innate antiviral response. Finally, LUC7L2 was induced following HSV-1 infection. Our findings reveal a feedback negative post-transcriptional regulatory mechanism for regulation of MITA-mediated innate immune response to viral and aberrant cellular DNA.

Absence of Non-Canonical, Inhibitory MYD88 Splice Variants in B Cell Lymphomas Correlates With Sustained NF-κB Signaling

Gain-of-function mutations of the TLR adaptor and oncoprotein MyD88 drive B cell lymphomagenesis via sustained NF-κB activation. In myeloid cells, both short and sustained TLR activation and NF-κB activation lead to the induction of inhibitory MYD88 splice variants that restrain prolonged NF-κB activation. We therefore sought to investigate whether such a negative feedback loop exists in B cells. Analyzing MYD88 splice variants in normal B cells and different primary B cell malignancies, we observed that MYD88 splice variants in transformed B cells are dominated by the canonical, strongly NF-κB-activating isoform of MYD88 and contain at least three novel, so far uncharacterized signaling-competent splice isoforms. Sustained TLR stimulation in B cells unexpectedly reinforces splicing of NF-κB-promoting, canonical isoforms rather than the 'MyD88s', a negative regulatory isoform reported to be typically induced by TLRs in myeloid cells. This suggests that an essential negative feedback loop restricting TLR signaling in myeloid cells at the level of alternative splicing, is missing in B cells when they undergo proliferation, rendering B cells vulnerable to sustained NF-κB activation and eventual lymphomagenesis. Our results uncover MYD88 alternative splicing as an unappreciated promoter of B cell lymphomagenesis and provide a rationale why oncogenic MYD88 mutations are exclusively found in B cells.

Integrative functional genomic analysis of intron retention in human and mouse brain with Alzheimer's disease

Intron retention (IR) has been implicated in the pathogenesis of complex diseases such as cancers; its association with Alzheimer's disease (AD) remains unexplored. We performed genome-wide analysis of IR through integrating genetic, transcriptomic, and proteomic data of AD subjects and mouse models from the Accelerating Medicines Partnership-Alzheimer's Disease project. We identified 4535 and 4086 IR events in 2173 human and 1736 mouse genes, respectively. Quantitation of IR enabled the identification of differentially expressed genes that conventional exon-level approaches did not reveal. There were significant correlations of intron expression within innate immune genes, like HMBOX1, with AD in humans. Peptides with a high probability of translation from intron-retained mRNAs were identified using mass spectrometry. Further, we established AD-specific intron expression Quantitative Trait Loci, and identified splicing-related genes that may regulate IR. Our analysis provides a novel resource for the search for new AD biomarkers and pathological mechanisms.

Reactive astrocytes in ALS display diminished intron retention

Reactive astrocytes are implicated in amyotrophic lateral sclerosis (ALS), although the mechanisms controlling reactive transformation are unknown. We show that decreased intron retention (IR) is common to human-induced pluripotent stem cell (hiPSC)-derived astrocytes carrying ALS-causing mutations in VCP, SOD1 and C9orf72. Notably, transcripts with decreased IR and increased expression are overrepresented in reactivity processes including cell adhesion, stress response and immune activation. This was recapitulated in public-datasets for (i) hiPSC-derived astrocytes stimulated with cytokines to undergo reactive transformation and (ii) in vivo astrocytes following selective deletion of TDP-43. We also re-examined public translatome sequencing (TRAP-seq) of astrocytes from a SOD1 mouse model, which revealed that transcripts upregulated in translation significantly overlap with transcripts exhibiting decreased IR. Using nucleocytoplasmic fractionation of VCP mutant astrocytes coupled with mRNA sequencing and proteomics, we identify that decreased IR in nuclear transcripts is associated with enhanced nonsense mediated decay and increased cytoplasmic expression of transcripts and proteins regulating reactive transformation. These findings are consistent with a molecular model for reactive transformation in astrocytes whereby poised nuclear reactivity-related IR transcripts are spliced, undergo nuclear-to-cytoplasmic translocation and translation. Our study therefore provides new insights into the molecular regulation of reactive transformation in astrocytes.

Identification of Tumor Microenvironment-Related Alternative Splicing Events to Predict the Prognosis of Endometrial Cancer

Background

Endometrial cancer (EC) is one of the most common female malignant tumors. The immunity is believed to be associated with EC patients’ survival, and growing studies have shown that aberrant alternative splicing (AS) might contribute to the progression of cancers.

Methods

We downloaded the clinical information and mRNA expression profiles of 542 tumor tissues and 23 normal tissues from The Cancer Genome Atlas (TCGA) database. ESTIMATE algorithm was carried out on each EC sample, and the OS-related different expressed AS (DEAS) events were identified by comparing the high and low stromal/immune scores groups. Next, we constructed a risk score model to predict the prognosis of EC patients. Finally, we used unsupervised cluster analysis to compare the relationship between prognosis and tumor immune microenvironment.

Results

The prognostic risk score model was constructed based on 16 OS-related DEAS events finally identified, and then we found that compared with high-risk group the OS in the low-risk group was notably better. Furthermore, according to the results of unsupervised cluster analysis, we found that the better the prognosis, the higher the patient’s ESTIMATE score and the higher the infiltration of immune cells.

Conclusions

We used bioinformatics to construct a gene signature to predict the prognosis of patients with EC. The gene signature was combined with tumor microenvironment (TME) and AS events, which allowed a deeper understanding of the immune status of EC patients, and also provided new insights for clinical patients with EC.

Nonsense-mediated RNA decay and its bipolar function in cancer

Nonsense-mediated decay (NMD) was first described as a quality-control mechanism that targets and rapidly degrades aberrant mRNAs carrying premature termination codons (PTCs). However, it was found that NMD also degrades a significant number of normal transcripts, thus arising as a mechanism of gene expression regulation. Based on these important functions, NMD regulates several biological processes and is involved in the pathophysiology of a plethora of human genetic diseases, including cancer. The present review aims to discuss the paradoxical, pro- and anti-tumorigenic roles of NMD, and how cancer cells have exploited both functions to potentiate the disease. Considering recent genetic and bioinformatic studies, we also provide a comprehensive overview of the present knowledge of the advantages and disadvantages of different NMD modulation-based approaches in cancer therapy, reflecting on the challenges imposed by the complexity of this disease. Furthermore, we discuss significant advances in the recent years providing new perspectives on the implications of aberrant NMD-escaping frameshifted transcripts in personalized immunotherapy design and predictive biomarker optimization. A better understanding of how NMD differentially impacts tumor cells according to their own genetic identity will certainly allow for the application of novel and more effective personalized treatments in the near future.

Identification of prognostic and metastasis-related alternative splicing signatures in hepatocellular carcinoma

As the most common neoplasm in digestive system, hepatocellular carcinoma (HCC) is one of the most important leading cause of cancer deaths worldwide. Its high-frequency metastasis and relapse rate lead to the poor survival of HCC patients. However, the mechanism of HCC metastasis is still unclear. Alternative splicing events (ASEs) have a great effect in cancer development, progression and metastasis. We downloaded RNA sequencing and seven types of ASEs data of HCC samples, in order to explore the mechanism of ASEs underlying tumorigenesis and metastasis of HCC. The data were taken from the The Cancer Genome Atlas (TCGA) and TCGASpliceSeq databases. Univariate Cox regression analysis was used to determine a total of 3197 overall survival-related ASEs (OS-SEs). And based on five OS-SEs screened by Lasso regression, we constructed a prediction model with the Area Under Curve of 0.765. With a good reliability of the model, the risk score was also proved to be an independent predictor. Among identified 390 candidate SFs, Y-box protein 3 (YBX3) was significantly correlated with OS and metastasis. Among 177 ASEs, ATP-binding cassette subfamily A member 6 (ABCA6)-43162-AT and PLIN5-46808-AT were identified both associated with OS, bone metastasis and co-expressed with SFs. Then we identified primary bile acid biosynthesis as survival-related (KEGG) pathway by Gene Set Variation Analysis (GSVA) and univariate regression analysis, which was correlated with ABCA6-43162-AT and PLIN5-46808-AT. Finally, we proposed that ABCA6-43162-AT and PLIN5-46808-AT may contribute to HCC poor prognosis and metastasis under the regulation of aberrant YBX3 through the pathway of primary bile acid biosynthesis.

Signature of gene aberrant alternative splicing events in pancreatic adenocarcinoma prognosis

Alternative splicing (AS), as an effective and universal mechanism of transcriptional regulation, is involved in the development and progression of cancer. Therefore, systematic analysis of alternative splicing in pancreatic adenocarcinoma (PAAD) is warranted. The corresponding clinical information of the RNA-Seq data and PAAD cohort was downloaded from the TCGA data portal. Then, a java application, SpliceSeq, was used to evaluate the RNA splicing pattern and calculate the splicing percentage index (PSI). Differentially expressed AS events (DEAS) were identified based on PSI values between PAAD cancer samples and normal samples of adjacent tissues. Kaplan-Meier and Cox regression analyses were used to assess the association between DEAS and patient clinical characteristics. Unsupervised cluster analysis used to reveal four clusters with different survival patterns. At the same time, GEO and TCGA combined with GTEx to verify the differential expression of AS gene and splicing factor. After rigorous filtering, a total of 45,313 AS events were identified, 1,546 of which were differentially expressed AS events. Nineteen DEAS were found to be associated with OS with a five-year overall survival rate of 0.946. And the subtype clusters results indicate that there are differences in the nature of individual AS that affect clinical outcomes. Results also identified 15 splicing factors associated with the prognosis of PAAD. And the splicing factors ESRP1 and RBM5 played an important role in the PAAD-associated AS events. The PAAD-associated AS events, splicing networks, and clusters identified in this study are valuable for deciphering the underlying mechanisms of AS in PAAD and may facilitate the establishment of therapeutic goals for further validation.

Development of alternative splicing signature in lung squamous cell carcinoma

Increasing evidence demonstrated that alternative splicing (AS) plays a vital role in tumorigenesis and clinical outcome of patient. However, systematical analysis of AS in lung squamous cell carcinoma (LUSC) is lacking and greatly necessary. Thus, this study was to systematically estimate the function of AS events served as prognostic indicators in LUSC. Among 31,345 mRNA AS events in 9633 genes, we detected 1996 AS in 1409 genes which have significant connection with overall survival (OS) of LUSC patients. Then, prognostic model based on seven types of AS events was established and we further constructed a combined prognostic model. The Kaplan–Meier curve results suggested that seven types of AS signatures and the combined prognostic model could exhibit robust performance in predicting prognosis. Patients in the high-risk group had significantly shorter OS than those in the low-risk group. The ROC showed all prognostic models had high accuracy and powerful predictive performance with different AUC ranging from 0.837 to 0.978. Moreover, the combined prognostic model had highest performance in risk stratification and predictive accuracy than single prognostic models and had higher accuracy than other mRNA model. Finally, a significant correlation network between survival-related AS genes and prognostic splicing factors (SFs) was established. In conclusion, our study provided several potential prognostic AS models and constructed splicing network between AS and SFs in LUSC, which could be used as potential indicators and treatment targets for LUSC patients.

Direct Nanopore Sequencing of mRNA Reveals Landscape of Transcript Isoforms in Apicomplexan Parasites

Alternative splicing is a widespread phenomenon in metazoans by which single genes are able to produce multiple isoforms of the gene product. However, this has been poorly characterized in apicomplexans, a major phylum of some of the most important global parasites. Efforts have been hampered by atypical transcriptomic features, such as the high AU content of Plasmodium RNA, but also the limitations of short-read sequencing in deciphering complex splicing events. In this study, we utilized the long read direct RNA sequencing platform developed by Oxford Nanopore Technologies to survey the alternative splicing landscape of Toxoplasma gondii and Plasmodium falciparum. We find that while native RNA sequencing has a reduced throughput, it allows us to obtain full-length or nearly full-length transcripts with comparable quantification to Illumina sequencing. By comparing these data with available gene models, we find widespread alternative splicing, particularly intron retention, in these parasites. Most of these transcripts contain premature stop codons, suggesting that in these parasites, alternative splicing represents a pathway to transcriptomic diversity, rather than expanding proteomic diversity. Moreover, alternative splicing rates are comparable between parasites, suggesting a shared splicing machinery, despite notable transcriptomic differences between the parasites. This study highlights a strategy in using long-read sequencing to understand splicing events at the whole-transcript level and has implications in the future interpretation of transcriptome sequencing studies.

IMPORTANCE We have used a novel nanopore sequencing technology to directly analyze parasite transcriptomes. The very long reads of this technology reveal the full-length genes of the parasites that cause malaria and toxoplasmosis. Gene transcripts must be processed in a process called splicing before they can be translated to protein. Our analysis reveals that these parasites very frequently only partially process their gene products, in a manner that departs dramatically from their human hosts.

RUNX1/RUNX1T1 mediates alternative splicing and reorganises the transcriptional landscape in leukemia

The fusion oncogene RUNX1/RUNX1T1 encodes an aberrant transcription factor, which plays a key role in the initiation and maintenance of acute myeloid leukemia. Here we show that the RUNX1/RUNX1T1 oncogene is a regulator of alternative RNA splicing in leukemic cells. The comprehensive analysis of RUNX1/RUNX1T1-associated splicing events identifies two principal mechanisms that underlie the differential production of RNA isoforms: (i) RUNX1/RUNX1T1-mediated regulation of alternative transcription start site selection, and (ii) direct or indirect control of the expression of genes encoding splicing factors. The first mechanism leads to the expression of RNA isoforms with alternative structure of the 5'-UTR regions. The second mechanism generates alternative transcripts with new junctions between internal cassettes and constitutive exons. We also show that RUNX1/RUNX1T1-mediated differential splicing affects several functional groups of genes and produces proteins with unique conserved domain structures. In summary, this study reveals alternative splicing as an important component of transcriptome re-organization in leukemia by an aberrant transcriptional regulator.

Recent trends in cancer therapy: A review on the current state of gene delivery

Cancer treatment has been always considered one of the most critical and vital themes of clinical issues. Many approaches have been developed, depending on the type and the stage of tumor. Gene therapy has the potential to revolutionize different cancer therapy. With the advent of recent bioinformatics technologies and genetic science, it become possible to identify, diagnose and determine the potential treatment using the technology of gene delivery. Several approaches have been developed and experimented in vitro and vivo for cancer therapy including: naked nucleic acids based therapy, targeting micro RNAs, oncolytic virotherapy, suicide gene based therapy, targeting telomerase, cell mediated gene therapy, and CRISPR/Cas9 based therapy. In this review, we present a straightforward introduction to cancer biology and occurrence, highlighting different viral and non-viral gene delivery systems for gene therapy and critically discussed the current and various strategies for cancer gene therapy.

Alternative splicing in aging and Alzheimer's disease: Highlighting the role of tau and estrogen receptor α isoforms in the hypothalamus

Human genes show the highest efficacy of alternative splicing (AS) in the brain as compared to other tissues. Within the brain, a remarkably rich diversity of AS events was identified in the hypothalamus. The AS frequency is increased in the aging brain. Such AS events, as intron retention and accumulation of circular RNAs, were acknowledged as some of the main hallmarks of the aging brain. In Alzheimer's disease (AD) pivotal (tau gene, in particular), risk, candidate and other genes show significant alterations in AS. Therefore AD has been suggested to be a disease of dysregulated AS. One of the reported risk factors for AD is estrogen deficiency that may interfere with the extension of neurobrillary tangles. Mounting evidence suggests that estrogens may decrease hyperphosphorylated tau deposition in the brain. Furthermore, AS of estrogen receptor α (ERα) mRNA is decreased in AD brain areas with the highest tau load. These potential interactions among tau, estrogens, and ERα AS may be important for the development of therapeutic and preventive strategies for AD. The intriguing point is that the amount of splice variants of ERα in the hypothalamus and the hippocampus is increased in aging and decreased in AD, while ERα is one of the regulators of AS and is subject to AS itself.

Perspective in Alternative Splicing Coupled to Nonsense-Mediated mRNA Decay

Alternative splicing (AS) of precursor mRNA (pre-mRNA) is a cellular post-transcriptional process that generates protein isoform diversity. Nonsense-mediated RNA decay (NMD) is an mRNA surveillance pathway that recognizes and selectively degrades transcripts containing premature translation-termination codons (PTCs), thereby preventing the production of truncated proteins. Nevertheless, NMD also fine-tunes the gene expression of physiological mRNAs encoding full-length proteins. Interestingly, around one third of all AS events results in PTC-containing transcripts that undergo NMD. Numerous studies have reported a coordinated action between AS and NMD, in order to regulate the expression of several genes, especially those coding for RNA-binding proteins (RBPs). This coupling of AS to NMD (AS-NMD) is considered a gene expression tool that controls the ratio of productive to unproductive mRNA isoforms, ultimately degrading PTC-containing non-functional mRNAs. In this review, we focus on the mechanisms underlying AS-NMD, and how this regulatory process is able to control the homeostatic expression of numerous RBPs, including splicing factors, through auto- and cross-regulatory feedback loops. Furthermore, we discuss the importance of AS-NMD in the regulation of biological processes, such as cell differentiation. Finally, we analyze interesting recent data on the relevance of AS-NMD to human health, covering its potential roles in cancer and other disorders.

Exploration of predictive and prognostic alternative splicing signatures in lung adenocarcinoma using machine learning methods

Background

Alternative splicing (AS) plays critical roles in generating protein diversity and complexity. Dysregulation of AS underlies the initiation and progression of tumors. Machine learning approaches have emerged as efficient tools to identify promising biomarkers. It is meaningful to explore pivotal AS events (ASEs) to deepen understanding and improve prognostic assessments of lung adenocarcinoma (LUAD) via machine learning algorithms.

Method

RNA sequencing data and AS data were extracted from The Cancer Genome Atlas (TCGA) database and TCGA SpliceSeq database. Using several machine learning methods, we identified 24 pairs of LUAD-related ASEs implicated in splicing switches and a random forest-based classifiers for identifying lymph node metastasis (LNM) consisting of 12 ASEs. Furthermore, we identified key prognosis-related ASEs and established a 16-ASE-based prognostic model to predict overall survival for LUAD patients using Cox regression model, random survival forest analysis, and forward selection model. Bioinformatics analyses were also applied to identify underlying mechanisms and associated upstream splicing factors (SFs).

Results

Each pair of ASEs was spliced from the same parent gene, and exhibited perfect inverse intrapair correlation (correlation coefficient = − 1). The 12-ASE-based classifier showed robust ability to evaluate LNM status of LUAD patients with the area under the receiver operating characteristic (ROC) curve (AUC) more than 0.7 in fivefold cross-validation. The prognostic model performed well at 1, 3, 5, and 10 years in both the training cohort and internal test cohort. Univariate and multivariate Cox regression indicated the prognostic model could be used as an independent prognostic factor for patients with LUAD. Further analysis revealed correlations between the prognostic model and American Joint Committee on Cancer stage, T stage, N stage, and living status. The splicing network constructed of survival-related SFs and ASEs depicts regulatory relationships between them.

Conclusion

In summary, our study provides insight into LUAD researches and managements based on these AS biomarkers.

Yellow-fruited phenotype is caused by 573 bp insertion at 5' UTR of YFT1 allele in yft1 mutant tomato

The yft1 tomato mutant has a yellow-fruited phenotype controlled by a recessive gene of YFT1 allele, which has been shown by map-based cloning to be a homolog of ETHYLENE INSENSITIVE 2 (EIN2). Genetic lesion of YFT1 allele in yft1 is attributed to a 573 bp DNA fragment (IF₅₇₃) insertion at 1,200 bp downstream of the transcription start site. Transcriptomic analysis revealed that YFT1 lesion resulted in 5,053 differentially expressed genes (DEGs) in yft1 pericarp compared with the M82 wild type cultivar. These were annotated as being involved in ethylene synthesis, chromoplast development, and carotenoid synthesis. The YFT1 lesion caused a reduction in its own transcript levels in yft1 and impaired ethylene emission and signal transduction, delayed chromoplast development and decreased carotenoid accumulation. The molecular mechanism underlying the downregulated YFT1 allele in yft1 was examined at both RNA and DNA levels. The IF₅₇₃ event appeared to introduce two negative regulatory sequences located at -272 to -173 bp and -172 to -73 bp in the YFT1 allele promoter, causing alterative splicing due to introduction of aberrant splicing sites, and breaking upstream open reading frames (uORF) structure in the 5'-UTR. Those results a new provided insight into molecular regulation of color formation in tomato fruit.

Widespread intron retention impairs protein homeostasis in C9orf72 ALS brains

The GGGGCC hexanucleotide expansion in C9orf72 (C9) is the most frequent known cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD), yet a clear understanding of how C9 fits into the broader context of ALS/FTD pathology has remained lacking. The repetitive RNA derived from the C9 repeat is known to sequester hnRNPH, a splicing regulator, into insoluble aggregates, resulting in aberrant alternative splicing. Furthermore, hnRNPH insolubility and altered splicing of a robust set of targets have been observed to correlate in C9 and sporadic ALS/FTD patients alike, suggesting that changes along this axis are a core feature of disease pathogenesis. Here, we characterize previously uncategorized RNA splicing defects involving widespread intron retention affecting almost 2000 transcripts in C9ALS/FTD brains exhibiting a high amount of sequestered, insoluble hnRNPH. These intron retention events appear not to alter overall expression levels of the affected transcripts but rather the protein-coding regions. These retained introns affect transcripts in multiple cellular pathways predicted to be involved in C9 as well as sporadic ALS/FTD etiology, including the proteasomal and autophagy systems. The retained intron pre-mRNAs display a number of characteristics, including enrichment of hnRNPH-bound splicing enhancer motifs and a propensity for G-quadruplex (G-Q) formation, linking the defective splicing directly to high amounts of sequestered hnRNPH. Together, our results reveal previously undetected splicing defects in high insoluble hnRNPH-associated C9ALS brains, suggesting a feedback between effective RNA-binding protein dosage and protein quality control in C9, and perhaps all, ALS/FTD.

Holding on to Junk Bonds: Intron Retention in Cancer and Therapy

Intron retention (IR) in cancer was for a long time overlooked by the scientific community, as it was previously considered to be an artifact of a dysfunctional spliceosome. Technological advancements made in the last decade offer unique opportunities to explore the role of IR as a widespread phenomenon that contributes to the transcriptional diversity of many cancers. Numerous studies in cancer have shed light on dysregulation of cellular mechanisms that lead to aberrant and pathologic IR. IR is not merely a mechanism of gene regulation, but rather it can mediate cancer pathogenesis and therapeutic resistance in various human diseases. The burden of IR in cancer is governed by perturbations to mechanisms known to regulate this phenomenon and include epigenetic variation, mutations within the gene body, and splicing factor dysregulation. This review summarizes possible causes for aberrant IR and discusses the role of IR in therapy or as a consequence of disease treatment. As neoepitopes originating from retained introns can be presented on the cancer cell surface, the development of personalized cancer vaccines based on IR-derived neoepitopes should be considered. Ultimately, a deeper comprehension about the origins and consequences of aberrant IR may aid in the development of such personalized cancer vaccines.

Identification of Survival-Associated Alternative Splicing Signatures in Lung Squamous Cell Carcinoma

Purpose: Alternative splicing (AS) is a post-transcriptional process that plays a significant role in enhancing the diversity of transcription and protein. Accumulating evidences have demonstrated that dysregulation of AS is associated with oncogenic processes. However, AS signature specifically in lung squamous cell carcinoma (LUSC) remains unknown. This study aimed to evaluate the prognostic values of AS events in LUSC patients.

Methods: The RNA-seq data, AS events data and corresponding clinical information were obtained from The Cancer Genome Atlas (TCGA) database. Univariate Cox regression analysis was performed to identify survival-related AS events and survival-related parent genes were subjected to Gene Ontology enrichment analysis and gene network analysis. The least absolute shrinkage and selection operator (LASSO) method and multivariate Cox regression analysis were used to construct prognostic prediction models, and their predictive values were assessed by Kaplan-Meier analysis and receiver operating characteristic (ROC) curves. Then a nomogram was established to predict the survival of LUSC patients. And the interaction network of splicing factors (SFs) and survival-related AS events was constructed by Spearman correlation analysis and visualized by Cytoscape.

Results: Totally, 467 LUSC patients were included in this study and 1,991 survival-related AS events within 1,433 genes were identified. SMAD4, FOS, POLR2L, and RNPS1 were the hub genes in the gene interaction network. Eight prognostic prediction models (seven types of AS and all AS) were constructed and all exhibited high efficiency in distinguishing good or poor survival of LUSC patients. The final integrated prediction model including all types of AS events exhibited the best prognostic power with the maximum AUC values of 0.778, 0.816, 0.814 in 1, 3, 5 years ROC curves, respectively. Meanwhile, the nomogram performed well in predicting the 1-, 3-, and 5-year survival of LUSC patients. In addition, the SF-AS regulatory network uncovered a significant correlation between SFs and survival-related AS events.

Conclusion: This is the first comprehensive study to analyze the role of AS events in LUSC specifically, which improves our understanding of the prognostic value of survival-related AS events for LUSC. And these survival-related AS events might serve as novel prognostic biomarkers and drug therapeutic targets for LUSC.

Splicing to Keep Cycling: The Importance of Pre-mRNA Splicing during the Cell Cycle

Pre-mRNA splicing is a fundamental process in mammalian gene expression, and alternative splicing plays an extensive role in generating protein diversity. Because the majority of genes undergo pre-mRNA splicing, most cellular processes depend on proper spliceosome function. We focus on the cell cycle and describe its dependence on pre-mRNA splicing and accurate alternative splicing. We outline the key cell-cycle factors and their known alternative splicing isoforms. We discuss different levels of pre-mRNA splicing regulation such as post-translational modifications and changes in the expression of splicing factors. We describe the effect of chromatin dynamics on pre-mRNA splicing during the cell cycle. In addition, we focus on spliceosome component SF3B1, which is mutated in many types of cancer, and describe the link between SF3B1 and its inhibitors and the cell cycle.

A new mechanism for a familiar mutation - bovine DGAT1 K232A modulates gene expression through multi-junction exon splice enhancement

BACKGROUND

The DGAT1 gene encodes an enzyme responsible for catalysing the terminal reaction in mammary triglyceride synthesis, and underpins a well-known pleiotropic quantitative trait locus (QTL) with a large influence on milk composition phenotypes. Since first described over 15 years ago, a protein-coding variant K232A has been assumed as the causative variant underlying these effects, following in-vitro studies that demonstrated differing levels of triglyceride synthesis between the two protein isoforms.

RESULTS

We used a large RNAseq dataset to re-examine the underlying mechanisms of this large milk production QTL, and hereby report novel expression-based functions of the chr14 g.1802265AA > GC variant that encodes the DGAT1 K232A substitution. Using expression QTL (eQTL) mapping, we demonstrate a highly-significant mammary eQTL for DGAT1, where the K232A mutation appears as one of the top associated variants for this effect. By conducting in vitro expression and splicing experiments in bovine mammary cell culture, we further show modulation of splicing efficiency by this mutation, likely through disruption of an exon splice enhancer as a consequence of the allele encoding the 232A variant.

CONCLUSIONS

The relative contributions of the enzymatic and transcription-based mechanisms now attributed to K232A remain unclear; however, these results suggest that transcriptional impacts contribute to the diversity of lactation effects observed at the DGAT1 locus.

Features of alternative splicing in stomach adenocarcinoma and their clinical implication: a research based on massive sequencing data

BACKGROUND

Alternative splicing (AS) offers a main mechanism to form protein polymorphism. A growing body of evidence indicates the correlation between splicing disorders and carcinoma. Nevertheless, an overall analysis of AS signatures in stomach adenocarcinoma (STAD) is absent and urgently needed.

RESULTS

2042 splicing events were confirmed as prognostic molecular events. Furthermore, the final prognostic signature constructed by 10 AS events gave good result with an area under the curve (AUC) of receiver operating characteristic (ROC) curve up to 0.902 for 5 years, showing high potency in predicting patient outcome. We built the splicing regulatory network to show the internal regulation mechanism of splicing events in STAD. QKI may play a significant part in the prognosis induced by splicing events.

CONCLUSIONS

In our study, a high-efficiency prognostic prediction model was built for STAD patients, and the results showed that AS events could become potential prognostic biomarkers for STAD. Meanwhile, QKI may become an important target for drug design in the future.

Alternative splicing and cancer: insights, opportunities, and challenges from an expanding view of the transcriptome

Over the past decade there has been increased awareness of the potential role of alternative splicing in the etiology of cancer. In particular, advances in RNA-Sequencing technology and analysis has led to a wave of discoveries in the last few years regarding the causes and functional relevance of alternative splicing in cancer. Here we discuss the current understanding of the connections between splicing and cancer, with a focus on the most recent findings. We also discuss remaining questions and challenges that must be addressed in order to use our knowledge of splicing to guide the diagnosis and treatment of cancer.

Development of a Nomogram With Alternative Splicing Signatures for Predicting the Prognosis of Glioblastoma: A Study Based on Large-Scale Sequencing Data

Purpose: Alternative splicing (AS) was reported to play a vital role in development and progression of glioblastoma (GBM), the most common and fatal brain tumor. Systematic analysis of survival-associated AS event profiles and prognostic prediction model based on multiple AS events in GBM was needed. Methods: Genome-wide AS and RNA sequencing profiles were generated in 152 patients with GBM in the cancer genome atlas (TCGA). Prognosis-associated AS events were screened by integrated Cox regression analysis to construct the prognostic risk score model in the training cohort (n = 101). The AS-based signature and clinicopathologic parameters were applied to construct a prognostic nomogram for 0.5-, 1-, and 3-year OS prediction. Finally, the regulatory networks between prognostic AS events and splicing factors (SFs) were constructed. Results: A total of 1,598 prognosis-related AS events from 1,183 source genes were determined. Eight prognostic risk score model based on integrated AS events and 7 AS types were established, respectively. Concordance index (C-index) and receiver operating characteristic (ROC) curve analysis demonstrated powerful ability in distinguishing patients' outcomes. Only Alternate Donor site (AD) and Exon Skip (ES) signature out of the eight types of AS signature were identified as independent prognostic factors for GBM, which was validated in the internal validation cohort. The nomogram with age, new event, pharmaceutical therapy, radiation therapy, AD signature and ES signature were constructed, with C-index of 0.892 (95% CI, 0.853-0.931; P = 5.13 × 10^-15). Calibration plots, ROC, and decision curve analysis suggested excellent predictive performance for the nomogram in both TCGA training cohort and validation cohort. Splicing network indicated distinguished correlations between prognostic AS events and SFs in GBM patients. Conclusions: AS-based prediction model could serve as a promising prognostic predictor and potential therapeutic target for GBM, facilitating better treatment strategies in clinical practice.

Systematic Profiling of Alternative Splicing for Sarcoma Patients Reveals Novel Prognostic Biomarkers Associated with Tumor Microenvironment and Immune Cells

Background

Alternative splicing (AS) events is a novel biomarker of tumor prognosis, but the role of AS events in sarcoma patients remains unclear.

Material/Methods

RNA-seq and clinicopathologic data of the sarcoma cohort were extracted from the TCGA database and data on AS events were downloaded from the TCGASpliceSeq database. Univariate Cox analysis, LASSO regression analysis, and multivariate Cox analysis were performed to determine the overall survival (OS)- and disease-free survival (DFS)-related AS events. Two nomograms were developed based on the independent variables, and subgroup analysis was performed. The area under the curve (AUC), calibration curve, and decision curve analysis (DCA) were used to evaluate the nomograms. Then, we used the CIBERSORT and ESTIMATE package to determine the immune cell proportion and tumor microenvironment (TME) score, respectively. The associations between AS events-based clusters and TME and immune cells were studied.

Results

We identified 1945 and 1831 AS events as OS- and DFS-related AS events, respectively. Two nomograms based on the AS events and clinical data were established and the AUCs of nomograms ranged from 0.807 to 0.894. The calibration curve and DCA showed excellent performance of nomograms. In addition, the results indicated the distinct relationships between AS events-based clusters and OS, DFS, immune score, stromal score, and 10 immune cells.

Conclusions

Our study indicated that AS events are novel prognostic biomarkers for sarcoma patients that may be associated with the TME and immune cells.

Genome-wide analyses of the prognosis-related mRNA alternative splicing landscape and novel splicing factors based on large-scale low grade glioma cohort

Alternative splicing (AS) changes are considered to be critical in predicting treatment response. Our study aimed to investigate differential splicing patterns and to elucidate the role of splicing factor (SF) as prognostic markers of low-grade glioma (LGG). We downloaded RNA-seq data from a cohort of 516 LGG tumors in The Cancer Genome Atlas and analyzed independent prognostic factors using LASSO regression and Cox proportional regression to build a network based on the correlation between SF-related survival AS events. We collected 100 patients from our center for immunohistochemistry and analyzed survival using χ2 test and Cox and Kaplan-Meier analyses. A total of 9,616 AS events related to LGG were screened and identified as well as established related models. Through analyzing specific splicing patterns in LGG, we screened 16 genes to construct a prognostic model to stratify the risk of LGG patients. Validation revealed that the expression level of the prognostic model in LGG tissue was increased, and patients with high expression showed worse prognosis. In summary, we demonstrated the role of SFs and AS events in the progression of LGG, which may provide insights into the clinical significance and aid the future exploration of LGG-associated AS.

Prognostic Value and Potential Role of Alternative mRNA Splicing Events in Cervical Cancer

Background

Increasing evidence suggests that aberrant alternative splicing (AS) events are associated with progression of cancer. This study evaluated the prognostic value and clarify the role of AS events in cervical cancer (CC).

Methods

Based on RNA-seq AS event data and clinical information of CC patients in The Cancer Genome Atlas (TCGA) database, we sought to identify prognosis-related AS events in this setting. We selected several survival-associated AS events to construct a prognostic predictor for CC through the least absolute shrinkage and selection operator (LASSO) and multivariate Cox regression. Moreover, Kyoto Encyclopedia of Genes and Genomes and Gene Ontology analyses were performed on genes with prognosis-related AS events and constructed an AS-splicing factors (SFs) regulatory network.

Results

2770 AS events were significantly correlated with overall survival (OS). The area under the curve (AUC) values of receiver-operator characteristic curve (ROC) for the final prognostic predictor were 0.926, 0.946 and 0.902 at 3, 5, and 10 years, respectively. These values indicated efficiency in prognostic risk stratification for patients with CC. The final prognostic predictor was an independent predictor of OS (HR: 1.24; 95% CI: 1.020–1.504; P < 0.05). The AS-SFs correlation network may reveal an underlying regulatory mechanism of AS events.

Conclusion

AS events are essential participants in the prognosis of CC and hold great potentials for the prognostic stratification and development of treatment strategy.

Intron Retention as a Mode for RNA-Seq Data Analysis

Intron retention (IR) is an alternative splicing mode whereby introns, rather than being spliced out as usual, are retained in mature mRNAs. It was previously considered a consequence of mis-splicing and received very limited attention. Only recently has IR become of interest for transcriptomic data analysis owing to its recognized roles in gene expression regulation and associations with complex diseases. In this article, we first review the function of IR in regulating gene expression in a number of biological processes, such as neuron differentiation and activation of CD4⁺ T cells. Next, we briefly review its association with diseases, such as Alzheimer's disease and cancers. Then, we describe state-of-the-art methods for IR detection, including RNA-seq analysis tools IRFinder and iREAD, highlighting their underlying principles and discussing their advantages and limitations. Finally, we discuss the challenges for IR detection and potential ways in which IR detection methods could be improved.

Somatic SF3B1 hotspot mutation in prolactinomas

The genetic basis and corresponding clinical relevance of prolactinomas remain poorly understood. Here, we perform whole genome sequencing (WGS) on 21 patients with prolactinomas to detect somatic mutations and then validate the mutations with digital polymerase chain reaction (PCR) analysis of tissue samples from 227 prolactinomas. We identify the same hotspot somatic mutation in splicing factor 3 subunit B1 (SF3B1R625H) in 19.8% of prolactinomas. These patients with mutant prolactinomas display higher prolactin (PRL) levels (p = 0.02) and shorter progression-free survival (PFS) (p = 0.02) compared to patients without the mutation. Moreover, we identify that the SF3B1R625H mutation causes aberrant splicing of estrogen related receptor gamma (ESRRG), which results in stronger binding of pituitary-specific positive transcription factor 1 (Pit-1), leading to excessive PRL secretion. Thus our study validates an important mutation and elucidates a potential mechanism underlying the pathogenesis of prolactinomas that may lead to the development of targeted therapeutics.

The Complex Relationship between HTLV-1 and Nonsense-Mediated mRNA Decay (NMD)

Before the establishment of an adaptive immune response, retroviruses can be targeted by several cellular host factors at different stages of the viral replication cycle. This intrinsic immunity relies on a large diversity of antiviral processes. In the case of HTLV-1 infection, these active innate host defense mechanisms are debated. Among these mechanisms, we focused on an RNA decay pathway called nonsense-mediated mRNA decay (NMD), which can target multiple viral RNAs, including HTLV-1 unspliced RNA, as has been recently demonstrated. NMD is a co-translational process that depends on the RNA helicase UPF1 and regulates the expression of multiple types of host mRNAs. RNA sensitivity to NMD depends on mRNA organization and the ribonucleoprotein (mRNP) composition. HTLV-1 has evolved several means to evade the NMD threat, leading to NMD inhibition. In the early steps of infection, NMD inhibition favours the production of HTLV-1 infectious particles, which may contribute to the survival of the fittest clones despite genome instability; however, its direct long-term impact remains to be investigated.

Stress-mediated convergence of splicing landscapes in male and female rock doves

BACKGROUND

The process of alternative splicing provides a unique mechanism by which eukaryotes are able to produce numerous protein products from the same gene. Heightened variability in the proteome has been thought to potentiate increased behavioral complexity and response flexibility to environmental stimuli, thus contributing to more refined traits on which natural and sexual selection can act. While it has been long known that various forms of environmental stress can negatively affect sexual behavior and reproduction, we know little of how stress can affect the alternative splicing associated with these events, and less still about how splicing may differ between sexes. Using the model of the rock dove (Columba livia), our team previously uncovered sexual dimorphism in the basal and stress-responsive gene transcription of a biological system necessary for facilitating sexual behavior and reproduction, the hypothalamic-pituitary-gonadal (HPG) axis. In this study, we delve further into understanding the mechanistic underpinnings of how changes in the environment can affect reproduction by testing the alternative splicing response of the HPG axis to an external stressor in both sexes.

RESULTS

This study reveals dramatic baseline differences in HPG alternative splicing between males and females. However, after subjecting subjects to a restraint stress paradigm, we found a significant reduction in these differences between the sexes. In both stress and control treatments, we identified a higher incidence of splicing activity in the pituitary in both sexes as compared to other tissues. Of these splicing events, the core exon event is the most abundant form of splicing and more frequently occurs in the coding regions of the gene. Overall, we observed less splicing activity in the 3'UTR (untranslated region) end of transcripts than the 5'UTR or coding regions.

CONCLUSIONS

Our results provide vital new insight into sex-specific aspects of the stress response on the HPG axis at an unprecedented proximate level. Males and females uniquely respond to stress, yet exhibit splicing patterns suggesting a convergent, optimal splicing landscape for stress response. This information has the potential to inform evolutionary theory as well as the development of highly-specific drug targets for stress-induced reproductive dysfunction.

Challenges in detecting and quantifying intron retention from next generation sequencing data

Intron retention (IR) occurs when an intron is transcribed into pre-mRNA and remains in the final mRNA. An increasing body of literature has demonstrated a major role for IR in numerous biological functions and in disease. Here we give an overview of the different computational approaches for detecting IR events from sequencing data. We show that these are based on different biological and computational assumptions that may lead to dramatically different results. We describe the various approaches for mitigating errors in detecting intron retention and for discovering IR signatures between different conditions.

iREAD: a tool for intron retention detection from RNA-seq data

BACKGROUND

Intron retention (IR) has been traditionally overlooked as 'noise' and received negligible attention in the field of gene expression analysis. In recent years, IR has become an emerging field for interrogating transcriptomes because it has been recognized to carry out important biological functions such as gene expression regulation and it has been found to be associated with complex diseases such as cancers. However, methods for detecting IR today are limited. Thus, there is a need to develop novel methods to improve IR detection.

RESULTS

Here we present iREAD (intron REtention Analysis and Detector), a tool to detect IR events genome-wide from high-throughput RNA-seq data. The command line interface for iREAD is implemented in Python. iREAD takes as input a BAM file, representing the transcriptome, and a text file containing the intron coordinates of a genome. It then 1) counts all reads that overlap intron regions, 2) detects IR events by analyzing the features of reads such as depth and distribution patterns, and 3) outputs a list of retained introns into a tab-delimited text file. iREAD provides significant added value in detecting IR compared with output from IRFinder with a higher AUC on all datasets tested. Both methods showed low false positive rates and high false negative rates in different regimes, indicating that use together is generally beneficial. The output from iREAD can be directly used for further exploratory analysis such as differential intron expression and functional enrichment. The software is freely available at https://github.com/genemine/iread.

CONCLUSION

Being complementary to existing tools, iREAD provides a new and generic tool to interrogate poly-A enriched transcriptomic data of intron regions. Intron retention analysis provides a complementary approach for understanding transcriptome.

Systematic analysis of survival-associated alternative splicing signatures in clear cell renal cell carcinoma

Alternative splicing (AS) constitutes a major reason for messenger RNA (mRNA) and protein diversity. Increasing studies have shown a link to splicing dysfunction associated with malignant neoplasia. Systematic analysis of AS events in kidney cancer remains poorly reported. Therefore, we generated AS profiles in 533 kidney renal clear cell carcinoma (KIRC) patients in The Cancer Genome Atlas (TCGA) database using RNA-seq data. Then, prognostic models were developed in a primary cohort (N = 351) and validated in a validation cohort (N = 182). In addition, splicing networks were built by integrating bioinformatics analyses. A total of 11 268 and 8083 AS variants were significantly associated with patient overall survival time in the primary and validation KIRC cohorts, respectively, including STAT1, DAZAP1, IDS, NUDT7, and KLHDC4. The AS events in the primary KIRC cohorts served as candidate AS events to screen the independent risk factors associated with survival in the primary cohort and to develop prognostic models. The area under the curve of the receiver-operator characteristic curve for prognostic prediction in the primary and validation KIRC cohorts was 0.84 and 0.82 at 2500 days of overall survival, respectively. In addition, splicing correlation networks revealed key splicing factors (SFs) in KIRC, such as HNRNPH1, HNRNPU, KHDBS1, KHDBS3, SRSF9, RBMX, SFQ, SRP54, HNRNPA0, and SRSF6. In this study, we analyzed the AS landscape in the TCGA KIRC cohort and detected predictors (prognostic) based on AS variants with high performance for risk stratification of the KIRC cohort and revealed key SFs in splicing networks, which could act as underlying mechanisms.

Food-Derived Compounds Apigenin and Luteolin Modulate mRNA Splicing of Introns with Weak Splice Sites

Cancer cells often exhibit extreme sensitivity to splicing inhibitors. We identified food-derived flavonoids, apigenin and luteolin, as compounds that modulate mRNA splicing at the genome-wide level, followed by proliferation inhibition. They bind to mRNA splicing-related proteins to induce a widespread change of splicing patterns in treated cells. Their inhibitory activity on splicing is relatively moderate, and introns with weak splice sites tend to be sensitive to them. Such introns remain unspliced, and the resulting intron-containing mRNAs are retained in the nucleus, resulting in the nuclear accumulation of poly(A)⁺ RNAs in these flavonoid-treated cells. Tumorigenic cells are more susceptible to these flavonoids than nontumorigenic cells, both for the nuclear poly(A)⁺ RNA-accumulating phenotype and cell viability. This study illustrates the possible mechanism of these flavonoids to suppress tumor progression in vivo that were demonstrated by previous studies and provides the potential of daily intake of moderate splicing inhibitors to prevent cancer development.

•

Food-derived compounds, apigenin and luteolin, modulate mRNA splicing

•

The treatment of these flavonoids causes numerous alternative splicing events

•

Splicing of introns with weak splice sites tend to be inhibited by these flavonoids

•

Tumorigenic cells are more sensitive to these flavonoids than non-tumorigenic cells

A missense mutation in SNRPE linked to non-syndromal microcephaly interferes with U snRNP assembly and pre-mRNA splicing

Malfunction of pre-mRNA processing factors are linked to several human diseases including cancer and neurodegeneration. Here we report the identification of a de novo heterozygous missense mutation in the SNRPE gene (c.65T>C (p.Phe22Ser)) in a patient with non-syndromal primary (congenital) microcephaly and intellectual disability. SNRPE encodes SmE, a basal component of pre-mRNA processing U snRNPs. We show that the microcephaly-linked SmE variant is unable to interact with the SMN complex and as a consequence fails to assemble into U snRNPs. This results in widespread mRNA splicing alterations in fibroblast cells derived from this patient. Similar alterations were observed in HEK293 cells upon SmE depletion that could be rescued by the expression of wild type but not mutant SmE. Importantly, the depletion of SmE in zebrafish causes aberrant mRNA splicing alterations and reduced brain size, reminiscent of the patient microcephaly phenotype. We identify the EMX2 mRNA, which encodes a protein required for proper brain development, as a major mis-spliced down stream target. Together, our study links defects in the SNRPE gene to microcephaly and suggests that alterations of cellular splicing of specific mRNAs such as EMX2 results in the neurological phenotype of the disease.

In higher eukaryotes, the protein coding genes are first transcribed as precursor mRNAs (pre-mRNAs) and further processed by the spliceosome to form the mature mRNA for translation. Malfunction of pre-mRNA processing factors are linked to several human diseases including cancer and neurodegeneration. Here we report the identification of a de novo heterozygous missense mutation in the SNRPE/SmE gene in a patient with non-syndromal primary (congenital) microcephaly and intellectual disability. The effect of identified de novo mutation on SNRPE/SmE was characterized in vitro. The zebrafish was used as in vivo model to further dissect the physiological consequence and pathomechanism. Finally, the EMX2 gene was identified as one of the major down stream target genes responsible for the phenotype. Our study links defects in the SNRPE/SmE gene to microcephaly and provides the new pathogenic mechanism for microcephaly.

Introns: Good Day Junk Is Bad Day Treasure

Introns are ubiquitous in eukaryotic transcripts. They are often viewed as junk RNA but the huge energetic burden of transcribing, removing, and degrading them suggests a significant evolutionary advantage. Ostensibly, an intron functions within the host pre-mRNA to regulate its splicing, transport, and degradation. However, recent studies have revealed an entirely new class of trans-acting functions where the presence of intronic RNA in the cell impacts the expression of other genes in trans. Here, we review possible new mechanisms of intron functions, with a focus on the role of yeast introns in regulating the cell growth response to starvation.

Identification of prognosis-related alternative splicing events in kidney renal clear cell carcinoma

Alternative splicing (AS) contributes to protein diversity by modifying most gene transcriptions. Cancer generation and progression are associated with specific splicing events. However, AS signature in kidney renal clear cell carcinoma (KIRC) remains unknown. In this study, genome‐wide AS profiles were generated in 537 patients with KIRC in the cancer genome atlas. With a total of 42 522 mRNA AS events in 10 600 genes acquired, 8164 AS events were significantly associated with the survival of patients with KIRC. Logistic regression analysis of the least absolute shrinkage and selection operator was conducted to identify an optimized multivariate prognostic predicting mode containing four predictors. In this model, the receptor‐operator characteristic curves of the training set were built, and the areas under the curves (AUCs) at different times were >0.88, thus indicating a stable and powerful ability in distinguishing patients' outcome. Similarly, the AUCs of the test set at different times were >0.73, verifying the results of the training set. Correlation and gene ontology analyses revealed some potential functions of prognostic AS events. This study provided an optimized survival‐predicting model and promising data resources for future in‐depth studies on AS mechanisms in KIRC.

Pan-cancer analysis of clinical relevance of alternative splicing events in 31 human cancers

Alternative splicing represents a critical posttranscriptional regulation of gene expression, which contributes to the protein complexity and mRNA processing. Defects of alternative splicing including genetic alteration and/or altered expression of both pre-mRNA and trans-acting factors give rise to many cancers. By integrally analyzing clinical data and splicing data from TCGA and SpliceSeq databases, a number of splicing events were found clinically relevant in tumor samples. Alternative splicing of KLK2 (KLK2_51239) was found as a potential inducement of nonsense-mediated mRNA decay and associated with poor survival in prostate cancer. Consensus K-means clustering analysis indicated that alternative splicing events could be potentially used for molecular subtype classification of cancers. By random forest survival algorithm, prognostic prediction signatures with well performances were constructed for 31 cancers by using survival-associated alternative splicing events. Furthermore, an online tool for visualization of Kaplan-Meier plots of splicing events in 31 cancers was explored. Briefly, alternative splicing was found of significant clinical relevance with cancers.

Distinct types of short open reading frames are translated in plant cells

Genomes contain millions of short (<100 codons) open reading frames (sORFs), which are usually dismissed during gene annotation. Nevertheless, peptides encoded by such sORFs can play important biological roles, and their impact on cellular processes has long been underestimated. Here, we analyzed approximately 70,000 transcribed sORFs in the model plant Physcomitrella patens (moss). Several distinct classes of sORFs that differ in terms of their position on transcripts and the level of evolutionary conservation are present in the moss genome. Over 5000 sORFs were conserved in at least one of 10 plant species examined. Mass spectrometry analysis of proteomic and peptidomic data sets suggested that tens of sORFs located on distinct parts of mRNAs and long noncoding RNAs (lncRNAs) are translated, including conserved sORFs. Translational analysis of the sORFs and main ORFs at a single locus suggested the existence of genes that code for multiple proteins and peptides with tissue-specific expression. Functional analysis of four lncRNA-encoded peptides showed that sORFs-encoded peptides are involved in regulation of growth and differentiation in moss. Knocking out lncRNA-encoded peptides resulted in a decrease of moss growth. In contrast, the overexpression of these peptides resulted in a diverse range of phenotypic effects. Our results thus open new avenues for discovering novel, biologically active peptides in the plant kingdom.

Alternative splicing events are prognostic in hepatocellular carcinoma

Alternative splicing events (ASEs) play a role in cancer development and progression. We investigated whether ASEs are prognostic for overall survival (OS) in hepatocellular carcinoma (HCC). RNA sequencing data was obtained for 343 patients included in The Cancer Genome Atlas. Matched splicing event data for these patients was then obtained from the TCGASpliceSeq database, which includes data for seven types of ASEs. Univariate and multivariate Cox regression analysis demonstrated that 3,814 OS-associated splicing events (OS-SEs) were correlated with OS. Prognostic indices were developed based on the most significant OS-SEs. The prognostic index based on all seven types of ASEs (PI-ALL) demonstrated superior efficacy in predicting OS of HCC patients at 2,000 days compared to those based on single ASE types. Patients were stratified into two risk groups (high and low) based on the median prognostic index. Kaplan-Meier survival analysis demonstrated that PI-ALL had the greatest capacity to distinguish between patients with favorable vs. poor outcomes. Finally, univariate Cox regression analysis demonstrated that the expression of 23 splicing factors was correlated with OS-SEs in the HCC cohort. Our data indicate that a prognostic index based on ASEs is prognostic for OS in HCC.