Multiple major increases and decreases in mitochondrial substitution rates in the plant family Geraniaceae

Background  

Rates of synonymous nucleotide substitutions are, in general, exceptionally low in plant mitochondrial genomes, several times lower than in chloroplast genomes, 10–20 times lower than in plant nuclear genomes, and 50–100 times lower than in many animal mitochondrial genomes. Several cases of moderate variation in mitochondrial substitution rates have been reported in plants, but these mostly involve correlated changes in chloroplast and/or nuclear substitution rates and are therefore thought to reflect whole-organism forces rather than ones impinging directly on the mitochondrial mutation rate. Only a single case of extensive, mitochondrial-specific rate changes has been described, in the angiosperm genus Plantago.     

Evolutionary rates of commonly used nuclear and organelle markers of Arabidopsis relatives (Brassicaceae)

Recovering the genetic divergence between species is one of the major interests in the evolutionary biology. It requires accurate estimation of the neutral substitution rates. Arabidopsis thaliana, the first whole-genome sequenced plant, and its out-crossing relatives provide an ideal model for examining the split between sister species. In the study, rates of molecular evolution at markers frequently used for systematics and population genetics, including 14 nuclear genes spanning most chromosomes, three noncoding regions of chloroplast genome, and one intron of mitochondrial genome, between A. thaliana and four relatives were estimated. No deviation from neutrality was detected in the genes examined. Based on the known divergence between A. thaliana and its sisters about 8.0-17.6 MYA, evolutionary rates of the eighteen genes were estimated. Accordingly, the ratio of rates of synonymous substitutions among mitochondrial, chloroplast and nuclear genes was calculated with an average and 95% confidence interval of 1 (0.25-1.75): 15.77 (7.48-114.09): 74.79 (36.27-534.61). Molecular evolutionary rates of nuclear genes varied, with a range of 0.383-0.856×10(-8) for synonymous substitutions per site per year and 0.036-0.081×10(-9) for nonsynonymous substitutions per site per year. Compared with orthologs in Populus, a long life-span tree, genes in Arabidopsis evolved faster in an order of magnitude at the gene level, agreeing with a generation time hypothesis. The estimated substitution rates of these genes can be used as a reference for molecular dating.     

Synonymous substitution rates in Drosophila: Mitochondrial versus nuclear genes

Synonymous substitution rates in mitochondrial and nuclear genes of Drosophila were compared. To make accurate comparisons, we considered the following: (1) relative synonymous rates, which do not require divergence time estimates, should be used; (2) methods estimating divergence should take into account base composition; (3) only very closely related species should be used to avoid effects of saturation; (4) the heterogeneity of rates should be examined. We modified the methods estimating synonymous substitution numbers to account for base composition bias. By using these methods, we found that mitochondrial genes have 1.7–3.4 times higher synonymous substitution rates than the fastest nuclear genes or 4.5–9.0 times higher rates than the average nuclear genes. The average rate of synonymous transversions was 2.7 (estimated from the melanogaster species subgroup) or 2.9 (estimated from the obscura group) times higher in mitochondrial genes than in nuclear genes. Synonymous transversions in mitochondrial genes occurred at an approximately equivalent rate to those in the fastest nuclear genes. This last result is not consistent with the hypothesis that the difference in turnover rates between mitochondrial and nuclear genomes is the major factor determining higher synonymous substitution rates in mtDNA. We conclude that the difference in synonymous substitution rates is due to a combination of two factors: a higher transitional mutation rate in mtDNA and constraints on nuclear genes due to selection for codon usage. Received: 27 November 1996 / Accepted: 8 May 1997     

Correlated rates of synonymous site evolution across plant genomes

Synonymous substitution rates have been shown to vary among evolutionary lineages of both nuclear and organellar genes across a broad range of taxonomic groups. In animals, rate heterogeneity does not appear to be correlated across nuclear and mitochondrial genes. In this paper, we contrast substitution rates in two plant groups and show that grasses evolve more rapidly than palms at synonymous sites in a mitochondrial, a nuclear, and a plastid gene. Furthermore, we show that the relative rates of synonymous substitution between grasses and palms are similar at the three loci. The correlation in synonymous substitution rates across genes is particularly striking because the three genes evolve at very different absolute rates. In contrast, relative rates of nonsynonymous substitution are not conserved among the three genes.      

