Computational Identification of Tissue-Specific Splicing Regulatory Elements in Human Genes from RNA-Seq Data

mRNA Transcript Variants Expressed in Mammalian Cells

Gene expression or gene regulation studies often assume one gene expresses one mRNA. However, contrary to the conventional idea, a single gene in mammalian cells can express multiple transcript variants translated into several different proteins. The transcript variants are generated through transcription from alternative start sites and alternative post-transcriptional processing of the precursor mRNA (pre-mRNA). In addition, gene mutations and RNA editing further enhance the diversity of the transcript variants. The transcript variants can encode proteins with various domains, expanding the functional repertoire of a single gene. Some transcript variants may not encode proteins but function as non-coding RNAs and regulate gene expression. The expression level of the transcript variants may vary between cell types or within the same cells under different biological conditions. Transcript variants are characteristic of cell differentiation in a particular tissue, and the variants may play a key role in normal development and aging. Studies also reported that some transcript variants may have roles in disease pathogenesis. The biological significances urge studying the complexity of gene expression at the transcript level. This article updates the molecular basis of transcript variants in mammalian cells, including the formation mechanisms and potential roles in host biology. Gaining insight into the transcript variants will not only identify novel mechanisms of gene regulation but also unravel the role of the variants in health and disease.

Alternative splicing of modulatory immune receptors in T lymphocytes: a newly identified and targetable mechanism for anticancer immunotherapy

Alternative splicing (AS) is a mechanism that generates translational diversity within a genome. Equally important is the dynamic adaptability of the splicing machinery, which can give preference to one isoform over others encoded by a single gene. These isoform preferences change in response to the cell's state and function. Particularly significant is the impact of physiological alternative splicing in T lymphocytes, where specific isoforms can enhance or reduce the cells' reactivity to stimuli. This process makes splicing isoforms defining features of cell states, exemplified by CD45 splice isoforms, which characterize the transition from naïve to memory states. Two developments have accelerated the use of AS dynamics for therapeutic interventions: advancements in long-read RNA sequencing and progress in nucleic acid chemical modifications. Improved oligonucleotide stability has enabled their use in directing splicing to specific sites or modifying sequences to enhance or silence particular splicing events. This review highlights immune regulatory splicing patterns with potential significance for enhancing anticancer immunotherapy.

Transcriptomic Alterations in Spliceosome Components in Advanced Heart Failure: Status of Cardiac-Specific Alternative Splicing Factors

Heart failure (HF) is associated with global changes in gene expression. Alternative mRNA splicing (AS) is a key regulatory mechanism underlying these changes. However, the whole status of molecules involved in the splicing process in human HF is unknown. Therefore, we analysed the spliceosome transcriptome in cardiac tissue (n = 36) from control subjects and HF patients (with ischaemic (ICM) and dilated (DCM) cardiomyopathies) using RNA-seq. We found greater deregulation of spliceosome machinery in ICM. Specifically, we showed widespread upregulation of the E and C complex components, highlighting an increase in SNRPD2 (FC = 1.35, p < 0.05) and DHX35 (FC = 1.34, p < 0.001) mRNA levels. In contrast, we observed generalised downregulation of the A complex and cardiac-specific AS factors, such as the multifunctional protein PCBP2 (FC = -1.29, p < 0.001) and the RNA binding proteins QKI (FC = -1.35, p < 0.01). In addition, we found a relationship between SNPRD2 (an E complex component) and the left ventricular mass index in ICM patients (r = 0.779; p < 0.01). On the other hand, we observed the specific underexpression of DDX46 (FC = -1.29), RBM17 (FC = -1.33), SDE2 (FC = -1.35) and RBFOX1 (FC = -1.33), p < 0.05, in DCM patients. Therefore, these aetiology-related alterations may indicate the differential involvement of the splicing process in the development of ICM and DCM.

The identification of protein and RNA interactors of the splicing factor Caper in the adult Drosophila nervous system

