The complete nucleotide sequence of the tobacco chloroplast genome: its gene organization and expression

Insights from the complete chloroplast genome into the evolution of Sesamum indicum L

Sesame (Sesamum indicum L.) is one of the oldest oilseed crops. In order to investigate the evolutionary characters according to the Sesame Genome Project, apart from sequencing its nuclear genome, we sequenced the complete chloroplast genome of S. indicum cv. Yuzhi 11 (white seeded) using Illumina and 454 sequencing. Comparisons of chloroplast genomes between S. indicum and the 18 other higher plants were then analyzed. The chloroplast genome of cv. Yuzhi 11 contains 153,338 bp and a total of 114 unique genes (KC569603). The number of chloroplast genes in sesame is the same as that in Nicotiana tabacum, Vitis vinifera and Platanus occidentalis. The variation in the length of the large single-copy (LSC) regions and inverted repeats (IR) in sesame compared to 18 other higher plant species was the main contributor to size variation in the cp genome in these species. The 77 functional chloroplast genes, except for ycf1 and ycf2, were highly conserved. The deletion of the cp ycf1 gene sequence in cp genomes may be due either to its transfer to the nuclear genome, as has occurred in sesame, or direct deletion, as has occurred in Panax ginseng and Cucumis sativus. The sesame ycf2 gene is only 5,721 bp in length and has lost about 1,179 bp. Nucleotides 1–585 of ycf2 when queried in BLAST had hits in the sesame draft genome. Five repeats (R10, R12, R13, R14 and R17) were unique to the sesame chloroplast genome. We also found that IR contraction/expansion in the cp genome alters its rate of evolution. Chloroplast genes and repeats display the signature of convergent evolution in sesame and other species. These findings provide a foundation for further investigation of cp genome evolution in Sesamum and other higher plants.

Complete Arabis alpina chloroplast genome sequence and insight into its polymorphism

The alpine plant Arabis alpina (alpine rock-cress) is a thoroughly studied species in the fields of perennial plant flowering regulation, phylogeography, and adaptation to harsh alpine climatic conditions. We report the complete A. alpina chloroplast genome sequence obtained through de novo assembly of Illumina paired-end reads produced by total DNA sequencing. The A. alpina cp circular genome is 152,866 bp in length and built of two inverted repeats of 26,933 bp separated by unique regions: a large single copy of 82,338 bp and a small single copy of 17,938 bp. The genome contains 131 genes, some of them being duplicated in the inverted repeats. Seventy-nine unique protein-coding genes were annotated, together with 29 tRNA genes and 4 ribosomal RNA genes. Sequencing and mapping of 23 additional A. alpina DNA samples enabled to gain insight into the intraspecies polymorphism of the sequenced cp genome. Genetic variability among genomes was detected as 44 indels, most of them being located in noncoding regions, and 130 single-nucleotide polymorphisms, 37 of them corresponding to mutations in coding regions. A. alpina chloroplast genome sequence will be helpful in population studies or investigations of chloroplast functions of this alpine plant species.

•

We report the complete A. alpina chloroplast genome sequence through de novo assembly.

•

Arabis alpina plastome is 152,866 in length and harbors 131 genes.

•

Intraspecies polymorphism analysis was investigated among 24 samples.

•

We detected 44 indels and 130 single nucleotide polymorphisms.

Complete chloroplast genome of Sedum sarmentosum and chloroplast genome evolution in Saxifragales

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

Complete sequencing of five araliaceae chloroplast genomes and the phylogenetic implications

Background

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

Methodology/Principal Findings

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

Conclusion

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

Mechanisms of functional and physical genome reduction in photosynthetic and nonphotosynthetic parasitic plants of the broomrape family

Nonphotosynthetic plants possess strongly reconfigured plastomes attributable to convergent losses of photosynthesis and housekeeping genes, making them excellent systems for studying genome evolution under relaxed selective pressures. We report the complete plastomes of 10 photosynthetic and nonphotosynthetic parasites plus their nonparasitic sister from the broomrape family (Orobanchaceae). By reconstructing the history of gene losses and genome reconfigurations, we find that the establishment of obligate parasitism triggers the relaxation of selective constraints. Partly because of independent losses of one inverted repeat region, Orobanchaceae plastomes vary 3.5-fold in size, with 45 kb in American squawroot (Conopholis americana) representing the smallest plastome reported from land plants. Of the 42 to 74 retained unique genes, only 16 protein genes, 15 tRNAs, and four rRNAs are commonly found. Several holoparasites retain ATP synthase genes with intact open reading frames, suggesting a prolonged function in these plants. The loss of photosynthesis alters the chromosomal architecture in that recombinogenic factors accumulate, fostering large-scale chromosomal rearrangements as functional reduction proceeds. The retention of DNA fragments is strongly influenced by both their proximity to genes under selection and the co-occurrence with those in operons, indicating complex constraints beyond gene function that determine the evolutionary survival time of plastid regions in nonphotosynthetic plants.

Reverse genetics in complex multigene operons by co-transformation of the plastid genome and its application to the open reading frame previously designated psbN

Reverse genetics approaches have contributed enormously to the elucidation of gene functions in plastid genomes and the determination of structure-function relationships in chloroplast multiprotein complexes. Gene knock-outs are usually performed by disrupting the reading frame of interest with a selectable marker cassette. Site-directed mutagenesis is done by placing the marker into the adjacent intergenic spacer and relying on co-integration of the desired mutation by homologous recombination. These strategies are not applicable to genes residing in large multigene operons or other gene-dense genomic regions, because insertion of the marker cassette into an operon-internal gene or into the nearest intergenic spacer is likely to interfere with expression of adjacent genes in the operon or disrupt cis-elements for the expression of neighboring genes and operons. Here we have explored the possibility of using a co-transformation strategy to mutate a small gene of unknown function (psbN) that is embedded in a complex multigene operon. Although inactivation of psbN resulted in strong impairment of photosynthesis, homoplasmic knock-out lines were readily recovered by co-transformation with a selectable marker integrating >38 kb away from the targeted psbN. Our results suggest co-transformation as a suitable strategy for the functional analysis of plastid genes and operons, which allows the recovery of unselected homoplasmic mutants even if the introduced mutations entail a significant selective disadvantage. Moreover, our data provide evidence for involvement of the psbN gene product in the biogenesis of both photosystem I and photosystem II. We therefore propose to rename the gene product 'photosystem biogenesis factor 1' and the gene pbf1.

Variations in efficiency of plastidial RNA editing within ndh transcripts of perennial ryegrass (Lolium perenne) are not linked to differences in drought tolerance

Projected climate change is likely to subject key temperate grassland species, such as perennial ryegrass (Lolium perenne) to drought stress. Previous studies have shown that the NADH dehydrogenase complex (NDH) is involved with countering oxidative stress during environmental stresses like drought. We studied RNA editing within plastidial transcripts of the NDH complex in relation to the drought response of several accessions of perennial ryegrass. We found dramatic and reproducible differences in RNA editing efficiency between accessions, but efficiency was not influenced by imposition of drought stress, and a direct relationship between editing behaviour and drought response was not detected.

