A chloroplast DNA mutational hotspot and gene conversion in a noncoding region near rbcL in the grass family (Poaceae)

Complete Chloroplast Genomes of 14 Subspecies of D. glomerata: Phylogenetic and Comparative Genomic Analyses

Orchardgrass (Dactylis glomerata L.) is a species in the Gramineae family that is highly important economically and valued for its role in ecology. However, the phylogeny and taxonomy of D. glomerata are still controversial based on current morphological and molecular evidence. The study of chloroplast (cp) genomes has developed into a powerful tool to develop molecular markers for related species and reveal the relationships between plant evolution and phylogenetics. In this study, we conducted comparative genomic analyses and phylogenetic inferences on 14 cp genomes of D. glomerata originating from the Mediterranean and Eurasia. The genome size ranged from 134,375 bp to 134,993 bp and exhibited synteny of gene organization and order. A total of 129–131 genes were identified, including 85–87 protein coding genes, 38 tRNA genes and 8 rRNA genes. The cp sequences were highly conserved, and key sequence variations were detected at the junctions of inverted repeats (IRs)/small single–copy (SSC) regions. Moreover, nine highly variable regions were identified among the subspecies based on a sequence divergence analysis. A total of 285 RNA editing sites were detected that were relevant to 52 genes, where rpoB exhibited the most abundant RNA editing sites. The phylogenetic analysis revealed that all Dactylis subspecies clustered into a monophyletic group and most branches provided a high support bootstrap. The main divergence time of D. glomerata was dated to the Miocene era, and this could have been due to changes in the climate. These findings will provide useful insights for further studies on phylogeny, the identification of subspecies and the development of hypotheses for the evolutionary history of the genus Dactylis and of the Gramineae family.

Comparative analysis of two Korean irises (Iris ruthenica and I. uniflora, Iridaceae) based on plastome sequencing and micromorphology

Iris ruthenica Ker Gawl. and I. uniflora Pall. ex Link, which are rare and endangered species in Korea, possess considerable horticultural and medicinal value among Korean irises. However, discrimination of the species is hindered by extensive morphological similarity. Thus, the aim of the present study was to identify discriminating features by comparing the species’ complete plastid genome (i.e., plastome) sequences and micromorphological features, including leaf margins, stomatal complex distribution (hypostomatic vs. amphistomatic leaves), anther stomata density, and tepal epidermal cell patterns. Plastome comparison revealed slightly divergent regions within intergenic spacer regions, and the most variable sequences, which were distributed in non-coding regions, could be used as molecular markers for the discrimination of I. ruthenica and I. uniflora. Phylogenetic analysis of the Iris species revealed that I. ruthenica and I. uniflora formed a well-supported clade. The comparison of plastomes and micromorphological features performed in this study provides useful information for elucidating taxonomic, phylogenetic, and evolutionary relationships in Iridaceae. Further studies, including those based on molecular cytogenetic approaches using species specific markers, will offer insights into species delimitation of the two closely related Iris species.

Context-Dependent Mutation Dynamics, Not Selection, Explains the Codon Usage Bias of Most Angiosperm Chloroplast Genes

Two competing proposals about the degree to which selection affects codon usage of angiosperm chloroplast genes are examined. The first, based on observations that codon usage does not match expectations under the naïve assumption that base composition will be identical at all neutral sites, is that selection plays a significant role. The second is that codon usage is determined almost solely by mutation bias and drift, with selection influencing only one or two highly expressed genes, in particular psbA. First it is shown that, as a result of an influence of neighboring base composition on mutation dynamics, compositional biases are expected to be widely divergent at different sites in the absence of selection. The observed mutation properties are then used to predict expected neutral codon usage biases and to show that observed deviations from the naïve expectations are in fact expected given the context-dependent mutational dynamics. It is also shown that there is a match between the observed and expected codon usage when context effects are taken into consideration, with psbA being a notable exception. Overall, the data support the model that selection is not a widespread factor affecting the codon usage of angiosperm chloroplast genes and highlight the need to have an accurate model of mutational dynamics.

Plastid Phylogenomic Data Offers Novel Insights Into the Taxonomic Status of the Trichosanthes kirilowii Complex (Cucurbitaceae) in South Korea

Trichosanthes is a genus in Cucurbitaceae comprising 90-100 species. Trichosanthes species are valuable as herbaceous medicinal ingredients. The fruits, seeds, and roots of species such as T. kirilowii and T. rosthornii are used in Korean traditional herbal medicines. T. rosthornii is only found in China, whereas in South Korea two varieties, T. kirilowii var. kirilowii and T. kirilowii var. japonica, are distributed. T. kirilowii var. kirilowii and T. kirilowii var. japonica have different fruit and leaf shapes but are recognized as belonging to the same species. Furthermore, although its members have herbal medicine applications, genomic information of the genus is still limited. The broad goals of this study were (i) to evaluate the taxonomy of Trichosanthes using plastid phylogenomic data and (ii) provide molecular markers specific for T. kirilowii var. kirilowii and T. kirilowii var. japonica, as these have differences in their pharmacological effectiveness and thus should not be confused and adulterated. Comparison of five Trichosanthes plastid genomes revealed locally divergent regions, mainly within intergenic spacer regions (trnT-UGU-trnL-UAA: marker name Tri, rrn4.5-rrn5: TRr, trnE-UUC-trnT-GGU: TRtt). Using these three markers as DNA-barcodes for important herbal medicine species in Trichosanthes, the identity of Trichosanthes material in commercial medicinal products in South Korea could be successfully determined. Phylogenetic analysis of the five Trichosanthes species revealed that the species are clustered within tribe Sicyoeae. T. kirilowii var. kirilowii and T. rosthornii formed a clade with T. kirilowii var. japonica as their sister group. As T. kirilowii in its current circumscription is paraphyletic and as the two varieties can be readily distinguished morphologically (e.g., in leaf shape), T. kirilowii var. japonica should be treated (again) as an independent species, T. japonica.

Gene Loss, Pseudogenization in Plastomes of Genus Allium (Amaryllidaceae), and Putative Selection for Adaptation to Environmental Conditions

Amaryllidaceae is a large family with more than 1,600 species, belonging to 75 genera. The largest genus—Allium—is vast, comprising about a thousand species. Allium species (as well as other members of the Amaryllidaceae) are widespread and diversified, they are adapted to a wide range of habitats from shady forests to open habitats like meadows, steppes, and deserts. The genes present in chloroplast genomes (plastomes) play fundamental roles for the photosynthetic plants. Plastome traits could thus be associated with geophysical abiotic characteristics of habitats. Most chloroplast genes are highly conserved and are used as phylogenetic markers for many families of vascular plants. Nevertheless, some studies revealed signatures of positive selection in chloroplast genes of many plant families including Amaryllidaceae. We have sequenced plastomes of the following nine Allium (tribe Allieae of Allioideae) species: A. zebdanense, A. moly, A. victorialis, A. macleanii, A. nutans, A. obliquum, A. schoenoprasum, A. pskemense, A. platyspathum, A. fistulosum, A. semenovii, and Nothoscordum bivalve (tribe Leucocoryneae of Allioideae). We compared our data with previously published plastomes and provided our interpretation of Allium plastome genes’ annotations because we found some noteworthy inconsistencies with annotations previously reported. For Allium species we estimated the integral evolutionary rate, counted SNPs and indels per nucleotide position as well as compared pseudogenization events in species of three main phylogenetic lines of genus Allium to estimate whether they are potentially important for plant physiology or just follow the phylogenetic pattern. During examination of the 38 species of Allium and the 11 of other Amaryllidaceae species we found that rps16, rps2, infA, ccsA genes have lost their functionality multiple times in different species (regularly evolutionary events), while the pseudogenization of other genes was stochastic events. We found that the “normal” or “pseudo” state of rps16, rps2, infA, ccsA genes correlates well with the evolutionary line of genus the species belongs to. The positive selection in various NADH dehydrogenase (ndh) genes as well as in matK, accD, and some others were found. Taking into account known mechanisms of coping with excessive light by cyclic electron transport, we can hypothesize that adaptive evolution in genes, coding subunits of NADH-plastoquinone oxidoreductase could be driven by abiotic factors of alpine habitats, especially by intensive light and UV radiation.

