Molecular evolution and phylogenetic utility of the petD group II intron: a case study in basal angiosperms

Sequences of spacers and group I introns in plant chloroplast genomes have recently been shown to be very effective in phylogenetic reconstruction at higher taxonomic levels and not only for inferring relationships among species. Group II introns, being more frequent in those genomes than group I introns, may be further promising markers. Because group II introns are structurally constrained, we assumed that sequences of a group II intron should be alignable across seed plants. We designed universal amplification primers for the petD intron and sequenced this intron in a representative selection of 47 angiosperms and three gymnosperms. Our sampling of taxa is the most representative of major seed plant lineages to date for group II introns. Through differential analysis of structural partitions, we studied patterns of molecular evolution and their contribution to phylogenetic signal. Nonpairing stretches (loops, bulges, and interhelical nucleotides) were considerably more variable in both substitutions and indels than in helical elements. Differences among the domains are basically a function of their structural composition. After the exclusion of four mutational hotspots accounting for less than 18% of sequence length, which are located in loops of domains I and IV, all sequences could be aligned unambiguously across seed plants. Microstructural changes predominantly occurred in loop regions and are mostly simple sequence repeats. An indel matrix comprising 241 characters revealed microstructural changes to be of lower homoplasy than are substitutions. In showing Amborella first branching and providing support for a magnoliid clade through a synapomorphic indel, the petD data set proved effective in testing between alternative hypotheses on the basal nodes of the angiosperm tree. Within angiosperms, group II introns offer phylogenetic signal that is intermediate in information content between that of spacers and group I introns on the one hand and coding sequences on the other.     

Relationship of viroids and certain other plant pathogenic nucleic acids to group I and II introns

Summary The nucleotide sequences of viroids contain features believed to be essential for the splicing of group I introns. Common sequence elements include a 16-nucleotide consensus sequence and three pairs of short sequences arranged in the same sequential order in both types of RNAs. The calculated probability of finding sequences resembling the 16-nucleotide consensus sequence in random nucleotide chains showed that at low fidelity (up to 5 mismatched nucleotides), the number of such sequences in viroids, plant viral satellite RNAs, plant viral RNAs and one plant viral DNA, group I introns and flanking exons does not significantly differ from the number expected at random. As the degree of fidelity is increased, the number in both introns and viroids, but not in exons or the other plant pathogens examined, greatly exceeds that expected in random chains. These findings suggest that viroids may have evolved from group I introns and/or that processing of viroid oligomers to monomers may have structural requirements similar to those of group I introns. The nucleotide sequences of viroids do not show close homology with two conserved regions of group II introns, the 14-base pair consensus region and the 5 terminal segment. However, close homology does exist between the conserved sequence of the 3 terminal segment of group II introns and viroids thus suggesting a possible evolutionary or functional relationship.     

Origin and evolution of the chloroplast trnK (matK) intron: a model for evolution of group II intron RNA structures

The trnK intron of plants encodes the matK open reading frame (ORF), which has been used extensively as a phylogenetic marker for classification of plants. Here we examined the evolution of the trnK intron itself as a model for group II intron evolution in plants. Representative trnK intron sequences were compiled from species spanning algae to angiosperms, and four introns were newly sequenced. Phylogenetic analyses showed that the matK ORFs belong to the ML (mitochondrial-like) subclass of group II intron ORFs, indicating that they were derived from a mobile group II intron of the class. RNA structures of the introns were folded and analyzed, which revealed progressive RNA structural deviations and degenerations throughout plant evolution. The data support a model in which plant organellar group II introns were derived from bacterial-like introns that had "standard" RNA structures and were competent for self-splicing and mobility and that subsequently the ribozyme structures degenerated to ultimately become dependent upon host-splicing factors. We propose that the patterns of RNA structure evolution seen for the trnK intron will apply to the other group II introns in plants.     

The dispersal of five group II introns among natural populations of Escherichia coli

