Module-intron correlation and intron sliding in family F/10 xylanase genes

Xylanases are classified into two families, numbered F/10 and G/11 according to the similarity of amino acid sequences of their catalytic domain (Henrissat, B., Bairoch, A., 1993. New families in the classification of glycosyl hydrolases based on amino acid sequence similarities. Biochem. J. 293, 781-788). Three-dimensional structure of the catalytic domain of the family F/10 xylanase was reported (White, A., Withers, S.G., Gilkes, N.R., Rose, D.R., 1994. Crystal structure of the catalytic domain of the beta-1,4-glycanase Cex from Cellulomonas fimi. Biochemistry 33, 12546-12552). The domain was decomposed into 22 modules by centripetal profiles (Go, M., Nosaka, M., 1987. Protein architecture and the origin of introns. Cold Spring Harbor Symp. Quant. Biol. 52, 915-924; Noguti, T., Sakakibara, H., Go, M., 1993. Localization of hydrogen-bonds within modules in barnase. Proteins 16, 357-363). A module is a contiguous polypeptide segment of amino acid residues having a compact conformation within a globular domain. Collected 31 intron sites of the family F/10 xylanase genes from fungus were found to be correlated to module boundaries with considerable statistical force (p values <0.001). The relationship between the intron locations and protein structures provides supporting evidence for the ancient origin of introns, because such a relationship cannot be expected by random insertion of introns into eukaryotic genes, but it rather suggests pre-existence of introns in the ancestral genes of prokaryotes and eukaryotes. A phylogenetic tree of the fungal and bacterial xylanase sequences made two clusters; one includes both the bacterial and fungal genes, but the other consists of only fungal genes. The mixed cluster of bacterial genes without introns and the fungal genes with introns further supports the ancient origin of introns. Comparison of the conserved base sequences of introns indicates that sliding of a splice site occurred in Aspergillus kawachii gene by one base from the ancestral position. Substrate-binding sites of xylanase are localized on eight modules, and introns are found at both termini of six out of these functional modules. This result suggests that introns might play a functional role in shuffling the exons encoding the substrate-binding modules.     

Phylogenetic and Exon–Intron Structure Analysis of Fungal Subtilisins: Support for a Mixed Model of Intron Evolution

Phylogenetic and exon–intron structure analyses of intra- and interspecific fungal subtilisins in this study provided support for a mixed model of intron evolution: a synthetic theory of introns-early and introns-late speculations. Intraspecifically, there were three phase zero introns in Pr1A and its introns 1 and 2 located at the highly conserved positions were phylogentically congruent with coding region, which is in favor of the view of introns-early speculation, while intron 3 had two different sizes and was evolutionarily incongruent with coding region, the evidence for introns-late speculation. Noticeably, the subtilisin Pr1J gene from different strains of M. ansiopliae contained different number of introns, the strong evidence in support of introns-late theory. Interspecifically, phylogenetic analysis of 60 retrievable fungal subtilisins provided a clear relationship between amino acid sequence and gene exon–intron structure that the homogeneous sequences usually have a similar exon–infron structure. There were 10 intron positions inserted by highly biased phase zero introns across examined fungal subtilisin genes, half of these positions were highly conserved, while the others were species-specific, appearing to be of recent origins due to intron insertion, in favor of the introns-late theory. High conservations of positions 1 and 2 inserted by the high percentage of phase zero introns as well as the evidence of phylogenetic congruence between the evolutionary histories of intron sequences and coding region suggested that the introns at these two positions were primordial.Reviewing Editor:Dr. Manyuan Long     

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     

Do introns favor or avoid regions of amino acid conservation?

Are intron positions correlated with regions of high amino acid conservation? For a set of ancient conserved proteins, with intronless prokaryotic but intron-containing eukaryotic homologs, multiple sequence alignments identified residues invariant throughout evolution. Intron positions between codons show no preferences. However, introns lying after the first base of a codon prefer conserved regions, markedly in glycines. Because glycines are in excess in conserved regions, this behavior could reflect phase-one introns entering glycine residues randomly in the ancestral sequences. Examination of intron positions within codons of evolutionarily invariable amino acids showed that roughly 50% of these introns are bordered by guanines at both 5'- and 3'-ends, 25% have a G only before the intron, and 5% have a G only after the intron, whereas about 20% are bordered by nonguanine bases.     

Structural context of exons in protein domains: implications for protein modelling and design