Post-transcriptional gene regulation is a fundamental mechanism that helps regulate the development and healthy aging of the nervous system. Mutations that disrupt the function of RNA-binding proteins (RBPs), which regulate post-transcriptional gene regulation, have increasingly been implicated in neurological disorders including amyotrophic lateral sclerosis, Fragile X Syndrome, and spinal muscular atrophy. Interestingly, although the majority of RBPs are expressed widely within diverse tissue types, the nervous system is often particularly sensitive to their dysfunction. It is therefore critical to elucidate how aberrant RNA regulation that results from the dysfunction of ubiquitously expressed RBPs leads to tissue specific pathologies that underlie neurological diseases. The highly conserved RBP and alternative splicing factor Caper is widely expressed throughout development and is required for the development of Drosophila sensory and motor neurons. Furthermore, caper dysfunction results in larval and adult locomotor deficits. Nonetheless, little is known about which proteins interact with Caper, and which RNAs are regulated by Caper. Here we identify proteins that interact with Caper in both neural and muscle tissue, along with neural specific Caper target RNAs. Furthermore, we show that a subset of these Caper-interacting proteins and RNAs genetically interact with caper to regulate Drosophila gravitaxis behavior.

Polymorphic differences within HLA-C alleles contribute to alternatively spliced transcripts lacking exon 5

The HLA genes are amongst the most polymorphic in the human genome. Alternative splicing could add an extra layer of complexity, but has not been studied extensively. Here, we applied an RNA based approach to study the influence of allele polymorphism on alternative splicing of HLA‐C in peripheral blood. RNA was isolated from these peripheral cells, converted into cDNA and amplified specifically for 12 common HLA‐C allele groups. Through subsequent sequencing of HLA‐C, we observed alternative splicing variants of HLA‐C*04 and *16 that resulted in exon 5 skipping and were co‐expressed with the mature transcript. Investigation of intron 4 sequences of HLA‐C*04 and *16 compared with other HLA‐C alleles demonstrated no effect on predicted splice sites and branch point. To further investigate if the unique polymorphic positions in exon 5 of HLA‐C*04 or *16 may facilitate alternative splicing by acting on splicing regulatory elements (SRE), in‐silico splicing analysis was performed. While the HLA‐C*04 specific SNP in exon 5 had no effect on predicted exonic SRE, the HLA‐C*16 specific exon 5 SNP did alter exonic SRE. Our findings provide experimental and theoretical support for the concept that polymorphisms within the HLA‐C alleles influence the alternative splicing of HLA‐C.

Transcript isoforms of Reep6 have distinct functions in the retina

Much of the complexity of the eukaryotic cell transcriptome is due to the alternative splicing of mRNA. However, knowledge on how transcriptome complexity is translated into functional complexity remains limited. For example, although different isoforms of a gene may show distinct temporal and spatial expression patterns, it is largely unknown whether these isoforms encode proteins with distinct functions matching their expression pattern. In this report, we investigated the function and relationship of the two isoforms of Reep6, namely Reep6.1 and Reep6.2, in rod photoreceptor cells. These two isoforms result from the alternative splicing of exon 5 and show mutually exclusive expression patterns. Reep6.2 is the canonical isoform that is expressed in non-retinal tissues while Reep6.1 is the only expressed isoform in the adult retina. The Reep6.1 isoform-specific knockout mouse, Reep6E5/E5, is generated by deleting exon 5 and a homozygous deletion phenotypically displayed a rod degeneration phenotype comparable to a Reep6 full knockout mouse, indicating that the Reep6.1 isoform is essential for the rod photoreceptor cell survival. Consistent with the results obtained from a loss-of-function experiment, overexpression of Reep6.2 failed to rescue the rod degeneration phenotype of Reep6 knockout mice while overexpression of Reep6.1 does lead to rescue. These results demonstrate that, consistent with the expression pattern of the isoform, Reep6.1 has rod-specific functions that cannot be substituted by its canonical isoform. Our findings suggested that a strict regulation of splicing is required for the maintenance of photoreceptor cells.

MTSplice predicts effects of genetic variants on tissue-specific splicing

