Alternative splicing as a regulator of development and tissue identity

Thermosensitive alternative splicing senses and mediates temperature adaptation in Drosophila

Circadian rhythms are generated by the cyclic transcription, translation, and degradation of clock gene products, including timeless (tim), but how the circadian clock senses and adapts to temperature changes is not completely understood. Here, we show that temperature dramatically changes the splicing pattern of tim in Drosophila. We found that at 18°C, TIM levels are low because of the induction of two cold-specific isoforms: tim-cold and tim-short and cold. At 29°C, another isoform, tim-medium, is upregulated. Isoform switching regulates the levels and activity of TIM as each isoform has a specific function. We found that tim-short and cold encodes a protein that rescues the behavioral defects of tim⁰¹ mutants, and that flies in which tim-short and cold is abrogated have abnormal locomotor activity. In addition, miRNA-mediated control limits the expression of some of these isoforms. Finally, data that we obtained using minigenes suggest that tim alternative splicing might act as a thermometer for the circadian clock.

Cassava AGPase genes and their encoded proteins are different from those of other plants

Cassava AGPase and AGPase genes have some unique characteristics. ADP-glucose pyrophosphorylase (AGPase) is a rate-limiting enzyme for starch synthesis. In this study, cassava AGPase genes (MeAGP) were analyzed based on six cultivars and one wild species. A total of seven MeAGPs was identified, including four encoding AGPase large subunits (MeAGPLs 1, 2, 3 and 4) and three encoding AGPase small subunits (MeAGPSs 1, 2 and 3). The copy number of MeAGPs varied in cassava germplasm materials. There were 14 introns for MeAGPLs 1, 2 and 3, 13 introns for MeAGPL4, and 8 introns for other three MeAGPSs. Multiple conservative amino acid sequence motifs were found in the MeAGPs. There were differences in amino acids at binding sites of substrates and regulators among different MeAGP subunits and between MeAGPs and a potato AGPase small subunit (1YP2:B). MeAGPs were all located in chloroplasts. MeAGP expression was not only associated with gene copy number and types/combinations, regions and levels of the DNA methylation but was also affected by environmental factors with the involvement of various transcription factors in multiple regulation networks and in various cis-elements in the gene promoter regions. The MeAGP activity also changed with environmental conditions and had potential differences among the subunits. Taken together, MeAGPs differ in number from those of Arabidopsis, potato, maize, banana, sweet potato, and tomato.

Intronic RNA: Ad'junk' mediator of post-transcriptional gene regulation

RNA splicing, the process through which intervening segments of noncoding RNA (introns) are excised from pre-mRNAs to allow for the formation of a mature mRNA product, has long been appreciated for its capacity to add complexity to eukaryotic proteomes. However, evidence suggests that the utility of this process extends beyond protein output and provides cells with a dynamic tool for gene regulation. In this review, we aim to highlight the role that intronic RNA plays in mediating specific splicing outcomes in pre-mRNA processing, as well as explore an emerging class of stable intronic sequences that have been observed to act in gene expression control. Building from underlying flexibility in both sequence and structure, intronic RNA provides mechanisms for post-transcriptional gene regulation that are amenable to the tissue and condition specific needs of eukaryotic cells. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

Identification of a novel RIG-I isoform and its truncating variant in Japanese eel, Anguilla japonica

Retinoic acid-inducible gene-I (RIG-I) is a cytoplasmic viral RNA sensor that triggers the production of type I interferons (IFNs) and proinflammatory cytokines during viral infection. RIG-I gene has been identified previously in Japanese eel, Anguilla japonica. In the present study, we have characterized a novel isoform of RIG-I (designated as AjRIG-Ib) and its truncated variant (AjRIG-Ibv). The AjRIG-Ib encodes 940 amino acids (aa) consisting of two N-terminal caspase activation and recruitment domains (CARDs), a DEX(D/H) box RNA helicase domain, and a C-terminal regulatory domain (CTD). The AjRIG-Ibv encodes a protein of 843 aa, that shares similar structural organization with AjRIG-Ib, but lacking CTD. The gene expression analyses showed that AjRIG-Ib and AjRIG-Ibv were detectable in all tissues/organs examined, and AjRIG-Ib was the predominant form. The mRNA level of AjRIG-Ibv was upregulated rapidly at 8 h after the Poly I:C injection, and the significant increase of AjRIG-Ib was observed at 16 and 24 h post-injection (hpi). Laser confocal microscopy showed that AjRIG-Ib and AjRIG-Ibv were both located in cytoplasm. In addition, the overexpression of AjRIG-Ib or AjRIG-Ibv led to the increased activity of IFN promoter in transient transfection assay. Taken together, our results indicated that AjRIG-Ib and AjRIG-Ibv may play cooperative or somewhat complementary roles in coordinating the antiviral response in fish.

TCF3 mutually exclusive alternative splicing is controlled by long-range cooperative actions between hnRNPH1 and PTBP1

TCF3, also known as E2A, is a well-studied transcription factor that plays an important role in stem cell maintenance and hematopoietic development. The TCF3 gene encodes two related proteins, E12 and E47, which arise from mutually exclusive alternative splicing (MEAS). Since these two proteins have different DNA binding and dimerization domains, this AS event must be strictly regulated to ensure proper isoform ratios. Previously, we found that heterogeneous nuclear ribonucleoprotein (hnRNP) H1/F regulates TCF3 AS by binding to exonic splicing silencers (ESSs) in exon 18b. Here, we identify conserved intronic splicing silencers (ISSs) located between, and far from, the two mutually exclusive exons, and show that they are essential for MEAS. Further, we demonstrate that the hnRNP PTBP1 binds the ISS and is a regulator of TCF3 AS. We also demonstrate that hnRNP H1 and PTBP1 regulate TCF3 AS reciprocally, and that position-dependent interactions between these factors are essential for proper TCF3 MEAS. Our study provides a new model in which MEAS is regulated by cooperative actions of distinct hnRNPs bound to ISSs and ESSs.

The UGTome: The expanding diversity of UDP glycosyltransferases and its impact on small molecule metabolism

The UDP glycosyltransferase (UGT) superfamily of enzymes is responsible for the metabolism and clearance of thousands of lipophilic chemicals including drugs, toxins and endogenous signaling molecules. They provide a protective interface between the organism and its chemical-rich environment, as well as controlling critical signaling pathways to maintain healthy tissue function. UGTs are associated with drug responses and interactions, as well as a wide range of diseases including cancer. The human genome contains 22 UGT genes; however as befitting their exceptionally diverse substrate ranges and biological activities, the output of these UGT genes is functionally diversified by multiple processes including alternative splicing, post-translational modification, homo- and hetero-oligomerization, and interactions with other proteins. All UGT genes are subject to extensive alternative splicing generating variant/truncated UGT proteins with altered functions including the capacity to dominantly modulate/inhibit cognate full-length forms. Heterotypic oligomerization of different UGTs can alter kinetic properties relative to monotypic complexes, and potentially produce novel substrate specificities. Moreover, the recently profiled interactions of UGTs with non-UGT proteins may facilitate coordination between different metabolic processes, as well as providing opportunities for UGTs to engage in novel 'moonlighting' functions. Herein we provide a detailed and comprehensive review of all known modes of UGT functional diversification and propose a UGTome model to describe the resulting expansion of metabolic capacity and its potential to modulate drug/xenobiotic responses and cell behaviours in normal and disease contexts.

