The determinants of alternative RNA splicing in human cells

Neuron and astrocyte specific 5mC and 5hmC signatures of BDNF's receptor, TrkB

Brain derived neurotrophic factor (BDNF) is the most studied trophic factor in the central nervous system (CNS), and its role in the maturation of neurons, including synapse development and maintenance has been investigated intensely for over three decades. The primary receptor for BDNF is the tropomyosin receptor kinase B (TrkB), which is broadly expressed as two primary isoforms in the brain; the full length TrkB (TrkB.FL) receptor, expressed mainly in neurons and the truncated TrkB (TrkB.T1) receptor. We recently demonstrated that TrkB.T1 is predominately expressed in astrocytes, and appears critical for astrocyte morphological maturation. Given the critical role of BDNF/TrkB pathway in healthy brain development and mature CNS function, we aimed to identify molecular underpinnings of cell-type specific expression of each TrkB isoform. Using Nanopore sequencing which enables direct, long read sequencing of native DNA, we profiled DNA methylation patterns of the entire TrkB gene, Ntrk2, in both neurons and astrocytes. Here, we identified robust differences in cell-type specific isoform expression associated with significantly different methylation patterns of the Ntrk2 gene in each cell type. Notably, astrocytes demonstrated lower 5mC methylation, and higher 5hmC across the entire gene when compared to neurons, including differentially methylated sites (DMSs) found in regions flanking the unique TrkB.T1 protein coding sequence (CDS). These data suggest DNA methylation patterns may provide instruction for isoform specific TrkB expression across unique CNS cell types.

Decoding the Epigenetic Landscape: Insights into 5mC and 5hmC Patterns in Mouse Cortical Cell Types

