Characterization of chloroplast DNA microsatellites from Saccharum spp and related species

Microsatellites, or simple sequence repeats (SSRs), and their flanking regions in chloroplast genomes (plastomes) of some species of the family Poaceae were analyzed in silico to look for DNA sequence variations. Comparison of the complete chloroplast DNA sequences (cpDNAs) of sugarcane (Saccharum hybrid cv. SP-80-3280 and S. officinarum cv. NCo310) and related species, Agrostis stolonifera, Brachypodium distachyon, Hordeum vulgare subsp vulgare, Lolium perenne, Oryza nivara, O. sativa subsp indica, O. sativa subsp japonica, Sorghum bicolor, Triticum aestivum, Zea mays, and Z. mays cv. B73, allowed us to examine the organization of chloroplast SSRs (cpSSRs) in genic and intergenic regions. We identified 204 cpSSRs in the sugarcane cpDNA; 22.5% were in genic regions. The ndh, rps, trn, and rpl gene clusters of the chloroplasts had the most repeats. Mononucleotide repeats were the most abundant cpSSRs in these species; however, di-, tri-, tetra-, penta-, and hexanucleotide repeats were also identified. Many base substitutions and deletions/insertions were identified in the cpSSR loci and their flanking regions. Multiple alignments of all cpSSR sequences of Poaceae species made identification of nucleotide variability possible; repeat motifs are not uniformly distributed across the Poaceae plastomes, but are mostly confined to intergenic regions. Phylogeny was determined by maximum parsimony and neighbor-joining inference methods. The cpSSRs of these species were found to be polymorphic. It appears that individual cpSSRs in the Poaceae are stable, at least over short periods of evolutionary time. We conclude that the plastome database can be exploited for phylogenetic analysis and biotechnological development.     

Complete nucleotide sequence of the sugarcane (Saccharum officinarum) chloroplast genome: a comparative analysis of four monocot chloroplast genomes.

The complete nucleotide sequence of the chloroplast genome of sugarcane (Saccharum officinarum) has been determined. It is a circular double-stranded DNA molecule, 141,182 bp in size, and is composed of a large single copy of 83,048 bp, a small single copy of 12,544 bp, and a pair of inverted repeat regions of 22,795 bp each. A comparative analysis among monocots showed that the sugarcane chloroplast genome was very similar to maize but not to rice or wheat. Between sugarcane and maize at the rps16-trnQ (UUG) region, however, a length polymorphism was identified. With regard to insertions/deletions equal to or longer than 5 bp, a total of 53 insertion and 31 deletion events were identified in the sugarcane chloroplast genome. Of the 84 loci identified, a pair of direct repeat sequences was located side by side in a tandem fashion in 47 loci (56.0%). A recombination event during plant evolution is discussed at two sites between the sugarcane and tobacco chloroplast genomes.     

Complete chloroplast genomes of Saccharum spontaneum,Saccharum officinarum and Miscanthus floridulus (Panicoideae: Andropogoneae) reveal the plastid view on sugarcane origins

Sugarcane (Saccharum hybrid cultivar) ranks among the world's top 10 food crops and annually provides 60–70% of the sugar produced worldwide. Despite its economic importance there has been no large-scale systematics study of genus Saccharum and the existing model of sugarcane origins has remained largely unchallenged for almost 50 years. For the first time, we have assembled the complete plastid genomes of Miscanthus floridulus (first report for this genus), Saccharum spontaneum and Saccharum officinarum allowing us to elucidate the phylogenetic origins of Saccharum s.s. species. We demonstrate that Saccharum s.s. is divided into four species, with S. spontaneum diverging from the remainder of the genus about 1.5 million years ago and S. robustum diverging 750,000 years ago. Two separate lineages, one leading to S. officinarum and the other leading to modern hybrid cultivars diverged from S. robustum 640,000 years ago. These findings overturn all previous hypotheses on sugarcane origins, demonstrating that sugarcane's antecedents could not have arisen by human action. All modern cultivars share a common Polynesian origin, whereas Old World canes, S. barberi and S. sinense, cluster as a distinct S. officinarum lineage. This makes modern cultivars a distinct species of genus Saccharum, and we formally propose the name Saccharum cultum for the ancestor of all lineages currently classified as Saccharum hybrid cultivars.     

