Pre-mRNA splicing during transcription in the mammalian system

Functional impact of splicing variants in the elaboration of complex traits in cattle

GWAS conducted directly on imputed whole genome sequence have led to the identification of numerous genetic variants associated with agronomic traits in cattle. However, such variants are often simply markers in linkage disequilibrium with the actual causal variants, which is a limiting factor for the development of accurate genomic predictions. It is possible to identify causal variants by integrating information on how variants impact gene expression into GWAS output. RNA splicing plays a major role in regulating gene expression. Thus, assessing the effect of variants on RNA splicing may explain their function. Here, we use a high-throughput strategy to functionally analyse putative splice-disrupting variants in the bovine genome. Using GWAS, massively parallel reporter assay and deep learning algorithms designed to predict splice-disrupting variants, we identify 38 splice-disrupting variants associated with complex traits in cattle, three of which could be classified as causal. Our results indicate that splice-disrupting variants are widely found in the quantitative trait loci related to these phenotypes. Using our combined approach, we also assess the validity of splicing predictors originally developed to analyse human variants in the context of the bovine genome.

Regulation of age-dependent expression patterns of five transcription factors in Larix kaempferi

To reveal the regulatory mechanisms underlying age-dependent expression patterns, we characterized seven age-related genes, LaDAL1, LaAGL2-2, LaAGL2-3, LaAGL11, LaSOC1-1, LaAP2-1, and LaAP2-2 in terms of transcription and intron splicing in Larix kaempferi. Based on the exon-intron structures, we quantified the pre-mRNA levels and mature mRNA levels of these seven genes using quantitative reverse transcription polymerase chain reaction experiments. We found that the pre-mRNA levels manifested age-related patterns, indicating that their transcription was primarily regulated by age. By comparing the increasing or decreasing rates of the pre-mRNA and spliced mRNA levels, we found that their splicing efficiencies also changed with age. These results clearly show that both pre-mRNA transcription and splicing of five age-related genes are regulated by age, indicating that age-dependent expression patterns are controlled at both transcriptional and post-transcriptional levels, and unveiling the underlying regulatory molecular mechanisms should focus on the transcription factors, epigenetic regulation, and RNA splicing. These data provide new insights into the age-mediated regulation of gene expression in woody perennials in terms of longevity.

Isoginkgetin leads to decreased protein synthesis and activates an ATF4-dependent transcriptional response

Isoginkgetin (IGG) is a small molecule inhibitor of pre-mRNA splicing. Failure to accurately remove introns could lead to the production of aberrant mRNAs and proteins. The cellular responses to splicing stress are not well defined. Here, we used oligonucleotide microarrays to assess genome wide changes in gene expression associated with exposure to IGG. Two of the 3 enriched pathways identified using PANTHER analysis of differentially expressed transcripts are linked to the ATF4 transcription factor. We confirmed that ATF4 was selectively translated and upregulated in response IGG despite an almost complete block to total protein synthesis. Importantly, partial disruption of the ATF4 gene using CRISPR-mediated gene editing prevented IGG-induced changes in gene expression. Remarkably, another spliceosome inhibitor, pladienolide B, did not inhibit translation, activate ATF4 or increase ATF4-dependent gene expression. Taken together, IGG activates ATF4 and an ATF4-dependent transcriptional response but these effects are not common to all spliceosome inhibitors.

The atypical RNA-binding protein Taf15 regulates dorsoanterior neural development through diverse mechanisms in Xenopus tropicalis

The FET family of atypical RNA-binding proteins includes Fused in sarcoma (FUS), Ewing's sarcoma (EWS) and the TATA-binding protein-associate factor 15 (TAF15). FET proteins are highly conserved, suggesting specialized requirements for each protein. Fus regulates splicing of transcripts required for mesoderm differentiation and cell adhesion in Xenopus, but the roles of Ews and Taf15 remain unknown. Here, we analyze the roles of maternally deposited and zygotically transcribed Taf15, which is essential for the correct development of dorsoanterior neural tissues. By measuring changes in exon usage and transcript abundance from Taf15-depleted embryos, we found that Taf15 may regulate dorsoanterior neural development through fgfr4 and ventx2.1. Taf15 uses distinct mechanisms to downregulate Fgfr4 expression, namely retention of a single intron within fgfr4 when maternal and zygotic Taf15 is depleted, and reduction in the total fgfr4 transcript when zygotic Taf15 alone is depleted. The two mechanisms of gene regulation (post-transcriptional versus transcriptional) suggest that Taf15-mediated gene regulation is target and co-factor dependent, contingent on the milieu of factors that are present at different stages of development.

