Heteroduplex formation and S1 digestion for mapping alternative splicing sites

The identification of alternatively spliced transcripts has contributed to a better comprehension of developmental mechanisms, tissue-specific physiological processes and human diseases. Polymerase chain reaction amplification of alternatively spliced variants commonly leads to the formation of heteroduplexes as a result of base pairing involving exons common between the two variants. S1 nuclease cleaves single-stranded loops of heteroduplexes and also nicks the opposite DNA strand. In order to establish a strategy for mapping alternative splice-prone sites in the whole transcriptome, we developed a method combining the formation of heteroduplexes between 2 distinct splicing variants and S1 nuclease digestion. For 20 consensuses identified here using this methodology, 5 revealed a conserved splice site after inspection of the cDNA alignment against the human genome (exact splice sites). For 8 other consensuses, conserved splice sites were mapped at 2 to 30 bp from the border, called proximal splice sites; for the other 7 consensuses, conserved splice sites were mapped at 40 to 800 bp, called distal splice sites. These latter cases showed a nonspecific activity of S1 nuclease in digesting double-strand DNA. From the 20 consensuses identified here, 5 were selected for reverse transcription-polymerase chain reaction validation, confirming the splice sites. These data showed the potential of the strategy in mapping splice sites. However, the lack of specificity of the S1 nuclease enzyme is a significant obstacle that impedes the use of this strategy in large-scale studies.

Alternative splicing enriched cDNA libraries identify breast cancer-associated transcripts

Background: Alternative splicing (AS) is a central mechanism in the generation of genomic complexity and is a major contributor to transcriptome and proteome diversity. Alterations of the splicing process can lead to deregulation of crucial cellular processes and have been associated with a large spectrum of human diseases. Cancer-associated transcripts are potential molecular markers and may contribute to the development of more accurate diagnostic and prognostic methods and also serve as therapeutic targets. Alternative splicing-enriched cDNA libraries have been used to explore the variability generated by alternative splicing. In this study, by combining the use of trapping heteroduplexes and RNA amplification, we developed a powerful approach that enables transcriptome-wide exploration of the AS repertoire for identifying AS variants associated with breast tumor cells modulated by ERBB2 (HER-2/neu) oncogene expression. Results: The human breast cell line (C5.2) and a pool of 5 ERBB2 over-expressing breast tumor samples were used independently for the construction of two AS-enriched libraries. In total, 2,048 partial cDNA sequences were obtained, revealing 214 alternative splicing sequence-enriched tags (ASSETs). A subset with 79 multiple exon ASSETs was compared to public databases and reported 138 different AS events. A high success rate of RT-PCR validation (94.5%) was obtained, and 2 novel AS events were identified. The influence of ERBB2-mediated expression on AS regulation was evaluated by capillary electrophoresis and probe-ligation approaches in two mammary cell lines (Hb4a and C5.2) expressing different levels of ERBB2. The relative expression balance between AS variants from 3 genes was differentially modulated by ERBB2 in this model system. Conclusions: In this study, we presented a method for exploring AS from any RNA source in a transcriptome-wide format, which can be directly easily adapted to next generation sequencers. We identified AS transcripts that were differently modulated by ERBB2-mediated expression and that can be tested as molecular markers for breast cancer. Such a methodology will be useful for completely deciphering the cancer cell transcriptome diversity resulting from AS and for finding more precise molecular markers.

Genetic analysis of complement C1s deficiency associated with systemic lupus erythernatosus highlights alternative splicing of normal C1s gene

Deficiencies of complement proteins of the classical pathway are strongly associated with the development of autoimmune diseases. Deficiency of Clr has been observed to occur concomitantly with deficiency in Cls and 9 out of 15 reported cases presented systemic lupus erythernatosus (SLE). Here, we describe a family in which all four children are deficient in Cls but only two of them developed SLE. Hemolytic activity mediated by the alternative and the lectin pathways were normal, but classical pathway activation was absent in all children`s sera. Cls was undetectable, while in the parents` sera it was lower than in the normal controls. The levels of Clr observed in the siblings and parents sera were lower than in the control, while the concentrations of other complement proteins (C3, C4, MBL and MASP-2) were normal in all family members. Impairment of Cls synthesis was observed in the patients` fibroblasts when analyzed by confocal microscopy. We show that all four siblings are homozygous for a mutation at position 938 in exon 6 of the Cls cDNA that creates a premature stop codon. Our investigations led us to reveal the presence of previously uncharacterized splice variants of Cls mRNA transcripts in normal human cells. These variants are derived from the skipping of exon 3 and from the use of an alternative 3` splice site within intron I which increases the size of exon 2 by 87 nucleotides. (c) 2007 Elsevier Ltd. All rights reserved.