三种厚朴叶绿体基因组的比较研究

为了深入发掘日本厚朴、厚朴、凹叶厚朴叶绿体基因组差异,筛选厚朴优良性状候选基因,开展三种厚朴的分子遗传研究,该文利用Illumina HiSeq高通量测序平台首次对日本厚朴叶绿体进行测序、组装,并与已有的厚朴、凹叶厚朴叶绿体基因组共同注释,获得三个物种叶绿体基因图谱,筛选出三个基因组中的差异基因,又与同科中11个亲缘物种进行叶绿体基因组比对,构建NJ遗传树。结果表明：(1)日本厚朴叶绿体基因组的Clean Reads为19 791 019,Q30为91.33%,组装后基因组全长160 051 bp, GC含量为39.2%,含tRNA 37个,rRNA 8个。(2)比对分析发现三种厚朴具有相似的IR、LSC和SSC结构,以及GC含量和tRNA数量,但编码基因种类和数量、内含子和外显子的数量和结构等存在差异。(3)日本厚朴的功能基因数目较厚朴、凹叶厚朴分别多6个和4个,主要分布于LSC区和IR区,涉及核糖体大亚基、核糖体小亚基和未知功能基因类群。(4)系统发育分析结果进一步显示日本厚朴与凹叶厚朴亲缘关系较近,其次是厚朴。该研究表明日本厚朴具有更丰富的叶绿体基因组结构、组成和变异特征,是其适...     

Genetic analyses of Oryza species by molecular markers for chloroplast genomes

Summary Relationships in a wide range of Oryza species (13 species) were analyzed using the large subunits (LS) of Fraction I protein (Rubisco) and the Bam HI restriction patterns of chloroplast DNA (ctDNA) as molecular markers. Four types of LS were detected by isoelectrofocusing with and without S-carboxymethylation. The close relation between AA and CCDD genome species was suggested by analyses of LS and ctDNA. Intraspecific variation in O. latifolia was detected at the levels of both LS and ctDNA. The LS of the BB, BBCC, and CC genomes and FF (O. brachyantha) were not distinguishable, although the native Rubisco of the latter was slightly different from those of the first three. It was also shown that O. australiensis, the only EE genome species, might have evolved differently than the other Oryza species.     

Comparative genome sequence and phylogenetic analysis of chloroplast for evolutionary relationship among Pinus species

Genus Pinus is a widely dispersed genus of conifer plants in the Northern Hemisphere. However, the inadequate accessibility of genomic knowledge limits our understanding of molecular phylogeny and evolution of Pinus species. In this study, the evolutionary features of complete plastid genome and the phylogeny of the Pinus genus were studied. A total of thirteen divergent hotspot regions (trnk-UUU, matK, trnQ-UUG, atpF, atpH, rpoC1, rpoC2, rpoB, ycf2, ycf1, trnD-GUC, trnY-GUA, and trnH-GUG) were identified that would be utilized as possible genetic markers for determination of phylogeny and population genetics analysis of Pinus species. Furthermore, seven genes (petD, psaI, psaM, matK, rps18, ycf1, and ycf2) with positive selection site in Pinus species were identified. Based on the whole genome this phylogenetic study showed that twenty-four Pinus species form a significant genealogical clade. Divergence time showed that the Pinus species originated about 100 million years ago (MYA) (95% HPD, 101.76.35–109.79 MYA), in lateral stages of Cretaceous. Moreover, two of the subgenera are consequently originated in 85.05 MYA (95% HPD, 81.04–88.02 MYA). This study provides a phylogenetic relationship and a chronological framework for the future study of the molecular evolution of the Pinus species.     

Microsatellites and microsynteny in the chloroplast genomes of Oryza and eight other Gramineae species