Cis-regulatory elements in conserved non-coding sequences of nuclear receptor genes indicate for crosstalk between endocrine systems

Nuclear receptors (NRs) are ligand-activated transcription factors that regulate gene expression when bound to specific DNA sequences. Crosstalk between steroid NR systems has been studied for understanding the development of hormone-driven cancers but not to an extent at a genetic level. This study aimed to investigate crosstalk between steroid NRs in conserved intron and exon sequences, with a focus on steroid NRs involved in prostate cancer etiology. For this purpose, we evaluated conserved intron and exon sequences among all 49 members of the NR Superfamily (NRS) and their relevance as regulatory sequences and NR-binding sequences. Sequence conservation was found to be higher in the first intron (35%), when compared with downstream introns. Seventy-nine percent of the conserved regions in the NRS contained putative transcription factor binding sites (TFBS) and a large fraction of these sequences contained splicing sites (SS). Analysis of transcription factors binding to putative intronic and exonic TFBS revealed that 5 and 16%, respectively, were NRs. The present study suggests crosstalk between steroid NRs, e.g., vitamin D, estrogen, progesterone, and retinoic acid endocrine systems, through cis-regulatory elements in conserved sequences of introns and exons. This investigation gives evidence for crosstalk between steroid hormones and contributes to novel targets for steroid NR regulation.

Biological Aspects of Cancer Health Disparities

Racial and ethnic disparities span the continuum of cancer care and are driven by a complex interplay among social, psychosocial, lifestyle, environmental, health system, and biological determinants of health. Research is needed to identify these determinants of cancer health disparities and to develop interventions to achieve cancer health equity. Herein, we focus on the overall burden of ancestry-related molecular alterations, the functional significance of the alterations in hallmarks of cancer, and the implications of the alterations for precision oncology and immuno-oncology. In conclusion, we reflect on the importance of estimating ancestry, improving diverse racial and ethnic participation in cancer clinical trials, and examining the intersection among determinants of cancer health disparities.

Intermittent leucine pulses during continuous feeding alters novel components involved in skeletal muscle growth of neonatal pigs

When neonatal pigs continuously fed formula are supplemented with leucine pulses, muscle protein synthesis and body weight gain are enhanced. To identify the responsible mechanisms, we combined plasma metabolomic analysis with transcriptome expression of the transcriptome and protein catabolic pathways in skeletal muscle. Piglets (n = 23, 7-day-old) were fed continuously a milk replacement formula via orogastric tube for 21 days with an additional parenteral infusion (800 μmol kg^-1 h^-1) of either leucine (LEU) or alanine (CON) for 1 h every 4 h. Plasma metabolites were measured by liquid chromatography-mass spectrometry. Gene and protein expression analyses of longissimus dorsi muscle were performed by RNA-seq and Western blot, respectively. Compared with CON, LEU pigs had increased plasma levels of leucine-derived metabolites, including 4-methyl-2-oxopentanoate, beta-hydroxyisovalerate, β-hydroxyisovalerylcarnitine, and 3-methylglutaconate (P ≤ 0.05). Leucine pulses downregulated transcripts enriched in the Kyoto Encyclopedia of Genes and Genomes terms "spliceosome," "GAP junction," "endocytosis," "ECM-receptor interaction," and "DNA replication". Significant correlations were identified between metabolites derived from leucine catabolism and muscle genes involved in protein degradation, transcription and translation, and muscle maintenance and development (P ≤ 0.05). Further, leucine pulses decreased protein expression of autophagic markers and serine/threonine kinase 4, involved in muscle atrophy (P ≤ 0.01). In conclusion, results from our studies support the notion that leucine pulses during continuous enteral feeding enhance muscle mass gain in neonatal pigs by increasing protein synthetic activity and downregulating protein catabolic pathways through concerted responses in the transcriptome and metabolome.

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.

