Alternative splicing of human genes: more the rule than the exception?

Molecular Basis of Influence of A501X Mutations in Penicillin-Binding Protein 2 of Neisseria gonorrhoeae Strain 35/02 on Ceftriaxone Resistance

The increase in the resistance of mutant strains of Neisseria gonorrhoeae to the antibiotic ceftriaxone is pronounced in the decrease in the second-order acylation rate constant, k2/KS, by penicillin-binding protein 2 (PBP2). These changes can be caused by both the decrease in the acylation rate constant, k2, and the weakening of the binding affinity, i.e., an increase in the substrate constant, KS. A501X mutations in PBP2 affect second-order acylation rate constants. The PBP2A501V variant exhibits a higher k2/KS value, whereas for PBP2A501R and PBP2A501P variants, these values are lower. We performed molecular dynamic simulations with both classical and QM/MM potentials to model both acylation energy profiles and conformational dynamics of four PBP2 variants to explain the origin of k2/KS changes. The acylation reaction occurs in two elementary steps, specifically, a nucleophilic attack by the oxygen atom of the Ser310 residue and C-N bond cleavage in the β-lactam ring accompanied by the elimination of the leaving group of ceftriaxone. The energy barrier of the first step increases for PBP2 variants with a decrease in the observed k2/KS value. Submicrosecond classic molecular dynamic trajectories with subsequent cluster analysis reveal that the conformation of the β3-β4 loop switches from open to closed and its flexibility decreases for PBP2 variants with a lower k2/KS value. Thus, the experimentally observed decrease in the k2/KS in A501X variants of PBP2 occurs due to both the decrease in the acylation rate constant, k2, and the increase in KS.

Alternative splicing isoform in succinate dehydrogenase complex, subunit C causes downregulation of succinate-coenzyme Q oxidoreductase activity in mitochondria

Mitochondrial succinate dehydrogenase (SDH) is localized to the inner mitochondrial membrane and is responsible for the redox of succinic acid. SDH is a tetrameric iron-sulfur flavoprotein of the tricarboxylic acid cycle and respiratory chain. The SDH complex, subunit C (SDHC) transcript has deletion-type alternative splicing sites. Generally, alternative splicing produces variant proteins and expression patterns, as products of different genes. In certain cases, specific alternative splicing variants (ASVs) have been associated with human disease. Due to a frameshift mutation causing loss of the heme binding region, the SDHC Δ5 isoform (lacking exon 5) exhibits no SDHC activity. To investigate whether the SDHC splicing variants can function as dominant-negative inhibitors, SDHC ASVs were overexpressed in HCT-15 human colorectal cancer cells. Using real-time reverse transcription-polymerase chain reaction, a dominant-negative effect of the Δ5 isoform on SDHC mRNA was shown. In addition, Δ5 overexpression increased the levels of reactive oxygen species. Furthermore, in the Δ5 isoform-overexpressing cells, SDH activity was reduced. SDHC activation is a significant event during the electron transport chain, and the function of the SDHC Δ5 variant may be significant for the differentiation of tumor cells.

An intron-retaining splice variant of human cyclin A2, expressed in adult differentiated tissues, induces a G1/S cell cycle arrest in vitro

Background

Human cyclin A2 is a key regulator of S phase progression and entry into mitosis. Alternative splice variants of the G1 and mitotic cyclins have been shown to interfere with full-length cyclin functions to modulate cell cycle progression and are therefore likely to play a role in differentiation or oncogenesis. The alternative splicing of human cyclin A2 has not yet been studied.

Methodology/Principal Findings

Sequence-specific primers were designed to amplify various exon–intron regions of cyclin A2 mRNA in cell lines and human tissues. Intron retaining PCR products were cloned and sequenced and then overexpressed in HeLa cells. The subcellular localization of the splice variants was studied using confocal and time-lapse microscopy, and their impact on the cell cycle by flow cytometry, immunoblotting and histone H1 kinase activity. We found a splice variant of cyclin A2 mRNA called A2V6 that partly retains Intron 6. The gene expression pattern of A2V6 mRNA in human tissues was noticeably different from that of wild-type cyclin A2 (A2WT) mRNA. It was lower in proliferating fetal tissues and stronger in some differentiated adult tissues, especially, heart. In transfected HeLa cells, A2V6 localized exclusively in the cytoplasm whereas A2WT accumulated in the nucleus. We show that A2V6 induced a clear G1/S cell cycle arrest associated with a p21 and p27 upregulation and an inhibition of retinoblastoma protein phosphorylation. Like A2WT, A2V6 bound CDK2, but the A2V6/CDK2 complex did not phosphorylate histone H1.

Conclusion/Significance

This study has revealed that some highly differentiated human tissues express an intron-retaining cyclin A2 mRNA that induced a G1/S block in vitro. Contrary to full-length cyclin A2, which regulates cell proliferation, the A2V6 splice variant might play a role in regulating nondividing cell states such as terminal differentiation or senescence.

BEAT: Bioinformatics Exon Array Tool to store, analyze and visualize Affymetrix GeneChip Human Exon Array data from disease experiments

BACKGROUND

It is known from recent studies that more than 90% of human multi-exon genes are subject to Alternative Splicing (AS), a key molecular mechanism in which multiple transcripts may be generated from a single gene. It is widely recognized that a breakdown in AS mechanisms plays an important role in cellular differentiation and pathologies. Polymerase Chain Reactions, microarrays and sequencing technologies have been applied to the study of transcript diversity arising from alternative expression. Last generation Affymetrix GeneChip Human Exon 1.0 ST Arrays offer a more detailed view of the gene expression profile providing information on the AS patterns. The exon array technology, with more than five million data points, can detect approximately one million exons, and it allows performing analyses at both gene and exon level. In this paper we describe BEAT, an integrated user-friendly bioinformatics framework to store, analyze and visualize exon arrays datasets. It combines a data warehouse approach with some rigorous statistical methods for assessing the AS of genes involved in diseases. Meta statistics are proposed as a novel approach to explore the analysis results. BEAT is available at http://beat.ba.itb.cnr.it.

RESULTS

BEAT is a web tool which allows uploading and analyzing exon array datasets using standard statistical methods and an easy-to-use graphical web front-end. BEAT has been tested on a dataset with 173 samples and tuned using new datasets of exon array experiments from 28 colorectal cancer and 26 renal cell cancer samples produced at the Medical Genetics Unit of IRCCS Casa Sollievo della Sofferenza.To highlight all possible AS events, alternative names, accession Ids, Gene Ontology terms and biochemical pathways annotations are integrated with exon and gene level expression plots. The user can customize the results choosing custom thresholds for the statistical parameters and exploiting the available clinical data of the samples for a multivariate AS analysis.

CONCLUSIONS