Primer pairs flanking ten chloroplast microsatellite loci, originally identified in Oryza sativa cv Nipponbare, were evaluated for amplification and allelic diversity using a panel of 13 diverse cultivars of rice (O. sativa), 19 accessions of wild rice (three O. officinalis, five O. latifolia, five O. minuta, four O. australiensis, one O. brachyantha and one O. ridleyi) and eight other Gramineae species (maize, teosinte, wheat, oat, barley, pearl millet, sorghum and sugarcane). Amplified products were obtained for all samples at nine out of ten loci. Among the rice cultivars, the number of alleles per locus ranged from one to four, with monomorphic patterns observed at five loci. The average polymorphism information content (PIC) value at the other five (polymorphic) loci was 0.54 among the 13 cultivars. When wild rice and the other Gramineae species were compared based on the proportion of shared alleles, their phylogenetic relationships were in agreement with previous studies using different types of markers; however, the magnitude of the differences based on chloroplast microsatellites underestimated the genetic distance separating these divergent species and genera. A sequence-based comparison of homologous regions of the rice and maize chloroplast genomes revealed that, while a high level of microsynteny is evident, the occurrence of actively evolving microsatellite motifs in specific regions of the rice chloroplast genome appears to be mainly a species or genome-specific phenomenon. Thus the chloroplast primer pairs used in this study bracketed mutationally active microsatellite motifs in rice but degenerate, interrupted motifs or highly conserved, mutationally inert motifs in distantly related genera. Received: 17 March 1999 / Accepted: 11 November 1999     

The complete chloroplast genomes of three Hamamelidaceae species: Comparative and phylogenetic analyses

Hamamelidaceae is an important group that represents the origin and early evolution of angiosperms. Its plants have many uses, such as timber, medical, spice, and ornamental uses. In this study, the complete chloroplast genomes of Loropetalum chinense (R. Br.) Oliver, Corylopsis glandulifera Hemsl., and Corylopsis velutina Hand.‐Mazz. were sequenced using the Illumina NovaSeq 6000 platform. The sizes of the three chloroplast genomes were 159,402 bp (C. glandulifera), 159,414 bp (C. velutina), and 159,444 bp (L. chinense), respectively. These chloroplast genomes contained typical quadripartite structures with a pair of inverted repeat (IR) regions (26,283, 26,283, and 26,257 bp), a large single‐copy (LSC) region (88,134, 88,146, and 88,160 bp), and a small single‐copy (SSC) region (18,702, 18,702, and 18,770 bp). The chloroplast genomes encoded 132–133 genes, including 85–87 protein‐coding genes, 37–38 tRNA genes, and 8 rRNA genes. The coding regions were composed of 26,797, 26,574, and 26,415 codons, respectively, most of which ended in A/U. A total of 37–43 long repeats and 175–178 simple sequence repeats (SSRs) were identified, and the SSRs contained a higher number of A + T than G + C bases. The genome comparison showed that the IR regions were more conserved than the LSC or SSC regions, while the noncoding regions contained higher variability than the gene coding regions. Phylogenetic analyses revealed that species in the same genus tended to cluster together. Chunia Hung T. Chang, Mytilaria Lecomte, and Disanthus Maxim. may have diverged early and Corylopsis Siebold & Zucc. was closely related to Loropetalum R. Br. This study provides valuable information for further species identification, evolution, and phylogenetic studies of Hamamelidaceae plants.     

A comparison of rice chloroplast genomes

Using high quality sequence reads extracted from our whole genome shotgun repository, we assembled two chloroplast genome sequences from two rice (Oryza sativa) varieties, one from 93-11 (a typical indica variety) and the other from PA64S (an indica-like variety with maternal origin of japonica), which are both parental varieties of the super-hybrid rice, LYP9. Based on the patterns of high sequence coverage, we partitioned chloroplast sequence variations into two classes, intravarietal and intersubspecific polymorphisms. Intravarietal polymorphisms refer to variations within 93-11 or PA64S. Intersubspecific polymorphisms were identified by comparing the major genotypes of the two subspecies represented by 93-11 and PA64S, respectively. Some of the minor genotypes occurring as intravarietal polymorphisms in one variety existed as major genotypes in the other subspecific variety, thus giving rise to intersubspecific polymorphisms. In our study, we found that the intersubspecific variations of 93-11 (indica) and PA64S (japonica) chloroplast genomes consisted of 72 single nucleotide polymorphisms and 27 insertions or deletions. The intersubspecific polymorphism rates between 93-11 and PA64S were 0.05% for single nucleotide polymorphisms and 0.02% for insertions or deletions, nearly 8 and 10 times lower than their respective nuclear genomes. Based on the total number of nucleotide substitutions between the two chloroplast genomes, we dated the divergence of indica and japonica chloroplast genomes as occurring approximately 86,000 to 200,000 years ago.     