The Role of Alternative RNA Splicing in the Regulation of hTERT, Telomerase, and Telomeres: Implications for Cancer Therapeutics

Alternative RNA splicing impacts the majority (>90%) of eukaryotic multi-exon genes, expanding the coding capacity and regulating the abundance of gene isoforms. Telomerase (hTERT) is a key example of a gene that is alternatively spliced during human fetal development and becomes dysregulated in nearly all cancers. Approximately 90% of human tumors use telomerase to synthesize de novo telomere repeats and obtain telomere-dependent cellular immortality. Paradigm shifting data indicates that hTERT alternative splicing, in addition to transcription, plays an important role in the regulation of active telomerase in cells. Our group and others are pursuing the basic science studies to progress this emerging area of telomerase biology. Recent evidence demonstrates that switching splicing of hTERT from the telomerase activity producing full-length hTERT isoform to alternatively spliced, non-coding isoforms may be a novel telomerase inhibition strategy to prevent cancer growth and survival. Thus, the goals of this review are to detail the general roles of telomerase in cancer development, explore the emerging regulatory mechanisms of alternative RNA splicing of the hTERT gene in various somatic and cancer cell types, define the known and potential roles of hTERT splice isoforms in cancer cell biology, and provide insight into new treatment strategies targeting hTERT in telomerase-positive cancers.

Small non-coding RNA within the endogenous spliceosome and alternative splicing regulation

Splicing and alternative splicing (AS), which occur in the endogenous spliceosome, play major roles in regulating gene expression, and defects in them are involved in numerous human diseases including cancer. Although the mechanism of the splicing reaction is well understood, the regulation of AS remains to be elucidated. A group of essential regulatory factors in gene expression are small non-coding RNAs (sncRNA): e.g. microRNA, mainly known for their inhibitory role in translation in the cytoplasm; and small nucleolar RNA, known for their role in methylating non-coding RNA in the nucleolus. Here I highlight a new aspect of sncRNAs found within the endogenous spliceosome. Assembled in non-canonical complexes and through different base pairing than their canonical ones, spliceosomal sncRNAs can potentially target different RNAs. Examples of spliceosomal sncRNAs regulating AS, regulating gene expression, and acting in a quality control of AS are reviewed, suggesting novel functions for spliceosomal sncRNAs. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

Shorter Phosphorodiamidate Morpholino Splice-Switching Oligonucleotides May Increase Exon-Skipping Efficacy in DMD

Duchenne muscular dystrophy is a fatal muscle disease, caused by mutations in DMD, leading to loss of dystrophin expression. Phosphorodiamidate morpholino splice-switching oligonucleotides (PMO-SSOs) have been used to elicit the restoration of a partially functional truncated dystrophin by excluding disruptive exons from the DMD messenger. The 30-mer PMO eteplirsen (EXONDYS51) developed for exon 51 skipping is the first dystrophin-restoring, conditionally FDA-approved drug in history. Clinical trials had shown a dose-dependent variable and patchy dystrophin restoration. The main obstacle for efficient dystrophin restoration is the inadequate uptake of PMOs into skeletal muscle fibers at low doses. The excessive cost of longer PMOs has limited the utilization of higher dosing. We designed shorter 25-mer PMOs directed to the same eteplirsen-targeted region of exon 51 and compared their efficacies in vitro and in vivo in the mdx52 murine model. Our results showed that skipped-dystrophin induction was comparable between the 30-mer PMO sequence of eteplirsen and one of the shorter PMOs, while the other 25-mer PMOs showed lower exon-skipping efficacies. Shorter PMOs would make higher doses economically feasible, and high dosing would result in better drug uptake into muscle, induce higher levels of dystrophin restoration in DMD muscle, and, ultimately, increase the clinical efficacy.

Splicing dysregulation as a driver of breast cancer

