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 organization of the 5' region of the thyroglobulin gene. Evidence for intron loss and "exonization" during evolution

More than one third of thyroglobulin (1190 residues out of 2750) is made of one peptide motif repeated ten times in tandem. Segments unrelated to the motif interrupt this structure at various places. The corresponding gene region, which extends over 40 x 10(3) bases, was studied in detail. All exon borders and exon/intron junctions were localized precisely and sequenced, and their positions were correlated with the repetitive organization of the protein. When intron positions were compiled on a consensus sequence of all repeats, three categories of introns were observed. Except between repeats numbers 5 and 6, an intron was invariably found within the Cys codon making the limit of each motif. This category of intron most probably reflects the serial duplication events responsible for the evolution of this region of the gene. All other introns, except no. 2, are found at positions were the repetitive structure is disrupted by "inserted" peptides. We present the hypothesis that this second category of introns was already present in the original unit before the first duplication. Thereafter, they would have experienced either complete loss (some units do not contain any intron) or partial or total exonization, resulting in the slipping of intronic material into coding sequence. Intron no. 2, finally, separates motif no. 1 at a position on the boundary between two segments presenting sequence homology. This last type of intron probably reflects an initial duplication event at the origin of a primordial thyroglobulin gene motif. With all these characteristics, the thyroglobulin gene is presented as a paradigm for the analysis of the fate of introns in gene evolution.     

Secondary structural features of modules M2 and M3 of barnase in solution by NMR experiment and distance geometry calculation

Proteins consist of structural units such as globular domains, secondary structures, and modules. Modules were originally defined by partitioning a globular domain into compact regions, each of which is a contiguous polypeptide segment having a compact conformation. Since modules show close correlations with the intron positions of genes, they are regarded as primordial polypeptide pieces encoded by exons and shuffled, leading to yield new combination of them in early biological evolution. Do modules maintain their native conformations in solution when they are excised at their boundaries? In order to find answers to this question, we have synthesized modules of barnase, one of the bacterial RNases, and studied the solution structures of modules M2 (amino acid residues 24–52) and M3 (52–73) by 2D NMR studies. Some local secondary structures, α-helix, and β-turns in M2 and β-turns in M3, were observed in the modules at the similar positions to those in the intact barnase but the overall state seems to be in a mixture of random and native conformations. The present result shows that the excised modules have propensity to form similar secondary structures to those of the intact barnase. © 1993 Wiley-Liss, Inc.     

Evidence that human genes of modular proteins have retained significantly more ancestral introns than their fly or worm orthologues

Comparison of the exon-intron structures of human, fly and worm orthologues of mosaic genes assembled from class 1-1 modules by exon-shuffling has revealed that human genes retained significantly more of the original inter-module introns than their protostome orthologues. It is suggested that the much higher rate of intron loss in the worm- and insect lineages than in the chordate lineage reflects their greater tendency for genome compaction.     

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     

Distribution of rRNA introns in the three-dimensional structure of the ribosome