Acute intermittent porphyria: alternative splicing of hydroxymethylbilane synthase mRNA excludes exons 3 and 12

The hydroxymethylbilane synthase (HMBS) mRNAs from 44 control individuals and 30 patients suffering from acute intermittent porphyria (AIP), were screened for length differences by reverse transcriptase polymerase chain reaction (RT-PCR) and any abnormalities were characterized by direct sequencing. Examination of the mRNAs extracted from the peripheral blood lymphocytes of the samples revealed varying degrees of alternative splicing, involving the removal of exons 3 and 12. Approximately 10-50% of the mRNA molecules were affected, despite the absence of genomic splice site mutations or any major deviance from consensus splice sequence values. The preliminary data obtained from this study suggest that this event is a normal occurrence in peripheral blood lymphocytes, and may not be associated with the molecular pathology responsible for AIP. (C) 1998 Academic Press Limited.

Insights into the biological significance of alternative splicing in Arabidopsis: functional characterization of a dual-targeted E3 ligase and the SCL30a SR protein

Dissertation presented to obtain the Ph.D degree in Molecular Biology

RT-PCR diagnosis of patients with acute nonlymphocytic leukemia and inv(16)(p13q22) and identification of new alternative splicing in CBFB-MYH11 transcripts.

As acute nonlymphocytic leukemia (ANLL) with inv(16) (p13q22) or t(16;16)(p13;q22) has been shown to result from the fusion of transcription factor subunit core binding factor (CBFB) to a myosin heavy chain (MYH11), we sought to design methods to detect this rearrangement using reverse transcriptase-polymerase chain reaction (RT-PCR). In all of 27 inv(16)(p13q22) and four t(16;16)(p13;q22) cases tested, a chimeric CBFB-MYH11 transcript coding for an in-frame fusion protein was detected. In a more extensive RT-PCR analysis with different primer pairs, we detected a second new chimeric CBFB-MYH11 transcript in 10 of 11 patients tested. The CBFB-MYH11 reading frame of the second transcript was maintained in one patient but not in the others. We show that the different CBFB-MYH11 transcripts in one patient arise from alternative splicing. Translation of the transcript in which the CBFB-MYH11 reading frame is not maintained leads to a slightly truncated CBFB protein.

PRPF31 alternative splicing and expression in human retina

PURPOSE: To provide a mechanistic link between mutations in PRPF31, and essential and ubiquitously expressed gene, and retinitis pigmentosa, a disorder restricted to the eye. METHODS: We investigated the existence of retina-specific PRPF31 isoforms and the expression of this gene in human retina and other tissues, as well as in cultured human cell lines. PRPF31 transcripts were examined by RT-PCR, quantitative PCR, cloning and sequencing. RESULTS: Database searching revealed the presence of a retina-specific PRPF31 isoform in mouse. However, this isoform could not be experimentally identified in transcripts from human retina or from a human whole eye. Nevertheless, four different PRPF31 isoforms, that were common to all analyzed tissues and cell lines, were isolated. Three of these harbored the full-length PRPF31 coding sequence, whereas the fourth was very short and probably non-coding. The amount of PRPF31 mRNA was previously found to be lower in patients with mutations in this gene than in healthy individuals, making it likely that retinal cells are more sensitive to variation in PRPF31 expression. However, quantitative PCR experiments revealed that PRPF31 mRNA levels in human retina were comparable to those detected in other tissues. CONCLUSIONS: Our results show that the retina-restricted phenotype caused by PRPF31 mutations cannot be explained by the presence of tissue-specific isoforms, or by differential expression of PRPF31 in the retina. As a consequence, the etiology of PRPF31-associated retinitis pigmentosa likely relies on other, probably more subtle molecular mechanisms.

A general definition and nomenclature for alternative splicing events

Understanding the molecular mechanisms responsible for the regulation of the transcriptome present in eukaryotic cells isone of the most challenging tasks in the postgenomic era. In this regard, alternative splicing (AS) is a key phenomenoncontributing to the production of different mature transcripts from the same primary RNA sequence. As a plethora ofdifferent transcript forms is available in databases, a first step to uncover the biology that drives AS is to identify thedifferent types of reflected splicing variation. In this work, we present a general definition of the AS event along with anotation system that involves the relative positions of the splice sites. This nomenclature univocally and dynamically assignsa specific ‘‘AS code’’ to every possible pattern of splicing variation. On the basis of this definition and the correspondingcodes, we have developed a computational tool (AStalavista) that automatically characterizes the complete landscape of ASevents in a given transcript annotation of a genome, thus providing a platform to investigate the transcriptome diversityacross genes, chromosomes, and species. Our analysis reveals that a substantial part—in human more than a quarter—ofthe observed splicing variations are ignored in common classification pipelines. We have used AStalavista to investigate andto compare the AS landscape of different reference annotation sets in human and in other metazoan species and found thatproportions of AS events change substantially depending on the annotation protocol, species-specific attributes, andcoding constraints acting on the transcripts. The AStalavista system therefore provides a general framework to conductspecific studies investigating the occurrence, impact, and regulation of AS.