The DNA modifications, 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC), represent powerful epigenetic regulators of temporal and spatial gene expression. Yet, how the cooperation of these genome-wide, epigenetic marks determine unique transcriptional signatures across different brain cell populations is unclear. Here we applied Nanopore sequencing of native DNA to obtain a complete, genome-wide, single-base resolution atlas of 5mC and 5hmC modifications in neurons, astrocytes and microglia in the mouse cortex (99% genome coverage, 40 million CpG sites). In tandem with RNA sequencing, analysis of 5mC and 5hmC patterns across cell types reveals astrocytes drive uniquely high brain 5hmC levels and support two decades of research regarding methylation patterns, gene expression and alternative splicing, benchmarking this resource. As such, we provide the most comprehensive DNA methylation data in mouse brain as an interactive, online tool (NAM-Me, https://olsenlab.shinyapps.io/NAMME/) to serve as a resource dataset for those interested in the methylome landscape.

CaTCH: Calculating transcript complexity of human genes

The findings based on whole transcriptome sequencing suggest that alternative splicing occurs in approximately 95% of human multi-exon genes, thus, playing a crucial role in promoting proteome diversity. According to the latest GENCODE annotations, most genes have less than four transcripts, positively correlating with the number of exons. Thus, it is more accurate to measure the splice variant efficiency of a gene with respect to the number of exons, which is a measure of Transcript Complexity (TC). In addition to that, the theoretical number of transcripts is substantially higher than the actual number of transcripts produced by Alternative Splicing Events, and the features restricting this phenomenon need to be explored. In this method, we have extracted the data of various features contributing to TC from different databases. Linear regression is used to identify the determinant features and to train and test the model of TC. The results indicate that exon length is the determining feature of TC, followed by coding potential, presence of chromatin signature, and 5' splice site dinucleotide, all of which negatively affect a gene's TC, except exon length. To further classify the genes based on TC, random forest is used to identify the determinant features.•The splicing efficiency of a gene can be inferred by the transcript complexity, which is the number of transcripts per exon.•CaTCH is a linear regression-based model to calculate the transcript complexity of human genes, which can be calculated from the exon length, coding potentiality, presence of chromatin signature/s, and 5' splice site dinucleotide.

The Many Roads from Alternative Splicing to Cancer: Molecular Mechanisms Involving Driver Genes

Cancer driver genes are either oncogenes or tumour suppressor genes that are classically activated or inactivated, respectively, by driver mutations. Alternative splicing-which produces various mature mRNAs and, eventually, protein variants from a single gene-may also result in driving neoplastic transformation because of the different and often opposed functions of the variants of driver genes. The present review analyses the different alternative splicing events that result in driving neoplastic transformation, with an emphasis on their molecular mechanisms. To do this, we collected a list of 568 gene drivers of cancer and revised the literature to select those involved in the alternative splicing of other genes as well as those in which its pre-mRNA is subject to alternative splicing, with the result, in both cases, of producing an oncogenic isoform. Thirty-one genes fall into the first category, which includes splicing factors and components of the spliceosome and splicing regulators. In the second category, namely that comprising driver genes in which alternative splicing produces the oncogenic isoform, 168 genes were found. Then, we grouped them according to the molecular mechanisms responsible for alternative splicing yielding oncogenic isoforms, namely, mutations in cis splicing-determining elements, other causes involving non-mutated cis elements, changes in splicing factors, and epigenetic and chromatin-related changes. The data given in the present review substantiate the idea that aberrant splicing may regulate the activation of proto-oncogenes or inactivation of tumour suppressor genes and details on the mechanisms involved are given for more than 40 driver genes.

Long pentraxin 3 and vitamin D receptor mRNA expression pattern of cumulus granulosa cells isolated from PCOS oocytes at different stages of nuclear maturation

BACKGROUND

A fine-tuned pro-inflammatory and anti-inflammatory balance in the follicular unit is essential for cumulus expansion and successful ovulation. While the long pentraxin 3 (PTX3) gene is required for the expansion of cumulus cells (CCs), ovulation, resumption of meiosis and fertilization, the vitamin D receptor gene (VDR-X2) is required for intra-follicle redox balance. This study was planned to determine the expression pattern of VDR-X2 and PTX3 mRNA in CCs isolated from germinal vesicle (GV), metaphase I (MI), and metaphase II (MII) oocytes of PCOS patients with ovulatory dysfunction.

METHODS

The relative expression of CC-PTX3 and CC-VDR-X2 mRNA were evaluated using qRT-PCR in a total of 79 CC samples collected from individual cumulus-oocyte complex of 40 infertile patients (20 PCOS and 20 non-PCOS normal responders) who underwent ovarian stimulation with the GnRH antagonist protocol.

RESULTS

Relative PTX3 mRNA expressions of CCMI-control and CCMII-control showed 3- and 9-fold significant upregulation compared to CCGV-control, respectively. The relative PTX3 mRNA expression of CCMII-control increased approximately three fold compared to CCMI-control. Compared to CCGV-pcos, a 3-fold increase was noted in the relative PTX3 mRNA expression of CCMI-pcos and an approximately 4-fold increase in the PTX3 mRNA expression of CCMII-pcos. Relative PTX3 mRNA expression values of CCMII-pcos and CCMI-pcos were similar. A 6-fold upregulation of relative PTX3 mRNA and a 4-fold upregulation of VDR-X2 mRNA were detected in CCMII-control compared to CCMII-pcos. CC-VDR-X2 expression patterns of the PCOS and control groups overlapped with the CC-PTX3 pattern. Fertilization rates of the PCOS group exhibiting failed transcript expression were similar to normal responders.

CONCLUSION

The fact that relative CC-PTX3 and CC-VDR mRNA expression does not increase during the transition from MI to MII stage in PCOS as in normal responders suggests that PTX3 and VDR expression may be defective in cumulus cells of PCOS patients with ovulatory dysfunction.

Analysis of alternative splicing in chicken macrophages infected with avian pathogenic E. coli (APEC)

Colibacillosis is a complex disease that caused by avian pathogenic Escherichia coli (APEC), resulting in huge economic loss to the global poultry industry and threatening to human health. Alternative splicing (AS) is a universal post-transcriptional regulatory mechanism, which can simultaneously produce many proteins from a single gene to involve in various diseases and individual development. Herein, we characterized genome-wide AS events in wild type macrophages (WT) and APEC infected macrophages (APEC) by high-throughput RNA sequencing technology. A total of 751 differentially expressed (DE) AS genes were identified in the comparison of APEC vs. WT, including 587 of SE, 114 of MXE, 25 of RI, 17 of A3 and 8 of A5 event. Functional analysis showed that these identified DE AS genes were involved in 'Endocytosis', 'p53 signaling pathway', 'MAPK signaling pathway', 'NOD-like receptor signaling pathway', 'Ubiquitin mediated proteolysis' and 'Focal adhesion' immune related pathways. In summary, we comprehensively investigate AS events during APEC infection. This study has expanded our understanding of the process of APEC infection and provided new insights for further treatment options for APEC infection.

The DNA binding high mobility group box protein family functionally binds RNA

Nucleic acid binding proteins regulate transcription, splicing, RNA stability, RNA localization, and translation, together tailoring gene expression in response to stimuli. Upon discovery, these proteins are typically classified as either DNA or RNA binding as defined by their in vivo functions; however, recent evidence suggests dual DNA and RNA binding by many of these proteins. High mobility group box (HMGB) proteins have a DNA binding HMGB domain, act as transcription factors and chromatin remodeling proteins, and are increasingly understood to interact with RNA as means to regulate gene expression. Herein, multiple layers of evidence that the HMGB family are dual DNA and RNA binding proteins is comprehensively reviewed. For example, HMGB proteins directly interact with RNA in vitro and in vivo, are localized to RNP granules involved in RNA processing, and their protein interactors are enriched in RNA binding proteins involved in RNA metabolism. Importantly, in cell-based systems, HMGB-RNA interactions facilitate protein-protein interactions, impact splicing outcomes, and modify HMGB protein genomic or cellular localization. Misregulation of these HMGB-RNA interactions are also likely involved in human disease. This review brings to light that as a family, HMGB proteins are likely to bind RNA which is essential to HMGB protein biology. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications.

Proteomics as a tool to study the osteoimmunomodulatory role of metallic ions in a sol-gel coating

The success of bone implants depends on the osteoimmunomodulatory (OIM) activity of the biomaterials in the interactions with the periimplantary tissues. Many in vitro tests have been conducted to evaluate the osteoimmunology effects of biomaterials. However, results of these tests have often been inconclusive. This study examines the properties of newly developed sol-gel coatings doped with two metal ions associated with bone regeneration, Ca and Zn. The study uses both proteomic methods and traditional in vitro assays. The results demonstrate that proteomics is an effective tool to scrutinize the OIM properties of the materials. Moreover, sol-gel coatings offer excellent base materials to evaluate the effects of metal ions on these properties. The obtained data highlight the highly tunable nature of sol-gel materials; studying the materials with different doping levels supplies valuable information on the interactions between the immune and bone-forming processes.

A wave of deep intronic mutations in X-linked Alport Syndrome

X-linked Alport syndrome (XLAS) is an inherited kidney disease caused exclusively by pathogenic variants in the COL4A5 gene. In 10-20% of cases, DNA sequencing of COL4A5 exons or flanking regions cannot identify molecular causes. Here, our objective was to use a transcriptomic approach to identify causative events in a group of 19 patients with XLAS without identified mutation by Alport gene panel sequencing. Bulk RNAseq and/or targeted RNAseq using a capture panel of kidney genes was performed. Alternative splicing events were compared to those of 15 controls by a developed bioinformatic score. When using targeted RNAseq, COL4A5 coverage was found to be 23-fold higher than with bulk RNASeq and revealed 30 significant alternative splicing events in 17 of the 19 patients. After computational scoring, a pathogenic transcript was found in all patients. A causative variant affecting COL4A5 splicing and absent in the general population was identified in all cases. Altogether, we developed a simple and robust method for identification of aberrant transcripts due to pathogenic deep-intronic COL4A5 variants. Thus, these variants, potentially targetable by specific antisense oligonucleotide therapies, were found in a high percentage of patients with XLAS in whom pathogenic variants were missed by conventional DNA sequencing.

Deciphering the role of alternative splicing as modulators of defense response in the MYMIV- Vigna mungo pathosystem

Alternative splicing (AS) is a crucial regulatory mechanism that impacts transcriptome and proteome complexity under stressful situations. Although its role in abiotic stresses is somewhat understood, our understanding of the mechanistic regulation of pre-mRNA splicing in plant-pathogen interaction is meagre. To comprehend this unexplored immune reprogramming mechanism, transcriptome profiles of Mungbean Yellow Mosaic India Virus (MYMIV)-resistant and susceptible Vigna mungo genotypes were analysed for AS genes that may underlie the resistance mechanism. Results revealed a repertoire of AS-isoforms accumulated during pathogenic infestation, with intron retention being the most common AS mechanism. Identification of 688 differential alternatively spliced (DAS) genes in the resistant host elucidates its robust antiviral response, whereas 322 DAS genes were identified in the susceptible host. Enrichment analyses confirmed DAS transcripts pertaining to stress, signalling, and immune system pathways have undergone maximal perturbations. Additionally, a strong regulation of the splicing factors has been observed both at transcriptional and post-transcriptional levels. qPCR validation of candidate DAS transcripts with induced expression upon MYMIV-infection demonstrated a competent immune response in the resistant background. The AS-impacted genes resulted either in partial/complete loss of functional domains or altered sensitivity to miRNA-mediated gene silencing. A complex regulatory module, miR7517-ATAF2, has been identified in an aberrantly spliced ATAF2 isoform that exposes an intronic miR7517 binding site, thereby suppressing the negative regulator to enhance defense reaction. The present study establishes AS as a non-canonical immune reprogramming mechanism that operates in parallel, thereby offering an alternative strategy for developing yellow mosaic-resistant V. mungo cultivars. This article is protected by copyright. All rights reserved.

Inefficient splicing of long non-coding RNAs is associated with higher transcript complexity in human and mouse

Recent reports show that long non-coding RNAs (lncRNAs) have inefficient splicing and fewer alternative splice variants than mRNAs. Here, we have explored the efficiency of lncRNAs and mRNAs in producing various splice variants, given the number of exons in humans and mice. Intriguingly, lncRNAs produce more splice variants per exon, referred to as Transcript Complexity, than mRNAs. Most lncRNA splice variants are the product of the alternative last exon and exon skipping. LncRNAs and mRNAs with higher transcript complexity have shorter intron lengths. Longer exon length and GC/AG at 5'/3' splice sites are associated with higher transcript complexity in lncRNAs. Lastly, our results indicate that inefficient splicing of lncRNAs may facilitate multiple introns splicing and, thus, more spliced products per exon.

Recruitment of a splicing factor to the nuclear lamina for its inactivation

Precursor messenger RNA splicing is a highly regulated process, mediated by a complex RNA-protein machinery, the spliceosome, that encompasses several hundred proteins and five small nuclear RNAs in humans. Emerging evidence suggests that the spatial organization of splicing factors and their spatio-temporal dynamics participate in the regulation of splicing. So far, methods to manipulate the spatial distribution of splicing factors in a temporally defined manner in living cells are missing. Here, we describe such an approach that takes advantage of a reversible chemical dimerizer, and outline the requirements for efficient, reversible re-localization of splicing factors to selected sub-nuclear compartments. In a proof-of-principle study, the partial re-localization of the PRPF38A protein to the nuclear lamina in HEK293T cells induced a moderate increase in intron retention. Our approach allows fast and reversible re-localization of splicing factors, has few side effects and can be applied to many splicing factors by fusion of a protein tag through genome engineering. Apart from the systematic analysis of the spatio-temporal aspects of splicing regulation, the approach has a large potential for the fast induction and reversal of splicing switches and can reveal mechanisms of splicing regulation in native nuclear environments.

Through the use of a reversible chemical dimerizer, the splicing factor PRPF38A is re-localized to the nuclear lamina, paving the way for a systematic analysis of spatio-temporal splicing regulation.

Post-Transcriptional Dynamics is Involved in Rapid Adaptation to Hypergravity in Jurkat T Cells

The transcriptome of human immune cells rapidly reacts to altered gravity in a highly dynamic way. We could show in previous experiments that transcriptional patterns show profound adaption after seconds to minutes of altered gravity. To gain further insight into these transcriptional alteration and adaption dynamics, we conducted a highly standardized RNA-Seq experiment with human Jurkat T cells exposed to 9xg hypergravity for 3 and 15 min, respectively. We investigated the frequency with which individual exons were used during transcription and discovered that differential exon usage broadly appeared after 3 min and became less pronounced after 15 min. Additionally, we observed a shift in the transcript pool from coding towards non-coding transcripts. Thus, adaption of gravity-sensitive differentially expressed genes followed a dynamic transcriptional rebound effect. The general dynamics were compatible with previous studies on the transcriptional effects of short hypergravity on human immune cells and suggest that initial up-regulatory changes mostly result from increased elongation rates. The shift correlated with a general downregulation of the affected genes. All chromosome bands carried homogenous numbers of gravity-sensitive genes but showed a specific tendency towards up- or downregulation. Altered gravity affected transcriptional regulation throughout the entire genome, whereby the direction of differential expression was strongly dependent on the structural location in the genome. A correlation analysis with potential mediators of the early transcriptional response identified a link between initially upregulated genes with certain transcription factors. Based on these findings, we have been able to further develop our model of the transcriptional response to altered gravity.

ErbB4 alternative splicing mediates fetal mouse alveolar type II cell differentiation in vitro

BACKGROUND

Alternative splicing (AS) creates different protein isoforms, an important mechanism regulating cell-specific function. Little is known about AS in lung development, particularly in alveolar type II (ATII) cells. ErbB4 receptor isoforms Jma and Jmb have significant and opposing functions in the brain, heart, and lung development and/or disease. However, the regulators of ErbB4 AS are unknown. ErbB4 AS regulators in fetal mouse ATII cells control its function in ATII cell maturation.

METHODS

Candidate ErbB4 AS regulators were found using in silico analysis. Their developmental expression was studied in fetal mouse ATII cells. The effects of splice factor downregulation and upregulation on ATII cell maturation were analyzed.

RESULTS

ErbB4-Jma increased significantly in ATII cells after gestation E16.5. In silico analysis found four candidate splice factors: FOX2, CUG/CELF1, TIAR, and HUB. Fetal ATII cells expressed these factors in distinct developmental profiles. HUB downregulation in E17.5 ATII cells increased Jma isoform levels and Sftpb gene expression and decreased Jmb. HUB overexpression decreased Jma and Sftpb.

CONCLUSIONS

ErbB4 AS is developmentally controlled by HUB in fetal ATII cells, promoting ATII differentiation. Regulated AS expression during ATII cell differentiation suggests novel therapeutic strategies to approach human disease.

IMPACT

Alternative splicing (AS) of the ErbB4 receptor, involving mutually exclusive exon inclusion, creates Jma and Jmb isoforms with distinct differences in receptor processing and function. The Jma isoform of ErbB4 promotes differentiation of fetal lung alveolar type II cells. The AS is mediated in part by the RNA-binding protein HUB. The molecular mechanism of AS for ErbB4 has not been previously described. The regulation of ErbB4 AS has important implications in the development of organs, such as the lung, brain, and heart, and for disease, including cancer.

Cryptic splicing events result in unexpected protein products from calpain-10 (CAPN10) cDNA

Calpain-10 (CAPN10) belongs to the calpain superfamily. Genetic polymorphisms of the CAPN10 gene are associated with susceptibility to develop type 2 diabetes mellitus. Although the role of CAPN10 in the pathophysiology of diabetes has been extensively investigated, its biochemical properties are largely unknown. In this report, we made the surprising discovery that CAPN10 cDNA transcripts are subject to cryptic splicing and unexpected protein products were expressed. The same set of splicing products was reproducibly detected in four types of cultured cells including the primary culture of mouse myoblast. At least, one of the products was identical to a natural splicing variant. Sequence analysis of the splicing potential of CAPN10 cDNA, together with mutagenesis studies, resulted in the identification of a powerful splicing acceptor site at the junction of the sequences encoded by exons 9 and 10. We successfully extended the analysis to create expression construct resistant to splicing for both human and mouse CAPN10. The construct allowed us to analyze two major CAPN10 isoforms and reveal their difference in substrate proteolysis and potential cell functions. These results demonstrate that proteins produced from cDNA do not necessarily reflect the original nucleotide sequence. We provide insight into the property of recombinantly expressed CAPN10 proteins in cultured cells circumventing unexpected protein products.

Modulation of RNA Splicing by Oligonucleotides: Mechanisms of Action and Therapeutic Implications

Dysregulation of RNA splicing causes many diseases and disorders. Several therapeutic approaches have been developed to correct aberrant alternative splicing events for the treatment of cancers and hereditary diseases, including gene therapy and redirecting splicing, using small molecules or splice switching oligonucleotides (SSO). Significant advances in the chemistry and pharmacology of nucleic acid have led to the development of clinically approved SSO drugs for the treatment of spinal muscular dystrophy and Duchenne muscular dystrophy (DMD). In this review, we discuss the mechanisms of SSO action with emphasis on "less common" approaches to modulate alternative splicing, including bipartite and bifunctional SSO, oligonucleotide decoys for splice factors and SSO-mediated mRNA degradation via AS-NMD and NGD pathways. We briefly discuss the current progress and future perspectives of SSO therapy for rare and ultrarare diseases.

Staufen1 unwinds the secondary structure and facilitates the translation of fatty acid binding protein 4 mRNA during adipogenesis

Adipogenesis is regulated by genetic interactions, in which post-transcriptional regulation plays an important role. Staufen double-stranded RNA binding protein 1 (Staufen1 or STAU1) plays diverse roles in RNA processing and adipogenesis. Previously, we found that the downregulation of STAU1 affects the expression of fatty acid-binding protein 4 (FABP4) at the protein level but not at the mRNA level. This study aimed to determine the mechanism underlying the regulation of FABP4 expression by STAU1, explaining the inconsistency between FABP4 mRNA and protein levels. We used RNA interference, photoactivatable ribonucleoside enhanced cross-linking and immunoprecipitation, and an adeno-associated virus to examine the functions of STAU1 in adipogenesis. Our results indicate that STAU1 binds to the coding sequences of FABP4, thereby regulating the translation of FABP4 mRNA by unwinding the double-stranded structure. Furthermore, STAU1 mediates adipogenesis by regulating the secretion of free fatty acids. However, STAU1 knockdown decreases the fat weight/body weight ratio but does not affect the plasma triglyceride levels. These findings describe the mechanisms involved in STAU1-mediated regulation of FABP4 expression at the translational level during adipogenesis.

Modulation of alternative splicing during early infection of human primary B lymphocytes with Epstein-Barr virus (EBV): a novel function for the viral EBNA-LP protein

Epstein-Barr virus (EBV) is a human herpesvirus associated with human cancers worldwide. Ex vivo, the virus efficiently infects resting human B lymphocytes and induces their continuous proliferation. This process is accompanied by a global reprogramming of cellular gene transcription. However, very little is known on the impact of EBV infection on the regulation of alternative splicing, a pivotal mechanism that plays an essential role in cell fate determination and is often deregulated in cancer. In this study, we have developed a systematic time-resolved analysis of cellular mRNA splice variant expression during EBV infection of resting B lymphocytes. Our results reveal that major modifications of alternative splice variant expression appear as early as day 1 post-infection and suggest that splicing regulation provides—besides transcription—an additional mechanism of gene expression regulation at the onset of B cell activation and proliferation. We also report a role for the viral proteins, EBNA2 and EBNA-LP, in the modulation of specific alternative splicing events and reveal a previously unknown function for EBNA-LP—together with the RBM4 splicing factor—in the alternative splicing regulation of two important modulators of cell proliferation and apoptosis respectively, NUMB and BCL-X.

Rapid structure-function insights via hairpin-centric analysis of big RNA structure probing datasets

The functions of RNA are often tied to its structure, hence analyzing structure is of significant interest when studying cellular processes. Recently, large-scale structure probing (SP) studies have enabled assessment of global structure-function relationships via standard data summarizations or local folding. Here, we approach structure quantification from a hairpin-centric perspective where putative hairpins are identified in SP datasets and used as a means to capture local structural effects. This has the advantage of rapid processing of big (e.g. transcriptome-wide) data as RNA folding is circumvented, yet it captures more information than simple data summarizations. We reformulate a statistical learning algorithm we previously developed to significantly improve precision of hairpin detection, then introduce a novel nucleotide-wise measure, termed the hairpin-derived structure level (HDSL), which captures local structuredness by accounting for the presence of likely hairpin elements. Applying HDSL to data from recent studies recapitulates, strengthens and expands on their findings which were obtained by more comprehensive folding algorithms, yet our analyses are orders of magnitude faster. These results demonstrate that hairpin detection is a promising avenue for global and rapid structure-function analysis, furthering our understanding of RNA biology and the principal features which drive biological insights from SP data.

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

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

CircPLCE1 facilitates the malignant progression of colorectal cancer by repressing the SRSF2-dependent PLCE1 pre-RNA splicing

Studies have demonstrated that circular RNAs (circRNAs) play important roles in various types of cancer; however, the mechanisms of circRNAs located in the nucleus have rarely been explored. Here, we report a novel circular RNA circPLCE1 (hsa_circ_0019230) that facilitates the malignant progression of colorectal cancer (CRC) by repressing serine/arginine‐rich splicing factor 2 (SRSF2)‐dependent phospholipase C epsilon 1 (PLCE1) pre‐RNA splicing. Quantitative real‐time polymerase chain reaction was used to determine the expression of circPLCE1 in CRC tissues and cells. Cell Counting Kit‐8, Transwell and flow cytometric assays were used to assess the role of circPLE1 in CRC cell proliferation, migration and apoptosis, respectively. An animal study was conducted to test the role of circPLCE1 in vivo. Furthermore, catRAPID and RPISeq were used to predict the possible binding proteins of circPLCE1. RNA fractionation and RNA immunoprecipitation assays were used to confirm the RNA‐protein interaction. In this study, we found that circPLCE1 was more significantly down‐regulated in CRC tissues compared with that in adjacent normal tissues. However, circPLCE1 knockdown suppressed CRC cell proliferation, migration and invasion and increased apoptosis. Nude mouse experiments showed that ectopic expression of circPLCE1 dramatically increased tumour growth in vivo. Mechanistically, circPLCE1 directly bound to the SRSF2 protein, repressing SRSF2‐dependent PLCE1 pre‐RNA splicing, resulting in the progression of CRC. Individually mutating the binding sites of circPLCE1 abolished the inhibition of PLCE1 mRNA production. Our study revealed a novel molecular mechanism in the regulation of PLCE1 and suggested a new function of circular RNA.

A Missense Mutation rs781536408 (c.2395G>A) of TYK2 Affects Splicing and Causes Skipping of Exon18 in vivo

TYK2 variants can impact disease onset or progression. In our previous study, we identified abnormal splicing that happened near rs781536408 in the TYK2 gene. The purpose of this research was to examine the effect of the mutation on alternative splicing in vivo and in vitro. Whole exome sequencing was performed to identify the mutations followed by bidirectional Sanger sequencing. Then the minigene analysis was carried out based on HeLa and HEK293T cell lines. The results showed that rs781536408 (c.2395G>A, p.G799R) was homozygous in the patient, but heterozygous in parents. PCR amplification confirmed the abnormal splicing in the somatic cells of the patients, but not in the parents. Sanger sequencing results showed that there was a skipping of exon18 near the mutation. For minigene analysis, there was no difference between the wild-type and the mutant type in the two minigene construction strategies, indicating that mutation c.2395G>A had no effect on splicing in vitro. Combining the results of in vivo, we speculated that the effect of the mutation on splicing was not absolute, but rather in degree.

A Homozygous Deep Intronic Mutation Alters the Splicing of Nebulin Gene in a Patient With Nemaline Myopathy

Nemaline myopathy is a rare disorder affecting the muscle sarcomere. Mutations in nebulin gene (NEB) are known to be responsible for about 50% of nemaline myopathy cases. Nebulin is a giant protein which is formed integrally with the sarcomeric thin filament. This complex gene is under extensive alternative splicing giving rise to multiple isoforms. In this study, we report a 6-year-old boy presenting with general muscular weaknesses. Identification of rod-shaped structures in the patient' biopsy raised doubt about the presence of a nemaline myopathy. Next-generation sequencing was used to identify a causative mutation for the patient syndrome. A homozygous deep intronic substitution was found in the intron 144 of the NEB. The variant was predicted by in silico tools to create a new donor splice site. Molecular analysis has shown that the mutation could alter splicing events of the nebulin gene leading to a significant decrease of isoforms level. This change in the expression level of nebulin could give rise to functional consequences in the sarcomere. These results are consistent with the phenotypes observed in the patient. Such a discovery of variants in this gene will allow a better understanding of the involvement of nebulin in neuromuscular diseases and help find new treatments for the nemaline myopathy.

Identification of prognostic and bone metastasis​-related alternative splicing signatures in mesothelioma

Mesothelioma (MESO) is an infrequent tumor derived from mesothelial cells of pleura, peritoneum, pericardium, and tunica vaginalis testis. Despite advancement in technologies and better understanding of tumor progression mechanism, the prognosis of MESO remains poor. The role of alternative splicing events (ASEs) in the oncogenesis, tumor metastasis and drug resistance has been widely discussed in multiple cancers. But the prognosis and potential therapeutic value of ASEs in MESO were not clearly studied by now. We constructed a prognostic model using RNA sequencing data and matched ASE data of MESO patients obtained from the TCGA and TCGASpliceSeq database. A total of 3,993 ASEs were identified associated with overall survival using Cox regression analysis. Eight of them were finally figured out to institute the model by lasso regression analysis. The risk score of the model can predict the prognosis independently. Among the identified 390 splicing factors (SF), HSPA1A and DDX3Y was significantly associated with 43 OS-SEs. Among these OS-SEs, SNX5-58744-AT (p = 0.048) and SNX5-58745-AT (p = 0.048) were significantly associated with bone metastasis. Co-expression analysis of signal pathways and SNX5-58744-AT, SNX5-58745-AT was also depicted using GSVA. Finally, we proposed that splicing factor (SF) HSPA1A could regulate SNX5-58744-AT (R = -0.414) and SNX5-58745-AT (R = 0.414) through the pathway "Class I MHC mediated antigen processing and presentation" (R = 0.400). In this way, tumorigenesis and bone metastasis of MESO were controlled.

Characterization of Alternative Splicing (AS) Events during Chicken (Gallus gallus) Male Germ-Line Stem Cell Differentiation with Single-Cell RNA-seq

Simple Summary

Studies have shown that alternative splicing (AS) has been utilized in a wide variety of life processes. However, there are very few studies on AS during germ cell development. In this study, we preliminarily investigated the variation of variable shear events during the formation of chicken germ cells through the RNA-seq data analysis of embryonic stem cells (ESCs), gonad PGCs (gPGCs), and spermatogonia stem cells (SSCs), and the critical AS mode for several crucial stage-specific genes, which were identified during germ cell development. The results of this study lay a theoretical foundation for further analysis of the regulation mechanism of key genes involved in germ cell formation.

Abstract

Alternative splicing (AS) is a ubiquitous, co-transcriptional, and post-transcriptional regulation mechanism during certain developmental processes, such as germ cell differentiation. A thorough understanding of germ cell differentiation will help us to open new avenues for avian reproduction, stem cell biology, and advances in medicines for human consumption. Here, based on single-cell RNA-seq, we characterized genome-wide AS events in manifold chicken male germ cells: embryonic stem cells (ESCs), gonad primordial germ cells (gPGCs), and spermatogonia stem cells (SSCs). A total of 38,494 AS events from 15,338 genes were detected in ESCs, with a total of 48,955 events from 14,783 genes and 49,900 events from 15,089 genes observed in gPGCs and SSCs, respectively. Moreover, this distribution of AS events suggests the diverse splicing feature of ESCs, gPGCs, and SSCs. Finally, several crucial stage-specific genes, such as NANOG, POU5F3, LIN28B, BMP4, STRA8, and LHX9, were identified in AS events that were transmitted in ESCs, gPGCs, and SSCs. The gene expression results of the RNA-seq data were validated by qRT-PCR. In summary, we provided a comprehensive atlas of the genome-wide scale of the AS event landscape in male chicken germ-line cells and presented its distribution for the first time. This research may someday improve treatment options for men suffering from male infertility.

Synonymous variants that disrupt messenger RNA structure are significantly constrained in the human population

BACKGROUND

The role of synonymous single-nucleotide variants in human health and disease is poorly understood, yet evidence suggests that this class of "silent" genetic variation plays multiple regulatory roles in both transcription and translation. One mechanism by which synonymous codons direct and modulate the translational process is through alteration of the elaborate structure formed by single-stranded mRNA molecules. While tools to computationally predict the effect of non-synonymous variants on protein structure are plentiful, analogous tools to systematically assess how synonymous variants might disrupt mRNA structure are lacking.

RESULTS

We developed novel software using a parallel processing framework for large-scale generation of secondary RNA structures and folding statistics for the transcriptome of any species. Focusing our analysis on the human transcriptome, we calculated 5 billion RNA-folding statistics for 469 million single-nucleotide variants in 45,800 transcripts. By considering the impact of all possible synonymous variants globally, we discover that synonymous variants predicted to disrupt mRNA structure have significantly lower rates of incidence in the human population.

CONCLUSIONS

These findings support the hypothesis that synonymous variants may play a role in genetic disorders due to their effects on mRNA structure. To evaluate the potential pathogenic impact of synonymous variants, we provide RNA stability, edge distance, and diversity metrics for every nucleotide in the human transcriptome and introduce a "Structural Predictivity Index" (SPI) to quantify structural constraint operating on any synonymous variant. Because no single RNA-folding metric can capture the diversity of mechanisms by which a variant could alter secondary mRNA structure, we generated a SUmmarized RNA Folding (SURF) metric to provide a single measurement to predict the impact of secondary structure altering variants in human genetic studies.

Consequences for Pancreatic β-Cell Identity and Function of Unregulated Transcript Processing

Mounting evidence suggests a role for alternative splicing (AS) of transcripts in the normal physiology and pathophysiology of the pancreatic β-cell. In the apparent absence of RNA repair systems, RNA decay pathways are likely to play an important role in controlling the stability, distribution and diversity of transcript isoforms in these cells. Around 35% of alternatively spliced transcripts in human cells contain premature termination codons (PTCs) and are targeted for degradation via nonsense-mediated decay (NMD), a vital quality control process. Inflammatory cytokines, whose levels are increased in both type 1 (T1D) and type 2 (T2D) diabetes, stimulate alternative splicing events and the expression of NMD components, and may or may not be associated with the activation of the NMD pathway. It is, however, now possible to infer that NMD plays a crucial role in regulating transcript processing in normal and stress conditions in pancreatic β-cells. In this review, we describe the possible role of Regulated Unproductive Splicing and Translation (RUST), a molecular mechanism embracing NMD activity in relationship to AS and translation of damaged transcript isoforms in these cells. This process substantially reduces the abundance of non-functional transcript isoforms, and its dysregulation may be involved in pancreatic β-cell failure in diabetes.

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

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

A novel missense mutation of RPGR identified from retinitis pigmentosa affects splicing of the ORF15 region and causes loss of transcript heterogeneity

Mutations in the retinitis pigmentosa GTPase regulator (RPGR) gene, are the major cause of X-linked retinitis pigmentosa (RP), in which exon open reading frame 15 (ORF15) of RPGR has been implicated to play a substantial role. We identified a novel hemizygous missense mutation E585K of RPGR from whole-exome sequencing of RP. RNA-Seq analysis and functional study were conducted to investigate the underlying pathogenic mechanism of the mutation. Our results showed that the mutation actually affected RPGR ORF15 splicing. RNA-Seq analysis of the human retina followed by validation in cells revealed a complex splicing pattern near the 3' boundary of RPGR exon 14 in the ORF15 region, resulting from a variety of alternative splicing events (ASEs). The wildtype RPGR mini-gene expressed in human 293T cells confirmed these ASEs in vitro. In contrast, without new RNA species detected, the mutant mini-gene disrupted the splicing pattern of the ORF15 region, and caused loss of RPGR transcript heterogeneity. The RNA species derived from the mutant mini-gene were predominated by a minor out-of-frame transcript that was also observed in wildtype RPGR, resulting from an upstream alternative 5' splice site in exon 14. Our findings therefore provide insights into the influence of RPGR exonic mutations on alternative splicing of the ORF15 region, and the underlying molecular mechanism of RP.

Structural disruption of exonic stem-loops immediately upstream of the intron regulates mammalian splicing

Recognition of highly degenerate mammalian splice sites by the core spliceosomal machinery is regulated by several protein factors that predominantly bind exonic splicing motifs. These are postulated to be single-stranded in order to be functional, yet knowledge of secondary structural features that regulate the exposure of exonic splicing motifs across the transcriptome is not currently available. Using transcriptome-wide RNA structural information we show that retained introns in mouse are commonly flanked by a short (≲70 nucleotide), highly base-paired segment upstream and a predominantly single-stranded exonic segment downstream. Splicing assays with select pre-mRNA substrates demonstrate that loops immediately upstream of the introns contain pre-mRNA-specific splicing enhancers, the substitution or hybridization of which impedes splicing. Additionally, the exonic segments flanking the retained introns appeared to be more enriched in a previously identified set of hexameric exonic splicing enhancer (ESE) sequences compared to their spliced counterparts, suggesting that base-pairing in the exonic segments upstream of retained introns could be a means for occlusion of ESEs. The upstream exonic loops of the test substrate promoted recruitment of splicing factors and consequent pre-mRNA structural remodeling, leading up to assembly of the early spliceosome. These results suggest that disruption of exonic stem-loop structures immediately upstream (but not downstream) of the introns regulate alternative splicing events, likely through modulating accessibility of splicing factors.

Intrinsic Regulatory Role of RNA Structural Arrangement in Alternative Splicing Control

Alternative splicing is a highly sophisticated process, playing a significant role in posttranscriptional gene expression and underlying the diversity and complexity of organisms. Its regulation is multilayered, including an intrinsic role of RNA structural arrangement which undergoes time- and tissue-specific alterations. In this review, we describe the principles of RNA structural arrangement and briefly decipher its cis- and trans-acting cellular modulators which serve as crucial determinants of biological functionality of the RNA structure. Subsequently, we engage in a discussion about the RNA structure-mediated mechanisms of alternative splicing regulation. On one hand, the impairment of formation of optimal RNA structures may have critical consequences for the splicing outcome and further contribute to understanding the pathomechanism of severe disorders. On the other hand, the structural aspects of RNA became significant features taken into consideration in the endeavor of finding potential therapeutic treatments. Both aspects have been addressed by us emphasizing the importance of ongoing studies in both fields.

A combined RNA-seq and whole genome sequencing approach for identification of non-coding pathogenic variants in single families

Inherited retinal degenerations (IRDs) are at the focus of current genetic therapeutic advancements. For a genetic treatment such as gene therapy to be successful, an accurate genetic diagnostic is required. Genetic diagnostics relies on the assessment of the probability that a given DNA variant is pathogenic. Non-coding variants present a unique challenge for such assessments as compared to coding variants. For one, non-coding variants are present at much higher number in the genome than coding variants. In addition, our understanding of the rules that govern the non-coding regions of the genome is less complete than our understanding of the coding regions. Methods that allow for both the identification of candidate non-coding pathogenic variants and their functional validation may help overcome these caveats allowing for a greater number of patients to benefit from advancements in genetic therapeutics. We present here an unbiased approach combining whole genome sequencing (WGS) with patient-induced pluripotent stem cell (iPSC)-derived retinal organoids (ROs) transcriptome analysis. With this approach, we identified and functionally validated a novel pathogenic non-coding variant in a small family with a previously unresolved genetic diagnosis.

Interactions of the GABRG2 polymorphisms and childhood trauma on suicide attempt and related traits in depressed patients

BACKGROUND

Previously, we reported that the longest variant of the GABA A receptor γ2 subunit (GABRG2) was associated with suicidal behavior. The present study therefore aimed to determine whether polymorphisms near the alternatively spliced exon of GABRG2 are associated with suicide attempt (SA) and its related traits, and how these variants might interact with reported childhood trauma (CT) in their association with suicidal behavior.

METHODS

We examined 5 single nucleotide polymorphisms (SNPs) of GABRG2. Subjects were suicide Attempters (N = 94), non-suicide attempters (N = 168) with MDD or Bipolar depression, and healthy volunteers (N = 100). Data on demographics, depression severity and suicide attempts were collected. Participants also completed a set of instruments assessing CT, and lifetime aggression and impulsivity.. GABRG2 polymorphisms were genotyped using Sanger sequencing.

RESULTS

Allele A of rs211034 was a protective factor for SA (OR = 0.50 (0.32, 0.80), p = 0.003), and had an interaction effect with emotional neglect (OR = 0.89 (0.82, 0.97), p = 0.006) on depression. One haploblock (consisting of rs211035 and rs211034) was identified within these SNPs, and subjects with haplotype GA (frequency = 7.3%), had lower rate of SA (OR=0.26(0.10, 0.67), p = 0.006). Cognitive impulsivity (OR=1.38)1.24,1.55), p < 0.001), non-planning impulsivity (OR = 1.18 (1.10,1.25), p < 0.001), anger (OR = 1.13 (1.07,1.19), p < 0.001), impulsivity total score (OR = 1.10(1.06,1.15), p < 0.001), hostility (OR = 1.10 (1.04, 1.15), p < 0.001), aggression total score (OR = 1.05 (1.03,1.07), p < 0.001) were associated with depression, meanwhile, hopelessness (OR = 2.18 (1.56, 3.04), p < 0.001) and impulsivity total score (OR = 1.05 (1.02,1.08), p < 0.001) were associated with the risk of SA, adjusted by age and gender. There was no mediation effect in the relationship among CT, gene polymorphisms and SA or depression through increased impulsivity or aggression.