Noncoding sequences from the slowly evolving chloroplast inverted repeat in addition to rbcL data do not support gnetalean affinities of angiosperms

We developed PCR primers against highly conserved regions of the rRNA operon located within the inverted repeat of the chloroplast genome and used these to amplify the region spanning from the 3' terminus of the 23S rRNA gene to the 5' terminus of the 5S rRNA gene. The sequence of this roughly 500-bp region, which includes the 4.5S rRNA gene and two chloroplast intergenic transcribed spacer regions (cpITS2 and cpITS3), was determined from 20 angiosperms, 7 gymnosperms, and 16 ferns (21,700 bp). Sequences for the large subunit of ribulose bisphosphate carboxylase/oxygenase (rbcL) from the same or confamilial genera were analyzed in both separate and combined data sets. Due to the low substitution rate in the inverted repeat region, noncoding sequences in the cpITS region are not saturated with substitutions, in contrast to synonymous sites in rbcL, which are shown to evolve roughly six times faster than noncoding cpITS sequences. Several length polymorphisms with very clear phylogenetic distributions were detected in the data set. Results of phylogenetic analyses provide very strong bootstrap support for monophyly of both spermatophytes and angiosperms. No support for a sister group relationship between Gnetales and angiosperms in either cpITS or rbcL data was found. Rather, weak bootstrap support for monophyly of gymnosperms studied and for a basal position for the aquatic angiosperm Nymphaea among angiosperms studied was observed. Noncoding sequences from the inverted repeat region of chloroplast DNA appear suitable for study of land plant evolution.      

Dramatically elevated rate of mitochondrial substitution in lice (Insecta: Phthiraptera)

Few estimates of relative substitution rates, and the underlying mutation rates, exist between mitochondrial and nuclear genes in insects. Previous estimates for insects indicate a 2-9 times faster substitution rate in mitochondrial genes relative to nuclear genes. Here we use novel methods for estimating relative rates of substitution, which incorporate multiple substitutions, and apply these methods to a group of insects (lice, Order: Phthiraptera). First, we use a modification of copath analysis (branch length regression) to construct independent comparisons of rates, consisting of each branch in a phylogenetic tree. The branch length comparisons use maximum likelihood models to correct for multiple substitution. In addition, we estimate codon-specific rates under maximum likelihood for the different genes and compare these values. Estimates of the relative synonymous substitution rates between a mitochondrial (COI) and nuclear (EF-1alpha) gene in lice indicate a relative rate of several 100 to 1. This rapid relative mitochondrial rate (>100 times) is at least an order of magnitude faster than previous estimates for any group of organisms. Comparisons using the same methods for another group of insects (aphids) reveals that this extreme relative rate estimate is not simply attributable to the methods we used, because estimates from aphids are substantially lower. Taxon sampling affects the relative rate estimate, with comparisons involving more closely related taxa resulting in a higher estimate. Relative rate estimates also increase with model complexity, indicating that methods accounting for more multiple substitution estimate higher relative rates.     

Genes for two mitochondrial ribosomal proteins in flowering plants are derived from their chloroplast or cytosolic counterparts

Often during flowering plant evolution, ribosomal protein genes have been lost from the mitochondrion and transferred to the nucleus. Here, we show that substitution by a duplicated, divergent gene originally encoding the chloroplast or cytosolic ribosomal protein counterpart accounts for two missing mitochondrial genes in diverse angiosperms. The rps13 gene is missing from the mitochondrial genome of many rosids, and a transferred copy of this gene is not evident in the nucleus of Arabidopsis, soybean, or cotton. Instead, these rosids contain a divergent nuclear copy of an rps13 gene of chloroplast origin. The product of this gene from all three rosids was shown to be imported into isolated mitochondria but not into chloroplasts. The rps8 gene is missing from the mitochondrion and nucleus of all angiosperms examined. A divergent copy of the gene encoding its cytosolic counterpart (rps15A) was identified in the nucleus of four angiosperms and one gymnosperm. The product of this gene from Arabidopsis and tomato was imported successfully into mitochondria. We infer that rps13 was lost from the mitochondrial genome and substituted by a duplicated nuclear gene of chloroplast origin early in rosid evolution, whereas rps8 loss and substitution by a gene of nuclear/cytosolic origin occurred much earlier, in a common ancestor of angiosperms and gymnosperms.     