Deciphering 3'ss selection in the yeast genome reveals an RNA thermosensor that mediates alternative splicing

Poor understanding of the spliceosomal mechanisms to select intronic 3' ends (3'ss) is a major obstacle to deciphering eukaryotic genomes. Here, we discern the rules for global 3'ss selection in yeast. We show that, in contrast to the uniformity of yeast splicing, the spliceosome uses all available 3'ss within a distance window from the intronic branch site (BS), and that in 70% of all possible 3'ss this is likely to be mediated by pre-mRNA structures. Our results reveal that one of these RNA folds acts as an RNA thermosensor, modulating alternative splicing in response to heat shock by controlling alternate 3'ss availability. Thus, our data point to a deeper role for the pre-mRNA in the control of its own fate, and to a simple mechanism for some alternative splicing.

The TPTE gene family: cellular expression, subcellular localization and alternative splicing.

The human TPTE (Transmembrane Phosphatase with TEnsin homology) gene family encodes a PTEN-related tyrosine phosphatase with four potential transmembrane domains. Chromosomal mapping revealed multiple copies of the TPTE gene on chromosomes 13, 15, 21, 22 and Y. Human chromosomes 13 and 21 copies encode two functional proteins, TPIP (TPTE and PTEN homologous Inositol lipid Phosphatase) and TPTE, respectively, whereas only one copy of the gene exists in the mouse genome. In the present study, we show that TPTE and TPIP proteins are expressed in secondary spermatocytes and/or prespermatids. In addition, we report the existence of several novel alternatively spliced isoforms of these two proteins with variable number of transmembrane domains. The latter has no influence on the subcellular localization of these different peptides as shown by co-immunofluorescence experiments. Finally, we identify another expressed TPTE copy, mapping to human chromosome 22, whose transcription appears to be under the control of the LTR of human endogenous retrovirus RTVL-H3.

The implications of alternative splicing in the ENCODE protein complement.

Alternative premessenger RNA splicing enables genes to generate more than one gene product. Splicing events that occur within protein coding regions have the potential to alter the biological function of the expressed protein and even to create new protein functions. Alternative splicing has been suggested as one explanation for the discrepancy between the number of human genes and functional complexity. Here, we carry out a detailed study of the alternatively spliced gene products annotated in the ENCODE pilot project. We find that alternative splicing in human genes is more frequent than has commonly been suggested, and we demonstrate that many of the potential alternative gene products will have markedly different structure and function from their constitutively spliced counterparts. For the vast majority of these alternative isoforms, little evidence exists to suggest they have a role as functional proteins, and it seems unlikely that the spectrum of conventional enzymatic or structural functions can be substantially extended through alternative splicing.

The origins, evolution, and functional potential of alternative splicing in vertebrates.

Alternative splicing (AS) has the potential to greatly expand the functional repertoire of mammalian transcriptomes. However, few variant transcripts have been characterized functionally, making it difficult to assess the contribution of AS to the generation of phenotypic complexity and to study the evolution of splicing patterns. We have compared the AS of 309 protein-coding genes in the human ENCODE pilot regions against their mouse orthologs in unprecedented detail, utilizing traditional transcriptomic and RNAseq data. The conservation status of every transcript has been investigated, and each functionally categorized as coding (separated into coding sequence [CDS] or nonsense-mediated decay [NMD] linked) or noncoding. In total, 36.7% of human and 19.3% of mouse coding transcripts are species specific, and we observe a 3.6 times excess of human NMD transcripts compared with mouse; in contrast to previous studies, the majority of species-specific AS is unlinked to transposable elements. We observe one conserved CDS variant and one conserved NMD variant per 2.3 and 11.4 genes, respectively. Subsequently, we identify and characterize equivalent AS patterns for 22.9% of these CDS or NMD-linked events in nonmammalian vertebrate genomes, and our data indicate that functional NMD-linked AS is more widespread and ancient than previously thought. Furthermore, although we observe an association between conserved AS and elevated sequence conservation, as previously reported, we emphasize that 30% of conserved AS exons display sequence conservation below the average score for constitutive exons. In conclusion, we demonstrate the value of detailed comparative annotation in generating a comprehensive set of AS transcripts, increasing our understanding of AS evolution in vertebrates. Our data supports a model whereby the acquisition of functional AS has occurred throughout vertebrate evolution and is considered alongside amino acid change as a key mechanism in gene evolution.

