Distribution of SR protein exonic splicing enhancer motifs in human protein-coding genes

Interferon-Regulated Expression of Cellular Splicing Factors Modulates Multiple Levels of HIV-1 Gene Expression and Replication

Type I interferons (IFN-Is) are pivotal in innate immunity against human immunodeficiency virus I (HIV-1) by eliciting the expression of IFN-stimulated genes (ISGs), which encompass potent host restriction factors. While ISGs restrict the viral replication within the host cell by targeting various stages of the viral life cycle, the lesser-known IFN-repressed genes (IRepGs), including RNA-binding proteins (RBPs), affect the viral replication by altering the expression of the host dependency factors that are essential for efficient HIV-1 gene expression. Both the host restriction and dependency factors determine the viral replication efficiency; however, the understanding of the IRepGs implicated in HIV-1 infection remains greatly limited at present. This review provides a comprehensive overview of the current understanding regarding the impact of the RNA-binding protein families, specifically the two families of splicing-associated proteins SRSF and hnRNP, on HIV-1 gene expression and viral replication. Since the recent findings show specifically that SRSF1 and hnRNP A0 are regulated by IFN-I in various cell lines and primary cells, including intestinal lamina propria mononuclear cells (LPMCs) and peripheral blood mononuclear cells (PBMCs), we particularly discuss their role in the context of the innate immunity affecting HIV-1 replication.

Development and validation of AI/ML derived splice-switching oligonucleotides

Splice-switching oligonucleotides (SSOs) are antisense compounds that act directly on pre-mRNA to modulate alternative splicing (AS). This study demonstrates the value that artificial intelligence/machine learning (AI/ML) provides for the identification of functional, verifiable, and therapeutic SSOs. We trained XGboost tree models using splicing factor (SF) pre-mRNA binding profiles and spliceosome assembly information to identify modulatory SSO binding sites on pre-mRNA. Using Shapley and out-of-bag analyses we also predicted the identity of specific SFs whose binding to pre-mRNA is blocked by SSOs. This step adds considerable transparency to AI/ML-driven drug discovery and informs biological insights useful in further validation steps. We applied this approach to previously established functional SSOs to retrospectively identify the SFs likely to regulate those events. We then took a prospective validation approach using a novel target in triple negative breast cancer (TNBC), NEDD4L exon 13 (NEDD4Le13). Targeting NEDD4Le13 with an AI/ML-designed SSO decreased the proliferative and migratory behavior of TNBC cells via downregulation of the TGFβ pathway. Overall, this study illustrates the ability of AI/ML to extract actionable insights from RNA-seq data.

Regulation of micro- and small-exon retention and other splicing processes by GRP20 for flower development

Pre-mRNA splicing is crucial for gene expression and depends on the spliceosome and splicing factors. Plant exons have an average size of ~180 nucleotides and typically contain motifs for interactions with spliceosome and splicing factors. Micro exons (<51 nucleotides) are found widely in eukaryotes and in genes for plant development and environmental responses. However, little is known about transcript-specific regulation of splicing in plants and about the regulators for micro exon splicing. Here we report that glycine-rich protein 20 (GRP20) is an RNA-binding protein and required for splicing of ~2,100 genes including those functioning in flower development and/or environmental responses. Specifically, GRP20 is required for micro-exon retention in transcripts of floral homeotic genes; these micro exons are conserved across angiosperms. GRP20 is also important for small-exon (51-100 nucleotides) splicing. In addition, GRP20 is required for flower development. Furthermore, GRP20 binds to poly-purine motifs in micro and small exons and a spliceosome component; both RNA binding and spliceosome interaction are important for flower development and micro-exon retention. Our results provide new insights into the mechanisms of micro-exon retention in flower development.

Splicing defects in rare diseases: transcriptomics and machine learning strategies towards genetic diagnosis

Genomic variants affecting pre-messenger RNA splicing and its regulation are known to underlie many rare genetic diseases. However, common workflows for genetic diagnosis and clinical variant interpretation frequently overlook splice-altering variants. To better serve patient populations and advance biomedical knowledge, it has become increasingly important to develop and refine approaches for detecting and interpreting pathogenic splicing variants. In this review, we will summarize a few recent developments and challenges in using RNA sequencing technologies for rare disease investigation. Moreover, we will discuss how recent computational splicing prediction tools have emerged as complementary approaches for revealing disease-causing variants underlying splicing defects. We speculate that continuous improvements to sequencing technologies and predictive modeling will not only expand our understanding of splicing regulation but also bring us closer to filling the diagnostic gap for rare disease patients.

Insights into established and emerging roles of SR protein family in plants and animals

Splicing of pre-mRNA is an essential part of eukaryotic gene expression. Serine-/arginine-rich (SR) proteins are highly conserved RNA-binding proteins present in all metazoans and plants. SR proteins are involved in constitutive and alternative splicing, thereby regulating the transcriptome and proteome diversity in the organism. In addition to their role in splicing, SR proteins are also involved in mRNA export, nonsense-mediated mRNA decay, mRNA stability, and translation. Due to their pivotal roles in mRNA metabolism, SR proteins play essential roles in normal growth and development. Hence, any misregulation of this set of proteins causes developmental defects in both plants and animals. SR proteins from the animal kingdom are extensively studied for their canonical and noncanonical functions. Compared with the animal kingdom, plant genomes harbor more SR protein-encoding genes and greater diversity of SR proteins, which are probably evolved for plant-specific functions. Evidence from both plants and animals confirms the essential role of SR proteins as regulators of gene expression influencing cellular processes, developmental stages, and disease conditions. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Processing > Splicing Regulation/Alternative Splicing.

Deciphering hERG Mutation in Long QT Syndrome Type 2 Using Antisense Oligonucleotide-Mediated Techniques: Lessons from Cystic Fibrosis

Long QT syndrome type 2 (LQT2) is a genetic disorder caused by mutations in the KCNH2 gene, also known as the human ether-a-go-go-related gene (hERG). Over 30% of hERG mutations result in a premature termination codon (PTC) that triggers a process called nonsense-mediated mRNA decay (NMD), where the mRNA transcript is degraded. NMD is a quality control mechanism that removes faulty mRNA to prevent the translation of truncated proteins. Recent advances in antisense oligonucleotide (ASO) technology in the field of cystic fibrosis (CF) have yielded significant progress, including the ASO-mediated comprehensive characterization of key NMD factors and exon-skipping therapy. These advances have contributed to our understanding of the role of PTC-containing mutations in disease phenotypes and have also led to the development of potentially useful therapeutic strategies. Historically, studies of CF have provided valuable insights for the research on LQT2, particularly concerning increasing the expression of hERG. In this article, we outline the current state of knowledge regarding ASO, NMD, and hERG and discuss the introduction of ASO technology in the CF to elucidate the pathogenic mechanisms through targeting NMD. We also discuss the potential clinical therapeutic benefits and limitations of ASO for the management of LQT2. By drawing on lessons learned from CF research, we explore the potential translational values of these advances into LQT2 studies.

Investigation of the regulatory effects of synthesized antisense oligonucleotides on androgen receptor (AR) exon 3 splicing in prostate cancer cells

The Androgen Receptor (AR) gene plays a key role in castration-resistant prostate cancer (CRPC). Controlling the progression of CRPC by inhibiting AR gene expression is one of the core directions for prostate cancer (Pca) drug development. A 23-amino acids retention, named exon 3a, into the DNA binding domain of the splice variant AR23 has been shown to prevent AR from entering the nucleus and restore the sensitivity of cancer cells to related therapies. In this study, we conducted a preliminary investigation of the splicing modulation of the AR gene in order to develop a splice-switching therapy for Pca by promoting exon 3a inclusion. Using mutagenesis-coupled RT-PCR with AR minigene and over-expression of certain splicing factors, we found that serine/arginine-rich (SR) proteins are key factors facilitating the recognition of the 3' splice site of exon 3a (L-3' SS), while the deletion or blocking of the polypyrimidine tract (PPT) region of the original 3' splice site of exon 3 (S-3' SS) could strongly enhance exon 3a splicing without affecting the function of any SR protein. Furthermore, we designed a series of antisense oligonucleotides (ASOs) to screen drug candidates, and ASOs targeting S-3' SS and its PPT region or the exonic region of exon 3 turned out to be most effective in rescuing exon 3a splicing. A dose-response test indicated ASO12 as the lead candidate drug significantly promoting the inclusion of exon 3a to more than 85%. MTT assay confirmed that the cell proliferation was significantly inhibited after ASO treatment. Our results provide the first glance to AR splicing regulation. With several promising therapeutic ASO candidates obtained here, further development of ASO drugs to treat CRPC is strongly encouraged.