Multi-strategic RNA-seq analysis reveals a high-resolution transcriptional landscape in cotton

Cotton is an important natural fiber crop, however, its comprehensive and high-resolution gene map is lacking. Here we integrate four complementary high-throughput techniques, including Pacbio long read Iso-seq, strand-specific RNA-seq, CAGE-seq, and PolyA-seq, to systematically explore the transcription landscape across 16 tissues or different organ types in Gossypium arboreum. We devise a computational pipeline, named IGIA, to reconstruct accurate gene structures from the integrated data. Our results reveal a dynamic and diverse transcriptional map in cotton: tissue-specific gene expression, alternative usage of TSSs and polyadenylation sites, hotspot of alternative splicing, and transcriptional read-through. These regulated events affect many genes in various aspects such as gain or loss of functional RNA motifs and protein domains, fine-tuning of DNA binding activity, and co-regulation for genes in the same complex or pathway. The methods and findings provide valuable resources for further functional genomic studies such as understanding natural SNP variations for plant community.

Genome-wide profiling of prognosis-related alternative splicing signatures in sarcoma

BACKGROUND

Sarcomas (SARCs) are rare malignant tumors with poor prognosis. Increasing evidence has suggested that aberrant alternative splicing (AS) is strongly associated with tumor initiation and progression. We considered whether survival-related AS events might serve as prognosis predictors and underlying targeted molecules in SARC treatment.

METHODS

RNA-Seq data of the SARC cohort were downloaded from The Cancer Genome Atlas (TCGA) database. Survival-related AS events were selected by univariate and multivariate Cox regression analyses. Metascape was used for constructing a gene interaction network and performing functional enrichment analysis. Then, prognosis predictors were established based on statistically significant survival-related AS events and evaluated by receiver operator characteristic (ROC) curve analysis. Finally, the potential regulatory network was analyzed via Pearson's correlation between survival-related AS events and splicing factors (SFs).

RESULTS

A total of 3,610 AS events and 2,291 genes were found to be prognosis-related in 261 SARC samples. The focal adhesion pathway was identified as the most critical molecular mechanism corresponding to poor prognosis. Notably, several prognosis predictors based on survival-related AS events showed excellent performance in prognosis prediction. The area under the curve of the ROC of the risk score was 0.85 in the integrated predictor. The splicing network proved complicated regulation between prognosis-related SFs and AS events. Also, driver gene mutations were significantly associated with AS in SARC patients.

CONCLUSIONS

Survival-related AS events may become ideal indictors for the prognosis prediction of SARCs. Corresponding splicing regulatory mechanisms are worth further exploration.

Systematic profiling of a novel prognostic alternative splicing signature in hepatocellular carcinoma

Alternative splicing (AS) is a pervasive and vital mechanism involved in the progression of cancer by expanding genomic encoding capacity and increasing protein complexity. However, the systematic analysis of AS in hepatocellular carcinoma (HCC) is lacking and urgently required. In the present study, genome‑wide AS events with corresponding clinical information were profiled in 290 patients with HCC from the Cancer Genome Atlas and SpliceSeq software. Functional enrichment analyses revealed the pivotal biological process of AS regulation. Univariate Cox regression analyses were performed, followed by stepwise forward multivariate analysis to develop the prognostic signatures. Spearman's correlation analyses were also used to construct potential regulatory network between the AS events and aberrant splicing factors. A total of 34,163 AS events were detected, among which 1,805 AS events from 1,314 parent genes were significantly associated with the overall survival (OS) of patients with HCC, and their parent genes serve crucial roles in HCC‑related oncogenic processes, including the p53 signaling pathway, AMPK signaling pathway and HIF‑1 signaling pathway. A prognostic AS signature was established that was found to be an independent prognostic factor for OS in stratified cohorts, harboring a noteworthy ability to distinguish between the distinct prognoses of patients with HCC (high‑risk vs. low‑risk, 827 vs. 3,125 days, P<2e‑16). Time‑dependent receiver‑-operator characteristic curves confirmed its robustness and clinical efficacy, with the area under the curves maintained >0.9 for short‑term and long‑term prognosis prediction. The splicing correlation network suggested a trend in the interactions between splicing factors and prognostic AS events, further revealing the underlying mechanism of AS in the oncogenesis of HCC. In conclusion, the present study provides a comprehensive portrait of global splicing alterations involved in the progression and HCC in addition to valuable prognostic factors for patients, which may represent as underappreciated hallmark and provide novel clues of therapeutic targets in HCC.

Targeting mRNA processing as an anticancer strategy

Discoveries in the past decade have highlighted the potential of mRNA as a therapeutic target for cancer. Specifically, RNA sequencing revealed that, in addition to gene mutations, alterations in mRNA can contribute to the initiation and progression of cancer. Indeed, precursor mRNA processing, which includes the removal of introns by splicing and the formation of 3' ends by cleavage and polyadenylation, is frequently altered in tumours. These alterations result in numerous cancer-specific mRNAs that generate altered levels of normal proteins or proteins with new functions, leading to the activation of oncogenes or the inactivation of tumour-suppressor genes. Abnormally spliced and polyadenylated mRNAs are also associated with resistance to cancer treatment and, unexpectedly, certain cancers are highly sensitive to the pharmacological inhibition of splicing. This Review summarizes recent progress in our understanding of how splicing and polyadenylation are altered in cancer and highlights how this knowledge has been translated for drug discovery, resulting in the production of small molecules and oligonucleotides that modulate the spliceosome and are in clinical trials for the treatment of cancer.

Systematic Profiling of Alternative mRNA Splicing Signature for Predicting Glioblastoma Prognosis

Emerging evidence suggests that alternative splicing (AS) is modified in cancer and is associated with cancer progression. Systematic analysis of AS signature in glioblastoma (GBM) is lacking and is greatly needed. We profiled genome-wide AS events in 498 GBM patients in TCGA using RNA-seq data, and splicing network and prognostic predictor were built by integrated bioinformatics analysis. Among 45,610 AS events in 10,434 genes, we detected 1,829 AS events in 1,311 genes, and 1,667 AS events in 1,146 genes that were significantly associated with overall survival and disease-free survival of GBM patients, respectively. Five potential feature genes, S100A4, ECE2, CAST, ASPH, and LY6K, were discovered after network mining as well as correlation analysis between AS and gene expression, most of which were related to carcinogenesis and development. Multivariate survival model analysis indicated that these five feature genes could classify the prognosis at AS event and gene expression level. This report opens up a new avenue for exploration of the pathogenesis of GBM through AS, thus more precisely guiding clinical treatment and prognosis judgment.