Plastid Genomes of Flowering Plants: Essential Principles

The plastid genome (plastome ) has proved a valuable source of data for evaluating evolutionary relationships among angiosperms. Through basic and applied approaches, plastid transformation technology offers the potential to understand and improve plant productivity, providing food, fiber, energy, and medicines to meet the needs of a burgeoning global population. The growing genomic resources available to both phylogenetic and biotechnological investigations is allowing novel insights and expanding the scope of plastome research to encompass new species. In this chapter, we present an overview of some of the seminal and contemporary research that has contributed to our current understanding of plastome evolution and attempt to highlight the relationship between evolutionary mechanisms and the tools of plastid genetic engineering.

SPInDel Analysis of the Non-Coding Regions of cpDNA as a More Useful Tool for the Identification of Rye (Poaceae: Secale) Species

Secale is a small but very diverse genus from the tribe Triticeae (family Poaceae), which includes annual, perennial, self-pollinating and open-pollinating, cultivated, weedy and wild species of various phenotypes. Despite its high economic importance, classification of this genus, comprising 3–8 species, is inconsistent. This has resulted in significantly reduced progress in the breeding of rye which could be enriched with functional traits derived from wild rye species. Our previous research has suggested the utility of non-coding sequences of chloroplast and mitochondrial DNA in studies on closely related species of the genus Secale. Here we applied the SPInDel (Species Identification by Insertions/Deletions) approach, which targets hypervariable genomic regions containing multiple insertions/deletions (indels) and exhibiting extensive length variability. We analysed a total of 140 and 210 non-coding sequences from cpDNA and mtDNA, respectively. The resulting data highlight regions which may represent useful molecular markers with respect to closely related species of the genus Secale, however, we found the chloroplast genome to be more informative. These molecular markers include non-coding regions of chloroplast DNA: atpB-rbcL and trnT-trnL and non-coding regions of mitochondrial DNA: nad1B-nad1C and rrn5/rrn18. Our results demonstrate the utility of the SPInDel concept for the characterisation of Secale species.

Dissection for Floral Micromorphology and Plastid Genome of Valuable Medicinal Borages Arnebia and Lithospermum (Boraginaceae)

The genera Arnebia and Lithospermum (Lithospermeae-Boraginaceae) comprise 25-30 and 50-60 species, respectively. Some of them are economically valuable, as their roots frequently contain a purple-red dye used in the cosmetic industry. Furthermore, dried roots of Arnebia euchroma, A. guttata, and Lithospermum erythrorhizon, which have been designated Lithospermi Radix, are used as traditional Korean herbal medicine. This study is the first report on the floral micromorphology and complete chloroplast (cp) genome sequences of A. guttata (including A. tibetana), A. euchroma, and L. erythrorhizon. We reveal great diversity in floral epidermal cell patterns, gynoecium, and structure of trichomes. The cp genomes were 149,361-150,465 bp in length, with conserved quadripartite structures. In total, 112 genes were identified, including 78 protein-coding regions, 30 tRNA genes, and four rRNA genes. Gene order, content, and orientation were highly conserved and were consistent with the general structure of angiosperm cp genomes. Comparison of the four cp genomes revealed locally divergent regions, mainly within intergenic spacer regions (atpH-atpI, petN-psbM, rbcL-psaI, ycf4-cemA, ndhF-rpl32, and ndhC-trnV-UAC). To facilitate species identification, we developed molecular markers psaA- ycf3 (PSY), trnI-CAU- ycf2 (TCY), and ndhC-trnV-UAC (NCTV) based on divergence hotspots. High-resolution phylogenetic analysis revealed clear clustering and a close relationship of Arnebia to its Lithospermum sister group, which was supported by strong bootstrap values and posterior probabilities. Overall, gynoecium characteristics and genetic distance of cp genomes suggest that A. tibetana, might be recognized as an independent species rather than a synonym of A. guttata. The present morphological and cp genomic results provide useful information for future studies, such as taxonomic, phylogenetic, and evolutionary analysis of Boraginaceae.

Complete plastome sequencing of Allium paradoxum reveals unusual rearrangements and the loss of the ndh genes as compared to Allium ursinum and other onions

In this work the complete chloroplast DNAs of Allium paradoxum and Allium ursinum, two edible species of Allium subg. Amerallium (the first lineage), were sequenced, assembled, annotated, and compared with complete Allium plastomes of the second and third evolutionary lines from GenBank database. The A. ursinum plastome contains 90 predicted genes (81 unique) including 5 pseudogenes, while A. paradoxum has 88 predicted genes (79 unique) including 19 pseudogenes. The comparative analysis has revealed that the A. paradoxum plastome differs markedly from those of other species. Due to many deletions, the A. paradoxum plastome is the shortest of known for Allium species, being only 145,819 bp long. The most prominent distinctions are (1) a 4825 bp long local inversion that spans from the ndhE to the rpl32 gene in the small single copy region and (2) pseudogenization, or the loss of all NADH-genes. In contrast, the plastome of A. ursinum - a species from the first evolutionary line (as well as A. paradoxum) - resembles the Allium species of the second and third evolutionary lines, showing no large rearrangements or discrepancies in gene content. It is unclear yet whether only A. paradoxum was affected by some evolutionary events or its close relatives from both sect. Briseis and other sections of Amerallium were altered as well. We speculate the sunlight-intolerant, shade-loving nature of A. paradoxum and the impairment of the ndh genes in its plastome could be interrelated phenomena.

The rpl23 gene and pseudogene are hotspots of illegitimate recombination in barley chloroplast mutator seedlings

Previously, through a TILLING (Targeting Induced Local Lesions in Genomes) approach applied on barley chloroplast mutator (cpm) seedlings a high frequency of polymorphisms in the rpl23 gene was detected. All the polymorphisms corresponded to five differences already known to exist in nature between the rpl23 gene located in the inverted repeats (IRs) and the rpl23 pseudogene located in the large single copy region (LSC). In this investigation, polymorphisms in the rpl23 gene were verified and besides, a similar situation was found for the pseudogene in cpm seedlings. On the other hand, no polymorphisms were found in any of those loci in 40 wild type barley seedlings. Those facts and the independent occurrence of polymorphisms in the gene and pseudogene in individual seedlings suggest that the detected polymorphisms initially arose from gene conversion between gene and pseudogene. Moreover, an additional recombination process involving small recombinant segments seems to occur between the two gene copies as a consequence of their location in the IRs. These and previous results support the hypothesis that the CPM protein is a component of the plastome mismatch repair (MMR) system, whose failure of the anti-recombination activity results in increased illegitimate recombination between the rpl23 gene and pseudogene.

Genetic diversity and relationship between cultivated, weedy and wild rye species as revealed by chloroplast and mitochondrial DNA non-coding regions analysis