LIMITATIONS

The main limitation of this study is its modest sample size. More genetic variants as well as epigenetic markers should be examined in future studies.

CONCLUSIONS

These findings add to evidence for the involvement of GABRG2 and impulsivity and hopelessness in SA independent from their association with depression. More research is needed on possible mediators of the relationship between GABA-related gene and SA.

mRNA Metabolism in Cardiac Development and Disease: Life After Transcription

The central dogma of molecular biology illustrates the importance of mRNAs as critical mediators between genetic information encoded at the DNA level and proteomes/metabolomes that determine the diverse functional outcome at the cellular and organ levels. Although the total number of protein-producing (coding) genes in the mammalian genome is ~20,000, it is evident that the intricate processes of cardiac development and the highly regulated physiological regulation in the normal heart, as well as the complex manifestation of pathological remodeling in a diseased heart, would require a much higher degree of complexity at the transcriptome level and beyond. Indeed, in addition to an extensive regulatory scheme implemented at the level of transcription, the complexity of transcript processing following transcription is dramatically increased. RNA processing includes post-transcriptional modification, alternative splicing, editing and transportation, ribosomal loading, and degradation. While transcriptional control of cardiac genes has been a major focus of investigation in recent decades, a great deal of progress has recently been made in our understanding of how post-transcriptional regulation of mRNA contributes to transcriptome complexity. In this review, we highlight some of the key molecular processes and major players in RNA maturation and post-transcriptional regulation. In addition, we provide an update to the recent progress made in the discovery of RNA processing regulators implicated in cardiac development and disease. While post-transcriptional modulation is a complex and challenging problem to study, recent technological advancements are paving the way for a new era of exciting discoveries and potential clinical translation in the context of cardiac biology and heart disease.