Intron boundaries were extracted from genomic data and mapped onto single-domain human and murine protein structures taken from the Protein Data Bank. A first analysis of this set of proteins shows that intron boundaries prefer to be in non-regular secondary structure elements, while avoiding alpha-helices and beta-strands. This fact alone suggests an evolutionary model in which introns are constrained by protein structure, particularly by tertiary structure contacts. In addition, in silico recombination experiments of a subset of these proteins together with their homologues, including those in different species, show that introns have a tendency to occur away from artificial crossover hot spots. Altogether, these findings support a model in which genes can preferentially harbour introns in less constrained regions of the protein fold they code for. In the light of these findings, we discuss some implications for protein modelling and design.     

Molecular evolution of eukaryotic genomes: hemiascomycetous yeast spliceosomal introns

As part of the exploratory sequencing program Génolevures, visual scrutinisation and bioinformatic tools were used to detect spliceosomal introns in seven hemiascomycetous yeast species. A total of 153 putative novel introns were identified. Introns are rare in yeast nuclear genes (<5% have an intron), mainly located at the 5′ end of ORFs, and not highly conserved in sequence. They all share a clear non-random vocabulary: conserved splice sites and conserved nucleotide contexts around splice sites. Homologues of metazoan snRNAs and putative homologues of SR splicing factors were identified, confirming that the spliceosomal machinery is highly conserved in eukaryotes. Several introns’ features were tested as possible markers for phylogenetic analysis. We found that intron sizes vary widely within each genome, and according to the phylogenetic position of the yeast species. The evolutionary origin of spliceosomal introns was examined by analysing the degree of conservation of intron positions in homologous yeast genes. Most introns appeared to exist in the last common ancestor of present day yeast species, and then to have been differentially lost during speciation. However, in some cases, it is difficult to exclude a possible sliding event affecting a pre-existing intron or a gain of a novel intron. Taken together, our results indicate that the origin of spliceosomal introns is complex within a given genome, and that present day introns may have resulted from a dynamic flux between intron conservation, intron loss and intron gain during the evolution of hemiascomycetous yeasts.     

Support for the Coelomata clade of animals from a rigorous analysis of the pattern of intron conservation

Many intron positions are conserved in varying subsets of eukaryotic genomes and, consequently, comprise a potentially informative class of phylogenetic characters. Roy and Gilbert developed a method of phylogenetic reconstruction using the patterns of intron presence-absence in eukaryotic genes and, applying this method to the analysis of animal phylogeny, obtained support for an Ecdysozoa clade (Roy SW, Gilbert W. 2005. Resolution of a deep animal divergence by the pattern of intron conservation. Proc Natl Acad Sci USA. 102:4403-4408). The critical assumption in the method was the independence of intron loss in different branches of the phylogenetic tree. Here, this assumption is refuted by showing that the branch-specific intron loss rates are strongly correlated. We show that different tree topologies are obtained, in each case with a significant statistical support, when different subsets of intron positions are analyzed. The analysis of the conserved intron positions supports the Coelomata topology, that is, a clade comprised of arthropods and chordates, whereas the analysis of more variable intron positions favors the Ecdysozoa topology, that is, a clade of arthropods and nematodes. We show, however, that the support for Ecdysozoa is fully explained by parallel loss of introns in nematodes and arthropods, a factor that does not contribute to the analysis of the conserved introns. The developed procedure for the identification and analysis of conserved introns and other characters with minimal or no homoplasy is expected to be useful for resolving many hard phylogenetic problems.     

Phylogenetic mapping of intron positions: a case study of translation initiation factor eIF2gamma

Eukaryotic translation initiation factor 2 (eIF2) is a G protein that delivers the methionyl initiator tRNA to the small ribosomal subunit and releases it upon GTP hydrolysis after the recognition of the initiation codon. eIF2 is composed of three subunits, alpha, beta, and gamma. Subunit gamma shows the strongest conservation, and it confers both tRNA and GTP/GDP binding. Using intron positioning and protein sequence alignment, here we show that eIF2gamma is a suitable phylogenetic marker for eukaryotes. We determined or completed the sequences of 13 arthropod eIF2gamma genes. Analyzing the phylogenetic distribution of 52 different intron positions in 55 distantly related eIF2gamma genes, we identified ancient ones and shared derived introns in our data set. Obviously, intron positioning in eIF2gamma is evolutionarily conserved. However, there were episodes of complete and partial intron losses followed by intron gains. We identified 17 clusters of intron positions based on their distribution. The evolution of these clusters appears to be connected with preferred exon length and can be used to estimate the relative timing of intron gain because nearby precursor introns had to be erased from the gene before the new introns could be inserted. Moreover, we identified a putative case of intron sliding that constitutes a synapomorphic character state supporting monophyly of Coleoptera, Lepidoptera, and Diptera excluding Hymenoptera. We also performed tree reconstructions using the eIF2gamma protein sequences and intron positioning as phylogenetic information. Our results support the monophyly of Viridoplantae, Ascomycota, Homobasidiomyceta, and Apicomplexa.     

