Relative rates of synonymous substitutions in the mitochondrial, chloroplast and nuclear genomes of seed plants

Previous studies have estimated that, in angiosperms, the synonymous substitution rate of chloroplast genes is three times higher than that of mitochondrial genes and that of nuclear genes is twelve times higher than that of mitochondrial genes. Here we used 12 genes in 27 seed plant species to investigate whether these relative rates of substitutions are common to diverse seed plant groups. We find that the overall relative rate of synonymous substitutions of mitochondrial, chloroplast and nuclear genes of all seed plants is 1:3:10, that these ratios are 1:2:4 in gymnosperms but 1:3:16 in angiosperms and that they go up to 1:3:20 in basal angiosperms. Our results show that the mitochondrial, chloroplast and nuclear genomes of seed plant groups have different synonymous substitutions rates, that these rates are different in different seed plant groups and that gymnosperms have smaller ratios than angiosperms.     

蕨类植物叶绿体rps12基因的分子进化研究

以68种蕨类植物和2种石松类植物的rps12基因为对象,在系统发育背景下,结合最大似然法,使用HyPhy和PAML软件对该基因进行进化速率和适应性进化研究。结果显示:位于IR区的外显子2~3,其替换率明显降低,rps12基因编码序列的替换率也随之降低,且rps12基因密码子第3位的GC含量明显升高;在蕨类植物的进化过程中,3′-rps12更倾向定位于IR区,以保持较低的替换率;rps12基因编码的123个氨基酸位点中,共检测到4个正选择位点和116个负选择位点。研究结果表明基因序列进入到IR区后,显示出降低的替换率;强烈的负选择压力表明RPS12蛋白的高度保守性以及rps12基因的功能和结构已经趋于稳定。     

Complete nucleotide sequence of the chloroplast genome from a leptosporangiate fern, Adiantum capillus-veneris L.

We determined the complete nucleotide sequence of the chloroplast genome of the leptosporangiate fern, Adiantum capillus-veneris L. (Pteridaceae). The circular genome is 150,568 bp, with a large single-copy region (LSC) of 82,282 bp, a small-single copy region (SSC) of 21,392 bp and inverted repeats (IR) of 23,447 bp each. We compared the sequence to other published chloroplast genomes to infer the location of putative genes. When the IR is considered only once, we assigned 118 genes, of which 85 encode proteins, 29 encode tRNAs and 4 encode rRNAs. Four protein-coding genes, all four rRNA genes and six tRNA genes occur in the IR. Most (57) putative protein-coding genes appear to start with an ATG codon, but we also detected five other possible start codons, some of which suggest tRNA editing. We also found 26 apparent stop codons in 18 putative genes, also suggestive of RNA editing. We found all but one of the tRNA genes necessary to encode the complete repertoire required for translation. The missing trnK gene appears to have been disrupted by a large inversion, relative to other published chloroplast genomes. We detected several structural rearrangements that may provide useful information for phylogenetic studies.     

Complete Sequencing of Five Araliaceae Chloroplast Genomes and the Phylogenetic Implications

Background

The ginseng family (Araliaceae) includes a number of economically important plant species. Previously phylogenetic studies circumscribed three major clades within the core ginseng plant family, yet the internal relationships of each major group have been poorly resolved perhaps due to rapid radiation of these lineages. Recent studies have shown that phyogenomics based on chloroplast genomes provides a viable way to resolve complex relationships.

Methodology/Principal Findings

We report the complete nucleotide sequences of five Araliaceae chloroplast genomes using next-generation sequencing technology. The five chloroplast genomes are 156,333–156,459 bp in length including a pair of inverted repeats (25,551–26,108 bp) separated by the large single-copy (86,028–86,566 bp) and small single-copy (18,021–19,117 bp) regions. Each chloroplast genome contains the same 114 unique genes consisting of 30 transfer RNA genes, four ribosomal RNA genes, and 80 protein coding genes. Gene size, content, and order, AT content, and IR/SC boundary structure are similar among all Araliaceae chloroplast genomes. A total of 140 repeats were identified in the five chloroplast genomes with palindromic repeat as the most common type. Phylogenomic analyses using parsimony, likelihood, and Bayesian inference based on the complete chloroplast genomes strongly supported the monophyly of the Asian Palmate group and the Aralia-Panax group. Furthermore, the relationships among the sampled taxa within the Asian Palmate group were well resolved. Twenty-six DNA markers with the percentage of variable sites higher than 5% were identified, which may be useful for phylogenetic studies of Araliaceae.

