The chloroplast genome sequence of Chara vulgaris sheds new light into the closest green algal relatives of land plants

The phylum Streptophyta comprises all land plants and six monophyletic groups of charophycean green algae (Mesostigmatales, Chlorokybales, Klebsormidiales, Zygnematales, Coleochaetales, and Charales). Phylogenetic analyses of four genes encoded in three cellular compartments suggest that the Charales are sister to land plants and that charophycean green algae evolved progressively toward an increasing cellular complexity. To validate this phylogenetic hypothesis and to understand how and when the highly conservative pattern displayed by land plant chloroplast DNAs (cpDNAs) originated in the Streptophyta, we have determined the complete chloroplast genome sequence (184,933 bp) of a representative of the Charales, Chara vulgaris, and compared this genome to those of Mesostigma (Mesostigmatales), Chlorokybus (Chlorokybales), Staurastrum and Zygnema (Zygnematales), Chaetosphaeridium (Coleochaetales), and selected land plants. The phylogenies we inferred from 76 cpDNA-encoded proteins and genes using various methods favor the hypothesis that the Charales diverged before the Coleochaetales and Zygnematales. The Zygnematales were identified as sister to land plants in the best tree topology (T1), whereas Chaetosphaeridium (T2) or a clade uniting the Zygnematales and Chaetosphaeridium (T3) occupied this position in alternative topologies. Chara remained at the same basal position in trees including more land plant taxa and inferred from 56 proteins/genes. Phylogenetic inference from gene order data yielded two most parsimonious trees displaying the T1 and T3 topologies. Analyses of additional structural cpDNA features (gene order, gene content, intron content, and indels in coding regions) provided better support for T1 than for the topology of the above-mentioned four-gene tree. Our structural analyses also revealed that many of the features conserved in land plant cpDNAs were inherited from their green algal ancestors. The intron content data predicted that at least 15 of the 21 land plant group II introns were gained early during the evolution of streptophytes and that a single intron was acquired during the transition from charophycean green algae to land plants. Analyses of genome rearrangements based on inversions predicted no alteration in gene order during the transition from charophycean green algae to land plants.     

Polymorphism and structural variation of rps16 group-II intron in the Solanum species

Nucleotide sequences of the self-splicing group-II intron of rps16 have first been determined in nine species of the Solanum genus. It was found that the observed variations in the intron length (855–864 bp) was associated with indels of 1 to 9 bp. Altogether, five indels and 50 nucleotide substitutions were detected, which were used to identify six Solanum haplotypes. Although the intron sequence was in general fairly well conserved, the distribution of the described mutations among its structural elements corresponding to six pre-RNA domains was qualitatively and quantitatively nonuniform. The highest polymorphism levels were observed in domains I, II, and IV. The sequence of domain V was absolutely invariable, which is in agreement with its functional significance. The chloroplast rpS16 intron sequences have been characterized in nine Solanum species. The intron length ranged from 855 bp to 864 bp, which is associated with 1–9-nucleotide indels. In total five indels and 50 nucleotide substitutions have been detected and six Solanum haplotypes have been revealed. Solanum rpS16 introns has been characterized by mutation rate heterogeneity between structure regions of all six domains its pre-RNA. Intron domains I, II, IV are shown to be more variable. Sequences of the domain V are invariant, that agrees with its functional significance.     

Structural characteristics of the chloroplast <Emphasis Type="Italic">rpS16</Emphasis> intron in <Emphasis Type="Italic">Allium sativum</Emphasis> and related <Emphasis Type="Italic">Allium</Emphasis> species

The intron sequence of chloroplast rpS16 and the secondary structure of its pre-mRNA were characterized for the first time in 26 Allium sativum accessions of different ecologo-geographical origins and seven related Allium species. The boundaries and main stem-loop consensus sequences were identified for all six domains of the intron. Polymorphism was estimated for the total intron and its regions. The structural regions of the rpS16 intron proved to be heterogeneous for mutation rate and spectrum. Mutations were most abundant in domains II and IV, and transition predominated in domains I, III, V, and VI. In addition to structural elements and motifs typical for group IIB introns, several Allium-specific micro- and macrostructural mutations were revealed. A 290-bp deletion involving domains III and IV and part of domain V was observed in A. altaicum, A. fistulosum, and A. schoenoprasum. Several indels and nucleotide substitutions were found to cause a deviation of the pre-mRNA secondary structure from the consensus model of group II introns.     