Group II introns are self-splicing RNAs that also act as retroelements in bacteria, mitochondria, and chloroplasts. Group II introns were identified in Escherichia coli in 1994, but have not been characterized since, and, instead, other bacterial group II introns have been studied for splicing and mobility properties. Despite their apparent intractability, at least five distinct group II introns exist naturally in E. coli strains. To illuminate their function and learn how the introns have dispersed in their natural host, we have investigated their distribution in the ECOR reference collection. Two introns were cloned and sequenced to complete their partial sequences. Unexpectedly, southern blots showed all ECOR strains to contain fragments and/or full-length copies of group II introns, with some strains containing up to 15 intron copies. One intron, E.c.14, has two natural homing sites in IS629 and IS911 elements, and the intron can be present in one, both, or neither homing site in a given strain. Nearly all strains that contain full-length introns also contain unfilled homing sites, suggesting either that mobility is highly inefficient or that most full-length copies are nonfunctional. The data indicate independent mobility of the introns, as well as mobility via the host DNA elements, and overall, the pattern of intron distribution resembles that of IS elements.     

A mixed group II/group III twintron in the Euglena gracilis chloroplast ribosomal protein S3 gene: evidence for intron insertion during gene evolution.

The splicing of a 409 nucleotide intron from the Euglena gracilis chloroplast ribosomal protein S3 gene (rps3) was examined by cDNA cloning and sequencing, and northern hybridization. Based on the characterization of a partially spliced pre-mRNA, the intron was characterized as a 'mixed' twintron, composed of a 311 nucleotide group II intron internal to a 98 nucleotide group III intron. Twintron excision is via a 2-step sequential splicing pathway, with removal of the internal group II intron preceding excision of the external group III intron. Based on secondary structural analysis of the twintron, we propose that group III introns may represent highly degenerate versions of group II introns. The existence of twintrons is interpreted as evidence that group II introns were inserted during the evolution of Euglena chloroplast genes from a common ancestor with eubacteria, archaebacteria, cyanobacteria, and other chloroplasts.     

Group I introns in the liverwort mitochondrial genome: the gene coding for subunit 1 of cytochrome oxidase shares five intron positions with its fungal counterparts.

The complete nucleotide sequence of the mitochondrial DNA (mtDNA) from a liverwort, Marchantia polymorpha, contains thirty-two introns. Twenty-five of these introns possess the characteristic secondary structures and consensus sequences of group II introns. The remaining seven are group I introns, six of which happen to interrupt the gene coding for subunit 1 of cytochrome oxidase (cox1). Interestingly, the insertion sites of one group II and four group I introns in the cox1 gene coincide with those of the respective fungal mitochondrial interns. Moreover, comparison of the four group I introns with their fungal counterparts shows that group I introns inserted at identical genomic sites in different organisms are indeed related to one another, in terms of the peptide sequences generated from the complete or fragmental ORFs encoded by these introns. At the same time, the liverwort introns turned out to be more divergent from their fungal cognates than the latter are from one another. We therefore conclude that vertical transmission from a common ancestor organism is the simplest explanation for the presence of cognate introns in liverwort and fungal mitochondrial genomes.     

Variation in sequence and RNA editing within core domains of mitochondrial group II introns among plants

The 3' regions of several group II introns within the mitochondrial genes nad1 and nad7 show unexpected sequence divergence among flowering plants, and the core domains 5 and 6 are predicted to have weaker helical structure than those in self-splicing group II introns. To assess whether RNA editing improves helical stability by the conversion of A-C mispairs to A-U pairs, we sequenced RT-PCR amplification products derived from excised intron RNAs or partially spliced precursors. Only in some cases was editing observed to strengthen the predicted helices. Moreover, the editing status within nad1 intron 1 and nad7 intron 4 was seen to differ among plant species, so that homologous intron sequences shared lower similarity at the RNA level than at the DNA level. Plant-specific variation was also seen in the length of the linker joining domains 5 and 6 of nad7 intron 3; it ranged from 4 nt in wheat to 11 nt in soybean, in contrast to the 2-4 nt length typical of classical group II introns. However, this intron is excised as a lariat structure with a domain 6 branchpoint adenosine. Our observations suggest that the core structures and sequences of these plant mitochondrial introns are subject to less stringent evolutionary constraints than conventional group II introns.     