Despite exon array chips being widely used for transcriptomics studies, there is a lack of analysis tools offering advanced statistical features and requiring no programming knowledge. BEAT provides a user-friendly platform for a comprehensive study of AS events in human diseases, displaying the analysis results with easily interpretable and interactive tables and graphics.

COMPARATIVE ANALYSIS OF DISEASE-ASSOCIATED MUTATIONS IN THE CODING REGIONS OF ALTERNATIVELY AND CONSTITUTIVELY SPLICED HUMAN GENES

To better understand the evolutionary and molecular mechanisms of alternative splicing causing human diseases, we have systematically compared the pattern, the distribution and the density of disease-associated mutations as well as the influence of codon usage bias on the single mutation between alternatively and constitutively spliced genes through analysis of the large datasets from human disease genes. The results indicated that:

1. The most common pattern of single mutation in alternatively and constitutively spliced genes are, respectively, C/T (25.17%), (22.81%) and G/A (21.54%), (22.73%), suggesting that the two types of disease genes are prone to C → T and G → A mutations.

2. There is an overall preponderance for transitions over transversions in alternatively (62.88% versus 37.12%) and constitutively (64.41% versus 35.59%) spliced disease genes.

3. For the second base of codons, there exist significant differences in transitions and transversions between the two types of genes.

4. Our data indicated that the single mutation tends to occur preferentially when the upstream neighboring-nucleotide is C or G in human disease genes.

5. Codon usage bias and synonymous codon usage have great influence on the single mutation in both alternatively and constitutively spliced genes. The GC content and gene length also have very evident influence on such mutations.

Our results seem to imply that disease-associated mutations within the coding regions of alternatively spliced human disease genes have different mechanisms from constitutively spliced genes. Such findings may facilitate understanding the molecular mechanism of alternative splicing causing human diseases, and the development of gene therapies for such diseases.

Alternative RNA splicing in expression of the glutathione synthetase gene in human cells

Ubiquitous free radical production occurs continuously in cells and tissues. Glutathione is the most abundant mammalian antioxidant, and is synthesized by glutathione synthetase (GSS). Therefore, GSS plays an important role in defending the cell against reactive oxygen species. The expression of GSS has been studied in human cells; however, sequence information about alternative splicing variants of GSS mRNA has not been reported. In the present study, we identified a novel alternative splicing variant (ASV) of the GSS gene in 10 human normal tissues and five human cancer cell lines. The deleted transcript of GSS was characterized by an in-frame deletion of 333 bp, corresponding to the complete loss of exons 4 and 5. Thus this GSS ASV causes protein truncation. We quantified the mRNA of GSS ASV in human normal tissues using real-time PCR. The ASV was detected in colon, kidney, lung, liver, placenta, peripheral blood and uterus, but not in heart, skeletal muscle and spleen tissue. Our results provide a basis for more detailed studies on the regulation of GSS, and for further evaluation of this and other possible roles of GSS. Understanding the regulation of GSS expression is very important for the development of new strategies for controlling the development of GSH-based redox homeostasis.

Expressed sequence tags (ESTs) and single nucleotide polymorphisms (SNPs): What large-scale sequencing projects can tell us about ADME

To date over 800 complete genomes have been sequenced, with many more partially complete. Coupled with the large amount of mRNA transcript sequence data being produced from expression studies, there is now a daunting amount of information available to the research scientist. This review examines how this information may be best used, focusing on examples from sequences encoding absorption, distribution, metabolism and excretion (ADME)-related proteins in particular. Through the use of phylogenetic, splice variant and single nucleotide polymorphism (SNP) analysis, the review examines not only how insights into species-specific responses to drug exposure may be gained, but also how best to utilize this information to predict both individual human responses and the impact of population variance in response.

SAF-3, a novel splice variant of the SAF-1/MAZ/Pur-1 family, is expressed during inflammation

The Cys2His2-type zinc finger transcription factor serum amyloid A activating factor 1 [SAF-1, also known as MAZ (myc-associated zinc finger protein) or Pur-1 (purine binding factor-1)] plays an important role in regulation of a variety of inflammation-responsive genes. An SAF-2 splice variant acting as a negative regulator of SAF-1 was identified previously, and the present study reports the identification of a novel SAF-3 splice variant that is expressed during inflammation. SAF-3 mRNA, isolated from a cDNA library produced from IL-1beta-induced cells, originates from a previously unknown first coding exon, and thereby contains a unique N-terminal domain but shares the same six zinc finger DNA-binding domains as present in SAF-1. In addition, a negatively functioning domain present at the N-terminus of SAF-1 and SAF-2 is spliced out in SAF-3. The expression of SAF-3 is very low in normal tissues and in cells grown under normal conditions. However, RT-PCR analysis of mRNAs from cytokine and growth factor-induced cells as well of mRNAs isolated from several diseased tissues revealed abundant expression of SAF-3. The transactivation potential of SAF-3 is much greater than that of the predominantly expressed splice variant SAF-1. These findings show that transcriptional regulation of downstream inflammation-responsive genes by SAF/MAZ/Pur-1 is likely to be more complex than previously assumed. In addition, we show that SAF-3 expression initiates from an upstream novel promoter. This is the first report of the existence of multiple promoters regulating expression of the SAF/MAZ/Pur-1 family of proteins.

BimL upregulation induced by BCR cross-linking in BL41 Burkitt's lymphoma results from a splicing mechanism of the BimEL mRNA

B lymphocyte receptor-mediated apoptosis is associated with increased expression of the BimL isoform of Bim. The mechanisms involved in the regulation of BimL protein expression are still unknown. We report that BimL expression following BCR activation is not associated with a specific increase of BimL mRNA but rather to the intron retention structure of the BimEL mRNA. Indeed, expression of a BimEL cDNA leads in Hela cells leads to the production of both BimEL and BimL proteins. Mutation of the intron-splicing GT sequence present in the exon 3 results in the production of only BimEL protein. Ectopic expression of BimEL cDNA resulted in a large increase of BimL expression upon BCR-stimulation, whereas cells transfected with the GT/AA mutated form of BimEL only produced BimEL proteins upon BCR-activation. These data showed that BimL expression induced by BCR activation may result from the splicing of BimEL mRNA independently of Bim promoter regulation.

ING1a expression increases during replicative senescence and induces a senescent phenotype