Molecular systematics of Apiaceae subfamily Apioideae: phylogenetic analyses of nuclear ribosomal DNA internal transcribed spacer and plastid RPO C1 intron sequences

Evolutionary relationships among representatives of Apiaceae (Umbelliferae) subfamily Apioideae have been inferred from phylogenetic analyses of nuclear ribosomal DNA internal transcribed spacer (ITS 1 and ITS 2) and plastid rpoC1 intron sequences. High levels of nucleotide sequence variation preclude the use of the ITS region for examining relationships across subfamilial boundaries in Apiaceae, whereas the rpoC1 intron is more suitably conserved for family-wide phylogenetic study but is too conserved for examining relationships among closely related taxa. In total, 126 ITS sequences from subfamily Apioideae and 100 rpoC1 intron sequences from Apiaceae (all three subfamilies) and outgroups Araliaceae and Pittosporaceae were examined. Phylogenies estimated using parsimony, neighbor-joining, and maximum likelihood methods reveal that: (1) Apiaceae subfamily Apioideae is monophyletic and is sister group to Apiaceae subfamily Saniculoideae; (2) Apiaceae subfamily Hydrocotyloideae is not monophyletic, with some members strongly allied to Araliaceae and others to Apioideae + Saniculoideae; and (3) Apiaceae subfamily Apioideae comprises several well-supported subclades, but none of these coincide with previously recognized tribal divisions based largely on morphological and anatomical characters of the fruit. Four major clades in Apioideae are provisionally recognized and provide the framework for future lower level phylogenetic analyses. A putative secondary structure model of the Daucus carota (carrot) rpoC1 group II intron is presented. Of its six major structural domains, domains II and III are the most, and domains V and VI the least, variable.     

Bindin from a sea star

SUMMARY The genetic basis for the evolution of development includes genes that encode proteins expressed on the surfaces of sperm and eggs. Previous studies of the sperm acrosomal protein bindin have helped to characterize the adaptive evolution of gamete compatibility and speciation in sea urchins. The absence of evidence for bindin expression in taxa other than the Echinoidea has limited such studies to sea urchins, and led to the suggestion that bindin might be a sea urchin-specific molecule. Here we characterize the gene that encodes bindin in a broadcast-spawning asterinid sea star ( Patiria miniata ). We describe the sequence and domain structure of a full-length bindin cDNA and its single intron. In comparison with sea urchins, P. miniata bindin is larger but the two molecules share several general features of their domain structure and some sequence features of two domains. Our results extend the known evolutionary history of bindin from the Mesozoic (among the crown group sea urchins) into the early Paleozoic (and the common ancestor of eleutherozoans), and present new opportunities for understanding the role of bindin molecular evolution in sexual selection, life history evolution, and speciation among sea stars.     

Evolution of the mitochondrial rps3 intron in perennial and annual angiosperms and homology to nad5 intron 1.

The plant mitochondrial rps3 intron was analyzed for substitution and indel rate variation among 15 monocot and dicot angiosperms from 10 genera, including perennial and annual taxa. Overall, the intron sequence was very conserved among angiosperms. Based on length polymorphism, 10 different alleles were identified among the 10 genera. These allelic differences were mainly attributable to large indels. An insertion of 133 nucleotides, observed in the Alnus intron was partially or completely absent in the other lineages of the family Betulaceae. This insertion was located within domain IV of the secondary-structure model of this group IIA intron. A mobile element of 47 nucleotides that showed homology to sequences located in rice rps3 intron and in intergenic plant mitochondrial genomes was found within this insertion. Both substitution and indel rates were low among the Betulaceae sequences, but substitution rates were increasingly larger than indel rates in comparisons involving more distantly related taxa. From a secondary-structure model, regions involved in helical structures were shown to be well preserved from indels as compared to substitutions, but compensatory changes were not observed among the angiosperm sequences analyzed. Using approximate divergence times based on the fossil record, substitution and indel rate heterogeneity was observed between different pairs of annual and perennial taxa. In particular, the annual petunia and primrose evolved more than 15 and 10 times faster, for substitution and indel rates respectively, than the perennial birch and alder. This is the first demonstration of an evolutionary rate difference between perennial and annual forms in noncoding DNA, lending support to neutral causes such as the generation time, population size, and speciation rate effects to explain such rate heterogeneity. Surprisingly, the sequence from the rps3 intron had a high identity with the sequence of intron 1 from the angiosperm mitochondrial nad5 gene, suggesting a common origin of these two group IIA introns.     

Domains 2 and 3 interact to form critical elements of the group II intron active site