Conclusion

The chloroplast genomes of Araliaceae are highly conserved in all aspects of genome features. The large-scale phylogenomic data based on the complete chloroplast DNA sequences is shown to be effective for the phylogenetic reconstruction of Araliaceae.     

Complete Chloroplast Genome of Sedum sarmentosum and Chloroplast Genome Evolution in Saxifragales

Comparative chloroplast genome analyses are mostly carried out at lower taxonomic levels, such as the family and genus levels. At higher taxonomic levels, chloroplast genomes are generally used to reconstruct phylogenies. However, little attention has been paid to chloroplast genome evolution within orders. Here, we present the chloroplast genome of Sedum sarmentosum and take advantage of several available (or elucidated) chloroplast genomes to examine the evolution of chloroplast genomes in Saxifragales. The chloroplast genome of S. sarmentosum is 150,448 bp long and includes 82,212 bp of a large single-copy (LSC) region, 16.670 bp of a small single-copy (SSC) region, and a pair of 25,783 bp sequences of inverted repeats (IRs).The genome contains 131 unique genes, 18 of which are duplicated within the IRs. Based on a comparative analysis of chloroplast genomes from four representative Saxifragales families, we observed two gene losses and two pseudogenes in Paeonia obovata, and the loss of an intron was detected in the rps16 gene of Penthorum chinense. Comparisons among the 72 common protein-coding genes confirmed that the chloroplast genomes of S. sarmentosum and Paeonia obovata exhibit accelerated sequence evolution. Furthermore, a strong correlation was observed between the rates of genome evolution and genome size. The detected genome size variations are predominantly caused by the length of intergenic spacers, rather than losses of genes and introns, gene pseudogenization or IR expansion or contraction. The genome sizes of these species are negatively correlated with nucleotide substitution rates. Species with shorter duration of the life cycle tend to exhibit shorter chloroplast genomes than those with longer life cycles.     

Distinctive Architecture of the Chloroplast Genome in the Chlorodendrophycean Green Algae Scherffelia dubia and Tetraselmis sp. CCMP 881

The Chlorodendrophyceae is a small class of green algae belonging to the core Chlorophyta, an assemblage that also comprises the Pedinophyceae, Trebouxiophyceae, Ulvophyceae and Chlorophyceae. Here we describe for the first time the chloroplast genomes of chlorodendrophycean algae (Scherffelia dubia, 137,161 bp; Tetraselmis sp. CCMP 881, 100,264 bp). Characterized by a very small single-copy (SSC) region devoid of any gene and an unusually large inverted repeat (IR), the quadripartite structures of the Scherffelia and Tetraselmis genomes are unique among all core chlorophytes examined thus far. The lack of genes in the SSC region is offset by the rich and atypical gene complement of the IR, which includes genes from the SSC and large single-copy regions of prasinophyte and streptophyte chloroplast genomes having retained an ancestral quadripartite structure. Remarkably, seven of the atypical IR-encoded genes have also been observed in the IRs of pedinophycean and trebouxiophycean chloroplast genomes, suggesting that they were already present in the IR of the common ancestor of all core chlorophytes. Considering that the relationships among the main lineages of the core Chlorophyta are still unresolved, we evaluated the impact of including the Chlorodendrophyceae in chloroplast phylogenomic analyses. The trees we inferred using data sets of 79 and 108 genes from 71 chlorophytes indicate that the Chlorodendrophyceae is a deep-diverging lineage of the core Chlorophyta, although the placement of this class relative to the Pedinophyceae remains ambiguous. Interestingly, some of our phylogenomic trees together with our comparative analysis of gene order data support the monophyly of the Trebouxiophyceae, thus offering further evidence that the previously observed affiliation between the Chlorellales and Pedinophyceae is the result of systematic errors in phylogenetic reconstruction.     

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     

Migration of endpoints of two genes relative to boundaries between regions of the plastid genome in the grass family (Poaceae)

