Impact of alternative initiation, splicing, and termination on the diversity of the mRNA transcripts encoded by the mouse transcriptome

We analyzed the FANTOM2 clone set of 60,770 RIKEN full-length mouse cDNA sequences and 44,122 public mRNA sequences. We developed a new computational procedure to identify and classify the forms of splice variation evident in this data set and organized the results into a publicly accessible database that can be used for future expression array construction, structural genomics, and analyses of the mechanism and regulation of alternative splicing. Statistical analysis shows that at least 41% and possibly as much as 60% of multiexon genes in mouse have multiple splice forms. Of the transcription units with multiple splice forms, 49% contain transcripts in which the apparent use of an alternative transcription start (stop) is accompanied by alternative splicing of the initial (terminal) exon. This implies that alternative transcription may frequently induce alternative splicing. The fact that 73% of all exons with splice variation fall within the annotated coding region indicates that most splice variation is likely to affect the protein form. Finally, we compared the set of constitutive (present in all transcripts) exons with the set of cryptic (present only in some transcripts) exons and found statistically significant differences in their length distributions, the nucleoticle distributions around their splice junctions, and the frequencies of occurrence of several short sequence motifs.

Molecular cloning of growth hormone receptor (GHR) from common carp (Cyprinus carpio L.) and identification of its two forms of mRNA transcripts

The cDNA of growth hormone receptor (GHR) was cloned from the liver of 2-year common carp (Cyprinus carpio L.) by reverse transcription-polymerase chain reaction (RT-PCR) and rapid amplification of cDNA end (RACE). Its open reading frame (ORF) of 1806 nucleotides is translated into a putative peptide of 602 amino acids, including an extracellular ligand-binding domain of 244 amino acids (aa), a single transmembrane domain of 24 aa and an intracellular signal-transduction domain of 334 aa. Sequence analysis indicated that common carp GHR is highly homologous to goldfish (Carassius auratus) GHR at both gene and protein levels. Using a pair of gene-specific primers, a GHR fragment was amplified from the cDNA of 2-year common carp, a 224 bp product was identified in liver and a 321 bp product in other tissues. The sequencing of the products and the partial genomic DNA indicated that the difference in product size was the result of a 97 bp intron that alternatively spliced. In addition, the 321 bp fragment could be amplified from all the tissues of 4-month common carp including liver, demonstrating the occurrence of the alternative splicing of this intron during the development of common carp. Moreover, a semi-quantitative RT-PCR was performed to analyze the expression level of GHR in tissues of 2-year common carp and 4-month common carp. The result revealed that in the tissues of gill, thymus and brain, the expression level of GHR in 2-year common carp was significantly tower than that of 4-month common carp.

Aberrant mRNA Transcripts and the Nonsense-Mediated Decay Proteins UPF2 and UPF3 Are Enriched in the Arabidopsis Nucleolus

The eukaryotic nucleolus is multifunctional and involved in the metabolism and assembly of many different RNAs and ribonucleoprotein particles as well as in cellular functions, such as cell division and transcriptional silencing in plants. We previously showed that Arabidopsis thaliana exon junction complex proteins associate with the nucleolus, suggesting a role for the nucleolus in mRNA production. Here, we report that the plant nucleolus contains mRNAs, including fully spliced, aberrantly spliced, and single exon gene transcripts. Aberrant mRNAs are much more abundant in nucleolar fractions, while fully spliced products are more abundant in nucleoplasmic fractions. The majority of the aberrant transcripts contain premature termination codons and have characteristics of nonsense-mediated decay (NMD) substrates. A direct link between NMD and the nucleolus is shown by increased levels of the same aberrant transcripts in both the nucleolus and in Up-frameshift (upf) mutants impaired in NMD. In addition, the NMD factors UPF3 and UPF2 localize to the nucleolus, suggesting that the Arabidopsis nucleolus is therefore involved in identifying aberrant mRNAs and NMD.

Granulosal and thecal expression of bone morphogenetic protein- and activin-binding protein mRNA transcripts during bovine follicle development and factors modulating their expression in vitro

Evidence supports local roles for TGFβ superfamily members including activins and bone morphogenetic proteins (BMP) in follicle development. Access of these ligands to signaling receptors is likely modulated by extracellular binding proteins (BP). In this study we compared expression of four BPs (chordin, gremlin, noggin, follistatin) in granulosal (GC) and theca interna (TC) compartments of developing bovine antral follicles (1-18mm). Effects of FSH and IGF on BMP and BP expression by cultured GC, and effects of LH and BMPs on BP expression by cultured TC were also examined. Follicular expression of all four BP transcripts was higher in GC than TC compartments (P<0.001) a finding confirmed by immunohistochemistry. Follicle category affected (P<0.01) gremlin and follistatin mRNA abundance, with a significant cell-type x follicle category interaction for chordin, follistatin and noggin. Noggin transcript abundance was lower (P<0.05) in GC of large 'E-active' than 'E-inactive' follicles while follistatin mRNA level was higher (P<0.01). FSH enhanced CYP19, FSHR, INHBA and follistatin by GC without affecting BMP or BMP-BP expression. IGF increased CYP19 and follistatin, reduced BMP4, noggin and gremlin but did not affect chordin or FSHR mRNA levels. LH increased TC androgen secretion but had no effect on BMP or BP expression. BMPs uniformly suppressed TC androgen production whilst increasing chordin, noggin, and gremlin mRNA levels up to 20-fold (P<0.01). These findings support the hypothesis that extracellular BP, mostly from GC, contribute to the regulation of intrafollicular BMP/activin signaling. Enhancement of thecal BP expression by BMP implies an autoregulatory feedback role to prevent excessive signaling.

