Pre-mRNA splicing in plants: characterization of Ser/Arg splicing factors.

The fact that animal introns are not spliced out in plants suggests that recognition of pre-mRNA splice sites differs between the two kingdoms. In plants, little is known about proteins required for splicing, as no plant in vitro splicing system is available. Several essential splicing factors from animals, such as SF2/ASF and SC-35, belong to a family of highly conserved proteins consisting of one or two RNA binding domain(s) (RRM) and a C-terminal Ser/Arg-rich (SR or RS) domain. These animal SR proteins are required for splice site recognition and spliceosome assembly. We have screened for similar proteins in plants by using monoclonal antibodies specific for a phosphoserine epitope of the SR proteins (mAb1O4) or for SF2/ASF. These experiments demonstrate that plants do possess SR proteins, including SF2/ASF-like proteins. Similar to the animal SR proteins, this group of proteins can be isolated by two salt precipitations. However, compared to the animal SR proteins, which are highly conserved in size and number, SR proteins from Arabidopsis, carrot, and tobacco exhibit a complex pattern of intra- and interspecific variants. These plant SR proteins are able to complement inactive HeLa cell cytoplasmic S1OO extracts that are deficient in SR proteins, yielding functional splicing extracts. In addition, plant SR proteins were active in a heterologous alternative splicing assay. Thus, these plant SR proteins are authentic plant splicing factors.

An RNA structure involved in feedback regulation of splicing and of translation is critical for biological fitness.

While studies of the regulation of gene expression have generally concerned qualitative changes in the selection or the level of expression of a gene, much of the regulation that occurs within a cell involves the continuous subtle optimization of the levels of proteins used in macromolecular complexes. An example is the biosynthesis of the ribosome, in which equimolar amounts of nearly 80 ribosomal proteins must be supplied by the cytoplasm to the nucleolus. We have found that the transcript of one of the ribosomal protein genes of Saccharomyces cerevisiae, RPL32, participates in such fine tuning. Sequences from exon I of the RPL32 transcript interact with nucleotides from the intron to form a structure that binds L32 to regulate splicing. In the spliced transcript, the same sequences interact with nucleotides from exon II to form a structure that binds L32 to regulate translation, thus providing two levels of autoregulation. We now show, by using a sensitive cocultivation assay, that these RNA structures and their interaction with L32 play a role in the fitness of the cell. The change of a single nucleotide within the 5' leader of the RPL32 transcript, which abolishes the site for L32 binding, leads to detectably slower growth and to eventual loss of the mutant strain from the culture. Experiments designed to assess independently the regulation of splicing and the regulation of translation are presented. These observations demonstrate that, in evolutionary terms, subtle regulatory compensations can be critical. The change in structure of an RNA, due to alteration of just one noncoding nucleotide, can spell the difference between biological success and failure.

Structural organization of mouse peroxisome proliferator-activated receptor gamma (mPPAR gamma) gene: alternative promoter use and different splicing yield two mPPAR gamma isoforms.

To gain insight into the regulation of expression of peroxisome proliferator-activated receptor (PPAR) isoforms, we have determined the structural organization of the mouse PPAR gamma (mPPAR gamma) gene. This gene extends > 105 kb and gives rise to two mRNAs (mPPAR gamma 1 and mPPAR gamma 2) that differ at their 5' ends. The mPPAR gamma 2 cDNA encodes an additional 30 amino acids N-terminal to the first ATG codon of mPPAR gamma 1 and reveals a different 5' untranslated sequence. We show that mPPAR gamma 1 mRNA is encoded by eight exons, whereas the mPPAR gamma 2 mRNA is encoded by seven exons. Most of the 5' untranslated sequence of mPPAR gamma 1 mRNA is encoded by two exons, whereas the 5' untranslated sequence and the extra 30 N-terminal amino acids of mPPAR gamma 2 are encoded by one exon, which is located between the second and third exons coding for mPPAR gamma 1. The last six exons of mPPAR gamma gene code for identical sequences in mPPAR gamma 1 and mPPAR gamma 2 isoforms. The mPPAR gamma 1 and mPPAR gamma 2 isoforms are transcribed from different promoters. The mPPAR gamma gene has been mapped to chromosome 6 E3-F1 by in situ hybridization using a biotin-labeled probe. These results establish that at least one of the PPAR genes yields more than one protein product, similar to that encountered with retinoid X receptor and retinoic acid receptor genes. The existence of multiple PPAR isoforms transcribed from different promoters could increase the diversity of ligand and tissue-specific transcriptional responses.

