Genes with a large intronic burden show greater evolutionary conservation on the protein level
© Gorlova et al.; licensee BioMed Central Ltd. 2014
Received: 30 July 2013
Accepted: 11 March 2014
Published: 16 March 2014
The existence of introns in eukaryotic genes is believed to provide an evolutionary advantage by increasing protein diversity through exon shuffling and alternative splicing. However, this eukaryotic feature is associated with the necessity of exclusion of intronic sequences, which requires considerable energy expenditure and can lead to splicing errors. The relationship between intronic burden and evolution is poorly understood. The goal of this study was to analyze the relationship between the intronic burden and the level of evolutionary conservation of the gene.
We found a positive correlation between the level of evolutionary conservation of a gene and its intronic burden. The level of evolutionary conservation was estimated using the conservation index (CI). The CI value was determined on the basis of the most distant ortholog of the human protein sequence and ranged from 0 (the gene was unique to the human genome) to 9 (an ortholog of the human gene was detected in plants). In multivariable model, both the number of introns and total intron size remained significant predictors of CI. We also found that the number of alternative splice variants was positively correlated with CI.
The expression level of a gene was negatively correlated with the number of introns and total size of intronic region. Genes with a greater intronic burden had lower density of missense and nonsense mutations in the coding regions of the gene, which suggests that they are under a stronger pressure from purifying selection.
We identified a positive association between intronic burden and CI. One of the possible explanations of this is the idea of a cost-benefits balance. Evolutionarily conserved (functionally important) genes can “afford” the negative consequences of maintaining multiple introns because these consequences are outweighed by the benefit of maintaining the gene. Evolutionarily conserved and functionally important genes may use introns to create novel splice variants to tune the gene function to developmental stage and tissue type.
KeywordsExon/intron structure Intronic burden Evolutionary conservation
The division of genes into introns and exons is a hallmark of eukaryotic evolution. This division is believed to be evolutionarily beneficial because it allows the production of multiple proteins from the same gene through alternative splicing and may accelerate the creation of novel proteins through exon shuffling [1–4].
However, little is known about the forces that influence the exon/intron structure of genes [5–9]. Several biologically important characteristics correlated with intronic burden have been identified. For example, highly expressed genes tend to have shorter introns . Similarly, a study of 391 genes from 19 eukaryotic species conducted by Carmel et al.  demonstrated that the probability of intronic gains is positively correlated with the level of evolutionary conservation of the gene. However, a whole-genome assessment of the association between the number and length of introns and the level of evolutionary conservation has not yet been conducted. We sought to determine whether evolutionary conservation was correlated with intronic load in a whole-genome analysis.
Results and discussion
Variation in number of introns among human genes
Association between the conservation index (CI), number of introns, and total intron size
We also observed a significant positive correlation between the number of introns in a gene and its CI (Spearman rank R = 0.16, N = 16,249, P < 10−6). Figure 2(c) and (d) show the relationship between the number of introns in a gene and the gene’s CI. On the basis of segmented regression analysis, the curve of this relationship was divided into two segments: the first segment of genes with 0 to 10 introns showed a linear positive correlation, and the second segment of genes with ≥ 11 introns showed a plateau in terms of correlation.
We used a multiple linear regression model to estimate the independent effects of the total intron length and the number of introns as predictors of CI. We analyzed only the genes that satisfied two conditions: 10 or fewer introns and total intron size between 3 and 30 kb (too few genes had a total intron size less than 3 kb). For this subset of the data the assumption of linearity is justified. We also checked the subset for independence of the errors, constant variance of the errors and normality of the distribution of the errors. Because all the conditions were met, multiple regression was applied to the data to assess whether the number of introns and the total intronic length are independent predictors of CI. In the multiple regression model including number of introns and total intron size both regression coefficients were significant: bnoi = 0.11, N = 1,420, P < 0.00003; btil = 0.21, N = 1,420, P < 10−6. These results suggest that both the number of introns and the total intron size are independent predictors of CI for a subset of the data we have used.
Intronic burden and gene expression
The cost-benefit hypothesis predicts that genes with a large intronic burden are more functionally important than genes with a small intronic burden. A gene’s level of evolutionary conservation reflects its functional significance, and conserved genes are more likely to be involved in basic biological functions [19, 20]. However, evolutionary conservation of a gene reflects not only its functional significance but also its evolutionary history, which is usually unknown. To further assess the functional significance of the genes, we used data on the density of functional polymorphisms. Genes with low functional importance accumulate functional (protein-changing) polymorphisms more easily [21, 22]. Using dbSNP data, we found that the number of synonymous single nucleotide polymorphisms (SNPs) per codon did not correlate with CI (R = 0.01, N = 16,194, P = 0.18), whereas the number of non-synonymous SNPs and stop-gained SNPs was negatively correlated with CI (R = −0.1, n = 16,194, P < 10−6 and R = −0.06, n = 16,194, P < 10−6, respectively). These findings are consistent with the idea that the density of potentially functional polymorphisms can be used as a measure of the functional importance of a gene. We used the ratio of missense to synonymous substitutions (MIS/SYN) and the ratio of nonsense to synonymous substitutions (NON/SYN) to assess the strength of purifying selection on a gene. Strength of purifying selection can be used as a measure of functional importance of the gene under the assumption that more important genes are less tolerant to missense and nonsense substitutions. We noted a significant negative correlation between MIS/SYN and the number of introns (R = −0.05, n = 3,363, P = 0.006). The correlation between the number of introns and NON/SYN was also negative and significant (R = −0.26, n = 3,363, P < 10−6). The correlation between MIS/SYN and the total intron length was negative but did not reach statistical significance (R = −0.03, n = 3,363, P = 0.11), whereas the correlation between NON/SYN and the total intron length was significant (R = −0.13, n = 3,363, P < 10−6). These results support the idea that genes with a large total intron size and/or multiple introns tend to be more functionally important than the genes with smaller intronic loads as manifested in lower density of missense and nonsense substitutions.
We assume that functionally important genes may have more and larger introns because they can better “afford” a larger intronic burden than can less important genes. The other possible explanation for the maintenance of introns in functionally important genes may be related to alternative splicing. Up to 90% of human genes undergo alternative splicing . Even if the majority of rare splice variants are the products of splicing errors , many of the alternative splice variants are functional [24–26]. Furthermore, exon boundaries often correspond to functional domains [27, 28]. Therefore, it is conceivable that a larger number of introns enables functionally important genes to use alternative splicing to adjust and modify their functions on the basis of developmental stage or tissue type. We found a significant positive correlation between the number of splice variants and the CI of genes (R = 0.07, N = 15.819, P < 10−6). The association remained significant after controlling for the number of introns or the the total intron length (P = 0.0001) and suggests that functionally important genes are more likely to undergo alternative splicing.
Conversely, genes with a small number of introns (0–2) and a small total intron size (<3 kb) may be different from other genes. We found that the smallest genes in the human genome have the lowest CIs and low expression levels. These genes also have a higher density of non-synonymous and stop-gained SNPs than other genes (t-test values of 2.1 [P = 0.03] and 2.02 [P = 0.02] respectively), but they do not have a higher density of synonymous SNPs (t-test 0.2, df = 16,192, P = 0.87). This observation is consistent with the results of the analysis conducted by Krylov at al. . In their analysis of 7 complete genomes, the authors demonstrated that the propensity of gene loss was positively correlated with the rate of accumulation of nucleotide substitutions and negatively correlated with the gene’s expression level. The smallest genes may be enriched with young genes that have yet to develop an important biological function and thus cannot accumulate or maintain multiple introns. An alternative possibility is that the smallest genes are mostly dying genes that for some reason (e.g., environmental changes) became less functional, lost their introns, and started to accumulate non-synonymous and stop-gained mutations.
Heterogeneity of genetic composition across intronic groups
Database for Annotation, Visualisation, and Integrated Discovery (DAVID) http://david.abcc.ncifcrf.gov/home.jsp was used to check if the top and bottom 5% of the genes in terms of intronic burden were enriched by gene categories. We found that the bottom 5% of the genes with lowest number of introns are enriched by G-protein coupled receptor genes (GPCRs) (Benjamini P = 6.8E-54). This enrichment, however, is unlikely to drive low CI in the group. We found no significant difference in CI between GPCRs and other genes: average CI was 4.01 ± 0.11 for GPCRs and 4.19 ± 0.09 for other genes, Mann–Whitney U Test, adjusted Z = 0.34, P = 0.73. We also did not notice a significant difference in expression levels between GPRCs and other genes (nonparametric Mann Whitney test = 0.56, P = 0.24). The top 5% of genes with largest number of introns are enriched by ATP-binding genes (P value for enrichment is 1.1E-40) and cytoskeleton-associated genes (P value for enrichment is 2.2E-19).
We believe that heterogeneity of genetic composition is unlikely to be a major driving force behind the association between the intronic burden and CI. Firstly, the curve describing the relationship between intronic burden and CI is monotonic (Figure 2) suggesting that the intronic number rather than effect of the gene composition drives the association. For the middle part of the distribution gene enrichment is relatively low but we still see strong positive association between intronic burden and CI. Also for several gene families showing enrichment we found no evidence that the enrichment significantly contributes to the association.
Our analysis demonstrated that evolutionarily conserved genes have a greater intronic burden. Previous research found a positive association between the level of evolutionary conservation and the size of the intronic region of a gene for a fraction of the human genes. We confirmed that association for the whole human genome. The results of our analysis also suggest that not only the total intronic length but also the number of introns is an independent predictor of the level of evolutionary conservation of a gene.
Using the latest and most precise estimates of gene expression levels, we demonstrated a negative association between the intronic burden and the level of expression. Similar to the above analysis, this investigation showed that both the total intronic length and the number of introns independently predict expression level. We found that the genes with the lowest intronic burden are different from most other human genes, suggesting that they could be evolutionarily young (which is why they tend to have a lower conservation index) and have yet to acquire an indispensable biological function.
In conclusion, the problem of evolutionary advantages (or disadvantages) of introns is complex. Having multiple introns is obviously associated with a burden in terms of energy and resources: a cell needs first to transcribe all introns, then to remove them. It has been shown that splicing is often associated with errors that produce abnormal product and may have a negative effect on cell survival . On the other hand, having multiple introns provides an opportunity for alternative splicing which can be associated with distinct and important biological functions . It is not clear if by alternative splicing the cell tries to make use of something already available, or alternative splicing is a driver of intronic additions. It is also possible that introns are gained more often in important genes because they have higher expression levels and more open chromatin structure. Additionally it is possible that the disadvantage of an intron gain is less in important compared to not-so-important genes. We demonstrated that both the total size of intronic sequences and the total number of introns are independent predictors of the level of evolutionary conservation. Our findings raise an interesting possibility that intronic burden could be used as a predictor of the functional importance of a gene.
Conservation index (CI) scale
An example of the most distant species with detectable human ortholog
Divergence time (million years)
Unique to humans
Data on the exon/intron structure of genes, including the number of alternatively spliced isoforms, were obtained from the Exon-Intron Database. This database was created and is curated by one of the authors (A.F.) , is based on NCBI Gene Bank data, and allows large-scale computational examination of exon/intron structure. We assessed gene expression levels using data from the RNA-Seq Atlas database . This database provides estimates of gene expression levels across 10 human tissues (colon, heart, hypothalamus, kidney, liver, lung, ovary, skeletal muscle, spleen, and testes) and is based on RNA sequencing, which provides a less biased estimate of gene expression than probe-based technologies. As previously published, Z-scores were used as quantitative measures of gene expression . In brief, Z-score is a measure of the expression of an individual gene to the gene’s relative to the expression distribution in a reference population. The reference distribution is based on the distribution of the expression of all genes in all tissues. The expression value of a given gene in given tissue (Z-score) represents the number of standard deviations away from the mean of expression in the reference population.
To assess the correlation between the density of functional polymorphisms (missense and nonsense mutations) and intronic burden, we used the dbSNP database. The density of functional polymorphisms in a given gene was computed by dividing the total number of reported missense and nonsense mutations by the size of the coding region. Only validated SNPs were used for the analysis. We used the number of SNPs per codon as a measure of the SNP density. The nonparametric Spearman correlation coefficient was used to assess the correlation between the size and number of introns and the CI. To identify segments of linearity we used R2-based method . Statistical analysis was done using STATA software (version 10, StataCorp LP, College Station, TX).
This research was supported in part by the David H. Koch Center for Applied Research of Genitourinary Cancers and the NIH grants 5 R01 CA149462 to OYG and 5 P30 CA016672 (MD Anderson’s Cancer Center Support Grant).
- Lynch M, Conery JS: The origins of genome complexity. Science. 2003, 302 (5649): 1401-1404. 10.1126/science.1089370.PubMedView ArticleGoogle Scholar
- Koonin EV, Csuros M, Rogozin IB: Whence genes in pieces: reconstruction of the exon-intron gene structures of the last eukaryotic common ancestor and other ancestral eukaryotes. Wiley Interdiscip Rev RNA. 2013, 4 (1): 93-105.PubMedView ArticleGoogle Scholar
- Rogozin IB, Carmel L, Csuros M, Koonin EV: Origin and evolution of spliceosomal introns. Biol Direct. 2012, 7: 11-10.1186/1745-6150-7-11.PubMedPubMed CentralView ArticleGoogle Scholar
- Fedorova L, Fedorov A: Introns in gene evolution. Genetica. 2003, 118 (2–3): 123-131.PubMedView ArticleGoogle Scholar
- Roy SW, Gilbert W: The evolution of spliceosomal introns: patterns, puzzles and progress. Nat Rev Genet. 2006, 7 (3): 211-221.PubMedGoogle Scholar
- Koralewski TE, Krutovsky KV: Evolution of exon-intron structure and alternative splicing. PLoS One. 2011, 6 (3): e18055-10.1371/journal.pone.0018055.PubMedPubMed CentralView ArticleGoogle Scholar
- Rogozin IB, Sverdlov AV, Babenko VN, Koonin EV: Analysis of evolution of exon-intron structure of eukaryotic genes. Brief Bioinform. 2005, 6 (2): 118-134. 10.1093/bib/6.2.118.PubMedView ArticleGoogle Scholar
- Roy M, Kim N, Xing Y, Lee C: The effect of intron length on exon creation ratios during the evolution of mammalian genomes. Rna. 2008, 14 (11): 2261-2273. 10.1261/rna.1024908.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang C, Li WH, Krainer AR, Zhang MQ: RNA landscape of evolution for optimal exon and intron discrimination. Proc Natl Acad Sci USA. 2008, 105 (15): 5797-5802. 10.1073/pnas.0801692105.PubMedPubMed CentralView ArticleGoogle Scholar
- Castillo-Davis CI, Mekhedov SL, Hartl DL, Koonin EV, Kondrashov FA: Selection for short introns in highly expressed genes. Nat Genet. 2002, 31 (4): 415-418.PubMedGoogle Scholar
- Carmel L, Rogozin IB, Wolf YI, Koonin EV: Evolutionarily conserved genes preferentially accumulate introns. Genome Res. 2007, 17 (7): 1045-1050. 10.1101/gr.5978207.PubMedPubMed CentralView ArticleGoogle Scholar
- Louhichi A, Fourati A, Rebai A: IGD: a resource for intronless genes in the human genome. Gene. 2011, 488 (1–2): 35-40.PubMedView ArticleGoogle Scholar
- Wagner AK, Soumerai SB, Zhang F, Ross-Degnan D: Segmented regression analysis of interrupted time series studies in medication use research. J Clin Pharm Ther. 2002, 27 (4): 299-309. 10.1046/j.1365-2710.2002.00430.x.PubMedView ArticleGoogle Scholar
- Krupp M, Marquardt JU, Sahin U, Galle PR, Castle J, Teufel A: RNA-Seq Atlas–a reference database for gene expression profiling in normal tissue by next-generation sequencing. Bioinformatics. 2012, 28 (8): 1184-1185. 10.1093/bioinformatics/bts084.PubMedView ArticleGoogle Scholar
- Lehninger AL, Nelson DL, Cox MM: Principles of Biochemistry. 1982, New York: WorthGoogle Scholar
- Hsu SN, Hertel KJ: Spliceosomes walk the line: splicing errors and their impact on cellular function. RNA Biol. 2009, 6 (5): 526-530. 10.4161/rna.6.5.9860.PubMedPubMed CentralView ArticleGoogle Scholar
- Hertel KJ: Combinatorial control of exon recognition. J Biol Chem. 2008, 283 (3): 1211-1215. 10.1074/jbc.R700035200.PubMedView ArticleGoogle Scholar
- Kandul NP, Noor MA: Large introns in relation to alternative splicing and gene evolution: a case study of Drosophila bruno-3. BMC Genet. 2009, 10: 67-PubMedPubMed CentralView ArticleGoogle Scholar
- Schena M: The evolutionary conservation of eukaryotic gene transcription. Experientia. 1989, 45 (10): 972-983. 10.1007/BF01953055.PubMedView ArticleGoogle Scholar
- Yuan J: Evolutionary conservation of a genetic pathway of programmed cell death. J Cell Biochem. 1996, 60 (1): 4-11. 10.1002/(SICI)1097-4644(19960101)60:1<4::AID-JCB2>3.0.CO;2-1.PubMedView ArticleGoogle Scholar
- Akashi H, Osada N, Ohta T: Weak selection and protein evolution. Genetics. 2012, 192 (1): 15-31. 10.1534/genetics.112.140178.PubMedPubMed CentralView ArticleGoogle Scholar
- Sleator RD: An overview of the processes shaping protein evolution. Sci Prog. 2010, 93 (Pt 1): 1-6.PubMedView ArticleGoogle Scholar
- Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, Mayr C, Kingsmore SF, Schroth GP, Burge CB: Alternative isoform regulation in human tissue transcriptomes. Nature. 2008, 456 (7221): 470-476. 10.1038/nature07509.PubMedPubMed CentralView ArticleGoogle Scholar
- Irimia M, Blencowe BJ: Alternative splicing: decoding an expansive regulatory layer. Curr Opin Cell Biol. 2012, 24 (3): 323-332. 10.1016/j.ceb.2012.03.005.PubMedView ArticleGoogle Scholar
- Kelemen O, Convertini P, Zhang Z, Wen Y, Shen M, Falaleeva M, Stamm S: Function of alternative splicing. Gene. 2012, 514 (1): 1-30.PubMedView ArticleGoogle Scholar
- Sanchez SE, Petrillo E, Kornblihtt AR, Yanovsky MJ: Alternative splicing at the right time. RNA Biol. 2011, 8 (6): 954-959. 10.4161/rna.8.6.17336.PubMedPubMed CentralView ArticleGoogle Scholar
- Cho YJ, Hitchcock-DeGregori SE: Relationship between alternatively spliced exons and functional domains in tropomyosin. Proc Natl Acad Sci USA. 1991, 88 (22): 10153-10157. 10.1073/pnas.88.22.10153.PubMedPubMed CentralView ArticleGoogle Scholar
- Schmidt EE, Davies CJ: The origins of polypeptide domains. Bioessays. 2007, 29 (3): 262-270. 10.1002/bies.20546.PubMedPubMed CentralView ArticleGoogle Scholar
- Krylov DM, Wolf YI, Rogozin IB, Koonin EV: Gene loss, protein sequence divergence, gene dispensability, expression level, and interactivity are correlated in eukaryotic evolution. Genome Res. 2003, 13 (10): 2229-2235. 10.1101/gr.1589103.PubMedPubMed CentralView ArticleGoogle Scholar
- Sayers EW, Barrett T, Benson DA, Bolton E, Bryant SH, Canese K, Chetvernin V, Church DM, Dicuccio M, Federhen S, Feolo M, Fingerman IM, Geer LY, Helmberg W, Kapustin Y, Krasnov S, Landsman D, Lipman DJ, Lu Z, Madden TL, Madej T, Maglott DR, Marchler-Bauer A, Miller V, Karsch-Mizrachi I, Ostell J, Panchenko A, Phan L, Pruitt KD, Schuler GD, et al: Database resources of the national center for biotechnology information. Nucleic Acids Res. 2012, 40 (Database issue): D13-D25.PubMedPubMed CentralView ArticleGoogle Scholar
- Kumar S, Filipski A, Swarna V, Walker A, Hedges SB: Placing confidence limits on the molecular age of the human-chimpanzee divergence. Proc Natl Acad Sci USA. 2005, 102 (52): 18842-18847. 10.1073/pnas.0509585102.PubMedPubMed CentralView ArticleGoogle Scholar
- Domazet-Loso T, Tautz D: An ancient evolutionary origin of genes associated with human genetic diseases. Mol Biol Evol. 2008, 25 (12): 2699-2707. 10.1093/molbev/msn214.PubMedPubMed CentralView ArticleGoogle Scholar
- Gibney G, Baxevanis AD: Searching NCBI databases using entrez. Curr Protoc Hum Genet. 2011, Chapter 6: Unit6 10-PubMedGoogle Scholar
- Kumar S, Hedges SB: A molecular timescale for vertebrate evolution. Nature. 1998, 392 (6679): 917-920. 10.1038/31927.PubMedView ArticleGoogle Scholar
- Chernikova D, Motamedi S, Csuros M, Koonin EV, Rogozin IB: A late origin of the extant eukaryotic diversity: divergence time estimates using rare genomic changes. Biol Direct. 2011, 6: 26-10.1186/1745-6150-6-26.PubMedPubMed CentralView ArticleGoogle Scholar
- Shepelev V, Fedorov A: Advances in the Exon-Intron Database (EID). Brief Bioinform. 2006, 7 (2): 178-185. 10.1093/bib/bbl003.PubMedView ArticleGoogle Scholar
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 credited.