Introns in chloroplast protein-coding genes of land plants

Several protein-coding genes from land plant chloroplasts have been shown to contain introns. The majority of these introns resemble the fungal mitochondrial group II introns due to considerable nucleotide sequence homology at their 5 and 3 ends and they can readily be folded to form six hairpins characteristic of the predicted secondary structure of the mitochondrial group II introns. Recently it has been demonstrated that some mitochondrial group II introns are capable of self-splicing in vitro in the absence of protein co-factors. However evidence presented in this overview suggests that this is probably not the case for chloroplast introns and that trans-acting factors are almost certainly involved in their processing reactions.Abbreviations kop kilobase pairs - ORF Open Reading Frame - pre-RNA precursor ribonucleic acid     

A CRM domain protein functions dually in group I and group II intron splicing in land plant chloroplasts

The CRM domain is a recently recognized RNA binding domain found in three group II intron splicing factors in chloroplasts, in a bacterial protein that associates with ribosome precursors, and in a family of uncharacterized proteins in plants. To elucidate the functional repertoire of proteins with CRM domains, we studied CFM2 (for CRM Family Member 2), which harbors four CRM domains. RNA coimmunoprecipitation assays showed that CFM2 in maize (Zea mays) chloroplasts is associated with the group I intron in pre-trnL-UAA and group II introns in the ndhA and ycf3 pre-mRNAs. T-DNA insertions in the Arabidopsis thaliana ortholog condition a defective-seed phenotype (strong allele) or chlorophyll-deficient seedlings with impaired splicing of the trnL group I intron and the ndhA, ycf3-int1, and clpP-int2 group II introns (weak alleles). CFM2 and two previously described CRM proteins are bound simultaneously to the ndhA and ycf3-int1 introns and act in a nonredundant fashion to promote their splicing. With these findings, CRM domain proteins are implicated in the activities of three classes of catalytic RNA: group I introns, group II introns, and 23S rRNA.     

Bacteria and Archaea Group II introns: additional mobile genetic elements in the environment

Self-splicing group II introns are present in the organelles of lower eukaryotes, plants and Bacteria and have been found recently in Archaea. It is generally accepted that group II introns originated in bacteria before spreading to mitochondria and chloroplasts. These introns are thought to be related to the progenitors of spliceosomal introns. Group II introns are also mobile genetic elements. In bacteria, they appear to spread using either other mobile genetic elements or low-expression regions as target sites. Bacteria and Archaea genome sequence annotations have revealed the diversity of group II intron classes and that they are involved in vertical and horizontal inheritance.     

The Euglena gracilis chloroplast ribulose-1,5-bisphosphate carboxylase gene. I. Complete DNA sequence and analysis of the nine intervening sequences

The nucleotide sequence of 6225 base pairs (bp) of Euglena gracilis chloroplast DNA including the complete DNA sequence of the chloroplast-encoded ribulose-1,5-bisphosphate carboxylase large subunit gene along with the flanking DNA sequences is presented. The gene is greater than 5.5 kilobase pairs in length and is organized as 10 exons coding for 475 amino acids, separated by 9 introns. The exons range in size from 45 to 438 bp, while the introns range in size from 382 to 568 bp. The introns have highly conserved boundary sequences with the consensus, 5'-N GTGTGGATTT...(intron)...TTAATTTTAT N-3'. The introns are 82-85 mol% AT, with a pronounced T greater than A greater than G greater than C base bias in the RNA-like strand. They do not appear to encode any polypeptides. In addition, the introns have a conserved sequence 30-50 bp from their 3'-ends with the consensus, 5'-TACAGTTTGAAAATGA-3'. The 5'-TACA sequence bears some homology to the 5'-end of the TACTAACA sequence found in a similar location in yeast nuclear mRNA introns. The conserved sequences of the Euglena rbcL introns may be indicative of a splicing mechanism similar to that of eucaryotic nuclear mRNA introns and group II mitochondrial introns.     