Variation within and among the chloroplast genomes ofMelaleuca alternifolia andM. linariifolia (Myrtaceae)

Melaleuca alternifolia andM. linariifolia are commercially important Australian species harvested for their essential oils. Both species have relatively narrow and disjunct distributions on the central coast of eastern Australia. Variation in the chloroplast genome was assessed for eight individuals from each of twelve populations, representing the species' geographic range. Low nucleotide diversity withinM. alternifolia contrasted with high nucleotide diversity inM. linariifolia. CpDNA data are consistent with the southern population ofM. alternifolia being a hybrid population withM. linariifolia. The two species are sympatric in this region. Variation inM. linariifolia was geographically structured, with northern populations differing from southern populations by seven restriction site mutations, five length mutations and an inversion. There was no evidence of hybridisation of the cp genome of northernM. linariifolia with the partially sympatric speciesM. trichostachya. Intra- and interspecific variation in the chloroplast genomes ofM. alternifolia, M. linariifolia, andM. trichostachya indicate considerable potential for the use of intraspecific cpDNA studies in examining phylogenetic relationships in melaleucas.     

The chloroplast genomes of azuki bean and its close relatives: a deletion mutation found in weed azuki bean

A physical map of the azuki bean (Vigna angularis cv. Erimo-shozu) chloroplast DNA (cpDNA) was constructed by localising the cleavage sites of PstI, SalI, SmaI, SacI, KpnI, PvuII and XhoI. The resulting map is more similar to the cpDNA-maps of two Vigna species (mung bean, V. radiata, and V. nakashimae) than to common bean (Phaseolus vulgaris) cpDNA map. Azuki bean was originally classified in the genus Phaseolus, and the inclusion of this taxon in the genus Vigna is a recent taxonomic decision. Our result is thus in favour of the taxonomic placement of azuki bean in the same genus as V. nakashimae and mung bean. We also found that a weed-form accession of azuki bean has a 96-bp deletion relative to the cultivar Erimo-shozu. The 96-bp deletion is located between the trnS-UGA and psbC genes in the large single-copy region of the chloroplast genome. This deletion is flanked by imperfect 9-bp direct repeats, suggesting that the deletion was a result of intra-molecular recombination mediated by these direct repeats.     

Fine structural features of the chloroplast genome: comparison of the sequenced chloroplast genomes.

The entire nucleotide sequences of the rice, tobacco and liverwort chloroplast genomes have been determined. We compared all the chloroplast genes, open reading frames and spacer regions in the plastid genomes of these three species in order to elucidate general structural features of the chloroplast genome. Analyses of homology, GC content and codon usage of the genes enabled us to classify them into two groups: photosynthesis genes and genetic system genes. Based on comparisons of homology, GC content and codon usage, unidentified ORFs can also be assigned to each of these groups such that it is possible to speculate about the functions of products which may be produced by these ORFs. The spacer regions and intron sequences were compared and found to have no obvious homology between rice and liverwort or between tobacco and liverwort.     

木芙蓉三个品种及近缘种的叶绿体基因组比较分析

木芙蓉(Hibiscus mutabilis)栽培历史悠久,是原产中国的古老园林树种和药用植物。为了探讨木芙蓉品种及近缘种的进化特征,厘清木芙蓉品种间及其与近缘种间的亲缘关系,以及探究木芙蓉叶绿体基因组(chloroplast DNA, cpDNA)的遗传方式,该文选择了一个杂交组合中的3个木芙蓉栽培品种(‘单瓣白’‘金秋颂’‘牡丹粉’),用高通量测序平台Illumina NovaSeq对其cpDNA进行首次测序。经组装注释后得到3条完整的cpDNA序列,结合该团队已经完成的近缘种台湾芙蓉(H.taiwanensis)和来自基因库的木槿、朱槿的cpDNA,对木槿属4种及木芙蓉种下的3个品种进行了cpDNA组成和结构特征的比较分析,并完成了其系统发育树重建。结果表明：(1)‘单瓣白’‘金秋颂’‘牡丹粉’的cpDNA序列长度分别为160 880、160 879、160 920 bp,基因数目均为130个,其中蛋白编码基因85个、核糖体RNA 8个和转运RNA 37个。(2)比较分析结果显示,木芙蓉的种下3个品种及其近缘种台湾芙蓉在cpDNA上高度保守,反向重复区(IR)均为26 300 b...     