Isolation of a differentially regulated splicing isoform of human NF-E2.

The transcription factor NF-E2 (nuclear factor erythroid 2), interacting via DNA motifs within regulatory regions of several hematopoietic genes, is thought to mediate the enhancer activity of the globin locus control regions. By screening a human fetal liver cDNA library with probes derived from mouse NF-E2, we have isolated a splicing variant of the NF-E2 gene (fNF-E2) that differs in the 5' untranslated region from the previously reported cDNA (aNF-E2). The fNF-E2 isoform is transcribed from an alternative promoter located in the 3' end of the first intron and joined by alternative splicing to the second and third exons, which are shared by both RNA isoforms. Although the two forms produce the same protein, they are expressed in different ratios during development. fNF-E2 is more abundant in the fetal liver and less abundant in the adult bone marrow compared to the previously described form. Their distribution apparently follows the differential expression of fetal and adult hemoglobins.

RNA-Seq expression profile of genes related to neurodegenerative disorders affecting the human retina

Human neurodegenerative diseases, such as Parkinson’s disease (PD) and the neuromuscular disorders called dystroglycanopathies (DGPs), cause retinal impairments. We have used RNA-Seq technology to catalog all known genes linked to PD and DGPs expressed in the human retina and quantitate their mRNA levels in terms of FPKM. We have also characterized their expression profiles in the retina by determining their exonic, intronic and exon-intron junction expression levels, as well as the alternative splicing pattern of particular genes. We believe these data could pave the way toward understanding the molecular bases of sight deficiencies associated with neurodegenerative disorders.

Alternative polyadenylation and splicing of mRNAs transcribed from the human Sin1 gene

Limited but significant sequence similarity has been observed between an uncharacterized human protein, SIN1, and the S. pombe SIN1, Dictyostelium RIP3 and S. cerevisiae AVO1 proteins. The human Sin1 gene has been automatically predicted (MAPKAP1; GenBank accession number NM_024117); however, this sequence appears to be incomplete. In this study, we have cloned and characterized the full-length human Sin1 mRNA and identified a highly conserved domain that defines the family of SIN1 orthologues, members of which are widely distributed in the fungal and metazoan kingdoms. We demonstrate that Sin1 transcripts can use alternative polyadenylation signals and describe a number of Sin1 splice variants that potentially encode functionally different isoforms. (C) 2004 Elsevier B.V. All rights reserved.

Molecular recognition of an RNA trafficking element by heterogeneous nuclear ribonucleoprotein A2

Heterogeneous nuclear ribonucleoprotein (hnRNP) A2 is a multitasking protein involved in RNA packaging, alternative splicing of pre-mRNA. telomere maintenance, cytoplasmic RNA trafficking, and translation. It binds short segments of single-stranded nucleic acids, including the A2RE11 RNA element that is necessary and sufficient for cytoplasmic transport of a subset of rnRNAs in oligodendrocytes and neurons. We have explored the structures of hnRNP A2, its RNA recognition motifs (RRMs) and Gly-rich module, and the RRM complexes with A2RE11. Circular dichroism spectroscopy showed that the secondary structure of the first 189 residues of hnRNP A2 parallels that of the tandem beta alpha beta beta alpha beta RRMs of its paralogue, hnRNP A1, previously deduced from X-ray diffraction studies. The unusual GRD was shown to have substantial beta-sheet and beta-turn structure. Sedimentation equilibrium and circular dichroism results were consistent with the tandem RRM region being monomeric and supported earlier evidence for the binding of two A2RE11 oligoribonucleotides to this domain, in contrast to the protein dimer formed by the complex of hnRNP A1 with the telomeric ssDNA repeat. A three-dimensional structure for the N-terminal, two-RRM-containing segment of hnRNP A2 was derived by homology modeling. This structure was used to derive a model for the complex with A2RE11 using the previously described interaction of pairs of stacked nucleotides with aromatic residues on the RRM beta-sheet platforms, conserved in other RRM-RNA complexes, together with biochemical data and molecular dynamics-based observations of inter-RRM mobility.

The role of splicing factors and small nuclear RNAs in spliceosomal formation