The ING family of tumor suppressor proteins affects cell growth, apoptosis and response to DNA damage by modulating chromatin structure through association with different HAT and HDAC complexes. The major splicing isoforms of the ING1 locus are ING1a and INGlb. While INGlb plays a role in inducing apoptosis, the function of ING1a is currently unknown. Here we show that alternative splicing of the ING1 message alters the INGla:INGlb ratio by approximately 30-fold in senescent compared to low passage primary fibroblasts. INGla antagonizes INGlb function in apoptosis, induces the formation of structures resembling senescence-associated heterochromatic foci containing heterochromatin protein 1 gamma, the accumulation of senescence-associated beta-galactosidase activity and promotes senescent cell morphology and cell cycle arrest. Phenotypic effects may result from differential effects on gene expression since ING1a increases levels of both retinoblastoma and the p16 cyclin-dependent kinase inhibitor and ING1a and ING1b have opposite effects on the expression of proliferating nuclear cell antigen (PCNA), which is required for cell growth. Gene expression appears to be altered by targeting of HDAC complexes to gene promoters since INGla associates with several-fold higher levels of HDAC1 in senescent, compared to replication-competent cells and ING1 is found on the PCNA promoter by chromatin immunoprecipitation analysis. These data demonstrate a novel role for the ING1 proteins in differentially regulating senescence-associated chromatin remodeling vs. apoptosis and support the idea that altered ratios of the ING1 splicing isoforms may contribute to establishing the senescent phenotype through HDAC and HAT complex-mediated effects on chromatin structure.

Regulation of expression of two LY-6 family genes by intron retention and transcription induced chimerism

Background

Regulation of the expression of particular genes can rely on mechanisms that are different from classical transcriptional and translational control. The LY6G5B and LY6G6D genes encode LY-6 domain proteins, whose expression seems to be regulated in an original fashion, consisting of an intron retention event which generates, through an early premature stop codon, a non-coding transcript, preventing expression in most cell lines and tissues.

Results

The MHC LY-6 non-coding transcripts have shown to be stable and very abundant in the cell, and not subject to Nonsense Mediated Decay (NMD). This retention event appears not to be solely dependent on intron features, because in the case of LY6G5B, when the intron is inserted in the artificial context of a luciferase expression plasmid, it is fully spliced but strongly stabilises the resulting luciferase transcript. In addition, by quantitative PCR we found that the retained and spliced forms are differentially expressed in tissues indicating an active regulation of the non-coding transcript. EST database analysis revealed that these genes have an alternative expression pathway with the formation of Transcription Induced Chimeras (TIC). This data was confirmed by RT-PCR, revealing the presence of different transcripts that would encode the chimeric proteins CSNKβ-LY6G5B and G6F-LY6G6D, in which the LY-6 domain would join to a kinase domain and an Ig-like domain, respectively.

Conclusion

In conclusion, the LY6G5B and LY6G6D intron-retained transcripts are not subjected to NMD and are more abundant than the properly spliced forms. In addition, these genes form chimeric transcripts with their neighbouring same orientation 5' genes. Of interest is the fact that the 5' genes (CSNKβ or G6F) undergo differential splicing only in the context of the chimera (CSNKβ-LY6G5B or G6F-LY6G6C) and not on their own.

Alternative splicing and expression profile analysis of expressed sequence tags in domestic pig

Domestic pig (Sus scrofa domestica) is one of the most important mammals to humans. Alternative splicing is a cellular mechanism in eukaryotes that greatly increases the diversity of gene products. Expression sequence tags (ESTs) have been widely used for gene discovery, expression profile analysis, and alternative splicing detection. In this study, a total of 712,905 ESTs extracted from 101 different non-normalized EST libraries of the domestic pig were analyzed. These EST libraries cover the nervous system, digestive system, immune system, and meat production related tissues from embryo, newborn, and adult pigs, making contributions to the analysis of alternative splicing variants as well as expression profiles in various stages of tissues. A modified approach was designed to cluster and assemble large EST datasets, aiming to detect alternative splicing together with EST abundance of each splicing variant. Much efforts were made to classify alternative splicing into different types and apply different filters to each type to get more reliable results. Finally, a total of 1,223 genes with average 2.8 splicing variants were detected among 16,540 unique genes. The overview of expression profiles would change when we take alternative splicing into account.

The rice Mybleu transcription factor increases tolerance to oxygen deprivation in Arabidopsis plants

Mybleu is a natural incomplete transcription factor of rice (Oryza sativa), consisting of a partial Myb repeat followed by a short leucine zipper. We previously showed its localization to the apical region of rice roots and coleoptiles. Specifically, in coleoptiles, Mybleu is expressed under both aerobic and anaerobic conditions, whereas in roots, it is expressed only under aerobic conditions. Mybleu is able to dimerize with canonical leucine zippers and to activate transcription selectively. To investigate Mybleu function in vivo, we transformed Arabidopsis thaliana and evaluated several morphological, physiological and biochemical parameters. In agreement with a hypothesized role of Mybleu in cell elongation in the differentiation zone, we found that the constitutive expression of this transcription factor in Arabidopsis induced elongation in the primary roots and in the internodal region of the floral stem; we also observed a modification of the root apex morphology in transformed lines. Based on the high expression of Mybleu in anaerobic rice coleoptiles, we studied the role of this transcription factor in transgenic plants grown under low-oxygen conditions. We found that overexpression of this transcription factor increased tolerance to oxygen deficit. In transgenic plants, this effect may depend both on the maintenance of a higher metabolism during stress and on the higher expression levels of certain genes involved in the anaerobic response.

Analysis of chameleon sequences and their implications in biological processes

Chameleon sequences have been implicated in amyloid related diseases. Here we report an analysis of two types of chameleon sequences, chameleon-HS (Helix vs. Strand) and chameleon-HE (Helix vs. Sheet), based on known structures in Protein Data Bank. Our survey shows that the longest chameleon-HS is eight residues while the longest chameleon-HE is seven residues. We have done a detailed analysis on the local and global environment that might contribute to the unique conformation of a chameleon sequence. We found that the existence of chameleon sequences does not present a problem for secondary structure prediction programs, including the first generation prediction programs, such as Chou-Fasman algorithm, and the third generation prediction programs that utilize evolution information. We have also investigated the possible implication of chameleon sequences in structural conservation and functional diversity of alternatively spliced protein isoforms.

Alternative splicing in cancer: Noise, functional, or systematic?

Pre-messenger RNA splicing is a fine-tuned process that generates multiple functional variants from individual genes. Various cell types and developmental stages regulate alternative splicing patterns differently in their generation of specific gene functions. In cancers, splicing is significantly altered, and understanding the underlying mechanisms and patterns in cancer will shed new light onto cancer biology. Cancer-specific transcript variants are promising biomarkers and targets for diagnostic, prognostic, and treatment purposes. In this review, we explore how alternative splicing cannot simply be considered as noise or an innocent bystander, but is actively regulated or deregulated in cancers. A special focus will be on aspects of cell biology and biochemistry of alternative splicing in cancer cells, addressing differences in splicing mechanisms between normal and malignant cells. The systems biology of splicing is only now applied to the field of cancer research. We explore functional annotations for some of the most intensely spliced gene classes, and provide a literature mining and clustering that reflects the most intensely investigated genes. A few well-established cancer-specific splice events, such as the CD44 antigen, are used to illustrate the potential behind the exploration of the mechanisms of their regulation. Accordingly, we describe the functional connection between the regulatory machinery (i.e., the spliceosome and its accessory proteins) and their global impact on qualitative transcript variation that are only now emerging from the use of genomic technologies such as microarrays. These studies are expected to open an entirely new level of genetic information that is currently still poorly understood.