A translational repression reporter assay for the analysis of RNA-binding protein consensus sites

RNA-binding proteins are essential regulators of RNA processing and function. Translational repression assays can be used to study how they interact with specific RNA sequences by insertion of such a consensus sequence into the 5' untranslated region of a reporter mRNA and measuring reporter protein translation. The straightforward set-up of these translational repression assays avoids the need for the isolation of the protein or the RNA providing speed, robustness and a low-cost method. Here, we report the optimization of the assay to function with linear RNA sequences instead of the previously reported hairpin type sequences to allow the study of a wider variety of RNA-binding proteins. Multiplication of a consensus sequence strongly improves the signal allowing analysis by both fluorescence intensity measurements and flow cytometry.

Alternative splicing of HSPA12A pre-RNA by SRSF11 contributes to metastasis potential of colorectal cancer

BACKGROUND

Dysregulation of alternative splicing (AS) induced by serine/arginine-rich proteins has recently been linked to cancer metastasis. Nonetheless, as a member of the serine/arginine-rich protein family, the involvement of SRSF11 in colorectal cancer (CRC) is unknown.

METHODS

The TCGA dataset and clinical samples were used to assess SRSF11 expression levels in CRC. For SRSF11, functional experiments were conducted both in vitro and in vivo. RNA-seq technology was used to analyze and screen SRSF11-triggered AS events, which were then confirmed by in vivo UV crosslinking and immunoprecipitation (CLIP) and mini-gene reporter assays. Jalview software was used to determine the preferential binding motif with relation to exon skipping (ES) events. Furthermore, coimmunoprecipitation (Co-IP) and Phospho-tag SDS-PAGE experiments were used to investigate PAK5-mediated phosphorylation regulation on SRSF11, and in vitro kinase experiments validated the interaction.

RESULTS

In CRC, SRSF11 was discovered to be overexpressed and associated with a poor prognosis. And SRSF11 played a pro-metastatic role in vitro and in vivo. By screening SRSF11-regulated AS events, we identified the binding motif of SRSF11-triggered splicing-switching of HSPA12A AS, which specifically regulated HSPA12A AS by directly binding to a motif in exon 2. Mechanistically, the HSPA12A transcript with exon 2 retention increased N-cadherin expression by promoting RNA stability. Furthermore, the oncogenic kinase PAK5 phosphorylated SRSF11 at serine 287, protecting it from ubiquitination degradation.

CONCLUSIONS

SRSF11 exerts pro-metastatic effects in CRC by inhibiting the AS of HSPA12A pre-RNA. Our findings point to SRSF11-regulated HSPA12A splicing as a novel relationship between SRSF11-regulated splicing and CRC metastasis and suggest a PAK5/SRSF11/HSPA12A axis as a potential therapeutic target and prognostic biomarker in CRC.

Code inside the codon: The role of synonymous mutations in regulating splicing machinery and its impact on disease

In eukaryotes, precise pre-mRNA processing, including alternative splicing, is essential to carry out the intricate protein translation process. Both point mutations (that alter the translated protein sequence) and synonymous mutations (that do not alter the translated protein sequence) are capable of affecting the splicing process. Synonymous mutations are known to affect gene expression via altering mRNA stability, mRNA secondary structure, splicing processes, and translational kinetics. In higher eukaryotes, precise splicing is regulated by three weakly conserved cis-elements, 5' and 3' splice sites and the branch site. Many other cis-acting elements (exonic/intronic splicing enhancers and silencers) and trans-acting splicing factors (serine and arginine-rich proteins and heterogeneous nuclear ribonucleoproteins) have also been found to enhance or suppress the splicing process. The appearance of synonymous mutations in cis-acting elements can alter the splicing process by changing the binding pattern of splicing factors to exonic splicing enhancers or silencer motifs. This results in exon skipping, intron retention, and various other forms of alternative splicing, eventually leading to the emergence of a wide range of diseases. The focus of this review is to elucidate the role of synonymous mutations and their impact on abnormal splicing mechanisms. Further, this study highlights the function of synonymous mutation in mediating abnormal splicing in cancer and development of X-linked, and autosomal inherited diseases.

Pseudoexon activation in disease by non-splice site deep intronic sequence variation - wild type pseudoexons constitute high-risk sites in the human genome

Accuracy of pre-messenger RNA (pre-mRNA) splicing is crucial for normal gene expression. Complex regulation supports the spliceosomal distinction between authentic exons and the many seemingly functional splice sites delimiting pseudoexons. Pseudoexons are nonfunctional intronic sequences that can be activated for aberrant inclusion in mRNA, which may cause disease. Pseudoexon activation is very challenging to predict, in particular when activation occurs by sequence variants that alter the splicing regulatory environment without directly affecting splice sites. As pseudoexon inclusion often evades detection due to activation of nonsense-mediated mRNA decay, and because conventional diagnostic procedures miss deep intronic sequence variation, pseudoexon activation is a heavily underreported disease mechanism. Pseudoexon characteristics have mainly been studied based on in silico predicted sequences. Moreover, because recognition of sequence variants that create or strengthen splice sites is possible by comparison with well-established consensus sequences, this type of pseudoexon activation is by far the most frequently reported. Here we review all known human disease-associated pseudoexons that carry functional splice sites and are activated by deep intronic sequence variants located outside splice site sequences. We delineate common characteristics that make this type of wild type pseudoexons distinct high-risk sites in the human genome.

Alternative Splicing Regulation of Low-Frequency Genetic Variants in Exon 2 of TREM2 in Alzheimer's Disease by Splicing-Based Aggregation

TREM2 is among the most well-known Alzheimer’s disease (AD) risk genes; however, the functional roles of its AD-associated variants remain to be elucidated, and most known risk alleles are low-frequency variants whose investigation is challenging. Here, we utilized a splicing-guided aggregation method in which multiple low-frequency TREM2 variants were bundled together to investigate the functional impact of those variants on alternative splicing in AD. We analyzed whole genome sequencing (WGS) and RNA-seq data generated from cognitively normal elderly controls (CN) and AD patients in two independent cohorts, representing three regions in the frontal lobe of the human brain: the dorsolateral prefrontal cortex (CN = 213 and AD = 376), frontal pole (CN = 72 and AD = 175), and inferior frontal (CN = 63 and AD = 157). We observed an exon skipping event in the second exon of TREM2, with that exon tending to be more frequently skipped (p = 0.0012) in individuals having at least one low-frequency variant that caused loss-of-function for a splicing regulatory element. In addition, genes differentially expressed between AD patients with high vs. low skipping of the second exon (i.e., loss of a TREM2 functional domain) were significantly enriched in immune-related pathways. Our splicing-guided aggregation method thus provides new insight into the regulation of alternative splicing of the second exon of TREM2 by low-frequency variants and could be a useful tool for further exploring the potential molecular mechanisms of multiple, disease-associated, low-frequency variants.

Discovery of a pre-mRNA structural scaffold as a contributor to the mammalian splicing code

The specific recognition of splice signals at or near exon-intron junctions is not explained by their weak conservation and instead is postulated to require a multitude of features embedded in the pre-mRNA strand. We explored the possibility of 3D structural scaffold of AdML-a model pre-mRNA substrate-guiding early spliceosomal components to the splice signal sequences. We find that mutations in the non-cognate splice signal sequences impede recruitment of early spliceosomal components due to disruption of the global structure of the pre-mRNA. We further find that the pre-mRNA segments potentially interacting with the early spliceosomal component U1 snRNP are distributed across the intron, that there is a spatial proximity of 5' and 3' splice sites within the pre-mRNA scaffold, and that an interplay exists between the structural scaffold and splicing regulatory elements in recruiting early spliceosomal components. These results suggest that early spliceosomal components can recognize a 3D structural scaffold beyond the short splice signal sequences, and that in our model pre-mRNA, this scaffold is formed across the intron involving the major splice signals. This provides a conceptual basis to analyze the contribution of recognizable 3D structural scaffolds to the splicing code across the mammalian transcriptome.