(-)-Epicatechin metabolites promote vascular health through epigenetic reprogramming of endothelial-immune cell signaling and reversing systemic low-grade inflammation

Ingestion of (-)-epicatechin flavanols reverses endothelial dysfunction by increasing flow mediated dilation and by reducing vascular inflammation and oxidative stress, monocyte-endothelial cell adhesion and transendothelial monocyte migration in vitro and in vivo. This involves multiple changes in gene expression and epigenetic DNA methylation by poorly understood mechanisms. By in silico docking and molecular modeling we demonstrate favorable binding of different glucuronidated, sulfated or methylated (-)-epicatechin metabolites to different DNA methyltransferases (DNMT1/DNMT3A). In favor of this model, genome-wide DNA methylation profiling of endothelial cells treated with TNF and different (-)-epicatechin metabolites revealed specific DNA methylation changes in gene networks controlling cell adhesion-extravasation endothelial hyperpermeability as well as gamma-aminobutyric acid, renin-angiotensin and nitric oxide hypertension pathways. Remarkably, blood epigenetic profiles of an 8 weeks intervention with monomeric and oligomeric flavanols (MOF) including (-)-epicatechin in male smokers revealed individual epigenetic gene changes targeting similar pathways as the in vitro exposure experiments in endothelial cells. Furthermore, epigenetic changes following MOF diet intervention oppose atherosclerosis associated epigenetic changes. In line with biological data, the individual epigenetic response to a MOF diet is associated with different vascular health parameters (glutathione peroxidase 1 and endothelin-1 expression, acetylcholine-mediated microvascular response), in part involving systemic shifts in blood immune cell types which reduce the neutrophil-lymphocyte ratio (NLR). Altogether, our study suggests that different (-)-epicatechin metabolites promote vascular health in part via epigenetic reprogramming of endothelial-immune cell signaling and reversing systemic low-grade inflammation.

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.