Overlapping genes occur widely in microorganisms and in some plastid genomes, but unique properties are observed when such genes span the boundaries between single-copy and repeat regions. The termini of ndhH and ndhF, situated near opposite ends of the small single-copy region (SSC) in the plastid genomes of grasses (Poaceae), have migrated repeatedly into and out of the adjacent inverted-repeat regions (IR). The two genes are transcribed in the same direction, and the 5' terminus of ndhH extends into the IR in some species, while the 3' terminus of ndhF extends into the IR in others. When both genes extend into the IR, portions of the genes overlap and are encoded by the same nucleotide positions. Fine-scale mapping of the SSC-IR junctions across a sample of 92 grasses and outgroups, integrated into a phylogenetic analysis, indicates that the earliest grasses resembled the related taxa Joinvillea (Joinvilleaceae) and Ecdeiocolea (Ecdeiocoleaceae), with ca. 180 nucleotides of ndhH extending into the IR, and with ndhF confined to the SSC. This structure is maintained in early-diverging grass lineages and in most species of the BEP clade. In the PACMAD clade, ndhH lies completely or nearly completely within the SSC, and ca. 20 nucleotides of ndhF extend into the IR. The nucleotide substitution rate has increased in the PACMAD clade in the portion of ndhH that has migrated into the SSC.     

Comparative genomics of four Liliales families inferred from the complete chloroplast genome sequence of Veratrum patulum O. Loes. (Melanthiaceae)

The sequence of the chloroplast genome, which is inherited maternally, contains useful information for many scientific fields such as plant systematics, biogeography and biotechnology because its characteristics are highly conserved among species. There is an increase in chloroplast genomes of angiosperms that have been sequenced in recent years. In this study, the nucleotide sequence of the chloroplast genome (cpDNA) of Veratrum patulum Loes. (Melanthiaceae, Liliales) was analyzed completely. The circular double-stranded DNA of 153,699 bp consists of two inverted repeat (IR) regions of 26,360 bp each, a large single copy of 83,372 bp, and a small single copy of 17,607 bp. This plastome contains 81 protein-coding genes, 30 distinct tRNA and four genes of rRNA. In addition, there are six hypothetical coding regions (ycf1, ycf2, ycf3, ycf4, ycf15 and ycf68) and two open reading frames (ORF42 and ORF56), which are also found in the chloroplast genomes of the other species. The gene orders and gene contents of the V. patulum plastid genome are similar to that of Smilax china, Lilium longiflorum and Alstroemeria aurea, members of the Smilacaceae, Liliaceae and Alstroemeriaceae (Liliales), respectively. However, the loss rps16 exon 2 in V. patulum results in the difference in the large single copy regions in comparison with other species. The base substitution rate is quite similar among genes of these species. Additionally, the base substitution rate of inverted repeat region was smaller than that of single copy regions in all observed species of Liliales. The IR regions were expanded to trnH_GUG in V. patulum, a part of rps19 in L. longiflorum and A. aurea, and whole sequence of rps19 in S. china. Furthermore, the IGS lengths of rbcL-accD-psaI region were variable among Liliales species, suggesting that this region might be a hotspot of indel events and the informative site for phylogenetic studies in Liliales. In general, the whole chloroplast genome of V. patulum, a potential medicinal plant, will contribute to research on the genetic applications of this genus.     

桃儿七叶绿体比较基因组学分析

为探究桃儿七(Sinopodophyllum hexandrum)不同叶绿体基因组特征,本研究以桃儿七5个叶绿体基因组为研究对象,借助生物信息学工具进行基因组图谱构建、重复序列分析、密码子偏好分析、反向重复序列区(inverted repeat, IR)/单拷贝区(single-copy, SC)边界分析、基因组序列比较分析及系统发育分析。结果表明:桃儿七5个叶绿体基因组全长为157 203–157 940 bp,为典型的叶绿体四分体结构,共注释出133–137个基因,说明桃儿七叶绿体基因组具有多样性。桃儿七不同叶绿体基因组简单重复序列(simple sequence repeat, SSR)位点不同,单核苷酸A/T占主要优势,散在重复序列包括正向重复、回文重复和反向重复3类。密码子偏好分析显示有效密码子(effective number of codon,ENc)值为51.14–51.17,密码子偏好性弱,GC与GC3s所占比例小于50%,密码子偏向使用A和U碱基并且以A和U碱基结尾。桃儿七5个叶绿体基因组IR/SC边界和基因组序列均比较保守。系统发育分析结果表明桃儿七和北美桃儿七亲...     

