Intron Gain Research Articles

Characterization of newly gained introns in Daphnia populations.

As one of the few known species in an active phase of intron proliferation, the microcrustacean Daphnia pulex is an especially attractive system for interrogating the gain and loss of introns in natural populations. In this study, we used a comparative population-genomic approach to identify and characterize 90 recently gained introns in this species. Molecular clock analyses indicate that these introns arose between 3.9 × 105 and 1.45 × 104 years ago, with a spike in intron proliferation approximately 5.2 × 104 to 1.22 × 105 years ago. Parallel gains at homologous positions contribute to 47.8% (43/90) of discovered new introns. A disproportionally large number of new introns were found in historically isolated populations in Oregon. Nonetheless, derived, intron-bearing alleles were also identified in a wide range of geographic locations, suggesting intron gain and, to a lesser degree, intron loss are important sources of genetic variation in natural populations of Daphnia. A majority (55/90 or 61.1%) of the identified neointrons have associated internal direct repeats with lengths and compositions that are unlikely to occur by chance, suggesting repeated bouts of staggered double-strand breaks (DSBs) during their evolution. Accordingly, internal, staggered DSBs may contribute to a passive trend toward increased length and sequence diversity in nascent introns.

Handling tRNA introns, archaeal way and eukaryotic way.

Introns are found in various tRNA genes in all the three kingdoms of life. Especially, archaeal and eukaryotic genomes are good sources of tRNA introns that are removed by proteinaceous splicing machinery. Most intron-containing tRNA genes both in archaea and eukaryotes possess an intron at a so-called canonical position, one nucleotide 3′ to their anticodon, while recent bioinformatics have revealed unusual types of tRNA introns and their derivatives especially in archaeal genomes. Gain and loss of tRNA introns during various stages of evolution are obvious both in archaea and eukaryotes from analyses of comparative genomics. The splicing of tRNA molecules has been studied extensively from biochemical and cell biological points of view, and such analyses of eukaryotic systems provided interesting findings in the past years. Here, I summarize recent progresses in the analyses of tRNA introns and the splicing process, and try to clarify new and old questions to be solved in the next stages.

Genome-wide analysis of the NADK gene family in plants.

BackgroundNAD(H) kinase (NADK) is the key enzyme that catalyzes de novo synthesis of NADP(H) from NAD(H) for NADP(H)-based metabolic pathways. In plants, NADKs form functional subfamilies. Studies of these families in Arabidopsis thaliana indicate that they have undergone considerable evolutionary selection; however, the detailed evolutionary history and functions of the various NADKs in plants are not clearly understood.Principal FindingsWe performed a comparative genomic analysis that identified 74 NADK gene homologs from 24 species representing the eight major plant lineages within the supergroup Plantae: glaucophytes, rhodophytes, chlorophytes, bryophytes, lycophytes, gymnosperms, monocots and eudicots. Phylogenetic and structural analysis classified these NADK genes into four well-conserved subfamilies with considerable variety in the domain organization and gene structure among subfamily members. In addition to the typical NAD_kinase domain, additional domains, such as adenylate kinase, dual-specificity phosphatase, and protein tyrosine phosphatase catalytic domains, were found in subfamily II. Interestingly, NADKs in subfamily III exhibited low sequence similarity (∼30%) in the kinase domain within the subfamily and with the other subfamilies. These observations suggest that gene fusion and exon shuffling may have occurred after gene duplication, leading to specific domain organization seen in subfamilies II and III, respectively. Further analysis of the exon/intron structures showed that single intron loss and gain had occurred, yielding the diversified gene structures, during the process of structural evolution of NADK family genes. Finally, both available global microarray data analysis and qRT-RCR experiments revealed that the NADK genes in Arabidopsis and Oryza sativa show different expression patterns in different developmental stages and under several different abiotic/biotic stresses and hormone treatments, underscoring the functional diversity and functional divergence of the NADK family in plants.ConclusionsThese findings will facilitate further studies of the NADK family and provide valuable information for functional validation of this family in plants.