Protein coding genes are comprised of protein-coding exons and non-protein-coding introns. The process of splicing involves removal of the introns and joining of the exons to form a mature messenger RNA, which subsequently undergoes translation into polypeptide. The spliceosome is a large, RNA/protein assembly of five small nuclear RNAs as well as over 300 proteins, which catalyzes intron removal and exon ligation. The selection of specific exons for inclusion in the mature messenger RNA is spatiotemporally regulated and results in production of an enormous diversity of polypeptides from a single gene locus. This phenomenon, known as alternative splicing, is regulated, in part, by protein splicing factors, which target the spliceosome to exon/intron boundaries. The first part of my dissertation (Chapters II and III) focuses on the discovery and characterization of the 45 kilodalton FK506 binding protein (FKBP45), which I discovered in the silk moth, Bombyx mori, as a U1 small nuclear RNA binding protein. This protein family binds the immunosuppressants FK506 and rapamycin and contains peptidyl-prolyl cis-trans isomerase activity, which converts polypeptides from cis to trans about a proline residue. This is the first time that an FKBP has been identified in the spliceosome. The second section of my dissertation (Chapters IV, V, VI and VII) is an investigation of the potential role of small nuclear RNA sequence variants in the control of splicing. I identified 46 copies of small nuclear RNAs in the 6X whole genome shotgun of the Bombyx mori p50T strain. These variants may play a role in differential binding of specific proteins that mediate alternative splicing. Along these lines, further investigation of U2 snRNA sequence variants in Bombyx mori demonstrated that some U2 snRNAs preferentially assemble into high molecular weight spliceosomal complexes over others. Expression of snRNA variants may represent another mechanism by which the cell is able to fine tune the splicing process.

Identification of Conserved Splicing Motifs in Mutually Exclusive Exons of 15 Insect Species

Background: During alternative splicing, the inclusion of an exon in the final mRNA molecule is determined by nuclear proteins that bind cis-regulatory sequences in a target pre-mRNA molecule. A recent study suggested that the regulatory codes of individual RNA-binding proteins may be nearly immutable between very diverse species such as mammals and insects. The model system Drosophila melanogaster therefore presents an excellent opportunity for the study of alternative splicing due to the availability of quality EST annotations in FlyBase. Methods: In this paper, we describe an in silico analysis pipeline to extract putative exonic splicing regulatory sequences from a multiple alignment of 15 species of insects. Our method, ESTs-to-ESRs (E2E), uses graph analysis of EST splicing graphs to identify mutually exclusive (ME) exons and combines phylogenetic measures, a sliding window approach along the multiple alignment and the Welch’s t statistic to extract conserved ESR motifs. Results: The most frequent 100% conserved word of length 5 bp in different insect exons was “ATGGA”. We identified 799 statistically significant “spike” hexamers, 218 motifs with either a left or right FDR corrected spike magnitude p-value < 0.05 and 83 with both left and right uncorrected p < 0.01. 11 genes were identified with highly significant motifs in one ME exon but not in the other, suggesting regulation of ME exon splicing through these highly conserved hexamers. The majority of these genes have been shown to have regulated spatiotemporal expression. 10 elements were found to match three mammalian splicing regulator databases. A putative ESR motif, GATGCAG, was identified in the ME-13b but not in the ME-13a of Drosophila N-Cadherin, a gene that has been shown to have a distinct spatiotemporal expression pattern of spliced isoforms in a recent study. Conclusions: Analysis of phylogenetic relationships and variability of sequence conservation as implemented in the E2E spikes method may lead to improved identification of ESRs. We found that approximately half of the putative ESRs in common between insects and mammals have a high statistical support (p < 0.01). Several Drosophila genes with spatiotemporal expression patterns were identified to contain putative ESRs located in one exon of the ME exon pairs but not in the other.

The Role of Splicing Factors and Small Nuclear RNAS in Spliceosomal Formation

