Complete sequence of the maize chloroplast genome: gene content, hotspots of divergence and fine tuning of genetic information by transcript editing

The First Complete Plastid Genome from Joinvilleaceae (J. ascendens; Poales) Shows Unique and Unpredicted Rearrangements

Joinvilleaceae is a family of tropical grass-like monocots that comprises only the genus Joinvillea. Previous studies have placed Joinvilleaceae in close phylogenetic proximity to the well-studied grass family. A full plastome sequence was determined and characterized for J. ascendens. The plastome was sequenced with next generation methods, fully assembled de novo and annotated. The assembly revealed two novel inversions specific to the Joinvilleaceae lineage and at least one novel plastid inversion in the Joinvilleaceae-Poaceae lineage. Two previously documented inversions in the Joinvilleaceae-Poaceae lineage and one previously documented inversion in the Poaceae lineage were also verified. Inversion events were identified visually and verified computationally by simulation mutations. Additionally, the loss and subsequent degradation of the accD gene in order Poales was explored extensively in Poaceae and J. ascendens. The two novel inversions along with changes in gene composition between families better delimited lineages in the Poales. The presence of large inversions and subsequent reversals in this small family suggested a high potential for large-scale rearrangements to occur in plastid genomes.

Comparative Analysis of the Complete Chloroplast Genomes of Five Quercus Species

Quercus is considered economically and ecologically one of the most important genera in the Northern Hemisphere. Oaks are taxonomically perplexing because of shared interspecific morphological traits and intraspecific morphological variation, which are mainly attributed to hybridization. Universal plastid markers cannot provide a sufficient number of variable sites to explore the phylogeny of this genus, and chloroplast genome-scale data have proven to be useful in resolving intractable phylogenetic relationships. In this study, the complete chloroplast genomes of four Quercus species were sequenced, and one published chloroplast genome of Quercus baronii was retrieved for comparative analyses. The five chloroplast genomes ranged from 161,072 bp (Q. baronii) to 161,237 bp (Q. dolicholepis) in length, and their gene organization and order, and GC content, were similar to those of other Fagaceae species. We analyzed nucleotide substitutions, indels, and repeats in the chloroplast genomes, and found 19 relatively highly variable regions that will potentially provide plastid markers for further taxonomic and phylogenetic studies within Quercus. We observed that four genes (ndhA, ndhK, petA, and ycf1) were subject to positive selection. The phylogenetic relationships of the Quercus species inferred from the chloroplast genomes obtained moderate-to-high support, indicating that chloroplast genome data may be useful in resolving relationships in this genus.

Chloroplast genomes: diversity, evolution, and applications in genetic engineering

Chloroplasts play a crucial role in sustaining life on earth. The availability of over 800 sequenced chloroplast genomes from a variety of land plants has enhanced our understanding of chloroplast biology, intracellular gene transfer, conservation, diversity, and the genetic basis by which chloroplast transgenes can be engineered to enhance plant agronomic traits or to produce high-value agricultural or biomedical products. In this review, we discuss the impact of chloroplast genome sequences on understanding the origins of economically important cultivated species and changes that have taken place during domestication. We also discuss the potential biotechnological applications of chloroplast genomes.

Evolutionary relationships in Panicoid grasses based on plastome phylogenomics (Panicoideae; Poaceae)

BACKGROUND

Panicoideae are the second largest subfamily in Poaceae (grass family), with 212 genera and approximately 3316 species. Previous studies have begun to reveal relationships within the subfamily, but largely lack resolution and/or robust support for certain tribal and subtribal groups. This study aims to resolve these relationships, as well as characterize a putative mitochondrial insert in one linage.

RESULTS

35 newly sequenced Panicoideae plastomes were combined in a phylogenomic study with 37 other species: 15 Panicoideae and 22 from outgroups. A robust Panicoideae topology largely congruent with previous studies was obtained, but with some incongruences with previously reported subtribal relationships. A mitochondrial DNA (mtDNA) to plastid DNA (ptDNA) transfer was discovered in the Paspalum lineage.

CONCLUSIONS

The phylogenomic analysis returned a topology that largely supports previous studies. Five previously recognized subtribes appear on the topology to be non-monophyletic. Additionally, evidence for mtDNA to ptDNA transfer was identified in both Paspalum fimbriatum and P. dilatatum, and suggests a single rare event that took place in a common progenitor. Finally, the framework from this study can guide larger whole plastome sampling to discern the relationships in Cyperochloeae, Steyermarkochloeae, Gynerieae, and other incertae sedis taxa that are weakly supported or unresolved.

The Complete Chloroplast Genome of the Hare's Ear Root, Bupleurum falcatum: Its Molecular Features

Bupleurum falcatum, which belongs to the family Apiaceae, has long been applied for curative treatments, especially as a liver tonic, in herbal medicine. The chloroplast (cp) genome has been an ideal model to perform the evolutionary and comparative studies because of its highly conserved features and simple structure. The Apiaceae family is taxonomically close to the Araliaceae family and there have been numerous complete chloroplast genome sequences reported in the Araliaceae family, while little is known about the Apiaceae family. In this study, the complete sequence of the B. falcatum chloroplast genome was obtained. The full-length of the cp genome is 155,989 nucleotides with a 37.66% overall guanine-cytosine (GC) content and shows a quadripartite structure composed of three nomenclatural regions: a large single-copy (LSC) region, a small single-copy (SSC) region, and a pair of inverted repeat (IR) regions. The genome occupancy is 85,912-bp, 17,517-bp, and 26,280-bp for LSC, SSC, and IR, respectively. B. falcatum was shown to contain 111 unique genes (78 for protein-coding, 29 for tRNAs, and four for rRNAs, respectively) on its chloroplast genome. Genic comparison found that B. falcatum has no pseudogenes and has two gene losses, accD in the LSC and ycf15 in the IRs. A total of 55 unique tandem repeat sequences were detected in the B. falcatum cp genome. This report is the first to describe the complete chloroplast genome sequence in B. falcatum and will open up further avenues of research to understand the evolutionary panorama and the chloroplast genome conformation in related plant species.

PBR1 selectively controls biogenesis of photosynthetic complexes by modulating translation of the large chloroplast gene Ycf1 in Arabidopsis