Maintenance of healthy grasslands is essential for efficient livestock production, yet projected climate change is likely to place a heavy drought stress burden on key grassland species, such as perennial ryegrass (Lolium perenne). It is therefore important to gather an in-depth knowledge of the underlying plant response to this stress. The present study is focused on RNA editing (post-transcriptional nucleotide modifications resulting in altered transcripts) within plastidial transcripts of the NADH:ubiquinone oxidoreductase (NDH) complex (NADH dehydrogenase complex) in relation to the drought response of several accessions of perennial ryegrass. Previous studies have shown that the NDH complex is involved in countering oxidative stress during environmental stresses like drought. Owing to the nature of RNA editing within this complex, the RNA editing machinery could play a potential role in regulating the activity of the NDH complex. The investigation revealed dramatic and reproducible differences in RNA editing efficiency between accessions, but efficiency was not influenced by imposition of drought stress, and a direct relationship between editing behaviour and drought response was not detected.

The complete chloroplast genome sequence of Mahonia bealei (Berberidaceae) reveals a significant expansion of the inverted repeat and phylogenetic relationship with other angiosperms

Mahonia bealei (Berberidaceae) is a frequently-used traditional Chinese medicinal plant with efficient anti-inflammatory ability. This plant is one of the sources of berberine, a new cholesterol-lowering drug with anti-diabetic activity. We have sequenced the complete nucleotide sequence of the chloroplast (cp) genome of M. bealei. The complete cp genome of M. bealei is 164,792 bp in length, and has a typical structure with large (LSC 73,052 bp) and small (SSC 18,591 bp) single-copy regions separated by a pair of inverted repeats (IRs 36,501 bp) of large size. The Mahonia cp genome contains 111 unique genes and 39 genes are duplicated in the IR regions. The gene order and content of M. bealei are almost unarranged which is consistent with the hypothesis that large IRs stabilize cp genome and reduce gene loss-and-gain probabilities during evolutionary process. A large IR expansion of over 12 kb has occurred in M. bealei, 15 genes (rps19, rpl22, rps3, rpl16, rpl14, rps8, infA, rpl36, rps11, petD, petB, psbH, psbN, psbT and psbB) have expanded to have an additional copy in the IRs. The IR expansion rearrangement occurred via a double-strand DNA break and subsequence repair, which is different from the ordinary gene conversion mechanism. Repeat analysis identified 39 direct/inverted repeats 30 bp or longer with a sequence identity ≥ 90%. Analysis also revealed 75 simple sequence repeat (SSR) loci and almost all are composed of A or T, contributing to a distinct bias in base composition. Comparison of protein-coding sequences with ESTs reveals 9 putative RNA edits and 5 of them resulted in non-synonymous modifications in rpoC1, rps2, rps19 and ycf1. Phylogenetic analysis using maximum parsimony (MP) and maximum likelihood (ML) was performed on a dataset composed of 65 protein-coding genes from 25 taxa, which yields an identical tree topology as previous plastid-based trees, and provides strong support for the sister relationship between Ranunculaceae and Berberidaceae. Molecular dating analyses suggest that Ranunculaceae and Berberidaceae diverged between 90 and 84 mya, which is congruent with the fossil records and with recent estimates of the divergence time of these two taxa.

Evolutionary rewiring: a modified prokaryotic gene-regulatory pathway in chloroplasts

Photosynthetic electron transport regulates chloroplast gene transcription through the action of a bacterial-type sensor kinase known as chloroplast sensor kinase (CSK). CSK represses photosystem I (PS I) gene transcription in PS I light and thus initiates photosystem stoichiometry adjustment. In cyanobacteria and in non-green algae, CSK homologues co-exist with their response regulator partners in canonical bacterial two-component systems. In green algae and plants, however, no response regulator partner of CSK is found. Yeast two-hybrid analysis has revealed interaction of CSK with sigma factor 1 (SIG1) of chloroplast RNA polymerase. Here we present further evidence for the interaction between CSK and SIG1. We also show that CSK interacts with quinone. Arabidopsis SIG1 becomes phosphorylated in PS I light, which then specifically represses transcription of PS I genes. In view of the identical signalling properties of CSK and SIG1 and of their interactions, we suggest that CSK is a SIG1 kinase. We propose that the selective repression of PS I genes arises from the operation of a gene-regulatory phosphoswitch in SIG1. The CSK-SIG1 system represents a novel, rewired chloroplast-signalling pathway created by evolutionary tinkering. This regulatory system supports a proposal for the selection pressure behind the evolutionary stasis of chloroplast genes.

A rapid ribosome profiling method elucidates chloroplast ribosome behavior in vivo

The profiling of ribosome footprints by deep sequencing has revolutionized the analysis of translation by mapping ribosomes with high resolution on a genome-wide scale. We present a variation on this approach that offers a rapid and cost-effective alternative for the genome-wide profiling of chloroplast ribosomes. Ribosome footprints from leaf tissue are hybridized to oligonucleotide tiling microarrays of the plastid ORFeome and report the abundance and translational status of every chloroplast mRNA. Each assay replaces several time-consuming traditional methods while also providing information that was previously inaccessible. To illustrate the utility of the approach, we show that it detects known defects in chloroplast gene expression in several nuclear mutants of maize (Zea mays) and that it reveals previously unsuspected defects. Furthermore, it provided firm answers to several lingering questions in chloroplast gene expression: (1) the overlapping atpB/atpE open reading frames, whose translation had been proposed to be coupled, are translated independently in vivo; (2) splicing is not a prerequisite for translation initiation on an intron-containing chloroplast RNA; and (3) a feedback control mechanism that links the synthesis of ATP synthase subunits in Chlamydomonas reinhardtii does not exist in maize. An analogous approach is likely to be useful for studies of mitochondrial gene expression.

Complete plastid genome sequencing of Trochodendraceae reveals a significant expansion of the inverted repeat and suggests a Paleogene divergence between the two extant species

The early-diverging eudicot order Trochodendrales contains only two monospecific genera, Tetracentron and Trochodendron. Although an extensive fossil record indicates that the clade is perhaps 100 million years old and was widespread throughout the Northern Hemisphere during the Paleogene and Neogene, the two extant genera are both narrowly distributed in eastern Asia. Recent phylogenetic analyses strongly support a clade of Trochodendrales, Buxales, and Gunneridae (core eudicots), but complete plastome analyses do not resolve the relationships among these groups with strong support. However, plastid phylogenomic analyses have not included data for Tetracentron. To better resolve basal eudicot relationships and to clarify when the two extant genera of Trochodendrales diverged, we sequenced the complete plastid genome of Tetracentron sinense using Illumina technology. The Tetracentron and Trochodendron plastomes possess the typical gene content and arrangement that characterize most angiosperm plastid genomes, but both genomes have the same unusual ∼4 kb expansion of the inverted repeat region to include five genes (rpl22, rps3, rpl16, rpl14, and rps8) that are normally found in the large single-copy region. Maximum likelihood analyses of an 83-gene, 88 taxon angiosperm data set yield an identical tree topology as previous plastid-based trees, and moderately support the sister relationship between Buxaceae and Gunneridae. Molecular dating analyses suggest that Tetracentron and Trochodendron diverged between 44-30 million years ago, which is congruent with the fossil record of Trochodendrales and with previous estimates of the divergence time of these two taxa. We also characterize 154 simple sequence repeat loci from the Tetracentron sinense and Trochodendron aralioides plastomes that will be useful in future studies of population genetic structure for these relict species, both of which are of conservation concern.

Additional pathway to translate the downstream ndhK cistron in partially overlapping ndhC-ndhK mRNAs in chloroplasts