Divergent Histories of rDNA Group I Introns in the Lichen Family Physciaceae

The wide but sporadic distribution of group I introns in protists, plants, and fungi, as well as in eubacteria, likely resulted from extensive lateral transfer followed by differential loss. The extent of horizontal transfer of group I introns can potentially be determined by examining closely related species or genera. We used a phylogenetic approach with a large data set (including 62 novel large subunit [LSU] rRNA group I introns) to study intron movement within the monophyletic lichen family Physciaceae. Our results show five cases of horizontal transfer into homologous sites between species but do not support transposition into ectopic sites. This is in contrast to previous work with Physciaceae small subunit (SSU) rDNA group I introns where strong support was found for multiple ectopic transpositions. This difference in the apparent number of ectopic intron movements between SSU and LSU rDNA genes may in part be explained by a larger number of positions in the SSU rRNA, which can support the insertion and/or retention of group I introns. In contrast, we suggest that the LSU rRNA may have fewer acceptable positions and therefore intron spread is limited in this gene. Reviewing Editor: Dr. W. Ford Doolittle     

An analysis of intron positions in relation to nucleotides, amino acids, and protein secondary structure

We present an analysis of intron positions in relation to nucleotides, amino acid residues, and protein secondary structure. Previous work has shown that intron sites in proteins are not randomly distributed with respect to secondary structures. Here we show that this preference can be almost totally explained by the nucleotide bias of splice site machinery, and may well not relate to protein stability or conformation at all. Each intron phase is preferentially associated with its own set of residues: phase 0 introns with lysine, glutamine, and glutamic acid before the intron, and valine after; phase 1 introns with glycine, alanine, valine, aspartic acid, and glutamic acid; and phase 2 introns with arginine, serine, lysine, and tryptophan. These preferences can be explained principally on the basis of nucleotide bias at intron locations, which is in accordance with previous literature. Although this work does not prove that introns are inserted into genomes at specific proto-splice sites, it shows that the nucleotide bias surrounding introns, however it originally occurred, explains the observed correlations between introns and protein secondary structure.     

Nucleotide sequence of the triosephosphate isomerase gene from Aspergillus nidulans: implications for a differential loss of introns

A functional cDNA from Aspergillus nidulans encoding triosephosphate isomerase (TPI) was isolated by its ability to complement a tpi1 mutation in Saccharomyces cerevisiae. This cDNA was used to obtain the corresponding gene, tpiA. Alignment of the cDNA and genomic DNA nucleotide sequences indicated that tpiA contains five introns. The intron positions in the tpiA gene were compared with those in the TPI genes of human, chicken, and maize. One intron is present at an identical position in all four organisms, two other introns are located in similar positions in A. nidulans and maize, and the remaining two introns are unique to A. nidulans. These Aspergillus-specific introns are located in regions of the protein that were predicted to be interrupted by introns based on analysis of a Go plot of chicken TPI. These comparisons are discussed in relation to the evolution of introns within TPI genes.     

USING RDNA GENES TO UNDERSTAND THE ORIGIN AND EVOLUTION OF SPLICEOSOMAL AND GROUP I INTRONS

Ribosomal DNA genes in lichen algae and lichen fungi are astonishingly rich in spliceosomal and group I introns. We use phylogenetic, secondary structure, and biochemical analyses to understand the evolution of these introns. Despite the widespread distribution of spliceosomal introns in nuclear pre-mRNA genes, their general mechanism of origin remains an open question because few proven cases of recent and pervasive intron origin have been documented. The lichen introns are valuable in this respect because they are undoubtedly of a "recent" origin and limited to the Euascomycetes. Our analyses suggest that rDNA spliceosomal introns have arisen through aberrant reverse-splicing (in trans) of free pre-mRNA introns into r RNAs. We propose that the spliceosome itself (and not an external agent; e.g. transposable elements, group II introns) has given rise to the introns. The rDNA introns are found most often between the flanking sequence G (78%) - intron-G (72%), and their clustered positions on secondary structures suggest that particular r RNA regions are preferred sites (i.e., proto-splice sites) for insertion. Mapping of intron positions on the newly available tertiary structures show that they are found most often in exposed regions of the ribosomes. This again is consistent with an intron origin through reverse-splicing. Remarkably, the distribution and phylogenetic relationships of most group I introns in nuclear rDNA genes are also consistent with a reverse-splicing origin. These data underline the value of lichens as a model system for understanding intron origin and stress the importance of RNA-level processes in the spread of these sequences in nuclear coding regions.     

