Genome-wide detection of tissue-specific alternative splicing in the human transcriptome

A family of splice variants of CstF-64 expressed in vertebrate nervous systems

Background

Alternative splicing and polyadenylation are important mechanisms for creating the proteomic diversity necessary for the nervous system to fulfill its specialized functions. The contribution of alternative splicing to proteomic diversity in the nervous system has been well documented, whereas the role of alternative polyadenylation in this process is less well understood. Since the CstF-64 polyadenylation protein is known to be an important regulator of tissue-specific polyadenylation, we examined its expression in brain and other organs.

Results

We discovered several closely related splice variants of CstF-64 – collectively called βCstF-64 – that could potentially contribute to proteomic diversity in the nervous system. The βCstF-64 splice variants are found predominantly in the brains of several vertebrate species including mice and humans. The major βCstF-64 variant mRNA is generated by inclusion of two alternate exons (that we call exons 8.1 and 8.2) found between exons 8 and 9 of the CstF-64 gene, and contains an additional 147 nucleotides, encoding 49 additional amino acids. Some variants of βCstF-64 contain only the first alternate exon (exon 8.1) while other variants contain both alternate exons (8.1 and 8.2). In mice, the predominant form of βCstF-64 also contains a deletion of 78 nucleotides from exon 9, although that variant is not seen in any other species examined, including rats. Immunoblot and 2D-PAGE analyses of mouse nuclear extracts indicate that a protein corresponding to βCstF-64 is expressed in brain at approximately equal levels to CstF-64. Since βCstF-64 splice variant family members were found in the brains of all vertebrate species examined (including turtles and fish), this suggests that βCstF-64 has an evolutionarily conserved function in these animals. βCstF-64 was present in both pre- and post-natal mice and in different regions of the nervous system, suggesting an important role for βCstF-64 in neural gene expression throughout development. Finally, experiments in representative cell lines suggest that βCstF-64 is expressed in neurons but not glia.

Conclusion

This is the first report of a family of splice variants encoding a key polyadenylation protein that is expressed in a nervous system-specific manner. We propose that βCstF-64 contributes to proteomic diversity by regulating alternative polyadenylation of neural mRNAs.

Genomic landscape of developing male germ cells

Spermatogenesis is a highly orchestrated developmental process by which spermatogonia develop into mature spermatozoa. This process involves many testis- or male germ cell-specific gene products whose expressions are strictly regulated. In the past decade the advent of high-throughput gene expression analytical techniques has made functional genomic studies of this process, particularly in model animals such as mice and rats, feasible and practical. These studies have just begun to reveal the complexity of the genomic landscape of the developing male germ cells. Over 50% of the mouse and rat genome are expressed during testicular development. Among transcripts present in germ cells, 40% - 60% are uncharacterized. A number of genes, and consequently their associated biological pathways, are differentially expressed at different stages of spermatogenesis. Developing male germ cells present a rich repertoire of genetic processes. Tissue-specific as well as spermatogenesis stage-specific alternative splicing of genes exemplifies the complexity of genome expression. In addition to this layer of control, discoveries of abundant presence of antisense transcripts, expressed psuedogenes, non-coding RNAs (ncRNA) including long ncRNAs, microRNAs (miRNAs) and Piwi-interacting RNAs (piRNAs), and retrogenes all point to the presence of multiple layers of expression and functional regulation in male germ cells. It is anticipated that application of systems biology approaches will further our understanding of the regulatory mechanism of spermatogenesis.

The PRMT1 gene expression pattern in colon cancer

The methylation of arginine has been implicated in many cellular processes, such as regulation of transcription, mRNA splicing, RNA metabolism and transport. The enzymes responsible for this modification are the protein arginine methyltransferases. The most abundant methyltransferase in human cells is protein arginine methyltransferase 1. Methylation processes appear to interfere in the emergence of several diseases, including cancer. During our study, we examined the expression pattern of protein arginine methyltransferase 1 gene in colon cancer patients. The emerging results showed that the expression of one of the gene variants is associated with statistical significant probability to clinical and histological parameters, such as nodal status and stage. This is a first attempt to acquire an insight on the possible relation of the expression pattern of protein arginine methyltransferase 1 and colon cancer progression.

Alternative isoform regulation in human tissue transcriptomes

Through alternative processing of pre-mRNAs, individual mammalian genes often produce multiple mRNA and protein isoforms that may have related, distinct or even opposing functions. Here we report an in-depth analysis of 15 diverse human tissue and cell line transcriptomes based on deep sequencing of cDNA fragments, yielding a digital inventory of gene and mRNA isoform expression. Analysis of mappings of sequence reads to exon-exon junctions indicated that 92-94% of human genes undergo alternative splicing (AS), ∼86% with a minor isoform frequency of 15% or more. Differences in isoform-specific read densities indicated that a majority of AS and of alternative cleavage and polyadenylation (APA) events vary between tissues, while variation between individuals was ∼2- to 3-fold less common. Extreme or ‘switch-like’ regulation of splicing between tissues was associated with increased sequence conservation in regulatory regions and with generation of full-length open reading frames. Patterns of AS and APA were strongly correlated across tissues, suggesting coordinated regulation of these processes, and sequence conservation of a subset of known regulatory motifs in both alternative introns and 3′ UTRs suggested common involvement of specific factors in tissue-level regulation of both splicing and polyadenylation.

Human and mouse switch-like genes share common transcriptional regulatory mechanisms for bimodality

Background

Gene expression is controlled over a wide range at the transcript level through complex interplay between DNA and regulatory proteins, resulting in profiles of gene expression that can be represented as normal, graded, and bimodal (switch-like) distributions. We have previously performed genome-scale identification and annotation of genes with switch-like expression at the transcript level in mouse, using large microarray datasets for healthy tissue, in order to study the cellular pathways and regulatory mechanisms involving this class of genes. We showed that a large population of bimodal mouse genes encoding for cell membrane and extracellular matrix proteins is involved in communication pathways. This study expands on previous results by annotating human bimodal genes, investigating their correspondence to bimodality in mouse orthologs and exploring possible regulatory mechanisms that contribute to bimodality in gene expression in human and mouse.

Results

Fourteen percent of the human genes on the HGU133A array (1847 out of 13076) were identified as bimodal or switch-like. More than 40% were found to have bimodal mouse orthologs. KEGG pathways enriched for bimodal genes included ECM-receptor interaction, focal adhesion, and tight junction, showing strong similarity to the results obtained in mouse. Tissue-specific modes of expression of bimodal genes among brain, heart, and skeletal muscle were common between human and mouse. Promoter analysis revealed a higher than average number of transcription start sites per gene within the set of bimodal genes. Moreover, the bimodal gene set had differentially methylated histones compared to the set of the remaining genes in the genome.

Conclusion