Compact genes are highly expressed in the moss Physcomitrella patens

Compact genes contain short and few introns, and they are highly expressed in different animal genomes. Recently, it has been shown that in Oryza sativa and Arabidopsis thaliana, highly expressed genes tend to be least compact, containing long and many introns. It has been suggested that selection on genome organization may have acted differently in plants compared with animals. Gene expression can be estimated as the number of hits when comparing a gene sequence with publicly available expressed sequence tags. Here it is shown that in the haploid moss Physcomitrella pates, highly expressed genes contain shorter introns than genes with low expression levels. This study therefore supports the hypothesis that selection may strongly favour transcriptional efficiency at least in the haploid phase of plant life cycles. It is concluded that plants do not necessarily respond to other selection pressures than animals regarding genome structuring.

A subgroup of MYB transcription factor genes undergoes highly conserved alternative splicing in Arabidopsis and rice

MYB transcription factor genes play important roles in many developmental processes and in various defence responses of plants. Two Arabidopsis R2R3-type MYB genes, AtMYB59 and AtMYB48, were found to undergo similar alternative splicing. Both genes have four distinctively spliced transcripts that encode either MYB-related proteins or R2R3-MYB proteins. An extensive BLAST search of the GenBank database resulted in finding and cloning two rice homologues, both of which were also found to share a similar alternative splicing pattern. In a semi-quantitative study, the expression of one splice variant of AtMYB59 was found to be differentially regulated in treatments with different phytohormones and stresses. GFP fusion protein analysis revealed that both of the two predicted nuclear localization signals (NLSs) in the R3 domain are required for localizing to the nucleus. Promoter-GUS analysis in transgenic plants showed that 5'-UTR is sufficient for the translation initiation of type 3 transcripts (encoding R2R3-MYB proteins), but not for type 2 transcripts (encoding MYB-related proteins). Moreover, a new type of non-canonical intron, with the same nucleotide repeats at the 5' and 3' splice sites, was identified. Thirty-eight Arabidopsis and rice genes were found to have this type of non-canonical intron, most of which undergo alternative splicing. These data suggest that this subgroup of transcription factor genes may be involved in multiple biological processes and may be transcriptionally regulated by alternative splicing.

Alternatively and constitutively spliced exons are subject to different evolutionary forces

There has been a controversy on whether alternatively spliced exons (ASEs) evolve faster than constitutively spliced exons (CSEs). Although it has been noted that ASEs are subject to weaker selective constraints than CSEs, so they evolve faster, there have also been studies that indicated slower evolution in ASEs than in CSEs. In this study, we retrieve more than 5,000 human-mouse orthologous exons and calculate the synonymous (KS) and nonsynonymous (KA) substitution rates in these exons. Our results show that ASEs have higher KA values and higher KA/KS ratios than CSEs, indicating faster amino acid-level evolution in ASEs. The faster evolution may be in part due to weaker selective constraints. It is also possible that the faster rate is in part due to faster functional evolution in ASEs. On the other hand, the majority of ASEs have lower KS values than CSEs. With reference to the substitution rate in introns, we show that the KS values in ASEs are close to the neutral substitution rate, whereas the synonymous substitution rate in CSEs has likely been accelerated. The elevated synonymous rate in CSEs is not related to CpG dinucleotides or low-complexity regions of protein but may be weakly related to codon usage bias. The overall trends of higher KA and lower KS in ASEs than in CSEs are also observed in human-rat and mouse-rat comparisons. Therefore, our observations hold for mammals of different molecular clocks.

High frequency of alternative first exons in erythroid genes suggests a critical role in regulating gene function

The human genome uses alternative pre-mRNA splicing as an important mechanism to encode a complex proteome from a relatively small number of genes. An unknown number of these genes also possess multiple transcriptional promoters and alternative first exons that contribute another layer of complexity to gene expression mechanisms. Using a collection of more than 100 erythroid-expressed genes as a test group, we used genome browser tools and genetic databases to assess the frequency of alternative first exons in the genome. Remarkably, 35% of these erythroid genes show evidence of alternative first exons. The majority of the candidate first exons are situated upstream of the coding exons, whereas a few are located internally within the gene. Computational analyses predict transcriptional promoters closely associated with many of the candidate first exons, supporting their authenticity. Importantly, the frequent presence of consensus translation initiation sites among the alternative first exons suggests that many proteins have alternative N-terminal structures whose expression can be coupled to promoter choice. These findings indicate that alternative promoters and first exons are more widespread in the human genome than previously appreciated and that they may play a major role in regulating expression of selected protein isoforms in a tissue-specific manner.

Evolutionarily conserved coupling of transcription and alternative splicing in the EPB41 (protein 4.1R) and EPB41L3 (protein 4.1B) genes

Recent studies have shown that transcription and alternative splicing can be mechanistically coupled. In the EPB41 (protein 4.1R) and EPB41L3 (protein 4.1B) genes, we showed previously that promoter/alternative first exon choice is coupled to downstream splicing events in exon 2. Here we demonstrate that this coupling is conserved among several vertebrate classes from fish to mammals. The EPB41 and EPB41L3 genes from fish, bird, amphibian, and mammal genomes exhibit shared features including alternative first exons and differential splice acceptors in exon 2. In all cases, the 5'-most exon (exon 1A) splices exclusively to a weaker internal acceptor site in exon 2, skipping a fragment designated as exon 2'. Conversely, alternative first exons 1B and 1C always splice to the stronger first acceptor site, retaining exon 2'. These correlations are independent of cell type or species of origin. Since exon 2' contains a translation initiation site, splice variants generate protein isoforms with distinct N-termini. We propose that these genes represent a physiologically relevant model system for mechanistic analysis of transcription-coupled alternative splicing.

Structural and functional diversity generated by alternative mRNA splicing

Alternative mRNA splicing is becoming increasingly recognized as an important mechanism for the generation of structural and functional diversity in proteins. Recent estimations predict that approximately 50% of all eukaryotic proteins can be alternatively spliced. Several lines of evidence suggest that alternative mRNA splicing results in small changes in protein structure and is likely to fine-tune the function and specificity of the affected protein. However, knowledge of how alternative splicing regulates cellular processes on the molecular level is still limited. It is only recently that structures of alternatively spliced proteins have been solved. These studies have shown that alternative splicing affects the structure not only in the vicinity of the splice site but also at long distance.