Rates of synonymous substitution in plant nuclear genes

Summary The rate of synonymous nucleotide substitution in nuclear genes of higher plants has been estimated. The rate varies among genes by a factor of up to two, in a manner that is not immediately explicable in terms of base composition or codon usage bias. The average rate, in both monocots and dicots, is about four times higher than that in chloroplast genes. This leads to an estimated absolute silent substitution rate of 6 × 10^–9 substitutions per site per year that falls within the range of average rates (2–8 × 10^–9) seen in different mammalian nuclear genomes.     

多花海棠叶绿体基因组分析

多花海棠（Malus floribunda Siebold.）是世界范围内广泛栽培的苹果属物种,具有较高的观赏价值和育种意义。对其进行叶绿体基因组比较分析,有利于完善苹果属系统进化以及种质利用的研究内容。基于全基因组测序数据,组装获得一个完整的具有四分体结构的多花海棠叶绿体基因组。该基因组包括大单拷贝区（88 142 bp）、反向重复区B （26 353 bp）、小单拷贝区（19 189 bp）与反向重复区A （26 353 bp）,共计160 037 bp。多花海棠叶绿体全基因组共注释到111个基因,包括78个蛋白编码基因、29个tRNA基因和4个rRNA基因。此外,在其基因组中识别到大量的重复序列,与三叶海棠和变叶海棠略有差异。通过计算相对同义密码子使用度,发现其高频密码子共30种,并且密码子具有偏向A/T结尾的使用模式。种间序列比对、边界分析的结果表明,大单拷贝区序列变异较大,8种苹果属植物SC区与IR区扩张收缩情况整体上较为相似。基于叶绿体基因组序列的系统进化分析,将多花海棠、湖北海棠和变叶海棠聚为一类。多花海棠叶绿体基因组的研究可为今后遗传标记开发与种质资源利用等提供数据支持。     

木犀科植物叶绿体基因组结构特征和系统发育关系

木犀科11属19个种叶绿体基因组的一般特征和变异特征的比较分析显示, 结果表明, 该科叶绿体基因组大小为154-165 kb, 其差异主要是大单拷贝(LSC)长度的差异所致。Jasminum属3个物种的叶绿体基因组长度与其余物种有较大差异, 该属clpP基因内含子和accD基因丢失。共线性分析表明, Jasminum属3个物种多个基因出现基因重排现象, 倒位可能是重排的主要原因。Jasminum属在IRb/SSC和SSC/IRa边界的基因均与其它物种不同; 重复序列与SSR数量检测结果表明, Jasminum属与其余物种在数量及重复长度上差异较大。基于CDS数据构建的系统发育树表明, Abeliophyllum distichum和Forsythia suspensa为木犀科中较早分化的类群。     

Whole chloroplast genome comparison of rice,maize, and wheat: implications for chloroplast gene diversification and phylogeny of cereals

The fully sequenced chloroplast genomes of maize (subfamily Panicoideae), rice (subfamily Bambusoideae), and wheat (subfamily Pooideae) provide the unique opportunity to investigate the evolution of chloroplast genes and genomes in the grass family (Poaceae) by whole-genome comparison. Analyses of nucleotide sequence variations in 106 cereal chloroplast genes with tobacco sequences as the outgroup suggested that (1) most of the genic regions of the chloroplast genomes of maize, rice, and wheat have evolved at similar rates; (2) RNA genes have highly conservative evolutionary rates relative to the other genes; (3) photosynthetic genes have been under strong purifying selection; (4) between the three cereals, 14 genes which account for about 28% of the genic region have evolved with heterogeneous nucleotide substitution rates; and (5) rice genes tend to have evolved more slowly than the others at loci where rate heterogeneity exists. Although the mechanism that underlies chloroplast gene diversification is complex, our analyses identified variation in nonsynonymous substitution rates as a genetic force that generates heterogeneity, which is evidence of selection in chloroplast gene diversification at the intrafamilial level. Phylogenetic trees constructed with the variable nucleotide sites of the chloroplast genes place maize basal to the rice-wheat clade, revealing a close relationship between the Bambusoideae and Pooideae.     

梅花草属叶绿体基因组进化分析