A 54-kDa short variant of DHX33 functions in regulating mRNA translation

DEAH box protein DHX33 has been found to be necessary for cell proliferation and early development of multicellular organisms. It plays diverse roles in regulating gene transcription, ribosome RNA synthesis, and protein translation. Dysregulation of DHX33 has been observed in various human cancers. In this study, we identified a short DHX33 variant in cells. The short DHX33 (hereafter referred to as DHX33-2) has only 534 amino acids, which completely matches the C-terminal helicase domain of full-length DHX33 (DHX33-1). Different from DHX33-1, which mainly localizes to the nucleus, DHX33-2 preferentially localizes to the cytoplasm. Through protein immunoprecipitation and RNA- immunoprecipitation analysis, we found that DHX33-2 interacts with DDX3, eIF3, hnRNPs, poly (A) binding protein, and a subset of mRNAs. Further RNA sequencing analysis showed that DHX33 binds to a subset of mRNAs important in cell proliferation. DHX33-2 stimulates the translation for specific mRNAs. Our study for the first time demonstrates the function of a short DHX33 variant in protein translation.

Alternative Splicing in the Regulation of Plant-Microbe Interactions

As sessile organisms, plants are continuously exposed to a wide range of biotic interactions. While some biotic interactions are beneficial or even essential for the plant (e.g. rhizobia and mycorrhiza), others such as pathogens are detrimental and require fast adaptation. Plants partially achieve this growth and developmental plasticity by modulating the repertoire of genes they express. In the past few years, high-throughput transcriptome sequencing have revealed that, in addition to transcriptional control of gene expression, post-transcriptional processes, notably alternative splicing (AS), emerged as a key mechanism for gene regulation during plant adaptation to the environment. AS not only can increase proteome diversity by generating multiple transcripts from a single gene but also can reduce gene expression by yielding isoforms degraded by mechanisms such as nonsense-mediated mRNA decay. In this review, we will summarize recent discoveries detailing the contribution of AS to the regulation of plant-microbe interactions, with an emphasis on the modulation of immunity receptor function and other components of the signaling pathways that deal with pathogen responses. We will also discuss emerging evidences that AS could contribute to dynamic reprogramming of the plant transcriptome during beneficial interactions, such as the legume-symbiotic interaction.

SRSF3-Regulated RNA Alternative Splicing Promotes Glioblastoma Tumorigenicity by Affecting Multiple Cellular Processes

Misregulated alternative RNA splicing (AS) contributes to the tumorigenesis and progression of human cancers, including glioblastoma (GBM). Here, we showed that a major splicing factor, serine and arginine rich splicing factor 3 (SRSF3), was frequently upregulated in clinical glioma specimens and that elevated SRSF3 was associated with tumor progression and a poor prognosis for patients with glioma. In patient-derived glioma stem-like cells (GSC), SRSF3 expression promoted cell proliferation, self-renewal, and tumorigenesis. Transcriptomic profiling identified more than 1,000 SRSF3-affected AS events, with a preference for exon skipping in genes involved with cell mitosis. Motif analysis identified the sequence of CA(G/C/A)CC(C/A) as a potential exonic splicing enhancer for these SRSF3-regulated exons. To evaluate the biological impact of SRSF3-affected AS events, four candidates were selected whose AS correlated with SRSF3 expression in glioma tissues, and their splicing pattern was modified using a CRISPR/Cas9 approach. Two functionally validated AS candidates were further investigated for the mechanisms underlying their isoform-specific functions. Specifically, following knockout of SRSF3, transcription factor ETS variant 1 (ETV1) gene showed exon skipping at exon 7, while nudE neurodevelopment protein 1 (NDE1) gene showed replacement of terminal exon 9 with a mutually exclusive exon 9'. SRSF3-regulated AS of these two genes markedly increased their oncogenic activity in GSCs. Taken together, our data demonstrate that SRSF3 is a key regulator of AS in GBM and that understanding mechanisms of misregulated AS could provide critical insights for developing effective therapeutic strategies against GBMs. SIGNIFICANCE: SRSF3 is a significant regulator of glioma-associated alternative splicing, implicating SRSF3 as an oncogenic factor that contributes to the tumor biology of GBM.

A high resolution A-to-I editing map in the mouse identifies editing events controlled by pre-mRNA splicing

Pre-mRNA-splicing and adenosine to inosine (A-to-I) RNA-editing occur mostly cotranscriptionally. During A-to-I editing, a genomically encoded adenosine is deaminated to inosine by adenosine deaminases acting on RNA (ADARs). Editing-competent stems are frequently formed between exons and introns. Consistently, studies using reporter assays have shown that splicing efficiency can affect editing levels. Here, we use Nascent-seq and identify ∼90,000 novel A-to-I editing events in the mouse brain transcriptome. Most novel sites are located in intronic regions. Unlike previously assumed, we show that both ADAR (ADAR1) and ADARB1 (ADAR2) can edit repeat elements and regular transcripts to the same extent. We find that inhibition of splicing primarily increases editing levels at hundreds of sites, suggesting that reduced splicing efficiency extends the exposure of intronic and exonic sequences to ADAR enzymes. Lack of splicing factors NOVA1 or NOVA2 changes global editing levels, demonstrating that alternative splicing factors can modulate RNA editing. Finally, we show that intron retention rates correlate with editing levels across different brain tissues. We therefore demonstrate that splicing efficiency is a major factor controlling tissue-specific differences in editing levels.

A Computational Analysis of Alternative Splicing across Mammalian Tissues Reveals Circadian and Ultradian Rhythms in Splicing Events

Mounting evidence points to a role of the circadian clock in the temporal regulation of post-transcriptional processes in mammals, including alternative splicing (AS). In this study, we carried out a computational analysis of circadian and ultradian rhythms on the transcriptome level to characterise the landscape of rhythmic AS events in published datasets covering 76 tissues from mouse and olive baboon. Splicing-related genes with 24-h rhythmic expression patterns showed a bimodal distribution of peak phases across tissues and species, indicating that they might be controlled by the circadian clock. On the output level, we identified putative oscillating AS events in murine microarray data and pairs of differentially rhythmic splice isoforms of the same gene in baboon RNA-seq data that peaked at opposing times of the day and included oncogenes and tumour suppressors. We further explored these findings using a new circadian RNA-seq dataset of human colorectal cancer cell lines. Rhythmic isoform expression patterns differed between the primary tumour and the metastatic cell line and were associated with cancer-related biological processes, indicating a functional role of rhythmic AS that might be implicated in tumour progression. Our data shows that rhythmic AS events are widespread across mammalian tissues and might contribute to a temporal diversification of the proteome.

The RNA-binding protein QKI controls alternative splicing in vascular cells, producing an effective model for therapy