[Alternative mRNA splicing, pathology and molecular therapeutics]

Pre-mRNA splicing operates towards at least 95 % of the transcript pool. It is subjected to a large number of variations, collectively regrouped under the term of alternative mRNA splicing, which occurs, on average, 6 to 8 times per pre-mRNA molecule. Consequently, many more proteins may be encoded from a single gene, which may satisfy a physiological need, or mark a pathological adaptation. The identification of mutations in sequences required for splicing, both constitutive and alternative, or for their control, has permitted to determine the causes of qualitative or quantitative variations in transcript levels associated with inherited diseases or cancer development. A number of molecular approaches have been undertaken to try to compensate for the effect of deleterious splicing mutations and to restore, at least in part, sufficient amounts of either the normal or a surrogate transcript. These include overexpression of splicing proteins, improvement of their activity by post-translational modification, splice-site increased or decreased usage, and RNA-mediated trans-splicing. Using such approaches, phenotypic improvements have been obtained in animal models, carrying new hopes for the development of therapeutic strategies aimed at correcting both inherited and acquired diseases that involve pre-mRNA splicing defects.

Expression of c-Fes Protein Isoforms Correlates with Differentiation in Myeloid Leukemias

The cellular fes gene encodes a 93-kilodalton protein-tyrosine kinase (p93) that is expressed in both normal and neoplastic myeloid cells. Increased c-Fes expression is associated with differentiation in normal myeloid cells and cell lines. Our hypothesis was that primary leukemia cells would show a similar pattern of increased expression in more differentiated cells. Therefore, we compared c-Fes expression in cells with an undifferentiated, blast phenotype (acute myelogenous leukemia--AML) to cells with a differentiated phenotype (chronic myelogenous leukemia--CML). Instead of differences in p93 expression levels, we found complex patterns of c-Fes immunoreactive proteins that corresponded with differentiation in normal and leukemic myeloid cells. The "blast" pattern consisted of c-Fes immunoreactive proteins p93, p74, and p70; the "differentiated" pattern showed two additional c-Fes immunoreactive proteins, p67 and p62. Using mRNA from mouse and human cell lines, we found deletion of one or more exons in the c-fes mRNA. Those deletions predicted truncation of conserved domains (CDC15/FCH and SH2) involved in protein-protein interactions. No deletions were found, however, within the kinase domain. We infer that alternative splicing generates a family of c-Fes proteins. This may be a mechanism to direct the c-Fes kinase domain to different subcellular locations and/or substrates at specific stages of myeloid cell differentiation.

Molecular mechanisms regulating the tumor-targeting potential of splice-activated gene expression

Previous studies have suggested that differences in the ability of normal and malignant cells to process certain alternatively spliced pre-mRNA transcripts can be exploited as a potentially powerful means of targeting the expression of therapeutic genes to tumor cells in vivo and in vitro. Specifically, it was shown that efficient processing of minigene constructs containing the alternatively spliced CD44 exons v9 and v10 only occurs in tumor cells that express CD44 isoforms that incorporate these exons (e.g. CD44R1). In the present study, efforts were made to define the molecular mechanisms that underlie the apparent specificity of this process. RT-PCR analysis and DNA sequencing were used to characterize the various splicing events that occur between CD44 exons v8, v9 and v10 following transfection of minigene constructs containing these various exons into CD44R1-positive (PC3) and CD44R1-negative (T24) cell lines. The results obtained confirm that although the v8-v9 intron is efficiently removed in both CD44R1-positive and CD44R1-negative cells, the corresponding v9-v10 intron is accurately spliced and the exons appropriately joined only in lines that express v10-containing CD44 isoforms (e.g. PC3). In CD44R1-negative cell lines (e.g. T24) alternative 5' and 3' splice sites located within the v9-v10 intron are preferentially used, resulting in various portions of the intron being retained within the final processed mRNA product. It is proposed that identification of these functionally important intronic sequence elements will facilitate the development of second generation "splice activated gene expression" vectors that may prove useful in various cancer gene therapy applications.

Insertion/deletion and nucleotide polymorphism data reveal constraints in Drosophila melanogaster introns and intergenic regions

Our study of nucleotide sequence and insertion/deletion polymorphism in Drosophila melanogaster noncoding DNA provides evidence for selective pressures in both intergenic regions and introns (of the large size class). Intronic and intergenic sequences show a similar polymorphic deletion bias. Insertions have smaller sizes and higher frequencies than deletions, supporting the hypothesis that insertions are selected to compensate for the loss of DNA caused by deletion bias. Analysis of a simple model of selective constraints suggests that the blocks of functional elements located in intergenic sequences are on average larger than those in introns, while the length distribution of relatively unconstrained sequences interspaced between these blocks is similar in intronic and intergenic regions.

Intron size correlates positively with recombination rate in Caenorhabditis elegans

A negative correlation between intron size and recombination rate has been reported for the Drosophila melanogaster and human genomes. Population-genetic models suggest that this pattern could be caused by an interaction between recombination rate and the efficacy of natural selection. To test this idea, we examined variation in intron size and recombination rate across the genome of the nematode Caenorhabditis elegans. Interestingly, we found that intron size correlated positively with recombination rate in this species.

Hypervariable and Highly Divergent Intron?Exon Organizations in the Chordate Oikopleura dioica

Oikopleura dioica is a pelagic tunicate with a very small genome and a very short life cycle. In order to investigate the intron-exon organizations in Oikopleura, we have isolated and characterized ribosomal protein EF-1alpha, Hox, and alpha-tubulin genes. Their intron positions have been compared with those of the same genes from various invertebrates and vertebrates, including four species with entirely sequenced genomes. Oikopleura genes, like Caenorhabditis genes, have introns at a large number of nonconserved positions, which must originate from late insertions or intron sliding of ancient insertions. Both species exhibit hypervariable intron-exon organization within their alpha-tubulin gene family. This is due to localization of most nonconserved intron positions in single members of this gene family. The hypervariability and divergence of intron positions in Oikopleura and Caenorhabditis may be related to the predominance of short introns, the processing of which is not very dependent upon the exonic environment compared to large introns. Also, both species have an undermethylated genome, and the control of methylation-induced point mutations imposes a control on exon size, at least in vertebrate genes. That introns placed at such variable positions in Oikopleura or C. elegans may serve a specific purpose is not easy to infer from our current knowledge and hypotheses on intron functions. We propose that new introns are retained in species with very short life cycles, because illegitimate exchanges including gene conversion are repressed. We also speculate that introns placed at gene-specific positions may contribute to suppressing these exchanges and thereby favor their own persistence.

Database and analyses of known alternatively spliced genes in plants