More than 1200 introns have been documented at over 150 unique sites in the small and large subunit ribosomal RNA genes (as of February 2002). Nearly all of these introns are assigned to one of four main types: group I, group II, archaeal and spliceosomal. This sequence information has been organized into a relational database that is accessible through the Comparative RNA Web Site (http://www.rna.icmb.utexas.edu/) While the rRNA introns are distributed across the entire tree of life, the majority of introns occur within a few phylogenetic groups. We analyzed the distributions of rRNA introns within the three-dimensional structures of the 30S and 50S ribosomes. Most sites in rRNA genes that contain introns contain only one type of intron. While the intron insertion sites occur at many different coordinates, the majority are clustered near conserved residues that form tRNA binding sites and the subunit interface. Contrary to our expectations, many of these positions are not accessible to solvent in the mature ribosome. The correlation between the frequency of intron insertions and proximity of the insertion site to functionally important residues suggests an association between intron evolution and rRNA function.     

Correlation of exons with structural domains in alcohol dehydrogenase.

The intron/exon arrangement in the gene sequence of maize alcohol dehydrogenase has been compared to the three dimensional structure of liver alcohol dehydrogenase. The co-enzyme binding domain is separated from the catalytic domain by introns four and nine. Intron seven separates the co-enzyme binding domain into two structurally similar mononucleotide binding units. The first of these units is divided by introns five and six into three structurally similar alpha beta modules. Implications of these results for protein evolution is discussed. All splice junctions map close to or at the surface of the domains, and several of these cannot be identified by distance maps.     

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.     

Eu-Chlorella large subunit rDNA sequences and group I introns in ribosomal DNA of the paramecian symbiotic alga NC64A

Although the examination of large subunit ribosomal RNA genes (LSU rDNA) is advanced in phylogenetic studies, no corresponding sequence data from trebouxiophytes have been published, with the exception of ‘Chlorella’ellipsoidea Gerneck. We determined the LSU rDNA sequence of Chlorella vulgaris Beijerinck and of the symbiotic alga of green paramecium, Chlorella sp. NC64A. A total of 59 nucleotide substitutions were found in the LSU rDNA of the two species, which are disproportionately distributed. Primarily, 65% of the substitutions were encountered in the first 800 bp of the alignment. This segment apparently has evolved eight times faster than the complete SSU rDNA sequence, making it a good candidate for a phylogenetic marker and giving a resolution level intermediate between small subunit (SSU) rDNA and internal transcribed spacers. Green algae are known as a group I intron‐rich group along with rhodophytes and fungi. NC64A is particularly rich in the introns; five introns were newly identified from the LSU rDNA sequence, which we named Cnc.L200, Cnc.L1688, Cnc.L1926, Cnc.L2184 and Cnc.L2437, following the insertion positions. In the present study we analyzed these introns with three others (Cnc.S943, Cnc.S1367 and Cnc.S1512) that had already been found in NC64A SSU rDNA. Secondary structure modeling placed these introns in the group I intron family, with four introns belonging to subgroup C1 and the other four introns belonging to subgroup E. Five of the intron insertion positions are unique to the paramecian symbiont, which may indicate relatively recent events of intron infections that includes transpositions. Intron phylogeny showed unprecedented relationships; four Cnc. IC1 introns made a clade with some green algal introns with insertions at nine different positions, whereas four Cnc. IE introns made a clade with the S651 intron (Chlorella sp. AN 1–3), which lay as a sister to the S516 insertion position subfamily.     

Analysis of ribosomal protein gene structures: implications for intron evolution

Many spliceosomal introns exist in the eukaryotic nuclear genome. Despite much research, the evolution of spliceosomal introns remains poorly understood. In this paper, we tried to gain insights into intron evolution from a novel perspective by comparing the gene structures of cytoplasmic ribosomal proteins (CRPs) and mitochondrial ribosomal proteins (MRPs), which are held to be of archaeal and bacterial origin, respectively. We analyzed 25 homologous pairs of CRP and MRP genes that together had a total of 527 intron positions. We found that all 12 of the intron positions shared by CRP and MRP genes resulted from parallel intron gains and none could be considered to be “conserved,” i.e., descendants of the same ancestor. This was supported further by the high frequency of proto-splice sites at these shared positions; proto-splice sites are proposed to be sites for intron insertion. Although we could not definitively disprove that spliceosomal introns were already present in the last universal common ancestor, our results lend more support to the idea that introns were gained late. At least, our results show that MRP genes were intronless at the time of endosymbiosis. The parallel intron gains between CRP and MRP genes accounted for 2.3% of total intron positions, which should provide a reliable estimate for future inferences of intron evolution.     

The structural organization of the ascidian, Halocynthia roretzi, calmodulin genes. The vicissitude of introns during the evolution of calmodulin genes

Two distinct calmodulin (CaM) genes are isolated from the ascidian, Halocynthia roretzi, (Hr-CaM A and Hr-CaM B) and those structures are determined. There are three nucleotide substitutions, producing two amino acid differences between Hr-CaM A and Hr-CaM B, and those are corresponding to two of the known eight variable residues among metazoan CaMs. Both Hr-CaM A and Hr-CaM B are constructed from six exons and five introns, and the positions of introns are identical. The positions of introns of Hr-CaMs are also identical with those of vertebrate CaMs, except third introns. The third introns of Hr-CaMs are inserted at 28bp upstream when compared with vertebrate CaMs. Thus, sliding of the third intron might have occurred in only the ascidian lineage prior to the gene duplication that also occurred only in that lineage. In addition, with the comparison of the intron positions, we attempt to investigate the vicissitude of introns during the evolution of metazoan CaMs.     

Centripetal modules and ancient introns

We have created an algorithm which instantiates the centripetal definition of modules, compact regions of protein structure, as introduced by Go and Nosaka (M. Go and M. Nosaka, 1987. Protein architecture and the origin of introns. Cold Spring Harbor Symp. Quant. Bio. 52, 915-924). That definition seeks the minima of a function that sums the squares of C-alpha carbon distances over a window around each amino acid residue in a three-dimensional protein structure and identifies such minima with module boundaries. We analyze a set of 44 ancient conserved proteins, with known three-dimensional structures, which have intronless homologues in bacteria and intron-containing homologues in the eukaryotes, with a corresponding set of 988 intron positions. We show that the phase zero intron positions are significantly correlated with the module boundaries (p = 0.0002), while the intron positions that lie within codons, in phase one and phase two, are not correlated with these 'centripetal' module boundaries. Furthermore, we analyze the phylogenetic distribution of intron positions and identify a subset of putatively 'ancient' intron positions: phase zero positions in one phylogenetic kingdom which have an associated intron either in an identical position or within three codons in another phylogenetic kingdom (a notion of intron sliding). This subset of 120 'ancient' introns lies closer to the module boundaries than does the full set of phase zero introns with high significance, a p-value of 0.008. We conclude that the behavior of this set of introns supports the prediction of a mixed theory: that some introns are very old and were used for exon shuffling in the progenote, while many introns have been lost and added since.     

The recent origins of introns

Accumulating evidence that introns are highly restricted in their phylogenetic distribution strongly supports the view that introns were inserted late in eukaryotic evolution into preformed genes and, hence, that exon-shuffling played no role in the assembly of primordial genes. Potential mechanisms of intron insertion and the possible evolution of nuclear introns and their splicing machinery from self-splicing group II introns are also discussed.     

Heterogeneous yet similar introns reside in identical positions of the rRNA genes in natural isolates of the archaeon Aeropyrum pernix

Some archaeal ribosomal DNA (rDNA) introns carry homing endonuclease-like genes and are therefore assumed to propagate by "intron homing". A previous study demonstrated that three introns are located within the rRNA operon (arnSL) of Aeropyrum pernix strain K1, two of which, Ialpha and Igamma, harbor open reading frames (ORFs) encoding putative LAGLIDADG-type endonucleases. In an effort to understand further the rDNA intron distribution in natural A. pernix populations, 11 A. pernix strains were isolated from marine hydrothermal biotopes, and comparative nucleotide sequence analysis of the arnSL alleles was performed. Of the 11 isolates, eight contained multiple introns, and three patterns of intron insertion were found. Three novel introns, Idelta (62 bp in length), Ivarepsilon (122 bp) and Izeta (57 bp) were identified. They were all ORF-less, but their predicted RNA secondary structure at the exon-intron junctions was consistent with the bulge-helix-bulge motif. The insertion positions and the terminal inverted repeat sequences of Idelta and Izeta were in agreement with those of Ialpha and Igamma, respectively. This suggests that these intron variants were generated by large indels (insertions/deletions) during their evolution.     

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.     

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.     

Polygalacturonase-inhibiting proteins--leucine-rich repeat proteins in plant defence

Plant polygalacturonase-inhibiting proteins (PGIPs) belong to the leucine-rich repeat (LRR) family and are known to prevent pathogen invasion by inhibiting the plant cell wall degrading enzyme, polygalacturonase. Our study reveals that these multigene-encoded defence proteins found in flowering plants only exhibit identical domain architecture with 10 tandemly-arranged LRRs. This implies that variations of PGIP inhibitory properties are not associated with the number of the repeats but with subtle changes in the sequence content of the repeats. The first and eighth repeat contain more mutations compared to the strict conservation of the plant-specific LRRs or any repeat at other positions. Each of these repeats forms a separate cluster in the phylogenetic tree, both within and across plant families, thus suggesting uniqueness with respect to their position. A study of the genes encoding PGIPs, shows the existence of two categories (i) single exon and hence no intron; and (ii) two exons with an intron in between. Analyses of the intron phase and correlation of the exon-intron structure with the compact structural modules in PGIPs support insertion of introns in the pre-existing single exon genes and thus the intron late model. Lack of conservation of phase across families and formation of individual clusters for each family in the phylogenetic tree drawn with the intron sequences illustrate the event of insertion that took place separately in each of these families.     

Intron Distribution in Ancient Paralogs Supports Random Insertion and Not Random Loss

The intron positions of ten different protein families were examined to determine (the statistical likelihood of) whether spliceosomal introns are the result of random insertion events into previously intronless genes, on the one hand, or the result of random loss from common ancestral introns, on the other. The number of expected matches for the alternative scenarios was calculated for a binomial distribution by considering currently observed introns relative to all possible locations for insertion or loss. Introns occurring at approximately the same location (hereafter called a ``match') were tallied for each of the paired proteins. Matches were identified by their positions in the multiple alignment and were defined as any two introns occurring within a window of 11 possible nucleotide positions, thereby allowing for possible alignment errors and ``intron sliding.' Matches were tallied from the raw data and compared with the expected number of matches for the two different scenarios. The results suggest that the distribution of introns in genes encoding proteins is due to random insertion and not random loss. Received: 8 September 1996 / Accepted: 24 January 1997     

The genomic structure of the scallop, Patinopecten yessoensis, troponin C gene: a hypothesis for the evolution of troponin C

Two cDNAs encoding troponin C (TnC) isoforms are isolated from the scallop, Patinopecten yessoensis, striated adductor muscle. The sequential differences between these isoforms, named TnC(long) and TnC(short), are restricted in several residues of the C-terminal region. TnC(long) is commonly expressed in both the striated and the smooth adductor muscle; however, TnC(short) is only in the striated adductor muscle. The TnC gene is a single copy gene in the scallop, thus they are expressed through the alternative splicing from the same gene. The scallop TnC gene is constructed from five exons and four introns, and positions of introns are identical with chordate TnC genes, although the scallop TnC possesses no corresponding intron to the fourth intron of chordates. The loss of this intron is also observed in Drosophila TnC; these may be remnants of their ancestor, namely the early metazoan TnC gene might be a five exons-four introns structure. In addition, the absence of the corresponding intron is also observed among protostomian calmodulins (CaMs), a molecule closely related to TnC. This suggests that the common ancestor gene of the TnC superfamily might also be a five exons-four introns structure. Assuming this to be true, the discordance of the fourth intron positions observed among members of the family is well explained by the evolutionary independent gain of the intron on each member's lineage.