Dysfunction of endothelial cells (ECs) and vascular smooth muscle cells (VSMCs) leads to ischaemia, the central pathology of cardiovascular disease. Stem cell technology will revolutionise regenerative medicine, but a need remains to understand key mechanisms of vascular differentiation. RNA-binding proteins have emerged as novel post-transcriptional regulators of alternative splicing and we have previously shown that the RNA-binding protein Quaking (QKI) plays roles in EC differentiation. In this study, we decipher the role of the alternative splicing isoform Quaking 6 (QKI-6) to induce VSMC differentiation from induced pluripotent stem cells (iPSCs). PDGF-BB stimulation induced QKI-6, which bound to HDAC7 intron 1 via the QKI-binding motif, promoting HDAC7 splicing and iPS-VSMC differentiation. Overexpression of QKI-6 transcriptionally activated SM22 (also known as TAGLN), while QKI-6 knockdown diminished differentiation capability. VSMCs overexpressing QKI-6 demonstrated greater contractile ability, and upon combination with iPS-ECs-overexpressing the alternative splicing isoform Quaking 5 (QKI-5), exhibited higher angiogenic potential in vivo than control cells alone. This study demonstrates that QKI-6 is critical for modulation of HDAC7 splicing, regulating phenotypically and functionally robust iPS-VSMCs. These findings also highlight that the QKI isoforms hold key roles in alternative splicing, giving rise to cells which can be used in vascular therapy or for disease modelling.This article has an associated First Person interview with the first author of the paper.

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.

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.

Protein-protein interaction analysis highlights the role of septins in membrane enclosed lumen and mRNA processing

Septins are a family of GTP-binding proteins that assemble into non-polar filaments which can be recruited to negatively charged membranes and serve as a scaffold to recruit cytosolic proteins and cytoskeletal elements such as microtubules and actin so that they can perform their important biological functions. Human septins consist of four groups, each with 13 members, and filaments formation usually involve members from each group in specific positions. However, little is known about the molecular mechanisms that drive the binding of septins to membranes and its importance to their biological functions. Here we have built a protein-protein interaction (PPI) network around human septins and highlighted the connections with 170 partners. Functional enrichment by inference of the network of septins and their partners revealed their participation in functions consistent with some of the roles described for septins, including cell cycle, cell division and cell shape, but we also identified septin partners in these functions that had not previously been described. Interestingly, we identified important and multiple connections between septins and mRNA processing and their export from the nucleus. Analysis of the enrichment of gene ontology cellular components highlighted some important interactions between molecules involved in the spliceosome with septin 2 and septin 7 in particular. RNA splicing regulates gene expression, and through it, cell fate, development and physiology. Mutations in components of the in the splicing machinery is linked to several diseases including cancer, thus taken together, the different analyses presented here open new perspectives to elucidate the pathobiological role of septins.

Alternative splicing is required for stage differentiation in malaria parasites

Background

In multicellular organisms, alternative splicing is central to tissue differentiation and identity. Unicellular protists lack multicellular tissue but differentiate into variable cell types during their life cycles. The role of alternative splicing in transitions between cell types and establishing cellular identity is currently unknown in any unicellular organism.

Results

To test whether alternative splicing in unicellular protists plays a role in cellular differentiation, we conduct RNA-seq to compare splicing in female and male sexual stages to asexual intraerythrocytic stages in the rodent malaria parasite Plasmodium berghei. We find extensive changes in alternative splicing between stages and a role for alternative splicing in sexual differentiation. Previously, general gametocyte differentiation was shown to be modulated by specific transcription factors. Here, we show that alternative splicing establishes a subsequent layer of regulation, controlling genes relating to consequent sex-specific differentiation of gametocytes.

Conclusions

We demonstrate that alternative splicing is reprogrammed during cellular differentiation of a unicellular protist. Disruption of an alternative splicing factor, PbSR-MG, perturbs sex-specific alternative splicing and decreases the ability of the parasites to differentiate into male gametes and oocysts, thereby reducing transmission between vertebrate and insect hosts. Our results reveal alternative splicing as an integral, stage-specific phenomenon in these protists and as a regulator of cellular differentiation that arose early in eukaryotic evolution.

Electronic supplementary material

The online version of this article (10.1186/s13059-019-1756-6) contains supplementary material, which is available to authorized users.

Alternative mRNA splicing in cancer immunotherapy

Immunotherapies are yielding effective treatments for several previously untreatable cancers. Still, the identification of suitable antigens specific to the tumour that can be targets for cancer vaccines and T cell therapies is a challenge. Alternative processing of mRNA, a phenomenon that has been shown to alter the proteomic diversity of many cancers, may offer the potential of a broadened target space. Here, we discuss the promise of analysing mRNA processing events in cancer cells, with an emphasis on mRNA splicing, for the identification of potential new targets for cancer immunotherapy. Further, we highlight the challenges that must be overcome for this new avenue to have clinical applicability.

A Novel Splice Variant of the Masculinizing Gene Masc with piRNA-Cleavage-Site Defect Functions in Female External Genital Development in the Silkworm, Bombyx mori

In the silkworm, the sex-determination primary signal Fem controls sex differentiation by specific binding of Fem-derived piRNA to the cleavage site in Masc mRNA, thus inhibiting Masc protein production in the female. In this study, we identified a novel splicing isoform of Masc, named Masc-S, which lacks the intact sequence of the cleavage site, encoding a C-terminal truncated protein. Results of RT-PCR showed that Masc-S was expressed in both sexes. Over-expression of Masc-S and Masc in female-specific cell lines showed that Masc-S could be translated against the Fem-piRNA cut. By RNA-protein pull-down, LC/MS/MS, and EMSA, we identified a protein BmEXU that specifically binds to an exclusive RNA sequence in Masc compared to Masc-S. Knockdown of Masc-S resulted in abnormal morphology in female external genital and increased expression of the Hox gene Abd-B, which similarly occurred by Bmexu RNAi. These results suggest that the splice variant Masc-S against Fem-piRNA plays an important role in female external genital development, of which function is opposite to that of full-length Masc. Our study provides new insights into the regulatory mechanism of sex determination in the silkworm.

CAGI 5 splicing challenge: Improved exon skipping and intron retention predictions with MMSplice

Pathogenic genetic variants often primarily affect splicing. However, it remains difficult to quantitatively predict whether and how genetic variants affect splicing. In 2018, the fifth edition of the Critical Assessment of Genome Interpretation proposed two splicing prediction challenges based on experimental perturbation assays: Vex-seq, assessing exon skipping, and MaPSy, assessing splicing efficiency. We developed a modular modeling framework, MMSplice, the performance of which was among the best on both challenges. Here we provide insights into the modeling assumptions of MMSplice and its individual modules. We furthermore illustrate how MMSplice can be applied in practice for individual genome interpretation, using the MMSplice VEP plugin and the Kipoi variant interpretation plugin, which are directly applicable to VCF files.

An engineered RNA binding protein with improved splicing regulation