The fact that bimodal genes were enriched within the cell membrane and extracellular environment make these genes as candidates for biomarkers for tissue specificity. The commonality of the important roles bimodal genes play in tissue differentiation in both the human and mouse indicates the potential value of mouse data in providing context for human tissue studies. The regulation motifs enriched in the bimodal gene set (TATA boxes, alternative promoters, methlyation) have known associations with complex diseases, such as cancer, providing further potential for the use of bimodal genes in studying the molecular basis of disease.

The effect of intron length on exon creation ratios during the evolution of mammalian genomes

Recent studies report that alternatively spliced exons tend to occur in longer introns, which is attributed to the length constraints for splice site pairing for the two major splicing mechanisms, intron definition versus exon definition. Using genome-wide studies of EST and microarray data from human and mouse, we have analyzed the distribution of various subsets of alternatively spliced exons, based on their inclusion level and evolutionary history, versus increasing intron length. Alternative exons may be included in either a major or minor fraction of all transcripts (known as major-form and minor-form exons, respectively). We find that major-form exons are seven- to eightfold more likely to be contained in short introns (<400 nt) than minor-form exons, which occur preferentially in longer introns. Since minor-form exons are more likely to be novel (approximately 75%), this implied that novel exons arise more frequently in longer introns. To test this hypothesis, we used whole genome alignments to classify exons according to their phylogenetic age. We find that older exons, i.e., exons that are conserved in all mammals, predominate at shorter intron lengths, for both major- and minor-form exons. In contrast, exons that arose recently during primate evolution are more prevalent at longer intron lengths (>1000 nt). This suggests that the observed correlation of longer intron lengths with alternatively spliced exons may be at least partly due to biases in the probability of exon creation, which is higher in long introns.

MADS: a new and improved method for analysis of differential alternative splicing by exon-tiling microarrays

We describe a method, microarray analysis of differential splicing (MADS), for discovery of differential alternative splicing from exon-tiling microarray data. MADS incorporates a series of low-level analysis algorithms motivated by the "probe-rich" design of exon arrays, including background correction, iterative probe selection, and removal of sequence-specific cross-hybridization to off-target transcripts. We used MADS to analyze Affymetrix Exon 1.0 array data on a mouse neuroblastoma cell line after shRNA-mediated knockdown of the splicing factor polypyrimidine tract binding protein (PTB). From a list of exons with predetermined inclusion/exclusion profiles in response to PTB depletion, MADS recognized all exons known to have large changes in transcript inclusion levels and offered improvement over Affymetrix's analysis procedure. We also identified numerous novel PTB-dependent splicing events. Thirty novel events were tested by RT-PCR and 27 were confirmed. This work demonstrates that the exon-tiling microarray design is an efficient and powerful approach for global, unbiased analysis of pre-mRNA splicing.

Solving the problem of Trans-Genomic Query with alignment tables

The trans-genomic query (TGQ) problem--enabling the free query of biological information, even across genomes--is a central challenge facing bioinformatics. Solutions to this problem can alter the nature of the field, moving it beyond the jungle of data integration and expanding the number and scope of questions that can be answered. An alignment table is a binary relationship on locations (sequence segments). An important special case of alignment tables are hit tables ? tables of pairs of highly similar segments produced by alignment tools like BLAST. However, alignment tables also include general binary relationships, and can represent any useful connection between sequence locations. They can be curated, and provide a high-quality queryable backbone of connections between biological information. Alignment tables thus can be a natural foundation for TGQ, as they permit a central part of the TGQ problem to be reduced to purely technical problems involving tables of locations.Key challenges in implementing alignment tables include efficient representation and indexing of sequence locations. We define a location datatype that can be incorporated naturally into common off-the-shelf database systems. We also describe an implementation of alignment tables in BLASTGRES, an extension of the open-source POSTGRESQL database system that provides indexing and operators on locations required for querying alignment tables. This paper also reviews several successful large-scale applications of alignment tables for Trans-Genomic Query. Tables with millions of alignments have been used in queries about alternative splicing, an area of genomic analysis concerning the way in which a single gene can yield multiple transcripts. Comparative genomics is a large potential application area for TGQ and alignment tables.

Resurrecting the role of transcription factor change in developmental evolution

A long-standing question in evolutionary and developmental biology concerns the relative contribution of cis-regulatory and protein changes to developmental evolution. Central to this argument is which mutations generate evolutionarily relevant phenotypic variation? A review of the growing body of evolutionary and developmental literature supports the notion that many developmentally relevant differences occur in the cis-regulatory regions of protein-coding genes, generally to the exclusion of changes in the protein-coding region of genes. However, accumulating experimental evidence demonstrates that many of the arguments against a role for proteins in the evolution of gene regulation, and the developmental evolution in general, are no longer supported and there is an increasing number of cases in which transcription factor protein changes have been demonstrated in evolution. Here, we review the evidence that cis-regulatory evolution is an important driver of phenotypic evolution and provide examples of protein-mediated developmental evolution. Finally, we present an argument that the evolution of proteins may play a more substantial, but thus far underestimated, role in developmental evolution.

Domains of WNK1 kinase in the regulation of ROMK1

WNK1 kinase belongs to a family of serine-threonine protein kinases with an atypical placement of the catalytic lysine. Increased expression of WNK1 causes hypertension and hyperkalemia in humans. WNK1 inhibits renal potassium channel ROMK1 by enhancing its endocytosis, likely contributing to hyperkalemia in affected patients. The domains of WNK1 involved in inhibition of ROMK1 have not been completely elucidated. Here, we reported that an NH2-terminal proline-rich domain (N-PRD; amino acids 1-119) is necessary and sufficient for WNK1 inhibition of ROMK1. A region (named "NL" for N-linker; amino acids 120-220) located between N-PRD and the kinase domain of WNK1 (amino acids 220-491) antagonized the inhibition of ROMK1 caused by N-PRD. The WNK1 kinase domain reversed the antagonism of NL on N-PRD. Mutagenesis studies revealed that charge-charge interactions between two conserved catalytic residues (Lys-233 and Asp-368) within the kinase domain (not the kinase activity) are critical for kinase domain to reverse the antagonism of NL domain. The WNK1 autoinhibitory domain (AID; amino acids 491-555) also affected ROMK, presumably by modulating the kinase domain conformation. Mutations of two conserved phenylalanine abolished the ability of AID to modulate ROMK1. Finally, the first coiled-coil domain (CC1; amino acids 555-640) of WNK1 alleviated the effect of AID domain toward kinase domain. Thus, multiple intra- and/or intermolecular interactions of WNK1 domains are at play for regulation of ROMK1 by WNK1.

Genome-wide survey of allele-specific splicing in humans

Background

Accurate mRNA splicing depends on multiple regulatory signals encoded in the transcribed RNA sequence. Many examples of mutations within human splice regulatory regions that alter splicing qualitatively or quantitatively have been reported and allelic differences in mRNA splicing are likely to be a common and important source of phenotypic diversity at the molecular level, in addition to their contribution to genetic disease susceptibility. However, because the effect of a mutation on the efficiency of mRNA splicing is often difficult to predict, many mutations that cause disease through an effect on splicing are likely to remain undiscovered.