The biogenesis of photosystem I (PSI), cytochrome b₆f (Cytb₆f) and NADH dehydrogenase (NDH) complexes relies on the spatially and temporally coordinated expression and translation of both nuclear and chloroplast genes. Here we report the identification of photosystem biogenesis regulator 1 (PBR1), a nuclear-encoded chloroplast RNA-binding protein that regulates the concerted biogenesis of NDH, PSI and Cytb₆f complexes. We identified Ycf1, one of the two largest chloroplast genome-encoded open reading frames as the direct downstream target protein of PBR1. Biochemical and molecular analyses reveal that PBR1 regulates Ycf1 translation by directly binding to its mRNA. Surprisingly, we further demonstrate that relocation of the chloroplast gene Ycf1 fused with a plastid-transit sequence to the nucleus bypasses the requirement of PBR1 for Ycf1 translation, which sufficiently complements the defects in biogenesis of NDH, PSI and Cytb₆f complexes in PBR1-deficient plants. Remarkably, the nuclear-encoded PBR1 tightly controls the expression of the chloroplast gene Ycf1 at the translational level, which is sufficient to sustain the coordinated biogenesis of NDH, PSI and Cytb₆f complexes as a whole. Our findings provide deep insights into better understanding of how a predominant nuclear-encoded factor can act as a migratory mediator and undergoes selective translational regulation of the target plastid gene in controlling biogenesis of photosynthetic complexes.

A PPR protein in the PLS subfamily stabilizes the 5'-end of processed rpl16 mRNAs in maize chloroplasts

Pentatricopeptide repeat (PPR) proteins are a large family of helical-repeat proteins that bind RNA in mitochondria and chloroplasts. Precise RNA targets and functions have been assigned to only a small fraction of the >400 members of the PPR family in plants. We used the amino acid code governing the specificity of RNA binding by PPR repeats to infer candidate-binding sites for the maize protein PPR103 and its ortholog Arabidopsis EMB175. Genetic and biochemical data confirmed a predicted binding site in the chloroplast rpl16 5′UTR to be a site of PPR103 action. This site maps to the 5′ end of transcripts that fail to accumulate in ppr103 mutants. A small RNA corresponding to the predicted PPR103 binding site accumulates in a PPR103-dependent fashion, as expected of PPR103's in vivo footprint. Recombinant PPR103 bound specifically to this sequence in vitro. These observations imply that PPR103 stabilizes rpl16 mRNA by impeding 5′→3′ RNA degradation. Previously described PPR proteins with this type of function consist of canonical PPR motifs. By contrast, PPR103 is a PLS-type protein, an architecture typically associated with proteins that specify sites of RNA editing. However, PPR103 is not required to specify editing sites in chloroplasts.

Evolutionary dynamics of the plastid inverted repeat: the effects of expansion, contraction, and loss on substitution rates

Rates of nucleotide substitution were previously shown to be several times slower in the plastid inverted repeat (IR) compared with single-copy (SC) regions, suggesting that the IR provides enhanced copy-correction activity. To examine the generality of this synonymous rate dependence on the IR, we compared plastomes from 69 pairs of closely related species representing 52 families of angiosperms, gymnosperms, and ferns. We explored the breadth of IR boundary shifts in land plants and demonstrate that synonymous substitution rates are, on average, 3.7 times slower in IR genes than in SC genes. In addition, genes moved from the SC into the IR exhibit lower synonymous rates consistent with other IR genes, while genes moved from the IR into the SC exhibit higher rates consistent with other SC genes. Surprisingly, however, several plastid genes from Pelargonium, Plantago, and Silene have highly accelerated synonymous rates despite their IR localization. Together, these results provide strong evidence that the duplicative nature of the IR reduces the substitution rate within this region. The anomalously fast-evolving genes in Pelargonium, Plantago, and Silene indicate localized hypermutation, potentially induced by a higher level of error-prone double-strand break repair in these regions, which generates substitutional rate variation.

Non-functional plastid ndh gene fragments are present in the nuclear genome of Norway spruce (Picea abies L. Karsch): insights from in silico analysis of nuclear and organellar genomes

Many genes have been lost from the prokaryote plastidial genome during the early events of endosymbiosis in eukaryotes. Some of them were definitively lost, but others were relocated and functionally integrated to the host nuclear genomes through serial events of gene transfer during plant evolution. In gymnosperms, plastid genome sequencing has revealed the loss of ndh genes from several species of Gnetales and Pinaceae, including Norway spruce (Picea abies). This study aims to trace the ndh genes in the nuclear and organellar Norway spruce genomes. The plastid genomes of higher plants contain 11 ndh genes which are homologues of mitochondrial genes encoding subunits of the proton-pumping NADH-dehydrogenase (nicotinamide adenine dinucleotide dehydrogenase) or complex I (electron transport chain). Ndh genes encode 11 NDH polypeptides forming the Ndh complex (analogous to complex I) which seems to be primarily involved in chloro-respiration processes. We considered ndh genes from the plastidial genome of four gymnosperms (Cryptomeria japonica, Cycas revoluta, Ginkgo biloba, Podocarpus totara) and a single angiosperm species (Arabidopsis thaliana) to trace putative homologs in the nuclear and organellar Norway spruce genomes using tBLASTn to assess the evolutionary fate of ndh genes in Norway spruce and to address their genomic location(s), structure, integrity and functionality. The results obtained from tBLASTn were subsequently analyzed by performing homology search for finding ndh specific conserved domains using conserved domain search. We report the presence of non-functional plastid ndh gene fragments, excepting ndhE and ndhG genes, in the nuclear genome of Norway spruce. Regulatory transcriptional elements like promoters, TATA boxes and enhancers were detected in the upstream regions of some ndh fragments. We also found transposable elements in the flanking regions of few ndh fragments suggesting nuclear rearrangements in those regions. These evidences support the hypothesis that, at least in Picea, ndh translocations from the plastid to the nuclear genome have occurred, and that there might have been a functional machinery at some time during evolution to accommodate them within a nuclear-encoded environment, or attempts to form it.

Chloroplast Genome Sequence of Pigeonpea (Cajanus cajan (L.) Millspaugh) and Cajanus scarabaeoides (L.) Thouars: Genome Organization and Comparison with Other Legumes

Pigeonpea (Cajanus cajan (L.) Millspaugh), a diploid (2n = 22) legume crop with a genome size of 852 Mbp, serves as an important source of human dietary protein especially in South East Asian and African regions. In this study, the draft chloroplast genomes of Cajanus cajan and Cajanus scarabaeoides (L.) Thouars were generated. Cajanus scarabaeoides is an important species of the Cajanus gene pool and has also been used for developing promising CMS system by different groups. A male sterile genotype harboring the C. scarabaeoides cytoplasm was used for sequencing the plastid genome. The cp genome of C. cajan is 152,242bp long, having a quadripartite structure with LSC of 83,455 bp and SSC of 17,871 bp separated by IRs of 25,398 bp. Similarly, the cp genome of C. scarabaeoides is 152,201bp long, having a quadripartite structure in which IRs of 25,402 bp length separates 83,423 bp of LSC and 17,854 bp of SSC. The pigeonpea cp genome contains 116 unique genes, including 30 tRNA, 4 rRNA, 78 predicted protein coding genes and 5 pseudogenes. A 50 kb inversion was observed in the LSC region of pigeonpea cp genome, consistent with other legumes. Comparison of cp genome with other legumes revealed the contraction of IR boundaries due to the absence of rps19 gene in the IR region. Chloroplast SSRs were mined and a total of 280 and 292 cpSSRs were identified in C. scarabaeoides and C. cajan respectively. RNA editing was observed at 37 sites in both C. scarabaeoides and C. cajan, with maximum occurrence in the ndh genes. The pigeonpea cp genome sequence would be beneficial in providing informative molecular markers which can be utilized for genetic diversity analysis and aid in understanding the plant systematics studies among major grain legumes.