Group II introns are self-splicing RNA molecules that also behave as mobile genetic elements. The secondary structure of group II intron RNAs is typically described as a series of six domains that project from a central wheel. Most structural and mechanistic analyses of the intron have focused on domains 1 and 5, which contain the residues essential for catalysis, and on domain 6, which contains the branch-point adenosine. Domains 2 and 3 (D2, D3) have been shown to make important contributions to intronic activity; however, information about their function is quite limited. To elucidate the role of D2 and D3 in group II ribozyme catalysis, we built a series of multi-piece ribozyme constructs based on the ai5gamma group II intron. These constructs are designed to shed light on the roles of D2 and D3 in some of the major reactions catalyzed by the intron: 5'-exon cleavage, branching, and substrate hydrolysis. Reactions with these constructs demonstrate that D3 stimulates the chemical rate constant of group II intron reactions, and that it behaves as a form of catalytic effector. However, D3 is unable to associate independently with the ribozyme core. Docking of D3 is mediated by a short duplex that is found at the base of D2. In addition to recruiting D3 into the core, the D2 stem directs the folding of the adjacent j(2/3) linker, which is among the most conserved elements in the group II intron active site. In turn, the D2 stem contributes to 5'-splice site docking and ribozyme conformational change. Nucleotide analog interference mapping suggests an interaction between the D2 stem and D3 that builds on the known theta-theta' interaction and extends it into D3. These results establish that D3 and the base of D2 are key elements of the group II intron core and they suggest a hierarchy for active-site assembly.     

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.     

Group II introns in wheat mitochondria have degenerate structural features and varied splicing pathways

Mitochondrial introns in flowering plant genes are virtually all classified as members of the group II ribozyme family although certain structural features have degenerated to varying degrees over evolutionary time. We are interested in the impact that unconventional intron architecture might have on splicing biochemistry in vivo and we have focused in particular on intronic domains V and VI, which for self-splicing introns provide a key component of the catalytic core and the bulged branchpoint adenosine, respectively. Notably, the two transesterification steps in classical group II splicing are the same as for nuclear spliceosomal introns and release the intron as a lariat. Using RT-PCR and circularized RT-PCR, we had previously demonstrated that several wheat mitochondrial introns which lack a branchpoint adenosine have atypical splicing pathways, and we have now extended this analysis to the full set of wheat introns, namely six trans-splicing and sixteen cis-splicing ones. A number of introns are excised using non-lariat pathways and interestingly, we find that several introns which do have a conventional domain VI also use pathways that appear to exploit other internal or external nucleophiles, with the lariat form being relatively minor. Somewhat surprisingly, several introns with weakly-structured domain V/VI helices still exhibit classical lariat splicing, suggesting that accessory factors aid in restoring a splicing-competent conformation. Our observations illustrate that the loss of conventional group II features during evolution is correlated with altered splicing biochemistry in an intron-distinctive manner.     

Molecular evolution of nitrate reductase genes

To understand the evolutionary mechanisms and relationships of nitrate reductases (NRs), the nucleotide sequences encoding 19 nitrate reductase (NR) genes from 16 species of fungi, algae, and higher plants were analyzed. The NR genes examined show substantial sequence similarity, particularly within functional domains, and large variations in GC content at the third codon position and intron number. The intron positions were different between the fungi and plants, but conserved within these groups. The overall and nonsynonymous substitution rates among fungi, algae, and higher plants were estimated to be 4.33 × 10⁻¹⁰ and 3.29 × 10⁻¹⁰ substitutions per site per year. The three functional domains of NR genes evolved at about one-third of the rate of the N-terminal and the two hinge regions connecting the functional domains. Relative rate tests suggested that the nonsynonymous substitution rates were constant among different lineages, while the overall nucleotide substitution rates varied between some lineages. The phylogenetic trees based on NR genes correspond well with the phylogeny of the organisms determined from systematics and other molecular studies. Based on the nonsynonymous substitution rate, the divergence time of monocots and dicots was estimated to be about 340 Myr when the fungi–plant or algae–higher plant divergence times were used as reference points and 191 Myr when the rice–barley divergence time was used as a reference point. These two estimates are consistent with other estimates of divergence times based on these reference points. The lack of consistency between these two values appears to be due to the uncertainty of the reference times. Received: 10 April 1995 / Accepted: 10 September 1995     

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.     

Intron evolution in a phylogenetic perspective: Divergent trends in the two copies of the duplicated def gene in Impatiens L. (Balsaminaceae)