Molecular phylogeny of C1 inhibitor depicts two immunoglobulin-like domains fusion in fishes and ray-finned fishes specific intron insertion after separation from zebrafish

C1 inhibitor (C1IN) is a multi-facet serine protease inhibitor in the plasma cascades, inhibiting several proteases, notably, regulates both complement and contact system activation. Despite huge advancements in the understanding of C1IN based on biochemical properties and its roles in the plasma cascades, the phylogenetic history of C1IN remains uncharacterized. To date, there is no comprehensive study illustrating the phylogenetic history of C1IN. Herein, we explored phylogenetic history of C1IN gene in vertebrates. Fishes have C1IN with two immunoglobulin like domains attached in the N-terminal region. The RCL regions of CIIN from fishes and tetrapod genomes have variations at the positions P2 and P1′. Gene structures of C1IN gene from selected ray-finned fishes varied in the Ig domain region with creation of novel intron splitting exon Im2 into Im2a and Im2b. This intron is limited to ray-finned fishes with genome size reduced below 1Gb. Hence, we suggest that genome compaction and associated double-strand break repairs are behind this intron gain. This study reveals the evolutionary history of C1IN and confirmed that this gene remains the same locus for ∼450 MY in 52 vertebrates analysed, but it is not found in frogs and lampreys.

Hitchhiking Selection Is Driving Intron Gain in a Pathogenic Fungus

The variability of intron density among eukaryotes is puzzling and still debated. Most previous studies have been limited because of the near absence of intron presence-absence polymorphism (IPAP) within species or because comparisons could be made only between distantly related species. We conducted population genetic analyses on eight loci showing IPAP to investigate the effect of natural selection on intron dynamics in a global collection of the panmictic fungal plant pathogen Zymoseptoria tritici and its very close relatives. Five of these loci likely represent recent intron gains because their absence is fixed among the closest relatives of Z. tritici, and three likely represent recent intron losses because their presence is fixed among the close relatives. We analyzed signatures of selection by comparing allele frequencies, nucleotide diversities, and rates of recombination and found compelling evidence that at least two out of the five intron-gain loci, a SWIM zinc-finger gene and a sugar transporter, are under directional selection favoring alleles that gained the intron. Our results suggest that the intron-present alleles of these loci are sweeping to fixation, providing a genetic hitchhiking mechanism to explain rapid intron gain in Z. tritici. The overall findings are consistent with the hypothesis that intron gains are more likely to be driven by selection while intron losses are more likely to be due to neutral processes such as genetic drift.

Evolution of Gene Structural Complexity: An Alternative-Splicing-Based Model Accounts for Intron-Containing Retrogenes

The structure of eukaryotic genes evolves extensively by intron loss or gain. Previous studies have revealed two models for gene structure evolution through the loss of introns: RNA-based gene conversion, dubbed the Fink model and retroposition model. However, retrogenes that experienced both intron loss and intron-retaining events have been ignored; evolutionary processes responsible for the variation in complex exon-intron structure were unknown. We detected hundreds of retroduplication-derived genes in human (Homo sapiens), fly (Drosophila melanogaster), rice (Oryza sativa), and Arabidopsis (Arabidopsis thaliana) and categorized them either as duplicated genes that have all introns lost or as duplicated genes that have at least lost one and retained one intron compared with the parental copy (intron-retaining [IR] type). Our new model attributes intron retention alternative splicing to the generation of these IR-type gene pairs. We presented 25 parental genes that have an intron retention isoform and have retained introns in the same locations in the IR-type duplicate genes, which directly support our hypothesis. Our alternative-splicing-based model in conjunction with the retroposition and Fink models can explain the IR-type gene observed. We discovered a greater percentage of IR-type genes in plants than in animals, which may be due to the abundance of intron retention cases in plants. Given the prevalence of intron retention in plants, this new model gives a support that plant genomes have very complex gene structures.