The muscleblind-like (MBNL) family of proteins are key developmental regulators of alternative splicing. Sequestration of MBNL proteins by expanded CUG/CCUG repeat RNA transcripts is a major pathogenic mechanism in the neuromuscular disorder myotonic dystrophy (DM). MBNL1 contains four zinc finger (ZF) motifs that form two tandem RNA binding domains (ZF1-2 and ZF3-4) which each bind YGCY RNA motifs. In an effort to determine the differences in function between these domains, we designed and characterized synthetic MBNL proteins with duplicate ZF1-2 or ZF3-4 domains, referred to as MBNL-AA and MBNL-BB, respectively. Analysis of splicing regulation revealed that MBNL-AA had up to 5-fold increased splicing activity while MBNL-BB had 4-fold decreased activity compared to a MBNL protein with the canonical arrangement of zinc finger domains. RNA binding analysis revealed that the variations in splicing activity are due to differences in RNA binding specificities between the two ZF domains rather than binding affinity. Our findings indicate that ZF1-2 drives splicing regulation via recognition of YGCY RNA motifs while ZF3-4 acts as a general RNA binding domain. Our studies suggest that synthetic MBNL proteins with improved or altered splicing activity have the potential to be used as both tools for investigating splicing regulation and protein therapeutics for DM and other microsatellite diseases.

Repeat-associated RNA structure and aberrant splicing

Over 30 hereditary disorders attributed to the expansion of microsatellite repeats have been identified. Despite variant nucleotide content, number of consecutive repeats, and different locations in the genome, many of these diseases have pathogenic RNA gain-of-function mechanisms. The repeat-containing RNAs can form structures in vitro predicted to contribute to the disease through assembly of intracellular RNA aggregates termed foci. The expanded repeat RNAs within these foci sequester RNA binding proteins (RBPs) with important roles in the regulation of RNA metabolism, most notably alternative splicing (AS). These deleterious interactions lead to downstream alterations in transcriptome-wide AS directly linked with disease symptoms. This review summarizes existing knowledge about the association between the repeat RNA structures and RBPs as well as the resulting aberrant AS patterns, specifically in the context of myotonic dystrophy. The connection between toxic, structured RNAs and dysregulation of AS in other repeat expansion diseases is also discussed. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

RNA processing errors triggered by cadmium and integrator complex disruption are signals for environmental stress

BACKGROUND

Adaptive responses to stress are essential for cell and organismal survival. In metazoans, little is known about the impact of environmental stress on RNA homeostasis.

RESULTS

By studying the regulation of a cadmium-induced gene named numr-1 in Caenorhabditis elegans, we discovered that disruption of RNA processing acts as a signal for environmental stress. We find that NUMR-1 contains motifs common to RNA splicing factors and influences RNA splicing in vivo. A genome-wide screen reveals that numr-1 is strongly and specifically induced by silencing of genes that function in basal RNA metabolism including subunits of the metazoan integrator complex. Human integrator processes snRNAs for functioning with splicing factors, and we find that silencing of C. elegans integrator subunits disrupts snRNA processing, causes aberrant pre-mRNA splicing, and induces the heat shock response. Cadmium, which also strongly induces numr-1, has similar effects on RNA and the heat shock response. Lastly, we find that heat shock factor-1 is required for full numr-1 induction by cadmium.

CONCLUSION

Our results are consistent with a model in which disruption of integrator processing of RNA acts as a molecular damage signal initiating an adaptive stress response mediated by heat shock factor-1. When numr-1 is induced via this pathway in C. elegans, its function in RNA metabolism may allow it to mitigate further damage and thereby promote tolerance to cadmium.

How RNA structure dictates the usage of a critical exon of spinal muscular atrophy gene

Role of RNA structure in pre-mRNA splicing has been implicated for several critical exons associated with genetic disorders. However, much of the structural studies linked to pre-mRNA splicing regulation are limited to terminal stem-loop structures (hairpins) sequestering splice sites. In few instances, role of long-distance interactions is implicated as the major determinant of splicing regulation. With the recent surge of reports of circular RNA (circRNAs) generated by backsplicing, role of Alu-associated RNA structures formed by long-range interactions are taking central stage. Humans contain two nearly identical copies of Survival Motor Neuron (SMN) genes, SMN1 and SMN2. Deletion or mutation of SMN1 coupled with the inability of SMN2 to compensate for the loss of SMN1 due to exon 7 skipping causes spinal muscular atrophy (SMA), one of the leading genetic diseases of children. In this review, we describe how structural elements formed by both local and long-distance interactions are being exploited to modulate SMN2 exon 7 splicing as a potential therapy for SMA. We also discuss how Alu-associated secondary structure modulates generation of a vast repertoire of SMN circRNAs. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

Conserved sequences in the Drosophila mod(mdg4) intron promote poly(A)-independent transcription termination and trans-splicing

Alternative splicing (AS) is a regulatory mechanism of gene expression that greatly expands the coding capacities of genomes by allowing the generation of multiple mRNAs from a single gene. In Drosophila, the mod(mdg4) locus is an extreme example of AS that produces more than 30 different mRNAs via trans-splicing that joins together the common exons and the 3′ variable exons generated from alternative promoters. To map the regions required for trans-splicing, we have developed an assay for measuring trans-splicing events and identified a 73-bp region in the last common intron that is critical for trans-splicing of three pre-mRNAs synthesized from different DNA strands. We have also found that conserved sequences in the distal part of the last common intron induce polyadenylation-independent transcription termination and are enriched by paused RNA polymerase II (RNAP II). These results suggest that all mod(mdg4) mRNAs are formed by joining in trans the 5′ splice site in the last common exon with the 3′ splice site in one of the alternative exons.

A Comprehensive Analysis of Ontogeny of Renal Drug Transporters: mRNA Analyses, Quantitative Proteomics, and Localization

Human renal membrane transporters play key roles in the disposition of renally cleared drugs and endogenous substrates, but their ontogeny is largely unknown. Using 184 human postmortem frozen renal cortical tissues (preterm newborns to adults) and a subset of 62 tissue samples, we measured the mRNA levels of 11 renal transporters and the transcription factor pregnane X receptor (PXR) with quantitative real‐time polymerase chain reaction, and protein abundance of nine transporters using liquid chromatography tandem mass spectrometry selective reaction monitoring, respectively. Expression levels of p‐glycoprotein, urate transporter 1, organic anion transporter 1, organic anion transporter 3, and organic cation transporter 2 increased with age. Protein levels of multidrug and toxin extrusion transporter 2‐K and breast cancer resistance protein showed no difference from newborns to adults, despite age‐related changes in mRNA expression. Multidrug and toxin extrusion transporter 1, glucose transporter 2, multidrug resistance‐associated protein 2, multidrug resistance‐associated protein 4 (MRP4), and PXR expression levels were stable. Using immunohistochemistry, we found that MRP4 localization in pediatric samples was similar to that in adult samples. Collectively, our study revealed that renal drug transporters exhibited different rates and patterns of maturation, suggesting that renal handling of substrates may change with age.

More than a messenger: Alternative splicing as a therapeutic target