Identification of novel prognostic alternative splicing signature in papillary renal cell carcinoma

Papillary renal cell carcinoma (pRCC) is a heterogeneous disease containing multifocal or solitary tumors with an aggressive phenotype. Increasing evidence has indicated the involvement of aberrant splicing variants in renal cell cancer, while systematic profiling of aberrant alternative splicing (AS) in pRCC was lacking and largely unknown. In the current study, comprehensive profiling of AS events were performed based on the integration of pRCC cohort from the Cancer Genome Atlas database and SpliceSeq software. With rigorous screening and univariate Cox analysis, a total of 2077 prognoses AS events from 1642 parent genes were identified. Then, stepwise least absolute shrinkage and selection operator method-penalized Cox regression analyses with 10-fold cross-validation followed by multivariate Cox regression were used to construct the prognostic AS signatures within each AS type. And a final 21 AS event-based signature was proposed which showed potent prognostic capability in stratifying patients into low- and high-risk subgroups (P < .0001). Furthermore, time-dependent receiver operating characteristics curves confirmed that the final AS signature was effective and robust in predicting overall survival for pRCC patients with the area under the curve above 0.9 from 1 to 5 years. In addition, splicing correlation network was built to uncover the potential regulatory pattern among prognostic splicing factors and candidate AS events. Besides, gene set enrichment analysis revealed the involvement of these candidates AS events in tumor-related pathways including extracellular matrix organization, oxidative phosphorylation, and P53 signaling pathways. Taken together, our results could contribute to elucidating the underlying mechanism of AS in the oncogenesis process and broaden the novel field of prognostic and clinical application of molecule-targeted approaches in pRCC.