Results

We have combined a genome-wide scan for sequence polymorphisms likely to affect mRNA splicing with analysis of publicly available Expressed Sequence Tag (EST) and exon array data. The genome-wide scan uses published tools and identified 30,977 SNPs located within donor and acceptor splice sites, branch points and exonic splicing enhancer elements. For 1,185 candidate splicing polymorphisms the difference in splicing between alternative alleles was corroborated by publicly available exon array data from 166 lymphoblastoid cell lines. We developed a novel probabilistic method to infer allele-specific splicing from EST data. The method uses SNPs and alternative mRNA isoforms mapped to EST sequences and models both regulated alternative splicing as well as allele-specific splicing. We have also estimated heritability of splicing and report that a greater proportion of genes show evidence of splicing heritability than show heritability of overall gene expression level. Our results provide an extensive resource that can be used to assess the possible effect on splicing of human polymorphisms in putative splice-regulatory sites.

Conclusion

We report a set of genes showing evidence of allele-specific splicing from an integrated analysis of genomic polymorphisms, EST data and exon array data, including several examples for which there is experimental evidence of polymorphisms affecting splicing in the literature. We also present a set of novel allele-specific splicing candidates and discuss the strengths and weaknesses of alternative technologies for inferring the effect of sequence variants on mRNA splicing.

Fine expression profiling of full-length transcripts using a size-unbiased cDNA library prepared with the vector-capping method

Recently, we have developed a vector-capping method for constructing a full-length cDNA library. In the present study, we performed in-depth analysis of the vector-capped cDNA library prepared from a single type of cell. As a result of single-pass sequencing analysis of 24 000 clones randomly isolated from the unamplified library, we identified 19 951 full-length cDNA clones whose intactness was confirmed by the presence of an additional G at their 5' end. The full-length cDNA content was >95%. Mapping these sequences to the human genome, we identified 4513 transcriptional units that include 36 antisense transcripts against known genes. Comparison of the frequencies of abundant clones showed that the expression profiles of different libraries, including the distribution of transcriptional start sites (TSSs), were reproducible. The analysis of long-sized cDNAs showed that this library contained many cDNAs with a long-sized insert up to 11 199 bp of golgin B, including multiple slicing variants for filamin A and filamin B. These results suggest that the size-unbiased full-length cDNA library constructed using the vector-capping method will be an ideal resource for fine expression profiling of transcriptional variants with alternative TSSs and alternative splicing.

Comparative analysis of sequence features involved in the recognition of tandem splice sites

Background

The splicing of pre-mRNAs is conspicuously often variable and produces multiple alternatively spliced (AS) isoforms that encode different messages from one gene locus. Computational studies uncovered a class of highly similar isoforms, which were related to tandem 5'-splice sites (5'ss) and 3'-splice sites (3'ss), yet with very sparse anecdotal evidence in experimental studies. To compare the types and levels of alternative tandem splice site exons occurring in different human organ systems and cell types, and to study known sequence features involved in the recognition and distinction of neighboring splice sites, we performed large-scale, stringent alignments of cDNA sequences and ESTs to the human and mouse genomes, followed by experimental validation.

Results

We analyzed alternative 5'ss exons (A5Es) and alternative 3'ss exons (A3Es), derived from transcript sequences that were aligned to assembled genome sequences to infer patterns of AS occurring in several thousands of genes. Comparing the levels of overlapping (tandem) and non-overlapping (competitive) A5Es and A3Es, a clear preference of isoforms was seen for tandem acceptors and donors, with four nucleotides and three to six nucleotides long exon extensions, respectively. A subset of inferred A5E tandem exons was selected and experimentally validated. With the focus on A5Es, we investigated their transcript coverage, sequence conservation and base-paring to U1 snRNA, proximal and distal splice site classification, candidate motifs for cis-regulatory activity, and compared A5Es with A3Es, constitutive and pseudo-exons, in H. sapiens and M. musculus. The results reveal a small but authentic enriched set of tandem splice site preference, with specific distances between proximal and distal 5'ss (3'ss), which showed a marked dichotomy between the levels of in- and out-of-frame splicing for A5Es and A3Es, respectively, identified a number of candidate NMD targets, and allowed a rough estimation of a number of undetected tandem donors based on splice site information.

Conclusion

This comparative study distinguishes tandem 5'ss and 3'ss, with three to six nucleotides long extensions, as having unusually high proportions of AS, experimentally validates tandem donors in a panel of different human tissues, highlights the dichotomy in the types of AS occurring at tandem splice sites, and elucidates that human alternative exons spliced at overlapping 5'ss posses features of typical splice variants that could well be beneficial for the cell.

Entropy measures quantify global splicing disorders in cancer

Most mammalian genes are able to express several splice variants in a phenomenon known as alternative splicing. Serious alterations of alternative splicing occur in cancer tissues, leading to expression of multiple aberrant splice forms. Most studies of alternative splicing defects have focused on the identification of cancer-specific splice variants as potential therapeutic targets. Here, we examine instead the bulk of non-specific transcript isoforms and analyze their level of disorder using a measure of uncertainty called Shannon's entropy. We compare isoform expression entropy in normal and cancer tissues from the same anatomical site for different classes of transcript variations: alternative splicing, polyadenylation, and transcription initiation. Whereas alternative initiation and polyadenylation show no significant gain or loss of entropy between normal and cancer tissues, alternative splicing shows highly significant entropy gains for 13 of the 27 cancers studied. This entropy gain is characterized by a flattening in the expression profile of normal isoforms and is correlated to the level of estimated cellular proliferation in the cancer tissue. Interestingly, the genes that present the highest entropy gain are enriched in splicing factors. We provide here the first quantitative estimate of splicing disruption in cancer. The expression of normal splice variants is widely and significantly disrupted in at least half of the cancers studied. We postulate that such splicing disorders may develop in part from splicing alteration in key splice factors, which in turn significantly impact multiple target genes.

Comparative analysis of transposed element insertion within human and mouse genomes reveals Alu's unique role in shaping the human transcriptome

Analysis of transposed elements in the human and mouse genomes reveals many effects on the transcriptomes, including a higher level of exonization of Alu elements than other elements.

Background

Transposed elements (TEs) have a substantial impact on mammalian evolution and are involved in numerous genetic diseases. We compared the impact of TEs on the human transcriptome and the mouse transcriptome.

Results