The chloroplast NAD(P)H dehydrogenase (NDH) C (ndhC) and ndhK genes partially overlap and are cotranscribed in many plants. We previously reported that the tobacco ndhC/K genes are translationally coupled but produce NdhC and NdhK, subunits of the NDH complex, in similar amounts. Generally, translation of the downstream cistron in overlapping mRNAs is very low. Hence, these findings suggested that the ndhK cistron is translated not only from the ndhC 5'UTR but also by an additional pathway. Using an in vitro translation system from tobacco chloroplasts, we report here that free ribosomes enter, with formylmethionyl-tRNA(fMet), at an internal AUG start codon that is located in frame in the middle of the upstream ndhC cistron, translate the 3' half of the ndhC cistron, reach the ndhK start codon, and that, at that point, some ribosomes resume ndhK translation. We detected a peptide corresponding to a 57-amino-acid product encoded by the sequence from the internal AUG to the ndhC stop codon. We propose a model in which the internal initiation site AUG is not designed for synthesizing a functional isoform but for delivering additional ribosomes to the ndhK cistron to produce NdhK in the amount required for the assembly of the NDH complex. This pathway is a unique type of translation to produce protein in the needed amount with the cost of peptide synthesis.

Essential nucleoid proteins in early chloroplast development

The plastid transcription machinery can be biochemically purified at different organisational levels as soluble RNA polymerase, transcriptionally active chromosome, or nucleoid. Recent proteomic studies have uncovered several novel proteins in these structures and functional genomic studies have indicated that a lack of many of these proteins results in chlorotic phenotypes of varying degree. The most severe cases exhibit complete albino phenotypes, which led to the conclusion that the proteins that were lacking had important regulatory roles in plastid gene expression and chloroplast development. In this opinion article, we propose an alternative model in which the structural establishment of a transcriptional subdomain within the nucleoid represents an early developmental bottleneck that leads to abortion of proper chloroplast biogenesis if disturbed.

Contradiction between plastid gene transcription and function due to complex posttranscriptional splicing: an exemplary study of ycf15 function and evolution in angiosperms

Plant chloroplast genes are usually co-transcribed while its posttranscriptional splicing is fairly complex and remains largely unsolved. On basis of sequencing the three complete Camellia (Theaceae) chloroplast genomes for the first time, we comprehensively analyzed the evolutionary patterns of ycf15, a plastid gene quite paradoxical in terms of its function and evolution, along the inferred angiosperm phylogeny. Although many species in separate lineages including the three species reported here contained an intact ycf15 gene in their chloroplast genomes, the phylogenetic mixture of both intact and obviously disabled ycf15 genes imply that they are all non-functional. Both intracellular gene transfer (IGT) and horizontal gene transfer (HGT) failed to explain such distributional anomalies. While, transcriptome analyses revealed that ycf15 was transcribed as precursor polycistronic transcript which contained ycf2, ycf15 and antisense trnL-CAA. The transcriptome assembly was surprisingly found to cover near the complete Camellia chloroplast genome. Many non-coding regions including pseudogenes were mapped by multiple transcripts, indicating the generality of pseudogene transcriptions. Our results suggest that plastid DNA posttranscriptional splicing may involve complex cleavage of non-functional genes.

Complete chloroplast genome sequences of Mongolia medicine Artemisia frigida and phylogenetic relationships with other plants

Background

Artemisia frigida Willd. is an important Mongolian traditional medicinal plant with pharmacological functions of stanch and detumescence. However, there is little sequence and genomic information available for Artemisia frigida, which makes phylogenetic identification, evolutionary studies, and genetic improvement of its value very difficult. We report the complete chloroplast genome sequence of Artemisia frigida based on 454 pyrosequencing.

Methodology/Principal Findings

The complete chloroplast genome of Artemisia frigida is 151,076 bp including a large single copy (LSC) region of 82,740 bp, a small single copy (SSC) region of 18,394 bp and a pair of inverted repeats (IRs) of 24,971 bp. The genome contains 114 unique genes and 18 duplicated genes. The chloroplast genome of Artemisia frigida contains a small 3.4 kb inversion within a large 23 kb inversion in the LSC region, a unique feature in Asteraceae. The gene order in the SSC region of Artemisia frigida is inverted compared with the other 6 Asteraceae species with the chloroplast genomes sequenced. This inversion is likely caused by an intramolecular recombination event only occurred in Artemisia frigida. The existence of rich SSR loci in the Artemisia frigida chloroplast genome provides a rare opportunity to study population genetics of this Mongolian medicinal plant. Phylogenetic analysis demonstrates a sister relationship between Artemisia frigida and four other species in Asteraceae, including Ageratina adenophora, Helianthus annuus, Guizotia abyssinica and Lactuca sativa, based on 61 protein-coding sequences. Furthermore, Artemisia frigida was placed in the tribe Anthemideae in the subfamily Asteroideae (Asteraceae) based on ndhF and trnL-F sequence comparisons.

Conclusion

The chloroplast genome sequence of Artemisia frigida was assembled and analyzed in this study, representing the first plastid genome sequenced in the Anthemideae tribe. This complete chloroplast genome sequence will be useful for molecular ecology and molecular phylogeny studies within Artemisia species and also within the Asteraceae family.

The complete chloroplast genome sequence of the medicinal plant Salvia miltiorrhiza

Salvia miltiorrhiza is an important medicinal plant with great economic and medicinal value. The complete chloroplast (cp) genome sequence of Salvia miltiorrhiza, the first sequenced member of the Lamiaceae family, is reported here. The genome is 151,328 bp in length and exhibits a typical quadripartite structure of the large (LSC, 82,695 bp) and small (SSC, 17,555 bp) single-copy regions, separated by a pair of inverted repeats (IRs, 25,539 bp). It contains 114 unique genes, including 80 protein-coding genes, 30 tRNAs and four rRNAs. The genome structure, gene order, GC content and codon usage are similar to the typical angiosperm cp genomes. Four forward, three inverted and seven tandem repeats were detected in the Salvia miltiorrhiza cp genome. Simple sequence repeat (SSR) analysis among the 30 asterid cp genomes revealed that most SSRs are AT-rich, which contribute to the overall AT richness of these cp genomes. Additionally, fewer SSRs are distributed in the protein-coding sequences compared to the non-coding regions, indicating an uneven distribution of SSRs within the cp genomes. Entire cp genome comparison of Salvia miltiorrhiza and three other Lamiales cp genomes showed a high degree of sequence similarity and a relatively high divergence of intergenic spacers. Sequence divergence analysis discovered the ten most divergent and ten most conserved genes as well as their length variation, which will be helpful for phylogenetic studies in asterids. Our analysis also supports that both regional and functional constraints affect gene sequence evolution. Further, phylogenetic analysis demonstrated a sister relationship between Salvia miltiorrhiza and Sesamum indicum. The complete cp genome sequence of Salvia miltiorrhiza reported in this paper will facilitate population, phylogenetic and cp genetic engineering studies of this medicinal plant.

Design of chimeric expression elements that confer high-level gene activity in chromoplasts