The linear plastid chromosomes of maize: terminal sequences, structures, and implications for DNA replication

The structure of a chromosomal DNA molecule may influence the way in which it is replicated and inherited. For decades plastid DNA (ptDNA) was believed to be circular, with breakage invoked to explain linear forms found upon extraction from the cell. Recent evidence indicates that ptDNA in vivo consists of linear molecules with discrete termini, although these ends were not characterized. We report the sequences of two terminal regions, End1 and End2, for maize (Zea mays L.) ptDNA. We describe structural features of these terminal regions and similarities found in other plant ptDNAs. The terminal sequences are within inverted repeat regions (leading to four genomic isomers) and adjacent to origins of replication. Conceptually, stem-loop structures may be formed following melting of the double-stranded DNA ends. Exonuclease digestion indicates that the ends in maize are unobstructed, but tobacco (Nicotiana tabacum L.) ends may have a 5'-protein. If the terminal structure of ptDNA molecules influences the retention of ptDNA, the unprotected molecular ends in mature leaves of maize may be more susceptible to degradation in vivo than the protected ends in tobacco. The terminal sequences and cumulative GC skew profiles are nearly identical for maize, wheat (Triticum aestivum L.) and rice (Oryza sativa L.), with less similarity among other plants. The linear structure is now confirmed for maize ptDNA and inferred for other plants and suggests a virus-like recombination-dependent replication mechanism for ptDNA. Plastid transformation vectors containing the terminal sequences may increase the chances of success in generating transplastomic cereals.

Embryo defective 14 encodes a plastid-targeted cGTPase essential for embryogenesis in maize

The embryo defective (emb) mutants in maize genetically define a unique class of loci that is required for embryogenesis but not endosperm development, allowing dissection of two developmental processes of seed formation. Through characterization of the emb14 mutant, we report here that Emb14 gene encodes a circular permuted, YqeH class GTPase protein that likely functions in 30S ribosome formation in plastids. Loss of Emb14 function in the null mutant arrests embryogenesis at the early transition stage. Emb14 was cloned by transposon tagging and was confirmed by analysis of four alleles. Subcellular localization indicated that the EMB14 is targeted to chloroplasts. Recombinant EMB14 is shown to hydrolyze GTP in vitro (Km  = 2.42 ± 0.3 μm). Emb14 was constitutively expressed in all tissues examined and high level of expression was found in transition stage embryos. Comparison of emb14 and WT indicated that loss of EMB14 function severely impairs accumulation of 16S rRNA and several plastid encoded ribosomal genes. We show that an EMB14 transgene complements the pale green, slow growth phenotype conditioned by mutations in AtNOA1, a closely related YqeH GTPase of Arabidopsis. Taken together, we propose that the EMB14/AtNOA1/YqeH class GTPases function in assembly of the 30S subunit of the chloroplast ribosome, and that this function is essential to embryogenesis in plants.

RNA Editing in Chloroplasts of Spirodela polyrhiza, an Aquatic Monocotelydonous Species

RNA editing is the post-transcriptional conversion from C to U before translation, providing a unique feature in the regulation of gene expression. Here, we used a robust and efficient method based on RNA-seq from non-ribosomal total RNA to simultaneously measure chloroplast-gene expression and RNA editing efficiency in the Greater Duckweed, Spirodela polyrhiza, a species that provides a new reference for the phylogenetic studies of monocotyledonous plants. We identified 66 editing sites at the genome-wide level, with an average editing efficiency of 76%. We found that the expression levels of chloroplast genes were relatively constant, but 11 RNA editing sites show significant changes in editing efficiency, when fronds turn into turions. Thus, RNA editing efficiency contributes more to the yield of translatable transcripts than steady state mRNA levels. Comparison of RNA editing sites in coconut, Spirodela, maize, and rice suggests that RNA editing originated from a common ancestor.

Analysis of the Complete Chloroplast Genome of a Medicinal Plant, Dianthus superbus var. longicalyncinus, from a Comparative Genomics Perspective

Dianthus superbus var. longicalycinus is an economically important traditional Chinese medicinal plant that is also used for ornamental purposes. In this study, D. superbus was compared to its closely related family of Caryophyllaceae chloroplast (cp) genomes such as Lychnis chalcedonica and Spinacia oleracea. D. superbus had the longest large single copy (LSC) region (82,805 bp), with some variations in the inverted repeat region A (IRA)/LSC regions. The IRs underwent both expansion and constriction during evolution of the Caryophyllaceae family; however, intense variations were not identified. The pseudogene ribosomal protein subunit S19 (rps19) was identified at the IRA/LSC junction, but was not present in the cp genome of other Caryophyllaceae family members. The translation initiation factor IF-1 (infA) and ribosomal protein subunit L23 (rpl23) genes were absent from the Dianthus cp genome. When the cp genome of Dianthus was compared with 31 other angiosperm lineages, the infA gene was found to have been lost in most members of rosids, solanales of asterids and Lychnis of Caryophyllales, whereas rpl23 gene loss or pseudogization had occurred exclusively in Caryophyllales. Nevertheless, the cp genome of Dianthus and Spinacia has two introns in the proteolytic subunit of ATP-dependent protease (clpP) gene, but Lychnis has lost introns from the clpP gene. Furthermore, phylogenetic analysis of individual protein-coding genes infA and rpl23 revealed that gene loss or pseudogenization occurred independently in the cp genome of Dianthus. Molecular phylogenetic analysis also demonstrated a sister relationship between Dianthus and Lychnis based on 78 protein-coding sequences. The results presented herein will contribute to studies of the evolution, molecular biology and genetic engineering of the medicinal and ornamental plant, D. superbus var. longicalycinus.

Evolutionary Stasis in Cycad Plastomes and the First Case of Plastome GC-Biased Gene Conversion