Alternative splicing is an important cellular mechanism that increases the diversity of gene products. The number of alternatively spliced genes reported so far in plants is much smaller than that in mammals, but is increasing as a result of the explosive growth of available EST and genomic sequences. We have searched for all alternatively spliced genes reported in GenBank and PubMed in all plant species under Viridiplantae. After careful merging and manual review of the search results, we obtained a comprehensive, high-quality collection of 168 genes reported to be alternatively spliced in plants, spanning 44 plant species (March 22, 2003 update). We developed a relational database with Web-based user interface to store and present the data, named the Plant Alternative Splicing Database (PASDB), freely available at http://pasdb.genomics.org.cn. We analyzed the functional categories that these genes belong to using the Gene Ontology. We also analyzed in detail the biological roles and gene structures of the four genes that are known to be alternatively spliced in more than one plant species. Finally, we studied the structural features of the splice sites in the alternatively spliced genes.

Detection and evaluation of intron retention events in the human transcriptome

Alternative splicing is a very frequent phenomenon in the human transcriptome. There are four major types of alternative splicing: exon skipping, alternative 3' splice site, alternative 5' splice site, and intron retention. Here we present a large-scale analysis of intron retention in a set of 21,106 known human genes. We observed that 14.8% of these genes showed evidence of at least one intron retention event. Most of the events are located within the untranslated regions (UTRs) of human transcripts. For those retained introns interrupting the coding region, the GC content, codon usage, and the frequency of stop codons suggest that these sequences are under selection for coding potential. Furthermore, 26% of the introns within the coding region participate in the coding of a protein domain. A comparison with mouse shows that at least 22% of all informative examples of retained introns in human are also present in the mouse transcriptome. We discuss that the data we present suggest that a significant fraction of the observed events is not spurious and might reflect biological significance. The analyses also allowed us to generate a reliable set of intron retention events that can be used for the identification of splicing regulatory elements.

Inferred functions of "novel" genes identified in fibroblasts chondroinduced by demineralized bone

Little is known about the cellular mechanisms that control postnatal chondrocyte differentiation. As a first step towards identifying those mechanisms, gene expression shifts were characterized in an in vitro model of chondroinduction. In previous studies, several functional classes of genes (cytoskeletal and matrix elements, cell adhesion proteins, peptide growth factors, and signal transduction proteins) were found to be altered in human dermal fibroblasts (hDFs) cultured in porous collagen sponges with chondroinductive demineralized bone powder (DBP) for 3 days. In addition, a number of "novel" sequences were identified. In this study, molecular techniques were combined with computational methods to characterize those sequences. Gene expression of all 10 novel sequences tested was found in hDFs by RT-PCR. The sequences were compared to the human genome, and their cellular functions were inferred from genes that mapped to the same chromosomal coordinates. Only one of the novel sequences contained a protein-coding region (kinesin superfamily protein 26B). The others contained 3' untranslated (osteonectin, alpha-V integrin, RAP2B) or other untranslated regions (PTPN21, GAS6) of mRNAs. The cellular functions of the DBP-regulated genes described in this study fall into similar categories as those previously identified. These results provide new details on the cellular response of hDFs exposed to DBP.

Alternative splicing and evolution

Alternative splicing is a critical post-transcriptional event leading to an increase in the transcriptome diversity. Recent bioinformatics studies revealed a high frequency of alternative splicing. Although the extent of AS conservation among mammals is still being discussed, it has been argued that major forms of alternatively spliced transcripts are much better conserved than minor forms. It suggests that alternative splicing plays a major role in genome evolution allowing new exons to evolve with less constraint.

ORESTES are enriched in rare exon usage variants affecting the encoded proteins

A significant fraction of the variability found in the human transcriptome is due to alternative splicing, including alternative exon usage (AEU), intron retention and use of cryptic splice sites. We present a comparison of a large-scale analysis of AEU in the human transcriptome through genome mapping of Open Reading Frame ESTs (ORESTES) and conventional ESTs. It is shown here that ORESTES probe low abundant messages more efficiently. In addition, most of the variants detected by ORESTES affect the structure of the corresponding proteins.

Isolation and characterization of human and mouse WDR19,a novel WD-repeat protein exhibiting androgen-regulated expression in prostate epithelium

Androgens regulate important processes involved in the normal development and function of the human and rodent prostate glands. Here we report the isolation and characterization of a new androgen-regulated gene, designated WDR19, that encodes repeating sequence motifs found in the WD-repeat family of proteins. The WD repeat is a conserved domain of approximately 40 amino acids that is typically bracketed by glycine-histidine and tryptophan-aspartic acid (WD) dipeptides. WD-repeat proteins are a large group of structurally related proteins that participate in a wide range of cellular functions, including transmembrane signaling, mRNA modification, vesicle formation, and vesicular trafficking. The WDR19 gene comprises 36 exons and is located on chromosome 4p15-4p11. The predicted protein contains six WD repeats, a clathrin heavy-chain repeat, and three transmembrane domains. Sequence analysis reveals that the WDR19 gene is conserved from Caenorhabditis elegans to human. WDR19 is expressed in normal and neoplastic prostate epithelium as demonstrated by RNA in situ hybridization and is regulated by androgenic hormones. WDR19 transcripts exhibit alternative splicing in which two isoforms appear to be prostate restricted, a property that could be exploited for designing diagnostic or therapeutic strategies for prostate carcinoma.

Expression of U1 small nuclear ribonucleoprotein 70K antisense transcript using APETALA3 promoter suppresses the development of sepals and petals

U1 small nuclear ribonucleoprotein (snRNP)-70K (U1-70K), a U1 snRNP-specific protein, is involved in the early stages of spliceosome formation. In non-plant systems, it is involved in constitutive and alternative splicing. It has been shown that U1snRNP is dispensable for in vitro splicing of some animal pre-mRNAs, and inactivation of U1-70K in yeast (Saccharomyces cerevisiae) is not lethal. As in yeast and humans (Homo sapiens), plant U1-70K is coded by a single gene. In this study, we blocked the expression of Arabidopsis U1-70K in petals and stamens by expressing U1-70K antisense transcript using the AP3 (APETALA3) promoter specific to these floral organs. Flowers of transgenic Arabidopsis plants expressing U1-70K antisense transcript showed partially developed stamens and petals that are arrested at different stages of development. In some transgenic lines, flowers have rudimentary petals and stamens and are male sterile. The severity of the phenotype is correlated with the level of the antisense transcript. Molecular analysis of transgenic plants has confirmed that the observed phenotype is not due to disruption of whorl-specific homeotic genes, AP3 or PISTILLATA, responsible for petal and stamen development. The AP3 transcript was not detected in transgenic flowers with severe phenotype. Flowers of Arabidopsis plants transformed with a reporter gene driven by the same promoter showed no abnormalities. These results show that U1-70K is necessary for the development of sepals and petals and is an essential gene in plants.

