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.     

ExInt: an Exon/Intron database

The Exon/Intron (ExInt) database incorporates information on the exon/intron structure of eukaryotic genes. Features in the database include: intron nucleotide sequence, amino acid sequence of the corresponding protein, position of the introns at the amino acid level and intron phase. From ExInt, we have also generated four additional databases each with ExInt entries containing predicted introns, introns experimentally defined, organelle introns or nuclear introns. ExInt is accessible through a retrieval system with pointers to GenBank. The database can be searched by keywords, locus name, NID, accession number or length of the protein. ExInt is freely accessible at http://intron.bic.nus.edu.sg/exint/exint.html     

A structural and phylogenetic analysis of the group IC1 introns in the order Bangiales (Rhodophyta)

Our previous study of the North American biogeography of Bangia revealed the presence of two introns inserted at positions 516 and 1506 in the nuclear-encoded SSU rRNA gene. We subsequently sequenced nuclear SSU rRNA in additional representatives of this genus and the sister genus Porphyra in order to examine the distribution, phylogeny, and structural characteristics of these group I introns. The lengths of these introns varied considerably, ranging from 467 to 997 nt for intron 516 and from 509 to 1,082 nt for intron 1506. The larger introns contained large insertions in the P2 domain of intron 516 and the P1 domain of intron 1506 that correspond to open reading frames (ORFs) with His-Cys box homing endonuclease motifs. These ORFs were found on the complementary strand of the 1506 intron in Porphyra fucicola and P. umbilicalis (HG), unlike the 516 intron in P. abbottae, P. kanakaensis, P. tenera (SK), Bangia fuscopurpurea (Helgoland), and B. fuscopurpurea (MA). Frameshifts were noted in the ORFs of the 516 introns in P. kanakaensis and B. fuscopurpurea (HL), and all ORFs terminated prematurely relative to the amino acid sequence for the homing endonuclease I-Ppo I. This raises the possibility that these sequences are pseudogenes. Phylogenies generated using sequences of both introns and the 18S rRNA gene were congruent, which indicated long-term immobility and vertical inheritance of the introns followed by subsequent loss in more derived lineages. The introns within the florideophyte species Hildenbrandia rubra (position 1506) were included to determine relationships with those in the Bangiales. The two sequences of intron 1506 analyzed in Hildenbrandia were positioned on a well-supported branch associated with members of the Bangiales, indicating possible common ancestry. Structural analysis of the intron sequences revealed a signature structural feature in the P5b domain of intron 516 that is unique to all Bangialean introns in this position and not seen in intron 1506 or other group IC1 introns.     

Distinct sites of phosphorothioate substitution interfere with folding and splicing of the Anabaena group I intron

Although the active site of group I introns is phylogenetically conserved, subclasses of introns have evolved different mechanisms of stabilizing the catalytic core. Large introns contain weakly conserved 'peripheral' domains that buttress the core through predicted interhelical contacts, while smaller introns use loop-helix interactions for stability. In all cases, specific and non-specific magnesium ion binding accompanies folding into the active structure. Whether similar RNA-RNA and RNA-magnesium ion contacts play related functional roles in different introns is not clear, particularly since it can be difficult to distinguish interactions directly involved in catalysis from those important for RNA folding. Using phosphorothioate interference with RNA activity and structure in the small (249 nt) group I intron from Anabaena, we used two independent assays to detect backbone phosphates important for catalysis and those involved in intron folding. Comparison of the interference sites identified in each assay shows that positions affecting catalysis cluster primarily in the conserved core of the intron, consistent with conservation of functionally important phosphates, many of which are magnesium ion binding sites, in diverse group I introns, including those from Azoarcus and Tetrahymena. However, unique sites of folding interference located outside the catalytic core imply that different group I introns, even within the same subclass, use distinct sets of tertiary interactions to stabilize the structure of the catalytic core.     

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.     

Mammalian Male Mutation Bias: Impacts of Generation Time and Regional Variation in Substitution Rates