Non-green plastids, such as chromoplasts, generally have much lower activity of gene expression than chloroplasts in photosynthetically active tissues. Suppression of plastid genes in non-green tissues occurs through a complex interplay of transcriptional and translational control, with the contribution of regulation of transcript abundance versus translational activity being highly variable between genes. Here, we have investigated whether the low expression of the plastid genome in chromoplasts results from inherent limitations in gene expression capacity, or can be overcome by designing appropriate combinations of promoters and translation initiation signals in the 5' untranslated region (5'-UTR). We constructed chimeric expression elements that combine promoters and 5'-UTRs from plastid genes, which are suppressed during chloroplast-to-chromoplast conversion in Solanum lycopersicum (tomato) fruit ripening, either just at the translational level or just at the level of mRNA accumulation. These chimeric expression elements were introduced into the tomato plastid genome by stable chloroplast transformation. We report the identification of promoter-UTR combinations that confer high-level gene expression in chromoplasts of ripe tomato fruits, resulting in the accumulation of reporter protein GFP to up to 1% of total cellular protein. Our work demonstrates that non-green plastids are capable of expressing genes to high levels. Moreover, the chimeric cis-elements for chromoplasts developed here are widely applicable in basic and applied research using transplastomic methods.

Complete chloroplast genome sequence of Magnolia grandiflora and comparative analysis with related species

Magnolia grandiflora is an important medicinal, ornamental and horticultural plant species. The chloroplast (cp) genome of M. grandiflora was sequenced using a 454 sequencing platform and the genome structure was compared with other related species. The complete cp genome of M. grandiflora was 159623 bp in length and contained a pair of inverted repeats (IR) of 26563 bp separated by large and small single copy (LSC, SSC) regions of 87757 and 18740 bp, respectively. A total of 129 genes were successfully annotated, 18 of which included introns. The identity, number and GC content of M. grandiflora cp genes were similar to those of other Magnoliaceae species genomes. Analysis revealed 218 simple sequence repeat (SSR) loci, most composed of A or T, contributing to a bias in base composition. The types and abundances of repeat units in Magnoliaceae species were relatively conserved and these loci will be useful for developing M. grandiflora cp genome vectors. In addition, results indicated that the cp genome size in Magnoliaceae species and the position of the IR border were closely related to the length of the ycf1 gene. Phylogenetic analyses based on 66 shared genes from 30 species using maximum parsimony (MP) and maximum likelihood (ML) methods provided strong support for the phylogenetic position of Magnolia. The availability of the complete cp genome sequence of M. grandiflora provides valuable information for breeding of desirable varieties, cp genetic engineering, developing useful molecular markers and phylogenetic analyses in Magnoliaceae.

Eukaryote Genomes

General overviews of eukaryote genomes are first discussed, including organelle genomes, introns, and junk DNAs. We then discuss the evolutionary features of eukaryote genomes, such as genome duplication, C-value paradox, and the relationship between genome size and mutation rates. Genomes of multicellular organisms, plants, fungi, and animals are then briefly discussed.

The complete chloroplast genome sequence of the medicinal plant Salvia miltiorrhiza

Salvia miltiorrhiza is an important medicinal plant with great economic and medicinal value. The complete chloroplast (cp) genome sequence of Salvia miltiorrhiza, the first sequenced member of the Lamiaceae family, is reported here. The genome is 151,328 bp in length and exhibits a typical quadripartite structure of the large (LSC, 82,695 bp) and small (SSC, 17,555 bp) single-copy regions, separated by a pair of inverted repeats (IRs, 25,539 bp). It contains 114 unique genes, including 80 protein-coding genes, 30 tRNAs and four rRNAs. The genome structure, gene order, GC content and codon usage are similar to the typical angiosperm cp genomes. Four forward, three inverted and seven tandem repeats were detected in the Salvia miltiorrhiza cp genome. Simple sequence repeat (SSR) analysis among the 30 asterid cp genomes revealed that most SSRs are AT-rich, which contribute to the overall AT richness of these cp genomes. Additionally, fewer SSRs are distributed in the protein-coding sequences compared to the non-coding regions, indicating an uneven distribution of SSRs within the cp genomes. Entire cp genome comparison of Salvia miltiorrhiza and three other Lamiales cp genomes showed a high degree of sequence similarity and a relatively high divergence of intergenic spacers. Sequence divergence analysis discovered the ten most divergent and ten most conserved genes as well as their length variation, which will be helpful for phylogenetic studies in asterids. Our analysis also supports that both regional and functional constraints affect gene sequence evolution. Further, phylogenetic analysis demonstrated a sister relationship between Salvia miltiorrhiza and Sesamum indicum. The complete cp genome sequence of Salvia miltiorrhiza reported in this paper will facilitate population, phylogenetic and cp genetic engineering studies of this medicinal plant.

Growth of transplastomic cells expressing D-amino acid oxidase in chloroplasts is tolerant to D-alanine and inhibited by D-valine

Dual-conditional positive/negative selection markers are versatile genetic tools for manipulating genomes. Plastid genomes are relatively small and conserved DNA molecules that can be manipulated precisely by homologous recombination. High-yield expression of recombinant products and maternal inheritance of plastid-encoded traits make plastids attractive sites for modification. Here, we describe the cloning and expression of a dao gene encoding D-amino acid oxidase from Schizosaccharomyces pombe in tobacco (Nicotiana tabacum) plastids. The results provide genetic evidence for the uptake of D-amino acids into plastids, which contain a target that is inhibited by D-alanine. Importantly, this nonantibiotic-based selection system allows the use of cheap and widely available D-amino acids, which are relatively nontoxic to animals and microbes, to either select against (D-valine) or for (D-alanine) cells containing transgenic plastids. Positive/negative selection with d-amino acids was effective in vitro and against transplastomic seedlings grown in soil. The dual functionality of dao is highly suited to the polyploid plastid compartment, where it can be used to provide tolerance against potential D-alanine-based herbicides, control the timing of recombination events such as marker excision, influence the segregation of transgenic plastid genomes, identify loci affecting dao function in mutant screens, and develop D-valine-based methods to manage the spread of transgenic plastids tagged with dao.

The contributions of wobbling and superwobbling to the reading of the genetic code

Reduced bacterial genomes and most genomes of cell organelles (chloroplasts and mitochondria) do not encode the full set of 32 tRNA species required to read all triplets of the genetic code according to the conventional wobble rules. Superwobbling, in which a single tRNA species that contains a uridine in the wobble position of the anticodon reads an entire four-fold degenerate codon box, has been suggested as a possible mechanism for how tRNA sets can be reduced. However, the general feasibility of superwobbling and its efficiency in the various codon boxes have remained unknown. Here we report a complete experimental assessment of the decoding rules in a typical prokaryotic genetic system, the plastid genome. By constructing a large set of transplastomic knock-out mutants for pairs of isoaccepting tRNA species, we show that superwobbling occurs in all codon boxes where it is theoretically possible. Phenotypic characterization of the transplastomic mutant plants revealed that the efficiency of superwobbling varies in a codon box-dependent manner, but--contrary to previous suggestions--it is independent of the number of hydrogen bonds engaged in codon-anticodon interaction. Finally, our data provide experimental evidence of the minimum tRNA set comprising 25 tRNA species, a number lower than previously suggested. Our results demonstrate that all triplets with pyrimidines in third codon position are dually decoded: by a tRNA species utilizing standard base pairing or wobbling and by a second tRNA species employing superwobbling. This has important implications for the interpretation of the genetic code and will aid the construction of synthetic genomes with a minimum-size translational apparatus.

Sequencing of chloroplast genome using whole cellular DNA and solexa sequencing technology