A novel mutation in intron 11 donor splice site, responsible of a rare genotype in thyroglobulin gene by altering the pre-mRNA splincing process. Cell expression and bioinformatic analysis

Thyroglobulin (TG) is a homodimeric glycoprotein synthesized by the thyroid gland. To date, two hundred twenty-seven variations of the TG gene have been identified in humans. Thyroid dyshormonogenesis due to TG gene mutations have an estimated incidence of approximately 1 in 100,000 newborns. The clinical spectrum ranges from euthyroid to mild or severe hypothyroidism. The purpose of the present study was to identify and characterize new variants in the TG gene. We report an Argentine patient with congenital hypothyroidism, enlarged thyroid gland and low levels of serum TG. Sequencing of DNA, expression of chimeric minigenes as well as bioinformatics analysis were performed. DNA sequencing identified the presence of compound heterozygous mutations in the TG gene: the maternal mutation consists of a c.3001+5G > A, whereas the paternal mutation consists of p.Arg296*. Minigen analysis of the variant c.3001+5A performed in HeLa, CV1 and Hek293T cell lines, showed a total lack of transcript expression. So, in order to validate that the loss of expression was caused by such variation, site-directed mutagenesis was performed on the mutated clone, which previously had a pSPL3 vector change, to give rise to a wild-type clone c.3001+5G, endorsing that the mutation c.3001+5G > A is the cause of the total lack of expression. In conclusion, we demonstrate that the c.3001+5G > A mutation causes a rare genotype, altering the splicing of the pre-mRNA. This work contributes to elucidating the molecular bases of TG defects associated with congenital hypothyroidism and expands our knowledge in relation to the pathologic roles of the position 5 in the donor splice site.

Structure of SRSF1 RRM1 bound to RNA reveals an unexpected bimodal mode of interaction and explains its involvement in SMN1 exon7 splicing

The human prototypical SR protein SRSF1 is an oncoprotein that contains two RRMs and plays a pivotal role in RNA metabolism. We determined the structure of the RRM1 bound to RNA and found that the domain binds preferentially to a CN motif (N is for any nucleotide). Based on this solution structure, we engineered a protein containing a single glutamate to asparagine mutation (E87N), which gains the ability to bind to uridines and thereby activates SMN exon7 inclusion, a strategy that is used to cure spinal muscular atrophy. Finally, we revealed that the flexible inter-RRM linker of SRSF1 allows RRM1 to bind RNA on both sides of RRM2 binding site. Besides revealing an unexpected bimodal mode of interaction of SRSF1 with RNA, which will be of interest to design new therapeutic strategies, this study brings a new perspective on the mode of action of SRSF1 in cells.

SRSF1 is an oncoprotein that plays important roles in RNA metabolism. We reveal the structure of the human SRSF1 RRM1 bound to RNA, and propose a bimodal mode of interaction of the protein with RNA. A single mutation in RRM1 changed SRSF1 specificity for RNA and made it active on SMN2 exon7 splicing.

Splicing analysis of SLC40A1 missense variations and contribution to hemochromatosis type 4 phenotypes

Hemochromatosis type 4, or ferroportin disease, is considered as the second leading cause of primary iron overload after HFE-related hemochromatosis. The disease, which is predominantly associated with missense variations in the SLC40A1 gene, is characterized by wide clinical heterogeneity. We tested the possibility that some of the reported missense mutations, despite their positions within exons, cause splicing defects. Fifty-eight genetic variants were selected from the literature based on two criteria: a precise description of the nucleotide change and individual evidence of iron overload. The selected variants were investigated by different in silico prediction tools and prioritized for midigene splicing assays. Of the 15 variations tested in vitro, only two were associated with splicing changes. We confirm that the c.1402G>A transition (p.Gly468Ser) disrupts the exon 7 donor site, leading to the use of an exonic cryptic splicing site and the generation of a truncated reading frame. We observed, for the first time, that the p.Gly468Ser substitution has no effect on the ferroportin iron export function. We demonstrate alternative splicing of exon 5 in different cell lines and show that the c.430A>G (p.Asn144Asp) variant promotes exon 5 inclusion. This could be part of a gain-of-function mechanism. We conclude that splicing mutations rarely contribute to hemochromatosis type 4 phenotypes. An in-depth investigation of exon 5 auxiliary splicing sequences may help to elucidate the mechanism by which splicing regulatory proteins regulate the production of the full length SLC40A1 transcript and to clarify its physiological importance.

The RNA binding protein SRSF1 is a master switch of gene expression and regulation in the immune system

Serine/Arginine splicing factor 1 (SRSF1) is an RNA binding protein abundantly expressed in most tissues. The pleiotropic functions of SRSF1 exert multiple roles in gene expression by regulating major steps in transcription, processing, export through the nuclear pores and translation of nascent RNA transcripts. The aim of this review is to highlight recent findings in the functions of this protein and to describe its role in immune system development, functions and regulation.

Mutations primarily alter the inclusion of alternatively spliced exons

Genetic analyses and systematic mutagenesis have revealed that synonymous, non-synonymous and intronic mutations frequently alter the inclusion levels of alternatively spliced exons, consistent with the concept that altered splicing might be a common mechanism by which mutations cause disease. However, most exons expressed in any cell are highly-included in mature mRNAs. Here, by performing deep mutagenesis of highly-included exons and by analysing the association between genome sequence variation and exon inclusion across the transcriptome, we report that mutations only very rarely alter the inclusion of highly-included exons. This is true for both exonic and intronic mutations as well as for perturbations in trans. Therefore, mutations that affect splicing are not evenly distributed across primary transcripts but are focussed in and around alternatively spliced exons with intermediate inclusion levels. These results provide a resource for prioritising synonymous and other variants as disease-causing mutations.

Detection of polymorphism within leptin gene in Egyptian river buffalo and predict its effects on different molecular levels

BACKGROUND

Leptin (LEP) regulates the glucose homeostasis directly and centrally by the regulation of the insulin levels or indirectly by alternation of the levels of the other glucose metabolism regulator hormones. The present investigation studied the polymorphism in LEP gene which is related to fertility in 81 female Egyptian river buffalo.

RESULTS

The PCR-RFLP pattern of the gene using the restriction enzyme Eco91I showed that all the animals had monomorphic pattern in the studied gene which consists of CC. A 511-bp fragment from LEP gene was amplified and sequenced. The homology between the amplified LEP gene fragment in buffalo and cattle, sheep, goat, human, and mouse on the nucleotides sequence level was 99, 97, 97, 87, and 79%, respectively, and on the translated amino acids sequence level was 100, 98, 98, 85, and 82%, respectively. Several SNPs were detected; among them, the T27C SNP disrupted an intronic splicing silencer. The A114G, A310G, G263A, and G379A SNPs disrupt exonic splicing enhancers, and the last two SNPs create new exonic splicing enhancers. The A114G, C163A, A211G, G288A, A310G, A322G, G330C, C348T, T360C, and G379A SNPs cause S71G, T87 N, N103S, E129K, E136G, Y140C, E143Q, R149W, S153P, and R159Q amino acids mutations. N103S, E129K, E136G, Y140C, E143Q, and S153P were classified as deleterious mutations. Y140, E143, N103, and R149 were the most conserved among the mutated amino acids. S71G only increased the stability of the leptin protein while the remaining mutations decreased it.

CONCLUSION

Four SNPs were revealed among the tested animals. Twenty-one SNPs were found between the sequenced amplicon and the buffalo records in the Genbank. Some SNPs were predicted to have several effects on different biological processes like mRNA splicing, protein stability, and the gene functions.

miR-1/206 downregulates splicing factor Srsf9 to promote C2C12 differentiation

Background