A purine-pyrimidine motif verifying an identical presence in almost all gene taxonomic groups

A statistical parameter identifies, with a high degree of significance, a motif which is present in protein-coding sequences of eukaryotes, prokaryotes, chloroplasts, mitochondria, viral introns, ribosomal RNA genes, and transfer RNA genes. The random probability of occurrence of such a situation is 10(-12). This motif has the following properties: (i) its significant presence in almost all present-day genes explains why it can be considered as primitive oligonucleotide, (ii) its nucleotide order is: YRY (N)6YRY, R being a purine base, Y a pyrimidine one and N any base, (iii) its length and its terminal trinucleotides YRY suggest a primordial function related to the spatial structure of the DNA sequences. This motif is found in some viral protein-coding genes, but not in eukaryotic introns.     

A contact-waiting-time metric and RNA folding rates

Metrics for indirectly predicting the folding rates of RNA sequences are of interest. In this letter, we introduce a simple metric of RNA structural complexity, which accounts for differences in the energetic contributions of RNA base contacts toward RNA structure formation. We apply the metric to RNA sequences whose folding rates were previously determined experimentally. We find that the metric has good correlation (correlation coefficient: −0.95, p?0.01) with the logarithmically transformed folding rates of those RNA sequences. This suggests that the metric can be useful for predicting RNA folding rates. We use the metric to predict the folding rates of bacterial and eukaryotic group II introns. Future applications of the metric (e.g., to predict structural RNAs) could prove fruitful.     

Group II introns as phylogenetic tools: structure, function, and evolutionary constraints

Group II introns comprise the majority of noncoding DNA in many plant chloroplast genomes and include the commonly sequenced regions trnK/matK, the rps16 intron, and the rpl16 intron. As demand increases for nucleotide characters at lower taxonomic levels, chloroplast introns may come to provide the bulk of plastome sequence data for assessment of evolutionary relationships in infrageneric, intergeneric, and interfamilial studies. Group II introns have many attractive properties for the molecular systematist: they are confined to organellar genomes in eukaryotes and the majority are single-copy; they share a well-defined and empirically tested secondary and tertiary structure; and many are easily amplified due to highly conserved sequence in flanking exons. However, structure-linked mutation patterns in group II intron sequences are more complex than generally supposed and have important implications for aligning nucleotides, assessing mutational biases in the data, and selecting appropriate models of character evolution for phylogenetic analysis. This paper presents a summary of group II intron function and structure, reviews the link between that structure and specific mutational constraints in group II intron sequences, and discusses strategies for accommodating the resulting complex mutational patterns in subsequent phylogenetic analyses.     

Survey of group I and group II introns in 29 sequenced genomes of the Bacillus cereus group: insights into their spread and evolution

Group I and group II introns are different catalytic self-splicing and mobile RNA elements that contribute to genome dynamics. In this study, we have analyzed their distribution and evolution in 29 sequenced genomes from the Bacillus cereus group of bacteria. Introns were of different structural classes and evolutionary origins, and a large number of nearly identical elements are shared between multiple strains of different sources, suggesting recent lateral transfers and/or that introns are under a strong selection pressure. Altogether, 73 group I introns were identified, inserted in essential genes from the chromosome or newly described prophages, including the first elements found within phages in bacterial plasmids. Notably, bacteriophages are an important source for spreading group I introns between strains. Furthermore, 77 group II introns were found within a diverse set of chromosomal and plasmidic genes. Unusual findings include elements located within conserved DNA metabolism and repair genes and one intron inserted within a novel retroelement. Group II introns are mainly disseminated via plasmids and can subsequently invade the host genome, in particular by coupling mobility with host cell replication. This study reveals a very high diversity and variability of mobile introns in B. cereus group strains.