Sequencing of the chloroplast (cp) genome using traditional sequencing methods has been difficult because of its size (>120 kb) and the complicated procedures required to prepare templates. To explore the feasibility of sequencing the cp genome using DNA extracted from whole cells and Solexa sequencing technology, we sequenced whole cellular DNA isolated from leaves of three Brassicarapa accessions with one lane per accession. In total, 246, 362, and 361 Mb sequence data were generated for the three accessions Chiifu-401-42, Z16, and FT, respectively. Micro-reads were assembled by reference-guided assembly using the cpDNA sequences of B. rapa, Arabidopsis thaliana, and Nicotiana tabacum. We achieved coverage of more than 99.96% of the cp genome in the three tested accessions using the B. rapa sequence as the reference. When A. thaliana or N. tabacum sequences were used as references, 99.7-99.8 or 95.5-99.7% of the B. rapa cp genome was covered, respectively. These results demonstrated that sequencing of whole cellular DNA isolated from young leaves using the Illumina Genome Analyzer is an efficient method for high-throughput sequencing of cp genome.

Composition and physiological function of the chloroplast NADH dehydrogenase-like complex in Marchantia polymorpha

The chloroplast NADH dehydrogenase-like (NDH) complex mediates cyclic electron transport and chloro-respiration and consists of five sub-omplexes, which in angiosperms further associate with photosystem I (PSI) to form a super-complex. In Marchantia polymorpha, 11 plastid-encoded subunits and all the nuclear-encoded subunits of the A, B, membrane and ferredoxin-binding sub-complexes are conserved. However, it is unlikely that the genome of this liverwort encodes Lhca5 and Lhca6, both of which mediate NDH-PSI super-complex formation. It is also unlikely that the subunits of the lumen sub-complex, PnsL1-L4, are encoded by the genome. Consistent with this in silico prediction, the results of blue-native gel electrophoresis showed that NDH subunits were detected in a protein complex with lower molecular mass in Marchantia than the NDH-PSI super-complex in Arabidopsis. Using the plastid transformation technique, we knocked out the ndhB gene in Marchantia. Although the wild-type genome copies were completely segregated out, the ΔndhB lines grew like the wild-type photoautotrophically. A post-illumination transient increase in chlorophyll fluorescence, which reflects NDH activity in vivo in angiosperms, was absent in the thalli of the ΔndhB lines. In ruptured chloroplasts, antimycin A-insensitive, and ferredoxin-dependent plastoquinone reduction was impaired, suggesting that chloroplast NDH mediates similar electron transport in Marchantia and Arabidopsis, despite its possible difference in structure. As in angiosperms, linear electron transport was not strongly affected in the ΔndhB lines. However, the plastoquinone pool was slightly more reduced at low light intensity, suggesting that chloroplast NDH functions in redox balancing of the inter system, especially under low light conditions.

Evolutionary constraints on the plastid tRNA set decoding methionine and isoleucine

The plastid (chloroplast) genomes of seed plants typically encode 30 tRNAs. Employing wobble and superwobble mechanisms, most codon boxes are read by only one or two tRNA species. The reduced set of plastid tRNAs follows the evolutionary trend of organellar genomes to shrink in size and coding capacity. A notable exception is the AUN codon box specifying methionine and isoleucine, which is decoded by four tRNA species in nearly all seed plants. However, three of these four tRNA genes were lost from the genomes of some parasitic plastid-containing lineages, possibly suggesting that less than four tRNA species could be sufficient to decode the triplets in the AUN box. To test this hypothesis, we have performed knockout experiments for the four AUN-decoding tRNAs in tobacco (Nicotiana tabacum) plastids. We find that all four tRNA genes are essential under both autotrophic and heterotrophic growth conditions, possibly suggesting tRNA import into plastids of parasitic plastid-bearing species. Phylogenetic analysis of the four plastid tRNA genes reveals striking conservation of all those bacterial features that are involved in discrimination between the different tRNA species containing CAU anticodons.

The comparative chloroplast genomic analysis of photosynthetic orchids and developing DNA markers to distinguish Phalaenopsis orchids

The chloroplast genome of Phalaenopsis equestris was determined and compared to those of Phalaenopsis aphrodite and Oncidium Gower Ramsey in Orchidaceae. The chloroplast genome of P. equestris is 148,959 bp, and a pair of inverted repeats (25,846 bp) separates the genome into large single-copy (85,967 bp) and small single-copy (11,300 bp) regions. The genome encodes 109 genes, including 4 rRNA, 30 tRNA and 75 protein-coding genes, but loses four ndh genes (ndhA, E, F and H) and seven other ndh genes are pseudogenes. The rate of inter-species variation between the two moth orchids was 0.74% (1107 sites) for single nucleotide substitution and 0.24% for insertions (161 sites; 1388 bp) and deletions (189 sites; 1393 bp). The IR regions have a lower rate of nucleotide substitution (3.5-5.8-fold) and indels (4.3-7.1-fold) than single-copy regions. The intergenic spacers are the most divergent, and based on the length variation of the three intergenic spacers, 11 native Phalaenopsis orchids could be successfully distinguished. The coding genes, IR junction and RNA editing sites are relatively more conserved between the two moth orchids than between those of Phalaenopsis and Oncidium spp.

Complete cpDNA genome sequence of Smilax china and phylogenetic placement of Liliales--influences of gene partitions and taxon sampling

The complete nucleotide sequence of the chloroplast genome (cpDNA) of Smilax china L. (Smilacaceae) is reported. It is the first complete cp genome sequence in Liliales. Genomic analyses were conducted to examine the rate and pattern of cpDNA genome evolution in Smilax relative to other major lineages of monocots. The cpDNA genomic sequences were combined with those available for Lilium to evaluate the phylogenetic position of Liliales and to investigate the influence of taxon sampling, gene sampling, gene function, natural selection, and substitution rate on phylogenetic inference in monocots. Phylogenetic analyses using sequence data of gene groups partitioned according to gene function, selection force, and total substitution rate demonstrated evident impacts of these factors on phylogenetic inference of monocots and the placement of Liliales, suggesting potential evolutionary convergence or adaptation of some cpDNA genes in monocots. Our study also demonstrated that reduced taxon sampling reduced the bootstrap support for the placement of Liliales in the cpDNA phylogenomic analysis. Analyses of sequences of 77 protein genes with some missing data and sequences of 81 genes (all protein genes plus the rRNA genes) support a sister relationship of Liliales to the commelinids-Asparagales clade, consistent with the APG III system. Analyses of 63 cpDNA protein genes for 32 taxa with few missing data, however, support a sister relationship of Liliales (represented by Smilax and Lilium) to Dioscoreales-Pandanales. Topology tests indicated that these two alignments do not significantly differ given any of these three cpDNA genomic sequence data sets. Furthermore, we found no saturation effect of the data, suggesting that the cpDNA genomic sequence data used in the study are appropriate for monocot phylogenetic study and long-branch attraction is unlikely to be the cause to explain the result of two well-supported, conflict placements of Liliales. Further analyses using sufficient nuclear data remain necessary to evaluate these two phylogenetic hypotheses regarding the position of Liliales and to address the causes of signal conflict among genes and partitions.

Complete chloroplast genome sequences of important oilseed crop Sesamum indicum L