We compiled a dataset of all TEs in the human and mouse genomes, identifying 3,932,058 and 3,122,416 TEs, respectively. We than extracted TEs located within human and mouse genes and, surprisingly, we found that 60% of TEs in both human and mouse are located in intronic sequences, even though introns comprise only 24% of the human genome. All TE families in both human and mouse can exonize. TE families that are shared between human and mouse exhibit the same percentage of TE exonization in the two species, but the exonization level of Alu, a primate-specific retroelement, is significantly greater than that of other TEs within the human genome, leading to a higher level of TE exonization in human than in mouse (1,824 exons compared with 506 exons, respectively). We detected a primate-specific mechanism for intron gain, in which Alu insertion into an exon creates a new intron located in the 3' untranslated region (termed 'intronization'). Finally, the insertion of TEs into the first and last exons of a gene is more frequent in human than in mouse, leading to longer exons in human.

Conclusion

Our findings reveal many effects of TEs on these two transcriptomes. These effects are substantially greater in human than in mouse, which is due to the presence of Alu elements in human.

Evidence of the modulation of mRNA splicing fidelity in humans by oxidative stress and p53

The majority of human genes generate mRNA splice variants and while there is little doubt that alternative splicing is an important biological phenomenon, it has also become apparent that some splice variants are associated with disease. To elucidate the molecular mechanisms responsible for generating aberrant splice variants, we have investigated alternative splicing of the human genes HPRT and POLB following oxidative stress in different genetic backgrounds. Our study revealed that splicing fidelity is sensitive to oxidative stress. Following treatment of cells with H2O2, the overall frequency of aberrant, unproductive splice variants increased in both loci. At least in POLB, splicing fidelity is p53 dependent. In the absence of p53, the frequency of POLB splice variants is elevated but oxidative stress does not further increase the frequency of splice variants. Our data indicate that mis-splicing following oxidative stress represents a novel and significant genotoxic outcome and that it is not simply DNA lesions induced by oxidative stress that lead to mis-splicing but changes in the alternative splicing machinery itself.

Bioinformatics detection of alternative splicing

In recent years, genome-wide detection of alternative splicing based on Expressed Sequence Tag (EST) sequence alignments with mRNA and genomic sequences has dramatically expanded our understanding of the role of alternative splicing in functional regulation. This chapter reviews the data, methodology, and technical challenges of these genome-wide analyses of alternative splicing, and briefly surveys some of the uses to which such alternative splicing databases have been put. For example, with proper alternative splicing database schema design, it is possible to query genome-wide for alternative splicing patterns that are specific to particular tissues, disease states (e.g., cancer), gender, or developmental stages. EST alignments can be used to estimate exon inclusion or exclusion level of alternatively spliced exons and evolutionary changes for various species can be inferred from exon inclusion level. Such databases can also help automate design of probes for RT-PCR and microarrays, enabling high throughput experimental measurement of alternative splicing.

SERpredict: detection of tissue- or tumor-specific isoforms generated through exonization of transposable elements

BACKGROUND

Transposed elements (TEs) are known to affect transcriptomes, because either new exons are generated from intronic transposed elements (this is called exonization), or the element inserts into the exon, leading to a new transcript. Several examples in the literature show that isoforms generated by an exonization are specific to a certain tissue (for example the heart muscle) or inflict a disease. Thus, exonizations can have negative effects for the transcriptome of an organism.

RESULTS

As we aimed at detecting other tissue- or tumor-specific isoforms in human and mouse genomes which were generated through exonization of a transposed element, we designed the automated analysis pipeline SERpredict (SER = Specific Exonized Retroelement) making use of Bayesian Statistics. With this pipeline, we found several genes in which a transposed element formed a tissue- or tumor-specific isoform.

CONCLUSION

Our results show that SERpredict produces relevant results, demonstrating the importance of transposed elements in shaping both the human and the mouse transcriptomes. The effect of transposed elements on the human transcriptome is several times higher than the effect on the mouse transcriptome, due to the contribution of the primate-specific Alu elements.

Discovery of tissue-specific exons using comprehensive human exon microarrays

BACKGROUND

Higher eukaryotes express a diverse population of messenger RNAs generated by alternative splicing. Large-scale methods for monitoring gene expression must adapt in order to accurately detect the transcript variation generated by this splicing.

RESULTS

We have designed a high-density oligonucleotide microarray with probesets for more than one million annotated and predicted exons in the human genome. Using these arrays and a simple algorithm that normalizes exon signal to signal from the gene as a whole, we have identified tissue-specific exons from a panel of 16 different normal adult tissues. RT-PCR validation confirms approximately 86% of the predicted tissue-enriched probesets. Pair-wise comparisons between the tissues suggest that as many as 73% of detected genes are differentially alternatively spliced. We also demonstrate how an inclusive exon microarray can be used to discover novel alternative splicing events. As examples, 17 new tissue-specific exons from 11 genes were validated by RT-PCR and sequencing.

CONCLUSION

In conjunction with a conceptually simple algorithm, comprehensive exon microarrays can detect tissue-specific alternative splicing events. Our data suggest significant expression outside of known exons and well annotated genes and a high frequency of alternative splicing events. In addition, we identified and validated a number of novel exons with tissue-specific splicing patterns. The tissue map data will likely serve as a valuable source of information on the regulation of alternative splicing.

Human MIP synthase splice variants in bipolar disorder

OBJECTIVES

Alternative splicing allows the production of multiple gene products with different functions from a given sequence, affecting cellular function control. Tissue-specific splicing is most prevalent in the brain. We therefore investigate whether splice variants contribute to complex psychiatric disorders. A database search suggested that the myo-inositol-1-phosphate (MIP) synthase gene, possibly involved in pathophysiology of bipolar disorder, has splice variants.

METHODS

Human RNA was purified from lymphocytes and postmortem brain. MIP synthase alternative splice variants were amplified using reverse transcription-polymerase chain reaction.

RESULTS

The bioinformatics finding was confirmed in both tissues. No difference in lymphocyte MIP synthase mRNA splice-variant levels was found between bipolar patients and controls. However, patients with family history of a major psychiatric disorder had significantly higher levels of the variant lacking exons 3 and 4 versus patients with no family history and controls.

CONCLUSIONS

As alternative splicing may be a mechanism by which the approximately 30,000 genes are amplified in mammalian brain, further studies with other candidate genes for psychiatric disorders are needed.

Global analysis of alternative splicing differences between humans and chimpanzees

Alternative splicing is a powerful mechanism affording extensive proteomic and regulatory diversity from a limited repertoire of genes. However, the extent to which alternative splicing has contributed to the evolution of primate species-specific characteristics has not been assessed previously. Using comparative genomics and quantitative microarray profiling, we performed the first global analysis of alternative splicing differences between humans and chimpanzees. Surprisingly, 6%-8% of profiled orthologous exons display pronounced splicing level differences in the corresponding tissues from the two species. Little overlap is observed between the genes associated with alternative splicing differences and the genes that display steady-state transcript level differences, indicating that these layers of regulation have evolved rapidly to affect distinct subsets of genes in humans and chimpanzees. The alternative splicing differences we detected are predicted to affect diverse functions including gene expression, signal transduction, cell death, immune defense, and susceptibility to diseases. Differences in expression at the protein level of the major splice variant of Glutathione S-transferase omega-2 (GSTO2), which functions in the protection against oxidative stress and is associated with human aging-related diseases, suggests that this enzyme is less active in human cells compared with chimpanzee cells. The results of this study thus support an important role for alternative splicing in establishing differences between humans and chimpanzees.