Introns and reading frames: correlation between splicing sites and their codon positions

Computer analyses of the entire GenBank database were conducted to examine correlation between splicing sites and codon positions in reading frames. Intron insertion patterns (i.e., splicing site locations with respect to codon positions) have been analyzed for all of the 74 codons of all the eukaryote taxonomic groups: primates, rodents mammals, vertebrates, invertebrates, and plants. We found that reading frames are interrupted by an intron at a codon boundary (as opposed to the middle of a codon) significantly more often than expected. This observation is consistent with the exon shuffling hypothesis, because exons that end at codon boundaries can be concatenated without causing a frame shift and thus are evolutionarily advantageous. On the other hand, when introns interrupt at the middles of codons, they exist in between the first and second bases much more frequently than between the second and third bases, despite the fact that boundaries between the first and second bases of codons are generally far more important than those between the second and third bases. The reason for this is not clear and yet to be explained. We also show that the length of an exon is a multiple of 3 more frequently than expected. Furthermore, the total length of two consecutive exons is also more frequently a multiple of 3. All the observations above are consistent with results recently published by Long, Rosenberg, and Gilbert (1995).      

Two group I introns with long internal open reading frames in the chloroplast psbA gene of Chlamydomonas moewusii.

We report the nucleotide sequence of the chloroplast psbA gene encoding the 32 kilodalton protein of photosystem II from Chlamydomonas moewusii. Like its land plant homologues, this green algal protein consists of 353 amino acids. The C. moewusii psbA gene is composed of three exons containing 252, 11 and 90 codons and of two group I introns containing 2363 and 1807 nucleotides. Each of the introns features an internal open reading frame (ORF) that potentially encodes a basic protein of more than 300 residues. The primary sequences of the putative intron-encoded proteins are unrelated and none of them shares conserved elements with any of the proteins predicted from the group I intron sequences published so far. The first C. moewusii intron is inserted at the same position as the fourth intron of the psbA gene from Chlamydomonas reinhardtii; the second intron lies at a novel site downstream of this position. On the basis of their RNA secondary structures, the C. moewusii introns 1 and 2 can be assigned to subgroups IA and IB, respectively. However, intron 1 is not typical of subgroup IA introns, its most unusual feature being the location of the ORF in the "loop L5" region. To our knowledge, this is the first time that an ORF is located in this region of the group I intron structure.     

Evolution of the intron-exon structure of eukaryotic genes

The origin and evolution of intron-exon structures continue to be controversial topics. Two alternative theories, the ‘exon theory of genes’ and the ‘insertional theory of introns’, debate the presence or absence of introns in primordial genes. Both sides of the argument have focused on the positions of introns with respect to protein and gene structures. A new approach has emerged in the study of the evolution of intron-exon structures: a population analysis of genes. One example is the statistical analysis of intron phases — the position of introns within or between codons. This analysis detected a significant signal of exon shuffling in the DNA sequence database containing both ancient and modern exon sequences: intron phase correlations, that is, the association together within genes of introns of the same phase. The results of this analysis suggest that exon shuffling played an important role in the origin of both ancient and modern genes.     

Generation of a database containing discordant intron positions in eukaryotic genes (MIDB)

MOTIVATION: Intron sliding is the relocation of intron-exon boundaries over short distances and is often also referred to as intron slippage or intron migration or intron drift. We have generated a database containing discordant intron positions in homologous genes (MIDB--Mismatched Intron DataBase). Discordant intron positions are those that are either closely located in homologous genes (within a window of 10 nucleotides) or an intron position that is present in one gene but not in any of its homologs. The MIDB database aims at systematically collecting information about mismatched introns in the genes from GenBank and organizing it into a form useful for understanding the genomics and dynamics of introns thereby helping understand the evolution of genes. RESULTS: Intron displacement or sliding is critically important for explaining the present distribution of introns among orthologous and paralogous genes. MIDB allows examining of intron movements and allows mapping of intron positions from homologous proteins onto a single sequence. The database is of potential use for molecular biologists in general and for researchers who are interested in gene evolution and eukaryotic gene structure. Partial analysis of this database allowed us to identify a few putative cases of intron sliding. AVAILABILITY: http://intron.bic.nus.edu.sg/midb/midb.html     

Evolutionary implications of intron–exon distribution and the properties and sequences of the RPL10A gene in eukaryotes