Alternative splicing of pre-mRNA is an essential post- and co-transcriptional mechanism of gene expression regulation that produces multiple mature mRNA transcripts from a single gene. Genetic mutations that affect splicing underlie numerous devastating diseases. The complexity of splicing regulation allows for multiple therapeutic approaches to correct disease-associated mis-splicing events. In this review, we first highlight recent findings from therapeutic strategies that have used splice switching antisense oligonucleotides and small molecules that bind directly to RNA. Second, we summarize different genetic and chemical approaches to target components of the spliceosome to correct splicing defects in pathological conditions. Finally, we present an overview of compounds that target kinases and accessory pathways that intersect with the splicing machinery. Advancements in the understanding of disease-specific defects caused by mis-regulation of alternative splicing will certainly increase the development of therapeutic options for the clinic. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

ALG-1 Influences Accurate mRNA Splicing Patterns in the Caenorhabditis elegans Intestine and Body Muscle Tissues by Modulating Splicing Factor Activities

MicroRNAs (miRNAs) are known to modulate gene expression, but their activity at the tissue-specific level remains largely uncharacterized. To study their contribution to tissue-specific gene expression, we developed novel tools to profile putative miRNA targets in the Caenorhabditis elegans intestine and body muscle. We validated many previously described interactions and identified ∼3500 novel targets. Many of the candidate miRNA targets curated are known to modulate the functions of their respective tissues. Within our data sets we observed a disparity in the use of miRNA-based gene regulation between the intestine and body muscle. The intestine contained significantly more putative miRNA targets than the body muscle highlighting its transcriptional complexity. We detected an unexpected enrichment of RNA-binding proteins targeted by miRNA in both tissues, with a notable abundance of RNA splicing factors. We developed in vivo genetic tools to validate and further study three RNA splicing factors identified as putative miRNA targets in our study (asd-2, hrp-2, and smu-2), and show that these factors indeed contain functional miRNA regulatory elements in their 3'UTRs that are able to repress their expression in the intestine. In addition, the alternative splicing pattern of their respective downstream targets (unc-60, unc-52, lin-10, and ret-1) is dysregulated when the miRNA pathway is disrupted. A reannotation of the transcriptome data in C. elegans strains that are deficient in the miRNA pathway from past studies supports and expands on our results. This study highlights an unexpected role for miRNAs in modulating tissue-specific gene isoforms, where post-transcriptional regulation of RNA splicing factors associates with tissue-specific alternative splicing.

Neural crest-specific deletion of Rbfox2 in mice leads to craniofacial abnormalities including cleft palate

Alternative splicing (AS) creates proteomic diversity from a limited size genome by generating numerous transcripts from a single protein-coding gene. Tissue-specific regulators of AS are essential components of the gene regulatory network, required for normal cellular function, tissue patterning, and embryonic development. However, their cell-autonomous function in neural crest development has not been explored. Here, we demonstrate that splicing factor Rbfox2 is expressed in the neural crest cells (NCCs), and deletion of Rbfox2 in NCCs leads to cleft palate and defects in craniofacial bone development. RNA-Seq analysis revealed that Rbfox2 regulates splicing and expression of numerous genes essential for neural crest/craniofacial development. We demonstrate that Rbfox2-TGF-β-Tak1 signaling axis is deregulated by Rbfox2 deletion. Furthermore, restoration of TGF-β signaling by Tak1 overexpression can rescue the proliferation defect seen in Rbfox2 mutants. We also identified a positive feedback loop in which TGF-β signaling promotes expression of Rbfox2 in NCCs.

Abnormalities affecting the head and face – such as a cleft lip or palate – are among the most common of all birth defects. These tissues normally develop from cells in the embryo known as the neural crest cells, and specifically a subset of these cells called the cranial neural crest cells. Most cases of cleft lip or palate are linked back to genes that affect the biology of this group of cells.

The list of genes implicated in the impaired development of cranial neural crest cells code for proteins with a wide range of different activities. Some encode transcription factors – proteins that switch genes on or off. Others code for chromatin remodeling factors, which control how the DNA is packed inside cells. However, the role of another group of proteins – the splicing factors – remains unclear and warrants further investigation.

When a gene is switched on its genetic code is first copied into a short-lived molecule called a transcript. These transcripts are then edited to form templates to build proteins. Splicing is one way that a transcript can be edited, which involves different pieces of the transcript being cut out and the remaining pieces being pasted together to form alternative versions of the final template. Splicing factors control this process.

Cibi et al. now show that neural crest cells from mice make a splicing factor called Rbfox2 and that deleting this gene for this protein from only these cells leads to mice with a cleft palate and defects in the bones of their head and face.

Further analysis helped to identify the transcripts that are spliced by Rbfox2, and the effects that these splicing events have on gene activity in mouse tissues that develop from cranial neural crest cells. Cibi et al. went on to find a signaling pathway that was impaired in the mutant cells that lacked Rbfox2. Forcing the mutant cells to over-produce one of the proteins involved in this signaling pathway (a protein named Tak1) was enough to compensate for the some of the defects caused by a lack of Rbfox2, suggesting it acts downstream of the splicing regulator.

Lastly, Cibi et al. showed that another protein in this signaling pathway, called TGF-β, acted to increase how much Rbfox2 was made by neural crest cells. Together these findings may be relevant in human disease studies, given that altered TGF-β signaling is a common feature in many birth defects seen in humans.

Predicting the change of exon splicing caused by genetic variant using support vector regression

Alternative splicing can be disrupted by genetic variants that are related to diseases like cancers. Discovering the influence of genetic variations on the alternative splicing will improve the understanding of the pathogenesis of variants. Here, we developed a new approach, PredPSI-SVR to predict the impact of variants on exon skipping events by using the support vector regression. From the sequence of a particular exon and its flanking regions, 42 comprehensive features related to splicing events were extracted. By using a greedy feature selection algorithm, we found eight features contributing most to the prediction. The trained model achieved a Pearson correlation coefficient (PCC) of 0.570 in the 10-fold cross-validation based on the training data set provided by the "vex-seq" challenge of the 5th Critical Assessment of Genome Interpretation. In the blind test also held by the challenge, our prediction ranked the 2nd with a PCC of 0.566 that demonstrates the robustness of our method. A further test indicated that the PredPSI-SVR is helpful in prioritizing deleterious synonymous mutations. The method is available on https://github.com/chenkenbio/PredPSI-SVR.

Alternative splicing events during adipogenesis from hMSCs

Adipogenesis, the developmental process of progenitor-cell differentiating into adipocytes, leads to fat metabolic disorders. Alternative splicing (AS), a ubiquitous regulatory mechanism of gene expression, allows the generation of more than one unique messenger RNA (mRNA) species from a single gene. Till now, alternative splicing events during adipogenesis from human mesenchymal stem cells (hMSCs) are not yet fully elucidated. We performed RNA-Seq coupled with bioinformatics analysis to identify the differentially expressed AS genes and events during adipogenesis from hMSCs. A global survey separately identified 1262, 1181, 1167, and 1227 ASE involved in the most common types of AS including cassette exon, alt3, and alt5, especially with cassette exon the most prevalent, at 7, 14, 21, and 28 days during adipogenesis. Interestingly, 122 differentially expressed ASE referred to 118 genes, and the three genes including ACTN1 (alt3 and cassette), LRP1 (alt3 and alt5), and LTBP4 (cassette, cassette_multi, and unknown), appeared in multiple AS types of ASE during adipogenesis. Except for all the identified ASE of LRP1 occurred in the extracellular topological domain, alt3 (84) in transmembrane domain significantly differentially expressed was the potential key event during adipogenesis. Overall, we have, for the first time, conducted the global transcriptional profiling during adipogenesis of hMSCs to identify differentially expressed ASE and ASE-related genes. This finding would provide extensive ASE as the regulator of adipogenesis and the potential targets for future molecular research into adipogenesis-related metabolic disorders.