Analysis of CACTA transposases reveals intron loss as major factor influencing their exon/intron structure in monocotyledonous and eudicotyledonous hosts.

BackgroundCACTA elements are DNA transposons and are found in numerous organisms. Despite their low activity, several thousand copies can be identified in many genomes. CACTA elements transpose using a ‘cut-and-paste’ mechanism, which is facilitated by a DDE transposase. DDE transposases from CACTA elements contain, despite their conserved function, different exon numbers among various CACTA families. While earlier studies analyzed the ancestral history of the DDE transposases, no studies have examined exon loss and gain with a view of mechanisms that could drive the changes.ResultsWe analyzed 64 transposases from different CACTA families among monocotyledonous and eudicotyledonous host species. The annotation of the exon/intron boundaries showed a range from one to six exons. A robust multiple sequence alignment of the 64 transposases based on their protein sequences was created and used for phylogenetic analysis, which revealed eight different clades. We observed that the exon numbers in CACTA transposases are not specific for a host genome. We found that ancient CACTA lineages diverged before the divergence of monocotyledons and eudicotyledons. Most exon/intron boundaries were found in three distinct regions among all the transposases, grouping 63 conserved intron/exon boundaries.ConclusionsWe propose a model for the ancestral CACTA transposase gene, which consists of four exons, that predates the divergence of the monocotyledons and eudicotyledons. Based on this model, we propose pathways of intron loss or gain to explain the observed variation in exon numbers. While intron loss appears to have prevailed, a putative case of intron gain was nevertheless observed.

Distribution of Conventional and Nonconventional Introns in Tubulin (α and β) Genes of Euglenids

The nuclear genomes of euglenids contain three types of introns: conventional spliceosomal introns, nonconventional introns for which a splicing mechanism is unknown (variable noncanonical borders, RNA secondary structure bringing together intron ends), and so-called intermediate introns, which combine features of conventional and nonconventional introns. Analysis of two genes, tubA and tubB, from 20 species of euglenids reveals contrasting distribution patterns of conventional and nonconventional introns—positions of conventional introns are conserved, whereas those of the nonconventional ones are unique to individual species or small groups of closely related taxa. Moreover, in the group of phototrophic euglenids, 11 events of conventional intron loss versus 15 events of nonconventional intron gain were identified. A comparison of all nonconventional intron sequences highlighted the most conserved elements in their sequence and secondary structure. Our results led us to put forward two hypotheses. 1) The first one posits that mutational changes in intron sequence could lead to a change in their excision mechanism—intermediate introns would then be a transitional form between the conventional and nonconventional introns. 2) The second hypothesis concerns the origin of nonconventional introns—because of the presence of inverted repeats near their ends, insertion of MITE-like transposon elements is proposed as a possible source of new introns.

The Complex Intron Landscape and Massive Intron Invasion in a Picoeukaryote Provides Insights into Intron Evolution

Genes in pieces and spliceosomal introns are a landmark of eukaryotes, with intron invasion usually assumed to have happened early on in evolution. Here, we analyze the intron landscape of Micromonas, a unicellular green alga in the Mamiellophyceae lineage, demonstrating the coexistence of several classes of introns and the occurrence of recent massive intron invasion. This study focuses on two strains, CCMP1545 and RCC299, and their related individuals from ocean samplings, showing that they not only harbor different classes of introns depending on their location in the genome, as for other Mamiellophyceae, but also uniquely carry several classes of repeat introns. These introns, dubbed introner elements (IEs), are found at novel positions in genes and have conserved sequences, contrary to canonical introns. This IE invasion has a huge impact on the genome, doubling the number of introns in the CCMP1545 strain. We hypothesize that each IE class originated from a single ancestral IE that has been colonizing the genome after strain divergence by inserting copies of itself into genes by intron transposition, likely involving reverse splicing. Along with similar cases recently observed in other organisms, our observations in Micromonas strains shed a new light on the evolution of introns, suggesting that intron gain is more widespread than previously thought.