Analysis of mRNA transcripts improves the success rate of molecular genetic testing in OTC deficiency

BACKGROUND: Ornithine transcarbamylase (OTC) deficiency is the most common inborn error of urea metabolism that can lead to hyperammonemic crises and orotic aciduria. To date, a total of 341 causative mutations within the OTC gene have been described. However, in about 20% of the patients with enzymatically confirmed OTC deficiency no mutation can be detected when sequencing of genomic DNA analyzing exons and adjacent intronic segments of the OTC gene is performed. METHODS: Standard genomic DNA analysis of the OTC gene in five consecutive patients from five families revealed no mutation. Hence, liver tissue was obtained by needle sampling or open biopsy and RNA extracted from liver was analyzed. RESULTS: Complex rearrangements of the OTC transcript (three insertions and two deletions) were found in all five patients. CONCLUSION: In patients with a strong suspicion of OTC deficiency despite normal results of sequencing exonic regions of the OTC gene, characterization of liver OTC mRNA is highly effective in resolving the genotype. Liver tissue sampling by needle aspiration allows for both enzymatic analysis and RNA based diagnostics of OTC deficiency.

Control of alternative pre-mRNA splicing by distributed pentameric repeats

Multiple copies of the hexamer TGCATG have been shown to regulate fibronectin pre-mRNA alternative splicing. GCATG repeats also are clustered near the regulated calcitonin-specific 3′ splice site in the rat calcitonin/CGRP gene. Specific mutagenesis of these repeats in calcitonin/CGRP pre-mRNA resulted in the loss of calcitonin-specific splicing, suggesting that the native repeats act to enhance alternative exon inclusion. Mutation of subsets of these elements implies that alternative splicing requires a minimum of two repeats, and that the combination of one intronic and one exonic repeat is necessary for optimal cell-specific splicing. However, multimerized intronic repeats inhibited calcitonin-specific splicing in both the wild-type context and in a transcript lacking endogenous repeats. These results suggest that both the number and distribution of repeats may be important features for the regulation of tissue-specific alternative splicing. Further, RNA containing a single repeat bound cell-specific protein complexes, but tissue-specific differences in protein binding were not detected by using multimerized repeats. Together, these data support a novel model for alternative splicing regulation that requires the cell-specific recognition of multiple, distributed sequence elements.

Expression of prostate-specific antigen variant MRNA transcripts in prostate cancer and benign prostatic hyperplasia

Prostate-specific antigen (PSA) is important in tumour detection, monitoring disease progression and tumour recurrence. however, PSA is not a cancerspecific marker as levels can also be elevated in benign prostatic disease. A number of different mRNA transcripts of PSA have also been identified in prostatic tissue, but have not been fully characterized (PSA 424, PSA 525, Schulz transcript). Tissue specimens from transurethral resection of the prostate (TURP) or radical prostatectomy were obtained from 17 men with BPH and 15 men with prostate cancer. Total RNA was extracted, and reverse-transcriptionpolymerase chain reaction (RT-PCR) and Southern analysis carried out using transcript-specific primers and probes to determine which mRNA PSA transcripts were expressed. Real-time PCR was performed to determine transcript levels between the two groups using transcript-specific primers and SYBR green fluorescence. Values obtained were normalized to a standard housekeeping gene, B2-microglobulin. Transcripts amplified by RT-PCR and real-time PCR were confirmed by DNA sequencing. Our results show that the transcripts were present in some, but not all, BPH and cancer samples indicating that they are not specific to either BPH or cancer. Analysis of real-time PCR normalized values using a Student’s t -test, shows that there is a significant difference between the two groups for PSA 424, but not wild-type PSA, PSA 525 or the Schulz transcript. Although a larger cohort of samples is needed to further confirm these results, these findings suggest that mRNA levels of PSA 424 may have some utility as a diagnostic or prognostic marker in prostate cancer detection.

EAAT2 splice variants in Alzheimer disease