The evolution of alternative splicing patterns and regulatory mechanisms

Alternative splicing produces multiple isoforms from the same gene, thus increasing the number of transcripts of the species. Alternative splicing is a virtually ubiquitous mechanism in eukaryotes, for example more than 90% of protein-coding genes in human are alternatively spliced. Recent evolutionary studies showed that alternative splicing is a fast evolving and highly species- specific mechanism. The rapid evolution of alternative splicing was considered as a contribution to the phenotypic diversity between species. However, the function of many isoforms produced by alternative splicing remains unclear and they might be the result of noisy splicing. Thus, the functional relevance of alternative splicing and the evolutionary mechanisms of its rapid divergence among species are still poorly understood. During my thesis, I performed a large-scale analysis of the regulatory mechanisms that drive the rapid evolution of alternative splicing. To study the evolution of alternative splicing regulatory mechanisms, I used an extensive RNA-sequencing dataset comprising 12 tetrapod species (human, chimpanzee and bonobo, gorilla, orangutan, macaque, marmoset, mouse, opossum, platypus, chicken and frog) and 8 tissues (cerebellum, brain, heart, kidney, liver, testis, placenta and ovary). To identify the catalogue of alternative splicing eis-acting regulatory elements in the different tetrapod species, I used a previously defined computational approach. This approach is a statistical analysis of exons/introns and splice sites composition and relies on a principle of compensation between splice sites strength and the presence of additional regulators. With an evolutionary comparative analysis of the exonic eis-acting regulators, I showed that these regulatory elements are generally shared among primates and more conserved than non-regulatory elements. In addition, I showed that the usage of these regulatory elements is also more conserved than expected by chance. In addition to the identification of species- specific eis-acting regulators, these results may explain the rapid evolution of alternative splicing. I also developed a new approach based on evolutionary sequence changes and corresponding alternative splicing changes to identify potential splicing eis-acting regulators in primates. The identification of lineage-specific substitutions and corresponding lineage-specific alternative splicing changes, allowed me to annotate the genomic sequences that might have played a role in the alternative splicing pattern differences among primates. Finally, I showed that the identified splicing eis-acting regulator datasets are enriched in human disease-causing mutations, thus confirming their biological relevance.

Alternative splicing of DNA polymerase beta in vertebrates

Alternative splicing (AS) is the predominant mechanism responsible for increasing eukaryotic transcriptome and proteome complexity. In this phenomenon, numerous mRNA transcripts are produced from a single pre-mRNA sequence. AS is reported to occur in 95% of human multi-exon genes; one specific gene that undergoes AS is DNA polymerase beta (POLB). POLB is the main DNA repair gene which performs short patch base excision repair (BER). In primate untransformed primary fibroblast cell lines, it was determined that the splice variant (SV) frequency of POLB correlates positively with species lifespan. To date, AS patterns of POLB have only been examined in mammals primarily through the use of cell lines. However, little attention has been devoted to investigating if such a relationship exists in non-mammals and whether cell lines reflect what is observed in vertebrate tissues. This idea was explored through cloning and characterization of 1,214 POLB transcripts from four non-mammalian species (Gallus gallus domesticus, Larus glaucescens, Xenopus laevis, and Pogona vitticeps) and two mammalian species (Sylvilagus floridanus and Homo sapiens) in two tissue types, liver and brain. POLB SV frequency occurred at low frequencies, < 3.2%, in non-mammalian tissues relative to mammalian (>20%). The highest POLB SV frequency was found in H. sapiens liver and brain tissues, occurring at 65.4% and 91.7%, respectively. Tissue specific AS of POLB was observed in L. glaucescens, P. vitticeps, and H. sapiens, but not G. gallus domesticus, X. laevis and S. floridanus.The AS patterns of a second gene, transient receptor potential cation channel subfamily V member 1 (TRPV1), were compared to those of POLB in liver and brain tissues of G. gallus domesticus, X. laevis and H. sapiens. This comparison was performed to investigate if any changes (either increase or decrease) observed in the AS of POLB were gene specific or if they were tissue specific, in which case similar changes in AS would be seen in POLB and TRPV1. Analysis did not reveal an increase or decrease in both the AS of POLB and TRPV1 in either the liver or brain tissues of G. gallus domesticus and H. sapiens. This result suggested that the AS patterns of POLB were not influenced by tissue specific rates of AS. Interestingly, an increase in the AS of both genes was only observed in X. laevis brain tissue. This result suggests that AS in general may be increased in the X. laevis brain as compared to liver tissue. No positive correlation between POLB SV frequency and species lifespan was found in non-mammalian tissues. The AS patterns of POLB in human primary untransformed fibroblast cell lines were representative of those seen in human liver tissue but not in brain tissue. Altogether, the AS patterns of POLB from vertebrate tissues and primate cell lines revealed a positive correlation between POLB SV frequency and lifespan in mammals, but not in non-mammals. It appears that this positive correlation does not exist in vertebrate species as a whole.