Novel alternatively spliced variant with a deletion of 52 BP in exon 6 of the progesterone receptor gene is observed frequently in breast cancer tissues

The human progesterone receptor (PR) is a ligand-activated nuclear transcription factor that mediates progesterone action in target tissues. We found a novel alternatively spliced variant (ASV) of the PR mRNA in breast cancer tissues. The deleted transcript was characterized by an out-of-frame deletion of 52 bp in exon 6 (PR delta6/2 ASV). The PR delta6/2 ASV mRNA results in a partial defect in the region of the ligand-binding domain of the hormone receptor, where conserved residues are missing from the core of the protein. To clarify the clinical significance of the PR delta6/2 ASV, we investigated the expression of this ASV in noncancerous and cancerous tissues from patients with breast cancer using RT-PCR. The novel PR delta6/2 mRNA was detected in 24 of 39 (61.5%) cancerous tissues and in 3 of 39 (7.7%) noncancerous tissues from patients with breast cancer. PR delta6/2 ASV mRNA was expressed more frequently in breast cancer tissues than in noncancerous tissues (p < 0.0001), which suggests a possible relationship between the expression of PR delta6/2 and breast cancer. Our observations may provide a novel strategy for the genetic diagnosis of breast cancer.

Expression profile of a gamma-deletion variant of the human telomerase reverse transcriptase gene

The human telomerase reverse transcriptase (hTERT) is an essential component of the holoenzyme complex that adds telomeric repeats to the ends of chromosomes. The hTERT transcript has been shown to have two deletion type alternative splicing sites. One deletion site induces the alpha-deletion variant, lacking 36 bp from exon 6, and the other induces the beta-deletion variant, lacking 182 bp from exons 7 and 8. Here, we identified a novel deletion variant of the hTERT transcript in hepatocellular carcinoma cell lines. The deleted transcript was characterized by an in-frame deletion of 189 bp, spanning nucleotides 2710 to 2898, corresponding to the complete loss of exon 11 (gamma-deletion). The region lacking in the gamma-deletion lies within RT motifs D and E, suggesting that it is missing conserved residues from the catalytic core of the protein. Both gamma- and alpha-deletion variants were occasionally detected, but the beta-deletion variant was frequently observed. Our results may provide important information for more detailed studies on the regulation of telomerase activity.

Characterization of the splice variant Gbeta3v of the human G-protein Gbeta3 subunit

A polymorphism (C825T) in the gene of the G-protein Gbeta3 (GNB3) has been the subject of numerous studies which have shown that the 825T-allele is associated with several cardiovascular and metabolic disorders. Furthermore, the T-allele is associated with the occurrence of the splice variant Gbeta3s which has been implicated in the pathogenesis of these disorders. Here, we characterise a novel splice variant of GNB3, termed Gbeta3v, which is generated by alternative splicing of parts from intron 9 as a novel exon 10a. Gbeta proteins belong to the superfamily of propeller proteins composed of seven regular WD-domains. In Gbeta3v, four of these WD-domains are retained but the protein has a novel C terminus. Gbeta3v forms dimers with Ggamma3 and Ggamma12 but these Gbetagamma complexes do not stimulate phospholipase C-beta2. Thus, a physiological role of Gbeta3v remains to be established. Gbeta3v transcripts are detectable in a wide variety of cells and tissues including fibroblasts, B lymphoblasts, retinoblastoma cells, retina, brain, umbilical cord and colon. However, there is no association with an allele of the GNB3 C825T polymorphism, which suggests that Gbeta3v does not contribute to the complex phenotype observed in association with the 825T-allele.

Increase of functional diversity by alternative splicing

A large-scale analysis of protein isoforms arising from alternative splicing shows that alternative splicing tends to insert or delete complete protein domains more frequently than expected by chance, whereas disruption of domains and other structural modules is less frequent. If domain regions are disrupted, the functional effect, as predicted from 3D structure, is frequently equivalent to removal of the entire domain. Also, short alternative splicing events within domains, which might preserve folded structure, target functional residues more frequently than expected. Thus, it seems that positive selection has had a major role in the evolution of alternative splicing.

Alternative core promoters regulate tissue-specific transcription from the autoimmune diabetes-related ICA1 (ICA69) gene locus

Islet cell autoantigen 69-kDa (ICA69), protein product of the human ICA1 gene, is one target of the immune processes defining the pathogenesis of Type 1 diabetes. We have characterized the genomic structure and functional promoters within the 5'-regulatory region of ICA1. 5'-RNA ligase-mediated rapid amplification of cDNA ends evaluation of ICA1 transcripts expressed in human islets, testis, heart, and cultured neuroblastoma cells reveals that three 5'-untranslated region exons are variably expressed from the ICA1 gene in a tissue-specific manner. Surrounding the transcription initiation sites are motifs characteristic of non-TATA, non-CAAT, GC-rich promoters, including consensus Sp1/GC boxes, an initiator element, cAMP-responsive element-binding protein (CREB) sites, and clusters of other putative transcription factor sites within a genomic CpG island. Luciferase reporter constructs demonstrate that the first two ICA1 exon promoters reciprocally stimulate luciferase expression within islet- (RIN 1046-38 cells) and brain-derived (NMB7) cells in culture; the exon A promoter exhibits greater activity in islet cells, whereas the exon B promoter more efficiently activates transcription in neuronal cells. Mutation of a CREB site within the ICA1 exon B promoter significantly enhances transcriptional activity in both cell lines. Our basic understanding of expression from the functional core promoter elements of ICA1 is an important advance that will not only add to our knowledge of the ICA69 autoantigen but will also facilitate a rational approach to discover the function of ICA69 and to identify relevant ICA1 promoter polymorphisms and their potential associations with disease.

SAF-2, a splice variant of SAF-1, acts as a negative regulator of transcription

Serum amyloid A-activating factor-1 (SAF-1), a Cys(2)His(2)-type zinc finger transcription factor, regulates inflammation-induced expression of serum amyloid A protein that is linked to the pathogenesis of reactive amyloidosis, rheumatoid arthritis, and atherosclerosis. Here we report the identification of a novel splice variant, SAF-2, of the SAF family bearing strong sequence similarity to SAF-1. The N-terminal 426 amino acids of both SAF-1 and SAF-2 are identical containing two polyalanine tracts, one proline-rich domain, and six zinc fingers. However, the C terminus of SAF-2 containing two additional zinc fingers is different from SAF-1, which indicates the capability of different biochemical function. We show that SAF-2 interacts more avidly with the SAF-binding element, but its transactivation potential is much lower than SAF-1. Furthermore, co-expression of SAF-2 markedly suppresses SAF-1-regulated promoter function. Finally, we show that the level of SAF-2 protein is reduced during many inflammatory conditions, whereas the SAF-1 protein level remains unchanged. Together, these data suggest that the relative abundance of SAF-2 plays a critical role in the fine tuned regulation of inflammation-responsive genes that are controlled by SAF-1.