Sesamum indicum is an important crop plant species for yielding oil. The complete chloroplast (cp) genome of S. indicum (GenBank acc no. JN637766) is 153,324 bp in length, and has a pair of inverted repeat (IR) regions consisting of 25,141 bp each. The lengths of the large single copy (LSC) and the small single copy (SSC) regions are 85,170 bp and 17,872 bp, respectively. Comparative cp DNA sequence analyses of S. indicum with other cp genomes reveal that the genome structure, gene order, gene and intron contents, AT contents, codon usage, and transcription units are similar to the typical angiosperm cp genomes. Nucleotide diversity of the IR region between Sesamum and three other cp genomes is much lower than that of the LSC and SSC regions in both the coding region and noncoding region. As a summary, the regional constraints strongly affect the sequence evolution of the cp genomes, while the functional constraints weakly affect the sequence evolution of cp genomes. Five short inversions associated with short palindromic sequences that form step-loop structures were observed in the chloroplast genome of S. indicum. Twenty-eight different simple sequence repeat loci have been detected in the chloroplast genome of S. indicum. Almost all of the SSR loci were composed of A or T, so this may also contribute to the A-T richness of the cp genome of S. indicum. Seven large repeated loci in the chloroplast genome of S. indicum were also identified and these loci are useful to developing S. indicum-specific cp genome vectors. The complete cp DNA sequences of S. indicum reported in this paper are prerequisite to modifying this important oilseed crop by cp genetic engineering techniques.

Complete chloroplast genome sequence of a major invasive species, crofton weed (Ageratina adenophora)

Background

Crofton weed (Ageratina adenophora) is one of the most hazardous invasive plant species, which causes serious economic losses and environmental damages worldwide. However, the sequence resource and genome information of A. adenophora are rather limited, making phylogenetic identification and evolutionary studies very difficult. Here, we report the complete sequence of the A. adenophora chloroplast (cp) genome based on Illumina sequencing.

Methodology/Principal Findings

The A. adenophora cp genome is 150, 689 bp in length including a small single-copy (SSC) region of 18, 358 bp and a large single-copy (LSC) region of 84, 815 bp separated by a pair of inverted repeats (IRs) of 23, 755 bp. The genome contains 130 unique genes and 18 duplicated in the IR regions, with the gene content and organization similar to other Asteraceae cp genomes. Comparative analysis identified five DNA regions (ndhD-ccsA, psbI-trnS, ndhF-ycf1, ndhI-ndhG and atpA-trnR) containing parsimony-informative characters higher than 2%, which may be potential informative markers for barcoding and phylogenetic analysis. Repeat structure, codon usage and contraction of the IR were also investigated to reveal the pattern of evolution. Phylogenetic analysis demonstrated a sister relationship between A. adenophora and Guizotia abyssinica and supported a monophyly of the Asterales.

Conclusion

We have assembled and analyzed the chloroplast genome of A. adenophora in this study, which was the first sequenced plastome in the Eupatorieae tribe. The complete chloroplast genome information is useful for plant phylogenetic and evolutionary studies within this invasive species and also within the Asteraceae family.

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.

Plastid genome evolution in mycoheterotrophic Ericaceae

Unlike parasitic plants, which are linked to their hosts directly through haustoria, mycoheterotrophic (MHT) plants derive all or part of their water and nutrients from autothrophs via fungal mycorrhizal intermediaries. Ericaceae, the heather family, are a large and diverse group of plants known to form elaborate symbiotic relationships with mycorrhizal fungi. Using PHYA sequence data, we first investigated relationships among mycoheterotrophic Ericaceae and their close autotrophic relatives. Phylogenetic results suggest a minimum of two independent origins of MHT within this family. Additionally, a comparative investigation of plastid genomes (plastomes) grounded within this phylogenetic framework was conducted using a slot-blot Southern hybridization approach. This survey encompassed numerous lineages of Ericaceae with different life histories and trophic levels, including multiple representatives from mixotrophic Pyroleae and fully heterotrophic Monotropeae and Pterosporeae. Fifty-four probes derived from all categories of protein coding genes typically found within the plastomes of flowering plants were used. Our results indicate that the holo-mycoheterotrophic Ericaceae exhibit extensive loss of genes relating to photosynthetic function and expression of the plastome but retain genes with possible functions outside photosynthesis. Mixotrophic taxa tend to retain most genes relating to photosynthetic functions but are varied regarding the plastid ndh gene content. This investigation extends previous inferences that the loss of the NDH complex occurs prior to becoming holo-heterotrophic and it shows that the pattern of gene losses among mycoheterotrophic Ericaceae is similar to that of haustorial parasites. Additionally, we identify the most desirable candidate species for entire plastome sequencing.

Translation of partially overlapping psbD-psbC mRNAs in chloroplasts: the role of 5'-processing and translational coupling

The chloroplast psbD and psbC genes encode the D2 and CP43 proteins of the photosystem II complex, and they are generally cotranscribed. We report studies on the basic translation process of tobacco psbD-psbC mRNAs using an in vitro translation system from tobacco chloroplasts. The primary transcript has an unusually long 5'-UTR (905 nt). We show that it is translatable. Processing of the 5'-UTR greatly enhances the translation efficiency of the psbD cistron. A striking feature is that psbD and psbC cistrons overlap by 14 nt. Removal of the psbD 5'-UTR plus the start codon and introduction of a premature termination codon in the psbD cistron considerably reduce the translation efficiency of the downstream psbC cistron. These results indicate that translation of the psbC cistron depends largely on that of the upstream psbD cistron and thus shows translational coupling; however, a portion is independently translated. These observations, together with the presence of monocistronic psbC mRNAs, suggest that the psbD and psbC cistrons are translated via multiple processes to produce necessary amounts of D2 and CP43 proteins.

DNA fingerprinting, DNA barcoding, and next generation sequencing technology in plants

DNA fingerprinting of plants has become an invaluable tool in forensic, scientific, and industrial laboratories all over the world. PCR has become part of virtually every variation of the plethora of approaches used for DNA fingerprinting today. DNA sequencing is increasingly used either in combination with or as a replacement for traditional DNA fingerprinting techniques. A prime example is the use of short, standardized regions of the genome as taxon barcodes for biological identification of plants. Rapid advances in "next generation sequencing" (NGS) technology are driving down the cost of sequencing and bringing large-scale sequencing projects into the reach of individual investigators. We present an overview of recent publications that demonstrate the use of "NGS" technology for DNA fingerprinting and DNA barcoding applications.

Plastid signalling under multiple conditions is accompanied by a common defect in RNA editing in plastids

Retrograde signalling from the plastid to the nucleus, also known as plastid signalling, plays a key role in coordinating nuclear gene expression with the functional state of plastids. Inhibitors that cause plastid dysfunction have been suggested to generate specific plastid signals related to their modes of action. However, the molecules involved in plastid signalling remain to be identified. Genetic studies indicate that the plastid-localized pentatricopeptide repeat protein GUN1 mediates signalling under several plastid signalling-related conditions. To elucidate further the nature of plastid signals, investigations were carried out to determine whether different plastid signal-inducing treatments had similar effects on plastids and on nuclear gene expression. It is demonstrated that norflurazon and lincomycin treatments and the plastid protein import2-2 (ppi2-2) mutation, which causes a defect in plastid protein import, all resulted in similar changes at the gene expression level. Furthermore, it was observed that these three treatments resulted in defective RNA editing in plastids. This defect in RNA editing was not a secondary effect of down-regulation of pentatricopeptide repeat protein gene expression in the nucleus. The results indicate that these three treatments, which are known to induce plastid signals, affect RNA editing in plastids, suggesting an unprecedented link between plastid signalling and RNA editing.

The primary transcriptome of barley chloroplasts: numerous noncoding RNAs and the dominating role of the plastid-encoded RNA polymerase