Phylogenetics and Gene Structure Dynamics of Polygalacturonase Genes in Aspergillus and Neurospora crassa.

Polygalacturonase (PG) gene is a typical gene family present in eukaryotes. Forty-nine PGs were mined from the genomes of Neurospora crassa and five Aspergillus species. The PGs were classified into 3 clades such as clade 1 for rhamno-PGs, clade 2 for exo-PGs and clade 3 for exo- and endo-PGs, which were further grouped into 13 sub-clades based on the polypeptide sequence similarity. In gene structure analysis, a total of 124 introns were present in 44 genes and five genes lacked introns to give an average of 2.5 introns per gene. Intron phase distribution was 64.5% for phase 0, 21.8% for phase 1, and 13.7% for phase 2, respectively. The introns varied in their sequences and their lengths ranged from 20 bp to 424 bp with an average of 65.9 bp, which is approximately half the size of introns in other fungal genes. There were 29 homologous intron blocks and 26 of those were sub-clade specific. Intron losses were counted in 18 introns in which no obvious phase preference for intron loss was observed. Eighteen introns were placed at novel positions, which is considerably higher than those of plant PGs. In an evolutionary sense both intron loss and gain must have taken place for shaping the current PGs in these fungi. Together with the small intron size, low conservation of homologous intron blocks and higher number of novel introns, PGs of fungal species seem to have recently undergone highly dynamic evolution.

Genome-wide search and comparative genomic analysis of the trypsin inhibitor-like cysteine-rich domain-containing peptides

It was shown that peptides containing trypsin inhibitor-like cysteine-rich (TIL) domain are able to inhibit proteinase activities, and thus play important roles in various biological processes, such as immune response and anticoagulation. However, only a limited number of the TIL peptides have been identified and characterized so far; and little has been known about the evolutionary relationships of the genes encoding the TIL peptides. BmKAPi is a TIL domain-containing peptide that was identified from Mesobuthus martensii Karsch. Here, we conducted genome-wide searches for new peptides that are homologous to BmKAPi or possess a cysteine pattern similar to that of BmKAPi. As a result, we identified a total of 80 different TIL peptides from 34 species of arthropods. We found that these peptides can be classified into seven evolutionarily distinct groups. Furthermore, we cloned the genomic sequence of BmKAPi; the genomic sequences of the majority of other TIL peptides were also identified from the GenBank database using bioinformatical approaches. Through phylogenetic and comparative genomic analysis, we found 26 cases of intron gain events occurred in the genes of the TIL peptides; however, no instances of intron loss were observed. Moreover, we found that alternative splicing contributes to the diversification of the TIL peptides. It is interesting to see that four genes of the TIL domain-containing peptides overlap in a DNA region located on the chromosome LG B15 of Bombus terretris. These data suggest that the evolution of the TIL peptide genes are dynamic, which was dominated by intron gain.

Gene invasion in distant eukaryotic lineages: discovery of mutually exclusive genetic elements reveals marine biodiversity

Inteins are rare, translated genetic parasites mainly found in bacteria and archaea, while spliceosomal introns are distinctly eukaryotic features abundant in most nuclear genomes. Using targeted metagenomics, we discovered an intein in an Atlantic population of the photosynthetic eukaryote, Bathycoccus, harbored by the essential spliceosomal protein PRP8 (processing factor 8 protein). Although previously thought exclusive to fungi, we also identified PRP8 inteins in parasitic (Capsaspora) and predatory (Salpingoeca) protists. Most new PRP8 inteins were at novel insertion sites that, surprisingly, were not in the most conserved regions of the gene. Evolutionarily, Dikarya fungal inteins at PRP8 insertion site a appeared more related to the Bathycoccus intein at a unique insertion site, than to other fungal and opisthokont inteins. Strikingly, independent analyses of Pacific and Atlantic samples revealed an intron at the same codon as the Bathycoccus PRP8 intein. The two elements are mutually exclusive and neither was found in cultured Bathycoccus or other picoprasinophyte genomes. Thus, wild Bathycoccus contain one of few non-fungal eukaryotic inteins known and a rare polymorphic intron. Our data indicate at least two Bathycoccus ecotypes exist, associated respectively with oceanic or mesotrophic environments. We hypothesize that intein propagation is facilitated by marine viruses; and, while intron gain is still poorly understood, presence of a spliceosomal intron where a locus lacks an intein raises the possibility of new, intein-primed mechanisms for intron gain. The discovery of nucleus-encoded inteins and associated sequence polymorphisms in uncultivated marine eukaryotes highlights their diversity and reveals potential sexual boundaries between populations indistinguishable by common marker genes.