Regional atrophy caused by neuronal loss is a characteristic of Alzheimer Disease (AD). Excitatory amino acid transporter-2 (EAAT2) is the major carrier responsible for clearing glutamate from the synaptic cleft in mammalian CNS. A localized attenuation of glutamate transport via reduced expression of functional forms of EAAT2 might contribute to regional excitotoxicity. The EAAT2 gene spans over 100 kb and encodes a 12-kb message. Several groups have identified alternative splice variants of EAAT2 in human brain tissue. These variants can affect transport by altering wild-type EAAT2 protein expression, localization, or transport efficiency. Alternative EAAT2 mRNA transcripts reportedly elicit a dominant-negative effect on glutamate uptake in cell culture. A 50% reduction in the expression in AD cortex of the truncated EAAT2 C-terminal isoform, EAAT2b, has been reported. We obtained cerebral cortex tissue, under informed written consent from the next of kin, from pathologically confirmed control, AD, and non-AD dementia cases. We aimed to determine the distribution and expression patterns of EAAT2 subtypes in susceptible and spared brain regions. We detected five alternate transcripts of EAAT2, two of which had not previously been reported. The relative contributions of novel variants, wild-type EAAT2, and previously discovered splice variants was investigated using Real-time PCR in AD, non-AD dementia, and age-matched control cortex. Our aim is to survey the relationship between these expression patterns and those of markers such as tau, GFAP, and b-amyloid, and to assess the correlation between variant-transporter expression and the level of cell loss.

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.

A Survey of hMLH1 mRNA Splicing Profiles in Human Cell Lines: Comparing Primary Cultured Fibroblasts Before and After Oxidative Stress and Transiently versus Stably Demethylated Colon Cancer Derived Cell Lines

Most human genes undergo alternative splicing and loss of splicing fidelity is associated with disease. Epigenetic silencing of hMLH 1 via promoter cytosine methylation is causally linked to a subset of sporadic non-polyposis colon cancer and is reversible by 5-aza-2' -deoxycytidine treatment. Here I investigated changes in hMLHI mRNA splicing profiles in normal fibroblasts and colon cancer-derived human cell lines. I established the types and frequencies of hMLHI mRNA transcripts generated under baseline conditions, after hydrogen peroxide induced oxidative stress, and in acutely 5-aza-2' -deoxycytidine-treated and stably derepressed cancer cell lines. I found that hMLHI is extensively spliced under all conditions including baseline (50% splice variants), the splice variant distribution changes in response to oxidative stress, and certain splice variants are sensitive to 5- aza-2' -deoxycytidine treatment: Splice variant diversity and frequency of exon 17 skipping correlates with the level of hMLHI promoter methylation suggesting a link between promoter methylation and mRNA splicing.

Heterogeneity in polyadenylation cleavage sites in mammalian mRNA sequences: implications for SAGE analysis

The analysis of a human thyroid serial analysis of gene expression (SAGE) library shows the presence of an abundant SAGE tag corresponding to the mRNA of thyroglobulin (TG). Additional, less abundant tags are present that can not be linked to any other known gene, but show considerable homology to the wild-type TG tag. To determine whether these tags represent TG mRNA molecules with alternative cleavage, 3′-RACE clones were sequenced. The results show that the three putative TG SAGE tags can be attributed to TG transcripts and reflect the use of alternative polyadenylation cleavage sites downstream of a single polyadenylation signal in vivo. By screening more than 300 000 sequences corresponding to human, mouse and rat transcripts for this phenomenon we show that a considerable percentage of mRNA transcripts (44% human, 22% mouse and 22% rat) show cleavage site heterogeneity. When analyzing SAGE-generated expression data, this phenomenon should be considered, since, according to our calculations, 2.8% of human transcripts show two or more different SAGE tags corresponding to a single gene because of alternative cleavage site selection. Both experimental and in silico data show that the selection of the specific cleavage site for poly(A) addition using a given polyadenylation signal is more variable than was previously thought.

Circular RNAs from transcripts of the rat cytochrome P450 2C24 gene: correlation with exon skipping.

The cytochrome P450 2C24 gene is characterized by the capability to generate, in rat kidney, a transcript containing exons 2 and 4 spliced at correct sites but having the donor site of exon 4 directly joined to the acceptor site of exon 2 (exon scrambling). By reverse transcriptase-PCR analysis, it is now shown that the only exons present in the scrambled transcript are exons 2, 3, and 4 and that this molecule lacks a poly(A)+ tail. Furthermore, the use of PCR primers in both orientations of either exon 2 or exon 4 revealed that the orders of the exons in the scrambled transcript are 2-3-4-2 and 4-2-3-4, respectively. These results, combined with the observation that P450 2C24 is a single-copy gene, with no duplication of the exon 2 to exon 4 segment, suggest that the scrambled transcript has properties consistent with that of a circular molecule. In line with this is the observation of an increased resistance of the transcript to phosphodiesterase I, a 3'-exonuclease. Moreover, an alternatively processed cytochrome P450 2C24 mRNA, lacking the three scrambled exons and having exon 1 directly joined to exon 5, has been identified in kidney and liver, tissues that express the scrambled transcript. This complete identity of the exons that are absent in the alternatively processed mRNA but present in the scrambled transcript is interpreted as indicative of the possibility that exon scrambling and exon skipping might be interrelated phenomena. It is therefore proposed that alternative pre-mRNA processing has the potential to generate not only mRNAs lacking one or more exons but also circular RNA molecules.