In mammals, males undergo a greater number of germline cell divisions compared with females. Thus, the male germline accumulates more DNA replication errors, which result in male mutation bias—a higher mutation rate for males than for females. The phenomenon of male mutation bias has been investigated mostly for rodents and primates, however, it has not been studied in detail for other mammalian orders. Here we sequenced and analyzed five introns of three genes (DBX/DBY, UTX/UTY, and ZFX/ZFY) homologous between X and Y chromosomes in several species of perissodactyls (horses and rhinos) and of primates. Male mutation bias was evident: substitution rate was higher for a Y chromosome intron than for its X chromosome homologue for all five intron pairs studied. Substitution rates varied regionally among introns sequenced on the same chromosome and this variation influenced male mutation bias inferred from each intron pair. Interestingly, we observed a positive correlation in substitution rates between homologous X and homologous Y introns as well as between orthologous primate and perissodactyl introns. The male-to-female mutation rate ratio estimated from concatenated sequences of five perissodactyl introns was 3.88 (95% CI = 2.90–6.07). Using the data generated here and estimates available in the literature, we compared male mutation bias among several mammalian orders. We conclude that male mutation bias is significantly higher for organisms with long generation times (primates, perissodactyls, and felids) than for organisms with short generation times (e.g., rodents) since the former undergo a greater number of male germline cell divisions. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Deborah Charlesworth]     

Remarkable sequence conservation of the last intron in the PKD1 gene

The last intron of the PKD1 gene (intron 45) was found to have exceptionally high sequence conservation across four mammalian species: human, mouse, rat, and dog. This conservation did not extend to the comparable intron in pufferfish. Pairwise comparisons for intron 45 showed 91% identity (human vs. dog) to 100% identity (mouse vs. rat) for an average for all four species of 94% identity. In contrast, introns 43 and 44 of the PKD1 gene had average pairwise identities of 57% and 54%, and exons 43, 44, and 45 and the coding region of exon 46 had average pairwise identities of 80%, 84%, 82%, and 80%. Intron 45 is 90 to 95 bp in length, with the major region of sequence divergence being in a central 4-bp to 9-bp variable region. RNA secondary structure analysis of intron 45 predicts a branching stem-loop structure in which the central variable region lies in one loop and the putative branch point sequence lies in another loop, suggesting that the intron adopts a specific stem-loop structure that may be important for its removal. Although intron 45 appears to conform to the class of small, G-triplet-containing introns that are spliced by a mechanism utilizing intron definition, its high sequence conservation may be a reflection of constraints imposed by a unique mechanism that coordinates splicing of this last PKD1 intron with polyadenylation.     

A self-splicing group I intron in the DNA polymerase gene of Bacillus subtilis bacteriophage SPO1

We report a self-splicing intron in bacteriophage SPO1, whose host is the gram-positive Bacillus subtilis. The intron contains all the conserved features of primary sequence and secondary structure previously described for the group IA introns of eukaryotic organelles and the gram-negative bacteriophage T4. The SPO1 intron contains an open reading frame of 522 nucleotides. As in the T4 introns, this open reading frame begins in a region that is looped out of the secondary structure, but ends in a highly conserved region of the intron core. The exons encode SPO1 DNA polymerase, which is highly similar to E. coli DNA polymerase I. The demonstration of self-splicing introns in viruses of both gram-positive and gram-negative eubacteria lends further evidence for their early origin in evolution.     

Self-splicing of group II introns in vitro: mapping of the branch point and mutational inhibition of lariat formation

Group II intron bl1 from yeast mitochondria can undergo self-splicing in vitro. Exons become correctly ligated, and the excised intron has a lariat structure similar to that of introns from nuclear mRNA. The branch point of the bl1 lariat is located eight or nine nucleotides upstream of the 3' end of the intron and is part of a hairpin structure that is well conserved among group II introns. Several mutations next to the branch point and in other parts of the core structure of group II introns are shown to affect lariat formation. One of them, carried by strain M4873, abolishes splicing in vivo and in vitro, apparently by changing the architecture of the hairpin structure containing the branch point. Similarities between group II introns and nuclear pre-mRNA introns are discussed in terms of evolutionary relatedness.     

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.     

The mitochondrial genome of the fission yeast Schizosaccharomyces pombe. The cytochrome b gene has an intron closely related to the first two introns in the Saccharomyces cerevisiae cox1 gene

The DNA sequence of the cob region of the Schizosaccharomyces pombe mitochondrial DNA has been determined. The cytochrome b structural gene is interrupted by an intron of 2526 base-pairs, which has an open reading frame of 2421 base-pairs in phase with the upstream exon. The position of the intron differs from those found in the cob genes of Saccharomyces cerevisiae, Aspergillus nidulans or Neurospora crassa. The Sch. pombe cob intron has the potential of assuming an RNA secondary structure almost identical to that proposed for the first two cox1 introns (group II) in S. cerevisiae and the p1-cox1 intron in Podospora anserina. It has most of the consensus nucleotides in the central core structure described for this group of introns and its comparison with other group II introns allows the identification of an additional conserved nucleotide stretch. A comparison of the predicted protein sequences of group II intronic coding regions reveals three highly conserved blocks showing pairwise amino acid identities of 34 to 53%. These regions comprise over 50% of the coding length of the intron but do not include the 5' region, which has strong secondary structural features. In addition to the potential intron folding, long helical structures involving repetitive sequences can be formed in the flanking cob exon regions. A comparison of the Sch. pombe cytochrome b sequence with those available from other organisms indicates that Sch. pombe is evolutionarily distant from both budding yeasts and filamentous fungi. As was seen for the Sch. pombe cox1 gene (Lang, 1984), the cob exons are translated using the universal genetic code and this distinguishes Sch. pombe mitochondria from all other fungal and animal mitochondrial systems.     