Alternative splicing by hnRNP L as a new modulator of hematopoietic cell differentiation, survival and migration

Les modifications post-transcriptionnelles de l’ARN messager (ARNm), comme l’épissage alternatif, jouent un rôle important dans la régulation du développement embryonnaire, de la fonction cellulaire et de l’immunité. De nouvelles évidences révèlent que l’épissage alternatif serait également impliqué dans la régulation de la maturation et de l’activation des cellules du système hématopoïétique. Le facteur hnRNP L a été identifié comme étant le principal régulateur de l’épissage alternatif du gène codant pour le récepteur CD45 in vitro. Le récepteur CD45 est une tyrosine phosphatase exprimée par toutes les cellules du système hématopoïétique qui contrôle le développement et l’activation des lymphocytes T. Dans un premier temps, nous avons étudié la fonction du facteur hnRNP L dans le développement des lymphocytes T et dans l’épissage de l’ARNm de CD45 in vivo en utilisant des souris dont le gène de hnRNP L a été supprimé spécifiquement dans les cellules T. La délétion de hnRNP L dans les thymocytes résulte en une expression aberrante des différents isoformes de CD45 avec une prédominance de l'isoforme CD45RA qui est généralement absent dans le thymus. Une conséquence de la délétion de hnRNP L est une diminution de la cellularité du thymus causée par un blocage partiel du développement des cellules pré-T au stade DN4. Cette réduction du nombre de cellules dans le thymus n’est pas liée à une hausse de la mort cellulaire. Les thymocytes déficients pour hnRNP L démontrent plutôt une prolifération augmentée comparée aux thymocytes sauvages due à une hyper-activation des kinases Lck, Erk1/2 et Akt. De plus, la délétion de hnRNP L dans le thymus cause une perte des cellules T en périphérie. Les résultats des expériences in vitro suggèrent que cette perte est principalement due à un défaut de migration des thymocytes déficients pour hnRNP L du thymus vers la périphérie en réponse aux chimiokines. L’épissage alternatif de CD45 ne peut expliquer ce phénotype mais l’identification de cibles par RNA-Seq a révélé un rôle de hnRNP L dans la régulation de l’épissage alternatif de facteurs impliqués dans la polymérisation de l’actine. Dans un second temps, nous avons étudié le rôle de hnRNP L dans l’hématopoïèse en utilisant des souris dont la délétion de hnRNP L était spécifique aux cellules hématopoïétiques dans les foies fœtaux et la moelle osseuse. L’ablation de hnRNP L réduit le nombre de cellules progénitrices incluant les cellules progénitrices lymphocytaires (CLPs), myéloïdes (CMPs, GMPs) et mégakaryocytes-érythrocytaires (MEPs) et une perte des cellules hématopoïétiques matures. À l’opposé des cellules progénitrices multipotentes (MPPs) qui sont affectées en absence de hnRNP L, la population de cellules souches hématopoïétiques (HSCs) n’est pas réduite et prolifère plus que les cellules contrôles. Cependant, les HSCs n’exprimant pas hnRNP L sont positives pour l'Annexin V et expriment CD95 ce qui suggère une mort cellulaire prononcée. Comme pour les thymocytes, une analyse par RNA-Seq des foies fœtaux a révélé différents gènes cibles de hnRNP L appartenant aux catégories reliées à la mort cellulaire, la réponse aux dommages à l’ADN et à l’adhésion cellulaire qui peuvent tous expliquer le phénotype des cellules n’exprimant pas le gène hnRNP L. Ces résultats suggèrent que hnRNP L et l’épissage alternatif sont essentiels pour maintenir le potentiel de différenciation des cellules souches hématopoïétiques et leur intégrité fonctionnelle. HnRNP L est aussi crucial pour le développement des cellules T par la régulation de l’épissage de CD45 ainsi que pour leur migration.