In angiosperms, gene conversion has been known to reduce the mutational load of plastid genomes (the plastomes). Particularly, more frequent gene conversions in inverted repeat (IR) than in single copy (SC) regions result in contrasting substitution rates between these two regions. However, little has been known about the effect of gene conversion in the evolution of gymnosperm plastomes. Cycads (Cycadophyta) are the second largest gymnosperm group. Evolutionary study of their plastomes is limited to the basal cycad genus, Cycas. In this study, we addressed three questions. 1) Do the plastomes of other cycad genera evolve slowly as previously observed in the plastome of Cycas taitungensis? 2) Do substitution rates differ between their SC and IR regions? And 3) Does gene conversion occur in the cycad plastomes? If yes, is it AT-biased or GC-biased? Plastomes of eight species from other eight genera of cycads were sequenced. These plastomes are highly conserved in genome organization. Excluding ginkgo, cycad plastomes have significantly lower synonymous and nonsynonymous substitution rates than other gymnosperms, reflecting their evolutionary stasis in nucleotide mutations. In the IRs of cycad plastomes, the reduced substitution rates and GC-biased mutations are associated with a GC-biased gene conversion (gBGC) mechanism. Further investigations suggest that in cycads, gBGC is able to rectify plastome-wide mutations. Therefore, this study is the first to uncover the plastomic gBGC in seed plants. We also propose a gBGC model to interpret the dissimilar evolutionary patterns as well as the compositionally biased mutations in the SC and IR regions of cycad plastomes.

Transcriptional Slippage and RNA Editing Increase the Diversity of Transcripts in Chloroplasts: Insight from Deep Sequencing of Vigna radiata Genome and Transcriptome

We performed deep sequencing of the nuclear and organellar genomes of three mungbean genotypes: Vigna radiata ssp. sublobata TC1966, V. radiata var. radiata NM92 and the recombinant inbred line RIL59 derived from a cross between TC1966 and NM92. Moreover, we performed deep sequencing of the RIL59 transcriptome to investigate transcript variability. The mungbean chloroplast genome has a quadripartite structure including a pair of inverted repeats separated by two single copy regions. A total of 213 simple sequence repeats were identified in the chloroplast genomes of NM92 and RIL59; 78 single nucleotide variants and nine indels were discovered in comparing the chloroplast genomes of TC1966 and NM92. Analysis of the mungbean chloroplast transcriptome revealed mRNAs that were affected by transcriptional slippage and RNA editing. Transcriptional slippage frequency was positively correlated with the length of simple sequence repeats of the mungbean chloroplast genome (R²=0.9911). In total, 41 C-to-U editing sites were found in 23 chloroplast genes and in one intergenic spacer. No editing site that swapped U to C was found. A combination of bioinformatics and experimental methods revealed that the plastid-encoded RNA polymerase-transcribed genes psbF and ndhA are affected by transcriptional slippage in mungbean and in main lineages of land plants, including three dicots (Glycine max, Brassica rapa, and Nicotiana tabacum), two monocots (Oryza sativa and Zea mays), two gymnosperms (Pinus taeda and Ginkgo biloba) and one moss (Physcomitrella patens). Transcript analysis of the rps2 gene showed that transcriptional slippage could affect transcripts at single sequence repeat regions with poly-A runs. It showed that transcriptional slippage together with incomplete RNA editing may cause sequence diversity of transcripts in chloroplasts of land plants.

Plastid primers for angiosperm phylogenetics and phylogeography

PREMISE OF THE STUDY

PCR primers are available for virtually every region of the plastid genome. Selection of which primer pairs to use is second only to selection of the genic region. This is particularly true for research at the species/population interface.

METHODS

Primer pairs for 130 regions of the chloroplast genome were evaluated in 12 species distributed across the angiosperms. Likelihood of amplification success was inferred based upon number and location of mismatches to target sequence. Intraspecific sequence variability was evaluated under three different criteria in four species.

RESULTS

Many published primer pairs should work across all taxa sampled, with the exception of failure due to genomic reorganization events. Universal barcoding primers were the least likely to work (65% success). The list of most variable regions for use within species has little in common with the lists identified in prior studies.

DISCUSSION

Published primer sequences should amplify a diversity of flowering plant DNAs, even those designed for specific taxonomic groups. "Universal" primers may have extremely limited utility. There was little consistency in likelihood of amplification success for any given publication across lineages or within lineage across publications.

Plastid phylogenomics of the cool-season grass subfamily: clarification of relationships among early-diverging tribes

Whole plastid genomes are being sequenced rapidly from across the green plant tree of life, and phylogenetic analyses of these are increasing resolution and support for relationships that have varied among or been unresolved in earlier single- and multi-gene studies. Pooideae, the cool-season grass lineage, is the largest of the 12 grass subfamilies and includes important temperate cereals, turf grasses and forage species. Although numerous studies of the phylogeny of the subfamily have been undertaken, relationships among some 'early-diverging' tribes conflict among studies, and some relationships among subtribes of Poeae have not yet been resolved. To address these issues, we newly sequenced 25 whole plastomes, which showed rearrangements typical of Poaceae. These plastomes represent 9 tribes and 11 subtribes of Pooideae, and were analysed with 20 existing plastomes for the subfamily. Maximum likelihood (ML), maximum parsimony (MP) and Bayesian inference (BI) robustly resolve most deep relationships in the subfamily. Complete plastome data provide increased nodal support compared with protein-coding data alone at nodes that are not maximally supported. Following the divergence of Brachyelytrum, Phaenospermateae, Brylkinieae-Meliceae and Ampelodesmeae-Stipeae are the successive sister groups of the rest of the subfamily. Ampelodesmeae are nested within Stipeae in the plastome trees, consistent with its hybrid origin between a phaenospermatoid and a stipoid grass (the maternal parent). The core Pooideae are strongly supported and include Brachypodieae, a Bromeae-Triticeae clade and Poeae. Within Poeae, a novel sister group relationship between Phalaridinae and Torreyochloinae is found, and the relative branching order of this clade and Aveninae, with respect to an Agrostidinae-Brizinae clade, are discordant between MP and ML/BI trees. Maximum likelihood and Bayesian analyses strongly support Airinae and Holcinae as the successive sister groups of a Dactylidinae-Loliinae clade.

Genome-wide analysis of thylakoid-bound ribosomes in maize reveals principles of cotranslational targeting to the thylakoid membrane

Chloroplast genomes encode ∼ 37 proteins that integrate into the thylakoid membrane. The mechanisms that target these proteins to the membrane are largely unexplored. We used ribosome profiling to provide a comprehensive, high-resolution map of ribosome positions on chloroplast mRNAs in separated membrane and soluble fractions in maize seedlings. The results show that translation invariably initiates off the thylakoid membrane and that ribosomes synthesizing a subset of membrane proteins subsequently become attached to the membrane in a nuclease-resistant fashion. The transition from soluble to membrane-attached ribosomes occurs shortly after the first transmembrane segment in the nascent peptide has emerged from the ribosome. Membrane proteins whose translation terminates before emergence of a transmembrane segment are translated in the stroma and targeted to the membrane posttranslationally. These results indicate that the first transmembrane segment generally comprises the signal that links ribosomes to thylakoid membranes for cotranslational integration. The sole exception is cytochrome f, whose cleavable N-terminal cpSecA-dependent signal sequence engages the thylakoid membrane cotranslationally. The distinct behavior of ribosomes synthesizing the inner envelope protein CemA indicates that sorting signals for the thylakoid and envelope membranes are distinguished cotranslationally. In addition, the fractionation behavior of ribosomes in polycistronic transcription units encoding both membrane and soluble proteins adds to the evidence that the removal of upstream ORFs by RNA processing is not typically required for the translation of internal genes in polycistronic chloroplast mRNAs.