Mechanisms of Intron Loss and Gain in the Fission Yeast Schizosaccharomyces

The fission yeast, Schizosaccharomyces pombe, is an important model species with a low intron density. Previous studies showed extensive intron losses during its evolution. To test the models of intron loss and gain in fission yeasts, we conducted a comparative genomic analysis in four Schizosaccharomyces species. Both intronization and de-intronization were observed, although both were at a low frequency. A de-intronization event was caused by a degenerative mutation in the branch site. Four cases of imprecise intron losses were identified, indicating that genomic deletion is not a negligible mechanism of intron loss. Most intron losses were precise deletions of introns, and were significantly biased to the 3′ sides of genes. Adjacent introns tended to be lost simultaneously. These observations indicated that the main force shaping the exon-intron structures of fission yeasts was precise intron losses mediated by reverse transcriptase. We found two cases of intron gains caused by tandem genomic duplication, but failed to identify the mechanisms for the majority of the intron gain events observed. In addition, we found that intron-lost and intron-gained genes had certain similar features, such as similar Gene Ontology categories and expression levels.

Genomic characterization of lytic Staphylococcus aureus phage GH15: providing new clues to intron shift in phages

Phage GH15 is a polyvalent phage that shows activity against a wide range of Staphylcoccus aureus strains. This study analysed the genome of GH15. The genome size of GH15 (139 806 bp) was found to be larger than that of the known staphylococcal phages, and the G+C content (30.23 mol%) of GH15 was lower than that of any other staphylococcal myovirus phages. By mass spectrometry, ten structural proteins were identified. Analysis revealed that GH15 was closely related to phages G1, ISP, A5W, Sb-1 and K, and was moderately related to Twort. In light of the variability in identity, coverage, G+C content and genome size, coupled with the large number of mosaicisms, there certainly were close evolutionary relationships from K to Sb-1, A5W, ISP, G1 and finally GH15. Interestingly, all the introns and inteins present in the above phages were absent in GH15 and there appeared to be intron loss in GH15 compared with the intron gain seen in other phages. A comparison of the intron- and intein-related genes demonstrated a clear distinction in the location of the insertion site between intron-containing and intron-free alleles, and this might lead to the establishment of a consensus sequence associated with the presence of an intron or intein. The comparative analysis of the GH15 genome sequence with other phages not only provides compelling evidence for the diversity of staphylococcal myovirus phages but also offers new clues to intron shift in phages.

Splicing-Related Features of Introns Serve to Propel Evolution

The role of spliceosomal intronic structures played in evolution has only begun to be elucidated. Comparative genomic analyses of fungal snoRNA sequences, which are often contained within introns and/or exons, revealed that about one-third of snoRNA-associated introns in three major snoRNA gene clusters manifested polymorphisms, likely resulting from intron loss and gain events during fungi evolution. Genomic deletions can clearly be observed as one mechanism underlying intron and exon loss, as well as generation of complex introns where several introns lie in juxtaposition without intercalating exons. Strikingly, by tracking conserved snoRNAs in introns, we found that some introns had moved from one position to another by excision from donor sites and insertion into target sties elsewhere in the genome without needing transposon structures. This study revealed the origin of many newly gained introns. Moreover, our analyses suggested that intron-containing sequences were more prone to sustainable structural changes than DNA sequences without introns due to intron's ability to jump within the genome via unknown mechanisms. We propose that splicing-related structural features of introns serve as an additional motor to propel evolution.