Six newly sequenced chloroplast genomes from prasinophyte green algae provide insights into the relationships among prasinophyte lineages and the diversity of streamlined genome architecture in picoplanktonic species

Background

Because they represent the earliest divergences of the Chlorophyta, the morphologically diverse unicellular green algae making up the prasinophytes hold the key to understanding the nature of the first viridiplants and the evolutionary patterns that accompanied the radiation of chlorophytes. Nuclear-encoded 18S rDNA phylogenies unveiled nine prasinophyte clades (clades I through IX) but their branching order is still uncertain. We present here the newly sequenced chloroplast genomes of Nephroselmis astigmatica (clade III) and of five picoplanktonic species from clade VI (Prasinococcus sp. CCMP 1194, Prasinophyceae sp. MBIC 106222 and Prasinoderma coloniale) and clade VII (Picocystis salinarum and Prasinophyceae sp. CCMP 1205). These chloroplast DNAs (cpDNAs) were compared with those of the six previously sampled prasinophytes (clades I, II, III and V) in order to gain information both on the relationships among prasinophyte lineages and on chloroplast genome evolution.

Results

Varying from 64.3 to 85.6 kb in size and encoding 100 to 115 conserved genes, the cpDNAs of the newly investigated picoplanktonic species are substantially smaller than those observed for larger-size prasinophytes, are economically packed and contain a reduced gene content. Although the Nephroselmis and Picocystis cpDNAs feature a large inverted repeat encoding the rRNA operon, gene partitioning among the single copy regions is remarkably different. Unexpectedly, we found that all three species from clade VI (Prasinococcales) harbor chloroplast genes not previously documented for chlorophytes (ndhJ, rbcR, rpl21, rps15, rps16 and ycf66) and that Picocystis contains a trans-spliced group II intron. The phylogenies inferred from cpDNA-encoded proteins are essentially congruent with 18S rDNA trees, resolving with robust support all six examined prasinophyte lineages, with the exception of the Pycnococcaceae.

Conclusions

Our results underscore the high variability in genome architecture among prasinophyte lineages, highlighting the strong pressure to maintain a small and compact chloroplast genome in picoplanktonic species. The unique set of six chloroplast genes found in the Prasinococcales supports the ancestral status of this lineage within the prasinophytes. The widely diverging traits uncovered for the clade-VII members (Picocystis and Prasinophyceae sp. CCMP 1205) are consistent with their resolution as separate lineages in the chloroplast phylogeny.     

Adaptive evolution of chloroplast genomes in ancestral grasses