Protein coding genes are comprised of protein-coding exons and non-protein-coding introns. The process of splicing involves removal of the introns and joining of the exons to form a mature messenger RNA, which subsequently undergoes translation into polypeptide. The spliceosome is a large, RNA/protein assembly of five small nuclear RNAs as well as over 300 proteins, which catalyzes intron removal and exon ligation. The selection of specific exons for inclusion in the mature messenger RNA is spatio-temporally regulated and results in production of an enormous diversity of polypeptides from a single gene locus. This phenomenon, known as alternative splicing, is regulated, in part, by protein splicing factors, which target the spliceosome to exon/intron boundaries. The first part of my dissertation (Chapters II and III) focuses on the discovery and characterization of the 45 kilodalton FK506 binding protein (FKBP45), which I discovered in the silk moth, Bombyx mori, as a U1 small nuclear RNA binding protein. This protein family binds the immunosuppressants FK506 and rapamycin and contains peptidyl-prolyl cis-trans isomerase activity, which converts polypeptides from cis to trans about a proline residue. This is the first time that an FKBP has been identified in the spliceosome. The second section of my dissertation (Chapters IV, V, VI and VII) is an investigation of the potential role of small nuclear RNA sequence variants in the control of splicing. I identified 46 copies of small nuclear RNAs in the 6X whole genome shotgun of the Bombyx mori p50T strain. These variants may play a role in differential binding of specific proteins that mediate alternative splicing. Along these lines, further investigation of U2 snRNA sequence variants in Bombyx mori demonstrated that some U2 snRNAs preferentially assemble into high molecular weight spliceosomal complexes over others. Expression of snRNA variants may represent another mechanism by which the cell is able to fine tune the splicing process.

Analysis of Histone Deacetylase 7 (HDAC7) Alternative Splicing and Its Role in Embryonic Stem Cell Differentiation Toward Smooth Muscle Lineage

Histone deacetylases (HDACs) have a central role in the regulation of gene expression, which undergoes alternative splicing during embryonic stem cell (ES) cell differentiation. Alternative splicing gives rise to vast diversity over gene information, arousing public concerns in the last decade. In this chapter, we describe a strategy to detect HDAC7 alternative splicing and analyze its function on ES cell differentiation.

Cryptic splicing events in the iron transporter ABCB7 and other key target genes in SF3B1-mutant myelodysplastic syndromes

The splicing factor SF3B1 is the most frequently mutated gene in myelodysplastic syndromes (MDS), and is strongly associated with the presence of ring sideroblasts (RS). We have performed a systematic analysis of cryptic splicing abnormalities from RNA sequencing data on hematopoietic stem cells (HSCs) of SF3B1-mutant MDS cases with RS. Aberrant splicing events in many downstream target genes were identified and cryptic 3' splice site usage was a frequent event in SF3B1-mutant MDS. The iron transporter ABCB7 is a well-recognized candidate gene showing marked downregulation in MDS with RS. Our analysis unveiled aberrant ABCB7 splicing, due to usage of an alternative 3' splice site in MDS patient samples, giving rise to a premature termination codon in the ABCB7 mRNA. Treatment of cultured SF3B1-mutant MDS erythroblasts and a CRISPR/Cas9-generated SF3B1-mutant cell line with the nonsense-mediated decay (NMD) inhibitor cycloheximide showed that the aberrantly spliced ABCB7 transcript is targeted by NMD. We describe cryptic splicing events in the HSCs of SF3B1-mutant MDS, and our data support a model in which NMD-induced downregulation of the iron exporter ABCB7 mRNA transcript resulting from aberrant splicing caused by mutant SF3B1 underlies the increased mitochondrial iron accumulation found in MDS patients with RS.Leukemia advance online publication, 17 June 2016; doi:10.1038/leu.2016.149.

Alternative splicing takes control of cytokine signaling

Thesis (Ph.D.)--University of Washington, 2016-08

Poly (A)(+) Transcriptome Assessment of ERBB2-Induced Alterations in Breast Cell Lines

We report the first quantitative and qualitative analysis of the poly (A)(+) transcriptome of two human mammary cell lines, differentially expressing (human epidermal growth factor receptor) an oncogene over-expressed in approximately 25% of human breast tumors. Full-length cDNA populations from the two cell lines were digested enzymatically, individually tagged according to a customized method for library construction, and simultaneously sequenced by the use of the Titanium 454-Roche-platform. Comprehensive bioinformatics analysis followed by experimental validation confirmed novel genes, splicing variants, single nucleotide polymorphisms, and gene fusions indicated by RNA-seq data from both samples. Moreover, comparative analysis showed enrichment in alternative events, especially in the exon usage category, in ERBB2 over-expressing cells, data indicating regulation of alternative splicing mediated by the oncogene. Alterations in expression levels of genes, such as LOX, ATP5L, GALNT3, and MME revealed by large-scale sequencing were confirmed between cell lines as well as in tumor specimens with different ERBB2 backgrounds. This approach was shown to be suitable for structural, quantitative, and qualitative assessment of complex transcriptomes and revealed new events mediated by ERBB2 overexpression, in addition to potential molecular targets for breast cancer that are driven by this oncogene.