The RPL10A gene encodes the RPL10 protein, required for joining 40S and 60S subunits into a functional 80S ribosome. This highly conserved gene, ubiquitous across all eukaryotic super-groups, is characterized by a variable number of spliceosomal introns, present in most organisms. These properties facilitate the recognition of orthologs among distant taxa and thus comparative studies of sequences as well as the distribution and properties of introns in taxonomically distant groups of eukaryotes. The present study examined the multiple ways in which RPL10A conservation vs. sequence changes in the gene over the course of evolution, including in exons, introns, and the encoded proteins, can be exploited for evolutionary analysis at different taxonomic levels. At least 25 different positions harboring introns within the RPL10A gene were determined in different taxa, including animals, plants, fungi, and alveolates. Generally, intron positions were found to be well conserved even across different kingdoms. However, certain introns seemed to be restricted to specific groups of organisms. Analyses of several properties of introns, including insertion site, phase, and length, along with exon and intron GC content and exon–intron boundaries, suggested biases within different groups of organisms. The use of a standard primer pair to analyze a portion of the intron-containing RPL10A gene in 12 genera of green algae within Chlorophyta is presented as a case study for evolutionary analyses of introns at intermediate and low taxonomic levels. Our study shows that phylogenetic reconstructions at different depths can be achieved using RPL10A nucleotide sequences from both exons and introns as well as the amino acid sequences of the encoded protein.     

<Emphasis Type="Italic">Oscheius tipulae</Emphasis> as an Example of eEF1A Gene Diversity in Nematodes

We characterized four eEF1A genes in the alternative rhabditid nematode model organism Oscheius tipulae. This is twice the copy number of eEF1A genes in C. elegans, C. briggsae, and, probably, many other free-living and parasitic nematodes. The introns show features remarkably different from those of other metazoan eEF1A genes. Most of the introns in the eEF1A genes are specific to O. tipulae and are not shared with any of the other genes described in metazoans. Most of the introns are phase 0 (inserted between two codons), and few are inserted in protosplice sites (introns inserted between the nucleotide sequence A/CAG and G/A). Two of these phase 0 introns are conserved in sequence in two or more of the four eEF1A gene copies, and are inserted in the same position in the genes. Neither of these characteristics has been detected in any of the nematode eEF1A genes characterized to date. The coding sequences were also compared with other eEF1A cDNAs from 11 different nematodes to determine the variability of these genes within the phylum Nematoda. Parsimony and distance trees yielded similar topologies, which were similar to those created using other molecular markers. The presence of more than one copy of the eEF1A gene with nearly identical coding regions makes it difficult to define the orthologous cDNAs. As shown by our data on O. tipulae, careful and extensive examination of intron positions in the eEF1A gene across the phylum is necessary to define their potential for use as valid phylogenetic markers.     

Association of intron phases with conservation at splice site sequences and evolution of spliceosomal introns.

How exon-intron structures of eukaryotic genes evolved under various evolutionary forces remains unknown. The phases of spliceosomal introns (the placement of introns with respect to reading frame) provide an opportunity to approach this question. When a large number of nuclear introns in protein-coding genes were analyzed, it was found that most introns were of phase 0, which keeps codons intact. We found that the phase distribution of spliceosomal introns is strongly correlated with the sequence conservation of splice signals in exons; the relatively underrepresented phase 2 introns are associated with the lowest conservation, the relatively overrepresented phase 0 introns display the highest conservation, and phase 1 introns are intermediate. Given the detrimental effect of mutations in exon sequences near splice sites as found in molecular experiments, the underrepresentation of phase 2 introns may be the result of deleterious-mutation-driven intron loss, suggesting a possible genetic mechanism for the evolution of intron-exon structures.     

Evidence for extensive recent intron transposition in closely related fungi

Though spliceosomal introns are a major structural component of most eukaryotic genes and intron density varies by more than three orders of magnitude among eukaryotes [1-3], the origins of introns are poorly understood, and only a few cases of unambiguous intron gain are known [4-8]. We utilized population genomic comparisons of three closely related fungi to identify crucial transitory phases of intron gain and loss. We found 74 intron positions showing intraspecific presence-absence polymorphisms (PAPs) for the entire intron. Population genetic analyses identified intron PAPs at different stages of fixation and showed that intron gain or loss was very recent. We found direct support for extensive intron transposition among unrelated genes. A substantial proportion of highly similar introns in the genome either were recently gained or showed a transient phase of intron PAP. We also identified an intron transfer among paralogous genes that created a new intron. Intron loss was due mainly to homologous recombination involving reverse-transcribed mRNA. The large number of intron positions in transient phases of either intron gain or loss shows that intron evolution is much faster than previously thought and provides an excellent model to study molecular mechanisms of intron gain.