Micro-processing events in mRNAs identified by DHPLC analysis

Post-transcriptional processes such as alternative splicing and RNA editing have a huge impact on the diversity of the proteome. Detecting alternatively spliced transcripts is difficult when they are rare. In addition, edited transcripts often differ from the genomic sequence by only a few nucleotides. Denaturing high performance liquid chromatography (DHPLC) is routinely used for single nucleotide polymorphism detection and we used this method to detect alternatively spliced or edited transcripts. As the sites of RNA editing appear to be conserved within gene families, we investigated whether editing sites are conserved in the murine homologue of the Drosophila cacophony transcript encoding the alpha1 subunit of a voltage-gated calcium channel that is edited at 10 independent positions. Although DHPLC analysis detects RNA editing at as low as 3% in control transcripts, no evidence of RNA editing was found in the analysed murine transcript. However an alternative exon was identified at the 3' end of the mouse Cacna1alpha transcript and an alternative micro-exon encoding only two amino acids (NP) was found in the extracellular loop before the IVS4 helix in the same transcript. In the homologous Drosophila transcript a micro-exon also encoding two amino acids was found at the same position before the IVS4 helix.

Hormonal genomics

The availability of the human genomic sequence is changing the way in which biological questions are addressed. Based on the prediction of genes from nucleotide sequences, homologies among their encoded amino acids can be analyzed and used to place them in distinct families. This serves as a first step in building hypotheses for testing the structural and functional properties of previously uncharacterized paralogous genes. As genomic information from more organisms becomes available, these hypotheses can be refined through comparative genomics and phylogenetic studies. Instead of the traditional single-gene approach in endocrine research, we are beginning to gain an understanding of entire mammalian genomes, thus providing the basis to reveal subfamilies and pathways for genes involved in ligand signaling. The present review provides selective examples of postgenomic approaches in the analysis of novel genes involved in hormonal signaling and their chromosomal locations, polymorphisms, splicing variants, differential expression, and physiological function. In the postgenomic era, scientists will be able to move from a gene-by-gene approach to a reconstructionistic one by reading the encyclopedia of life from a global perspective. Eventually, a community-based approach will yield new insights into the complexity of intercellular communications, thereby offering us an understanding of hormonal physiology and pathophysiology.

Alternative splicing: combinatorial output from the genome

Alternative splicing has emerged as a mechanism that can account for a large proportion of the disparity between the modest number of genes in the human genome and the much higher complexity of the expressed proteome. At least a third, and probably the majority, of human genes are alternatively spliced, and some genes can generate thousands of protein isoforms by complex alternative splicing events. Analysis of the transcriptome will therefore require the development of massively parallel technologies that are able to encompass the complexity arising from alternative splicing.

A draft sequence of the rice genome (Oryza sativa L. ssp. indica)

We have produced a draft sequence of the rice genome for the most widely cultivated subspecies in China, Oryza sativa L. ssp. indica, by whole-genome shotgun sequencing. The genome was 466 megabases in size, with an estimated 46,022 to 55,615 genes. Functional coverage in the assembled sequences was 92.0%. About 42.2% of the genome was in exact 20-nucleotide oligomer repeats, and most of the transposons were in the intergenic regions between genes. Although 80.6% of predicted Arabidopsis thaliana genes had a homolog in rice, only 49.4% of predicted rice genes had a homolog in A. thaliana. The large proportion of rice genes with no recognizable homologs is due to a gradient in the GC content of rice coding sequences.

Identification and ethnic distribution of major haplotypes in the gene GNB3 encoding the G-protein beta3 subunit

The 825T-allele of the gene GNB3 encoding the G protein beta3 subunit is associated with hypertension and obesity, and identifies individuals highly responsive to diuretic therapy. Gbeta3s, a Gbeta3 protein variant generated by alternative splicing in carriers of the 825T-allele, is linked to increased signal transduction and is a potential cause for the observed pathophysiology. Here, we searched the entire GNB3 gene for additional polymorphisms and analysed their prevalence in Caucasian, black African and Asian populations. We detected six novel single nucleotide polymorphisms which were termed according to their location as G76A, G1906T, G2906A, A3882C, G5177A, and G5249A. Furthermore, we found a CACA-insertion-deletion polymorphism at position 6496. Genotyping and association studies resulted in the definition of two major GNB3 haplotypes, termed 'C-haplotype' (alleles 825C, 3882A, 5249G, 6496CACA-) and 'T-haplotype' (alleles 825T, 3882C, 5249A, 6496CACA+). Molecular modelling studies revealed that the pre-mRNA structures of both haplotypes exhibit marked differences which may account for the alternative splicing predominantly observed with the T-haplotype. The prevalence of these haplotypes in major ethnic populations differs considerably. Furthermore, we detected additional frequent GNB3 polymorphisms. These variants were restricted to one or two major ethnic populations. Our results will aid future studies on population-specific effects of the GNB3 variants on risk and course of frequent diseases, including hypertension, obesity, stroke and myocardial infarction. Furthermore, they will contribute to the understanding of GNB3-related population-specific pharmacogenetic differences in the response to major drugs, as already shown for diuretics and antidepressants.

Intergenic mRNA molecules resulting from trans-splicing

Accumulated recent evidence is indicating that alternative splicing represents a generalized process that increases the complexity of human gene expression. Here we show that mRNA production may not necessarily be limited to single genes, as human liver also has the potential to produce a variety of hybrid cytochrome P450 3A mRNA molecules. The four known cytochrome P450 3A genes in humans, CYP3A4, CYP3A5, CYP3A7, and CYP3A43, share a high degree of similarity, consist of 13 exons with conserved exon-intron boundaries, and form a cluster on chromosome 7. The chimeric CYP3A mRNA molecules described herein are characterized by CYP3A43 exon 1 joined at canonical splice sites to distinct sets of CYP3A4 or CYP3A5 exons. Because the CYP3A43 gene is in a head-to-head orientation with the CYP3A4 and CYP3A5 genes, bypassing transcriptional termination can not account for the formation of hybrid CYP3A mRNAs. Thus, the mechanism generating these molecules has to be an RNA processing event that joins exons of independent pre-mRNA molecules, i.e. trans-splicing. Using quantitative real-time polymerase chain reaction, the ratio of one CYP3A43/3A4 intergenic combination was estimated to be approximately 0.15% that of the CYP3A43 mRNAs. Moreover, trans-splicing has been found not to interfere with polyadenylation. Heterologous expression of the chimeric species composed of CYP3A43 exon 1 joined to exons 2-13 of CYP3A4 revealed catalytic activity toward testosterone.