The genus Secale is small but very diverse. Despite the high economic importance, phylogenetic relationships of rye species have not been fully determined, and they are extremely important for the process of breeding of new cultivars that can be enriched with functional traits derived from wild rye species. The study analyzed the degree of relationship of 35 accessions of the genus Secale, representing 13 most often distinguished species and subspecies, originating from various seed collections in the world, based on the analysis of non-coding regions of the chloroplast (cpDNA) and mitochondrial genome (mtDNA), widely used in phylogenetic and population plant studies, because of a higher rate of evolution than the coding regions. There was no clear genetic structure between different species and subspecies, which may indicated the introgression between these taxa. The obtained data confirmed that S. vavilovii was very similar to S. cereale, which confirmed the assumption that they might share a common ancestor. The results also confirmed the divergence of S. sylvestre from other species and subspecies of rye. Areas that may be useful molecular markers in studies on closely related species of the genus Secale were also indicated.

Plastid super-barcodes as a tool for species discrimination in feather grasses (Poaceae: Stipa)

Present study was designed to verify which or if any of plastome loci is a hotspot region for mutations and hence might be useful for molecular species identification in feather grasses. 21 newly sequenced complete plastid genomes representing 19 taxa from the genus of Stipa were analyzed in search of the most variable and the most discriminative loci within Stipa. The results showed that the problem with selecting a good barcode locus for feather grasses lies in the very low level of genetic diversity within its plastome. None of the single chloroplast loci is polymorphic enough to play a role of a barcode or a phylogenetic marker for Stipa. The biggest number of taxa was successfully identified by the analysis of 600 bp long DNA fragment comprising a part of rbcL gene, the complete rbcL-rpl23 spacer and a part of rpl23 gene. The effectiveness of multi-locus barcode composed of six best-performing loci for Stipa (ndhH, rpl23, ndhF-rpl32, rpl32-ccsA, psbK-psbI and petA-psbJ) didn't reach 70% of analyzed taxa. The analysis of complete plastome sequences as a super-barcode for Stipa although much more effective, still didn't allow for discrimination of all the analyzed taxa of feather grasses.

Acetyl-CoA carboxylase overexpression in herbicide-resistant large crabgrass (Digitaria sanguinalis)

BACKGROUND

The occurrence of herbicide-resistant weed biotypes is increasing and this report of an acetyl-CoA carboxylase (ACCase) inhibitor-resistant Digitaria sanguinalis L. Scop. from southwestern Ontario is another example. The identified weed escaped control in an onion and carrot rotation in which graminicides were used for several consecutive years. Our goal was to characterize the level and mechanism of resistance of the biotype.

RESULTS

The biotype was resistant to all five ACCase inhibitor herbicides tested. Gene-expression profiling was performed because none of the mutations known to confer resistance in the ACCase gene were detected. RNASeq and quantitative reverse-transcriptase PCR (qRT-PCR) results indicated that transcription of ACCase was 3.4-9.3 times higher in the resistant biotype than the susceptible biotype. ACCase gene copy number was determined by qPCR to be five to seven times higher in the resistant compared with the susceptible biotype. ACCase gene overexpression was directly related to the increase of the ACCase gene copy number.

CONCLUSION

Our results are consistent with the hypothesis that overexpression of the herbicide target gene ACCase confers resistance to the herbicide. This is the first reported case of target gene duplication conferring resistance to a herbicide other than glyphosate. © 2017 Society of Chemical Industry See related Article.

Mitogenomic analyses support the recent division of the genus Orthotrichum (Orthotrichaceae, Bryophyta)

A recently presented taxonomical arrangement of the moss genus Orthotrichum Hedw. s.l. substantially changed the traditional view of the taxon that had been accepted throughout the twentieth century. This paper provides the results of mitogenomic studies that strongly support the new taxonomical concept. Comparative analyses presented in this study confirmed the stable structure of moss mitogenomes. Moreover, 17 complete mitogenome sequences were used to identify the major evolutionary groups, including 11 newly sequenced ones, for this study. The analysis of mitochondrial hotspots revealed intron 4 of the cox1 gene to be the most variable non-coding region. The most variable protein-coding genes in the tribe Orthotricheae were ccmFC and tatC. The intergenic and intronic hotspots of Orthotrichum s.l. identified in the present study do not correspond to those described in vascular plant mitogenomes.

Recombination-dependent replication and gene conversion homogenize repeat sequences and diversify plastid genome structure

PREMISE OF THE STUDY

There is a misinterpretation in the literature regarding the variable orientation of the small single copy region of plastid genomes (plastomes). The common phenomenon of small and large single copy inversion, hypothesized to occur through intramolecular recombination between inverted repeats (IR) in a circular, single unit-genome, in fact, more likely occurs through recombination-dependent replication (RDR) of linear plastome templates. If RDR can be primed through both intra- and intermolecular recombination, then this mechanism could not only create inversion isomers of so-called single copy regions, but also an array of alternative sequence arrangements.

METHODS

We used Illumina paired-end and PacBio single-molecule real-time (SMRT) sequences to characterize repeat structure in the plastome of Monsonia emarginata (Geraniaceae). We used OrgConv and inspected nucleotide alignments to infer ancestral nucleotides and identify gene conversion among repeats and mapped long (>1 kb) SMRT reads against the unit-genome assembly to identify alternative sequence arrangements.

RESULTS

Although M. emarginata lacks the canonical IR, we found that large repeats (>1 kilobase; kb) represent ∼22% of the plastome nucleotide content. Among the largest repeats (>2 kb), we identified GC-biased gene conversion and mapping filtered, long SMRT reads to the M. emarginata unit-genome assembly revealed alternative, substoichiometric sequence arrangements.

CONCLUSION

We offer a model based on RDR and gene conversion between long repeated sequences in the M. emarginata plastome and provide support that both intra-and intermolecular recombination between large repeats, particularly in repeat-rich plastomes, varies unit-genome structure while homogenizing the nucleotide sequence of repeats.

Chloroplast Genome Analysis of Resurrection Tertiary Relict Haberlea rhodopensis Highlights Genes Important for Desiccation Stress Response

Haberlea rhodopensis is a paleolithic tertiary relict species, best known as a resurrection plant with remarkable tolerance to desiccation. When exposed to severe drought stress, H. rhodopensis shows an ability to maintain the structural integrity of its photosynthetic apparatus, which re-activates easily upon rehydration. We present here the results from the assembly and annotation of the chloroplast (cp) genome of H. rhodopensis, which was further subjected to comparative analysis with the cp genomes of closely related species. H. rhodopensis showed a cp genome size of 153,099 bp, harboring a pair of inverted repeats (IR) of 25,415 bp separated by small and large copy regions (SSC and LSC) of 17,826 and 84,443 bp. The genome structure, gene order, GC content and codon usage are similar to those of the typical angiosperm cp genomes. The genome hosts 137 genes representing 70.66% of the plastome, which includes 86 protein-coding genes, 36 tRNAs, and 4 rRNAs. A comparative plastome analysis with other closely related Lamiales members revealed conserved gene order in the IR and LSC/SSC regions. A phylogenetic analysis based on protein-coding genes from 33 species defines this species as belonging to the Gesneriaceae family. From an evolutionary point of view, a site-specific selection analysis detected positively selected sites in 17 genes, most of which are involved in photosynthesis (e.g., rbcL, ndhF, accD, atpE, etc.). The observed codon substitutions may be interpreted as being a consequence of molecular adaptation to drought stress, which ensures an evolutionary advantage to H. rhodopensis.

Plastome Mutations and Recombination Events in Barley Chloroplast Mutator Seedlings