Breast cancer is known to be a heterogeneous disease driven by a large repertoire of molecular abnormalities, which contribute to its diverse clinical behaviour. Despite the success of targeted therapy approaches for breast cancer patient management, there is still a lack of the molecular understanding of aggressive forms of the disease and clinical management of these patients remains difficult. The advent of high-throughput sequencing technologies has paved the way for a more complete understanding of the molecular make-up of the breast cancer genome. As such, it is becoming apparent that disruption of canonical splicing within breast cancer governs its clinical progression. In this review, we discuss the role of dysregulation of spliceosomal component genes and associated factors in the progression of breast cancer, their role in therapy resistance and the use of quantitative isoform expression as potential prognostic and predictive biomarkers with a particular focus on oestrogen receptor-positive breast cancer.

Role of cleavage and polyadenylation specificity factor 100: anchoring poly(A) sites and modulating transcription termination

CPSF100 is a core component of the cleavage and polyadenylation specificity factor (CPSF) complex for 3'-end formation of mRNA, but it still has no clear functional assignment. CPSF100 was reported to play a role in RNA silencing and promote flowering in Arabidopsis. However, the molecular mechanisms underlying these phenomena are not fully understood. Our genetics analyses indicate that plants with a hypomorphic mutant of CPSF100 (esp5) show defects in embryogenesis, reduced seed production or altered root morphology. To unravel this puzzle, we employed a poly(A) tag sequencing protocol and uncovered a different poly(A) profile in esp5. This transcriptome-wide analysis revealed alternative polyadenylation of thousands of genes, most of which result in transcriptional read-through in protein-coding genes. AtCPSF100 also affects poly(A) signal recognition on the far-upstream elements; in particular it prefers less U-rich sequences. Importantly, AtCPSF100 was found to exert its functions through the change of poly(A) sites on genes encoding binding proteins, such as nucleotide-binding, RNA-binding and poly(U)-binding proteins. In addition, through its interaction with RNA Polymerase II C-terminal domain (CTD) and affecting the expression level of CTD phosphatase-like 3 (CPL3), AtCPSF100 is shown to potentially ensure transcriptional termination by dephosphorylation of Ser2 on the CTD. These data suggest a key role for CPSF100 in locating poly(A) sites and affecting transcription termination.

The nuts and bolts of the endogenous spliceosome

The complex life of pre-mRNA from transcription to the production of mRNA that can be exported from the nucleus to the cytoplasm to encode for proteins entails intricate coordination and regulation of a network of processing events. Coordination is required between transcription and splicing and between several processing events including 5' and 3' end processing, splicing, alternative splicing and editing that are major contributors to the diversity of the human proteome, and occur within a huge and dynamic macromolecular machine-the endogenous spliceosome. Detailed mechanistic insight of the splicing reaction was gained from studies of the in vitro spliceosome assembled on a single intron. Because most pre-mRNAs are multiintronic that undergo alternative splicing, the in vivo splicing machine requires additional elements to those of the in vitro machine, to account for all these diverse functions. Information about the endogenous spliceosome is emerging from imaging studies in intact and live cells that support the cotranscriptional commitment to splicing model and provide information about splicing kinetics in vivo. Another source comes from studies of the in vivo assembled spliceosome, isolated from cell nuclei under native conditions-the supraspliceosome-that individually package pre-mRNA transcripts of different sizes and number of introns into complexes of a unique structure, indicating their universal nature. Recent years have portrayed new players affecting alternative splicing and novel connections between splicing, transcription and chromatin. The challenge ahead is to elucidate the structure and function of the endogenous spliceosome and decipher the regulation and coordination of its network of processing activities. WIREs RNA 2017, 8:e1377. doi: 10.1002/wrna.1377 For further resources related to this article, please visit the WIREs website.

A Splicing Reporter Tuned to Non-AG Acceptor Sites Reveals that Luteolin Enhances the Recognition of Non-canonical Acceptor Sites

