Molecular evolution of calmodulin-like domain protein kinases (CDPKs) in plants and protists

Many genes for calmodulin-like domain protein kinases (CDPKs) have been identified in plants and Alveolate protists. To study the molecular evolution of the CDPK gene family, we performed a phylogenetic analysis of CDPK genomic sequences. Analysis of introns supports the phylogenetic analysis; CDPK genes with similar intron/exon structure are grouped together on the phylogenetic tree. Conserved introns support a monophyletic origin for plant CDPKs, CDPK-related kinases, and phosphoenolpyruvate carboxylase kinases. Plant CDPKs divide into two major branches. Plant CDPK genes on one branch share common intron positions with protist CDPK genes. The introns shared between protist and plant CDPKs presumably originated before the divergence of plants from Alveolates. Additionally, the calmodulin-like domains of protist CDPKs have intron positions in common with animal and fungal calmodulin genes. These results, together with the presence of a highly conserved phase zero intron located precisely at the beginning of the calmodulin-like domain, suggest that the ancestral CDPK gene could have originated from the fusion of protein kinase and calmodulin genes facilitated by recombination of ancient introns. Received: 11 July 2000 / Accepted: 18 April 2001     

Complex Evolutionary Patterns of tRNAUAALeu Group I Introns in Cyanobacterial Radiation

Based on the findings that plastids and cyanobacteria have similar group I introns inserted into tRNAUAALeu genes, these introns have been suggested to be immobile and of ancient origin. In contrast, recent evidence suggests lateral transfer of cyanobacterial group I introns located in tRNAUAALeu genes. In light of these new findings, we have readdressed the evolution and lateral transfer of tRNAUAALeu group I introns in cyanobacteral radiation. We determined the presence of introns in 38 different strains, representing the major cyanobacterial lineages, and characterized the introns in 22 of the strains. Notably, two of these strains have two tRNAUAALeu genes, with each of these genes interrupted by introns, while three of the strains have both interrupted and uninterrupted genes. Two evolutionary distinct clusters of tRNA genes, with the genes interrupted by introns belonging to two distinct intron clusters, were identified. We also compared 16S rDNA and intron evolution for both closely and distantly related strains. The distribution of the introns in the clustered groups, as defined from 16S rDNA analysis, indicates relatively recent gain and/or loss of the introns in some of these lineages. The comparative analysis also suggests differences in the phylogenetic trees for 16S rDNA and the tRNAUAALeu group I introns. Taken together, our results show that the evolution of the intron is considerably more complex than previous studies found to be the case. We discuss, based on our results, evolutionary models involving lateral intron transfer and models involving differential loss of the intron.     

Evolutionary convergence on highly-conserved 3' intron structures in intron-poor eukaryotes and insights into the ancestral eukaryotic genome

The presence of spliceosomal introns in eukaryotes raises a range of questions about genomic evolution. Along with the fundamental mysteries of introns' initial proliferation and persistence, the evolutionary forces acting on intron sequences remain largely mysterious. Intron number varies across species from a few introns per genome to several introns per gene, and the elements of intron sequences directly implicated in splicing vary from degenerate to strict consensus motifs. We report a 50-species comparative genomic study of intron sequences across most eukaryotic groups. We find two broad and striking patterns. First, we find that some highly intron-poor lineages have undergone evolutionary convergence to strong 3' consensus intron structures. This finding holds for both branch point sequence and distance between the branch point and the 3' splice site. Interestingly, this difference appears to exist within the genomes of green alga of the genus Ostreococcus, which exhibit highly constrained intron sequences through most of the intron-poor genome, but not in one much more intron-dense genomic region. Second, we find evidence that ancestral genomes contained highly variable branch point sequences, similar to more complex modern intron-rich eukaryotic lineages. In addition, ancestral structures are likely to have included polyT tails similar to those in metazoans and plants, which we found in a variety of protist lineages. Intriguingly, intron structure evolution appears to be quite different across lineages experiencing different types of genome reduction: whereas lineages with very few introns tend towards highly regular intronic sequences, lineages with very short introns tend towards highly degenerate sequences. Together, these results attest to the complex nature of ancestral eukaryotic splicing, the qualitatively different evolutionary forces acting on intron structures across modern lineages, and the impressive evolutionary malleability of eukaryotic gene structures.     

Exon junction sequences as cryptic splice sites: implications for intron origin