The grass family, Poaceae, is one of the most successful families among angiosperms. Although it has long been suggested that the chloroplast genomes of the Poaceae have undergone an elevated evolutionary rate compared to other angiosperms, little was known about the details of this phenomenon. By using chloroplast genome data from 31 seed plants species, we recently showed that episodic rate acceleration occurred in the common ancestral branch of the core Poaceae (a clade formed by rice Oryza sativa, wheat Triticum aestivum, maize Zea mays and their allies) accompanied by elevated non-synonymous/synonymous rate ratio, while the rate and the non-synonymous/synonymous rate ratio reverted to the low level typical of most monocot species in the terminal branches. It was further shown that positive selection or adaptive evolution operated in several chloroplast proteins during the evolution of ancestral grasses, and the amino acid sites which putatively experienced positive selection have been identified. These findings illustrate the importance of future works of structural biological research of chloroplasts to understand the background of the evolution of the successful group, Poaceae.Key words: rate acceleration, positive selection, non-synonymous/synonymous rate ratio, Poaceae, structural biologyThe grass family, Poaceae, is one of the largest plant families, comprising about 10,000 species including the most important agricultural plants, rice, wheat and maize, as well as grass-dominated ecosystems which comprise about one-third of Earth''s vegetative cover and support a vast number of animals.¹ The chloroplast genes of the grass family Poaceae are known to have undergone accelerations in their evolutionary rates,^2,3 yet little was known about the details of this acceleration. It has become increasingly feasible to estimate the phylogenetic tree of angiosperms and to clarify the tempo and model of molecular evolution by using chloroplast genome sequences.^4–6 By using chloroplast genome data from 13 monocot species and 18 species from dicots and gymnosperms (31 species in total), we recently examined the details of this phenomenon from several aspects.Figure 1 shows the Poales + Musa part of the chloroplast ML tree of the 31 species, and the elongated branches of Poaceae show the rate acceleration in that particular group. Moreover, longer distances of the Poaceae species from Musa than the Typha/Musa distance by more than two times both in terms of non-synonymous and synonymous substitutions (Fig. 1) indicate that both types of substitutions have undergone rate acceleration along the line leading to Poaceae. To explore the pattern of rate change during the course of grass evolution more in detail, we estimated the time-scale of Angiosperm phylogeny with a relaxed clock based on the Bayesian method implemented in MCMCTREE program of PAML.⁷ As is apparent from Figure 1, the molecular evolutionary rate (substitution rate) differs among different lineages, and therefore we used the relaxed clock method which takes account of the evolution of the evolutionary rate in estimating the divergence times and the pattern of rate change. Based on fossil evidence, we assumed the followings in calibrating the relaxed clock; (1) the Gymnosperm/Angiosperm divergence occurred at 280–310 Ma (million years ago),^5,6 (2) the divergence of Poales from other monocots occurred before 115 Ma,^8,9 and (3) the most basal divergence in eudicots occurred before 125 Ma.⁵Open in a separate windowFigure 1Poales + Musa part of the chloroplast genome tree from 31 seed plant species. The branch lengths are proportional to the estimated lengths by the ML with the codon-substitution model (CODEML in PAML⁷). Non-synonymous (d_N) and synonymous (d_S) distances of Poales from Musa and ω = d_N/d_S (along branches) were estimated by the same program.We further gave a constraint of >65 Ma for the Zea/Oryza divergence based on the recent finding of 65 Ma grass phytoliths in dinosaur coprolites which places the diversification of the grasses to the Cretaceous period.^10,11 As a result, it was found that the rate acceleration was limited to the common ancestral branch of Poaceae after they diverged from Musa and that the rate reverted to the slow rate typical of most monocot species in the terminal branches. Even when the constraint was removed, almost the same pattern of rate change was obtained, suggesting that our conclusion regarding accelerated rate in the ancestral grasses followed by the reverted slow rate in contemporary Poaceae is robust.Non-synonymous/synonymous rate ratio (ω = d_N/d_S) is widely used as an indicator of positive selection or adaptive evolution.¹² Figure 1 also indicates a pronounced increase of ω ratio in the common ancestral branch of Poaceae after their divergence from Typha, followed by reversion in the terminal branches to the lower level typical of basal lineages. The elevation of the ω ratio can be due either by adaptive evolution or by relaxation of selective constraints. An ω value higher than 1 is usually regarded as an evidence of adaptive evolution, but since the ω values shown in the figure averages over the entire protein-encoding genes, we would not obtain such a high value even if positive selection operated in some regions of some proteins. To identify positively selected sites, among 75 protein-encoding genes, we at first selected 14 genes, for which the model with higher ω in the ancestral grass branch than others is significantly better than the model with homogeneous ω, and by using the branch-site model,^13,14 we identified 5 genes (atpE, cemA, clpP, rpoB and rps11) which have p value of the branch-site likelihood ratio test less than 0.05 and contain positively selected sites. The amino acid sites and substitutions identified to have experienced positive selection are as follows; atpE (2T→K, 17S→C, 41A→N, 64M→W, 132V→W), cemA (55N→R, 76Y→K, 161W→F, 190I→F, 204I→C), clpP (26R→V, 48V→T, 86F→T, 112I→P, 134E→R, 182T→D), rpoB (90R→F, 338G→K, 1026G→N), rps11 (54V→P, 62A→S, 82A→R, 105L→S, 115R→A, 120L→R) where the numberings of amino acid sites are those of Zea mays.¹⁵ We anticipate that these amino acid substitutions might have relevance to the successful evolution of grasses. To clarify the implication of these findings, structural biological studies of chloroplast proteins on how the amino acid changes affect their functions are needed.Rates of molecular evolution can be potentially linked to life history of organisms. By comparing evolutionary rates of chloroplast, nuclear and mitochondrial genes across five groups of angiosperms, Smith and Donoghue¹⁶ found that the rates are generally low in trees/shrubs compared to related herbs. Our finding, however, suggests that the pattern of rate change during evolution is more complicated than has previously been anticipated, and highlights the need for distinguishing rates of internal branches and those of terminal branches rather than averaging along a lineage in addressing this complicated problem.As Theodosius Dobzhansky¹⁷ wrote, nothing in biology makes sense except in the light of evolution, and the functional background of the molecular machinery in chloroplasts should be interpreted in the light of evolution. We hope our molecular evolutionary analysis of the chloroplast genomes is the first step towards this goal, and hope collaboration of molecular evolutionists with structural biologists becomes fruitful in the future.     