Myogenesis is driven by specific changes in the transcriptome that occur during the different stages of muscle differentiation. In addition to controlled transcriptional transitions, several other post-transcriptional mechanisms direct muscle differentiation. Both alternative splicing and miRNA activity regulate gene expression and production of specialized protein isoforms. Importantly, disruption of either process often results in severe phenotypes as reported for several muscle diseases. Thus, broadening our understanding of the post-transcriptional pathways that operate in muscles will lay the foundation for future therapeutic interventions.

Methods

We employed bioinformatics analysis in concert with the well-established C2C12 cell system for predicting and validating novel miR-1 and miR-206 targets engaged in muscle differentiation. We used reporter gene assays to test direct miRNA targeting and studied C2C12 cells stably expressing one of the cDNA candidates fused to a heterologous, miRNA-resistant 3′ UTR. We monitored effects on differentiation by measuring fusion index, myotube area, and myogenic gene expression during time course differentiation experiments.

Results

Gene ontology analysis revealed a strongly enriched set of putative miR-1 and miR-206 targets associated with RNA metabolism. Notably, the expression levels of several candidates decreased during C2C12 differentiation. We discovered that the splicing factor Srsf9 is a direct target of both miRNAs during myogenesis. Persistent Srsf9 expression during differentiation impaired myotube formation and blunted induction of the early pro-differentiation factor myogenin as well as the late differentiation marker sarcomeric myosin, Myh8.

Conclusions

Our data uncover novel miR-1 and miR-206 cellular targets and establish a functional link between the splicing factor Srsf9 and myoblast differentiation. The finding that miRNA-mediated clearance of Srsf9 is a key myogenic event illustrates the coordinated and sophisticated interplay between the diverse components of the gene regulatory network.

Transcriptome analysis reveals significant differences between primary plasma cell leukemia and multiple myeloma even when sharing a similar genetic background

Primary plasma cell leukemia (pPCL) is a highly aggressive plasma cell dyscrasia characterised by short remissions and very poor survival. Although the 17p deletion is associated with poor outcome and extramedullary disease in MM, its presence does not confer the degree of aggressiveness observed in pPCL. The comprehensive exploration of isoform expression and RNA splicing events may provide novel information about biological differences between the two diseases. Transcriptomic studies were carried out in nine newly diagnosed pPCL and ten MM samples, all of which harbored the 17p deletion. Unsupervised cluster analysis clearly distinguished pPCL from MM samples. In total 3584 genes and 20033 isoforms were found to be deregulated between pPCL and MM. There were 2727 significantly deregulated isoforms of non-differentially expressed genes. Strangely enough, significant differences were observed in the expression of spliceosomal machinery components between pPCL and MM, in respect of the gene, isoform and the alternative splicing events expression. In summary, transcriptome analysis revealed significant differences in the relative abundance of isoforms between pPCL and MM, even when they both had the 17p deletion. The mRNA processing pathway including RNA splicing machinery emerged as one of the most remarkable mechanisms underlying the biological differences between the two entities.

Association of Kiss1 and GPR54 Gene Polymorphisms with Polycystic Ovary Syndrome among Sri Lankan Women

Polycystic ovary syndrome (PCOS) is the commonest endocrine disorder affecting women of reproductive age. Its aetiology, though yet unclear, is presumed to have an oligogenic basis interacting with environmental factors. Kisspeptins are peptide products of Kiss1 gene that control the hypothalamic pituitary (HPG) axis by acting via G protein-coupled receptor known as GPR54. There is paucity of data on the role of Kiss1 and GPR54 gene in PCOS. We aimed to identify the polymorphisms in Kiss1 and GPR54 genes and explore their association with serum kisspeptin levels among Sri Lankan women with well-characterized PCOS. Consecutive women with PCOS manifesting from adolescence (n=55) and adult controls (n=110) were recruited. Serum kisspeptin and testosterone levels were determined by ELISA method. Whole gene sequencing was performed to identify the polymorphisms in Kiss1 and GPR54 genes. Serum kisspeptin and testosterone concentrations were significantly higher in women with PCOS than controls: kisspeptin 4.873nmol/L versus 4.127nmol/L; testosterone 4.713nmol/L versus 3.415 nmol/L, p<0.05. Sequencing the GPR54 gene revealed 5 single nucleotide polymorphisms (SNPs), rs10407968, rs1250729403, rs350131, chr19:918686, and chr19:918735, with two novel SNPs (chr19:918686 and chr19:918735), while sequencing the Kiss1 gene revealed 2 SNPs, rs5780218 and rs4889. All identified SNPs showed no significant difference in frequency between patients and controls. GPR54 gene rs350131 polymorphism (G/T) was detected more frequently in our study population. The heterozygous allele (AG) of GPR54 gene novel polymorphism chr19:918686 showed a marginal association with serum kisspeptin levels (p=0.053). Genetic variations in GPR54 and Kiss1 genes are unlikely to be associated with PCOS among Sri Lankan women manifesting from adolescence. Meanwhile the heterozygous allele of chr19:918686 is probably associated with serum kisspeptin concentrations, which suggests a potential role in the aetiology of PCOS.

Quantitative Activity Profile and Context Dependence of All Human 5' Splice Sites

Pre-mRNA splicing is an essential step in the expression of most human genes. Mutations at the 5' splice site (5'ss) frequently cause defective splicing and disease due to interference with the initial recognition of the exon-intron boundary by U1 small nuclear ribonucleoprotein (snRNP), a component of the spliceosome. Here, we use a massively parallel splicing assay (MPSA) in human cells to quantify the activity of all 32,768 unique 5'ss sequences (NNN/GYNNNN) in three different gene contexts. Our results reveal that although splicing efficiency is mostly governed by the 5'ss sequence, there are substantial differences in this efficiency across gene contexts. Among other uses, these MPSA measurements facilitate the prediction of 5'ss sequence variants that are likely to cause aberrant splicing. This approach provides a framework to assess potential pathogenic variants in the human genome and streamline the development of splicing-corrective therapies.

Splicing heterogeneity: separating signal from noise

Single-cell analyses have revealed a tremendous variety among cells in the abundance and chemical composition of RNA. Much of this heterogeneity is due to alternative splicing by the spliceosome. Little is known about how many of the resulting isoforms are biologically functional or just provide noise with little to no impact. The dynamic nature of the spliceosome provides numerous opportunities for regulation but is also the source of stochastic fluctuations. We discuss possible origins of splicing stochasticity, the experimental approaches for studying heterogeneity in isoforms, and the potential biological significance of noisy splicing in development and disease.

Importance of Novel Sequence Alterations in the FHIT Gene on Formation of Breast Cancer

Aims and background

The character, role and impact of FHIT gene alterations, for which recent studies have shown that the gene has a role in the early stage of carcinogenesis in breast cancer, are still unclear. Thus, the current study evaluated FHIT gene mutations from breast tissue of women with malignant and benign breast disease and to elucidate the frequency and type of mutations in this gene.

Patients and methods

Mutations in exons 5–9 of the FHIT gene were screened using the intronic primer pairs in 83 breast (67 malignant and 16 benign) tissue samples by single-strand conformational polymorphism and sequencing analysis.

Results

FHIT mutations were detected in 13 of the 67 malignant cases (19.4%) and 2 of the 16 benign cases (12.5%). Four different sequence variants were determined: two novel frame shift mutations (codon 90 insA, codon 146 delT), one intronic novel mutation (IVS8 −17 insA), and one previously identified silent transition type alteration (codon 88 C to T). In addition, determination of this silent alteration caused formation of new exonic splicing enhancer (ESE) motifs on mutated sequences by using the ESEfinder program.

Conclusions

Our data contribute significantly to that currently known about the presence of FHIT gene mutations on the formation of breast cancer.

Advances in Transcriptomics of Plants