Gene expression in plastids of higher plants is dependent on two different transcription machineries, a plastid-encoded bacterial-type RNA polymerase (PEP) and a nuclear-encoded phage-type RNA polymerase (NEP), which recognize distinct types of promoters. The division of labor between PEP and NEP during plastid development and in mature chloroplasts is unclear due to a lack of comprehensive information on promoter usage. Here, we present a thorough investigation into the distribution of PEP and NEP promoters within the plastid genome of barley (Hordeum vulgare). Using a novel differential RNA sequencing approach, which discriminates between primary and processed transcripts, we obtained a genome-wide map of transcription start sites in plastids of mature first leaves. PEP-lacking plastids of the albostrians mutant allowed for the unambiguous identification of NEP promoters. We observed that the chloroplast genome contains many more promoters than genes. According to our data, most genes (including genes coding for photosynthesis proteins) have both PEP and NEP promoters. We also detected numerous transcription start sites within operons, indicating transcriptional uncoupling of genes in polycistronic gene clusters. Moreover, we mapped many transcription start sites in intergenic regions and opposite to annotated genes, demonstrating the existence of numerous noncoding RNA candidates.

Conservation of plastid sequences in the plant nuclear genome for millions of years facilitates endosymbiotic evolution

The nuclear genome of eukaryotes contains large amounts of cytoplasmic organelle DNA (nuclear integrants of organelle DNA [norgs]). The recent sequencing of many mitochondrial and chloroplast genomes has enabled investigation of the potential role of norgs in endosymbiotic evolution. In this article, we describe a new polymerase chain reaction-based method that allows the identification and evolutionary study of recent and older norgs in a range of eukaryotes. We tested this method in the genus Nicotiana and obtained sequences from seven nuclear integrants of plastid DNA (nupts) totaling 25 kb in length. These nupts were estimated to have been transferred 0.033 to 5.81 million years ago. The spectrum of mutations present in the potential protein-coding sequences compared with the noncoding sequences of each nupt revealed that nupts evolve in a nuclear-specific manner and are under neutral evolution. Indels were more frequent in noncoding regions than in potential coding sequences of former chloroplastic DNA, most probably due to the presence of a higher number of homopolymeric sequences. Unexpectedly, some potential protein-coding sequences within the nupts still contained intact open reading frames for up to 5.81 million years. These results suggest that chloroplast genes transferred to the nucleus have in some cases several millions of years to acquire nuclear regulatory elements and become functional. The different factors influencing this time frame and the potential role of nupts in endosymbiotic gene transfer are discussed.

The downstream atpE cistron is efficiently translated via its own cis-element in partially overlapping atpB-atpE dicistronic mRNAs in chloroplasts

The chloroplast atpB and atpE genes encode subunits β and ε of the ATP synthase, respectively. They are co-transcribed as dicistronic mRNAs in flowering plants. An unusual feature is an overlap (AUGA) of the atpB stop codon (UGA) with the atpE start codon (AUG). Hence, atpE translation has been believed to depend on atpB translation (i.e. translational coupling). Using an in vitro translation system from tobacco chloroplasts, we showed that both atpB and atpE cistrons are translated from the tobacco dicistronic mRNA, and that the efficiency of atpB translation is higher than that of atpE translation. When the atpB 5'-UTR was replaced with lower efficiency 5'-UTRs, atpE translation was higher than atpB translation. Removal of the entire atpB 5'-UTR arrested atpB translation but atpE translation still proceeded. Introduction of a premature stop codon in the atpB cistron did not abolish atpE translation. These results indicate that atpE translation is independent of atpB translation. Mutation analysis showed that the atpE cistron possesses its own cis-element(s) for translation, located ~25 nt upstream from the start codon.

Signalling by the global regulatory molecule ppGpp in bacteria and chloroplasts of land plants

The hyperphosphorylated guanine ribonucleotide ppGpp mediates the stringent response in bacteria. Biochemical and genetic studies of this response in Escherichia coli have shown that the biosynthesis of ppGpp is catalysed by two homologous enzymes, RelA and SpoT. RelA is activated in response to amino acid starvation, and SpoT responds to abiotic physical stress beside nutritional stress. All free-living bacteria, including Gram-positive firmicutes, contain RelA-SpoT homologues (RSH). Further, novel ppGpp biosynthetic enzymes, designated small alarmone synthetases (SASs), were recently identified in a subset of bacteria, including the Gram-positive organism Bacillus subtilis, and were shown to consist only of a ppGpp synthetase domain. Studies suggest that these SAS proteins contribute to ppGpp signalling in response to stressful conditions in a manner distinct from that of RelA-SpoT enzymes. SAS proteins currently appear to always occur in addition to RSH enzymes in various combinations but never alone. RSHs have also been identified in chloroplasts, organelles of photosynthetic eukaryotes that originated from endosymbiotic photosynthetic bacteria. These chloroplast RSHs are exclusively encoded in nuclear DNA and targeted into chloroplasts. The findings suggest that ppGpp may regulate chloroplast functions similar to those regulated in bacteria, including transcription and translation. In addition, a novel ppGpp synthetase that is regulated by Ca²⁺ as a result of the presence of two EF-hand motifs at its COOH terminus was recently identified in chloroplasts of land plants. This finding indicates the existence of a direct connection between eukaryotic Ca²⁺ signalling and prokaryotic ppGpp signalling in chloroplasts. The new observations with regard to ppGpp signalling in land plants suggest that such signalling contributes to the regulation of a wider range of cellular functions than previously anticipated.

Nonessential plastid-encoded ribosomal proteins in tobacco: a developmental role for plastid translation and implications for reductive genome evolution

Plastid genomes of higher plants contain a conserved set of ribosomal protein genes. Although plastid translational activity is essential for cell survival in tobacco (Nicotiana tabacum), individual plastid ribosomal proteins can be nonessential. Candidates for nonessential plastid ribosomal proteins are ribosomal proteins identified as nonessential in bacteria and those whose genes were lost from the highly reduced plastid genomes of nonphotosynthetic plastid-bearing lineages (parasitic plants, apicomplexan protozoa). Here we report the reverse genetic analysis of seven plastid-encoded ribosomal proteins that meet these criteria. We have introduced knockout alleles for the corresponding genes into the tobacco plastid genome. Five of the targeted genes (ribosomal protein of the large subunit22 [rpl22], rpl23, rpl32, ribosomal protein of the small subunit3 [rps3], and rps16) were shown to be essential even under heterotrophic conditions, despite their loss in at least some parasitic plastid-bearing lineages. This suggests that nonphotosynthetic plastids show elevated rates of gene transfer to the nuclear genome. Knockout of two ribosomal protein genes, rps15 and rpl36, yielded homoplasmic transplastomic mutants, thus indicating nonessentiality. Whereas Δrps15 plants showed only a mild phenotype, Δrpl36 plants were severely impaired in photosynthesis and growth and, moreover, displayed greatly altered leaf morphology. This finding provides strong genetic evidence that chloroplast translational activity influences leaf development, presumably via a retrograde signaling pathway.

Plastid expression of a double-pentameric vaccine candidate containing human papillomavirus-16 L1 antigen fused with LTB as adjuvant: transplastomic plants show pleiotropic phenotypes