The barley chloroplast mutator (cpm) is an allele of a nuclear gene that when homozygous induces several types of cytoplasmically inherited chlorophyll deficiencies. In this work, a plastome Targeting Induced Local Lesions in Genomes (TILLING) strategy based on mismatch digestion was used on families that carried the cpm genotype through many generations. Extensive scanning of 33 plastome genes and a few intergenic regions was conducted. Numerous polymorphisms were detected on both genic and intergenic regions. The detected polymorphisms can be accounted for by at least 61 independent mutational events. The vast majority of the polymorphisms originated in substitutions and small indels (insertions/deletions) in microsatellites. The rpl23 and the rps16 genes were the most polymorphic. Interestingly, the variation observed in the rpl23 gene consisted of several combinations of 5 different one nucleotide polymorphisms. Besides, 4 large indels that have direct repeats at both ends were also observed, which appear to be originated from recombinational events. The cpm mutation spectrum suggests that the CPM gene product is probably involved in plastome mismatch repair. The numerous subtle molecular changes that were localized in a wide range of plastome sites show the cpm as a valuable source of plastome variability for plant research and/or plant breeding. Moreover, the cpm mutant appears to be an interesting experimental material for investigating the mechanisms responsible for maintaining the stability of plant organelle DNA.

Genomic Resources of Three Pulsatilla Species Reveal Evolutionary Hotspots, Species-Specific Sites and Variable Plastid Structure in the Family Ranunculaceae

Background: The European continent is presently colonized by nine species of the genus Pulsatilla, five of which are encountered only in mountainous regions of southwest and south-central Europe. The remaining four species inhabit lowlands in the north-central and eastern parts of the continent. Most plants of the genus Pulsatilla are rare and endangered, which is why most research efforts focused on their biology, ecology and hybridization. The objective of this study was to develop genomic resources, including complete plastid genomes and nuclear rRNA clusters, for three sympatric Pulsatilla species that are most commonly found in Central Europe. The results will supply valuable information about genetic variation, which can be used in the process of designing primers for population studies and conservation genetics research. The complete plastid genomes together with the nuclear rRNA cluster can serve as a useful tool in hybridization studies. Methodology/principal findings: Six complete plastid genomes and nuclear rRNA clusters were sequenced from three species of Pulsatilla using the Illumina sequencing technology. Four junctions between single copy regions and inverted repeats and junctions between the identified locally-collinear blocks (LCB) were confirmed by Sanger sequencing. Pulsatilla genomes of 120 unique genes had a total length of approximately 161–162 kb, and 21 were duplicated in the inverted repeats (IR) region. Comparative plastid genomes of newly-sequenced Pulsatilla and the previously-identified plastomes of Aconitum and Ranunculus species belonging to the family Ranunculaceae revealed several variations in the structure of the genome, but the gene content remained constant. The nuclear rRNA cluster (18S-ITS1-5.8S-ITS2-26S) of studied Pulsatilla species is 5795 bp long. Among five analyzed regions of the rRNA cluster, only Internal Transcribed Spacer 2 (ITS2) enabled the molecular delimitation of closely-related Pulsatilla patens and Pulsatillavernalis. Conclusions/significance: The determination of complete plastid genome and nuclear rRNA cluster sequences in three species of the genus Pulsatilla is an important contribution to our knowledge of the evolution and phylogeography of those endangered taxa. The resulting data can be used to identify regions that are particularly useful for barcoding, phylogenetic and phylogeographic studies. The investigated taxa can be identified at each stage of development based on their species-specific SNPs. The nuclear and plastid genomic resources enable advanced studies on hybridization, including identification of parent species, including their roles in that process. The identified nonsynonymous mutations could play an important role in adaptations to changing environments. The results of the study will also provide valuable information about the evolution of the plastome structure in the family Ranunculaceae.

The plastid genomes of flowering plants

The plastid genome (plastome) has proved a valuable source of data for evaluating evolutionary relationships among angiosperms. Through basic and applied approaches, plastid transformation technology offers the potential to understand and improve plant productivity, providing food, fiber, energy and medicines to meet the needs of a burgeoning global population. The growing genomic resources available to both phylogenetic and biotechnological investigations are allowing novel insights and expanding the scope of plastome research to encompass new species. In this chapter we present an overview of some of the seminal and contemporary research that has contributed to our current understanding of plastome evolution and attempt to highlight the relationship between evolutionary mechanisms and tools of plastid genetic engineering.

The complete chloroplast DNA sequence of Eleutherococcus senticosus (Araliaceae); comparative evolutionary analyses with other three asterids

This study reports the complete chloroplast (cp) DNA sequence of Eleutherococcus senticosus (GenBank: JN 637765), an endangered endemic species. The genome is 156,768 bp in length, and contains a pair of inverted repeat (IR) regions of 25,930 bp each, a large single copy (LSC) region of 86,755 bp and a small single copy (SSC) region of 18,153 bp. The structural organization, gene and intron contents, gene order, AT content, codon usage, and transcription units of the E. senticosus chloroplast genome are similar to that of typical land plant cp DNA. We aligned and analyzed the sequences of 86 coding genes, 19 introns and 113 intergenic spacers (IGS) in three different taxonomic hierarchies; Eleutherococcus vs. Panax, Eleutherococcus vs. Daucus, and Eleutherococcus vs. Nicotiana. The distribution of indels, the number of polymorphic sites and nucleotide diversity indicate that positional constraint is more important than functional constraint for the evolution of cp genome sequences in Asterids. For example, the intron sequences in the LSC region exhibited base substitution rates 5-11-times higher than that of the IR regions, while the intron sequences in the SSC region evolved 7-14-times faster than those in the IR region. Furthermore, the Ka/Ks ratio of the gene coding sequences supports a stronger evolutionary constraint in the IR region than in the LSC or SSC regions. Therefore, our data suggest that selective sweeps by base collection mechanisms more frequently eliminate polymorphisms in the IR region than in other regions. Chloroplast genome regions that have high levels of base substitutions also show higher incidences of indels. Thirty-five simple sequence repeat (SSR) loci were identified in the Eleutherococcus chloroplast genome. Of these, 27 are homopolymers, while six are di-polymers and two are tri-polymers. In addition to the SSR loci, we also identified 18 medium size repeat units ranging from 22 to 79 bp, 11 of which are distributed in the IGS or intron regions. These medium size repeats may contribute to developing a cp genome-specific gene introduction vector because the region may use for specific recombination sites.

The evolution of chloroplast genome structure in ferns

The plastid genome (plastome) is a rich source of phylogenetic and other comparative data in plants. Most land plants possess a plastome of similar structure. However, in a major group of plants, the ferns, a unique plastome structure has evolved. The gene order in ferns has been explained by a series of genomic inversions relative to the plastome organization of seed plants. Here, we examine for the first time the structure of the plastome across fern phylogeny. We used a PCR-based strategy to map and partially sequence plastomes. We found that a pair of partially overlapping inversions in the region of the inverted repeat occurred in the common ancestor of most ferns. However, the ancestral (seed plant) structure is still found in early diverging branches leading to the osmundoid and filmy fern lineages. We found that a second pair of overlapping inversions occurred on a branch leading to the core leptosporangiates. We also found that the unique placement of the gene matK in ferns (lacking a flanking intron) is not a result of a large-scale inversion, as previously thought. This is because the intron loss maps to an earlier point on the phylogeny than the nearby inversion. We speculate on why inversions may occur in pairs and what this may mean for the dynamics of plastome evolution.

The chloroplast genome of Anomochloa marantoidea (Anomochlooideae; Poaceae) comprises a mixture of grass-like and unique features