The current global population of 7.3 billion is estimated to reach 9.7 billion in the year 2050. Rapid population growth is driving up global food demand. Additionally, global climate change, environmental degradation, drought, emerging diseases, and salty soils are the current threats to global food security. In order to mitigate the adverse effects of these diverse agricultural productivity constraints and enhance crop yield and stress-tolerance in plants, we need to go beyond traditional and molecular plant breeding. The powerful new tools for genome editing, Transcription Activator-Like Effector Nucleases (TALENs) and Clustered Regulatory Interspaced Short Palindromic Repeats (CRISPR)/Cas systems (CRISPR-Cas9), have been hailed as a quantum leap forward in the development of stress-resistant plants. Plant breeding techniques, however, have several drawbacks. Hence, identification of transcriptional regulatory elements and deciphering mechanisms underlying transcriptional regulation are crucial to avoiding unintended consequences in modified crop plants, which could ultimately have negative impacts on human health. RNA splicing as an essential regulated post-transcriptional process, alternative polyadenylation as an RNA-processing mechanism, along with non-coding RNAs (microRNAs, small interfering RNAs and long non-coding RNAs) have been identified as major players in gene regulation. In this chapter, we highlight new findings on the essential roles of alternative splicing and alternative polyadenylation in plant development and response to biotic and abiotic stresses. We also discuss biogenesis and the functions of microRNAs (miRNAs) and small interfering RNAs (siRNAs) in plants and recent advances in our knowledge of the roles of miRNAs and siRNAs in plant stress response. Graphical Abstract.

RNA splicing in human disease and in the clinic

Defects at the level of the pre-mRNA splicing process represent a major cause of human disease. Approximately 15-50% of all human disease mutations have been shown to alter functioning of basic and auxiliary splicing elements. These elements are required to ensure proper processing of pre-mRNA splicing molecules, with their disruption leading to misprocessing of the pre-mRNA molecule and disease. The splicing process is a complex process, with much still to be uncovered before we are able to accurately predict whether a reported genomic sequence variant (GV) represents a splicing-associated disease mutation or a harmless polymorphism. Furthermore, even when a mutation is correctly identified as affecting the splicing process, there still remains the difficulty of providing an exact evaluation of the potential impact on disease onset, severity and duration. In this review, we provide a brief overview of splicing diagnostic methodologies, from in silico bioinformatics approaches to wet lab in vitro and in vivo systems to evaluate splicing efficiencies. In particular, we provide an overview of how the latest developments in high-throughput sequencing can be applied to the clinic, and are already changing clinical approaches.

Evolutionarily conserved exon definition interactions with U11 snRNP mediate alternative splicing regulation on U11-48K and U11/U12-65K genes

Many splicing regulators bind to their own pre-mRNAs to induce alternative splicing that leads to formation of unstable mRNA isoforms. This provides an autoregulatory feedback mechanism that regulates the cellular homeostasis of these factors. We have described such an autoregulatory mechanism for two core protein components, U11-48K and U11/U12-65K, of the U12-dependent spliceosome. This regulatory system uses an atypical splicing enhancer element termed USSE (U11 snRNP-binding splicing enhancer), which contains two U12-type consensus 5' splice sites (5'ss). Evolutionary analysis of the USSE element from a large number of animal and plant species indicate that USSE sequence must be located 25-50 nt downstream from the target 3' splice site (3'ss). Together with functional evidence showing a loss of USSE activity when this distance is reduced and a requirement for RS-domain of U11-35K protein for 3'ss activation, our data suggests that U11 snRNP bound to USSE uses exon definition interactions for regulating alternative splicing. However, unlike standard exon definition where the 5'ss bound by U1 or U11 will be subsequently activated for splicing, the USSE element functions similarly as an exonic splicing enhancer and is involved only in upstream splice site activation but does not function as a splicing donor. Additionally, our evolutionary and functional data suggests that the function of the 5'ss duplication within the USSE elements is to allow binding of two U11/U12 di-snRNPs that stabilize each others' binding through putative mutual interactions.

Viral interactions with components of the splicing machinery

Eukaryotic genes are often interrupted by stretches of sequence with no protein coding potential or obvious function. After transcription, these interrupting sequences must be removed to give rise to the mature messenger RNA. This fundamental process is called RNA splicing and is achieved by complicated machinery made of protein and RNA that assembles around the RNA to be edited. Viruses also use RNA splicing to maximize their coding potential and economize on genetic space, and use clever strategies to manipulate the splicing machinery to their advantage. This article gives an overview of the splicing process and provides examples of viral strategies that make use of various components of the splicing system to promote their replicative cycle. Representative virus families have been selected to illustrate the interaction with various regulatory proteins and ribonucleoproteins. The unifying theme is fine regulation through protein-protein and protein-RNA interactions with the spliceosome components and associated factors to promote or prevent spliceosome assembly on given splice sites, in addition to a strong influence from cis-regulatory sequences on viral transcripts. Because there is an intimate coupling of splicing with the processes that direct mRNA biogenesis, a description of how these viruses couple the regulation of splicing with the retention or stability of mRNAs is also included. It seems that a unique balance of suppression and activation of splicing and nuclear export works optimally for each family of viruses.

Novel Intronic RNA Structures Contribute to Maintenance of Phenotype in Saccharomyces cerevisiae

The Saccharomyces cerevisiae genome has undergone extensive intron loss during its evolutionary history. It has been suggested that the few remaining introns (in only 5% of protein-coding genes) are retained because of their impact on function under stress conditions. Here, we explore the possibility that novel noncoding RNA structures (ncRNAs) are embedded within intronic sequences and are contributing to phenotype and intron retention in yeast. We employed de novo RNA structure prediction tools to screen intronic sequences in S. cerevisiae and 36 other fungi. We identified and validated 19 new intronic RNAs via RNA sequencing (RNA-seq) and RT-PCR. Contrary to the common belief that excised introns are rapidly degraded, we found that, in six cases, the excised introns were maintained intact in the cells. In another two cases we showed that the ncRNAs were further processed from their introns. RNA-seq analysis confirmed that introns in ribosomal protein genes are more highly expressed when they contain predicted RNA structures. We deleted the novel intronic RNA structure within the GLC7 intron and showed that this region, rather than the intron itself, is responsible for the cell’s ability to respond to salt stress. We also showed a direct association between the in cis presence of the intronic RNA and GLC7 expression. Overall, these data support the notion that some introns may have been maintained in the genome because they harbor functional RNA structures.

Compound heterozygous DUOX2 gene mutations (c.2335-1G>C/c.3264_3267delCAGC) associated with congenital hypothyroidism. Characterization of complex cryptic splice sites by minigene analysis

Iodide Organification defects (IOD) represent 10% of cases of congenital hypothyroidism (CH) being the main genes affected that of TPO (thyroid peroxidase) and DUOX2 (dual oxidasa 2). From a patient with clinical and biochemical criteria suggestive with CH associated with IOD, TPO and DUOX2 genes were analyzed by means of PCR-Single Strand Conformation Polymorphism analysis and sequencing. A novel heterozygous compound to the mutations c.2335-1G>C (paternal mutation, intron 17) and c.3264_3267delCAGC (maternal mutation, exon 24) was identified in the DUOX2 gene. Ex-vivo splicing assays and subsequent RT-PCR and sequencing analyses were performed on mRNA isolated from the HeLa cells transfected with wild-type and mutant pSPL3 expression vectors. The wild-type and c.2335-1G>C mutant alleles result in the complete inclusion or exclusion of exon 18, or in the activation of an exonic cryptic 5' ss with the consequent deletion of 169 bp at the end of this exon. However, we observed only a band of the expected size in normal thyroid tissue by RT-PCR. Additionally, the c.2335-1G>C mutation activates an unusual cryptic donor splice site in intron 17, located at position -14 of the authentic intron 17/exon 18 junction site, with an insertion of the last 14 nucleotides of the intron 17 in mutant transcripts with complete and partial inclusion of exon 18. The theoretical consequences of splice site mutation, predicted with the bioinformatics NNSplice, Fsplice, SPL, SPLM and MaxEntScan programs were investigated and evaluated in relation with the experimental evidence. These analyses confirm that c.2335-1G>C mutant allele would result in the abolition of the authentic splice acceptor site. The results suggest the coexistence in our patient of four putative truncated proteins of 786, 805, 806 and 1105 amino acids, with conservation of peroxidase-like domain and loss of gp91(phox)/NOX2-like domain. In conclusion a novel heterozygous compound was identified being responsible of IOD. Cryptic splicing sites have been characterized in DUOX2 gene for the first time. The use of molecular biology techniques is a valuable tool for understanding the molecular pathophysiology of this type of thyroid defects.