Beta chloroplast genomes: analysis of Fraction I protein and chloroplast DNA variation

Summary The interrelationships of Beta chloroplast genomes have been investigated on the basis of the analysis of Fraction I protein and chloroplast (ct) DNA. Three groups of the chloroplast genomes could be demonstrated by the difference in isoelectric points of the large subunit of Fraction I protein. Restriction enzyme analysis revealed inter- and intra-specific variations among the ctDNAs, which enabled us to detect seven distinct ctDNA types. In Vulgares and Corollinae species, the observed differences were physically mapped taking advantage of the restriction fragment map available for sugar beet (B. vulgaris) ctDNA. The DNA variations were found to result either from gains or losses of restriction sites or from small deletions/ insertions, and most of them were located in the large single-copy region of the genome. Moreover, the ctDNAs from Patellares species are more diverged from those of other Beta taxa. Our results also indicate that there is a close correlation between the chloroplast genome diversity and the accepted taxonomic classification of the species included in this survey.     

Complete chloroplast genomes of Liliaceae (s.l.) species: comparative genomic and phylogenetic analyses

We first report the complete chloroplast (cp) genome of Fritillaria taipaiensis and determine its characteristics, sequence divergence and phylogenetic relationships by comparing it with complete cp genomes of Liliaceae s.l. (including e.g. Nartheciaceae, Amaryllidaceae and Asparagaceae) species obtained from NCBI Genbank. We show that the ycf1, ycf15 and infA genes have become pseudogenes or are lost in some of the seventeen Liliaceae species, and that dispersed repeats are prevailing among the four types of repeats (dispersed, palindromic, complement and tandem repeats). The number of simple sequence repeats ranged from 53 to 84 in the seventeen species, with mononucleotide repeats being the most abundant, followed by dinucleotides. A total of nine genes with positive selection sites were identified (atpB, atpE, ndhF, ndhH, petB, rpl2, rpl20, rpl22 and ycf2). Furthermore, we examined 19 mutational hotspot regions, including three coding regions (rps16, infA and rpl22) and sixteen non-coding regions. A phylogenetic analysis of the complete cp genomes and protein-coding sequences showed that Fritillaria is most closely related to Lilium. Moreover, Asparagus and Polygonatum, Hosta and Yucca are closely related to the Liliaceae. These results will contribute to further study of evolutionary patterns and phylogenetic relationships in Liliaceae s.l.     

8. Biogeography of Vaucheria species from European freshwater/soil habitats: implications from chloroplast genomes

Infraspecific heterogeneity of chloroplast genomes was found in four Vaucheria species (V. bursata, V. cruciata, V. geminata, V. prolifera) collected from six European countries. The degree of sequence variability among strains of each of the four species, as demonstrated by restriction site analysis, exceeds that of higher plant species or even genera. Mainly single base substitutions and, to a much lesser extent, minor insertions/deletions account for such differences, whereas the linear gene arrangement remains unaffected. Chloroplast genotypes found to be identical among strains collected from different geographical localities are considered the common genotype of a given species. These findings are discussed with respect to evolution, biogeographical distribution and the species concept of this genus.