Removal of an intron requires precise recognition of the splice donor and acceptor sites located at the 5' and 3' termini of introns. Although the roles of these sequences differ, mutations in both sites easily block normal splicing and produce an aberrant mRNA. For example, many splice-site mutations occur in patients with inherited diseases. Several approaches have been evaluated to restore expression of a functional protein; however, because of the strict requirement for an AG dinucleotide at the 3' terminus of a U2-type intron, no method is available to correct splicing at a mutated sequence. To identify compounds that allow splicing at the non-AG acceptor site, in the present study we constructed a reporter gene with a modified polypyrimidine tract. However, the modified polypyrimidine tract mediated splicing at adjacent non-canonical acceptor sites, including the original mutated site. Further, we show that certain flavones such as luteolin and apigenin enhanced aberrant splicing at the non-canonical acceptor site of the reporter gene. These results suggest that the reporter gene and luteolin may be useful for further screening to identify molecules that correct aberrant splicing caused by a disease-associated splice acceptor site mutation.

Quantification of co-transcriptional splicing from RNA-Seq data

During gene expression, protein-coding transcripts are shaped by multiple processing events: 5' end capping, pre-mRNA splicing, RNA editing, and 3' end cleavage and polyadenylation. These events are required to produce mature mRNA, which can be subsequently translated. Nearly all of these RNA processing steps occur during transcription, while the nascent RNA is still attached to the DNA template by RNA polymerase II (i.e. co-transcriptionally). Polyadenylation occurs after 3' end cleavage or post-transcriptionally. Pre-mRNA splicing - the removal of introns and ligation of exons - can be initiated and concluded co-transcriptionally, although this is not strictly required. Recently, a number of studies using global methods have shown that the majority of splicing is co-transcriptional, yet not all published studies agree in their conclusions. Short read sequencing of RNA (RNA-Seq) is the prevailing approach to measuring splicing levels in nascent RNA, mRNA or total RNA. Here, we compare four different strategies for analyzing and quantifying co-transcriptional splicing. To do so, we reanalyze two nascent RNA-Seq datasets of the same species, but different cell type and RNA isolation procedure. Average co-transcriptional splicing values calculated on a per intron basis are similar, independent of the strategy used. We emphasize the technical requirements for identifying co-transcriptional splicing events with high confidence, e.g. how to calculate co-transcriptional splicing from nascent RNA- versus mRNA-Seq data, the number of biological replicates needed, depletion of polyA+RNA, and appropriate normalization. Finally, we present guidelines for planning a nascent RNA-Seq experiment.

Detained introns are a novel, widespread class of post-transcriptionally spliced introns

Deep sequencing of embryonic stem cell RNA revealed many specific internal introns that are significantly more abundant than the other introns within polyadenylated transcripts. Boutz et al. identified thousands of these “detained” introns (DIs) in human and mouse cell lines as well as the adult mouse liver. Drug inhibition of Clk, a stress-responsive kinase, triggered rapid splicing changes for a specific subset of DIs, altering transcript pools of >300 genes. Srsf4 regulates the splicing of some DIs, particularly in genes encoding RNA processing and splicing factors.

Deep sequencing of embryonic stem cell RNA revealed many specific internal introns that are significantly more abundant than the other introns within polyadenylated transcripts; we classified these as “detained” introns (DIs). We identified thousands of DIs, many of which are evolutionarily conserved, in human and mouse cell lines as well as the adult mouse liver. DIs can have half-lives of over an hour yet remain in the nucleus and are not subject to nonsense-mediated decay (NMD). Drug inhibition of Clk, a stress-responsive kinase, triggered rapid splicing changes for a specific subset of DIs; half showed increased splicing, and half showed increased intron detention, altering transcript pools of >300 genes. Srsf4, which undergoes a dramatic phosphorylation shift in response to Clk kinase inhibition, regulates the splicing of some DIs, particularly in genes encoding RNA processing and splicing factors. The splicing of some DIs—including those in Mdm4, a negative regulator of p53—was also altered following DNA damage. After 4 h of Clk inhibition, the expression of >400 genes changed significantly, and almost one-third of these are p53 transcriptional targets. These data suggest a widespread mechanism by which the rate of splicing of DIs contributes to the level of gene expression.

Identification of a new isoform of the murine Sh2d1a gene and its functional implications