Features in the complete plastome of Anomochloa marantoidea (Poaceae) were investigated. This species is one of four of Anomochlooideae, the crown node of which diverged before those of any other grass subfamily. The plastome was sequenced from overlapping amplicons using previously designed primers. The plastome of A. marantoidea is 138412 bp long with a typical gene content for Poaceae. Five regions were examined in detail because of prior surveys that identified structural alterations among graminoid Poales. Anomochloa marantoidea was found to have an intron in rpoC1, unlike other Poaceae. The insertion region of rpoC2 is unusually short in A. marantoidea compared with those of other grasses, but with atypically long subrepeats. Both ycf1 and ycf2 are nonfunctional as is typical in grasses, but A. marantoidea has a uniquely long ψycf1. Finally, the rbcL-psaI spacer in A. marantoidea is atypically short with no evidence of the ψrpl23 locus found in all other Poaceae. Some of these features are of noteworthy dissimilarity between A. marantoidea and those crown grasses for which entire plastomes have been sequenced. Complete plastome sequences of other Anomochlooideae and outgroups will further advance our understanding of the evolutionary events in the plastome that accompanied graminoid diversification.

Molecular analysis of Stemona plants in China based on sequences of four chloroplast DNA regions

Stemona sessilifolia, S. japonica and S. tuberosa are the three original sources of Stemonae Radix specified in the Chinese Pharmacopoeia (CP), and have been traditionally used for antitussive and insecticidal remedy. Significant variations in alkaloids composition and content, as well as different degrees of antitussive activities were found among them. In order to identify the genuine sources of Stemonae Radix accurately in genetic level, we determined the nucleotide sequences of chloroplast DNA trnL-trnF, trnH-psbA, petB-petD and trnK-rps16 regions of the species recorded in CP and S. parviflora, as well as the common counterfeits of Stemonae Radix, Asparagus species. The results revealed that the sequences of petB-petD and trnK-rps16 regions, showing relatively high substitution rate, were more informative than those of trnL-trnF and trnH-psbA regions. The sequences from all the four regions provided useful information to discriminate the three CP species from each other and from S. parviflora and the counterfeits. A phylogenetic tree reconstructed by the trnH-psbA sequences for 9 Stemona species distributed in China and Thailand showed that the three CP species belonged to the same clade, among which S. japonica and S. sessillifolia formed a sister group, showing closer relations to each other than to S. tuberosa. By contrast, S. parviflora was genetically far from the three CP species. Intra-species variations were observed in the three CP species. Especially, in S. tuberosa two types of petB-petD sequence and four types each of trnL-trnF, trnK-rps16 and trnH-psbA sequences resulted in 6 haplotypes; whereas, these differences had no relation with the different chemical types, but seemed to be consistent with geographical distribution.

Extensive rearrangements in the chloroplast genome of Trachelium caeruleum are associated with repeats and tRNA genes

Chloroplast genome organization, gene order, and content are highly conserved among land plants. We sequenced the chloroplast genome of Trachelium caeruleum L. (Campanulaceae), a member of an angiosperm family known for highly rearranged genomes. The total genome size is 162,321 bp, with an inverted repeat (IR) of 27,273 bp, large single-copy (LSC) region of 100,114 bp, and small single-copy (SSC) region of 7,661 bp. The genome encodes 112 different genes, with 17 duplicated in the IR, a tRNA gene (trnI-cau) duplicated once in the LSC region, and a protein-coding gene (psbJ) with two duplicate copies, for a total of 132 putatively intact genes. ndhK may be a pseudogene with internal stop codons, and clpP, ycf1, and ycf2 are so highly diverged that they also may be pseudogenes. ycf15, rpl23, infA, and accD are truncated and likely nonfunctional. The most conspicuous feature of the Trachelium genome is the presence of 18 internally unrearranged blocks of genes inverted or relocated within the genome relative to the ancestral gene order of angiosperm chloroplast genomes. Recombination between repeats or tRNA genes has been suggested as a mechanism of chloroplast genome rearrangements. The Trachelium chloroplast genome shares with Pelargonium and Jasminum both a higher number of repeats and larger repeated sequences in comparison to eight other angiosperm chloroplast genomes, and these are concentrated near rearrangement endpoints. Genes for tRNAs occur at many but not all inversion endpoints, so some combination of repeats and tRNA genes may have mediated these rearrangements.

Relationship between mRNA secondary structure and sequence variability in Chloroplast genes: possible life history implications

BACKGROUND

Synonymous sites are freer to vary because of redundancy in genetic code. Messenger RNA secondary structure restricts this freedom, as revealed by previous findings in mitochondrial genes that mutations at third codon position nucleotides in helices are more selected against than those in loops. This motivated us to explore the constraints imposed by mRNA secondary structure on evolutionary variability at all codon positions in general, in chloroplast systems.

RESULTS

We found that the evolutionary variability and intrinsic secondary structure stability of these sequences share an inverse relationship. Simulations of most likely single nucleotide evolution in Psilotum nudum and Nephroselmis olivacea mRNAs, indicate that helix-forming propensities of mutated mRNAs are greater than those of the natural mRNAs for short sequences and vice-versa for long sequences. Moreover, helix-forming propensity estimated by the percentage of total mRNA in helices increases gradually with mRNA length, saturating beyond 1000 nucleotides. Protection levels of functionally important sites vary across plants and proteins: r-strategists minimize mutation costs in large genes; K-strategists do the opposite.

CONCLUSION

Mrna length presumably predisposes shorter mRNAs to evolve under different constraints than longer mRNAs. The positive correlation between secondary structure protection and functional importance of sites suggests that some sites might be conserved due to packing-protection constraints at the nucleic acid level in addition to protein level constraints. Consequently, nucleic acid secondary structure a priori biases mutations. The converse (exposure of conserved sites) apparently occurs in a smaller number of cases, indicating a different evolutionary adaptive strategy in these plants. The differences between the protection levels of functionally important sites for r- and K-strategists reflect their respective molecular adaptive strategies. These converge with increasing domestication levels of K-strategists, perhaps because domestication increases reproductive output.

Plastidic trnFUUC pseudogenes in North American genus Boechera (Brassicaceae): mechanistic aspects of evolution

The origin and maintenance of a plastidic tandem repeat next to the TRNF (UUC) gene were analyzed in the genus BOECHERA in a phylogenetic context and were compared to published analogous examples that emerged in parallel in the Asteraceae and Juncaceae, respectively. Although we identified some features common to these taxonomic groups with respect to structure and origin of the region, obvious differences were encountered, which argue against a specific mechanism or evolutionary principle underlying the parallel origin and maintenance of the TRNF-tandem repeats in those families. In contrast to the situation in the Asteraceae, no reciprocal recombinant repeat types have been observed in the Brassicaceae. Forty copy types, classified into three groups, were isolated from 103 chloroplast haplotypes of BOECHERA and it was demonstrated that they are composed of four subregions of various origins. We discuss various mutation mechanisms such as DNA replication slippage, and inter- and intrachromosomal recombination which were reported to mediate variation in copy numbers and other types of observed sequence length polymorphism. It is shown that the observed molecular structure of the tandem repeat region did not fully fit the particular patterns expected under a scenario of evolution including any of the known mechanisms. Nevertheless, it appeared that intermolecular unequal crossing-over is most likely the driving force in the evolution of this tandem repeat. However, it remains to be explained, why no reciprocal recombinant copy types have been observed. The reconstructed phylogenetic relationships among copies reflected different evolutionary scenarios as follows: (1) A single and ancient origin of copies pre-dates the radiation of BOECHERA. (2) Parallel expansion and shortening of the tandem repeat within different BOECHERA lineages. (3) Conservation of the first copy, as it was the only one present in all chloroplast haplotypes.