We develop the free and open-source model Multi-tissue Splicing (MTSplice) to predict the effects of genetic variants on splicing of cassette exons in 56 human tissues. MTSplice combines MMSplice, which models constitutive regulatory sequences, with a new neural network that models tissue-specific regulatory sequences. MTSplice outperforms MMSplice on predicting tissue-specific variations associated with genetic variants in most tissues of the GTEx dataset, with largest improvements on brain tissues. Furthermore, MTSplice predicts that autism-associated de novo mutations are enriched for variants affecting splicing specifically in the brain. We foresee that MTSplice will aid interpreting variants associated with tissue-specific disorders.

Creation and validation of models to predict response to primary treatment in serous ovarian cancer

Nearly a third of patients with high-grade serous ovarian cancer (HGSC) do not respond to initial therapy and have an overall poor prognosis. However, there are no validated tools that accurately predict which patients will not respond. Our objective is to create and validate accurate models of prediction for treatment response in HGSC. This is a retrospective case–control study that integrates comprehensive clinical and genomic data from 88 patients with HGSC from a single institution. Responders were those patients with a progression-free survival of at least 6 months after treatment. Only patients with complete clinical information and frozen specimen at surgery were included. Gene, miRNA, exon, and long non-coding RNA (lncRNA) expression, gene copy number, genomic variation, and fusion-gene determination were extracted from RNA-sequencing data. DNA methylation analysis was performed. Initial selection of informative variables was performed with univariate ANOVA with cross-validation. Significant variables (p < 0.05) were included in multivariate lasso regression prediction models. Initial models included only one variable. Variables were then combined to create complex models. Model performance was measured with area under the curve (AUC). Validation of all models was performed using TCGA HGSC database. By integrating clinical and genomic variables, we achieved prediction performances of over 95% in AUC. Most performances in the validation set did not differ from the training set. Models with DNA methylation or lncRNA underperformed in the validation set. Integrating comprehensive clinical and genomic data from patients with HGSC results in accurate and robust prediction models of treatment response.

Novel insights into biological roles of inducible cAMP early repressor ICER

ICER corresponds to a group of alternatively spliced Inducible cAMP Early Repressors with high similarity, but multiple roles, including in circadian rhythm, and are involved in attenuation of cAMP-dependent gene expression. We present experimental and in silico data revealing biological differences between the isoforms with exon gamma (ICER) or without it (ICERγ). Both isoforms are expressed in the liver and the adrenal glands and can derive from differential splicing. In adrenals the expression is circadian, with maximum at ZT12 and higher amplitude of Icerγ. In the liver, the expression of Icerγ is lower than Icer in the 24 h time frame. Icer mRNA has a delayed early response to forskolin. The longer ICER protein binds to three DNA grooves of the Per1 promoter, while ICERγ only to two, as deduced by molecular modelling. This is in line with gel shift competition assays showing stronger binding of ICER to Per1 promotor. Only Icerγ siRNA provoked an increase of Per1 expression. In conclusion, we show that ICER and ICERγ have distinct biochemical properties in tissue expression, DNA binding, and response to forskolin. Data are in favour of ICERγ as the physiologically important form in hepatic cells where weaker binding of repressor might be preferred in guiding the cAMP-dependent response.

Microcephaly Modeling of Kinetochore Mutation Reveals a Brain-Specific Phenotype

Most genes mutated in microcephaly patients are expressed ubiquitously, and yet the brain is the only major organ compromised in most patients. Why the phenotype remains brain specific is poorly understood. In this study, we used in vitro differentiation of human embryonic stem cells to monitor the effect of a point mutation in kinetochore null protein 1 (KNL1;CASC5), identified in microcephaly patients, during in vitro brain development. We found that neural progenitors bearing a patient mutation showed reduced KNL1 levels, aneuploidy, and an abrogated spindle assembly checkpoint. By contrast, no reduction of KNL1 levels or abnormalities was observed in fibroblasts and neural crest cells. We established that the KNL1 patient mutation generates an exonic splicing silencer site, which mainly affects neural progenitors because of their higher levels of splicing proteins. Our results provide insight into the brain-specific phenomenon, consistent with microcephaly being the only major phenotype of patients bearing KNL1 mutation.