Identification and characterization of a new alpha-synuclein isoform and its role in Lewy body diseases

Alternative splicing is an important mechanism to generate a large number of mRNAs, thus increasing proteome diversity and tissue specificity. Three transcript variants of alpha-synuclein, a neuronal protein mainly involved in synapses, have been described so far. Whereas alpha-synuclein 140 is the whole and main transcript, alpha-synuclein 112 and 126 are short proteins that result from in-frame deletions of exons 3 and 5, respectively. Because the aforesaid alpha-synuclein isoforms show differential expression changes in Lewy body diseases (LBDs), in the present work, we searched for a fourth alpha-synuclein isoform and studied its expression levels in LBD brains. By using isoform-specific primers, isoform co-amplification and direct sequencing, we identified alpha-synuclein 98, which lacks exons 3 and 5. mRNA expression analyses in non-neuronal tissue revealed that alpha-synuclein 98 is a brain-specific splice variant with varying expression levels in different areas of fetal and adult brain. Additionally, we studied alpha-synuclein 98 expression levels by real-time semi-quantitative RT-PCR in the frontal cortices of LBD patients and compared them with those of Alzheimer disease (AD) patients and control subjects. Overexpression of alpha-synuclein 98 in LBD and AD brains would indicate its specific involvement in the pathogenesis of these neurodegenerative disorders.

Systematic analysis of alternative first exons in plant genomes

BACKGROUND

Alternative splicing (AS) contributes significantly to protein diversity, by selectively using different combinations of exons of the same gene under certain circumstances. One particular type of AS is the use of alternative first exons (AFEs), which can have consequences far beyond the fine-tuning of protein functions. For example, AFEs may change the N-termini of proteins and thereby direct them to different cellular compartments. When alternative first exons are distant, they are usually associated with alternative promoters, thereby conferring an extra level of gene expression regulation. However, only few studies have examined the patterns of AFEs, and these analyses were mainly focused on mammalian genomes. Recent studies have shown that AFEs exist in the rice genome, and are regulated in a tissue-specific manner. Our current understanding of AFEs in plants is still limited, including important issues such as their regulation, contribution to protein diversity, and evolutionary conservation.

RESULTS

We systematically identified 1,378 and 645 AFE-containing clusters in rice and Arabidopsis, respectively. From our data sets, we identified two types of AFEs according to their genomic organisation. In genes with type I AFEs, the first exons are mutually exclusive, while most of the downstream exons are shared among alternative transcripts. Conversely, in genes with type II AFEs, the first exon of one gene structure is an internal exon of an alternative gene structure. The functionality analysis indicated about half and approximately 19% of the AFEs in Arabidopsis and rice could alter N-terminal protein sequences, and approximately 5% of the functional alteration in type II AFEs involved protein domain addition/deletion in both genomes. Expression analysis indicated that 20-66% of rice AFE clusters were tissue- and/or development- specifically transcribed, which is consistent with previous observations; however, a much smaller percentage of Arabidopsis AFEs was regulated in this manner, which suggests different regulation mechanisms of AFEs between rice and Arabidopsis. Statistical analysis of some features of AFE clusters, such as splice-site strength and secondary structure formation further revealed differences between these two species. Orthologous search of AFE-containing gene pairs detected only 19 gene pairs conserved between rice and Arabidopsis, accounting only for a few percent of AFE-containing clusters.

CONCLUSION

Our analysis of AFE-containing genes in rice and Arabidopsis indicates that AFEs have multiple functions, from regulating gene expression to generating protein diversity. Comparisons of AFE clusters revealed different features in the two plant species, which indicates that AFEs may have evolved independently after the separation of rice (a model monocot) and Arabidopsis (a model dicot).

Identification of differentially regulated splice variants and novel exons in glial brain tumors using exon expression arrays

Aberrant splice variants are involved in the initiation and/or progression of glial brain tumors. We therefore set out to identify splice variants that are differentially expressed between histologic subgroups of gliomas. Splice variants were identified using a novel platform that profiles the expression of virtually all known and predicted exons present in the human genome. Exon-level expression profiling was done on 26 glioblastomas, 22 oligodendrogliomas, and 6 control brain samples. Our results show that Human Exon arrays can identify subgroups of gliomas based on their histologic appearance and genetic aberrations. We next used our expression data to identify differentially expressed splice variants. In two independent approaches, we identified 49 and up to 459 exons that are differentially spliced between glioblastomas and oligodendrogliomas, a subset of which (47% and 33%) were confirmed by reverse transcription-PCR (RT-PCR). In addition, exon level expression profiling also identified >700 novel exons. Expression of approximately 67% of these candidate novel exons was confirmed by RT-PCR. Our results indicate that exon level expression profiling can be used to molecularly classify brain tumor subgroups, can identify differentially regulated splice variants, and can identify novel exons. The splice variants identified by exon level expression profiling may help to detect the genetic changes that cause or maintain gliomas and may serve as novel treatment targets.

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.

Genome-wide detection of testis- and testicular cancer-specific alternative splicing

Alternative pre-messenger RNA (mRNA) splicing is a key molecular event that allows for protein diversity and plays important roles in development and disease. Alternative pre-mRNA splicing regulations during spermatogenesis and alternative pre-mRNA splicing etiology in testicular tumorigenesis are yet to be characterized. By genome-wide analysis, here we describe alternative splicing features that distinguish distinctive patterns of alternative pre-mRNA splicing among human testis, testicular cancer and mouse testis. Through computationally subtractive analysis, we detected 80 testis-specific transcript candidates in human testis, 175 in human testicular cancer and 262 in mouse testis, which were integrated into a database. Reverse transcription-polymerase chain reaction confirmed that most of these transcript candidates from mouse testis were testis specific. Around 40% of the transcripts were from unknown/hypothetical genes, which were useful for further functional analysis. These transcripts were not overlapped, indicating lack of evolutionary conservation. Further chromosome mapping showed distinct chromosomal preference of alternative pre-mRNA splicing events. Comparison analysis indicated that alternative pre-mRNA splicing in human testicular tumor shared some characters/trends with those in mouse testis. Moreover, human testicular tumor tended to use rare splice sites and there were also distinct sequences adjacent dominant splice sites between normal testis and testicular tumor. These special features of alternative pre-mRNA splicing in human testicular tumor suggested that testicular tumorigenesis was involved in multiple steps/levels of alternative splicing events. Using alternative splicing as a potential source for new clinical diagnostic, prognostic and therapeutic strategies for treatment of testicular tumors seems to have a bright prospect.