Molecular Mechanisms and Determinants of Innovative Correction Approaches in Coagulation Factor Deficiencies

Molecular strategies tailored to promote/correct the expression and/or processing of defective coagulation factors would represent innovative therapeutic approaches beyond standard substitutive therapy. Here, we focus on the molecular mechanisms and determinants underlying innovative approaches acting at DNA, mRNA and protein levels in inherited coagulation factor deficiencies, and in particular on: (i) gene editing approaches, which have permitted intervention at the DNA level through the specific recognition, cleavage, repair/correction or activation of target sequences, even in mutated gene contexts; (ii) the rescue of altered pre-mRNA processing through the engineering of key spliceosome components able to promote correct exon recognition and, in turn, the synthesis and secretion of functional factors, as well as the effects on the splicing of missense changes affecting exonic splicing elements; this section includes antisense oligonucleotide- or siRNA-mediated approaches to down-regulate target genes; (iii) the rescue of protein synthesis/function through the induction of ribosome readthrough targeting nonsense variants or the correction of folding defects caused by amino acid substitutions. Overall, these approaches have shown the ability to rescue the expression and/or function of potentially therapeutic levels of coagulation factors in different disease models, thus supporting further studies in the future aimed at evaluating the clinical translatability of these new strategies.

Differential isoform expression and alternative splicing in sex determination in mice