Genome-wide analysis of NBS-encoding disease resistance genes in Cucumis sativusand phylogenetic study of NBS-encoding genes in Cucurbitaceae crops

BackgroundPlant nucleotide-binding site (NBS)-leucine-rich repeat (LRR) proteins encoded by resistance genes play an important role in the responses of plants to various pathogens, including viruses, bacteria, fungi, and nematodes. In this study, a comprehensive analysis of NBS-encoding genes within the whole cucumber genome was performed, and the phylogenetic relationships of NBS-encoding resistance gene homologues (RGHs) belonging to six species in five genera of Cucurbitaceae crops were compared.ResultsCucumber has relatively few NBS-encoding genes. Nevertheless, cucumber maintains genes belonging to both Toll/interleukine-1 receptor (TIR) and CC (coiled-coil) families. Eight commonly conserved motifs have been established in these two families which support the grouping into TIR and CC families. Moreover, three additional conserved motifs, namely, CNBS-1, CNBS-2 and TNBS-1, have been identified in sequences from CC and TIR families. Analyses of exon/intron configurations revealed that some intron loss or gain events occurred during the structural evolution between the two families. Phylogenetic analyses revealed that gene duplication, sequence divergence, and gene loss were proposed as the major modes of evolution of NBS-encoding genes in Cucurbitaceae species. Compared with NBS-encoding sequences from the Arabidopsis thaliana genome, the remaining seven TIR familes of NBS proteins and RGHs from Cucurbitaceae species have been shown to be phylogenetically distinct from the TIR family of NBS-encoding genes in Arabidopsis, except for two subfamilies (TIR4 and TIR9). On the other hand, in the CC-NBS family, they grouped closely with the CC family of NBS-encoding genes in Arabidopsis. Thus, the NBS-encoding genes in Cucurbitaceae crops are shown to be ancient, and NBS-encoding gene expansions (especially the TIR family) may have occurred before the divergence of Cucurbitaceae and Arabidopsis.ConclusionThe results of this paper will provide a genomic framework for the further isolation of candidate disease resistance NBS-encoding genes in cucumber, and contribute to the understanding of the evolutionary mode of NBS-encoding genes in Cucurbitaceae crops.

Convergent intron gains in hymenopteran elongation factor-1α

The eukaryotic translation elongation factor-1α gene (eEF1A) has been used extensively in higher level phylogenetics of insects and other groups, despite being present in two or more copies in several taxa. Orthology assessment has relied heavily on the position of introns, but the basic assumption of low rates of intron loss and absence of convergent intron gains has not been tested thoroughly. Here, we study the evolution of eEF1A based on a broad sample of taxa in the insect order Hymenoptera. The gene is universally present in two copies – F1 and F2 – both of which apparently originated before the emergence of the order. An elevated ratio of non-synonymous versus synonymous substitutions and differences in rates of amino acid replacements between the copies suggest that they evolve independently, and phylogenetic methods clearly cluster the copies separately. The F2 copy appears to be ancient; it is orthologous with the copy known as F1 in Diptera, and is likely present in most insect orders. The hymenopteran F1 copy, which may or may not be unique to this order, apparently originated through retroposition and was originally intron free. During the evolution of the Hymenoptera, it has successively accumulated introns, at least three of which have appeared at the same position as introns in the F2 copy or in eEF1A copies in other insects. The sites of convergent intron gain are characterized by highly conserved nucleotides that strongly resemble specific intron-associated sequence motifs, so-called proto-splice sites. The significant rate of convergent intron gain renders intron–exon structure unreliable as an indicator of orthology in eEF1A, and probably also in other protein-coding genes.

Multiple intron gain and loss events occurred during the evolution of Cenp-A gene