Alternative splicing events identified in human embryonic stem cells and neural progenitors

Human embryonic stem cells (hESCs) and neural progenitor (NP) cells are excellent models for recapitulating early neuronal development in vitro, and are key to establishing strategies for the treatment of degenerative disorders. While much effort had been undertaken to analyze transcriptional and epigenetic differences during the transition of hESC to NP, very little work has been performed to understand post-transcriptional changes during neuronal differentiation. Alternative RNA splicing (AS), a major form of post-transcriptional gene regulation, is important in mammalian development and neuronal function. Human ESC, hESC-derived NP, and human central nervous system stem cells were compared using Affymetrix exon arrays. We introduced an outlier detection approach, REAP (Regression-based Exon Array Protocol), to identify 1,737 internal exons that are predicted to undergo AS in NP compared to hESC. Experimental validation of REAP-predicted AS events indicated a threshold-dependent sensitivity ranging from 56% to 69%, at a specificity of 77% to 96%. REAP predictions significantly overlapped sets of alternative events identified using expressed sequence tags and evolutionarily conserved AS events. Our results also reveal that focusing on differentially expressed genes between hESC and NP will overlook 14% of potential AS genes. In addition, we found that REAP predictions are enriched in genes encoding serine/threonine kinase and helicase activities. An example is a REAP-predicted alternative exon in the SLK (serine/threonine kinase 2) gene that is differentially included in hESC, but skipped in NP as well as in other differentiated tissues. Lastly, comparative sequence analysis revealed conserved intronic cis-regulatory elements such as the FOX1/2 binding site GCAUG as being proximal to candidate AS exons, suggesting that FOX1/2 may participate in the regulation of AS in NP and hESC. In summary, a new methodology for exon array analysis was introduced, leading to new insights into the complexity of AS in human embryonic stem cells and their transition to neural stem cells.

Mapping of oxidative stress responses of human tumor cells following photodynamic therapy using hexaminolevulinate

Background

Photodynamic therapy (PDT) involves systemic or topical administration of a lesion-localizing photosensitizer or its precursor, followed by irradiation of visible light to cause singlet oxygen-induced damage to the affected tissue. A number of mechanisms seem to be involved in the protective responses to PDT, including activation of transcription factors, heat shock proteins, antioxidant enzymes and apoptotic pathways.

Results

In this study, we address the effects of a destructive/lethal hexaminolevulinate (HAL) mediated PDT dose on the transcriptome by using transcriptional exon evidence oligo microarrays. Here, we confirm deviations in the steady state expression levels of previously identified early defence response genes and extend this to include unreported PDT inducible gene groups, most notably the metallothioneins and histones. HAL-PDT mediated stress also altered expression of genes encoded by mitochondrial DNA (mtDNA). Further, we report PDT stress induced alternative splicing. Specifically, the ATF3 alternative isoform (deltaZip2) was up-regulated, while the full-length variant was not changed by the treatment. Results were independently verified by two different technological microarray platforms. Good microarray, RT-PCR and Western immunoblotting correlation for selected genes support these findings.

Conclusion

Here, we report new insights into how destructive/lethal PDT alters the transcriptome not only at the transcriptional level but also at post-transcriptional level via alternative splicing.

WNK4-mediated regulation of renal ion transport proteins

Point mutations in WNK4 [for With No K (lysine)], a serine-threonine kinase that is expressed in the distal nephron of the kidney, are linked to familial hyperkalemic hypertension (FHH). The imbalanced electrolyte homeostasis in FHH has led to studies toward an understanding of WNK4-mediated regulation of ion transport proteins in the kidney. A growing number of ion transport proteins for Na(+), K(+), Ca(2+), and Cl(-), including ion channels and transporters in the transcellular pathway and claudins in the paracellular pathway, are shown to be regulated by WNK4 from studies using models ranging from Xenopus laevis oocytes to transgenic and knockin mice. WNK4 regulates these transport proteins in different directions and by different cellular mechanisms. The common theme of WNK4-mediated regulation is to alter the abundance of ion transport proteins at the plasma membrane, with the exception of claudins, which are phosphorylated in the presence of WNK4. The regulation of WNK4 can be blocked by the full-length WNK1, whose action is in turn antagonized by a kidney-specific WNK1 variant lacking the kinase domain. In addition, WNK4 also activates stress-related serine-threonine kinases to regulate members of the SLC12 family members of cation-chloride cotransporters. In many cases, the FHH-causing mutants of WNK4 exhibit differences from wild-type WNK4 in regulating ion transport proteins. These regulations well explain the clinical features of FHH and provide insights into the multilayered regulation of ion transport processes in the distal nephron.

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.

A new advance in alternative splicing databases: from catalogue to detailed analysis of regulation of expression and function of human alternative splicing variants

Background

Most human genes produce several transcripts with different exon contents by using alternative promoters, alternative polyadenylation sites and alternative splice sites. Much effort has been devoted to describing known gene transcripts through the development of numerous databases. Nevertheless, owing to the diversity of the transcriptome, there is a need for interactive databases that provide information about the potential function of each splicing variant, as well as its expression pattern.

Description

After setting up a database in which human and mouse splicing variants were compiled, we developed tools (1) to predict the production of protein isoforms from these transcripts, taking account of the presence of open reading frames and mechanisms that could potentially eliminate transcripts and/or inhibit their translation, i.e. nonsense-mediated mRNA decay and microRNAs; (2) to support studies of the regulation of transcript expression at multiple levels, including transcription and splicing, particularly in terms of tissue specificity; and (3) to assist in experimental analysis of the expression of splicing variants. Importantly, analyses of all features from transcript metabolism to functional protein domains were integrated in a highly interactive, user-friendly web interface that allows the functional and regulatory features of gene transcripts to be assessed rapidly and accurately.

Conclusion

In addition to identifying the transcripts produced by human and mouse genes, fast DB provides tools for analyzing the putative functions of these transcripts and the regulation of their expression. Therefore, fast DB has achieved an advance in alternative splicing databases by providing resources for the functional interpretation of splicing variants for the human and mouse genomes. Because gene expression studies are increasingly employed in clinical analyses, our web interface has been designed to be as user-friendly as possible and to be readily searchable and intelligible at a glance by the whole biomedical community.

Biological Cross-talk between WNK1 and the Transforming Growth Factor β-Smad Signaling Pathway

WNKs (with no lysine (K)), unique serine/threonine protein kinases, have been best studied in the context of cell volume regulation and ion homeostasis. Here we describe a biological link between WNKs and transforming growth factor (TGF) beta-Smad signaling. Both WNK1 and WNK4 directly bind to and phosphorylate Smad2. Knockdown of WNK1 in HeLa cells using small interfering RNA reduces Smad2 protein expression; this decrease is at least partially due to down-regulation of Smad2 transcription. In contrast, phosphorylated Smad2 significantly accumulated in the nucleus as a consequence of depletion of WNK1, resulting in Smad-mediated transcriptional responses. In addition, TGFbeta-induced target gene transcripts were increased in WNK1 small interfering RNA cells. These findings suggest WNK1 as a dual modulator of TGFbeta-Smad signaling pathways.