Angiosperm phylogeny inferred from sequences of four mitochondrial genes

An angiosperm phylogeny was reconstructed in a maximum likelihood analysis of sequences of four mitochondrial genes, atpl, matR, had5, and rps3, from 380 species that represent 376 genera and 296 families of seed plants. It is largely congruent with the phylogeny of angiosperms reconstructed from chloroplast genes atpB, matK, and rbcL, and nuclear 18S rDNA. The basalmost lineage consists of Amborella and Nymphaeales (including Hydatellaceae). Austrobaileyales follow this clade and are sister to the mesangiosperms, which include Chloranthaceae, Ceratophyllum, magnoliids, monocots, and eudicots. With the exception of Chloranthaceae being sister to Ceratophyllum, relationships among these five lineages are not well supported. In eudicots, Ranunculales, Sabiales, Proteales, Trochodendrales, Buxales, Gunnerales, Saxifragales, Vitales, Berberidopsidales, and Dilleniales form a basal grade of lines that diverged before the diversification of rosids and asterids. Within rosids, the COM (Celastrales-Oxalidales-Malpighiales) clade is sister to malvids (or rosid Ⅱ), instead of to the nitrogen-fixing clade as found in all previous large-scale molecular analyses of angiosperms. Santalales and Caryophyllales are members of an expanded asterid clade. This study shows that the mitochondrial genes are informative markers for resolving relationships among genera, families, or higher rank taxa across angiosperms. The low substitution rates and low homoplasy levels of the mitochondrial genes relative to the chloroplast genes, as found in this study, make them particularly useful for reconstructing ancient phylogenetic relationships. A mitochondrial gene-based angiosperm phylogeny provides an independent and essential reference for comparison with hypotheses of angiosperm phylogeny based on chloroplast genes, nuclear genes, and non-molecular data to reconstruct the underlying organismal phylogeny.     

Rapidly evolving mitochondrial genome and directional selection in mitochondrial genes in the parasitic wasp nasonia (hymenoptera: pteromalidae)

We sequenced the nearly complete mtDNA of 3 species of parasitic wasps, Nasonia vitripennis (2 strains), Nasonia giraulti, and Nasonia longicornis, including all 13 protein-coding genes and the 2 rRNAs, and found unusual patterns of mitochondrial evolution. The Nasonia mtDNA has a unique gene order compared with other insect mtDNAs due to multiple rearrangements. The mtDNAs of these wasps also show nucleotide substitution rates over 30 times faster than nuclear protein-coding genes, indicating among the highest substitution rates found in animal mitochondria (normally <10 times faster). A McDonald and Kreitman test shows that the between-species frequency of fixed replacement sites relative to silent sites is significantly higher compared with within-species polymorphisms in 2 mitochondrial genes of Nasonia, atp6 and atp8, indicating directional selection. Consistent with this interpretation, the Ka/Ks (nonsynonymous/synonymous substitution rates) ratios are higher between species than within species. In contrast, cox1 shows a signature of purifying selection for amino acid sequence conservation, although rates of amino acid substitutions are still higher than for comparable insects. The mitochondrial-encoded polypeptides atp6 and atp8 both occur in F0F1ATP synthase of the electron transport chain. Because malfunction in this fundamental protein severely affects fitness, we suggest that the accelerated accumulation of replacements is due to beneficial mutations necessary to compensate mild-deleterious mutations fixed by random genetic drift or Wolbachia sweeps in the fast evolving mitochondria of Nasonia. We further propose that relatively high rates of amino acid substitution in some mitochondrial genes can be driven by a "Compensation-Draft Feedback"; increased fixation of mildly deleterious mutations results in selection for compensatory mutations, which lead to fixation of additional deleterious mutations in nonrecombining mitochondrial genomes, thus accelerating the process of amino acid substitutions.     

Slow molecular evolution in 18S rDNA, rbcL and nad5 genes of mosses compared with higher plants