Genome-wide profiling of the alternative splicing provides insights into development in Plutella xylostella

Background

The diamondback moth (DBM), Plutella xylostella (L.), is a major pest of cruciferous crops worldwide. While the species has become a model for genomics, post-transcriptional mechanisms associated with development and sex determination have not been comprehensively studied and the lack of complete structure of mRNA transcripts limits further research.

Results

Here, we combined the methods of single-molecule long-read sequencing technology (IsoSeq) and RNA-seq to re-annotate the published DBM genome and present the genome-wide identification of alternative splicing (AS) associated with development and sex determination of DBM. In total, we identified ~ 13,900 genes (~ 77%) annotated in the DBM genome (version-2), resulting in the correction of 1586 wrongly annotated genes and identification of 78,000 previously unannotated transcripts. We also identified 1804 genes showing alternative splicing (AS) in each of the developmental stages and sexes, suggesting that AS events are ubiquitous in DBM. Comparative analyses showed that these AS events were rarely shared among developmental stages, indicating that they may play key specific roles in regulation of insect development. Further, we found 156 genes showing different AS events and expression patterns between males and females, linking them to potential functions in sex determination.

Conclusion

Overall, the P. xylostella transcriptome provides the significant information about regulatory alternative splicing events, which are shown to be involved in development and sex determination. Our work presents a solid foundation to better understand the mechanism of post-transcriptional regulation, and offers wider insights into insect development and sex determination.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5838-3) contains supplementary material, which is available to authorized users.

Alcohol exposure alters pre-mRNA splicing of antiapoptotic Mcl-1L isoform and induces apoptosis in neural progenitors and immature neurons

Alternative splicing and expression of splice variants of genes in the brain may lead to the modulation of protein functions, which may ultimately influence behaviors associated with alcohol dependence and neurotoxicity. We recently showed that ethanol exposure can lead to pre-mRNA missplicing of Mcl-1, a pro-survival member of the Bcl-2 family, by downregulating the expression levels of serine/arginine rich splicing factor 1 (SRSF1). Little is known about the physiological expression of these isoforms in neuronal cells and their role in toxicity induced by alcohol exposure during the developmental period. In order to investigate the impact of alcohol exposure on alternative splicing of Mcl-1 pre-mRNA and its role in neurotoxicity, we developed a unique primary human neuronal culture model where neurospheres (hNSPs), neural progenitors (hNPCs), immature neurons, and mature neurons were cultured from the matching donor fetal brain tissues. Our data suggest that neural progenitors and immature neurons are highly sensitive to the toxic effects of ethanol, while mature neuron cultures showed resistance to ethanol exposure. Further analysis of Mcl-1 pre-mRNA alternative splicing by semi-quantitative and quantitative analysis revealed that ethanol exposure causes a significant decrease in Mcl-1L/Mcl-1S ratio in a dose and time dependent manner in neural progenitors. Interestingly, ectopic expression of Mcl-1L isoform in neural progenitors was able to recover the viability loss and apoptosis induced by alcohol exposure. Altogether, these observations suggest that alternative splicing of Mcl-1 may play a crucial role in neurotoxicity associated with alcohol exposure in the developing fetal brain.

Increased intron retention is a post-transcriptional signature associated with progressive aging and Alzheimer's disease

Intron retention (IR) by alternative splicing is a conserved regulatory mechanism that can affect gene expression and protein function during adult development and age-onset diseases. However, it remains unclear whether IR undergoes spatial or temporal changes during different stages of aging or neurodegeneration like Alzheimer's disease (AD). By profiling the transcriptome of Drosophila head cells at different ages, we observed a significant increase in IR events for many genes during aging. Differential IR affects distinct biological functions at different ages and occurs at several AD-associated genes in older adults. The increased nucleosome occupancy at the differentially retained introns in young animals suggests that it may regulate the level of IR during aging. Notably, an increase in the number of IR events was also observed in healthy older mouse and human brain tissues, as well as in the cerebellum and frontal cortex from independent AD cohorts. Genes with differential IR shared many common features, including shorter intron length, no perturbation in their mRNA level, and enrichment for biological functions that are associated with mRNA processing and proteostasis. The differentially retained introns identified in AD frontal cortex have higher GC content, with many of their mRNA transcripts showing an altered level of protein expression compared to control samples. Taken together, our results suggest that an increased IR is an conserved signature that is associated with aging. By affecting pathways involved in mRNA and protein homeostasis, changes of IR pattern during aging may regulate the transition from healthy to pathological state in late-onset sporadic AD.

N6-methyladenosine modifications: interactions with novel RNA-binding proteins and roles in signal transduction

RNA epigenetics has received a great deal of attention in recent years, and the reversible N6-methyladenosine (m6A) modification on messenger RNAs (mRNAs) has emerged as a widespread phenomenon. The vital roles of m6A in diverse biological processes are dependent on many RNA-binding proteins (RBPs) with 'reader' or 'nonreader' functions. Moreover, m6A effector proteins affect cellular processes, such as stem cell differentiation, tumor development and the immune response by controlling signal transduction. This review provides an overview of the interactions of m6A with various RBPs, including the 'reader' proteins (excluding the YT521-B homology (YTH) domain proteins and the heterogeneous nuclear ribonucleoproteins (hnRNPs)), and the functional 'nonreader' proteins, and this review focuses on their specific RNA-binding domains and their associations with other m6A effectors. Furthermore, we summarize key m6A-marked targets in distinct signaling pathways, leading to a better understanding of the cellular m6A machinery.

Alternative splicing coupled to nonsense-mediated mRNA decay contributes to the high-altitude adaptation of maca (Lepidium meyenii)