Dynamic usage of transcription start sites within core promoters

BACKGROUND

Mammalian promoters do not initiate transcription at single, well defined base pairs, but rather at multiple, alternative start sites spread across a region. We previously characterized the static structures of transcription start site usage within promoters at the base pair level, based on large-scale sequencing of transcript 5' ends.

RESULTS

In the present study we begin to explore the internal dynamics of mammalian promoters, and demonstrate that start site selection within many mouse core promoters varies among tissues. We also show that this dynamic usage of start sites is associated with CpG islands, broad and multimodal promoter structures, and imprinting.

CONCLUSION

Our results reveal a new level of biologic complexity within promoters--fine-scale regulation of transcription starting events at the base pair level. These events are likely to be related to epigenetic transcriptional regulation.

SP1 transcription factors in male germ cell development and differentiation

Transcription factor SP1 is a zinc finger protein that has been implicated in regulating the expression of several genes involved in cellular differentiation and embryonic development. The zinc finger region of SP1 transcription factors binds to GC or GT-box elements present in the promoters of a number of male germ cell target genes that are developmentally expressed during spermatogenesis. The glutamine and serine/threonine-rich regions of the SP1 proteins recruit co-regulatory factors to the multi-protein preinitiation complex that are important for mediating transcriptional activation in male germ cells. Studies in our laboratory have identified several alternatively spliced transcripts encoding SP1 isoforms that display stage and cell-type-specific expression profiles in differentiating germ cells in the seminiferous epithelium of the testis. This review summarizes the expression patterns and functional significance of these SP1 transcription factor variants during spermatogenesis.

Virtual Northern analysis of the human genome

Background

We applied the Virtual Northern technique to human brain mRNA to systematically measure human mRNA transcript lengths on a genome-wide scale.

Methodology/Principal Findings

We used separation by gel electrophoresis followed by hybridization to cDNA microarrays to measure 8,774 mRNA transcript lengths representing at least 6,238 genes at high (>90%) confidence. By comparing these transcript lengths to the Refseq and H-Invitational full-length cDNA databases, we found that nearly half of our measurements appeared to represent novel transcript variants. Comparison of length measurements determined by hybridization to different cDNAs derived from the same gene identified clones that potentially correspond to alternative transcript variants. We observed a close linear relationship between ORF and mRNA lengths in human mRNAs, identical in form to the relationship we had previously identified in yeast. Some functional classes of protein are encoded by mRNAs whose untranslated regions (UTRs) tend to be longer or shorter than average; these functional classes were similar in both human and yeast.

Conclusions/Significance

Human transcript diversity is extensive and largely unannotated. Our length dataset can be used as a new criterion for judging the completeness of cDNAs and annotating mRNA sequences. Similar relationships between the lengths of the UTRs in human and yeast mRNAs and the functions of the proteins they encode suggest that UTR sequences serve an important regulatory role among eukaryotes.

Discovery and analysis of evolutionarily conserved intronic splicing regulatory elements

Knowledge of the functional cis-regulatory elements that regulate constitutive and alternative pre-mRNA splicing is fundamental for biology and medicine. Here we undertook a genome-wide comparative genomics approach using available mammalian genomes to identify conserved intronic splicing regulatory elements (ISREs). Our approach yielded 314 ISREs, and insertions of ~70 ISREs between competing splice sites demonstrated that 84% of ISREs altered 5′ and 94% altered 3′ splice site choice in human cells. Consistent with our experiments, comparisons of ISREs to known splicing regulatory elements revealed that 40%–45% of ISREs might have dual roles as exonic splicing silencers. Supporting a role for ISREs in alternative splicing, we found that 30%–50% of ISREs were enriched near alternatively spliced (AS) exons, and included almost all known binding sites of tissue-specific alternative splicing factors. Further, we observed that genes harboring ISRE-proximal exons have biases for tissue expression and molecular functions that are ISRE-specific. Finally, we discovered that for Nova1, neuronal PTB, hnRNP C, and FOX1, the most frequently occurring ISRE proximal to an alternative conserved exon in the splicing factor strongly resembled its own known RNA binding site, suggesting a novel application of ISRE density and the propensity for splicing factors to auto-regulate to associate RNA binding sites to splicing factors. Our results demonstrate that ISREs are crucial building blocks in understanding general and tissue-specific AS regulation and the biological pathways and functions regulated by these AS events.

During RNA splicing, sequences (introns) in a pre-mRNA are excised and discarded, and the remaining sequences (exons) are joined to form the mature RNA. Splicing is regulated not only by the binding of the basic splicing machinery to splice sites located at the exon–intron boundaries, but also by the combined effects of various other splicing factors that bind to a multitude of sequence elements located both in the exons as well as the flanking introns. Instances of alternative splicing, where usage of splice site(s) is incomplete or different between tissues, cell types, or lineages, can be created by the interaction of sequence elements and tissue, cell type, and stage-specific splicing factors. To better understand constitutive and alternative pre-mRNA splicing, the authors describe a comparative genomics approach, using available mammalian genomes, to systematically identify splicing regulatory elements located in the introns proximal to exons. A quarter of the elements were tested experimentally, and most of them altered splicing in human cells. The authors also showed that that the intronic elements are close to tissue-specific alternative exons and are more likely to be located in specific positions in the introns, suggestive of potential regulatory function. These elements are also frequently found in tissue-specific genes, suggesting a coupling between expression and alternative splicing of these genes. Finally, the authors propose a strategy using the elements to identify the binding sites of several splicing factors.

Plant SET domain-containing proteins: structure, function and regulation

Modification of the histone proteins that form the core around which chromosomal DNA is looped has profound epigenetic effects on the accessibility of the associated DNA for transcription, replication and repair. The SET domain is now recognized as generally having methyltransferase activity targeted to specific lysine residues of histone H3 or H4. There is considerable sequence conservation within the SET domain and within its flanking regions. Previous reviews have shown that SET proteins from Arabidopsis and maize fall into five classes according to their sequence and domain architectures. These classes generally reflect specificity for a particular substrate. SET proteins from rice were found to fall into similar groupings, strengthening the merit of the approach taken. Two additional classes, VI and VII, were established that include proteins with truncated/interrupted SET domains. Diverse mechanisms are involved in shaping the function and regulation of SET proteins. These include protein-protein interactions through both intra- and inter-molecular associations that are important in plant developmental processes, such as flowering time control and embryogenesis. Alternative splicing that can result in the generation of two to several different transcript isoforms is now known to be widespread. An exciting and tantalizing question is whether, or how, this alternative splicing affects gene function. For example, it is conceivable that one isoform may debilitate methyltransferase function whereas the other may enhance it, providing an opportunity for differential regulation. The review concludes with the speculation that modulation of SET protein function is mediated by antisense or sense-antisense RNA.