Chloroplast DNA sequence utility for the lowest phylogenetic and phylogeographic inferences in angiosperms: the tortoise and the hare IV

PREMISE OF THE STUDY

Noncoding chloroplast DNA (NC-cpDNA) sequences are the staple data source of low-level phylogeographic and phylogenetic studies of angiosperms. We followed up on previous papers (tortoise and hare II and III) that sought to identify the most consistently variable regions of NC-cpDNA. We used an exhaustive literature review and newly available whole plastome data to assess applicability of previous conclusions at low taxonomic levels.

METHODS

We aligned complete plastomes of 25 species pairs from across angiosperms, comparing the number of genetic differences found in 107 NC-cpDNA regions and matK. We surveyed Web of Science for the plant phylogeographic literature between 2007 and 2013 to assess how NC-cpDNA has been used at the intraspecific level.

KEY RESULTS

Several regions are consistently the most variable across angiosperm lineages: ndhF-rpl32, rpl32-trnL((UAG)), ndhC-trnV((UAC)), 5'rps16-trnQ((UUG)), psbE-petL, trnT((GGU))-psbD, petA-psbJ, and rpl16 intron. However, there is no universally best region. The average number of regions applied to low-level studies is ∼2.5, which may be too little to access the full discriminating power of this genome.

CONCLUSIONS

Plastome sequences have been used successfully at lower and lower taxonomic levels. Our findings corroborate earlier works, suggesting that there are regions that are most likely to be the most variable. However, while NC-cpDNA sequences are commonly used in plant phylogeographic studies, few of the most variable regions are applied in that context. Furthermore, it appears that in most studies too few NC-cpDNAs are used to access the discriminating power of the cpDNA genome.

Comparative chloroplast genomes of photosynthetic orchids: insights into evolution of the Orchidaceae and development of molecular markers for phylogenetic applications

The orchid family Orchidaceae is one of the largest angiosperm families, including many species of important economic value. While chloroplast genomes are very informative for systematics and species identification, there is very limited information available on chloroplast genomes in the Orchidaceae. Here, we report the complete chloroplast genomes of the medicinal plant Dendrobium officinale and the ornamental orchid Cypripedium macranthos, demonstrating their gene content and order and potential RNA editing sites. The chloroplast genomes of the above two species and five known photosynthetic orchids showed similarities in structure as well as gene order and content, but differences in the organization of the inverted repeat/small single-copy junction and ndh genes. The organization of the inverted repeat/small single-copy junctions in the chloroplast genomes of these orchids was classified into four types; we propose that inverted repeats flanking the small single-copy region underwent expansion or contraction among Orchidaceae. The AT-rich regions of the ycf1 gene in orchids could be linked to the recombination of inverted repeat/small single-copy junctions. Relative species in orchids displayed similar patterns of variation in ndh gene contents. Furthermore, fifteen highly divergent protein-coding genes were identified, which are useful for phylogenetic analyses in orchids. To test the efficiency of these genes serving as markers in phylogenetic analyses, coding regions of four genes (accD, ccsA, matK, and ycf1) were used as a case study to construct phylogenetic trees in the subfamily Epidendroideae. High support was obtained for placement of previously unlocated subtribes Collabiinae and Dendrobiinae in the subfamily Epidendroideae. Our findings expand understanding of the diversity of orchid chloroplast genomes and provide a reference for study of the molecular systematics of this family.

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.

A major role for the plastid-encoded RNA polymerase complex in the expression of plastid transfer RNAs

Chloroplast transcription in land plants relies on collaboration between a plastid-encoded RNA polymerase (PEP) of cyanobacterial ancestry and a nucleus-encoded RNA polymerase of phage ancestry. PEP associates with additional proteins that are unrelated to bacterial transcription factors, many of which have been shown to be important for PEP activity in Arabidopsis (Arabidopsis thaliana). However, the biochemical roles of these PEP-associated proteins are not known. We describe phenotypes conditioned by transposon insertions in genes encoding the maize (Zea mays) orthologs of five such proteins: ZmPTAC2, ZmMurE, ZmPTAC10, ZmPTAC12, and ZmPRIN2. These mutants have similar ivory/virescent pigmentation and similar reductions in plastid ribosomes and photosynthetic complexes. RNA gel-blot and microarray hybridizations revealed numerous changes in plastid transcript populations, many of which resemble those reported for the orthologous mutants in Arabidopsis. However, unanticipated reductions in the abundance of numerous transfer RNAs (tRNAs) dominated the microarray data and were validated on RNA gel blots. The magnitude of the deficiencies for several tRNAs was similar to that of the most severely affected messenger RNAs, with the loss of trnL-UAA being particularly severe. These findings suggest that PEP and its associated proteins are critical for the robust transcription of numerous plastid tRNAs and that this function is essential for the prodigious translation of plastid-encoded proteins that is required during the installation of the photosynthetic apparatus.

Patterns of genomic integration of nuclear chloroplast DNA fragments in plant species

The transfer of organelle DNA fragments to the nuclear genome is frequently observed in eukaryotes. These transfers are thought to play an important role in gene and genome evolution of eukaryotes. In plants, such transfers occur from plastid to nuclear [nuclear plastid DNAs (NUPTs)] and mitochondrial to nuclear (nuclear mitochondrial DNAs) genomes. The amount and genomic organization of organelle DNA fragments have been studied in model plant species, such as Arabidopsis thaliana and rice. At present, publicly available genomic data can be used to conduct such studies in non-model plants. In this study, we analysed the amount and genomic organization of NUPTs in 17 plant species for which genome sequences are available. The amount and distribution of NUPTs varied among the species. We also estimated the distribution of NUPTs according to the time of integration (relative age) by conducting sequence similarity analysis between NUPTs and the plastid genome. The age distributions suggested that the present genomic constitutions of NUPTs could be explained by the combination of the rapidly eliminated deleterious parts and few but constantly existing less deleterious parts.

Comparative chloroplast genomes of camellia species

BACKGROUND

Camellia, comprising more than 200 species, is a valuable economic commodity due to its enormously popular commercial products: tea leaves, flowers, and high-quality edible oils. It is the largest and most important genus in the family Theaceae. However, phylogenetic resolution of the species has proven to be difficult. Consequently, the interspecies relationships of the genus Camellia are still hotly debated. Phylogenomics is an attractive avenue that can be used to reconstruct the tree of life, especially at low taxonomic levels.