The history of MADS box genes is well-known in angiosperms. While duplication events and gene losses occur frequently, gene structure and intron positions are very conserved. We investigated all six introns in a duplicated MADS box gene (deficiens, def) in selected Impatiens taxa, thereby assessing intron features. For the first time, our study provides a comparison of molecular changes in all introns of a gene from a phylogenetic perspective. Interestingly, a uniform pattern of molecular evolution in the introns of each copy was not observed, but intron length increases, decreases, and size retention can be found in each copy. A tendency to accumulate long autapomorphic indels is also present, thus, a longer intron length does not reflect a higher number of parsimony-informative characters. Substitution rates vary between introns of each gene copy. While four of the six introns of def1 exhibit a change in their substitution rate, five of the six def2 introns maintain their rates throughout the genus albeit at different levels. In MADS box genes several regulatory sequences are found residing in introns. Thus, presence of putative regulatory motifs was investigated. Most of them are not conserved in position and usually present in only one of the gene copies. In addition, the potential for phylogenetic reconstruction of introns in both def copies is shortly discussed.     

Group II introns deleted for multiple substructures retain self-splicing activity.

Group II introns can be folded into highly conserved secondary structures with six major substructures or domains. Domains 1 and 5 are known to play key roles in self-splicing, while the roles of domains 2, 3, 4, and 6 are less clear. A trans assay for domain 5 function has been developed which indicates that domain 5 has a binding site on the precursor RNA that is not predicted from any secondary structure element. In this study, the self-splicing group II intron 5 gamma of the coxI gene of yeast mitochondrial DNA was deleted for various intron domains, singly and in combinations. Those mutant introns were characterized for self-splicing reactions in vitro as a means of locating the domain 5 binding site. A single deletion of domain 2, 3, 4, or 6 does not block in vitro reactions at either splice junction, though the deletion of domain 6 reduces the fidelity of 3' splice site selection somewhat. Even the triple deletion lacking domains 2, 4, and 6 retains some self-splicing activity. The deletion of domains 2, 3, 4, and 6 blocks the reaction at the 3' splice junction but not at the 5' junction. From these results, we conclude that the binding site for domain 5 is within domain 1 and that the complex of 5' exon, domain 1, and domain 5 (plus short connecting sequences) constitutes the essential catalytic core of this intron.     

Extensive structural conservation exists among several homologs of two Euglena chloroplast group II introns

82 of the 155 chloroplast introns in Euglena gracilis have been categorized as group II introns. Because they are shorter and more divergent than group II introns from other organisms, the assignment of these Euglena introns to the group II class has been questioned. In the current study, two homologs of E. gracilispetB intron 1 and four homologs of psbC intron 2 have been isolated from related species and characterized. Based on a comparative sequence analysis of intron homologs, the intron core and four of the six helical domains present in the canonical group II intron structural model are conserved in E. gracilispetB intron 1 and psbC intron 2 and all of their homologs. Distal portions of domain I, which are involved in most of the tertiary interactions, are less well conserved than the central core.     

The complete nucleotide sequence of the cassava (Manihot esculenta) chloroplast genome and the evolution of atpF in Malpighiales: RNA editing and multiple losses of a group II intron

The complete sequence of the chloroplast genome of cassava (Manihot esculenta, Euphorbiaceae) has been determined. The genome is 161,453 bp in length and includes a pair of inverted repeats (IR) of 26,954 bp. The genome includes 128 genes; 96 are single copy and 16 are duplicated in the IR. There are four rRNA genes and 30 distinct tRNAs, seven of which are duplicated in the IR. The infA gene is absent; expansion of IRb has duplicated 62 amino acids at the 3′ end of rps19 and a number of coding regions have large insertions or deletions, including insertions within the 23S rRNA gene. There are 17 intron-containing genes in cassava, 15 of which have a single intron while two (clpP, ycf3) have two introns. The usually conserved atpF group II intron is absent and this is the first report of its loss from land plant chloroplast genomes. The phylogenetic distribution of the atpF intron loss was determined by a PCR survey of 251 taxa representing 34 families of Malpighiales and 16 taxa from closely related rosids. The atpF intron is not only missing in cassava but also from closely related Euphorbiaceae and other Malpighiales, suggesting that there have been at least seven independent losses. In cassava and all other sequenced Malphigiales, atpF gene sequences showed a strong association between C-to-T substitutions at nucleotide position 92 and the loss of the intron, suggesting that recombination between an edited mRNA and the atpF gene may be a possible mechanism for the intron loss.