Diverse regulation of 3' splice site usage

The regulation of splice site (SS) usage is important for alternative pre-mRNA splicing and thus proper expression of protein isoforms in cells; its disruption causes diseases. In recent years, an increasing number of novel regulatory elements have been found within or nearby the 3'SS in mammalian genes. The diverse elements recruit a repertoire of trans-acting factors or form secondary structures to regulate 3'SS usage, mostly at the early steps of spliceosome assembly. Their mechanisms of action mainly include: (1) competition between the factors for RNA elements, (2) steric hindrance between the factors, (3) direct interaction between the factors, (4) competition between two splice sites, or (5) local RNA secondary structures or longer range loops, according to the mode of protein/RNA interactions. Beyond the 3'SS, chromatin remodeling/transcription, posttranslational modifications of trans-acting factors and upstream signaling provide further layers of regulation. Evolutionarily, some of the 3'SS elements seem to have emerged in mammalian ancestors. Moreover, other possibilities of regulation such as that by non-coding RNA remain to be explored. It is thus likely that there are more diverse elements/factors and mechanisms that influence the choice of an intron end. The diverse regulation likely contributes to a more complex but refined transcriptome and proteome in mammals.

A Quantitative Profiling Tool for Diverse Genomic Data Types Reveals Potential Associations between Chromatin and Pre-mRNA Processing

High-throughput sequencing, and genome-based datasets in general, are often represented as profiles centered at reference points to study the association of protein binding and other signals to particular regulatory mechanisms. Although these profiles often provide compelling evidence of these associations, they do not provide a quantitative assessment of the enrichment, which makes the comparison between signals and conditions difficult. In addition, a number of biases can confound profiles, but are rarely accounted for in the tools currently available. We present a novel computational method, ProfileSeq, for the quantitative assessment of biological profiles to provide an exact, nonparametric test that specific regions of the test profile have higher or lower signal densities than a control set. The method is applicable to high-throughput sequencing data (ChIP-Seq, GRO-Seq, CLIP-Seq, etc.) and to genome-based datasets (motifs, etc.). We validate ProfileSeq by recovering and providing a quantitative assessment of several results reported before in the literature using independent datasets. We show that input signal and mappability have confounding effects on the profile results, but that normalizing the signal by input reads can eliminate these biases while preserving the biological signal. Moreover, we apply ProfileSeq to ChIP-Seq data for transcription factors, as well as for motif and CLIP-Seq data for splicing factors. In all examples considered, the profiles were robust to biases in mappability of sequencing reads. Furthermore, analyses performed with ProfileSeq reveal a number of putative relationships between transcription factor binding to DNA and splicing factor binding to pre-mRNA, adding to the growing body of evidence relating chromatin and pre-mRNA processing. ProfileSeq provides a robust way to quantify genome-wide coordinate-based signal. Software and documentation are freely available for academic use at https://bitbucket.org/regulatorygenomicsupf/profileseq/.

Splicing factor SRSF1 negatively regulates alternative splicing of MDM2 under damage

Genotoxic stress induces alternative splicing of the oncogene MDM2 generating MDM2-ALT1, an isoform attributed with tumorigenic properties. However, the mechanisms underlying this event remain unclear. Here we explore MDM2 splicing regulation by utilizing a novel minigene that mimics endogenous MDM2 splicing in response to UV and cisplatinum-induced DNA damage. We report that exon 11 is necessary and sufficient for the damage-specific alternative splicing of the MDM2 minigene and that the splicing factor SRSF1 binds exon 11 at evolutionarily conserved sites. Interestingly, mutations disrupting this interaction proved sufficient to abolish the stress-induced alternative splicing of the MDM2 minigene. Furthermore, SRSF1 overexpression promoted exclusion of exon 11, while its siRNA-mediated knockdown prevented the stress-induced alternative splicing of endogenous MDM2. Additionally, we observed elevated SRSF1 levels under stress and in tumors correlating with the expression of MDM2-ALT1. Notably, we demonstrate that MDM2-ALT1 splicing can be blocked by targeting SRSF1 sites on exon 11 using antisense oligonucleotides. These results present conclusive evidence supporting a negative role for SRSF1 in MDM2 alternative splicing. Importantly, we define for the first time, a clear-cut mechanism for the regulation of damage-induced MDM2 splicing and present potential strategies for manipulating MDM2 expression via splicing modulation.

Novel compound heterozygous Thyroglobulin mutations c.745+1G>A/c.7036+2T>A associated with congenital goiter and hypothyroidism in a Vietnamese family. Identification of a new cryptic 5' splice site in the exon 6

Several patients were identified with dyshormonogenesis caused by mutations in the thyroglobulin (TG) gene. These defects are inherited in an autosomal recessive manner and affected individuals are either homozygous or compound heterozygous for the mutations. The aim of the present study was to identify new TG mutations in a patient of Vietnamese origin affected by congenital hypothyroidism, goiter and low levels of serum TG. DNA sequencing identified the presence of compound heterozygous mutations in the TG gene: the maternal mutation consists of a novel c.745+1G>A (g.IVS6 + 1G>A), whereas the hypothetical paternal mutation consists of a novel c.7036+2T>A (g.IVS40 + 2T>A). The father was not available for segregation analysis. Ex-vivo splicing assays and subsequent RT-PCR analyses were performed on mRNA isolated from the eukaryotic-cells transfected with normal and mutant expression vectors. Minigene analysis of the c.745+1G>A mutant showed that the exon 6 is skipped during pre-mRNA splicing or partially included by use of a cryptic 5' splice site located to 55 nucleotides upstream of the authentic exon 6/intron 6 junction site. The functional analysis of c.7036+2T>A mutation showed a complete skipping of exon 40. The theoretical consequences of splice site mutations, predicted with the bioinformatics tool NNSplice, Fsplice, SPL, SPLM and MaxEntScan programs were investigated and evaluated in relation with the experimental evidence. These analyses predicted that both mutant alleles would result in the abolition of the authentic splice donor sites. The c.745+1G>A mutation originates two putative truncated proteins of 200 and 1142 amino acids, whereas c.7036+2T>A mutation results in a putative truncated protein of 2277 amino acids. In conclusion, we show that the c.745+1G>A mutation promotes the activation of a new cryptic donor splice site in the exon 6 of the TG gene. The functional consequences of these mutations could be structural changes in the protein molecule that alter the biosynthesis of thyroid hormones.

Purifying selection on splice-related motifs, not expression level nor RNA folding, explains nearly all constraint on human lincRNAs

There are two strong and equally important predictors of rates of human protein evolution: The amount the gene is expressed and the proportion of exonic sequence devoted to control splicing, mediated largely by selection on exonic splice enhancer (ESE) motifs. Is the same true for noncoding RNAs, known to be under very weak purifying selection? Prior evidence suggests that selection at splice sites in long intergenic noncoding RNAs (lincRNAs) is important. We now report multiple lines of evidence indicating that the great majority of purifying selection operating on lincRNAs in humans is splice related. Splice-related parameters explain much of the between-gene variation in evolutionary rate in humans. Expression rate is not a relevant predictor, although expression breadth is weakly so. In contrast to protein-coding RNAs, we observe no relationship between evolutionary rate and lincRNA stability. As in protein-coding genes, ESEs are especially abundant near splice junctions and evolve slower than non-ESE sequence equidistant from boundaries. Nearly all constraint in lincRNAs is at exon ends (N.B. the same is not witnessed in Drosophila). Although we cannot definitely answer the question as to why splice-related selection is so important, we find no evidence that splicing might enable the nonsense-mediated decay pathway to capture transcripts incorrectly processed by ribosomes. We find evidence consistent with the notion that splicing modifies the underlying chromatin through recruitment of splice-coupled chromatin modifiers, such as CHD1, which in turn might modulate neighbor gene activity. We conclude that most selection on human lincRNAs is splice mediated and suggest that the possibility of splice-chromatin coupling is worthy of further scrutiny.

FTO-dependent demethylation of N6-methyladenosine regulates mRNA splicing and is required for adipogenesis