Using 3D neural spheroids, Javed et al. investigate a mutation in KNL1 that causes microcephaly. Their study shows that, despite ubiquitous mutant KNL1 expression, KNL1 mRNA processing is affected only in neural precursors due to difference in splicing protein levels, offering insights into why the phenotype remains brain specific in patients.

Sequence and Evolutionary Features for the Alternatively Spliced Exons of Eukaryotic Genes

Alternative splicing of pre-mRNAs is a crucial mechanism for maintaining protein diversity in eukaryotes without requiring a considerable increase of genes in the number. Due to rapid advances in high-throughput sequencing technologies and computational algorithms, it is anticipated that alternative splicing events will be more intensively studied to address different kinds of biological questions. The occurrences of alternative splicing mean that all exons could be classified to be either constitutively or alternatively spliced depending on whether they are virtually included into all mature mRNAs. From an evolutionary point of view, therefore, the alternatively spliced exons would have been associated with distinctive biological characteristics in comparison with constitutively spliced exons. In this paper, we first outline the representative types of alternative splicing events and exon classification, and then review sequence and evolutionary features for the alternatively spliced exons. The main purpose is to facilitate understanding of the biological implications of alternative splicing in eukaryotes. This knowledge is also helpful to establish computational approaches for predicting the splicing pattern of exons.

The LL-100 panel: 100 cell lines for blood cancer studies

For many years, immortalized cell lines have been used as model systems for cancer research. Cell line panels were established for basic research and drug development, but did not cover the full spectrum of leukemia and lymphoma. Therefore, we now developed a novel panel (LL-100), 100 cell lines covering 22 entities of human leukemia and lymphoma including T-cell, B-cell and myeloid malignancies. Importantly, all cell lines are unequivocally authenticated and assigned to the correct tissue. Cell line samples were proven to be free of mycoplasma and non-inherent virus contamination. Whole exome sequencing and RNA-sequencing of the 100 cell lines were conducted with a uniform methodology to complement existing data on these publicly available cell lines. We show that such comprehensive sequencing data can be used to find lymphoma-subtype-characteristic copy number aberrations, mRNA isoforms, transcription factor activities and expression patterns of NKL homeobox genes. These exemplary studies confirm that the novel LL-100 panel will be useful for understanding the function of oncogenes and tumor suppressor genes and to develop targeted therapies.

Combinatorial Pattern of Histone Modifications in Exon Skipping Event

Histone modifications are associated with alternative splicing. It has been suggested that histone modifications act in combinational patterns in gene expression regulation. However, how they interact with each other and what is their casual relationships in the process of RNA splicing remain unclear. In this study, the combinatorial patterns of 38 kinds of histone modifications in the exon skipping event of the CD4⁺ T cell were analyzed by constructing Bayesian networks. Distinct combinatorial patterns of histone modifications that illustrating their casual relationships were observed in excluded/included exons and the surrounding intronic regions. The Bayesian networks also indicate that some histone modifications directly correlate with RNA splicing. We anticipate that this work could provide novel insights into the effects of histone modifications on RNA splicing regulation.

Pathogenicity and selective constraint on variation near splice sites

Mutations that perturb normal pre-mRNA splicing are significant contributors to human disease. We used exome sequencing data from 7833 probands with developmental disorders (DDs) and their unaffected parents, as well as more than 60,000 aggregated exomes from the Exome Aggregation Consortium, to investigate selection around the splice sites and quantify the contribution of splicing mutations to DDs. Patterns of purifying selection, a deficit of variants in highly constrained genes in healthy subjects, and excess de novo mutations in patients highlighted particular positions within and around the consensus splice site of greater functional relevance. By using mutational burden analyses in this large cohort of proband–parent trios, we could estimate in an unbiased manner the relative contributions of mutations at canonical dinucleotides (73%) and flanking noncanonical positions (27%), and calculate the positive predictive value of pathogenicity for different classes of mutations. We identified 18 patients with likely diagnostic de novo mutations in dominant DD-associated genes at noncanonical positions in splice sites. We estimate 35%–40% of pathogenic variants in noncanonical splice site positions are missing from public databases.

Advances in Analyzing Virus-Induced Alterations of Host Cell Splicing