The evolutionary potential of bryophytes (mosses, liverworts and hornworts) has been debated for decades. Fossil record and biogeographical distribution patterns suggest very slow morphological evolution and the retainment of several ancient traits since the split with vascular plants some 450 million years ago. Many have argued that bryophytes may evolve as rapidly as higher plants on the molecular level, but this hypothesis has not been tested so far. Here, it is shown that mosses have experienced significantly lower rates of molecular evolution than higher plants within 18S rDNA (nuclear), rbcL (chloroplast) and nad5 (mitochondrial) genes. Mosses are on an average evolving 2-3 times slower than ferns, gymnosperms and angiosperms; and also green algae seem to be evolving faster than nonvascular plants. These results support the observation of a general correlation between morphological and molecular evolutionary rates in plants and also show that mosses are 'evolutionary sphinxes' regarding both morphological and molecular evolutionary potential.     

Transfer of chloroplast genomic DNA to mitochondrial genome occurred at least 300 MYA

With the completion of the first gymnosperm mitochondrial genome (mtDNA) from Cycas taitungensis and the availability of more mtDNA taxa in the past 5 years, we have conducted a systematic analysis of DNA transfer from chloroplast genomes (cpDNAs) to mtDNAs (mtpts) in 11 plants, including 2 algae, 1 liverwort, 1 moss, 1 gymnosperm, 3 monocots, and 3 eudicots. By using shared gene order and boundaries between different mtpts as the criterion, the timing of cpDNA transfer during plant evolution was estimated from the phylogenetic tree reconstructed independently from concatenated protein-coding genes of 11 available mtDNAs. Several interesting findings emerged. First, frequent DNA transfer from cpDNA to mtDNA occurred at least as far back as the common ancestor of extant gymnosperms and angiosperms, about 300 MYA. The oldest mtpt is trnV(uac)-trnM(cau)-atpE-atpB-rbcL. Three other mtpts--psaA-psaB, rps19-trnH(gug)-rpl2-rpl23, and psbE-psbF--were dated to the common ancestor of extant angiosperms, at least 150 MYA. However, all protein-coding genes of mtpts have degenerated since their first transfer. Therefore, mtpts contribute nothing to the functioning of mtDNA but junk sequences. We discovered that the cpDNA transfers have occurred randomly at any positions of the cpDNAs. We provide strong evidence that the cp-derived tRNA-trnM(cau) is the only mtpt (1 out of 3 cp-derived tRNA shared by seed plants) truly transferred from cpDNA to mtDNA since the time of the common ancestor of extant gymnosperms and angiosperms. Our observations support the proposition of Richly and Leister (2004) that "primary insertions of organellar DNAs are large and then diverge and fragment over evolutionary time."     

南昌市不同植物类群叶片氮磷浓度及其化学计量比

对南昌大学前湖校区89种主要植物叶片的N、P浓度及其化学计量比进行了研究,结果表明:乔灌、常绿、针叶、种子、裸子和单子叶植物类群的N浓度分别低于相对应的草本、落叶、阔叶、蕨类、被子和双子叶植物类群,而C3和C4植物差异不显著;乔灌、常绿和裸子植物类群的P浓度含量分别低于相对应的草本、落叶和被子植物类群,而针叶和阔叶、蕨类和种子、单子叶和双子叶、C3和C4植物类群间差异不显著;乔木、阔叶、被子和双子叶植物类群叶片N/P分别高于相对应的灌草、针叶、裸子和单子叶植物类群,而常绿和落叶、蕨类和种子、C3和C4植物类群之间差异不显著.可见,不同类型植物对N和P的吸收利用存在差异,且对不同养分供应采取不同的适应对策.结合研究区土壤养分现状,建议优先选择常绿、针叶、裸子和单子叶植物类群作为城市园林植物.     

植物三种不同遗传方式基因的地理家系谱理论与应用初探

将已知用于从地理空间上离散或连续分布群体随机抽取基因样本的基因家系谱理论推广到两性异交植物上。由于存在不同的群体间基因多率,地３种不同遗传方式的植物基因组（核、叶绿体和线粒体ＤＮＡ）分别进行了研究。理论上证明对于不同遗传方式的基因,通过相应调整有效群体大小和迁移率,现有的基因家系谱理论可直接应用于植物群体上,其中一个结论就是当从离散分布群体中随机抽取ｎ个基因样本时,亚群体间的花粉流和种子流的相对比     