The biology of intron gain and loss

Intron density in eukaryote genomes varies by more than three orders of magnitude, so there must have been extensive intron gain and/or intron loss during evolution. A favored and partial explanation for this range of intron densities has been that introns have accumulated stochastically in large eukaryote genomes during their evolution from an intron-poor ancestor. However, recent studies have shown that some eukaryotes lost many introns, whereas others accumulated and/or gained many introns. In this article, we discuss the growing evidence that these differences are subject to selection acting on introns depending on the biology of the organism and the gene involved.     

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.     

Long-term evolution of the S788 fungal nuclear small subunit rRNA group I introns

More than 1000 group I introns have been identified in fungal rDNA. Little is known, however, of the splicing and secondary structure evolution of these ribozymes. Here, we use a combination of comparative and biochemical methods to address the evolution and splicing of a vertically inherited group I intron found at position 788 in the fungal small subunit (S) rRNA. The ancestral state of the S788 intron contains a highly conserved core and an extended P5 domain typical of IC1 introns. In contrast, the more derived introns have lost most of P5, and have an accelerated divergence rate within the core region with three functionally important substitutions that unambiguously separate them from the ancestral pool. Of 14 S788 group I introns that were tested for splicing, five, all of the ancestral type, were able to self-splice and produced intron RNA circles in vitro. The more derived S788 introns did not self-splice, and potentially rely on fungal-specific factors to facilitate splicing. In summary, we demonstrate one possible fate of vertically inherited group I introns, the loss of secondary structure elements, lessened selective constraints in the intron core, and ultimately, dependence on host-mediated splicing.     

Widespread intron loss suggests retrotransposon activity in ancient apicomplexans

Several facets of spliceosomal intron in apicomplexans remain mysterious. First, intron numbers vary across species by 2 orders of magnitude, indicating massive intron loss and/or gain. Second, previous studies have shown very different evolutionary patterns over different timescales, with 100-fold higher rates of intron loss/gain between genera than within genera. Third, the timing and dynamics of nearly complete intron loss in Cryptosporidium species, as well as reasons for retention of the few remaining introns, remain unknown. We compared intron positions in 785 orthologous genes between 3 moderate to intron-rich apicomplexan species. We estimate that the Theileria-Plasmodium ancestor had 4.5 times as many introns as modern Plasmodium species and 38% more than modern Theileria species, and that subsequent intron losses have outnumbered intron gains by 5.8 to 1 in Theileria and by some 56 to 1 in Plasmodium. Several patterns suggest that these intron losses occurred by recombination with reverse-transcribed mRNAs. Intriguingly, this finding suggests significant retrotransposon activity in the lineages leading to both Theileria and Plasmodium, in contrast to the dearth of known retrotransposons and intron loss within modern species from both genera. We also compared genomes from Cryptosporidium parvum and C. hominis and found no evidence of ongoing intron loss, nor of intron gain. By contrast, Cryptosporidium introns are less evolutionary conserved with Toxoplasma than are introns from other apicomplexans; thus the few remaining introns are not simply indispensable ancestral introns.     

Genomic structure of the sandworm, Perinereis vancaurica tetradentata, troponin C

Troponin C (TnC) superfamily genes essentially possess five introns, the positions of all but the fourth being highly conserved. The fourth intron is frequently absent from protostomian invertebrate genes, such as calmodulin or TnC. We previously proposed that the common ancestor of TnC superfamily genes never possessed an intron corresponding to today's fourth introns, and that members of the superfamily independently gained a fourth intron in the evolutionary pathway of each lineage. In the present study, we isolated the TnC cDNA from the sandworm, Perinereis vancaurica tetradentata and determined its genomic structure. Sandworm TnC appears to exist as a single copy gene consisting of six exons and five introns. The positions of the first, second, third and fifth introns are identical to other TnCs, but that of the fourth intron is unique. This is in good agreement with the above-mentioned scheme, i.e. the gain of the fourth intron of sandworm TnC might have occurred within the annelid lineage after annelida/mollusca divergence.