Comparison of whole chloroplast genome sequences to choose noncoding regions for phylogenetic studies in angiosperms: the tortoise and the hare III

Although the chloroplast genome contains many noncoding regions, relatively few have been exploited for interspecific phylogenetic and intraspecific phylogeographic studies. In our recent evaluation of the phylogenetic utility of 21 noncoding chloroplast regions, we found the most widely used noncoding regions are among the least variable, but the more variable regions have rarely been employed. That study led us to conclude that there may be unexplored regions of the chloroplast genome that have even higher relative levels of variability. To explore the potential variability of previously unexplored regions, we compared three pairs of single-copy chloroplast genome sequences in three disparate angiosperm lineages: Atropa vs. Nicotiana (asterids); Lotus vs. Medicago (rosids); and Saccharum vs. Oryza (monocots). These three separate sequence alignments highlighted 13 mutational hotspots that may be more variable than the best regions of our former study. These 13 regions were then selected for a more detailed analysis. Here we show that nine of these newly explored regions (rpl32-trnL((UAG)), trnQ((UUG))-5'rps16, 3'trnV((UAC))-ndhC, ndhF-rpl32, psbD-trnT((GGU)), psbJ-petA, 3'rps16-5'trnK((UUU)), atpI-atpH, and petL-psbE) offer levels of variation better than the best regions identified in our earlier study and are therefore likely to be the best choices for molecular studies at low taxonomic levels.

Elimination of deleterious mutations in plastid genomes by gene conversion

Asexual reproduction is believed to be detrimental, mainly because of the accumulation of deleterious mutations over time, a hypothesis known as Muller's ratchet. In seed plants, most asexually reproducing genetic systems are polyploid, with apomictic species (plants forming seeds without fertilization) as well as plastids and mitochondria providing prominent examples. Whether or not polyploidy helps asexual genetic systems to escape Muller's ratchet is unknown. Gene conversion, particularly when slightly biased, represents a potential mechanism that could allow asexual genetic systems to reduce their mutation load in a genome copy number-dependent manner. However, direct experimental evidence for the operation of gene conversion between genome molecules to correct mutations is largely lacking. Here we describe an experimental system based on transgenic tobacco chloroplasts that allows us to analyze gene conversion events in higher plant plastid genomes. We provide evidence for gene conversion acting as a highly efficient mechanism by which the polyploid plastid genetic system can correct deleterious mutations and make one good genome out of two bad ones. Our finding that gene conversion can be biased may provide a molecular link between asexual reproduction, high genome copy numbers and low mutation rates.

Variations in a hotspot region of chloroplast DNAs among common wheat and Aegilops revealed by nucleotide sequence analysis

The second largest BamHI fragment (B2) of the chloroplast DNA in Triticum (wheat) and Aegilops contains a highly variable region (a hotspot), resulting in four types of B2 of different size, i.e. B2l (10.5kb), B2m (10.2kb), B2 (9.6kb) and B2s (9.4kb). In order to gain a better understanding of the molecular nature of the variations in length and explain unexpected identity among B2 of Ae. ovata, Ae. speltoides and common wheat (T. aestivum), the nucleotide sequence between a stop codon of rbcL and a HindIII site in cemA in the hotspot was determined for Ae. ovata, Ae. speltoides, Ae. caudata and Ae. mutica. The total number of nucleotides in the region was 2808, 2810, 3302, and 3594 bp, for Ae. speltoides, Ae. ovata, Ae. caudata and Ae. mutica, respectively, and the sequences were compared with the corresponding ones of Ae. crassa 4x, T. aestivum and Ae. squarrosa. Compared with the largest B2l fragment of Ae. mutica, a 791bp and a 793 bp deletion were found in Ae. speltoides and Ae. ovata, respectively, and the possible site of deletion in the two species is the same as that of T. aestivum. However, a deleted segment in Ae. ovata is 2 bp longer than that of Ae. speltoides (and T. aestivum), demonstrating that recurrent deletions had occurred in the chloroplast genomes of both species. Comparison of the sequences from Ae. caudata and Ae. crassa 4x with that of Ae. mutica revealed a 289 bp and a 61 bp deletion at the same site in Ae. caudata and Ae. crassa 4x, respectively. Sequence comparison using wild Aegilops plants showed that the large length variations in a hotspot are fixed to each species. A considerable number of polymorphisms are observed in a loop in the 3' of rbcL. The study reveals the relative importance of the large and small indels and minute inversions to account for variations in the chloroplast genomes among closely related species.

Inferring the pattern of spontaneous mutation from the pattern of substitution in unitary pseudogenes of Mycobacterium leprae and a comparison of mutation patterns among distantly related organisms

The pattern of spontaneous mutation can be inferred from the pattern of substitution in pseudogenes, which are known to be under very weak or no selective constraint. We modified an existing method (Gojobori T, et al., J Mol Evol 18:360, 1982) to infer the pattern of mutation in bacteria by using 569 pseudogenes from Mycobacterium leprae. In Gojobori et al.'s method, the pattern is inferred by using comparisons involving a pseudogene, a conspecific functional paralog, and an outgroup functional ortholog. Because pseudogenes in M. leprae are unitary, we replaced the missing paralogs by functional orthologs from M. tuberculosis. Functional orthologs from Streptomyces coelicolor served as outgroups. We compiled a database consisting of 69,378 inferred mutations. Transitional mutations were found to constitute more than 56% of all mutations. The transitional bias was mainly due to C-->T and G-->A, which were also the most frequent mutations on the leading strand and the only ones that were significantly more frequent than the random expectation. The least frequent mutations on the leading strand were A-->T and T-->A, each with a relative frequency of less than 3%. The mutation pattern was found to differ between the leading and the lagging strands. This asymmetry is thought to be the cause for the typical chirochoric structure of bacterial genomes. The physical distance of the pseudogene from the origin of replication (ori) was found to have almost no effect on the pattern of mutation. A surprising similarity was found between the mutation pattern in M. leprae and previously inferred patterns for such distant taxa as human and Drosophila. The mutation pattern on the leading strand of M. leprae was also found to share some common features with the pattern inferred for the heavy strand of the human mitochondrial genome. These findings indicate that taxon-specific factors may only play secondary roles in determining patterns of mutation.

Structural analysis of chloroplast DNA in Prunus (Rosaceae): evolution, genetic diversity and unequal mutations

In order to understand the evolutionary aspects of the chloroplast DNA (cpDNA) structures in Rosaceous plants, a physical map of peach (Prunus persica cv. Hakuhou) cpDNA was constructed. Fourteen lambda phage clones which covered the entire sequence of the peach cpDNA were digested by restriction enzymes (SalI, XhoI, BamHI, SacI, and PstI) used singly or in combination. The molecular size of peach cpDNA was estimated to be about 152 kb. The gene order and contents were revealed to be equivalent to those of standard type of angiosperms by the localization of 31 genes on the physical map. Eighteen accessions from 14 Prunus species (P. persica, P. mira, P. davidiana, P. cerasis, P. cerasifera, P. domestica, P. insititia, P. spinosa, P. salicina, P. maritima, P. armeniaca, P. mume, P. tomentosa, P. zippeliana, and P. salicifolia) and one interspecific hybrid were used for the structural analysis of cpDNAs. Seventeen mutations (16 recognition site changes and one length mutation) were found in the cpDNA of these 18 accessions by RFLP analysis allowing a classification into 11 genome types. Although the base substitution rate in the recognition site (100p = 0.72) of cpDNA in Prunus was similar to that of other plants, i.e., Triticum-Aegilops, Brassica, and Pisum, it differed from Pyrus (100p = 0.15) in Rosaceae. Seven mutations including one length mutation were densely located within a region of about 9.1 kb which includes psbA and atpA in the left border of a large single-copy region of Prunus cpDNAs. The length mutation was detected only in P. persica and consisted of a 277 bp deletion which occurred in a spacer region between the trnS and trnG genes within the 9.1 kb region. Additional fragment length mutations (insertion/deletion), which were not detected by RFLP analysis, were revealed by PCR and sequence analyses in P. zippeliana and P. salicifolia. All of these length mutations occurred within the 9.1 kb region between psbA and atpA. This region could be an intra-molecular recombinational hotspot in Prunus species.