Orengo C, Patthy L, Thornton JM, Tramontano A, Valencia A. 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.

Alternative isoform discrimination by the next generation of expression profiling microarrays

Microarray expression profiling has revolutionised the way that many therapeutic targets have been identified over the past 10 years. High-density microarrays have allowed scientists to simultaneously scrutinise the expression of all genes encoded on a given genome. Although the data collected from classically designed microarrays greatly enriched the pool of information available to help guide the selection and design of new therapeutic strategies, they were unable to tell the complete story. The major limitation with most array designs is that they can only produce a global expression value for all transcripts produced from a specific locus and cannot monitor each individual alternative isoform produced from the interrogated locus. Recently, new array designs have been described, and become commercially available, that can efficiently monitor individual alternatively spliced isoforms produced from a single locus, allowing the research community to get a more accurate picture of the biological landscape of the expressed transcripts.

Natural genome-editing competences of viruses

It is becoming increasingly evident that the driving forces of evolutionary novelty are not randomly derived chance mutations of the genetic text, but a precise genome editing by omnipresent viral agents. These competences integrate the whole toolbox of natural genetic engineering, replication, transcription, translation, genomic imprinting, genomic creativity, enzymatic inventions and all types of genetic repair patterns. Even the non-coding, repetitive DNA sequences which were interpreted as being ancient remnants of former evolutionary stages are now recognized as being of viral descent and crucial for higher-order regulatory and constitutional functions of protein structural vocabulary. In this article I argue that non-randomly derived natural genome editing can be envisioned as (a) combinatorial (syntactic), (b) context-specific (pragmatic) and (c) content-sensitive (semantic) competences of viral agents. These three-leveled biosemiotic competences could explain the emergence of complex new phenotypes in single evolutionary events. After short descriptions of the non-coding regulatory networks, major viral life strategies and pre-cellular viral life three of the major steps in evolution serve as examples: There is growing evidence that natural genome-editing competences of viruses are essential (1) for the evolution of the eukaryotic nucleus, (2) the adaptive immune system and (3) the placental mammals.

SpliceMiner: a high-throughput database implementation of the NCBI Evidence Viewer for microarray splice variant analysis

Background

There are many fewer genes in the human genome than there are expressed transcripts. Alternative splicing is the reason. Alternatively spliced transcripts are often specific to tissue type, developmental stage, environmental condition, or disease state. Accurate analysis of microarray expression data and design of new arrays for alternative splicing require assessment of probes at the sequence and exon levels.

Description

SpliceMiner is a web interface for querying Evidence Viewer Database (EVDB). EVDB is a comprehensive, non-redundant compendium of splice variant data for human genes. We constructed EVDB as a queryable implementation of the NCBI Evidence Viewer (EV). EVDB is based on data obtained from NCBI Entrez Gene and EV. The automated EVDB build process uses only complete coding sequences, which may or may not include partial or complete 5' and 3' UTRs, and filters redundant splice variants. Unlike EV, which supports only one-at-a-time queries, SpliceMiner supports high-throughput batch queries and provides results in an easily parsable format. SpliceMiner maps probes to splice variants, effectively delineating the variants identified by a probe.

Conclusion

EVDB can be queried by gene symbol, genomic coordinates, or probe sequence via a user-friendly web-based tool we call SpliceMiner (). The EVDB/SpliceMiner combination provides an interface with human splice variant information and, going beyond the very valuable NCBI Evidence Viewer, supports fluent, high-throughput analysis. Integration of EVDB information into microarray analysis and design pipelines has the potential to improve the analysis and bioinformatic interpretation of gene expression data, for both batch and interactive processing. For example, whenever a gene expression value is recognized as important or appears anomalous in a microarray experiment, the interactive mode of SpliceMiner can be used quickly and easily to check for possible splice variant issues.

The testis-specific human protein RBMY recognizes RNA through a novel mode of interaction

The RBMY (RNA-binding motif gene on Y chromosome) protein encoded by the human Y chromosome is important for normal sperm development. Although its precise molecular RNA targets are unknown at present, it is suggested that human RBMY (hRBMY) participates in splicing in the testis. Using systematic evolution of ligands by exponential enrichment, we found that RNA stem-loops capped by a C(A)/(U)CAA pentaloop are high-affinity binding targets for hRBMY. Subsequent nuclear magnetic resonance structural determination of the hRBMY RNA recognition motif (RRM) in complex with a high-affinity target showed two distinct modes of RNA recognition. First, the RRM beta-sheet surface binds to the RNA loop in a sequence-specific fashion. Second, the beta2-beta3 loop of the hRBMY inserts into the major groove of the RNA stem. The first binding mode might be conserved in the paralogous protein heterogeneous nuclear RNP G, whereas the second mode of binding is found only in hRBMY. This structural difference could be at the origin of the function of RBMY in spermatogenesis.

Relating alternative splicing to proteome complexity and genome evolution

Prior to genomics, studies of alternative splicing primarily focused on the function and mechanism of alternative splicing in individual genes and exons. This has changed dramatically since the late 1990s. High-throughput genomics technologies, such as EST sequencing and microarrays designed to detect changes in splicing, led to genome-wide discoveries and quantification of alternative splicing in a wide range of species from human to Arabidopsis. Consensus estimates of AS frequency in the human genome grew from less than 5% in mid-1990s to as high as 60-74% now. The rapid growth in sequence and microarray data for alternative splicing has made it possible to look into the global impact of alternative splicing on protein function and evolution of genomes. In this chapter, we review recent research on alternative splicing's impact on proteomic complexity and its role in genome evolution.

Mammalian mRNA splice-isoform selection is tightly controlled

Post-transcriptional RNA processing is an important regulatory control mechanism for determining the phenotype of eukaryotic cells. The processing of a transcribed RNA species into alternative splice isoforms yields products that can perform different functions. Each type of cell in a multi-cellular organism is presumed to actively control the relative quantities of alternative splice isoforms. In this study, the alternatively spliced isoforms of five mRNA transcription units were examined by quantitative reverse transcription-PCR amplification. We show that interindividual variation in splice-isoform selection is very highly constrained when measured in a large population of genetically diverse mice (i.e., full siblings; N = 150). Remarkably, splice-isoform ratios are among the most invariant phenotypes measured in this population and are confirmed in a second, genetically distinct population. In addition, the patterns of splice-isoform selection show tissue-specific and age-related changes. We propose that splice-isoform selection is exceptionally robust to genetic and environmental variability and may provide a control point for cellular homeostasis. As a consequence, splice-isoform ratios may be useful as a practical quantitative measure of the physiological status of cells and tissues.