Imprecise intron losses are less frequent than precise intron losses but are not rare in plants
Biology Direct volume 10, Article number: 24 (2015)
In this study, we identified 19 intron losses, including 11 precise intron losses (PILs), six imprecise intron losses (IILs), one de-exonization, and one exon deletion in tomato and potato, and 17 IILs in Arabidopsis thaliana. Comparative analysis of related genomes confirmed that all of the IILs have been fixed during evolution. Consistent with previous studies, our results indicate that PILs are a major type of intron loss. However, at least in plants, IILs are unlikely to be as rare as previously reported.
This article was reviewed by Jun Yu and Zhang Zhang. For complete reviews, see the Reviewers’ Reports section.
Theoretically, five different types of molecular events can inactivate introns or cause the deletion of an intron from a gene, thereby contributing to a decrease in intron abundance (Additional file 1). The first type of event is precise intron loss (PIL); in this case, intron losses do not affect the integrity of flanking exons. The second type of event is imprecise intron loss (IIL), which is accompanied by the insertion and/or deletion (indel) of nucleotides into/from flanking exons. The third type of event is termed de-intronization; in this case, sequences are not deleted from the genome, but rather an intronic sequence is converted into an exonic sequence by mutations that deactivate splicing signals. The fourth type of event is termed de-exonization, which is the conversion of an internal exon into an internal portion of an intron by mutations that deactivate splicing signals. This process leads to the fusion of an exon and its flanking two introns, which creates a larger intron and therefore decreases the intron number. Finally, the deletion of an internal exon also results in the fusion of two neighboring introns and consequently decreases the intron number by one. In this paper, we used the term “intron loss” in a broad sense to include all five of the above types of intron variations. Almost all previously observed intron losses are PILs; IILs and other types of intron losses appear only rarely [1-15]. There are three possibilities for the observed patterns. The first possibility is that they occur at quite different frequencies. For example, if intron losses are mediated by mRNA molecules, then all intron losses should be PILs . The second possibility is that intron losses that change coding sequences have essentially been eliminated by purifying selection. The third possibility is that there is a methodological bias toward the identification of PILs. It is possible for intron losses to introduce indels into coding sequences and therefore significantly reduce the similarities between flanking coding sequences and their orthologous regions. To be confident in identifying cases of intron loss, researchers generally discard poorly aligned regions [9-11].
IILs are less frequent than PILs but are not rare in Solanum or Arabidopsis thaliana
As the genomes of tomatoes (Solanum lycopersicum) and potatoes (Solanum tuberosum) diverge by less than 10% , we obtained reliable alignments for most of their orthologs. By surveying intron-exon structural changes in tomato and potato, we found 11 cases of PIL and six cases of IIL (Figure 1 and Additional file 2).
The species A. thaliana diverged from its relative Arabidopsis lyrata less than ten million years ago . In comparing the genomes of these species, we found 17 IILs in A. thaliana (Additional file 3). A close examination of 114 cases of intron loss from A. thaliana that were reported in a previous study  revealed that 104 of these cases, which occurred in 98 genes, were PILs, two were IILs, and eight lacked support due to insufficient numbers of RNA-Seq reads. The two cases of IILs were included in our 17-case dataset of IILs.
Transcriptome data showed that all of the variant genes in our study were still actively expressed (Additional files 2 and 3). A close examination of these genes did not reveal any premature stop codons that were introduced by intron loss mutations. The IIL genes are unlikely to be pseudogenized. Indels caused by IILs in coding sequences are more likely to be selected against. It is possible that the IILs that we observed were recent events that would soon be eliminated. This possibility could be excluded if the variations that we observed in one species were also found to exist in another. For this reason, we investigated whether tomato’s wild relative, Solanum pimpinellifolium, shares intron variations with tomato. In S. lycopersicum, there are five IILs and five PILs. We surveyed the orthologous genes in the genome of S. pimpinellifolium and found all ten of the variations. Therefore, all of the intron losses that we observed in S. lycopersicum have been fixed during evolution.
Taking advantage of the availability of genome sequences corresponding to multiple A. thaliana lines, we also tested whether the intron variations that were observed in A. thaliana have been fixed during evolution. By surveying 17 IIL genes in 180 lines of A. thaliana from Sweden , we found that all 17 cases of IILs had been fixed. Similarly, by surveying 104 PILs in 180 lines of A. thaliana, we found that 101 cases of PIL have been fixed during evolution and only three PILs remained polymorphic, in genes AT1G48420, AT3G23080, and AT4G00350. However, the allele frequencies of these three PILs were high: 97.7% for AT1G48420, 36.1% for AT3G23080 and 97.2% for AT4G00350.
IILs are expected to be under negative selection because of the indels that they cause in coding sequences; therefore, a majority of IILs might be eliminated. Nevertheless, the fixed IILs were still found to comprise an appreciable proportion of intron losses in Solanum and Arabidopsis.
The relative frequencies of IILs were underestimated when comparing distantly related genomes
The practice of filtering unreliable alignments may have led to underestimates of the frequencies of IILs in previous studies. If this were the case, a higher proportion of IILs than PILs would have been undetectable when we compared each of the two Solanum genomes with a distantly related species versus when we compared the two Solanum genomes to each other. To test this possibility, we surveyed for the presence of Solanum IL genes in the genome of a rice, Oryza sativa, that diverged from Solanum 163 million years ago . Among the 11 PIL genes and six IIL genes that we surveyed, we found O. sativa orthologs for all 11 PIL genes and only two IIL genes. In principle, when a very low identity is observed between two aligned sequences, it indicates that either a large number of mutations accumulated after divergence or that a low-quality alignment was produced. To maintain accuracy, these alignments are generally discarded. We calculated the identities of the coding sequence alignments of the sites flanking the intron losses. For each intron loss, 45-bp-long regions of coding sequence (not including gaps in the alignments) corresponding to positions on each side of the loss were included in the calculation. Using the first quartile of the identities of all aligned coding sequences that were generated between Solanum and O. sativa (0.53) as the threshold to filter unreliable alignments, one IIL and two PILs in Solanum were discarded. In summary, only one of the six IILs that occurred in Solanum could be detected when comparing its genome to rice. In contrast, nine of the 11 PILs could be detected using the same method. This difference is statistically significant (χ2 test, P = 0.03). Similar analyses have been carried out on the introns that have been lost from the genome of A. thaliana. Rice and A. thaliana also diverged 163 million years ago . Among 17 IIL genes and 98 PIL genes that were found in A. thaliana, we found four IILs and 47 PILs when comparing its genome against rice. Thus, a lower frequency of IILs than PILs was observed in A. thaliana when the reference genome was O. sativa (23.5% vs. 45.1%, respectively).
Similar to intron losses, the majority of previous studies on intron gains have been restricted to highly conserved orthologous genes. Among these genes, very few or no intron gains were found in humans, mice, and Arabidopsis thaliana [19,21-23]. This is in stark contrast to a study that specifically explored intron gains by evaluating segmental duplications, in which tens of intron gains had been revealed in each of these three species . In that study, intron gains that were accompanied by insertions and/or deletions of coding sequences were not excluded.
Identifying the true relative frequencies of PILs and IILs presents a dilemma: accurate IIL to PIL ratios can only be obtained when recently diverged genomes are compared. However, recent divergence means that only a limited amount of time has elapsed to enable the accumulation of intron loss variations. For these reasons, it would be helpful in the future to extend the current study to additional eukaryotic lineages.
De-exonization and exon deletion in Solanum
We identified one case of de-exonization, one case of exon deletion, two cases of intronization, and one case of exonization in tomato and potato (Additional file 2). In the tomato gene Solyc09g016940.2, an internal exon and the 5′ splicing signal of its downstream intron were both lost. In the potato gene PGSC0003DMG400004043, an internal exon has been converted into an internal region of a larger intron. In addition to GT-AG boundaries, there are many cis-acting sequence elements and trans-acting factors that facilitate intron recognition and splicing ; therefore, it is reasonable to hypothesize that some intron variations do not involve changes to these boundaries.
Ma et al. reports an initial look into precisely how intron loss has happened within a particular plant species, where two genome sequences – one domesticated and another wild – are available, and found 19 intron losses, which are supported by transcription evidence. They also took an addition look on the Arabidopsis genome, inspired by their finding from the Solanum species. Different from intron gain, intron loss should be rather rare event as purifying selection always prevents its loss-of-function effect, such as what may happen to IILs. In addition, the form of intron losses in a context of gene structure is of curiosity also, where functional consequences are complex for different forms, such as PILs vs. IILs; the latter may have more severe loss-of-function effect than the former that would not change protein coding sequences in theory. The results from Ma et al. are consistent to this speculation. A bit of concern is Figure 4, where results from diverse species were plotted into a trend that is not supported by adequate evidence across enough data from multiple species.
Authors’ response: This sentence has been deleted.
Change “IILs: not much less frequent than PILs in Solanum” to “IIL is (or IILs are) not much…”
Authors’ response: This problem has been corrected.
Replace “focused” with “enriched” in “Similar to intron losses, most previous explorations of intron gains were also focused on highly conserved orthologous genes”
Authors’ response: This has been corrected.
Quality of written English: Needs some language corrections before being published
Authors’ response: The language of this manuscript has been edited by a professional language-editing service.
Reviewer 2: Zhang Zhang, Beijing Institute of Genomics, Chinese Academy of Sciences
The manuscript by Ma et al. presented comprehensive investigations on intron loss by comparing multiple plant genome sequences, including tomato, potato, Arabidopsis and rice. Based on considerate filtration and exclusions of questionable data, the manuscript concluded that precise intron losses are the major type of intron loss and imprecise intron losses are not so rare as previously reported.
Although I am not an expert in this field, the manuscript is well-written and provides solid results. However, one of my major concerns is why plant species are used for studying intron loss and how about human, monkey, chimp, etc., if also used. As mentioned, it is due to low divergence (e.g., <10% between tomato and potato), but I feel it might be better to provide more background on a variety of species. Accordingly, the related concern is the title “Imprecise intron losses are less frequent than precise intron losses, but not so rare”. If it is unexplored or does not hold true in other species, it would be safe to add “in plants” in the title.
Authors’ response: We added “in plants” to the title of the revised manuscript and will monitor imprecise intron losses of other lineages.
The y-axis title in Figure 4 should be consistent with the word used in the main text. Also, how’s the correlation between PIL% and divergence time in different species. I think it is positively correlated. If available, plot together in Figure 4 and estimate the loss rate.
Coulombe-Huntington J, Majewski J. Intron loss and gain in Drosophila. Mol Biol Evol. 2007;24(12):2842–50.
Farlow A, Meduri E, Dolezal M, Hua L, Schlotterer C. Nonsense-mediated decay enables intron gain in Drosophila. PLoS Genet. 2010;6(1):e1000819.
Roy SW, Penny D. Large-scale intron conservation and order-of-magnitude variation in intron loss/gain rates in apicomplexan evolution. Genome Res. 2006;16(10):1270–5.
Roy SW, Penny D. Patterns of intron loss and gain in plants: Intron loss-dominated evolution and genome-wide comparison of O. sativa and A. thaliana. Mol Biol Evol. 2007;24(1):171–81.
Roy SW, Penny D. Widespread intron loss suggests retrotransposon activity in ancient apicomplexans. Mol Biol Evol. 2007;24(9):1926–33.
Roy SW, Hartl DL. Very little intron loss/gain in Plasmodium: Intron loss/gain mutation rates and intron number. Genome Res. 2006;16(6):750–6.
Da Lage JL, Binder M, Hua-Van A, Janecek S, Casane D. Gene make-up: rapid and massive intron gains after horizontal transfer of a bacterial alpha-amylase gene to Basidiomycetes. BMC Evol Biol. 2013;13:40.
Mitrovich QM, Tuch BB, De La Vega FM, Guthrie C, Johnson AD. Evolution of yeast noncoding RNAs reveals an alternative mechanism for widespread intron loss. Science. 2010;330(6005):838–41.
Irimia M, Rukov JL, Penny D, Vinther J, Garcia-Fernandez J, Roy SW. Origin of introns by ‘intronization’ of exonic sequences. Trends Genet. 2008;24(8):378–81.
Yenerall P, Krupa B, Zhou L. Mechanisms of intron gain and loss in Drosophila. BMC Evol Biol. 2011;11(1):364.
Fink GR. Pseudogenes in yeast? Cell. 1987;49(1):5–6.
The Tomato Genome Consortium. The tomato genome sequence provides insights into fleshy fruit evolution. Nature. 2012;485(7400):635–41.
Hedges SB, Dudley J, Kumar S. TimeTree: a public knowledge-base of divergence times among organisms. Bioinformatics. 2006;22(23):2971–2.
Yang YF, Zhu T, Niu DK. Association of intron loss with high mutation rate in Arabidopsis: implications for genome size evolution. Genome Biol Evol. 2013;5(4):723–33.
Long Q, Rabanal FA, Meng D, Huber CD, Farlow A, Platzer A, et al. Massive genomic variation and strong selection in Arabidopsis thaliana lines from Sweden. Nat Genet. 2013;45(8):884–90.
Roy SW, Fedorov A, Gilbert W. Large-scale comparison of intron positions in mammalian genes shows intron loss but no gain. Proc Natl Acad Sci U S A. 2003;100(12):7158–62.
Coulombe-Huntington J, Majewski J. Characterization of intron loss events in mammals. Genome Res. 2007;17(1):23–32.
Fawcett JA, Rouzé P, van de Peer Y. Higher intron loss rate in Arabidopsis thaliana than A. lyrata is consistent with stronger selection for a smaller genome. Mol Biol Evol. 2012;29(2):849–59.
Gao X, Lynch M. Ubiquitous internal gene duplication and intron creation in eukaryotes. Proc Natl Acad Sci U S A. 2009;49:20818–23.
Kornblihtt AR, Schor IE, Allo M, Blencowe BJ. When chromatin meets splicing. Nat Struct Mol Biol. 2009;16(9):902–3.
Ming-Yue Ma and Tao Zhu contributed equally to this work.
Illustration of different types of intron loss.
Details of intron variations in potato and tomato. A complete list of supplementary materials, methods, and figures is provided in this file.
Details of intron variations in Arabidopsis thaliana. A complete list of supplementary materials, methods, and figures is given in this file.
Cite this article
Ma, MY., Zhu, T., Li, XN. et al. Imprecise intron losses are less frequent than precise intron losses but are not rare in plants. Biol Direct 10, 24 (2015). https://doi.org/10.1186/s13062-015-0056-7