The role of Fat Mass and Obesity-associated protein (FTO) and its substrate N6-methyladenosine (m⁶A) in mRNA processing and adipogenesis remains largely unknown. We show that FTO expression and m⁶A levels are inversely correlated during adipogenesis. FTO depletion blocks differentiation and only catalytically active FTO restores adipogenesis. Transcriptome analyses in combination with m⁶A-seq revealed that gene expression and mRNA splicing of grouped genes are regulated by FTO. M⁶A is enriched in exonic regions flanking 5′- and 3′-splice sites, spatially overlapping with mRNA splicing regulatory serine/arginine-rich (SR) protein exonic splicing enhancer binding regions. Enhanced levels of m⁶A in response to FTO depletion promotes the RNA binding ability of SRSF2 protein, leading to increased inclusion of target exons. FTO controls exonic splicing of adipogenic regulatory factor RUNX1T1 by regulating m⁶A levels around splice sites and thereby modulates differentiation. These findings provide compelling evidence that FTO-dependent m⁶A demethylation functions as a novel regulatory mechanism of RNA processing and plays a critical role in the regulation of adipogenesis.

Genome-wide association study identifies phospholipase C zeta 1 (PLCz1) as a stallion fertility locus in Hanoverian warmblood horses

A consistently high level of stallion fertility plays an economically important role in modern horse breeding. We performed a genome-wide association study for estimated breeding values of the paternal component of the pregnancy rate per estrus cycle (EBV-PAT) in Hanoverian stallions. A total of 228 Hanoverian stallions were genotyped using the Equine SNP50 Beadchip. The most significant association was found on horse chromosome 6 for a single nucleotide polymorphism (SNP) within phospholipase C zeta 1 (PLCz1). In the close neighbourhood to PLCz1 is located CAPZA3 (capping protein (actin filament) muscle Z-line, alpha 3). The gene PLCz1 encodes a protein essential for spermatogenesis and oocyte activation through sperm induced Ca²⁺-oscillation during fertilization. We derived equine gene models for PLCz1 and CAPZA3 based on cDNA and genomic DNA sequences. The equine PLCz1 had four different transcripts of which two contained a premature termination codon. Sequencing all exons and their flanking sequences using genomic DNA samples from 19 Hanoverian stallions revealed 47 polymorphisms within PLCz1 and one SNP within CAPZA3. Validation of these 48 polymorphisms in 237 Hanoverian stallions identified three intronic SNPs within PLCz1 as significantly associated with EBV-PAT. Bioinformatic analysis suggested regulatory effects for these SNPs via transcription factor binding sites or microRNAs. In conclusion, non-coding polymorphisms within PLCz1 were identified as conferring stallion fertility and PLCz1 as candidate locus for male fertility in Hanoverian warmblood. CAPZA3 could be eliminated as candidate gene for fertility in Hanoverian stallions.

Elevated SRPK1 lessens apoptosis in breast cancer cells through RBM4-regulated splicing events

Imbalanced splicing of premessenger RNA is typical of tumorous malignancies, and the regulatory mechanisms involved in several tumorigenesis-associated splicing events are identified. Elevated expression of serine-arginine protein kinase 1 (SRPK1) may participate in the pathway responsible for the dysregulation of splicing events in malignant tumor cells. In this study, we observed a correlation between the cytoplasmic accumulation of RNA-binding motif protein 4 (RBM4) and up-regulated SRPK1 in breast cancer cells. The production of the IR-B and MCL-1S transcripts was induced separately by the overexpression of RBM4 and SRPK1 gene silencing. Overexpressed RBM4 simultaneously bound to the CU-rich elements within the MCL-1 exon2 and the downstream intron, which subsequently facilitated the exclusion of the regulated exon. Breast cancer cells are deprived of apoptotic resistance through the RBM4-mediated up-regulation of the IR-B and MCL-1S transcripts. These findings suggest that the splicing events regulated by the SRPK1-RMB4 network may contribute to tumorigenesis through altered sensitivity to apoptotic signals in breast cancer cells.

Slight variations in the SC35 ESE sequence motif among human chromosomes: a computational approach

Gene expression is initiated by the binding of transcription factors to cis-regulatory modules such as enhancer elements binding to the Serine/Arginine proteins. Recently, we noticed an increased ability to identify the location as well as the motifs of enhancers using genome-wide information on spliceosomal factor occupancy, cofactor recruitment and chromatin modifications. In this study, we have undertaken a large-scale genomic analysis in an attempt to uncover if the exonic splicing enhancer motif binding to the SC35 and the SRp40 SR proteins is conserved among several groups of human genes. For the SRp40, the results showed that the ESE consensus is conserved among human genes. Concerning the SC35 SR protein, results showed an ESE motif conserved among human tissues and between different levels of muscular cell differentiation and within the same chromosome. However, this motif displays subtle discrepancies between genes localized in different chromosomes. These results emphasize the presence of different translational isoforms of the SFRS2 gene encoding for the SC35, or different post-translational protein maturations in different chromosomes, confirming that chromatin structure is another layer of gene regulation. These links between chromatin pattern and splicing give further mechanistic support to functional interconnections between splicing, transcription and chromatin structure, and raise the intriguing possibility of the existence of a memory for splicing patterns to be inherited through epigenetic modifications.

Emerging functions of SRSF1, splicing factor and oncoprotein, in RNA metabolism and cancer

Serine/Arginine Splicing Factor 1 (SRSF1) is the archetype member of the SR protein family of splicing regulators. Since its discovery over two decades ago, SRSF1 has been repeatedly surprising and intriguing investigators by the plethora of complex biologic pathways it regulates. These include several key aspects of mRNA metabolism, such as mRNA splicing, stability, and translation, as well as other mRNA-independent processes, such as miRNA processing, protein sumoylation, and the nucleolar stress response. In this review, the structural features of SRSF1 are discussed as they relate to the intricate mechanism of splicing and the multiplicity of functions it performs. Similarly, a list of relevant alternatively spliced transcripts and SRSF1 interacting proteins is provided. Finally, emphasis is given to the deleterious consequences of overexpression of the SRSF1 proto-oncogene in human cancers, and the complex mechanisms and pathways underlying SRSF1-mediated transformation. The accumulated knowledge about SRSF1 provides critical insight into the integral role it plays in maintaining cellular homeostasis and suggests new targets for anticancer therapy. Mol Cancer Res; 12(9); 1195-204. ©2014 AACR.

Homologous SV40 RNA trans-splicing: a new mechanism for diversification of viral sequences and phenotypes

Simian Virus 40 (SV40) is a polyomavirus found in both monkeys and humans, which causes cancer in some animal models. In humans, SV40 has been reported to be associated with cancers but causality has not been proven yet. The transforming activity of SV40 is mainly due to its 94-kD large T antigen, which binds to the retinoblastoma (pRb) and p53 tumor suppressor proteins, and thereby perturbs their functions. Here we describe a 100 kD super T antigen harboring a duplication of the pRB binding domain that was associated with unusual high cell transformation activity and that was generated by a novel mechanism involving homologous RNA trans-splicing of SV40 early transcripts in transformed rodent cells. Enhanced trans-splice activity was observed in clones carrying a single point mutation in the large T antigen 5' donor splice site (ss). This mutation impaired cis-splicing in favor of an alternative trans-splice reaction via a cryptic 5'ss within a second cis-spliced SV40 pre-mRNA molecule and enabled detectable gene expression. Next to the cryptic 5'ss we identified additional trans-splice helper functions, including putative dimerization domains and a splice enhancer sequence. Our findings suggest RNA trans-splicing as a SV40-intrinsic mechanism that supports the diversification of viral RNA and phenotypes.

Cauli: a mouse strain with an Ift140 mutation that results in a skeletal ciliopathy modelling Jeune syndrome