Introns are flanked by a partially conserved coding sequence that forms the immediate exon junction sequence following intron removal from pre-mRNA. Phylogenetic evidence indicates that these sequences have been targeted by numerous intron insertions during evolution, but little is known about this process. Here, we test the prediction that exon junction sequences were functional splice sites that existed in the coding sequence of genes prior to the insertion of introns. To do this, we experimentally identified nine cryptic splice sites within the coding sequence of actin genes from humans, Arabidopsis, and Physarum by inactivating their normal intron splice sites. We found that seven of these cryptic splice sites correspond exactly to the positions of exon junctions in actin genes from other species. Because actin genes are highly conserved, we could conclude that at least seven actin introns are flanked by cryptic splice sites, and from the phylogenetic evidence, we could also conclude that actin introns were inserted into these cryptic splice sites during evolution. Furthermore, our results indicate that these insertion events were dependent upon the splicing machinery. Because most introns are flanked by similar sequences, our results are likely to be of general relevance.     

Evolutionary dynamics of introns in plastid-derived genes in plants: saturation nearly reached but slow intron gain continues

Some of the principal transitions in the evolution of eukaryotes are characterized by engulfment of prokaryotes by primitive eukaryotic cells. In particular, approximately 1.6 billion years ago, engulfment of a cyanobacterium that became the ancestor of chloroplasts and other plastids gave rise to Plantae, the major branch of eukaryotes comprised of glaucophytes, red algae, green algae, and green plants. After endosymbiosis, there was large-scale migration of genes from the endosymbiont to the nuclear genome of the host such that approximately 18% of the nuclear genes in Arabidopsis appear to be of chloroplast origin. To gain insights into the process of evolution of gene structure in these, originally, intronless genes, we compared the properties and the evolutionary dynamics of introns in genes of plastid origin and ancestral eukaryotic genes in Arabidopsis, poplar, and rice genomes. We found that intron densities in plastid-derived genes were slightly but significantly lower than those in ancestral eukaryotic genes. Although most of the introns in both categories of genes were conserved between monocots (rice) and dicots (Arabidopsis and poplar), lineage-specific intron gain was more pronounced in plastid-derived genes than in ancestral genes, whereas there was no significant difference in the intron loss rates between the 2 classes of genes. Thus, after the transfer to the nuclear genome, the plastid-derived genes have undergone a massive intron invasion that, by the time of the divergence of dicots and monocots (150-200 MYA), yielded intron densities only slightly lower than those in ancestral genes. Nevertheless, the accumulation of introns in plastid-derived genes appears not to have reached saturation and continues to this time, albeit at a low rate. The overall pattern of intron gain and loss in the plastid-derived genes is shaped by this continuing gain and the more general tendency for loss that is characteristic of the recent evolution of plant genes.     

Prevalence of intron gain over intron loss in the evolution of paralogous gene families

The mechanisms and evolutionary dynamics of intron insertion and loss in eukaryotic genes remain poorly understood. Reconstruction of parsimonious scenarios of gene structure evolution in paralogous gene families in animals and plants revealed numerous gains and losses of introns. In all analyzed lineages, the number of acquired new introns was substantially greater than the number of lost ancestral introns. This trend held even for lineages in which vertical evolution of genes involved more intron losses than gains, suggesting that gene duplication boosts intron insertion. However, dating gene duplications and the associated intron gains and losses based on the molecular clock assumption showed that very few, if any, introns were gained during the last ~100 million years of animal and plant evolution, in agreement with previous conclusions reached through analysis of orthologous gene sets. These results are generally compatible with the emerging notion of intensive insertion and loss of introns during transitional epochs in contrast to the relative quiet of the intervening evolutionary spans.     

Tempo and Mode of Spliceosomal Intron Evolution in Actin of Foraminifera

Spliceosomal introns are present in almost all eukaryotic genes, yet little is known about their origin and turnover in the majority of eukaryotic phyla. There is no agreement whether most introns are ancestral and have been lost in some lineage or have been gained recently. We addressed this question by analyzing the spatial and temporal distribution of introns in actins of foraminifera, a group of testate protists whose exceptionally rich fossil record permits the calibration of molecular phylogenies to date intron origins. We identified 24 introns dispersed along the sequence of two foraminiferan actin paralogues and actin deviating proteins, an unconventional type of fast-evolving actin found in some foraminifera. Comparison of intron positions indicates that 20 of 24 introns are specific to foraminifera. Four introns shared between foraminifera and other eukaryotes were interpreted as parallel gains because they have been found only in single species belonging to phylogenetically distinctive lineages. Moreover, additional recent intron gain due to the transfer between the actin paralogues was observed in two cultured species. Based on a relaxed molecular clock timescale, we conclude that intron gains in actin took place throughout the evolution of foraminifera, with the oldest introns inserted between 550 and 500 million years ago and the youngest ones acquired less than 100 million years ago. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Debashish Bhattacharya]     