Chloroplast genome (cpDNA) of Cycas taitungensis and 56 cp protein-coding genes of Gnetum parvifolium: insights into cpDNA evolution and phylogeny of extant seed plants

Phylogenetic relationships among the 5 groups of extant seed plants are presently unsettled. To reexamine this long-standing debate, we determine the complete chloroplast genome (cpDNA) of Cycas taitungensis and 56 protein-coding genes encoded in the cpDNA of Gnetum parvifolium. The cpDNA of Cycas is a circular molecule of 163,403 bp with 2 typical large inverted repeats (IRs) of 25,074 bp each. We inferred phylogenetic relationships among major seed plant lineages using concatenated 56 protein-coding genes in 37 land plants. Phylogenies, generated by the use of 3 independent methods, provide concordant and robust support for the monophylies of extant seed plants, gymnosperms, and angiosperms. Within the modern gymnosperms are 2 highly supported sister clades: Cycas-Ginkgo and Gnetum-Pinus. This result agrees with both the "gnetifer" and "gnepines" hypotheses. The sister relationships in Cycas-Ginkgo and Gnetum-Pinus clades are further reinforced by cpDNA structural evidence. Branch lengths of Cycas-Ginkgo and Gnetum were consistently the shortest and the longest, respectively, in all separate analyses. However, the Gnetum relative rate test revealed this tendency only for the 3rd codon positions and the transversional sites of the first 2 codon positions. A PsitufA located between psbE and petL genes is here first detected in Anthoceros (a hornwort), cycads, and Ginkgo. We demonstrate that the PsitufA is a footprint descended from the chloroplast tufA of green algae. The duplication of ycf2 genes and their shift into IRs should have taken place at least in the common ancestor of seed plants more than 300 MYA, and the tRNAPro-GGG gene was lost from the angiosperm lineage at least 150 MYA. Additionally, from cpDNA structural comparison, we propose an alternative model for the loss of large IR regions in black pine. More cpDNA data from non-Pinaceae conifers are necessary to justify whether the gnetifer or gnepines hypothesis is valid and to generate solid structural evidence for the monophyly of extant gymnosperms.     

Cenozoic extinctions account for the low diversity of extant gymnosperms compared with angiosperms

We test the widely held notion that living gymnosperms are 'ancient' and 'living fossils' by comparing them with their sister group, the angiosperms. This perception derives partly from the lack of gross morphological differences between some Mesozoic gymnosperm fossils and their living relatives (e.g. Ginkgo, cycads and dawn redwood), suggesting that the rate of evolution of gymnosperms has been slow. We estimated the ages and diversification rates of gymnosperm lineages using Bayesian relaxed molecular clock dating calibrated with 21 fossils, based on the phylogenetic analysis of alignments of matK chloroplast DNA (cpDNA) and 26S nuclear ribosomal DNA (nrDNA) sequences, and compared these with published estimates for angiosperms. Gymnosperm crown groups of Cenozoic age are significantly younger than their angiosperm counterparts (median age: 32 Ma vs 50 Ma) and have long unbranched stems, indicating major extinctions in the Cenozoic, in contrast with angiosperms. Surviving gymnosperm genera have diversified more slowly than angiosperms during the Neogene as a result of their higher extinction rate. Compared with angiosperms, living gymnosperm groups are not ancient. The fossil record also indicates that gymnosperms suffered major extinctions when climate changed in the Oligocene and Miocene. Extant gymnosperm groups occupy diverse habitats and some probably survived after making adaptive shifts.     

Synonymous and Nonsynonymous Substitution Rates in Diatoms: A Comparison Between Chloroplast and Nuclear Genes