METHODOLOGY/PRINCIPAL FINDINGS

Seven complete chloroplast (cp) genomes were sequenced from six species representing different subdivisions of the genus Camellia using Illumina sequencing technology. Four junctions between the single-copy segments and the inverted repeats were confirmed and genome assemblies were validated by PCR-based product sequencing using 123 pairs of primers covering preliminary cp genome assemblies. The length of the Camellia cp genome was found to be about 157kb, which contained 123 unique genes and 23 were duplicated in the IR regions. We determined that the complete Camellia cp genome was relatively well conserved, but contained enough genetic differences to provide useful phylogenetic information. Phylogenetic relationships were analyzed using seven complete cp genomes of six Camellia species. We also identified rapidly evolving regions of the cp genome that have the potential to be used for further species identification and phylogenetic resolution.

CONCLUSIONS/SIGNIFICANCE

In this study, we wanted to determine if analyzing completely sequenced cp genomes could help settle these controversies of interspecies relationships in Camellia. The results demonstrate that cp genome data are beneficial in resolving species definition because they indicate that organelle-based "barcodes", can be established for a species and then used to unmask interspecies phylogenetic relationships. It reveals that phylogenomics based on cp genomes is an effective approach for achieving phylogenetic resolution between Camellia species.

Differential regulation of Arabidopsis plastid gene expression and RNA editing in non-photosynthetic tissues

RNA editing is one of the post-transcriptional processes that commonly occur in plant plastids and mitochondria. In Arabidopsis, 34 C-to-U RNA editing events, affecting transcripts of 18 plastid genes, have been identified. Here, we examined the editing and expression of these transcripts in different organs, and in green and non-green seedlings (etiolated, cia5-2, ispF and ispG albino mutants, lincomycin-, and norflurazon-treated). The editing efficiency of Arabidopsis plastid transcripts varies from site to site, and may be specifically regulated in different tissues. Steady state levels of plastid transcripts are low or undetectable in etiolated seedlings, but most editing sites are edited with efficiencies similar to those observed in green seedlings. By contrast, the editing of some sites is completely lost or significantly reduced in other non-green tissues; for instance, the editing of ndhB-149, ndhB-1255, and ndhD-2 is completely lost in roots and in lincomycin-treated seedlings. The editing of ndhD-2 is also completely lost in albino mutants and norflurazon-treated seedlings. However, matK-640 is completely edited, and accD-794, atpF-92, psbE-214, psbF-77, psbZ-50, and rps14-50 are completely or highly edited in both green and non-green tissues. In addition, the expression of nucleus-encoded RNA polymerase dependent transcripts is specifically induced by lincomycin, and the splicing of ndhB transcripts is significantly reduced in the albino mutants and inhibitor-treated seedlings. Our results indicate that plastid gene expression, and the splicing and editing of plastid transcripts are specifically and differentially regulated in various types of non-green tissues.

The requirement of WHIRLY1 for embryogenesis is dependent on genetic background in maize

Plastid gene expression is essential to embryogenesis in higher plants, but the underlying mechanism is obscure. Through molecular characterization of an embryo defective 16 (emb16) locus, here we report that the requirement of plastid translation for embryogenesis is dependent on the genetic background in maize (Zea mays). The emb16 mutation arrests embryogenesis at transition stage and allows the endosperm to develop largely normally. Molecular cloning reveals that Emb16 encodes WHIRLY1 (WHY1), a DNA/RNA binding protein that is required for genome stability and ribosome formation in plastids. Interestingly, the previous why1 mutant alleles (why1-1 and why1-2) do not affect embryogenesis, only conditions albino seedlings. The emb16 allele of why1 mutation is in the W22 genetic background. Crosses between emb16 and why1-1 heterozygotes resulted in both defective embryos and albino seedlings in the F1 progeny. Introgression of the emb16 allele from W22 into A188, B73, Mo17, Oh51a and the why1-1 genetic backgrounds yielded both defective embryos and albino seedlings. Similar results were obtained with two other emb mutants (emb12 and emb14) that are impaired in plastid protein translation process. These results indicate that the requirement of plastid translation for embryogenesis is dependent on genetic backgrounds, implying a mechanism of embryo lethality suppression in maize.

Comparative genome analysis and phylogenetic relationship of order Liliales insight from the complete plastid genome sequences of two Lilies (Lilium longiflorum and Alstroemeria aurea)

Monocots are one of the most diverse, successful and economically important clades of angiosperms. We attempt to analyse the complete plastid genome sequences of two lilies and their lengths were 152,793bp in Lilium longiflorum (Liliaceae) and 155,510bp in Alstroemeria aurea (Alstroemeriaceae). Phylogenetic analyses were performed for 28 taxa including major lineages of monocots using the sequences of 79 plastid genes for clarifying the phylogenetic relationship of the order Liliales. The sister relationship of Liliales and Asparagales-commelinids was improved with high resolution. Comparative analyses of inter-familial and inter-specific sequence variation were also carried out among three families of Liliaceae, Smilacaceae, and Alstroemeriaceae, and between two Lilium species of L. longflorum and L. superbum. Gene content and order were conserved in the order Liliales except infA loss in Smilax and Alstroemeria. IR boundaries were similar in IRa, however, IRb showed different extension patterns as JLB of Smilax and JSB in Alstroemeria. Ka/Ks ratio was high in matK among the pair-wise comparison of three families and the most variable genes were psaJ, ycf1, rpl32, rpl22, matK, and ccsA among the three families and rps15, rpoA, matK, and ndhF between Lilium.

Embryo defective12 encodes the plastid initiation factor 3 and is essential for embryogenesis in maize

Embryo-specific mutants in maize define a unique class of genetic loci that affect embryogenesis without a significant deleterious impact on endosperm development. Here we report the characterization of an embryo specific12 (emb12) mutant in maize. Embryogenesis in the emb12 mutants is arrested at or before transition stage. The mutant embryo at an early stage exhibits abnormal cell structure with increased vacuoles and dramatically reduced internal membrane organelles. In contrast, the mutant endosperm appears normal in morphology, cell structure, starch, lipid and protein accumulation. The Emb12 locus was cloned by transposon tagging and predicts a protein with a high similarity to prokaryotic translation initiation factor 3 (IF3). EMB12-GFP fusion analysis indicates that EMB12 is localized in plastids. The RNA in situ hybridization and protein immunohistochemical analyses indicate that a high level of Emb12 expression localizes in the embryo proper at early developmental stages and in the embryo axis at later stages. Western analysis indicates that plastid protein synthesis is impaired. These results indicate that Emb12 encodes the plastid IF3 which is essential for embryogenesis but not for endosperm development in maize.