Human papillomavirus (HPV) causes cervical cancer in women worldwide, which is currently prevented by vaccines based on virus-like particles (VLPs). However, these vaccines have certain limitations in their availability to developing countries, largely due to elevated costs. Concerning the highest burden of disease in resource-poor countries, development of an improved mucosal and cost-effective vaccine is a necessity. As an alternative to VLPs, capsomeres have been shown to be highly immunogenic and can be used as vaccine candidate. Furthermore, coupling of an adjuvant like Escherichia coli heat-labile enterotoxin subunit B (LTB) to an antigen can increase its immunogenicity and reduce the costs related to separate co-administration of adjuvants. Our study demonstrates the expression of two pentameric proteins: the modified HPV-16 L1 (L1_2xCysM) and LTB as a fusion protein in tobacco chloroplasts. Homoplasmy of the transplastomic plants was confirmed by Southern blotting. Western blot analysis showed that the LTB-L1 fusion protein was properly expressed in the plastids and the recombinant protein was estimated to accumulate up to 2% of total soluble protein. Proper folding and display of conformational epitopes for both LTB and L1 in the fusion protein was confirmed by GM1-ganglioside binding assay and antigen capture ELISA, respectively. However, all transplastomic lines showed chlorosis, male sterility and growth retardation, which persisted in the ensuing four generations studied. Nevertheless, plants reached maturity and produced seeds by pollination with wild-type plants. Taken together, these results pave the way for the possible development of a low-cost adjuvant-coupled vaccine with potentially improved immunogenicity against cervical cancer.

Complete chloroplast genome sequence of Magnolia kwangsiensis (Magnoliaceae): implication for DNA barcoding and population genetics

Here, we report a completely sequenced plastome using Illumina/Solexa sequencing-by-synthesis (SBS) technology. The plastome of Magnolia kwangsiensis Figlar & Noot. is 159 667 bp in length with a typical quadripartite structure: 88 030 bp large single-copy (LSC) and 18 669 bp small single-copy (SSC) regions, separated by two 26 484 bp inverted repeat (IR) regions. The overall predicted gene number is 129, among which 17 genes are duplicated in IR regions. The plastome of M. kwangsiensis is identical in its gene order to previously published plastomes of magnoliids. Furthermore, the C-to-U type RNA editing frequency of 114 seed plants is positively correlated with plastome GC content and plastome length, whereas plastome length is not correlated with GC content. A total of 16 potential putative barcoding or low taxonomic level phylogenetic study markers in Magnoliaceae were detected by comparing the coding and noncoding regions of the plastome of M. kwangsiensis with that of Liriodendron tulipifera L. At least eight markers might be applied not only to Magnoliaceae but also to other taxa. The 86 mononucleotide cpSSRs that distributed in single-copy noncoding regions are highly valuable to study population genetics and conservation genetics of this endangered rare species.

The origin and characterization of new nuclear genes originating from a cytoplasmic organellar genome

Endosymbiotic transfer of DNA and functional genes from the cytoplasmic organelles (mitochondria and chloroplasts) to the nucleus has been a major factor driving the origin of new nuclear genes, a process central to eukaryote evolution. Although organelle DNA transfers very frequently to the nucleus, most is quickly deleted, decays, or is alternatively scrapped. However, a very small proportion of it gives rise, immediately or eventually, to functional genes. To simulate the process of functional transfer, we screened for nuclear activation of a chloroplast reporter gene aadA, which had been transferred from the chloroplast to independent nuclear loci in 16 different plant lines. Cryptic nuclear activity of the chloroplast promoter was revealed, which became conspicuous when present in multiple nuclear copies. We screened ∼50 million cells of each line and retrieved three plants in which aadA showed strong nuclear activation. Activation occurred by acquisition of the CaMV 35S nuclear promoter or by nuclear activation of the native chloroplast promoter. Two fortuitous sites within the 3' UTR of aadA mRNA both promoted polyadenylation without any sequence change. Complete characterization of one nuclear sequence before and after gene transfer demonstrated integration by nonhomologous end joining involving simultaneous insertion of multiple chloroplast DNA fragments. The real-time observation of three different means by which a chloroplast gene can become expressed in the nucleus suggests that the process, though rare, may be more readily achieved than previously envisaged.

The evolution of the plastid chromosome in land plants: gene content, gene order, gene function

This review bridges functional and evolutionary aspects of plastid chromosome architecture in land plants and their putative ancestors. We provide an overview on the structure and composition of the plastid genome of land plants as well as the functions of its genes in an explicit phylogenetic and evolutionary context. We will discuss the architecture of land plant plastid chromosomes, including gene content and synteny across land plants. Moreover, we will explore the functions and roles of plastid encoded genes in metabolism and their evolutionary importance regarding gene retention and conservation. We suggest that the slow mode at which the plastome typically evolves is likely to be influenced by a combination of different molecular mechanisms. These include the organization of plastid genes in operons, the usually uniparental mode of plastid inheritance, the activity of highly effective repair mechanisms as well as the rarity of plastid fusion. Nevertheless, structurally rearranged plastomes can be found in several unrelated lineages (e.g. ferns, Pinaceae, multiple angiosperm families). Rearrangements and gene losses seem to correlate with an unusual mode of plastid transmission, abundance of repeats, or a heterotrophic lifestyle (parasites or myco-heterotrophs). While only a few functional gene gains and more frequent gene losses have been inferred for land plants, the plastid Ndh complex is one example of multiple independent gene losses and will be discussed in detail. Patterns of ndh-gene loss and functional analyses indicate that these losses are usually found in plant groups with a certain degree of heterotrophy, might rendering plastid encoded Ndh1 subunits dispensable.

High-throughput sequencing of six bamboo chloroplast genomes: phylogenetic implications for temperate woody bamboos (Poaceae: Bambusoideae)

Background

Bambusoideae is the only subfamily that contains woody members in the grass family, Poaceae. In phylogenetic analyses, Bambusoideae, Pooideae and Ehrhartoideae formed the BEP clade, yet the internal relationships of this clade are controversial. The distinctive life history (infrequent flowering and predominance of asexual reproduction) of woody bamboos makes them an interesting but taxonomically difficult group. Phylogenetic analyses based on large DNA fragments could only provide a moderate resolution of woody bamboo relationships, although a robust phylogenetic tree is needed to elucidate their evolutionary history. Phylogenomics is an alternative choice for resolving difficult phylogenies.

Methodology/Principal Findings

Here we present the complete nucleotide sequences of six woody bamboo chloroplast (cp) genomes using Illumina sequencing. These genomes are similar to those of other grasses and rather conservative in evolution. We constructed a phylogeny of Poaceae from 24 complete cp genomes including 21 grass species. Within the BEP clade, we found strong support for a sister relationship between Bambusoideae and Pooideae. In a substantial improvement over prior studies, all six nodes within Bambusoideae were supported with ≥0.95 posterior probability from Bayesian inference and 5/6 nodes resolved with 100% bootstrap support in maximum parsimony and maximum likelihood analyses. We found that repeats in the cp genome could provide phylogenetic information, while caution is needed when using indels in phylogenetic analyses based on few selected genes. We also identified relatively rapidly evolving cp genome regions that have the potential to be used for further phylogenetic study in Bambusoideae.

Conclusions/Significance

The cp genome of Bambusoideae evolved slowly, and phylogenomics based on whole cp genome could be used to resolve major relationships within the subfamily. The difficulty in resolving the diversification among three clades of temperate woody bamboos, even with complete cp genome sequences, suggests that these lineages may have diverged very rapidly.