叶绿体基因组序列变异和基因组成等特征可有效反映植物类群间的系统发育和进化关系。本研究利用Illumina高通量测序平台对梅花草属（Parnassia）及其近缘属5种植物的叶绿体基因组进行测序和组装,同时基于已发表的近缘种叶绿体基因组信息,对梅花草属叶绿体基因组结构特征、序列遗传变异和蛋白编码基因密码子偏好性比对分析。结果显示：梅花草属叶绿体基因组整体结构较为保守,均为四分体结构;梅花草多个基因出现假基因化,而本属其他物种叶绿体基因组成一致,均编码115个基因;与近缘属物种相比,本属所有物种均丢失rpl16基因的内含子;蛋白质编码基因的非同义/同义替代率比值较低,叶绿体基因可能经历纯化选择作用;密码子偏好性聚类结果与蛋白编码序列重建的系统发育关系结果一致。本研究表明选择压力可能在梅花草属叶绿体基因组蛋白编码基因进化过程中发挥作用,有助于进一步理解梅花草属植物的进化和适应机制。     

Chloroplast genome characteristics and phylogenetic analysis of the medicinal plant Blumea balsamifera (L.) DC

Blumea balsamifera (L.) DC., a medicinal plant with high economic value in the Asteraceae family, is widely distributed in China and Southeast Asia. However, studies on the population structure or phylogenetic relationships with other related species are rare owing to the lack of genome information. In this study, through high-throughput sequencing, we found that the chloroplast genome of B. balsamifera was 151,170 bp in length, with a pair of inverted repeat regions (IRa and IRb) comprising 24,982 bp, a large single-copy (LSC) region comprising 82,740 bp, and a small single-copy (SSC) region comprising 18,466 bp. A total of 130 genes were identified in the chloroplast genome of B. balsamifera, including 85 protein-coding, 37 transfer RNA, and 8 ribosomal RNA genes; furthermore, sequence analysis identified 53 simple sequence repeats. Whole chloroplast genome comparison indicated that the inverted regions (IR) were more conserved than large single-copy and SSC regions. Phylogenetic analysis showed that B. balsamifera is closely related to Pluchea indica. Conclusively, the chloroplast genome of B. balsamifera was helpful for species identification and analysis of the genetic diversity and evolution in the genus Blumea and family Asteraceae.     

木犀科植物叶绿体基因组结构特征和系统发育关系

木犀科11属19个种叶绿体基因组的一般特征和变异特征的比较分析显示, 结果表明, 该科叶绿体基因组大小为154-165 kb, 其差异主要是大单拷贝(LSC)长度的差异所致。Jasminum属3个物种的叶绿体基因组长度与其余物种有较大差异, 该属clpP基因内含子和accD基因丢失。共线性分析表明, Jasminum属3个物种多个基因出现基因重排现象, 倒位可能是重排的主要原因。Jasminum属在IRb/SSC和SSC/IRa边界的基因均与其它物种不同; 重复序列与SSR数量检测结果表明, Jasminum属与其余物种在数量及重复长度上差异较大。基于CDS数据构建的系统发育树表明, Abeliophyllum distichum和Forsythia suspensa为木犀科中较早分化的类群。     

基于全基因组从头测序技术盐肤木叶绿体基因组的测序分析

盐肤木是一种重要的经济树种,可为医药和工业染料提供原料。盐肤木具有较强的抗旱、耐寒、耐盐,可在温带、暖温带和亚热带地区生长。本研究首次对盐肤木叶绿体基因组进行从头测序(de novo sequencing)组装研究。结果表明,盐肤木叶绿体基因组长度为159082 bp,具有典型的四部分结构,两个单拷贝区被一对反向重复区分隔。LSC和SSC的长度分别为85394 bp和18663 bp。叶绿体基因组总共编码126个基因,其中包括88个蛋白编码基因,8个rRNA基因,30个tRNA基因。在叶绿体基因组中,61.97%的序列为基因编码区。在盐肤木叶绿体基因组中,只有8个基因含有内含子,除ycf3基因(2个内含子)外,其余均含有1个内含子。盐肤木叶绿体基因组总共存在755个SSR位点。SSR主要由二核苷酸和单核苷酸组成,分别占60%(453)和28.74%(217)。聚类分析结果表明,漆树科与盐肤木最为接近,其次为槭树科和无患子科。本研究为盐肤木的分类提供了分子基础。本研究是关于盐肤木叶绿体基因组的首次报道,对了解其光合作用、进化和叶绿体转基因工程具有重要意义。