Nested Evolution of a tRNALeu(UAA) Group I Intron by both Horizontal Intron Transfer and Recombination of the Entire tRNA Locus

The origin and evolution of bacterial introns are still controversial issues. Here we present data on the distribution and evolution of a recently discovered divergent tRNA(Leu)(UAA) intron. The intron shows a higher sequence affiliation with introns in tRNA(Ile)(CAU) and tRNA(Arg)(CCU) genes in alpha- and beta-proteobacteria, respectively, than with other cyanobacterial tRNA(Leu)(UAA) group I introns. The divergent tRNA(Leu)(UAA) intron is sporadically distributed both within the Nostoc and the Microcystis radiations. The complete tRNA gene, including flanking regions and intron from Microcystis aeruginosa strain NIVA-CYA 57, was sequenced in order to elucidate the evolutionary pattern of this intron. Phylogenetic reconstruction gave statistical evidence for different phylogenies for the intron and exon sequences, supporting an evolutionary model involving horizontal intron transfer. The distribution of the tRNA gene, its flanking regions, and the introns were addressed by Southern hybridization and PCR amplification. The tRNA gene, including the flanking regions, were absent in the intronless stains but present in the intron-containing strains. This suggests that the sporadic distribution of this intron within the Microcystis genus cannot be attributed to intron mobility but rather to an instability of the entire tRNA(Leu)(UAA) intron-containing genome region. Taken together, the complete data set for the evolution of this intron can best be explained by a model involving a nested evolution of the intron, i.e., wherein the intron has been transferred horizontally (probably through a single or a few events) to a tRNA(Leu)(UAA) gene which is located within a unstable genome region.     

Conservation versus parallel gains in intron evolution

Orthologous genes from distant eukaryotic species, e.g. animals and plants, share up to 25–30% intron positions. However, the relative contributions of evolutionary conservation and parallel gain of new introns into this pattern remain unknown. Here, the extent of independent insertion of introns in the same sites (parallel gain) in orthologous genes from phylogenetically distant eukaryotes is assessed within the framework of the protosplice site model. It is shown that protosplice sites are no more conserved during evolution of eukaryotic gene sequences than random sites. Simulation of intron insertion into protosplice sites with the observed protosplice site frequencies and intron densities shows that parallel gain can account but for a small fraction (5–10%) of shared intron positions in distantly related species. Thus, the presence of numerous introns in the same positions in orthologous genes from distant eukaryotes, such as animals, fungi and plants, appears to reflect mostly bona fide evolutionary conservation.     

SPIn: model selection for phylogenetic mixtures via linear invariants

In phylogenetic inference, an evolutionary model describes the substitution processes along each edge of a phylogenetic tree. Misspecification of the model has important implications for the analysis of phylogenetic data. Conventionally, however, the selection of a suitable evolutionary model is based on heuristics or relies on the choice of an approximate input tree. We introduce a method for model Selection in Phylogenetics based on linear INvariants (SPIn), which uses recent insights on linear invariants to characterize a model of nucleotide evolution for phylogenetic mixtures on any number of components. Linear invariants are constraints among the joint probabilities of the bases in the operational taxonomic units that hold irrespective of the tree topologies appearing in the mixtures. SPIn therefore requires no input tree and is designed to deal with nonhomogeneous phylogenetic data consisting of multiple sequence alignments showing different patterns of evolution, for example, concatenated genes, exons, and/or introns. Here, we report on the results of the proposed method evaluated on multiple sequence alignments simulated under a variety of single-tree and mixture settings for both continuous- and discrete-time models. In the simulations, SPIn successfully recovers the underlying evolutionary model and is shown to perform better than existing approaches.     

Structure, chromosome location, and expression of the human smooth muscle (enteric type) gamma-actin gene: evolution of six human actin genes.