Alteration of host cell splicing is a common feature of many viral infections which is underappreciated because of the complexity and technical difficulty of studying alternative splicing (AS) regulation. Recent advances in RNA sequencing technologies revealed that up to several hundreds of host genes can show altered mRNA splicing upon viral infection. The observed changes in AS events can be either a direct consequence of viral manipulation of the host splicing machinery or result indirectly from the virus-induced innate immune response or cellular damage. Analysis at a higher resolution with single-cell RNAseq, and at a higher scale with the integration of multiple omics data sets in a systems biology perspective, will be needed to further comprehend this complex facet of virus-host interactions.

Transcriptional Profiling of Host Cell Responses to Virulent Haemophilus parasuis: New Insights into Pathogenesis

Haemophilus parasuis is the causative agent of Glässer’s disease in pigs. H. parasuis can cause vascular damage, although the mechanism remains unclear. In this study, we investigated the host cell responses involved in the molecular pathway interactions in porcine aortic vascular endothelial cells (PAVECs) induced by H. parasuis using RNA-Seq. The transcriptome results showed that when PAVECs were infected with H. parasuis for 24 h, 281 differentially expressed genes (DEGs) were identified; of which, 236 were upregulated and 45 downregulated. The 281 DEGs were involved in 136 KEGG signaling pathways that were organismal systems, environmental information processing, metabolism, cellular processes, and genetic information processing. The main pathways were the Rap1, FoxO, and PI3K/Akt signaling pathways, and the overexpressed genes were determined and verified by quantitative reverse transcription polymerase chain reaction. In addition, 252 genes were clustered into biological processes, molecular processes, and cellular components. Our study provides new insights for understanding the interaction between bacterial and host cells, and analyzed, in detail, the possible mechanisms that lead to vascular damage induced by H. parasuis. This may lead to development of novel therapeutic targets to control H. parasuis infection.

Computer Analysis of Glioma Transcriptome Profiling: Alternative Splicing Events

Here we present the analysis of alternative splicing events on an example of glioblastoma cell culture samples using a set of computer tools in combination with database integration. The gene expression profiles of glioblastoma were obtained from cell culture samples of primary glioblastoma which were isolated and processed for RNA extraction. Transcriptome profiling of normal brain samples and glioblastoma were done by Illumina sequencing. The significant differentially expressed exon-level probes and their corresponding genes were identified using a combination of the splicing index method. Previous studies indicated that tumor-specific alternative splicing is important in the regulation of gene expression and corresponding protein functions during cancer development. Multiple alternative splicing transcripts have been identified as progression markers, including generalized splicing abnormalities and tumor- and stage-specific events. We used a set of computer tools which were recently applied to analysis of gene expression in laboratory animals to study differential splicing events. We found 69 transcripts that are differentially alternatively spliced. Three cancer-associated genes were considered in detail, in particular: APP (amyloid beta precursor protein), CASC4 (cancer susceptibility candidate 4) and TP53. Such alternative splicing opens new perspectives for cancer research.

The determinants of alternative RNA splicing in human cells

Alternative splicing represents an important level of the regulation of gene function in eukaryotic organisms. It plays a critical role in virtually every biological process within an organism, including regulation of cell division and cell death, differentiation of tissues in the embryo and the adult organism, as well as in cellular response to diverse environmental factors. In turn, studies of the last decade have shown that alternative splicing itself is controlled by different mechanisms. Unfortunately, there is no clear understanding of how these diverse mechanisms, or determinants, regulate and constrain the set of alternative RNA species produced from any particular gene in every cell of the human body. Here, we provide a consolidated overview of alternative splicing determinants including RNA-protein interactions, epigenetic regulation via chromatin remodeling, coupling of transcription-to-alternative splicing, effect of secondary structures in pre-RNA, and function of the RNA quality control systems. We also extensively and critically discuss some mechanistic insights on coordinated inclusion/exclusion of exons during the formation of mature RNA molecules. We conclude that the final structure of RNA is pre-determined by a complex interplay between cis- and trans-acting factors. Altogether, currently available empirical data significantly expand our understanding of the functioning of the alternative splicing machinery of cells in normal and pathological conditions. On the other hand, there are still many blind spots that require further deep investigations.