- Open Access
Variation in alternative splicing across human tissues
© Yeo et al.; licensee BioMed Central Ltd. 2004
Received: 19 April 2004
Accepted: 27 July 2004
Published: 13 September 2004
Alternative pre-mRNA splicing (AS) is widely used by higher eukaryotes to generate different protein isoforms in specific cell or tissue types. To compare AS events across human tissues, we analyzed the splicing patterns of genomically aligned expressed sequence tags (ESTs) derived from libraries of cDNAs from different tissues.
Controlling for differences in EST coverage among tissues, we found that the brain and testis had the highest levels of exon skipping. The most pronounced differences between tissues were seen for the frequencies of alternative 3' splice site and alternative 5' splice site usage, which were about 50 to 100% higher in the liver than in any other human tissue studied. Quantifying differences in splice junction usage, the brain, pancreas, liver and the peripheral nervous system had the most distinctive patterns of AS. Analysis of available microarray expression data showed that the liver had the most divergent pattern of expression of serine-arginine protein and heterogeneous ribonucleoprotein genes compared to the other human tissues studied, possibly contributing to the unusually high frequency of alternative splice site usage seen in liver. Sequence motifs enriched in alternative exons in genes expressed in the brain, testis and liver suggest specific splicing factors that may be important in AS regulation in these tissues.
This study distinguishes the human brain, testis and liver as having unusually high levels of AS, highlights differences in the types of AS occurring commonly in different tissues, and identifies candidate cis-regulatory elements and trans-acting factors likely to have important roles in tissue-specific AS in human cells.
The differentiation of a small number of cells in the developing embryo into the hundreds of cell and tissue types present in a human adult is associated with a multitude of changes in gene expression. In addition to many differences between tissues in transcriptional and translational regulation of genes, alternative pre-mRNA splicing (AS) is also frequently used to regulate gene expression and to generate tissue-specific mRNA and protein isoforms [1–5]. Between one-third and two-thirds of human genes are estimated to undergo AS [6–11] and the disruption of specific AS events has been implicated in several human genetic diseases . The diverse and important biological roles of alternative splicing have led to significant interest in understanding its regulation.
Insights into the regulation of AS have come predominantly from the molecular dissection of individual genes (reviewed in [1, 12]). Prominent examples include the tissue-specific splicing of the c-src N1 exon , cancer-associated splicing of the CD44 gene  and the alternative splicing cascade involved in Drosophila melanogaster sex determination . Biochemical studies of these and other genes have described important classes of trans-acting splicing-regulatory factors, implicating members of the ubiquitously expressed serine/arginine-rich protein (SR protein) and heterogeneous nuclear ribonucleoprotein (hnRNP) families, and tissue-specific factors including members of the CELF  and NOVA  families of proteins, as well as other proteins and protein families, in control of specific splicing events. A number of cis-regulatory elements in exons or introns that play key regulatory roles have also been identified, using a variety of methods including site-directed mutagenesis, systematic evolution of ligands by exponential enrichment (SELEX) and computational approaches [18–22]. In addition, DNA microarrays and polymerase colony approaches have been developed for higher-throughput analysis of alternative mRNA isoforms [23–26] and a cross-linking/immunoprecipitation strategy (CLIP) has been developed for systematic detection of the RNAs bound by a given splicing factor . These new methods suggest a path towards increasingly parallel experimental analysis of splicing regulation.
From another direction, the accumulation of large databases of cDNA and expressed sequence tag (EST) sequences has enabled large-scale computational studies, which have assessed the scope of AS in the mammalian transcriptome [3, 8, 10, 28]. Other computational studies have analyzed the tissue specificity of AS events and identified sets of exons and genes that exhibit tissue-biased expression [29, 30]. However, a number of significant questions about tissue-specific alternative splicing have not yet been comprehensively addressed. Which tissues have the highest and lowest proportions of alternative splicing? Do tissues differ in their usage of different AS types, such as exon skipping, alternative 5' splice site choice or alternative 3' splice site choice? Which tissues are most distinct from other tissues in the spectrum of alternative mRNA isoforms they express? And to what extent do expression levels of known splicing factors explain AS patterns in different tissues?
Here, we describe an initial effort to answer these questions using a large-scale computational analysis of ESTs derived from about two dozen human tissues, which were aligned to the assembled human genome sequence to infer patterns of AS occurring in thousands of human genes. Our results distinguish specific tissues as having high levels and distinctive patterns of AS, identify pronounced differences between the proportions of alternative 5' splice site and alternative 3' splice site usage between tissues, and predict candidate cis-regulatory elements and trans-acting factors involved in tissue-specific AS.
Results and discussion
Variation in the levels of alternative splicing in different human tissues
To determine whether differences occur in the proportions of these three types of AS events across human tissues, we assessed the frequencies of genes containing skipped exons, alternative 3' splice site exons or alternative 5' splice site exons for 16 human tissues (see Figure 1 for the list of tissues) for which sufficiently large numbers of EST sequences were available. Because the availability of a larger number of ESTs derived from a gene increases the chance of observing alternative isoforms of that gene, the proportion of AS genes observed in a tissue will tend to increase with increasing EST coverage of genes [10, 31]. Since the number of EST sequences available differs quite substantially among human tissues (for example, the dbEST database contains about eight times more brain-derived ESTs than heart-derived ESTs), in order to compare the proportion of AS in different tissues in an unbiased way, we used a sampling strategy that ensured that all genes/tissues studied were represented by equal numbers of ESTs.
It is important to point out that our analysis does not make use of the concept of a canonical transcript for each gene because it is not clear that such a transcript could be chosen objectively or that this concept is biologically meaningful. Instead, AS events are defined only through pairwise comparison of ESTs.
Our objective was to control for differences in EST abundance across tissues while retaining sufficient power to detect a reasonable fraction of AS events. For each tissue we considered genes that had at least 20 aligned EST sequences derived from human cDNA libraries specific to that tissue ('tissue-derived' ESTs). For each such gene, a random sample of 20 of these ESTs was chosen (without replacement) to represent the splicing of the given gene in the given human tissue. For the gene and tissue combinations included in this analysis, the median number of EST sequences per gene was not dramatically different between tissues, ranging from 25 to 35 (see Additional data file 1). The sampled ESTs for each gene were then compared to each other to identify AS events occurring within the given tissue (see Materials and methods). The random sampling was repeated 20 times and the mean fraction of AS genes observed in these 20 trials was used to assess the fraction of AS genes for each tissue (Figure 1a). Different random subsets of a relatively large pool will have less overlap in the specific ESTs chosen (and therefore in the specific AS events detected) than for random subsets of a smaller pool of ESTs, and increased numbers of ESTs give greater coverage of exons. However, there is no reason that the expected number of AS events detected per randomly sampled subset should depend on the size of the pool the subset was chosen from. While the error (standard deviation) of the measured AS frequency per gene should be lower when restricting to genes with larger minimum pools of ESTs, such a restriction would not change the expected value. Unfortunately, the reduction in error of the estimated AS frequency per gene is offset by an increase in the expected error of the tissue-level AS frequency resulting from the use of fewer genes. The inclusion of all genes with at least 20 tissue-derived ESTs represents a reasonable trade-off between these factors.
The human brain had the highest fraction of AS genes in this analysis (Figure 1a), with more than 40% of genes exhibiting one or more AS events, followed by the liver and testis. Previous EST-based analyses have identified high proportions of splicing in human brain and testis tissues [29, 30, 32]. These studies did not specifically control for the highly unequal representation of ESTs from different human tissues. As larger numbers of ESTs increase the chance of observing a larger fraction of the expressed isoforms of a gene, the number of available ESTs has a direct impact on estimated proportions of AS, as seen previously in analyses comparing the levels of AS in different organisms . Thus, the results obtained in this study confirm that the human brain and testis possess an unusually high level of AS, even in the absence of EST-abundance advantages over other tissues. We also observe a high level of AS in the human liver, a tissue with much lower EST coverage, where higher levels of AS have been previously reported in cancerous cells [33, 34]. Human muscle, uterus, breast, stomach and pancreas had the lowest levels of AS genes in this analysis (less than 25% of genes). Lowering the minimum EST count for inclusion in this analysis from 20 to 10 ESTs, and sampling 10 (out of 10 or more) ESTs to represent each gene in each tissue, did not alter the results qualitatively (data not shown).
Differences in the levels of exon skipping in different tissues
Alternatively spliced genes in this analysis exhibited on average between one and two distinct AS exons. Analyzing the different types of AS events separately, we found that the human brain and testis had the highest levels of skipped exons, with more than 20% of genes containing SEs (Figure 1b). The high level of skipped exons observed in the brain is consistent with previous analyses [29, 30, 32]. At the other extreme, the human ovary, muscle, uterus and liver had the lowest levels of skipped exons (about 10% of genes).
Differences in the levels of alternative splice site usage in different tissues
Analyzing the proportions of AS events involving the usage of A5Es and A3Es revealed a very different pattern (Figure 1c,d). Notably, the fraction of genes containing A3Es was more than twice as high in the liver as in any other human tissue studied (Figure 1d), and the level of A5Es was also about 40-50% higher in the liver than in any other tissue (Figure 1c). The tissue with the second highest level of alternative usage for both 5' splice sites and 3' splice sites was the brain. Another group of human tissues including muscle, uterus, breast, pancreas and stomach - similar to the low SE frequency group above - had the lowest level of A5Es and A3Es (less than 5% of genes in each category). Thus, a picture emerges in which certain human tissues such as muscle, uterus, breast, pancreas and stomach, have low levels of AS of all types, whereas other tissues, such as the brain and testis, have relatively high levels of AS of all types and the liver has very high levels of A3Es and A5Es, but exhibits only a modest level of exon skipping. To our knowledge, this study represents the first systematic analysis of the proportions of different types of AS events occurring in different tissues. Repeating the analyses by removing ESTs from disease-associated tissue libraries, using available library classifications , gave qualitatively similar results (see Additional data files 2, 3, and 4). These data show that ESTs derived from diseased tissues show modestly higher frequencies of exon skipping, but the relative rankings of tissues remain similar. The fractions of genes containing A5Es and A3Es were not changed substantially when diseased-tissue ESTs were excluded.
Human genes expressed in the liver with alternative 3' splice site exons (A3Es) or alternative 5' splice site exons (A5Es)
Ensembl gene ID
Fold-change above median expression, HG-U95A
Fold-change above median expression, MG-U74A
Serotransferrin precursor, TF
8, 9; 4
Fibronectin precursor, FN1
Antithrombin-III precursor, SERPINC1
Fructose-bisphosphate aldolase, ALDOB
3, 8; 4
Haptoglobin-related protein precursor, HPR
Tryptophan 2,3-dioxygenase, TDO2
3, 5; 4
Ceruloplasmin precursor, CP
Inter-alpha-trypsin inhibitor heavy chain H1 precursor, ITIH1
Aldehyde dehydrogenase, ALDH2
Pigment epithelium-derived factor precursor, SERPINF1
Aldehyde oxidase, AOX1
Isocitrate dehydrogenase, IDH1
Betaine-homocysteine S-methyltransferase, BHMT
Cytochrome P450, CYP3A4
Glycine amidinotransferase, GATM
Heat shock protein HSP 90-alpha, HSPCA
Glutathione S-transferase, GSTA2
Protein precursor, AMBP
Telomerase-binding protein P23, TEBP
Hydroxymethylglutaryl-CoA synthase, HMGCS2
Cytochrome P450, CYP2C8
Alpha-2-HS-glycoprotein precursor, AHSG
Serum albumin precursor, ALB
Fibrinogen gamma chain precursor, FGG
Glutathione S-transferase, GSTA3
Sequences and splicing patterns for two of these genes for which orthologous mouse exons/genes and transcripts could be identified - the genes BHMT and CYP2C8 - are shown in detail in Figure 2b,c. In the event depicted for BHMT, the exons involved are highly conserved between the human and mouse orthologs (Figure 2b), consistent with the possibility that the splicing event may have a (conserved) regulatory role. This AS event preserves the reading frame of downstream exons, so the two isoforms are both likely to produce functional proteins, differing by the insertion/deletion of 23 amino acids. In the event depicted for CYP2C8, usage of an alternative 3' splice site removes 71 nucleotides, shifting the reading frame and leading to a premature termination codon in the exon (Figure 2c). In this case, the shorter alternative transcript is a potential substrate for nonsense-mediated decay [41, 42] and the AS event may be used to regulate the level of functional mRNA/protein produced.
Differences in splicing factor expression between tissues
While most of the tissues examined showed a very high degree of correlation in the expression levels of the 20 splicing factors studied (typically with r > 0.75; Figure 3), the human adult liver was clearly an outlier, with low concordance in splicing-factor expression to most other tissues (typically r < 0.6, and often much lower). The unusual splicing-factor expression in the human liver was seen consistently in data from two independent DNA microarray studies using different probe sets (compare the two halves of Figure 3). The low correlation observed between liver and other tissues in splicing factor expression is statistically significant even relative to arbitrary collections of 20 genes (see Additional data file 8). Examining the relative levels of specific splicing factors in the human adult liver versus other tissues, the relative level of SRp30c message was consistently higher in the liver and the relative levels of SRp40, hnRNP A2/B2 and Srp54 messages were consistently lower. A well established paradigm in the field of RNA splicing is that usage of alternative splice sites is often controlled by the relative concentrations of specific SR proteins and hnRNP proteins [49–52]. This functional antagonism between particular SR and hnRNP proteins is often due to competition for binding of nearby sites on pre-mRNAs [49, 53, 54]. Therefore, it seems likely that the unusual patterns of expression seen in the human adult liver for these families of splicing factors may contribute to the high level of alternative splice site usage seen in this tissue. It is also interesting that splicing-factor expression in the human fetal liver is highly concordant with most other tissues, but has low concordance with the adult liver (Figure 3). This observation suggests that substantial changes in splicing-factor expression may occur during human liver development, presumably leading to a host of changes in the splicing patterns of genes expressed in human liver. Currently available EST data were insufficient to allow systematic analysis of the patterns of AS in fetal relative to adult liver.
An important caveat to these results is that the DNA microarray data used in this analysis measure mRNA expression levels rather than protein levels or activities. The relation between the amount of mRNA expressed from a gene and the concentration of the corresponding protein has been examined previously in several studies in yeast as well as in human and mouse liver tissues [55–58]. These studies have generally found that mRNA expression levels correlate positively with protein concentrations, but with fairly wide divergences for a significant fraction of genes.
Over-represented motifs in alternative exons in the human brain, testis and liver
Sequence motifs enriched in skipped exons (SEs) and alternative 5' splice site exons (A5Es)
AS type /tissue (motif name)
Consensus (% of exons containing)
C CUCC C
G UAGC U
U UAGG G
Core A5E/liver (LI1)
Core A5E/liver (LI2)
A comparison of human testis-derived skipped exons to exons constitutively included in genes expressed in the testis identified only a single cluster of sequences, TE1, which share the tetramer UAGG. Enrichment of this motif, common to the brain-specific cluster BR7, suggests a role for regulation of exon skipping by hnRNP A1 - or a trans-acting factor with similar binding preferences - in the testis.
Alternative splice site usage gives rise to two types of exon segments - the 'core' portion common to both splice forms and the 'extended' portion that is present only in the longer isoform. Two clusters of sequence motifs enriched in the core sequences of A5Es in genes expressed in the liver relative to the core segments of A5Es resulting from alignments of non-liver-derived ESTs were identified - LI1 and LI2. Both are adenosine-rich, with consensus tetramers AAAC and UAAA, respectively. The former motif matches a candidate ESE motif identified previously using the computational/experimental RESCUE-ESE approach (motif 3F with consensus [AG]AA [AG]C) . The enrichment of a probable ESE motif in exons exhibiting alternative splice site usage in the liver is consistent with the model that such splicing events are often controlled by the relative levels of SR proteins (which bind many ESEs) and hnRNP proteins. Insufficient data were available for the analysis of motifs in the extended portions of liver A5Es (which tend to be significantly shorter than the core regions) or for the analysis of liver A3Es.
A measure of dissimilarity between mRNA isoforms
The SJD value can be generalized to compare splicing patterns between two sets of transcripts from a gene - for example, to compare the splicing patterns of the sets of ESTs derived from two different tissues. In this case, the SJD is defined by counting the number of splice junctions that differ between all pairs of transcripts (i, j), with transcript i coming from set 1 (for example, heart-derived ESTs), and transcript j coming from set 2 (for example, lung-derived ESTs), and dividing this number by the total number of splice junctions in all pairs of transcripts compared, again considering only those splice junctions that occur in regions of overlap between the transcript pairs considered. Note that this definition has the desirable property that pairs of transcripts that have larger numbers of overlapping splice junctions contribute more to the total than transcript pairs that overlap less. As an example of the splice junction difference between two sets of transcripts, consider the set S1, consisting of transcripts (1,2) from Figure 4, and set S2, consisting of transcripts (3,4) from Figure 4. Using the notation introduced in Figure 4, SJD(S1,S2) = d(S1,S2) / t(S1,S2) = [d(1,3) + d(1,4) + d(2,3) + d(2,4)]/ [t(1,3) +t(1,4) + t(2,3) + t(2,4)] = [3 + 4 + 2 + 3]/ [3 + 4 + 4 + 5] = 12/16 = 0.75, reflecting a high level of dissimilarity between the isoforms in these sets, whereas the SJD falls to 0.57 for the more similar sets S1 = transcripts (1,2) versus S3 = transcripts (2,3). Note that in cases where multiple similar/identical transcripts occur in a given set, the SJD measure effectively weights the isoforms by their abundance, reflecting an average dissimilarity when comparing randomly chosen pairs of transcripts from the two tissues. For example, the SJD computed for the set S4 = (1,2,2,2,2), that is, one transcript aligning as transcript 1 in Figure 4 and four transcripts aligning as transcript 2, and the set S5 = (2,2,2,2,3) is 23/95 = 0.24, substantially lower than the SJD value for sets S1 versus S3 above, reflecting the higher fraction of identically spliced transcripts between sets S4 and S5.
Global comparison of splicing patterns between tissues
To make a global comparison of patterns of splicing between two different human tissues, a tissue-level SJD value was computed by comparing the splicing patterns of ESTs from all genes for which at least one EST was available from cDNA libraries representing both tissues. The 'inter-tissue' SJD value is then defined as the ratio of the sum of d(SA,SB) values for all such genes, divided by the sum of t(SA,SB) values for all of these genes, where SA and SB refer to the set of ESTs for a gene derived from tissues A and B, respectively, and d(SA,SB) and t(SA,SB) are defined in terms of comparison of all pairs of ESTs from the two sets as described above. This analysis uses all available ESTs for each gene in each tissue (rather than samples of a fixed size). A large SJD value between a pair of tissues indicates that mRNA isoforms of genes expressed in the two tissues tend to be more dissimilar in their splicing patterns than is the case for two tissues with a smaller inter-tissue SJD value. This definition puts greater weight on those genes for which more ESTs are available.
The systematic analysis of transcripts generated from the human genome is just beginning, but promises to deepen our understanding of how changes in the program of gene expression contribute to development and differentiation. Here, we have observed pronounced differences between human tissues in the set of alternative mRNA isoforms that they express. Because our approach normalizes the EST coverage per gene in each tissue, there is higher confidence that these differences accurately reflect differences in splicing patterns between tissues. As human tissues are generally made up of a mixture of cell types, each of which may have its own unique pattern of gene expression and splicing, it will be important in the future to develop methods for systematic analysis of transcripts in different human cell types.
Understanding the mechanisms and regulatory consequences of AS will require experimental and computational analyses at many levels. At its core, AS involves the generation of alternative transcripts mediated by interactions between cis-regulatory elements in exons or introns and trans-acting splicing factors. The current study has integrated these three elements, inferring alternative transcripts from EST-genomic alignments, identifying candidate regulatory sequence motifs enriched in alternative exons from different tissues, and analyzing patterns of splicing-factor expression in different tissues. Our results emphasize differences in the frequencies of exon skipping versus alternative splice site usage in different tissues and highlight the liver, brain and testis as having particularly high levels of AS, supporting the idea that tissue-regulated AS plays important roles in the differentiation of these tissues. The high levels of alternative splice site usage in the liver may relate to the unusual patterns of splicing-factor expression observed in the adult liver, suggesting aspects of developmental regulation of AS at the tissue level. Obtaining a more comprehensive picture of AS will require the integration of additional types of data upstream and downstream of these core interactions. Upstream, splicing factors themselves may be differentially regulated in different tissues or in response to different stimuli at the level of transcription, splicing, or translation, and are frequently regulated by post-translational modifications such as phosphorylation, so systematic measurements of splicing factor levels and activities will be required. Downstream, AS may affect the stability of alternative transcripts (for example, in cases of messages subject to nonsense-mediated mRNA decay), and frequently alters functional properties of the encoded proteins, so systematic measurements of AS transcript and protein isoforms and functional assays will also be needed to fully understand the regulatory consequences of AS events. Ultimately, it will be important to place regulatory events involving AS into the context of regulatory networks involving control at the levels of transcription, translation and post-translational modifications.
Materials and methods
Data and resources
Chromosome assemblies of the human genome (hg13) were obtained from public databases . Transcript databases included approximately 94,000 human cDNA sequences obtained from GenBank (release 134.0, gbpri and gbhtc categories), and approximately 5 million human expressed sequence tags (ESTs) from dbEST (repository 02202003). Human ESTs were designated according to their cDNA library source (in total about 800) into different tissue types. Pertinent information about cDNA libraries and the corresponding human tissue or cell line was extracted from dbEST and subsequently integrated with library information retrieved from the Mammalian Gene Collection Initiative (MGC) , the Integrated Molecular Analysis of Gene Expression Consortium (IMAGE)  and the Cancer Genome Anatomy Project (CGAP) . Library information obtained from MGC, IMAGE and CGAP is provided in Additional data file 5.
Genome annotation by alignment of spliced transcripts
The GENOA genome annotation script  was used to align spliced cDNA and EST sequences to the human genome. GENOA uses BLASTN to detect significant blocks of identity between repeat-masked cDNA sequences and genomic DNA, and then aligns cDNAs to the genomic loci identified by BLASTN using the spliced-alignment algorithm MRNAVSGEN . This algorithm is similar in concept to SIM4  but was developed specifically to align high-quality cDNAs rather than ESTs and thus requires higher alignment quality (at least 93% identity) and consensus terminal dinucleotides at the ends of all introns (that is, GT..AG, GC..AG or AT..AC). EST sequences were aligned using SIM4 to those genomic regions that had aligned cDNAs. Stringent alignment criteria were imposed: ESTs were required to overlap cDNAs (so that all the genes studied were supported by at least one cDNA-genomic alignment); the first and last aligned segments of ESTs were required to be at least 30 nucleotides in length, with at least 90% sequence identity; and the entire EST sequence alignment was required to extend over at least 90% of the length of the EST with at least 90% sequence identity.
In total, GENOA aligned about 85,900 human cDNAs and about 890,300 ESTs to the human genome. The relatively low fraction of aligned ESTs (about 18%), and average aligned length of about 550 bases (the average lengths were not significantly different between different tissues, see Additional data file 6), reflect the stringent alignment-quality criteria that were imposed so as to be as confident as possible in the inferred splicing patterns. The aligned sequences yielded about 17,800 gene regions with more than one transcript aligned that exhibited a multi-exon structure. Of these, about 60% exhibited evidence of alternative splicing of internal exons. Our analysis did not examine differences in 3'-terminal and 5'-terminal exons, inclusion of which is frequently dictated by alternative polyadenylation or alternative transcription start sites and therefore does not represent 'pure' AS [73, 74]. The EST alignments were then used to categorize all internal exons as: constitutive exons; A3Es, A5Es, skipped exons, multiply alternatively spliced exons (for example, exons that exhibited both skipping and alternative 5' splice site usage), and exons that contained retained introns. An internal exon present in at least one transcript was identified as a skipped exon if it was precisely excluded in one or more other transcripts, such that the boundaries of both the 5' and 3' flanking exons were the same in the transcripts that included and skipped the exon (for example, exon E3 in Figure 1). Similarly, an internal exon present in at least one transcript was identified as an A3E or A5E if at least one other transcript contained an exon differing in length by the use of an alternative 3' splice site or 5' splice site. The 'core' of an A3E or A5E is defined as the exon portion that is common to all transcripts used to infer the AS event. The extension of an alternatively spliced exon is the exon portion added to the core region by the use of an alternative 3' splice site or 5' splice site) that is present in some, but not all transcripts used to infer the AS event. Pairs of inferred A3Es or A5Es differing by fewer than six nucleotides were excluded from further analysis, as in , because of the possibility that such small differences might sometimes result from EST sequencing or alignment errors. As the frequency of insertion-deletions errors greater than three bases using modern sequencing techniques is vanishingly small (P. Green, personal communication), a six-base cutoff should exclude the vast majority of such errors. Alternatively spliced exons/genes identified in specific tissues are available for download from the GENOA website .
Quantifying splice junction differences between alternative mRNA isoforms
To quantify the difference in splicing patterns between mRNAs or ESTs derived from a gene locus, the splice junction difference ratio (SJD) was calculated. For any pair of mRNAs/ESTs that have been aligned to overlapping portions of a genomic locus, the SJD is defined as the fraction of the splice junctions that occur in overlapping portions of the two transcripts that differ in one or both splice sites. A sample calculation is given in Figure 4. The SJD measure was calculated by taking the ratio of the number of 'valid' splice junctions that differ between two sequences over the total number of splice junctions, when comparing a pair of ESTs across all splice junctions present in overlapping portions of the two transcripts. A splice junction was considered valid if: the 5' splice site and the 3' splice site satisfied either the GT..AG or the GC..AG dinucleotide sequences at exon-intron junctions; and if the splice junction was observed at least twice in different transcripts.
Identification of candidate splicing regulatory motifs
Over-represented sequence motifs (k-mers) were identified by comparing the number of occurrences of k-mers (for k in the range of 4 to 6 bases) in a test set of alternative exons versus a control set. In this analysis, monomeric tandem repeats (for example, poly(A) sequences) were excluded. The enrichment score of candidate k-mers in the test set versus the control set was evaluated by computing χ2 (chi-squared) values with a Yates correction term , using an approach similar in spirit to that described by Brudno et al. . We randomly sampled 500 subsets of the same size as the test set from the control set. The enrichment scores for k-mers over-represented in the sampled subset versus the remainder of the control set were computed as above. The estimated p-value for observing the given enrichment score (χ2-value) associated with an over-represented sequence motif of length k was defined as the fraction of subsets that contained any k-mer with enrichment score (χ2-value) higher than the tested motif. Correcting for multiple testing is not required as the p-value was defined relative to the most enriched k-mer for each sampled set. For sets of skipped exons from human brain- and testis-derived EST sequences, the test sets comprised 1,265 and 517 exons skipped in brain- and testis-derived ESTs, respectively, and the control sets comprised 12,527 and 8,634 exons constitutively included in human brain- and testis-derived ESTs, respectively. Candidate sequence motifs in skipped exons from brain and testis-derived ESTs with associated p-values less than 0.002 were retained. For the set of A5E and A3E events from human liver-derived EST sequences, the test set comprised 44 A3Es and 45 A5Es, and the control set comprised 1,619 A3Es and 1,481 A5Es identified using ESTs from all tissues excluding liver. In this analysis, A3Es and A5Es with extension sequences of less than 25 bases were excluded and sequences longer than 150 bases were truncated to 150 bases, by retaining the exon sequence segment closest to the internal alternative splice junction. Over-represented sequence motifs in A3Es and A5Es from liver-derived EST sequences with associated p-values less than 0.01 were retained.
Gene-expression analysis of trans-acting splicing factors
SR proteins, SR-related proteins, and hnRNPs were derived from published proteomic analyses of the spliceosome [46–48]. Expression values for these genes were obtained from the 'gene expression atlas' using the HG-U95A DNA microarray  and from a similar set of expression data using the HG-U133A DNA microarray . Altogether, 20 splicing factors - ASF/SF2, SRm300, SC35, SRp40, SRp55, SRp30c, 9G8, SRp54, SFRS10, SRp20, hnRNPs A1, A2/B2, C, D, G, H1, K, L, M, and RALY - were studied in 26 different tissues present in both microarray experiments (Figure 5). The data from each gene chip - HG-U95A and HG-U133A - were analyzed separately. The average difference (AD) value of each probe was used as the indicator of expression level. In analyzing these microarray data, AD values smaller than 20 were standardized to 20, as in . When two or more probes mapped to a single gene, the values from those probes were averaged. Pearson (product-moment) correlation coefficients between 20-dimensional vectors for all tissue pairs were calculated, using data from each of the two DNA microarray studies separately.
Additional data files
Additional data files containing the following supplementary data, tables and figures are available with the online version of this paper and from the GENOA genome annotation website . The lists of GenBank accession numbers of human cDNAs and ESTs that were mapped to the human genome by the GENOA pipeline, GENOA gene locus identifiers, and gene loci with spliced alignments for the 22 human autosomes and two sex chromosomes are provided at our website . Sets of constitutive and alternative exons in genes expressed in the human brain, testis and liver, and control sets used are also provided . Additional data file 1 lists the average and median number of ESTs per gene and tissue, and the total number of genes per tissue using different minimum numbers of ESTs. Additional data file 2 lists the average total number of AS genes and AS genes containing SEs, A3Es and A5Es using ESTs derived from normal, non-diseased tissues. Additional data file 3 lists the number of constitutively spliced and AS genes, and AS genes containing SEs, A3Es and A5Es. Additional data file 4 shows the average fractions of AS genes and average fractions of AS genes containing SEs, A3Es and A5Es using ESTs derived from normal, non-disease-derived tissues. Additional data file 5 lists categories of cDNA libraries and designated tissues derived from the MGC, IMAGE and CGAP. Additional data file 6 shows the average lengths of ESTs that aligned to gene loci expressed in different tissues. Additional data file 7 lists human splicing factors of SR, SR-related and hnRNP genes, corresponding Ensembl gene numbers and Affymetrix microarray probe identification numbers. Additional data file 8 shows the distribution of the average Pearson correlation coefficient values across different tissues for expression levels of random sets of genes obtained from the Affymetrix microarray data.
We thank T. Poggio and P. Sharp for stimulating discussions and the anonymous reviewers for constructive suggestions. This work was supported by grants from the National Science Foundation and the National Institutes of Health (C.B.B.), and by a Lee Kuan Yew fellowship (G.Y.) and a Whiteman fellowship (G.K.).
- Black DL: Mechanisms of alternative pre-messenger RNA splicing. Annu Rev Biochem. 2003, 72: 291-336. 10.1146/annurev.biochem.72.121801.161720.PubMedView ArticleGoogle Scholar
- Cartegni L, Chew SL, Krainer AR: Listening to silence and understanding nonsense: exonic mutations that affect splicing. Nat Rev Genet. 2002, 3: 285-298. 10.1038/nrg775.PubMedView ArticleGoogle Scholar
- Graveley BR: Alternative splicing: increasing diversity in the proteomic world. Trends Genet. 2001, 17: 100-107. 10.1016/S0168-9525(00)02176-4.PubMedView ArticleGoogle Scholar
- Lopez AJ: Alternative splicing of pre-mRNA: developmental consequences and mechanisms of regulation. Annu Rev Genet. 1998, 32: 279-305. 10.1146/annurev.genet.32.1.279.PubMedView ArticleGoogle Scholar
- Grabowski PJ: Genetic evidence for a Nova regulator of alternative splicing in the brain. Neuron. 2000, 25: 254-256. 10.1016/S0896-6273(00)80888-0.PubMedView ArticleGoogle Scholar
- Zavolan M, Kondo S, Schonbach C, Adachi J, Hume DA, Hayashizaki Y, Gaasterland T, RIKEN GER Group, GSL members: Impact of alternative initiation, splicing, and termination on the diversity of the mRNA transcripts encoded by the mouse transcriptome. Genome Res. 2003, 13: 1290-1300. 10.1101/gr.1017303.PubMedPubMed CentralView ArticleGoogle Scholar
- Mironov AA, Fickett JW, Gelfand MS: Frequent alternative splicing of human genes. Genome Res. 1999, 9: 1288-1293. 10.1101/gr.9.12.1288.PubMedPubMed CentralView ArticleGoogle Scholar
- Modrek B, Resch A, Grasso C, Lee C: Genome-wide detection of alternative splicing in expressed sequences of human genes. Nucleic Acids Res. 2001, 29: 2850-2859. 10.1093/nar/29.13.2850.PubMedPubMed CentralView ArticleGoogle Scholar
- Modrek B, Lee C: A genomic view of alternative splicing. Nat Genet. 2002, 30: 13-19. 10.1038/ng0102-13.PubMedView ArticleGoogle Scholar
- Kan Z, States D, Gish W: Selecting for functional alternative splices in ESTs. Genome Res. 2002, 12: 1837-1845. 10.1101/gr.764102.PubMedPubMed CentralView ArticleGoogle Scholar
- Brett D, Hanke J, Lehmann G, Haase S, Delbruck S, Krueger S, Reich J, Bork P: EST comparison indicates 38% of human mRNAs contain possible alternative splice forms. FEBS Lett. 2000, 474: 83-86. 10.1016/S0014-5793(00)01581-7.PubMedView ArticleGoogle Scholar
- Faustino NA, Cooper TA: Pre-mRNA splicing and human disease. Genes Dev. 2003, 17: 419-437. 10.1101/gad.1048803.PubMedView ArticleGoogle Scholar
- Modafferi EF, Black DL: A complex intronic splicing enhancer from the c-src pre-mRNA activates inclusion of a heterologous exon. Mol Cell Biol. 1997, 17: 6537-6545.PubMedPubMed CentralView ArticleGoogle Scholar
- Naor D, Nedvetzki S, Golan I, Melnik L, Faitelson Y: CD44 in cancer. Crit Rev Clin Lab Sci. 2002, 39: 527-579.PubMedView ArticleGoogle Scholar
- MacDougall C, Harbison D, Bownes M: The developmental consequences of alternate splicing in sex determination and differentiation in Drosophila. Dev Biol. 1995, 172: 353-376. 10.1006/dbio.1995.8047.PubMedView ArticleGoogle Scholar
- Ladd AN, Charlet N, Cooper TA: The CELF family of RNA binding proteins is implicated in cell-specific and developmentally regulated alternative splicing. Mol Cell Biol. 2001, 21: 1285-1296. 10.1128/MCB.21.4.1285-1296.2001.PubMedPubMed CentralView ArticleGoogle Scholar
- Jensen KB, Dredge BK, Stefani G, Zhong R, Buckanovich RJ, Okano HJ, Yang YY, Darnell RB: Nova-1 regulates neuron-specific alternative splicing and is essential for neuronal viability. Neuron. 2000, 25: 359-371. 10.1016/S0896-6273(00)80900-9.PubMedView ArticleGoogle Scholar
- Lim LP, Burge CB: A computational analysis of sequence features involved in recognition of short introns. Proc Natl Acad Sci USA. 2001, 98: 11193-11198. 10.1073/pnas.201407298.PubMedPubMed CentralView ArticleGoogle Scholar
- Fairbrother WG, Yeh RF, Sharp PA, Burge CB: Predictive identification of exonic splicing enhancers in human genes. Science. 2002, 297: 1007-1013. 10.1126/science.1073774.PubMedView ArticleGoogle Scholar
- Liu HX, Zhang M, Krainer AR: Identification of functional exonic splicing enhancer motifs recognized by individual SR proteins. Genes Dev. 1998, 12: 1998-2012.PubMedPubMed CentralView ArticleGoogle Scholar
- Schaal TD, Maniatis T: Selection and characterization of pre-mRNA splicing enhancers: identification of novel SR protein-specific enhancer sequences. Mol Cell Biol. 1999, 19: 1705-1719.PubMedPubMed CentralView ArticleGoogle Scholar
- Tian H, Kole R: Strong RNA splicing enhancers identified by a modified method of cycled selection interact with SR protein. J Biol Chem. 2001, 276: 33833-33839. 10.1074/jbc.M102957200.PubMedView ArticleGoogle Scholar
- Zhu J, Shendure J, Mitra RD, Church GM: Single molecule profiling of alternative pre-mRNA splicing. Science. 2003, 301: 836-838. 10.1126/science.1085792.PubMedView ArticleGoogle Scholar
- Johnson JM, Castle J, Garrett-Engele P, Kan Z, Loerch PM, Armour CD, Santos R, Schadt EE, Stoughton R, Shoemaker DD: Genome-wide survey of human alternative pre-mRNA splicing with exon junction microarrays. Science. 2003, 302: 2141-2144. 10.1126/science.1090100.PubMedView ArticleGoogle Scholar
- Hu GK, Madore SJ, Moldover B, Jatkoe T, Balaban D, Thomas J, Wang Y: Predicting splice variant from DNA chip expression data. Genome Res. 2001, 11: 1237-1245. 10.1101/gr.165501.PubMedPubMed CentralView ArticleGoogle Scholar
- Clark TA, Sugnet CW, Ares M: Genome-wide analysis of mRNA processing in yeast using splicing-specific microarrays. Science. 2002, 296: 907-910. 10.1126/science.1069415.PubMedView ArticleGoogle Scholar
- Ule J, Jensen KB, Ruggiu M, Mele A, Ule A, Darnell RB: CLIP identifies Nova-regulated RNA networks in the brain. Science. 2003, 302: 1212-1215. 10.1126/science.1090095.PubMedView ArticleGoogle Scholar
- Clark F, Thanaraj TA: Categorization and characterization of transcript-confirmed constitutively and alternatively spliced introns and exons from human. Hum Mol Genet. 2002, 11: 451-464. 10.1093/hmg/11.4.451.PubMedView ArticleGoogle Scholar
- Stamm S, Zhu J, Nakai K, Stoilov P, Stoss O, Zhang MQ: An alternative-exon database and its statistical analysis. DNA Cell Biol. 2000, 19: 739-756. 10.1089/104454900750058107.PubMedView ArticleGoogle Scholar
- Xu Q, Modrek B, Lee C: Genome-wide detection of tissue-specific alternative splicing in the human transcriptome. Nucleic Acids Res. 2002, 30: 3754-3766. 10.1093/nar/gkf492.PubMedPubMed CentralView ArticleGoogle Scholar
- Brett D, Pospisil H, Valcarcel J, Reich J, Bork P: Alternative splicing and genome complexity. Nat Genet. 2002, 30: 29-30. 10.1038/ng803.PubMedView ArticleGoogle Scholar
- Lee CJ, Irizarry K: Alternative splicing in the nervous system: an emerging source of diversity and regulation. Biol Psychiatry. 2003, 54: 771-776. 10.1016/S0006-3223(03)00375-5.PubMedView ArticleGoogle Scholar
- Wang Z, Lo HS, Yang H, Gere S, Hu Y, Buetow KH, Lee MP: Computational analysis and experimental validation of tumor-associated alternative RNA splicing in human cancer. Cancer Res. 2003, 63: 655-657.PubMedGoogle Scholar
- Hui L, Zhang X, Wu X, Lin Z, Wang Q, Li Y, Hu G: Identification of alternatively spliced mRNA variants related to cancers by genome-wide ESTs alignment. Oncogene. 2004, 23: 3013-3023. 10.1038/sj.onc.1207362.PubMedView ArticleGoogle Scholar
- Kirkpatrick LL, McIlwain KA, Nelson DL: Alternative splicing in the murine and human FXR1 genes. Genomics. 1999, 59: 193-202. 10.1006/geno.1999.5868.PubMedView ArticleGoogle Scholar
- Kirkpatrick LL, McIlwain KA, Nelson DL: Comparative genomic sequence analysis of the FXR gene family: FMR1, FXR1, and FXR2. Genomics. 2001, 78: 169-177. 10.1006/geno.2001.6667.PubMedView ArticleGoogle Scholar
- Sugnet CW, Kent WJ, Ares M, Haussler D: Transcriptome and genome conservation of alternative splicing events in humans and mice. In Biocomputing 2004: Pac Symp Biocomput. Edited by: Altman RB, Dunker AK, Hunter L, Jung TA, Klein TE. 2004, Singapore: World Scientific, 66-77.Google Scholar
- Sorek R, Ast G: Intronic sequences flanking alternatively spliced exons are conserved between human and mouse. Genome Res. 2003, 13: 1631-1637. 10.1101/gr.1208803.PubMedPubMed CentralView ArticleGoogle Scholar
- Kaufmann D, Kenner O, Nurnberg P, Vogel W, Bartelt B: In NF1, CFTR, PER3, CARS and SYT7, alternatively included exons show higher conservation of surrounding intron sequences than constitutive exons. Eur J Hum Genet. 2004, 12: 139-149. 10.1038/sj.ejhg.5201098.PubMedView ArticleGoogle Scholar
- Megy K, Audic S, Claverie JM: Positional clustering of differentially expressed genes on human chromosomes 20, 21 and 22. Genome Biol. 2003, 4: P1-10.1186/gb-2003-4-2-p1.PubMedView ArticleGoogle Scholar
- Hillman RT, Green RE, Brenner SE: An unappreciated role for RNA surveillance. Genome Biol. 2004, 5: R8-10.1186/gb-2004-5-2-r8.PubMedPubMed CentralView ArticleGoogle Scholar
- Green RE, Lewis BP, Hillman RT, Blanchette M, Lareau LF, Garnett AT, Rio DC, Brenner SE: Widespread predicted nonsense-mediated mRNA decay of alternatively-spliced transcripts of human normal and disease genes. Bioinformatics. 2003, 19 (Suppl 1): I118-I121. 10.1093/bioinformatics/btg1015.PubMedView ArticleGoogle Scholar
- Su AI, Cooke MP, Ching KA, Hakak Y, Walker JR, Wiltshire T, Orth AP, Vega RG, Sapinoso LM, Moqrich A, et al: Large-scale analysis of the human and mouse transcriptomes. Proc Natl Acad Sci USA. 2002, 99: 4465-4470. 10.1073/pnas.012025199.PubMedPubMed CentralView ArticleGoogle Scholar
- Gene Expression Atlas. [http://expression.gnf.org]
- Su AI, Wiltshire T, Batalov S, Lapp H, Ching KA, Block D, Zhang J, Soden R, Hayakawa M, Kreiman G, et al: A gene atlas of the mouse and human protein-encoding transcriptomes. Proc Natl Acad Sci USA. 2004, 101: 6062-6067. 10.1073/pnas.0400782101.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhou Z, Licklider LJ, Gygi SP, Reed R: Comprehensive proteomic analysis of the human spliceosome. Nature. 2002, 419: 182-185. 10.1038/nature01031.PubMedView ArticleGoogle Scholar
- Jurica MS, Moore MJ: Pre-mRNA splicing: awash in a sea of proteins. Mol Cell. 2003, 12: 5-14. 10.1016/S1097-2765(03)00270-3.PubMedView ArticleGoogle Scholar
- Rappsilber J, Ryder U, Lamond AI, Mann M: Large-scale proteomic analysis of the human spliceosome. Genome Res. 2002, 12: 1231-1245. 10.1101/gr.473902.PubMedPubMed CentralView ArticleGoogle Scholar
- Hanamura A, Caceres JF, Mayeda A, Franza BR, Krainer AR: Regulated tissue-specific expression of antagonistic pre-mRNA splicing factors. RNA. 1998, 4: 430-444.PubMedPubMed CentralGoogle Scholar
- Bai Y, Lee D, Yu T, Chasin LA: Control of 3' splice site choice in vivo by ASF/SF2 and hnRNP A1. Nucleic Acids Res. 1999, 27: 1126-1134. 10.1093/nar/27.4.1126.PubMedPubMed CentralView ArticleGoogle Scholar
- Eperon IC, Makarova OV, Mayeda A, Munroe SH, Caceres JF, Hayward DG, Krainer AR: Selection of alternative 5' splice sites: role of U1 snRNP and models for the antagonistic effects of SF2/ASF and hnRNP A1. Mol Cell Biol. 2000, 20: 8303-8318. 10.1128/MCB.20.22.8303-8318.2000.PubMedPubMed CentralView ArticleGoogle Scholar
- Kamma H, Portman DS, Dreyfuss G: Cell type-specific expression of hnRNP proteins. Exp Cell Res. 1995, 221: 187-196. 10.1006/excr.1995.1366.PubMedView ArticleGoogle Scholar
- Caputi M, Mayeda A, Krainer AR, Zahler AM: hnRNP A/B proteins are required for inhibition of HIV-1 pre-mRNA splicing. EMBO J. 1999, 18: 4060-4067. 10.1093/emboj/18.14.4060.PubMedPubMed CentralView ArticleGoogle Scholar
- Caputi M, Zahler AM: SR proteins and hnRNP H regulate the splicing of the HIV-1 tev-specific exon 6D. EMBO J. 2002, 21: 845-855. 10.1093/emboj/21.4.845.PubMedPubMed CentralView ArticleGoogle Scholar
- Anderson L, Seilhamer J: A comparison of selected mRNA and protein abundances in human liver. Electrophoresis. 1997, 18: 533-537.PubMedView ArticleGoogle Scholar
- Futcher B, Latter GI, Monardo P, McLaughlin CS, Garrels JI: A sampling of the yeast proteome. Mol Cell Biol. 1999, 19: 7357-7368.PubMedPubMed CentralView ArticleGoogle Scholar
- Griffin TJ, Gygi SP, Ideker T, Rist B, Eng J, Hood L, Aebersold R: Complementary profiling of gene expression at the transcriptome and proteome levels in Saccharomyces cerevisiae. Mol Cell Proteomics. 2002, 1: 323-333. 10.1074/mcp.M200001-MCP200.PubMedView ArticleGoogle Scholar
- Kawamoto S, Matsumoto Y, Mizuno K, Okubo K, Matsubara K: Expression profiles of active genes in human and mouse livers. Gene. 1996, 174: 151-158. 10.1016/0378-1119(96)00512-4.PubMedView ArticleGoogle Scholar
- Brudno M, Gelfand MS, Spengler S, Zorn M, Dubchak I, Conboy JG: Computational analysis of candidate intron regulatory elements for tissue-specific alternative pre-mRNA splicing. Nucleic Acids Res. 2001, 29: 2338-2348. 10.1093/nar/29.11.2338.PubMedPubMed CentralView ArticleGoogle Scholar
- Chou MY, Underwood JG, Nikolic J, Luu MH, Black DL: Multisite RNA binding and release of polypyrimidine tract binding protein during the regulation of c-src neural-specific splicing. Mol Cell. 2000, 5: 949-957. 10.1016/S1097-2765(00)80260-9.PubMedView ArticleGoogle Scholar
- Chan RC, Black DL: The polypyrimidine tract binding protein binds upstream of neural cell-specific c-src exon N1 to repress the splicing of the intron downstream. Mol Cell Biol. 1997, 17: 4667-4676.PubMedPubMed CentralView ArticleGoogle Scholar
- Grabowski PJ: Splicing regulation in neurons: tinkering with cell-specific control. Cell. 1998, 92: 709-712. 10.1016/S0092-8674(00)81399-9.PubMedView ArticleGoogle Scholar
- Wollerton MC, Gooding C, Wagner EJ, Garcia-Blanco MA, Smith CW: Autoregulation of polypyrimidine tract binding protein by alternative splicing leading to nonsense-mediated decay. Mol Cell. 2004, 13: 91-100. 10.1016/S1097-2765(03)00502-1.PubMedView ArticleGoogle Scholar
- Burd CG, Dreyfuss G: RNA binding specificity of hnRNP A1: significance of hnRNP A1 high-affinity binding sites in pre-mRNA splicing. EMBO J. 1994, 13: 1197-1204.PubMedPubMed CentralGoogle Scholar
- Sirand-Pugnet P, Durosay P, Brody E, Marie J: An intronic (A/U)GGG repeat enhances the splicing of an alternative intron of the chicken beta-tropomyosin pre-mRNA. Nucleic Acids Res. 1995, 23: 3501-3507.PubMedPubMed CentralView ArticleGoogle Scholar
- McCullough AJ, Berget SM: G triplets located throughout a class of small vertebrate introns enforce intron borders and regulate splice site selection. Mol Cell Biol. 1997, 17: 4562-4571.PubMedPubMed CentralView ArticleGoogle Scholar
- UCSC Genome Browser. [http://genome.ucsc.edu]
- Mammalian Gene Collection (MGC) Initiative. [http://mgc.nci.nih.gov]
- The I.M.A.G.E. Consortium. [http://image.llnl.gov/image]
- The Cancer Gene Anatomy Project. [http://cgap.nci.nih.gov]
- GENOA file server. [http://genes.mit.edu/genoa]
- Florea L, Hartzell G, Zhang Z, Rubin GM, Miller W: A computer program for aligning a cDNA sequence with a genomic DNA sequence. Genome Res. 1998, 8: 967-974.PubMedPubMed CentralGoogle Scholar
- Berget SM: Exon recognition in vertebrate splicing. J Biol Chem. 1995, 270: 2411-2414.PubMedView ArticleGoogle Scholar
- Majewski J, Ott J: Distribution and characterization of regulatory elements in the human genome. Genome Res. 2002, 12: 1827-1836. 10.1101/gr.606402.PubMedPubMed CentralView ArticleGoogle Scholar
- Glantz SA: Primer of Biostatistics. 1997, New York: McGraw-Hill, 4Google Scholar
- Burge lab: Additional data directory 1. [http://genes.mit.edu/burgelab/Supplementary/yeo_holste04/Add_Datadir_1/]
- Burge lab: Additional data directory 2. [http://genes.mit.edu/burgelab/Supplementary/yeo_holste04/Add_Datadir_2/]
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.