Centromere protein A (CENP-A) is a histone H3 like protein, and it plays a very important role in chromosomal segregation during mitosis and meiosis. The analyses on the exon-intron organization of the Cenp-A gene in representative genomes revealed that multiple intron gain and loss events have occurred during the evolution of Cenp-A gene in opisthokonta (common ancestor of fungi and animals). Moreover, our results revealed that at least two positions were conserved in the intron gain and loss events during the evolution of the Cenp-A gene.

Complete plastid genomes from Ophioglossum californicum, Psilotum nudum, and Equisetum hyemale reveal an ancestral land plant genome structure and resolve the position of Equisetales among monilophytes.

BackgroundPlastid genome structure and content is remarkably conserved in land plants. This widespread conservation has facilitated taxon-rich phylogenetic analyses that have resolved organismal relationships among many land plant groups. However, the relationships among major fern lineages, especially the placement of Equisetales, remain enigmatic.ResultsIn order to understand the evolution of plastid genomes and to establish phylogenetic relationships among ferns, we sequenced the plastid genomes from three early diverging species: Equisetum hyemale (Equisetales), Ophioglossum californicum (Ophioglossales), and Psilotum nudum (Psilotales). A comparison of fern plastid genomes showed that some lineages have retained inverted repeat (IR) boundaries originating from the common ancestor of land plants, while other lineages have experienced multiple IR changes including expansions and inversions. Genome content has remained stable throughout ferns, except for a few lineage-specific losses of genes and introns. Notably, the losses of the rps16 gene and the rps12i346 intron are shared among Psilotales, Ophioglossales, and Equisetales, while the gain of a mitochondrial atp1 intron is shared between Marattiales and Polypodiopsida. These genomic structural changes support the placement of Equisetales as sister to Ophioglossales + Psilotales and Marattiales as sister to Polypodiopsida. This result is augmented by some molecular phylogenetic analyses that recover the same relationships, whereas others suggest a relationship between Equisetales and Polypodiopsida.ConclusionsAlthough molecular analyses were inconsistent with respect to the position of Marattiales and Equisetales, several genomic structural changes have for the first time provided a clear placement of these lineages within the ferns. These results further demonstrate the power of using rare genomic structural changes in cases where molecular data fail to provide strong phylogenetic resolution.

Modeling the evolution dynamics of exon-intron structure with a general random fragmentation process

BackgroundMost eukaryotic genes are interrupted by spliceosomal introns. The evolution of exon-intron structure remains mysterious despite rapid advance in genome sequencing technique. In this work, a novel approach is taken based on the assumptions that the evolution of exon-intron structure is a stochastic process, and that the characteristics of this process can be understood by examining its historical outcome, the present-day size distribution of internal translated exons (exon). Through the combination of simulation and modeling the size distribution of exons in different species, we propose a general random fragmentation process (GRFP) to characterize the evolution dynamics of exon-intron structure. This model accurately predicts the probability that an exon will be split by a new intron and the distribution of novel insertions along the length of the exon.ResultsAs the first observation from this model, we show that the chance for an exon to obtain an intron is proportional to its size to the 3rd power. We also show that such size dependence is nearly constant across gene, with the exception of the exons adjacent to the 5′ UTR. As the second conclusion from the model, we show that intron insertion loci follow a normal distribution with a mean of 0.5 (center of the exon) and a standard deviation of 0.11. Finally, we show that intron insertions within a gene are independent of each other for vertebrates, but are more negatively correlated for non-vertebrate. We use simulation to demonstrate that the negative correlation might result from significant intron loss during evolution, which could be explained by selection against multi-intron genes in these organisms.ConclusionsThe GRFP model suggests that intron gain is dynamic with a higher chance for longer exons; introns are inserted into exons randomly with the highest probability at the center of the exon. GRFP estimates that there are 78 introns in every 10 kb coding sequences for vertebrate genomes, agreeing with empirical observations. GRFP also estimates that there are significant intron losses in the evolution of non-vertebrate genomes, with extreme cases of around 57% intron loss in Drosophila melanogaster, 28% in Caenorhabditis elegans, and 24% in Oryza sativa.