Alpine plants remain the least studied plant communities in terrestrial ecosystems. However, how they adapt to high-altitude environments is far from clear. Here, we used RNA-seq to investigate a typical alpine plant maca (Lepidium meyenii) to understand its high-altitude adaptation at transcriptional and post-transcriptional level. At transcriptional level, we found that maca root significantly up-regulated plant immunity genes in day-time comparing to night-time, and up-regulated abiotic (cold/osmotic) stress response genes in Nov and Dec comparing to Oct. In addition, 17 positively selected genes were identified, which could be involved in mitochondrion. At post-transcriptional level, we found that maca had species-specific characterized alternative splicing (AS) profile which could be influenced by stress environments. For example, the alternative 3' splice site events (A3SS, 39.62%) were predominate AS events in maca, rather than intron retention (IR, 23.17%). Interestingly, besides serine/arginine-rich (SR) proteins and long non-coding RNAs (lncRNAs), a lot of components in nonsense-mediated mRNA decay (NMD) were identified under differential alternative splicing (DAS), supporting AS coupled to NMD as essential mechanisms for maca's stress responses and high-altitude adaptation. Taken together, we first attempted to unveil maca's high-altitude adaptation mechanisms based on transcriptome and post-transcriptome evidence. Our data provided valuable insights to understand the high-altitude adaptation of alpine plants.

Function and Evolution of Upstream ORFs in Eukaryotes

There is growing interest in the role of translational regulation in cellular homeostasis during organismal development. Translation initiation is the rate-limiting step in mRNA translation and is central to translational regulation. Upstream open reading frames (uORFs) are regulatory elements that are prevalent in eukaryotic mRNAs. uORFs modulate the translation initiation rate of downstream coding sequences (CDSs) by sequestering ribosomes. Over the past several years, genome-wide studies have revealed the widespread regulatory functions of uORFs in different species in different biological contexts. Here, we review the current understanding of uORF-mediated translational regulation from the perspective of functional and evolutionary genomics and address remaining gaps that deserve further study.

Analysis of NKX2-5 in 439 Chinese Patients with Sporadic Atrial Septal Defect

Background

The NKX2 gene family is made up of core transcription factors that are involved in the morphogenesis of the vertebrate heart. NKx2-5 plays a pivotal role in mouse cardiogenesis, and mutations in NKx2-5 result in an abnormal structure and function of the heart, including atrial septal defect and cardiac electrophysiological abnormalities.

Material/Methods

To investigate the genetic variation of NKX2-5 in Chinese patients with sporadic atrial septal defect, we sequenced the full length of the NKX2-5 gene in the participants of the study. Four hundred thirty-nine patients and 567 healthy unrelated individuals were recruited. Genomic DNA was extracted from the peripheral blood leukocytes of the participants. DNA samples from the participants were amplified by multiplex PCR and sequenced on an Illumina HiSeq platform. Variations were detected by comparison with a standard reference genome and annotation with a variant effect predictor.

Results

Thirty variations were detected in Chinese patients with sporadic atrial septal defect, and 6 single nucleotide polymorphisms (SNPs) had a frequency greater than 1%. Among the 30 variations, the SNPs rs2277923 and rs3729753 were extremely prominent, with a high frequency and odds ratio in patients.

Conclusions

Single nucleotide variations are the prominent genetic variations of NKX2-5 in Chinese patients with sporadic atrial septal defect. The SNPs rs2277923 and rs3729753 are prominent single nucleotide variations (SNVs) in Chinese patients with sporadic atrial septal defect.

IsoformSwitchAnalyzeR: analysis of changes in genome-wide patterns of alternative splicing and its functional consequences

Summary

Alternative splicing is an important mechanism involved in health and disease. Recent work highlights the importance of investigating genome-wide changes in splicing patterns and the subsequent functional consequences. Current computational methods only support such analysis on a gene-by-gene basis. Therefore, we extended IsoformSwitchAnalyzeR R library to enable analysis of genome-wide changes in specific types of alternative splicing and predicted functional consequences of the resulting isoform switches. As a case study, we analyzed RNA-seq data from The Cancer Genome Atlas and found systematic changes in alternative splicing and the consequences of the associated isoform switches.

Availability and implementation

Windows, Linux and Mac OS: http://bioconductor.org/packages/IsoformSwitchAnalyzeR.

Supplementary information

Supplementary data are available at Bioinformatics online.

Bioinformatics challenges and perspectives when studying the effect of epigenetic modifications on alternative splicing

It is widely known that epigenetic modifications are important in regulating transcription, but several have also been reported in alternative splicing. The regulation of pre-mRNA splicing is important to explain proteomic diversity and the misregulation of splicing has been implicated in many diseases. Here, we give a brief overview of the role of epigenetics in alternative splicing and disease. We then discuss the bioinformatics methods that can be used to model interactions between epigenetic marks and regulators of splicing. These models can be used to identify alternative splicing and epigenetic changes across different phenotypes.

This article is part of a discussion meeting issue ‘Frontiers in epigenetic chemical biology’.

An alternative CTCF isoform antagonizes canonical CTCF occupancy and changes chromatin architecture to promote apoptosis

CTCF plays key roles in gene regulation, chromatin insulation, imprinting, X chromosome inactivation and organizing the higher-order chromatin architecture of mammalian genomes. Previous studies have mainly focused on the roles of the canonical CTCF isoform. Here, we explore the functions of an alternatively spliced human CTCF isoform in which exons 3 and 4 are skipped, producing a shorter isoform (CTCF-s). Functionally, we find that CTCF-s competes with the genome binding of canonical CTCF and binds a similar DNA sequence. CTCF-s binding disrupts CTCF/cohesin binding, alters CTCF-mediated chromatin looping and promotes the activation of IFI6 that leads to apoptosis. This effect is caused by an abnormal long-range interaction at the IFI6 enhancer and promoter. Taken together, this study reveals a non-canonical function for CTCF-s that antagonizes the genomic binding of canonical CTCF and cohesin, and that modulates chromatin looping and causes apoptosis by stimulating IFI6 expression.

CTCF plays key roles in gene regulation, chromatin insulation and organizing the higher-order chromatin architecture of mammalian genomes. Here the authors investigate the function an alternatively spliced shorter CTCF isoform, finding that this isoform antagonizes canonical CTCF occupancy and changes chromatin architecture to promote apoptosis.

Production of an aberrant splice variant of CCL5 is not caused by genetic mutation in the mammary glands of mastitis‑infected Holstein cows

Genetic mutations, including single nucleotide polymorphisms (SNPs), result in aberrant alternatively splicing of gene and involves in susceptibility of inflammatory diseases, including bovine mastitis. C‑C motif chemokine ligand 5 (CCL5) is an immune‑associated gene, but its alternative splicing (AS) mechanism of gene expression has not yet been understood. The present study identified the splice variant of CCL5 and the compared differential expression of various transcripts between healthy and mastitic mammary gland tissue from cows. A novel transcript lacking exon 2 with a deletion of 112 bp (referred to as CCL5‑AS) was identified in the mammary gland. The expression of CCL5‑AS was lower compared with CCL5‑reference in the healthy and mastitic mammary tissues. A total of two novel SNPs (g.1647 C>T and g.1804 G>A) were identified in exon 2 of CCL5. Using the splicing mini‑gene reporter assay in bovine mammary epithelial cells (MAC‑T) and 293T cells, it was confirmed that the production of CCL5‑AS was not caused by the two SNPs. The present findings suggested that alternative splicing is one of the mechanisms of CCL5 expression regulation and is involved in mastitis infection, but that genetic mutation was not responsible for the generation of the abnormal transcript of CCL5 in cows.