Signaling lymphocytic activation molecule (SLAM)-associated protein (SAP) is a Src homology (SH) domain 2-containing intracellular adaptor protein that is predominantly expressed in the hematopoietic system by T lymphocytes and NK cells. SAP protein is encoded by the SH2D1A gene located on the X chromosome. Loss-of-function mutations in SAP cause the X-linked lymphoproliferative disease (XLP), a severe immunodeficiency characterized by heightened susceptibility to Epstein-Barr virus and impaired humoral immunity. Normal individuals express several functional and non-functional isoforms of SAP as a result of alternative splicing. In this study, we identify a cryptic exon in the murine Sh2d1a gene. At the mRNA level, the new isoform of SAP (SAP-2) that includes this new exon is widely expressed in lymphoid tissues by C57BL/6 and 129 strains of inbred mice. SAP-2 accounts for approximately 1%-3% of total SAP transcripts, and it is dynamically regulated during lymphocyte activation. At the protein level, the SAP-2 isoform is a 144 amino-acid protein. Compared to the dominant 126 aminoacid SAP-1 isoform, the additional 18 amino acids are inserted into a structural region that is critical for phosphotyrosine binding. Our functional analysis in vitro indicates that SAP-2 is a non-functional isoform due to decreased protein stability. Thus, both human and mouse have multiple SAP splice isoforms that may or may not function. Modulation of relative proportions of these isoforms is potentially a mechanism whereby cells can regulate SAP-mediated biological activities.

Splicing kinetics and transcript release from the chromatin compartment limit the rate of Lipid A-induced gene expression

The expression of eukaryotic mRNAs is achieved though an intricate series of molecular processes that provide many steps for regulating the production of a final gene product. However, the relationships between individual steps in mRNA biosynthesis and the rates at which they occur are poorly understood. By applying RNA-seq to chromatin-associated and soluble nucleoplasmic fractions of RNA from Lipid A-stimulated macrophages, we examined the timing of exon ligation and transcript release from chromatin relative to the induction of transcription. We find that for a subset of genes in the Lipid A response, the ligation of certain exon pairs is delayed relative to the synthesis of the complete transcript. In contrast, 3' end cleavage and polyadenylation occur rapidly once transcription extends through the cleavage site. Our data indicate that these transcripts with delayed splicing are not released from the chromatin fraction until all the introns have been excised. These unusual kinetics result in a chromatin-associated pool of completely transcribed and 3'-processed transcripts that are not yet fully spliced. We also find that long introns containing repressed exons that will be excluded from the final mRNA are excised particularly slowly relative to other introns in a transcript. These results indicate that the kinetics of splicing and transcript release contribute to the timing of expression for multiple genes of the inflammatory response.

Counting on co-transcriptional splicing

Splicing is the removal of intron sequences from pre-mRNA by the spliceosome. Researchers working in multiple model organisms – notably yeast, insects and mammalian cells – have shown that pre-mRNA can be spliced during the process of transcription (i.e. co-transcriptionally), as well as after transcription termination (i.e. post-transcriptionally). Co-transcriptional splicing does not assume that transcription and splicing machineries are mechanistically coupled, yet it raises this possibility. Early studies were based on a limited number of genes, which were often chosen because of their experimental accessibility. Since 2010, eight studies have used global datasets as counting tools, in order to quantify co-transcriptional intron removal. The consensus view, based on four organisms, is that the majority of splicing events take place co-transcriptionally in most cells and tissues. Here, we discuss the nature of the various global datasets and how bioinformatic analyses were conducted. Considering the broad differences in experimental approach and analysis, the level of agreement on the prevalence of co-transcriptional splicing is remarkable.

Neuronal regulation of pre-mRNA splicing by polypyrimidine tract binding proteins, PTBP1 and PTBP2

Alternative splicing patterns are regulated by RNA binding proteins that assemble onto each pre-mRNA to form a complex RNP structure. The polypyrimidine tract binding protein, PTB, has served as an informative model for understanding how RNA binding proteins affect spliceosome assembly and how changes in the expression of these proteins can control complex programs of splicing in tissues. In this review, we describe the mechanisms of splicing regulation by PTB and its function, along with its paralog PTBP2, in neuronal development.