Background

Alternative splicing (AS) may play an important role in gonadal sex determination (GSD) in mammals. The present study was designed to identify differentially expressed isoforms and AS modifications accompanying GSD in mice.

Results

Using deep RNA-sequencing, we performed a transcriptional analysis of XX and XY gonads during sex determination on embryonic days 11 (E11) and 12 (E12). Analysis of differentially expressed genes (DEG) identified hundreds of genes related to GSD and early sex differentiation that may represent good candidates for sex reversal. Expression at time point E11 in males was significantly enriched in RNA splicing and mRNA processing Gene Ontology terms. Differentially expressed isoform analysis identified hundreds of specific isoforms related to GSD, many of which showed no differences in the DEG analysis. Hundreds of AS events were identified as modified at E11 and E12. Female E11 gonads featured sex-biased upregulation of intron retention (in genes related to regulation of transcription, protein phosphorylation, protein transport and mRNA splicing) and exon skipping (in genes related to chromatin repression) suggesting AS as a post-transcription mechanism that controls sex determination of the bipotential fetal gonad.

Conclusion

Our data suggests an important role of splicing regulatory mechanisms for sex determination in mice.

Electronic supplementary material

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

Recurrence and Familial Inheritance of Intronic NIPBL Pathogenic Variant Associated With Mild CdLS

Splicing pathogenic variants account for a notable fraction of NIPBL alterations underlying Cornelia de Lange syndrome but are likely underrepresented, due to overlooking of non-canonical intronic variants by traditional and contemporary sequencing methods. We describe five subjects, belonging to three families, displaying a mild Cornelia de Lange syndrome phenotype who carry the NIPBL pathogenic variant c.5329–15A>G, affecting the IVS27 branch site, yet reported in a single case. By RNA analysis we evidenced two alternative transcripts: the exon 28 in frame skipped transcript, described in the published case and an out-of-frame transcript retaining 14 nucleotides of IVS27 3′end. Even if both aberrant transcripts are at negligible levels, their presence justifies the CdLS phenotype shared by our patients consisting of borderline-mild cognitive impairment and slight but typical facial dysmorphisms. Transmission of the pathogenic variant from pauci-symptomatic mother to her siblings emphasizes the need of molecular diagnosis extended to deep intronic regions in patients with subtle but recognizable CdLS phenotype.

Systematic evaluation of isoform function in literature reports of alternative splicing

BACKGROUND

Although most genes in mammalian genomes have multiple isoforms, an ongoing debate is whether these isoforms are all functional as well as the extent to which they increase the functional repertoire of the genome. To ground this debate in data, it would be helpful to have a corpus of experimentally-verified cases of genes which have functionally distinct splice isoforms (FDSIs).

RESULTS

We established a curation framework for evaluating experimental evidence of FDSIs, and analyzed over 700 human and mouse genes, strongly biased towards genes that are prominent in the alternative splicing literature. Despite this bias, we found experimental evidence meeting the classical definition for functionally distinct isoforms for ~ 5% of the curated genes. If we relax our criteria for inclusion to include weaker forms of evidence, the fraction of genes with evidence of FDSIs remains low (~ 13%). We provide evidence that this picture will not change substantially with further curation and conclude there is a large gap between the presumed impact of splicing on gene function and the experimental evidence. Furthermore, many functionally distinct isoforms were not traceable to a specific isoform in Ensembl, a database that forms the basis for much computational research.

CONCLUSIONS

We conclude that the claim that alternative splicing vastly increases the functional repertoire of the genome is an extrapolation from a limited number of empirically supported cases. We also conclude that more work is needed to integrate experimental evidence and genome annotation databases. Our work should help shape research around the role of splicing on gene function from presuming large general effects to acknowledging the need for stronger experimental evidence.

Detection of Splicing Abnormalities and Genotype-Phenotype Correlation in X-linked Alport Syndrome

BACKGROUND

X-linked Alport syndrome (XLAS) is a progressive hereditary nephropathy caused by mutations in the COL4A5 gene. Genotype-phenotype correlation in male XLAS is relatively well established; relative to truncating mutations, nontruncating mutations exhibit milder phenotypes. However, transcript comparison between XLAS cases with splicing abnormalities that result in a premature stop codon and those with nontruncating splicing abnormalities has not been reported, mainly because transcript analysis is not routinely conducted in patients with XLAS.