Complete chloroplast genome of the genus Cymbidium: lights into the species identification, phylogenetic implications and population genetic analyses

BACKGROUND

Cymbidium orchids, including some 50 species, are the famous flowers, and they possess high commercial value in the floricultural industry. Furthermore, the values of different orchids are great differences. However, species identification is very difficult. To a certain degree, chloroplast DNA sequence data are a versatile tool for species identification and phylogenetic implications in plants. Different chloroplast loci have been utilized for evaluating phylogenetic relationships at each classification level among plant species, including at the interspecies and intraspecies levels. However, there is no evidence that a short sequence can distinguish all plant species from each other in order to infer phylogenetic relationships. Molecular markers derived from the complete chloroplast genome can provide effective tools for species identification and phylogenetic resolution.

RESULTS

The complete nucleotide sequences of eight individuals from a total of five Cymbidium species' chloroplast (cp) genomes were determined using Illumina sequencing technology of the total DNA via a combination of de novo and reference-guided assembly. The length of the Cymbidium cp genome is about 155 kb. The cp genomes contain 123 unique genes, and the IR regions contain 24 duplicates. Although the genomes, including genome structure, gene order and orientation, are similar to those of other orchids, they are not evolutionarily conservative. The cp genome of Cymbidium evolved moderately with more than 3% sequence divergence, which could provide enough information for phylogeny. Rapidly evolving chloroplast genome regions were identified and 11 new divergence hotspot regions were disclosed for further phylogenetic study and species identification in Orchidaceae.

CONCLUSIONS

Phylogenomic analyses were conducted using 10 complete chloroplast genomes from seven orchid species. These data accurately identified the individuals and established the phylogenetic relationships between the species. The results reveal that phylogenomics based on organelle genome sequencing lights the species identification-organelle-scale "barcodes", and is also an effective approach for studying whole populations and phylogenetic characteristics of Cymbidium.

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.

Ferredoxin:thioredoxin reductase (FTR) links the regulation of oxygenic photosynthesis to deeply rooted bacteria

Uncovered in studies on photosynthesis 35 years ago, redox regulation has been extended to all types of living cells. We understand a great deal about the occurrence, function, and mechanism of action of this mode of regulation, but we know little about its origin and its evolution. To help fill this gap, we have taken advantage of available genome sequences that make it possible to trace the phylogenetic roots of members of the system that was originally described for chloroplasts-ferredoxin, ferredoxin:thioredoxin reductase (FTR), and thioredoxin as well as target enzymes. The results suggest that: (1) the catalytic subunit, FTRc, originated in deeply rooted microaerophilic, chemoautotrophic bacteria where it appears to function in regulating CO(2) fixation by the reverse citric acid cycle; (2) FTRc was incorporated into oxygenic photosynthetic organisms without significant structural change except for addition of a variable subunit (FTRv) seemingly to protect the Fe-S cluster against oxygen; (3) new Trxs and target enzymes were systematically added as evolution proceeded from bacteria through the different types of oxygenic photosynthetic organisms; (4) an oxygenic type of regulation preceded classical light-dark regulation in the regulation of enzymes of CO(2) fixation by the Calvin-Benson cycle; (5) FTR is not universally present in oxygenic photosynthetic organisms, and in certain early representatives is seemingly functionally replaced by NADP-thioredoxin reductase; and (6) FTRc underwent structural diversification to meet the ecological needs of a variety of bacteria and archaea.

Plastome sequences of two New World bamboos--Arundinaria gigantea and Cryptochloa strictiflora (Poaceae)--extend phylogenomic understanding of Bambusoideae

PREMISE OF THE STUDY

Two New World species of Bambusoideae, Arundinaria gigantea and Crytpochloa strictiflora, were investigated in a phylogenomic context. Complete plastome sequences have been previously determined and analyzed for nine bambusoid species that exclusively represent Old World lineages. The addition of New World species provides more complete information on relationships within Bambusoideae. •

METHODS

Plastomes from A. gigantea and C. strictiflora were sequenced using Sanger methods. Phylogenomic and divergence estimate analyses were conducted on both species with 23 other Poaceae. •

KEY RESULTS

Phylogenomic and divergence analyses suggested that A. gigantea diverged from within Arundinarieae between 1.94-3.92 mya and that C. strictiflora diverged as the sister to tropical woody species between 24.83 and 40.22 mya. These results are correlated with modern relative diversities in the two lineages. •

CONCLUSIONS

The two New World bamboos show unique plastome features accumulated and maintained in biogeographic isolation from Old World taxa. The overall evidence for A. gigantea is consistent with recent dispersal, and that for C. strictiflora is consistent with vicariance.

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.

Rogowsky P. PPR8522 encodes a chloroplast-targeted pentatricopeptide repeat protein necessary for maize embryogenesis and vegetative development

The pentatricopeptide repeat (PPR) domain is an RNA binding domain allowing members of the PPR superfamily to participate in post-transcriptional processing of organellar RNA. Loss of PPR8522 from maize (Zea mays) confers an embryo-specific (emb) phenotype. The emb8522 mutation was isolated in an active Mutator (Mu) population and co-segregation analysis revealed that it was tightly linked to a MuDR insertion in the first exon of PPR8522. Independent evidence that disruption of PPR8522 caused the emb phenotype was provided by fine mapping to a region of 116kb containing no other gene than PPR8522 and complementation of the emb8522 mutant by a PPR8522 cDNA. The deduced PPR8522 amino acid sequence of 832 amino acids contains 10 PPR repeats and a chloroplast target peptide, the function of which was experimentally demonstrated by transient expression in Nicotiana benthamiana. Whereas mutant endosperm is apparently normal, mutant embryos deviate from normal development as early as 3 days after pollination, are reduced in size, exhibit more or less severe morphological aberrations depending on the genetic background, and generally do not germinate. The emb8522 mutation is the first to associate the loss of a PPR gene with an embryo-lethal phenotype in maize. Analyses of mutant plantlets generated by embryo-rescue experiments indicate that emb8522 also affects vegetative plant growth and chloroplast development. The loss of chloroplast transcription dependent on plastid-encoded RNA polymerase is the likely cause for the lack of an organized thylakoid network and an albino, seedling-lethal phenotype.

Identification of cis-elements conferring high levels of gene expression in non-green plastids