Compositional properties of green-plant plastid genomes

We studied variation of GC contents among plastid (Pt) genomes of green plants. In the green plants, the GC contents of the whole Pt genomes range from 42.14 to 28.81%. These values are similar to those observed in the mitochondrial (Mt) genomes of the green plants, however, the GC contents in the Pt genomes are not related to those in the Mt genomes or the nuclear (Nc) genomes. In addition, some compositional properties of the three types of genomes are different. Thus, it is suggested that the GC contents of the Pt genomes are maintained independently of the other genomes within a cell. We found that the compositional bias toward AT is strong at the third codon position and in intergenic spacer (IGS) regions in the Pt genomes, and the GC contents (GC3 and GCIGS) at these sites are generally similar within each genome. Additionally, the GC3 and GCIGS are strongly related to the whole-genome GC content. Therefore, the interspecific variation of the GC contents in the Pt genomes is suggested to be mainly caused by the variation of the GC3 and GCIGS, both of which are considered to be under weak selective constraints. Using a maximum likelihood approach, we estimated equilibrium GC3 (eqGC(3)) of 12 genes in the land-plant Pt genomes. We found an increase in eqGC(3) after the divergence of liverworts. These results suggest that genome-wide factors such as GC mutational bias are important for the biased base composition in the Pt genomes.

Mitochondrial genome recombination in the zone of contact between two hybridizing conifers

Variation in mitochondrial DNA was surveyed at four gene loci in and around the zone of contact between two naturally hybridizing conifers, black spruce (Picea mariana) and red spruce (P. rubens) in northeastern North America. Most of the mtDNA diversity of these species was found in populations next to or into the zone of contact, where some individuals bore rare mitotypes intermediate between the common mitotypes observed in the allopatric areas of each species. Sequence analysis and tests for mtDNA recombination point to this phenomenon, rather than to recurrent mutation, as the most tenable hypothesis for the origin of these rare mitotypes. From the 10 mitotypes observed, at least 4 would be the product of recombination between 4 of the 5 putative ancestral mitotypes. Tests for cytonuclear disequilibrium and geographical structure of the putative recombinant mitotypes suggest that mtDNA recombination is not frequent and relatively recent on the geological time scale. mtDNA recombination would have been promoted by transient heteroplasmy due to leakage of paternal mtDNA since the Holocene secondary contact between the two species.

The tortoise and the hare II: relative utility of 21 noncoding chloroplast DNA sequences for phylogenetic analysis

Chloroplast DNA sequences are a primary source of data for plant molecular systematic studies. A few key papers have provided the molecular systematics community with universal primer pairs for noncoding regions that have dominated the field, namely trnL-trnF and trnK/matK. These two regions have provided adequate information to resolve species relationships in some taxa, but often provide little resolution at low taxonomic levels. To obtain better phylogenetic resolution, sequence data from these regions are often coupled with other sequence data. Choosing an appropriate cpDNA region for phylogenetic investigation is difficult because of the scarcity of information about the tempo of evolutionary rates among different noncoding cpDNA regions. The focus of this investigation was to determine whether there is any predictable rate heterogeneity among 21 noncoding cpDNA regions identified as phylogenetically useful at low levels. To test for rate heterogeneity among the different cpDNA regions, we used three species from each of 10 groups representing eight major phylogenetic lineages of phanerogams. The results of this study clearly show that a survey using as few as three representative taxa can be predictive of the amount of phylogenetic information offered by a cpDNA region and that rate heterogeneity exists among noncoding cpDNA regions.

The complete nucleotide sequence of wild rice (Oryza nivara) chloroplast genome: first genome wide comparative sequence analysis of wild and cultivated rice

We determined the complete nucleotide sequence of the chloroplast genome of wild rice, Oryza nivara and compared it with the corresponding published sequence of relative cultivated rice, Oryza sativa. The genome was 134,494 bp long with a large single-copy region of 80,544 bp, a small single-copy region of 12,346 bp and two inverted repeats of 20,802 bp each. The overall A+T content was 61.0%. The O. nivara chloroplast genome encoded identical functional genes to O. sativa in the same order along the genome. On the other hand, detailed analysis revealed 57 insertion, 61 deletion and 159 base substitution events in the entire chloroplast genome of O. nivara. Among substitutions, transversions were much higher than transitions with the former even more frequent than the latter in the coding region. Most of the insertions/deletions were single-base but a few large length mutations were also detected. The frequency of insertion/deletion events was more in the coding region within inverted repeats. In contrast, a very few substitution events were identified in the coding region. Polymorphism was observed among rice cultivars at loci of large insertion/deletion events. This is the first report describing comparative and genome wide chloroplast analysis between a wild and cultivated crop.

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.

Short dispersed repeats in the Chlamydomonas chloroplast genome are collocated with sites for mRNA 3' end formation

The Chlamydomonas reinhardtii chloroplast genome possesses thousands of small dispersed repeats (SDRs), which are of unknown function. Here, we used the petA gene as a model to investigate the role of SDRs in mRNA 3' end formation. In wild-type cells, petA mRNA accumulated as a major 1.3-kb transcript, whose 3' end was mapped to the distal end of a predicted stem-loop structure. To determine whether this stem-loop was required for petA mRNA stability, a series of deletions was constructed. These deletion strains accumulated a variety of petA mRNAs, for which approximate 3' ends were deduced. These 3' ends were found to flank stem-loop structures, many of which were formed partially or completely from inverted copies of SDRs. All strains accumulated wild-type levels of cytochrome f, demonstrating that alternative 3' termini are compatible with efficient translation. The ability to form alternative mRNA termini using SDRs lends additional flexibility to the chloroplast gene expression apparatus and thus could confer an evolutionary advantage.

Noncoding plastid trnT-trnF sequences reveal a well resolved phylogeny of basal angiosperms

Recent contributions from DNA sequences have revolutionized our concept of systematic relationships in angiosperms. However, parts of the angiosperm tree remain unclear. Previous studies have been based on coding or rDNA regions of relatively conserved genes. A phylogeny for basal angiosperms based on noncoding, fast-evolving sequences of the chloroplast genome region trnT-trnF is presented. The recognition of simple direct repeats allowed a robust alignment. Mutational hot spots appear to be confined to certain sectors, as in two stem-loop regions of the trnL intron secondary structure. Our highly resolved and well-supported phylogeny depicts the New Caledonian Amborella as the sister to all other angiosperms, followed by Nymphaeaceae and an Austrobaileya-Illicium-Schisandra clade. Ceratophyllum is substantiated as a close relative of monocots, as is a monophyletic eumagnoliid clade consisting of Piperales plus Winterales sister to Laurales plus Magnoliales. Possible reasons for the striking congruence between the trnT-trnF based phylogeny and phylogenies generated from combined multi-gene, multi-genome data are discussed.