Cilia are architecturally complex organelles that protrude from the cell membrane and have signalling, sensory and motility functions that are central to normal tissue development and homeostasis. There are two broad categories of cilia; motile and non-motile, or primary, cilia. The central role of primary cilia in health and disease has become prominent in the past decade with the recognition of a number of human syndromes that result from defects in the formation or function of primary cilia. This rapidly growing class of conditions, now known as ciliopathies, impact the development of a diverse range of tissues including the neural axis, craniofacial structures, skeleton, kidneys, eyes and lungs. The broad impact of cilia dysfunction on development reflects the pivotal position of the primary cilia within a signalling nexus involving a growing number of growth factor systems including Hedgehog, Pdgf, Fgf, Hippo, Notch and both canonical Wnt and planar cell polarity. We have identified a novel ENU mutant allele of Ift140, which causes a mid-gestation embryonic lethal phenotype in homozygous mutant mice. Mutant embryos exhibit a range of phenotypes including exencephaly and spina bifida, craniofacial dysmorphism, digit anomalies, cardiac anomalies and somite patterning defects. A number of these phenotypes can be attributed to alterations in Hedgehog signalling, although additional signalling systems are also likely to be involved. We also report the identification of a homozygous recessive mutation in IFT140 in a Jeune syndrome patient. This ENU-induced Jeune syndrome model will be useful in delineating the origins of dysmorphology in human ciliopathies.

Skeletal ciliopathies are an emerging field of human disease in which skeletal birth defects arise due to abnormal communication between cells. This failure in communication arises following mutation in components of the primary cilia, a hair-like structure present on every cell. The skeletal ciliopathies are debilitating and in severe cases lead to death in early infancy. However, the mechanisms by which these malformations come about remains unclear. Mouse models are often used to delineate the causes of human birth defects and we have identified a model that mimics one of these conditions known as Jeune syndrome. It is the first mouse model with a mutation in the Ift140 gene, and these mice exhibit phenotypes that are often seen in this set of human syndromes. We have complimented this study with the discovery of a patient that presents with Jeune Syndrome resulting from mutation of human IFT140. This model will allow us to explore the role of IFT140 and the primary cilia in normal human development and provide insight into the field of human skeletal ciliopathies.

Production of ACAT1 56-kDa isoform in human cells via trans-splicing involving the ampicillin resistance gene

Trans-splicing, a process involving the cleavage and joining of two separate transcripts, can expand the transcriptome and proteome in eukaryotes. Chimeric RNAs generated by trans-splicing are increasingly described in literatures. The widespread presence of antibiotic resistance genes in natural environments and human intestines is becoming an important challenge for public health. Certain antibiotic resistance genes, such as ampicillin resistance gene (Amp^r), are frequently used in recombinant plasmids. Until now, trans-splicing involving recombinant plasmid-derived exogenous transcripts and endogenous cellular RNAs has not been reported. Acyl-CoA:cholesterol acyltransferase 1 (ACAT1) is a key enzyme involved in cellular cholesterol homeostasis. The 4.3-kb human ACAT1 chimeric mRNA can produce 50-kDa and 56-kDa isoforms with different enzymatic activities. Here, we show that human ACAT1 56-kDa isoform is produced from an mRNA species generated through the trans-splicing of an exogenous transcript encoded by the antisense strand of Amp^r (asAmp) present in common Amp^r-plasmids and the 4.3-kb endogenous ACAT1 chimeric mRNA, which is presumably processed through a prior event of interchromosomal trans-splicing. Strikingly, DNA fragments containing the asAmp with an upstream recombined cryptic promoter and the corresponding exogenous asAmp transcripts have been detected in human cells. Our findings shed lights on the mechanism of human ACAT1 56-kDa isoform production, reveal an exogenous-endogenous trans-splicing system, in which recombinant plasmid-derived exogenous transcripts are linked with endogenous cellular RNAs in human cells, and suggest that exogenous DNA might affect human gene expression at both DNA and RNA levels.

Genome-wide analysis reveals SR protein cooperation and competition in regulated splicing

SR proteins are well-characterized RNA binding proteins that promote exon inclusion by binding to exonic splicing enhancers (ESEs). However, it has been unclear whether regulatory rules deduced on model genes apply generally to activities of SR proteins in the cell. Here, we report global analyses of two prototypical SR proteins, SRSF1 (SF2/ASF) and SRSF2 (SC35), using splicing-sensitive arrays and CLIP-seq on mouse embryo fibroblasts (MEFs). Unexpectedly, we find that these SR proteins promote both inclusion and skipping of exons in vivo, but their binding patterns do not explain such opposite responses. Further analyses reveal that loss of one SR protein is accompanied by coordinated loss or compensatory gain in the interaction of other SR proteins at the affected exons. Therefore, specific effects on regulated splicing by one SR protein actually depend on a complex set of relationships with multiple other SR proteins in mammalian genomes.

Combined computational-experimental analyses of CFTR exon strength uncover predictability of exon-skipping level

With the increased number of identified nucleotide sequence variations in genes, the current challenge is to classify them as disease causing or neutral. These variants of unknown clinical significance can alter multiple processes, from gene transcription to RNA splicing or protein function. Using an approach combining several in silico tools, we identified some exons presenting weaker splicing motifs than other exons in the Cystic Fibrosis Transmembrane conductance Regulator (CFTR) gene. These exons exhibit higher rates of basal skipping than exons harboring no identifiable weak splicing signals using minigene assays. We then screened 19 described mutations in three different exons, and identified exon-skipping substitutions. These substitutions induced higher skipping levels in exons having one or more weak splicing motifs. Indeed, this level remained under 2% for exons with strong splicing motifs and could reach 40% for exons having at least one weak motif. Further analysis revealed a functional exon splicing enhancer within exon 3 that was associated with the SR protein SF2/ASF and whose disruption induced exon skipping. Exon skipping was confirmed in vivo in two nasal epithelial cell brushing samples. Our approach, which point out exons with some splicing signals weaknesses, will help spot splicing mutations of clinical relevance.

Multiple tandem splicing silencer elements suppress aberrant splicing within the long exon 26 of the human Apolipoprotein B gene

BACKGROUND

Apolipoprotein B (APOB) is an integral component of the chylomicron and the atherogenic lipoproteins LDL and Lp(a). Exon 26 of the APOB pre-mRNA is unusually long at 7,572 nt and is constitutively spliced. It is also subject to RNA editing in the intestine, which generates a shortened isoform, APOB48, assembled exclusively into chylomicrons. Due to its length, exon 26 contains multiple pseudo splice sites which are not spliced, but which conform to the degenerate splice site consensus.

RESULTS

We demonstrate that these pseudo splice sites are repressed by multiple, tandem splicing silencers distributed along the length of exon 26. The distribution of these elements appears to be heterogeneous, with a greater frequency in the middle 4,800 nt of the exon.

CONCLUSION

Repression of these splice sites is key to maintaining the integrity of exon 26 during RNA splicing and therefore the correct expression of both isoforms of APOB.

Function of alternative splicing

Almost all polymerase II transcripts undergo alternative pre-mRNA splicing. Here, we review the functions of alternative splicing events that have been experimentally determined. The overall function of alternative splicing is to increase the diversity of mRNAs expressed from the genome. Alternative splicing changes proteins encoded by mRNAs, which has profound functional effects. Experimental analysis of these protein isoforms showed that alternative splicing regulates binding between proteins, between proteins and nucleic acids as well as between proteins and membranes. Alternative splicing regulates the localization of proteins, their enzymatic properties and their interaction with ligands. In most cases, changes caused by individual splicing isoforms are small. However, cells typically coordinate numerous changes in 'splicing programs', which can have strong effects on cell proliferation, cell survival and properties of the nervous system. Due to its widespread usage and molecular versatility, alternative splicing emerges as a central element in gene regulation that interferes with almost every biological function analyzed.

RNA helicases in splicing

In eukaryotic cells, introns are spliced from pre-mRNAs by the spliceosome. Both the composition and the structure of the spliceosome are highly dynamic, and eight DExD/H RNA helicases play essential roles in controlling conformational rearrangements. There is evidence that the various helicases are functionally and physically connected with each other and with many other factors in the spliceosome. Understanding the dynamics of those interactions is essential to comprehend the mechanism and regulation of normal as well as of pathological splicing. This review focuses on recent advances in the characterization of the splicing helicases and their interactions, and highlights the deep integration of splicing helicases in global mRNP biogenesis pathways.