Although our knowledge about the mechanisms of gene expression in chloroplasts has increased substantially over the past decades, next to nothing is known about the signals and factors that govern expression of the plastid genome in non-green tissues. Here we report the development of a quantitative method suitable for determining the activity of cis-acting elements for gene expression in non-green plastids. The in vivo assay is based on stable transformation of the plastid genome and the discovery that root length upon seedling growth in the presence of the plastid translational inhibitor kanamycin is directly proportional to the expression strength of the resistance gene nptII in transgenic tobacco plastids. By testing various combinations of promoters and translation initiation signals, we have used this experimental system to identify cis-elements that are highly active in non-green plastids. Surprisingly, heterologous expression elements from maize plastids were significantly more efficient in conferring high expression levels in root plastids than homologous expression elements from tobacco. Our work has established a quantitative method for characterization of gene expression in non-green plastid types, and has led to identification of cis-elements for efficient plastid transgene expression in non-green tissues, which are valuable tools for future transplastomic studies in basic and applied research.

The influence of identity elements on the aminoacylation of tRNA(Arg) by plant and Escherichia coli arginyl-tRNA synthetases

Identity elements determine the accurate recognition between tRNAs and aminoacyl-tRNA synthetases. The arginine system from yeast and Escherichia coli has been studied extensively in the past. However, information about the enzymes from higher eukaryotes is limited and plant aminoacyl-tRNA synthetases have been largely ignored in this respect. We have designed in vitro tRNA transcripts, based on the soybean tRNA(Arg) primary structure, aiming to investigate its specific aminoacylation by two recombinant plant arginyl-tRNA synthetases and to compare this with the enzyme from E. coli. Identity elements at positions 20 and 35 in plants parallel those previously established for bacteria. Cryptic identity elements in the plant system that are not revealed within a tRNA(Arg) consensus sequence compiled from isodecoders corresponding to nine distinct cytoplasmic, mitochondrial and plastid isoaccepting sequences are located in the acceptor stem. Additionally, it has been shown that U20a and A38 are essential for a fully efficient cognate E. coli arginylation, whereas, for the plant arginyl-tRNA synthetases, these bases can be replaced by G20a and C38 with full retention of activity. G10, a constituent of the 10:25:45 tertiary interaction, is essential for both plant and E. coli activity. Amino acid recognition in terms of discriminating between arginine and canavanine by the arginyl-tRNA synthetase from both kingdoms may be manipulated by changes at different sites within the tRNA structure.

Choosing the best plant for the job: a cost-effective assay to prescreen ancient plant remains destined for shotgun sequencing

DNA extracted from ancient plant remains almost always contains a mixture of endogenous (that is, derived from the plant) and exogenous (derived from other sources) DNA. The exogenous 'contaminant' DNA, chiefly derived from microorganisms, presents significant problems for shotgun sequencing. In some samples, more than 90% of the recovered sequences are exogenous, providing limited data relevant to the sample. However, other samples have far less contamination and subsequently yield much more useful data via shotgun sequencing. Given the investment required for high-throughput sequencing, whenever multiple samples are available, it is most economical to sequence the least contaminated sample. We present an assay based on quantitative real-time PCR which estimates the relative amounts of fungal and bacterial DNA in a sample in comparison to the endogenous plant DNA. Given a collection of contextually-similar ancient plant samples, this low cost assay aids in selecting the best sample for shotgun sequencing.

Activity, polypeptide and gene identification of thylakoid Ndh complex in trees: potential physiological relevance of fluorescence assays

Three evergreen (Laurus nobilis, Viburnum tinus and Thuja plicata) and two autumnal abscission deciduous trees (Cydonia oblonga and Prunus domestica) have been investigated for the presence (zymogram and immunodetection) and functionality (post-illumination chlorophyll fluorescence) of the thylakoid Ndh complex. The presence of encoding ndh genes has also been investigated in T. plicata. Western assays allowed tentative identification of zymogram NADH dehydrogenase bands corresponding to the Ndh complex after native electrophoresis of solubilized fractions from L. nobilis, V. tinus, C. oblonga and P. domestica leaves, but not in those of T. plicata. However, Ndh subunits were detected after SDS-PAGE of thylakoid solubilized proteins of T. plicata. The leaves of the five plants showed the post-illumination chlorophyll fluorescence increase dependent on the presence of active Ndh complex. The fluorescence increase was higher in autumn in deciduous, but not in evergreen trees, which suggests that the thylakoid Ndh complex could be involved in autumnal leaf senescence. Two ndhB genes were sequenced from T. plicata that differ at the 350 bp 3' end sequence. Comparison with the mRNA revealed that ndhB genes have a 707-bp type II intron between exons 1 (723 bp) and 2 (729 bp) and that the UCA 259th codon is edited to UUA in mRNA. Phylogenetically, the ndhB genes of T. plicata group close to those of Metasequoia, Cryptomeria, Taxodium, Juniperus and Widdringtonia in the cupresaceae branch and are 5' end shortened by 18 codons with respect to that of angiosperms.

Analysis of complete nucleotide sequences of 12 Gossypium chloroplast genomes: origin and evolution of allotetraploids

BACKGROUND

Cotton (Gossypium spp.) is a model system for the analysis of polyploidization. Although ascertaining the donor species of allotetraploid cotton has been intensively studied, sequence comparison of Gossypium chloroplast genomes is still of interest to understand the mechanisms underlining the evolution of Gossypium allotetraploids, while it is generally accepted that the parents were A- and D-genome containing species. Here we performed a comparative analysis of 13 Gossypium chloroplast genomes, twelve of which are presented here for the first time.

METHODOLOGY/PRINCIPAL FINDINGS

The size of 12 chloroplast genomes under study varied from 159,959 bp to 160,433 bp. The chromosomes were highly similar having >98% sequence identity. They encoded the same set of 112 unique genes which occurred in a uniform order with only slightly different boundary junctions. Divergence due to indels as well as substitutions was examined separately for genome, coding and noncoding sequences. The genome divergence was estimated as 0.374% to 0.583% between allotetraploid species and A-genome, and 0.159% to 0.454% within allotetraploids. Forty protein-coding genes were completely identical at the protein level, and 20 intergenic sequences were completely conserved. The 9 allotetraploids shared 5 insertions and 9 deletions in whole genome, and 7-bp substitutions in protein-coding genes. The phylogenetic tree confirmed a close relationship between allotetraploids and the ancestor of A-genome, and the allotetraploids were divided into four separate groups. Progenitor allotetraploid cotton originated 0.43-0.68 million years ago (MYA).

CONCLUSION

Despite high degree of conservation between the Gossypium chloroplast genomes, sequence variations among species could still be detected. Gossypium chloroplast genomes preferred for 5-bp indels and 1-3-bp indels are mainly attributed to the SSR polymorphisms. This study supports that the common ancestor of diploid A-genome species in Gossypium is the maternal source of extant allotetraploid species and allotetraploids have a monophyletic origin. G. hirsutum AD1 lineages have experienced more sequence variations than other allotetraploids in intergenic regions. The available complete nucleotide sequences of 12 Gossypium chloroplast genomes should facilitate studies to uncover the molecular mechanisms of compartmental co-evolution and speciation of Gossypium allotetraploids.

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 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.