Development and characterization of chloroplast microsatellite markers inMacaranga(Euphorbiaceae)

As part of our study on the phylogeography of the ant-plant genus Macaranga, we have screened for polymorphic regions in the chloroplast genome. Initially, ten universal PCR primer pairs targeted at chloroplast microsatellite loci were applied to a small set of specimens, covering various taxonomic levels from intrafamilial to intraspecific. Eight primer pairs produced PCR fragments that behaved as single and discrete bands on agarose gels. The five most promising candidate pairs were further analysed with an extended set of DNA templates, and PCR products were separated on sequencing gels. The number of size variants per locus varied from two to eight, combining into 17 haplotypes among 29 Macaranga accessions from 10 species. Comparative sequencing demonstrated that microsatellites were responsible for the observed size variation at three of five loci, whereas variation at the other loci was caused by larger insertions and (or) deletions (indels). In addition to poly(A) and poly(T) repeats, which are typically found in chloroplast DNA, we also identified a variable (CT)nrepeat, with n = 4 to n = 8. Sequencing revealed three examples of size homoplasy, one of which was caused by a single base substitution that raised the actual number of haplotypes to 18. Relationships between haplotypes were assessed by phenetic analyses of size variants and by constructing a parsimony network based on sequence variation. For both types of analysis, the distribution of haplotypes correlated with geographically circumscribed regions rather than with taxonomic boundaries.Key words: chloroplast microsatellites, ccmp primers, Macaranga, haplotype network, interspecific gene flow.

Characterization of proteins responsive to gibberellin in the leaf-sheath of rice (Oryza sativa L.) seedling using proteome analysis

In plants gibberellins (GAs) are responsible for triggering stem or internodal elongation. To comprehend the molecular basis of internodal elongation in rice, a proteomics approach using differentially displayed proteins on two-dimensional polyacrylamide gel electrophoresis (2-D PAGE) was carried out to identify the proteins expressed during the GA controlled leaf-sheath elongation response. Out of 352 protein spots detected on 2-D PAGE, 32 proteins showed modulation in the expression levels in GA3-treated leaf-sheath for 48 h as compared to control. These proteins were analyzed using protein sequencer and/or mass spectrometry in conjunction with the protein database to assign putative identities. The twin protein spots (LS079 and LS083), identified as calreticulin, showed different isoelectric points and expression level in GA3-treated leaf-sheath. The expression level of LS083 (pI 4.0) was down-regulated as compared to the up-regulation of LS079 (pI 4.3). In the presence of GA3 and growth inhibitor, uniconazole and abscisic acid, respectively, no elongation in leaf-sheath occurred and calreticulin did not shift from LS083 to LS079. Over-expression of calreticulin in rice inhibited the callus regeneration and seedling growth. These results suggest that calreticulin is an important component in the GA signaling pathway that regulates rice seedling leaf-sheath elongation.

The Chlamydomonas reinhardtii plastid chromosome: islands of genes in a sea of repeats

Chlamydomonas reinhardtii is a unicellular eukaryotic alga possessing a single chloroplast that is widely used as a model system for the study of photosynthetic processes. This report analyzes the surprising structural and evolutionary features of the completely sequenced 203,395-bp plastid chromosome. The genome is divided by 21.2-kb inverted repeats into two single-copy regions of approximately 80 kb and contains only 99 genes, including a full complement of tRNAs and atypical genes encoding the RNA polymerase. A remarkable feature is that >20% of the genome is repetitive DNA: the majority of intergenic regions consist of numerous classes of short dispersed repeats (SDRs), which may have structural or evolutionary significance. Among other sequenced chlorophyte plastid genomes, only that of the green alga Chlorella vulgaris appears to share this feature. The program MultiPipMaker was used to compare the genic complement of Chlamydomonas with those of other chloroplast genomes and to scan the genomes for sequence similarities and repetitive DNAs. Among the results was evidence that the SDRs were not derived from extant coding sequences, although some SDRs may have arisen from other genomic fragments. Phylogenetic reconstruction of changes in plastid genome content revealed that an accelerated rate of gene loss also characterized the Chlamydomonas/Chlorella lineage, a phenomenon that might be independent of the proliferation of SDRs. Together, our results reveal a dynamic and unusual plastid genome whose existence in a model organism will allow its features to be tested functionally.

Molecular analysis of the complete set of length mutations found in the plastomes of Triticum-Aegilops species

Precise location and nature of each of 14 length mutations detected among chloroplast DNAs of Triticum-Aegilops species by RFLP analysis were determined at the nucleotide sequence level. Each mutation was compared with at least three non-mutated wild-type plastomes as standards. These 14 length mutations were classified into 4 duplications and 10 deletions. One duplication occurred in the small single-copy region close to the border of the inverted repeat, and the remaining 13 length mutations took place in the large single-copy region. All length mutations occurred in the intergenic regions, suggesting that these length mutations do not affect plastid gene expression. Saltatory replication was the cause of all duplications, whereas intramolecular recombination mediated by short direct repeats played a substantial role in the deletions. Recurrent occurrences of certain deletion events were found in some AT-rich regions, which constituted hot spots for deletion. Out of four hypervariable regions detected among the grass plastomes, two (downstream of rbcL and a tRNA gene accumulated region) were still active after differentiation of Triticum and Aegilops complex.

Molecular and morphological evidence for an origin of the aberrant genus Milula within himalayan species of Allium (Alliacae)

Phylogenetic relationships between Allium and the monotypic Himalayan genus Milula were analyzed using sequences of the nuclear ribosomal DNA internal transcribed spacer (ITS) region and of the intergenic spacers from the chloroplast trnD(GUC)-trnT(GGU) region. Both marker systems unambiguously placed Milula spicata within Allium subgenus Rhizirideum, close to A. cyathophorum. Morphologically, the main difference between Allium and Milula is the conspicuous spicate inflorescence of Milula vs the mostly capitate or umbellate inflorescences in Allium. Anatomical investigations of leaf characters support a close relationship of Milula with A. cyathophorum and A. mairei, whereas root characters are distinctive from other species of section Cyathophora. To maintain Allium as monophyletic, Milula has been included as A. spicatum in Allium subgenus Rhizirideum.

Hairpins involving both inverted and direct repeats are associated with homoplasious indels in non-coding chloroplast DNA of Taraxacum (Lactuceae: Asteraceae)

Sequence variation in 2.2 kb of non-coding regions of the chloroplast genome of eight dandelions (Taraxacum: Lactuceae) from Asia and Europe is interpreted in the light of the phylogenetic signal of base substitutions vs. indels (insertions-deletions). The four non-coding regions displayed a total of approximately 30 structural mutations of which 9 are potentially phylogenetically informative. Insertions, deletions, and an inversion were found that involved consecutive stretches of up to 172 bases. When compared to phylogenetic relationships of the chloroplast genomes based on nucleotide substitutions only, many homoplasious indels (33%) were detected that differed considerably in length and did not comprise simple sequence repeats typically associated with replication slippage. Though many indels in the intergenic spacers were associated with direct repeats, frequently, the variable stretches participated in inverted repeat stabilized hairpins. In each intergenic spacer or intron examined, nucleotide stretches ranging from 30 to 60 bp were able to fold into stabilized secondary structures. When these indels were homoplasious, they always ranked among the most stabilized hairpins in the non-coding regions. The association of higher order structures that involve both classes of repeats and parallel structural mutations in hot spot regions of the chloroplast genome can be used to differentiate among mutations that differ in phylogenetic reliability.