Recombinant phages that carry the human smooth muscle (enteric type) gamma-actin gene were isolated from human genomic DNA libraries. The amino acid sequence deduced from the nucleotide sequence matches those of cDNAs but differs from the protein sequence previously reported at one amino acid position, codon 359. The gene containing one 5' untranslated exon and eight coding exons extends for 27 kb on human chromosome 2. The intron between codons 84 and 85 (site 3) is unique to the two smooth muscle actin genes. In the 5' flanking region, there are several CArG boxes and E boxes, which are regulatory elements in some muscle-specific genes. Hybridization with the 3' untranslated region, which is specific for the human smooth muscle gamma-actin gene, suggests the single gene in the human genome and specific expressions in enteric and aortic tissues. From characterized molecular structures of the six human actin isoform genes, we propose a hypothesis of evolutionary pathway of the actin gene family. A presumed ancestral actin gene had introns at least sites 1, 2, and 4 through 8. Cytoplasmic actin genes may have directly evolved from it through loss of introns at sites 5 and 6. However, through duplication of the ancestral actin gene with substitutions of many amino acids, a prototype of muscle actin genes had been created. Subsequently, striated muscle actin and smooth muscle actin genes may have evolved from this prototype by loss of an intron at site 4 and acquisition of a new intron at site 3, respectively.     

Analysis of evolution of exon-intron structure of eukaryotic genes

The availability of multiple, complete eukaryotic genome sequences allows one to address many fundamental evolutionary questions on genome scale. One such important, long-standing problem is evolution of exon-intron structure of eukaryotic genes. Analysis of orthologous genes from completely sequenced genomes revealed numerous shared intron positions in orthologous genes from animals and plants and even between animals, plants and protists. The data on shared and lineage-specific intron positions were used as the starting point for evolutionary reconstruction with parsimony and maximum-likelihood approaches. Parsimony methods produce reconstructions with intron-rich ancestors but also infer lineage-specific, in many cases, high levels of intron loss and gain. Different probabilistic models gave opposite results, apparently depending on model parameters and assumptions, from domination of intron loss, with extremely intron-rich ancestors, to dramatic excess of gains, to the point of denying any true conservation of intron positions among deep eukaryotic lineages. Development of models with adequate, realistic parameters and assumptions seems to be crucial for obtaining more definitive estimates of intron gain and loss in different eukaryotic lineages. Many shared intron positions were detected in ancestral eukaryotic paralogues which evolved by duplication prior to the divergence of extant eukaryotic lineages. These findings indicate that numerous introns were present in eukaryotic genes already at the earliest stages of evolution of eukaryotes and are compatible with the hypothesis that the original, catastrophic intron invasion accompanied the emergence of the eukaryotic cells. Comparison of various features of old and younger introns starts shedding light on probable mechanisms of intron insertion, indicating that propagation of old introns is unlikely to be a major mechanism for origin of new ones. The existence and structure of ancestral protosplice sites were addressed by examining the context of introns inserted within codons that encode amino acids conserved in all eukaryotes and, accordingly, are not subject to selection for splicing efficiency. It was shown that introns indeed predominantly insert into or are fixed in specific protosplice sites which have the consensus sequence (A/C)AG|Gt.     

Selective constraints on intron evolution in Drosophila

Intron sizes show an asymmetrical distribution in a number of organisms, with a large number of "short" introns clustered around a minimal intron length and a much broader distribution of longer introns. In Drosophila melanogaster, the short intron class is centered around 61 bp. The narrow length distribution suggests that natural selection may play a role in maintaining intron size. A comparison of 15 orthologous introns among species of the D. melanogaster subgroup indicates that, in general, short introns are not under greater DNA sequence or length constraints than long introns. There is a bias toward deletions in all introns (deletion/insertion ratio is 1.66), and the vast majority of indels are of short length (<10 bp). Indels occurring on the internal branches of the phylogenetic tree are significantly longer than those occurring on the terminal branches. These results are consistent with a compensatory model of intron length evolution in which slightly deleterious short deletions are frequently fixed within species by genetic drift, and relatively rare larger insertions that restore intron length are fixed by positive selection. A comparison of paralogous introns shared among duplicated genes suggests that length constraints differ between introns within the same gene. The janusA, janusB, and ocnus genes share two short introns derived from a common ancestor. The first of these introns shows significantly fewer indels than the second intron, although the two introns show a comparable number of substitutions. This indicates that intron-specific selective constraints have been maintained following gene duplication, which preceded the divergence of the D. melanogaster species subgroup.     

Analysis of nonuniformity in intron phase distribution.

The distribution of different intron groups with respect to phases has been analyzed. It has been established that group II introns and nuclear introns have a minimum frequency of phase 2 introns. Since the phase of introns is an extremely conservative measure the observed minimum reflects evolutionary processes. A sample of all known, group I introns was too small to provide a valid characteristic of their phase distribution. The findings observed for the unequal distribution of phases cannot be explained solely on the basis of the mobile properties of introns. One of the most likely explanations for this nonuniformity in the intron phase distribution is the process of exon shuffling. It is proposed that group II introns originated at the early stages of evolution and were involved in the process of exon shuffling.     