METHODS

We examined transcript expression for all patients with suspected splicing abnormalities who were treated at one hospital between January of 2006 and July of 2017. Additionally, we recruited 46 males from 29 families with splicing abnormalities to examine genotype-phenotype correlation in patients with truncating (n=21, from 14 families) and nontruncating (n=25, from 15 families) mutations at the transcript level.

RESULTS

We detected 41 XLAS families with abnormal splicing patterns and described novel XLAS atypical splicing patterns (n=14) other than exon skipping caused by point mutations in the splice consensus sequence. The median age for developing ESRD was 20 years (95% confidence interval, 14 to 23 years) among patients with truncating mutations and 29 years (95% confidence interval, 25 to 40 years) among patients with nontruncating mutations (P=0.001).

CONCLUSIONS

We report unpredictable atypical splicing in the COL4A5 gene in male patients with XLAS and reveal that renal prognosis differs significantly for patients with truncating versus nontruncating splicing abnormalities. Our results suggest that splicing modulation should be explored as a therapy for XLAS with truncating mutations.

Alternative Splicing of Differentiated Myeloid Cell Transcripts after Infection by Anaplasma phagocytophilum Impacts a Selective Group of Cellular Programs

Eukaryotic proteome diversity exceeds that encoded within individual genes, and results in part from alternative splicing events of pre-messenger RNA. The diversity of these splicing events can shape the outcome in development and differentiation of normal tissues, and is important in pathogenic circumstances such as cancer and some heritable conditions. A role for alternative splicing of eukaryotic genes in response to viral and intracellular bacterial infections has only recently been recognized, and plays an important role in providing fitness for microbial survival, while potentially enhancing pathogenicity. Anaplasma phagocytophilum survives within mammalian neutrophils by reshaping transcriptional programs that govern cellular functions. We applied next generation RNAseq to ATRA-differentiated HL-60 cells established to possess transcriptional and functional responses similar to A. phagocytophilum-infected human neutrophils. This demonstrated an increase in transcripts with infection and high proportion of alternatively spliced transcript events (ASEs) for which predicted gene ontology processes were in part distinct from those identified by evaluation of single transcripts or gene-level analyses alone. The alternative isoforms are not on average shorter, and no alternative splicing in genes encoding spliceosome components is noted. Although not evident at gene-level analyses, individual spliceosome transcripts that impact nearly all spliceosome components were significantly upregulated. How the distinct GO processes predicted by ASEs are regulated by infection and whether they are relevant to fitness or pathogenicity of A. phagocytophilum should be addressed in more detailed studies.

Haplotype-specific MAPT exon 3 expression regulated by common intronic polymorphisms associated with Parkinsonian disorders

BACKGROUND

Genome wide association studies have identified microtubule associated protein tau (MAPT) H1 haplotype single nucleotide polymorphisms (SNPs) as leading common risk variants for Parkinson's disease, progressive supranuclear palsy and corticobasal degeneration. The MAPT risk variants fall within a large 1.8 Mb region of high linkage disequilibrium, making it difficult to discern the functionally important risk variants. Here, we leverage the strong haplotype-specific expression of MAPT exon 3 to investigate the functionality of SNPs that fall within this H1 haplotype region of linkage disequilibrium.

METHODS

In this study, we dissect the molecular mechanisms by which haplotype-specific SNPs confer allele-specific effects on the alternative splicing of MAPT exon 3. Firstly, we use haplotype-hybrid whole-locus genomic MAPT vectors studies to identify functional SNPs. Next, we characterise the RNA-protein interactions at two loci by mass spectrometry. Lastly, we knockdown candidate splice factors to determine their effect on MAPT exon 3 using a novel allele-specific qPCR assay.

RESULTS

Using whole-locus genomic DNA expression vectors to express MAPT haplotype variants, we demonstrate that rs17651213 regulates exon 3 inclusion in a haplotype-specific manner. We further investigated the functionality of this region using RNA-electrophoretic mobility shift assays to show differential RNA-protein complex formation at the H1 and H2 sequence variants of SNP rs17651213 and rs1800547 and subsequently identified candidate trans-acting splicing factors interacting with these functional SNPs sequences by RNA-protein pull-down experiment and mass spectrometry. Finally, gene knockdown of candidate splice factors identified by mass spectrometry demonstrate a role for hnRNP F and hnRNP Q in the haplotype-specific regulation of exon 3 inclusion.

CONCLUSIONS

We identified common splice factors hnRNP F and hnRNP Q regulating the haplotype-specific splicing of MAPT exon 3 through intronic variants rs1800547 and rs17651213. This work demonstrates an integrated approach to characterise the functionality of risk variants in large regions of linkage disequilibrium.

Small molecule modulation of splicing factor expression is associated with rescue from cellular senescence

Background

Altered expression of mRNA splicing factors occurs with ageing in vivo and is thought to be an ageing mechanism. The accumulation of senescent cells also occurs in vivo with advancing age and causes much degenerative age-related pathology. However, the relationship between these two processes is opaque. Accordingly we developed a novel panel of small molecules based on resveratrol, previously suggested to alter mRNA splicing, to determine whether altered splicing factor expression had potential to influence features of replicative senescence.

Results

Treatment with resveralogues was associated with altered splicing factor expression and rescue of multiple features of senescence. This rescue was independent of cell cycle traverse and also independent of SIRT1, SASP modulation or senolysis. Under growth permissive conditions, cells demonstrating restored splicing factor expression also demonstrated increased telomere length, re-entered cell cycle and resumed proliferation. These phenomena were also influenced by ERK antagonists and agonists.

Conclusions

This is the first demonstration that moderation of splicing factor levels is associated with reversal of cellular senescence in human primary fibroblasts. Small molecule modulators of such targets may therefore represent promising novel anti-degenerative therapies.

Electronic supplementary material

The online version of this article (10.1186/s12860-017-0147-7) contains supplementary material, which is available to authorized users.

A Genetic Screen for Pre-mRNA Splicing Mutants of Arabidopsis thaliana Identifies Putative U1 snRNP Components RBM25 and PRP39a

In a genetic screen for mutants showing modified splicing of an alternatively spliced GFP reporter gene in Arabidopsis thaliana, we identified mutations in genes encoding the putative U1 small nuclear ribonucleoprotein (snRNP) factors RBM25 and PRP39a. The latter has not yet been studied for its role in pre-messenger RNA (pre-mRNA) splicing in plants. Both proteins contain predicted RNA-binding domains and have been implicated in 5′ splice site selection in yeast and metazoan cells. In rbm25 mutants, splicing efficiency of GFP pre-mRNA was reduced and GFP protein levels lowered relative to wild-type plants. By contrast, prp39a mutants exhibited preferential splicing of a U2-type AT-AC intron in GFP pre-mRNA and elevated levels of GFP protein. These opposing findings indicate that impaired function of either RBM25 or PRP39a can differentially affect the same pre-mRNA substrate. Given a prior genome-wide analysis of alternative splicing in rbm25 mutants, we focused on examining the alternative splicing landscape in prp39a mutants. RNA-seq experiments performed using two independent prp39a alleles revealed hundreds of common genes undergoing changes in alternative splicing, including PRP39a itself, a second putative U1 snRNP component PRP40b, and genes encoding a number of general transcription-related proteins. The prp39a mutants displayed somewhat delayed flowering, shorter stature, and reduced seed set but no other obvious common defects under normal conditions. Mutations in PRP39b, the paralog of PRP39a, did not visibly alter GFP expression, indicating the paralogs are not functionally equivalent in this system. Our study provides new information on the contribution of PRP39a to alternative splicing and expands knowledge of plant splicing factors.