Rates of synonymous and nonsynonymous nucleotide substitutions and codon usage bias (ENC) were estimated for a number of nuclear and chloroplast genes in a sample of centric and pennate diatoms. The results suggest that DNA evolution has taken place, on an average, at a slower rate in the chloroplast genes than in the nuclear genes: a rate variation pattern similar to that observed in land plants. Synonymous substitution rates in the chloroplast genes show a negative association with the degree of codon usage bias, suggesting that genes with a higher degree of codon usage bias have evolved at a slower rate. While this relationship has been shown in both prokaryotes and multicellular eukaryotes, it has not been demonstrated before in diatoms. Received: 3 June 1998 / Accepted: 11 August 1998     

Molecular evolution of chloroplast DNA sequences

Comparative data on the evolution of chloroplast genes are reviewed. The chloroplast genome has maintained a similar structural organization over most plant taxa so far examined. Comparisons of nucleotide sequence divergence among chloroplast genes reveals marked similarity across the plant kingdom and beyond to the cyanobacteria (blue-green algae). Estimates of rates of nucleotide substitution indicate a synonymous rate of 1.1 x 10(-9) substitutions per site per year. Noncoding regions also appear to be constrained in their evolution, although addition/deletion events are common. There have also been evolutionary changes in the distribution of introns in chloroplast encoded genes. Relative to mammalian mitochondrial DNA, the chloroplast genome evolves at a conservative rate.      

Abundance of plastid DNA insertions in nuclear genomes of rice and Arabidopsis

Pairwise comparison of whole plastid and draft nuclear genomic sequences of Arabidopsis thaliana and Oryza sativa L. ssp. indica shows that rice nuclear genomic sequences contain homologs of plastid DNA covering about 94 kb (83%) of plastid genome and including one or more full-length intact (without mutations resulting in premature stop codons) homologues of 26 known protein-coding (KPC) plastid genes. By contrast, only about 20 kb (16%) of chloroplast DNA, including a single intact plastid-derived KPC gene, is presented in the nucleus of A. thaliana. Sixteen rice plastid genes have at least one nuclear copy without any mutation or with only synonymous substitutions. Nuclear copies for other ten plastid genes contain both synonymous and non-synonymous substitutions. Multiple ESTs for 25 out of 26 KPC genes were also found, as well as putative promoters for some of them. The study of substitutions pattern shows that some of nuclear homologues of plastid genes may be functional and/or are under the pressure of the positive natural selection. The similar comparative analysis performed on rice chromosome 1 revealed 27 contigs containing plastid-derived sequences, totalling about 84 kb and covering two thirds of chloroplast DNA, with the intact nuclear copies of 26 different KPC genes. One of these contigs, AP003280, includes almost 57 kb (45%) of chloroplast genome with the intact copies of 22 KPC genes. At the same time, we observed that relative locations of homologues in plastid DNA and the nuclear genome are significantly different.     

Molecular Analysis of the Hot Spot Region Related to Length Mutations in Wheat Chloroplast Dnas. I. Nucleotide Divergence of Genes and Intergenic Spacer Regions Located in the Hot Spot Region

The nucleotide divergence of chloroplast DNAs around the hot spot region related to length mutation in Triticum (wheat) and Aegilops was analyzed. DNA sequences (ca. 4.5 kbp) of three chloroplast genome types of wheat complex were compared with one another and with the corresponding region of other grasses. The sequences region contained rbcL and psaI, two open reading frames, and a pseudogene, rpl23' (pseudogene for ribosomal protein L23) disrupted by AT-rich intergic spacer regions. The evolution of these genes in the closely related wheat complex is characterized by nonbiased nucleotide substitutions in terms of being synonymous/nonsynonymous, having A-T pressure transitions over transversions, and frequent changes at the third codon position, in contrast with the gene evolution among more distant plant groups where biased nucleotide substitutions have frequently occurred. The sequences of these genes had diverged almost in proportion to taxonomic distance. The sequence of the pseudogene rpl23' changed approximately two times faster than that of the coding region. Sequence comparison between the pseudogene and its protein-coding counterpart revealed different degrees of nucleotide homology in wheat, rice and maize, suggesting that the transposition timing of the pseudogene differed and/or that different rates of gene conversion operated on the pseudogene in the cpDNA of the three plant groups in Gramineae. The intergenic spacer regions diverged approximately ten times faster than the genes. The divergence of wheat from barley, and that from rice are estimated based on the nucleotide similarity to be 1.5, 10 and 40 million years, respectively.