Standardized phylogenetic tree: a reference to discover functional evolution

Functional evolution is often driven by positive natural selection. Although it is thought to be rare in evolution at the molecular level, its effects may be observed as the accelerated evolutionary rates. Therefore one of the effective ways to identify functional evolution is to identify accelerated evolution. Many methods have been developed to test the statistical significance of the accelerated evolutionary rate by comparison with the appropriate reference rate. The rates of synonymous substitution are one of the most useful and popular references, especially for large-scale analyses. On the other hand, these rates are applicable only to a limited evolutionary time period because they saturate quickly--i.e., multiple substitutions happen frequently because of the lower functional constraint. The relative rate test is an alternative method. This technique has an advantage in terms of the saturation effect but is not sufficiently powerful when the evolutionary rate differs considerably among phylogenetic lineages. For the aim to provide a universal reference tree, we propose a method to construct a standardized tree which serves as the reference for accelerated evolutionary rate. The method is based upon multiple molecular phylogenies of single genes with the aim of providing higher reliability. The tree has averaged and normalized branch lengths with standard deviations for statistical neutrality limits. The standard deviation also suggests the reliability level of the branch order. The resulting tree serves as a reference tree for the reliability level of the branch order and the test of evolutionary rate acceleration even when some of the species lineages show an accelerated evolutionary rate for most of their genes due to bottlenecking and other effects.     

Near intron positions are reliable phylogenetic markers: an application to holometabolous insects

Today, the reconstruction of the organismal evolutionary tree is based mainly on molecular sequence data. However, sequence data are sometimes insufficient to reliably resolve in particular deep branches. Thus, it is highly desirable to find novel, more reliable types of phylogenetic markers that can be derived from the wealth of genomic data. Here, we consider the gain of introns close to older preexisting ones. Because correct splicing is impeded by very small exons, nearby pairs of introns very rarely coexist, that is, the gain of the new intron is nearly always associated with the loss of the old intron. Both events may even be directly connected as in cases of intron migration. Therefore, it should be possible to identify one of the introns as ancient (plesiomorphic) and the other as novel (derived or apomorphic). To test the suitability of such near intron pairs (NIPs) as a marker class for phylogenetic analysis, we undertook an analysis of the evolutionary positions of bees and wasps (Hymenoptera) and beetles (Coleoptera) in relation to moths (Lepidoptera) and dipterans (Diptera) using recently completed genome project data. By scanning 758 putatively orthologous gene structures of Apis mellifera (Hymenoptera) and Tribolium castaneum (Coleoptera), we identified 189 pairs of introns, one from each species, which are located less than 50 nt from each other. A comparison with genes from 5 other holometabolan and 9 metazoan outgroup genomes resulted in 22 shared derived intron positions found in beetle as well as in butterflies and/or dipterans. This strongly supports a basal position of hymenopterans in the holometabolous insect tree. In addition, we found 31 and 12 intron positions apomorphic for A. mellifera and T. castaneum, respectively, which seem to represent changes inside these branches. Another 12 intron pairs indicate parallel intron gains or extraordinarily small exons. In conclusion, we show here that the analysis of phylogenetically nested, nearby intron pairs is suitable to identify evolutionarily younger intron positions and to determine their relative age, which should be of equal importance for the understanding of intron evolution and the reconstruction of the eukaryotic tree.     

Remarkable interkingdom conservation of intron positions and massive,lineage-specific intron loss and gain in eukaryotic evolution

Sequencing of eukaryotic genomes allows one to address major evolutionary problems, such as the evolution of gene structure. We compared the intron positions in 684 orthologous gene sets from 8 complete genomes of animals, plants, fungi, and protists and constructed parsimonious scenarios of evolution of the exon-intron structure for the respective genes. Approximately one-third of the introns in the malaria parasite Plasmodium falciparum are shared with at least one crown group eukaryote; this number indicates that these introns have been conserved through >1.5 billion years of evolution that separate Plasmodium from the crown group. Paradoxically, humans share many more introns with the plant Arabidopsis thaliana than with the fly or nematode. The inferred evolutionary scenario holds that the common ancestor of Plasmodium and the crown group and, especially, the common ancestor of animals, plants, and fungi had numerous introns. Most of these ancestral introns, which are retained in the genomes of vertebrates and plants, have been lost in fungi, nematodes, arthropods, and probably Plasmodium. In addition, numerous introns have been inserted into vertebrate and plant genes, whereas, in other lineages, intron gain was much less prominent.