The complete sequence of the rice (Oryza sativa) chloroplast genome: intermolecular recombination between distinct tRNA genes accounts for a major plastid DNA inversion during the evolution of the cereals

Chloroplast genome of Aconitum barbatum var. puberulum (Ranunculaceae) derived from CCS reads using the PacBio RS platform

The chloroplast genome (cp genome) of Aconitum barbatum var. puberulum was sequenced using the third-generation sequencing platform based on the single-molecule real-time (SMRT) sequencing approach. To our knowledge, this is the first reported complete cp genome of Aconitum, and we anticipate that it will have great value for phylogenetic studies of the Ranunculaceae family. In total, 23,498 CCS reads and 20,685,462 base pairs were generated, the mean read length was 880 bp, and the longest read was 2,261 bp. Genome coverage of 100% was achieved with a mean coverage of 132× and no gaps. The accuracy of the assembled genome is 99.973%; the assembly was validated using Sanger sequencing of six selected genes from the cp genome. The complete cp genome of A. barbatum var. puberulum is 156,749 bp in length, including a large single-copy region of 87,630 bp and a small single-copy region of 16,941 bp separated by two inverted repeats of 26,089 bp. The cp genome contains 130 genes, including 84 protein-coding genes, 34 tRNA genes and eight rRNA genes. Four forward, five inverted and eight tandem repeats were identified. According to the SSR analysis, the longest poly structure is a 20-T repeat. Our results presented in this paper will facilitate the phylogenetic studies and molecular authentication on Aconitum.

The rice ALS3 encoding a novel pentatricopeptide repeat protein is required for chloroplast development and seedling growth

BACKGROUND

Pentatricopeptide repeat (PPR) proteins play essential roles in modulating the expression of organelle genes and have expanded greatly in higher plants. However, molecular mechanisms of most rice PPR genes remain unclear.

RESULTS

In this study, a new rice PPR mutant, asl3 (albino seedling lethality3) exhibits an albino lethal phenotype at the seedling stage. This albino phenotype was associated with altered photosynthetic-pigment and chloroplast development. Map-based cloning showed that ASL3 encodes a novel rice PPR protein with 10 tandem PPR motifs, which localizes to the chloroplast. ASL3 showed tissue-specific expression, as it was highly expressed in the chlorenchyma, but expressed at much lower levels in roots and panicles. RNAi of ASL3 confirmed that ASL3 plays an essential role in the early development and chloroplast development in rice. Moreover, expression analysis revealed that the asl3 mutation severely affected the transcriptional levels of important genes associated with plastid translation machinery and photosynthesis, which may impair photosynthesis and finally led to the seedling death in asl3 mutant. These results evidenced the important role of ASL3 in the early development of rice, especially chloroplast development.

CONCLUSIONS

The ASL3 gene encoded a novel chloroplast-targeted PPR protein with 10 tandem PPR motifs in rice. Disruption of the ASL3 would lead to a defective chloroplast and seedling lethality, and affected expression levels of genes associated with chloroplast development and photosynthesis at early leaf stage of rice.

Processing of the 5'-UTR and existence of protein factors that regulate translation of tobacco chloroplast psbN mRNA

The chloroplast psbB operon includes five genes encoding photosystem II and cytochrome b 6 /f complex components. The psbN gene is located on the opposite strand. PsbN is localized in the thylakoid and is present even in the dark, although its level increases upon illumination and then decreases. However, the translation mechanism of the psbN mRNA remains unclear. Using an in vitro translation system from tobacco chloroplasts and a green fluorescent protein as a reporter protein, we show that translation occurs from a tobacco primary psbN 5'-UTR of 47 nucleotides (nt). Unlike many other chloroplast 5'-UTRs, the psbN 5'-UTR has two processing sites, at -39 and -24 upstream from the initiation site. Processing at -39 enhanced the translation rate fivefold. In contrast, processing at -24 did not affect the translation rate. These observations suggest that the two distinct processing events regulate, at least in part, the level of PsbN during development. The psbN 5'-UTR has no Shine-Dalgarno (SD)-like sequence. In vitro translation assays with excess amounts of the psbN 5'-UTR or with deleted psbN 5'-UTR sequences demonstrated that protein factors are required for translation and that their binding site is an 18 nt sequence in the 5'-UTR. Mobility shift assays using 10 other chloroplast 5'-UTRs suggested that common or similar proteins are involved in translation of a set of mRNAs lacking SD-like sequences.

The whole chloroplast genome of wild rice (Oryza australiensis)

The whole chloroplast genome of wild rice (Oryza australiensis) is characterized in this study. The genome size is 135,224  bp, exhibiting a typical circular structure including a pair of 25,776  bp inverted repeats (IRa,b) separated by a large single-copy region (LSC) of 82,212  bp and a small single-copy region (SSC) of 12,470  bp. The overall GC content of the genome is 38.95%. 110 unique genes were annotated, including 76 protein-coding genes, 4 ribosomal RNA genes, and 30t RNA genes. Among these, 18 are duplicated in the inverted repeat regions, 13 genes contain one intron, and 2 genes (rps12 and ycf3) have two introns.

The rice OsDG2 encoding a glycine-rich protein is involved in the regulation of chloroplast development during early seedling stage

OsDG2 gene encoded a novel chloroplast-targeted GRP in rice. Disruption of the OsDG2 would lead to delayed greening phenotype and affected expression levels of genes associated with chloroplast development at early leaf stage of rice. Glycine-rich proteins (GRPs) participate in various biological processes in plants. However, the evidence of GRPs involved in chloroplast development in plants is quite limited. In this study, we identified a rice GRP gene mutant named osdg2 (O ryza s ativa d elayed g reening 2), which exhibits delayed greening phenotype characterized as bright yellow leaves before the three-leaf stage and thereafter turns to normal green. Further study showed that the mutant phenotype was consistent with changes in chlorophyll content and chloroplast development. The rice OsDG2 gene, encoding a novel GRP protein, was located on chromosome 2 through map-based cloning method and confirmed by molecular complementation tests. Subcellular localization results showed that OsDG2 was targeted in chloroplasts. In addition, the OsDG2 transcripts were highly expressed in leaves and undetectable in other tissues, showing the tissue-specific expression. In osdg2 mutant, the expression levels of most genes associated with chloroplast development were severely decreased in the 3rd leaves, but almost recovered to wild-type level in the 4th leaves. Our findings indicated that the nuclear-encoded OsDG2 plays important roles in chloroplast development at early leaf stage of rice.

Chloroplast genome evolution in early diverged leptosporangiate ferns

In this study, the chloroplast (cp) genome sequences from three early diverged leptosporangiate ferns were completed and analyzed in order to understand the evolution of the genome of the fern lineages. The complete cp genome sequence of Osmunda cinnamomea (Osmundales) was 142,812 base pairs (bp). The cp genome structure was similar to that of eusporangiate ferns. The gene/intron losses that frequently occurred in the cp genome of leptosporangiate ferns were not found in the cp genome of O. cinnamomea. In addition, putative RNA editing sites in the cp genome were rare in O. cinnamomea, even though the sites were frequently predicted to be present in leptosporangiate ferns. The complete cp genome sequence of Diplopterygium glaucum (Gleicheniales) was 151,007 bp and has a 9.7 kb inversion between the trnL-CAA and trnVGCA genes when compared to O. cinnamomea. Several repeated sequences were detected around the inversion break points. The complete cp genome sequence of Lygodium japonicum (Schizaeales) was 157,142 bp and a deletion of the rpoC1 intron was detected. This intron loss was shared by all of the studied species of the genus Lygodium. The GC contents and the effective numbers of codons (ENCs) in ferns varied significantly when compared to seed plants. The ENC values of the early diverged leptosporangiate ferns showed intermediate levels between eusporangiate and core leptosporangiate ferns. However, our phylogenetic tree based on all of the cp gene sequences clearly indicated that the cp genome similarity between O. cinnamomea (Osmundales) and eusporangiate ferns are symplesiomorphies, rather than synapomorphies. Therefore, our data is in agreement with the view that Osmundales is a distinct early diverged lineage in the leptosporangiate ferns.

Comparative analysis of complete chloroplast genome sequence and inversion variation in Lasthenia burkei (Madieae, Asteraceae)

PREMISE OF THE STUDY

Complete chloroplast genome studies can help resolve relationships among large, complex plant lineages such as Asteraceae. We present the first whole plastome from the Madieae tribe and compare its sequence variation to other chloroplast genomes in Asteraceae.

METHODS

We used high throughput sequencing to obtain the Lasthenia burkei chloroplast genome. We compared sequence structure and rates of molecular evolution in the small single copy (SSC), large single copy (LSC), and inverted repeat (IR) regions to those for eight Asteraceae accessions and one Solanaceae accession.

KEY RESULTS

The chloroplast sequence of L. burkei is 150 746 bp and contains 81 unique protein coding genes and 4 coding ribosomal RNA sequences. We identified three major inversions in the L. burkei chloroplast, all of which have been found in other Asteraceae lineages, and a previously unreported inversion in Lactuca sativa. Regions flanking inversions contained tRNA sequences, but did not have particularly high G + C content. Substitution rates varied among the SSC, LSC, and IR regions, and rates of evolution within each region varied among species. Some observed differences in rates of molecular evolution may be explained by the relative proportion of coding to noncoding sequence within regions.

CONCLUSIONS

Rates of molecular evolution vary substantially within and among chloroplast genomes, and major inversion events may be promoted by the presence of tRNAs. Collectively, these results provide insight into different mechanisms that may promote intramolecular recombination and the inversion of large genomic regions in the plastome.

Chloroplast genes as genetic markers for inferring patterns of change, maternal ancestry and phylogenetic relationships among Eleusine species

Assessment of phylogenetic relationships is an important component of any successful crop improvement programme, as wild relatives of the crop species often carry agronomically beneficial traits. Since its domestication in East Africa, Eleusine coracana (2n = 4x = 36), a species belonging to the genus Eleusine (x = 8, 9, 10), has held a prominent place in the semi-arid regions of India, Nepal and Africa. The patterns of variation between the cultivated and wild species reported so far and the interpretations based upon them have been considered primarily in terms of nuclear events. We analysed, for the first time, the phylogenetic relationship between finger millet (E. coracana) and its wild relatives by species-specific chloroplast deoxyribonucleic acid (cpDNA) polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) and chloroplast simple sequence repeat (cpSSR) markers/sequences. Restriction fragment length polymorphism of the seven amplified chloroplast genes/intergenic spacers (trnK, psbD, psaA, trnH-trnK, trnL-trnF, 16S and trnS-psbC), nucleotide sequencing of the chloroplast trnK gene and chloroplast microsatellite polymorphism were analysed in all nine known species of Eleusine. The RFLP of all seven amplified chloroplast genes/intergenic spacers and trnK gene sequences in the diploid (2n = 16, 18, 20) and allotetraploid (2n = 36, 38) species resulted in well-resolved phylogenetic trees with high bootstrap values. Eleusine coracana, E. africana, E. tristachya, E. indica and E. kigeziensis did not show even a single change in restriction site. Eleusine intermedia and E. floccifolia were also shown to have identical cpDNA fragment patterns. The cpDNA diversity in Eleusine multiflora was found to be more extensive than that of the other eight species. The trnK gene sequence data complemented the results obtained by PCR-RFLP. The maternal lineage of all three allotetraploid species (AABB, AADD) was the same, with E. indica being the maternal diploid progenitor species. The markers specific to certain species were also identified.

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.

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.

Disruption of the rice plastid ribosomal protein s20 leads to chloroplast developmental defects and seedling lethality

Plastid ribosomal proteins (PRPs) are essential for ribosome biogenesis, plastid protein biosynthesis, chloroplast differentiation, and early chloroplast development. This study identifies the first rice PRP mutant, asl1 (albino seedling lethality1), which exhibits an albino lethal phenotype at the seedling stage. This albino phenotype was associated with altered chlorophyll (Chl) content and chloroplast development. Map-based cloning revealed that ASL1 encodes PRP S20 (PRPS20), which localizes to the chloroplast. ASL1 showed tissue-specific expression, as it was highly expressed in plumule and young seedlings but expressed at much lower levels in other tissues. In addition, ASL1 expression was regulated by light. The transcript levels of nuclear genes for Chl biosynthesis and chloroplast development were strongly affected in asl1 mutants; transcripts of some plastid genes for photosynthesis were undetectable. Our findings indicate that nuclear-encoded PRPS20 plays an important role in chloroplast development in rice.

Plastid trnF pseudogenes are present in Jaltomata, the sister genus of Solanum (Solanaceae): molecular evolution of tandemly repeated structural mutations

Extensive gene duplication arranged in a tandem array is rare in the plastome of embryophytes. Interestingly, we found pseudogene copies of the trnF gene in the genus Jaltomata, the sister genus of Solanum where such gene duplication has been previously reported. In each Jaltomata sequence available we found two pseudogene copies in close 5'-proximity to the original functional gene. The size of each pseudogene copy ranged between 17 and 48 bp and the anticodon domain was identified as the most conserved element. A common ATT(G)n motif is particularly interesting and its modifications were found to border the 3' of the duplicated regions. Other motifs were partial residues, or entire parts of the T- and D-domains, and both domains proved to be variable in length among the pseudogenes identified. The residues of the 3' and 5' acceptor stem were not found among the copies. We further compared the newly discovered copies of Jaltomata with those ones previously described from Solanum and inferred phylogenetic relationships of the copies aligned. The evolution of Solanum copies, in contrast to Jaltomata, is hard to explain as resulting only in parsimonious changes since reticulate evolutionary patterns were detected among the copies. The dynamic evolutionary patterns of Solanum might be explained by possible inter- or intrachromosomal recombination.

The complete chloroplast genome of banana (Musa acuminata, Zingiberales): insight into plastid monocotyledon evolution

BACKGROUND

Banana (genus Musa) is a crop of major economic importance worldwide. It is a monocotyledonous member of the Zingiberales, a sister group of the widely studied Poales. Most cultivated bananas are natural Musa inter-(sub-)specific triploid hybrids. A Musa acuminata reference nuclear genome sequence was recently produced based on sequencing of genomic DNA enriched in nucleus.

METHODOLOGY/PRINCIPAL FINDINGS

The Musa acuminata chloroplast genome was assembled with chloroplast reads extracted from whole-genome-shotgun sequence data. The Musa chloroplast genome is a circular molecule of 169,972 bp with a quadripartite structure containing two single copy regions, a Large Single Copy region (LSC, 88,338 bp) and a Small Single Copy region (SSC, 10,768 bp) separated by Inverted Repeat regions (IRs, 35,433 bp). Two forms of the chloroplast genome relative to the orientation of SSC versus LSC were found. The Musa chloroplast genome shows an extreme IR expansion at the IR/SSC boundary relative to the most common structures found in angiosperms. This expansion consists of the integration of three additional complete genes (rps15, ndhH and ycf1) and part of the ndhA gene. No such expansion has been observed in monocots so far. Simple Sequence Repeats were identified in the Musa chloroplast genome and a new set of Musa chloroplastic markers was designed.

CONCLUSION

The complete sequence of M. acuminata ssp malaccensis chloroplast we reported here is the first one for the Zingiberales order. As such it provides new insight in the evolution of the chloroplast of monocotyledons. In particular, it reinforces that IR/SSC expansion has occurred independently several times within monocotyledons. The discovery of new polymorphic markers within Musa chloroplast opens new perspectives to better understand the origin of cultivated triploid bananas.

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.

Diversification and genetic differentiation of cultivated melon inferred from sequence polymorphism in the chloroplast genome

Molecular analysis encouraged discovery of genetic diversity and relationships of cultivated melon (Cucumis melo L.). We sequenced nine inter- and intra-genic regions of the chloroplast genome, about 5500 bp, using 60 melon accessions and six reference accessions of wild species of Cucumis to show intra-specific variation of the chloroplast genome. Sequence polymorphisms were detected among melon accessions and other Cucumis species, indicating intra-specific diversification of the chloroplast genome. Melon accessions were classified into three subclusters by cytoplasm type and then into 12 subgroups. Geographical origin and seed size also differed between the three subclusters. Subcluster Ia contained small-seed melon from Southern Africa and South and East Asia and subcluster Ib mainly consisted of large-seed melon from northern Africa, Europe and USA. Melon accessions of subcluster Ic were only found in West, Central and Southern Africa. Our results indicated that European melon groups and Asian melon groups diversified independently and shared the same maternal lineage with northern African large-seed melon and Southern African small-seed melon, respectively. Cultivated melon of subcluster Ic may have been domesticated independently in Africa. The presence of 11 cytoplasm types in Africa strongly supported African origin of cultivated melon and indicated the importance of germplasm from Africa.

Plastome Sequence Determination and Comparative Analysis for Members of the Lolium-Festuca Grass Species Complex

Chloroplast genome sequences are of broad significance in plant biology, due to frequent use in molecular phylogenetics, comparative genomics, population genetics, and genetic modification studies. The present study used a second-generation sequencing approach to determine and assemble the plastid genomes (plastomes) of four representatives from the agriculturally important Lolium-Festuca species complex of pasture grasses (Lolium multiflorum, Festuca pratensis, Festuca altissima, and Festuca ovina). Total cellular DNA was extracted from either roots or leaves, was sequenced, and the output was filtered for plastome-related reads. A comparison between sources revealed fewer plastome-related reads from root-derived template but an increase in incidental bacterium-derived sequences. Plastome assembly and annotation indicated high levels of sequence identity and a conserved organization and gene content between species. However, frequent deletions within the F. ovina plastome appeared to contribute to a smaller plastid genome size. Comparative analysis with complete plastome sequences from other members of the Poaceae confirmed conservation of most grass-specific features. Detailed analysis of the rbcL-psaI intergenic region, however, revealed a "hot-spot" of variation characterized by independent deletion events. The evolutionary implications of this observation are discussed. The complete plastome sequences are anticipated to provide the basis for potential organelle-specific genetic modification of pasture grasses.

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.

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.

Chloroplast DNA polymorphism and evolutional relationships between Asian cultivated rice (Oryza sativa) and its wild relatives (O. rufipogon)

We analyzed chloroplast DNA (cpDNA) polymorphism and phylogenic relationships between 6 typical indica rice, 4 japonica rice, 8 javanica rice, and 12 Asian common wild rice (Oryza rufipogon) strains collected from different latitudes in China by comparing polymorphism at 9 highly variable regions. One hundred and forty-four polymorphic bases were detected. The O. rufipogon samples had 117 polymorphic bases, showing rich genetic diversity. One hundred and thirty-one bases at 13 sites were identified with indica/japonica characteristics; they showed differences between the indica and japonica subspecies at these sites. The javanica strains and japonica shared similar bases at these 131 polymorphic sites, suggesting that javanica is closely related to japonica. On the basis of length analyses of the open reading frame (ORF)100 and (ORF)29-tRNA-Cys(GCA) (TrnC(GCA)) fragments, the O. rufipogon strains were classified into indica/japonica subgroups, which was consistent with the results of the phylogenic tree assay based on concatenated datasets. These results indicated that differences in indica and japonica also exist in the cpDNA genome of the O. rufipogon strains. However, these differences demonstrated a certain degree of primitiveness and incompleteness, as an O. rufipogon line may show different indica/ japonica attributes at different sites. Consequently, O. rufipogon cannot be simply classified into the indica/japonica types as O. sativa. Our data support the hypothesis that Asian cultivated rice, O. indica and O. japonica, separately evolved from Asian common wild rice (O. rufipogon) strains, which have different indica-japonica differentiation trends.

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 organelle genomes of Hassawi rice (Oryza sativa L.) and its hybrid in saudi arabia: genome variation, rearrangement, and origins

Hassawi rice (Oryza sativa L.) is a landrace adapted to the climate of Saudi Arabia, characterized by its strong resistance to soil salinity and drought. Using high quality sequencing reads extracted from raw data of a whole genome sequencing project, we assembled both chloroplast (cp) and mitochondrial (mt) genomes of the wild-type Hassawi rice (Hassawi-1) and its dwarf hybrid (Hassawi-2). We discovered 16 InDels (insertions and deletions) but no SNP (single nucleotide polymorphism) is present between the two Hassawi cp genomes. We identified 48 InDels and 26 SNPs in the two Hassawi mt genomes and a new type of sequence variation, termed reverse complementary variation (RCV) in the rice cp genomes. There are two and four RCVs identified in Hassawi-1 when compared to 93–11 (indica) and Nipponbare (japonica), respectively. Microsatellite sequence analysis showed there are more SSRs in the genic regions of both cp and mt genomes in the Hassawi rice than in the other rice varieties. There are also large repeats in the Hassawi mt genomes, with the longest length of 96,168 bp and 96,165 bp in Hassawi-1 and Hassawi-2, respectively. We believe that frequent DNA rearrangement in the Hassawi mt and cp genomes indicate ongoing dynamic processes to reach genetic stability under strong environmental pressures. Based on sequence variation analysis and the breeding history, we suggest that both Hassawi-1 and Hassawi-2 originated from the Indonesian variety Peta since genetic diversity between the two Hassawi cultivars is very low albeit an unknown historic origin of the wild-type Hassawi rice.

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.

Evolution and biogeography of the slipper orchids: Eocene vicariance of the conduplicate genera in the Old and New World Tropics

Intercontinental disjunctions between tropical regions, which harbor two-thirds of the flowering plants, have drawn great interest from biologists and biogeographers. Most previous studies on these distribution patterns focused on woody plants, and paid little attention to herbs. The Orchidaceae is one of the largest families of angiosperms, with a herbaceous habit and a high species diversity in the Tropics. Here we investigate the evolutionary and biogeographical history of the slipper orchids, which represents a monophyletic subfamily (Cypripedioideae) of the orchid family and comprises five genera that are disjunctly distributed in tropical to temperate regions. A relatively well-resolved and highly supported phylogeny of slipper orchids was reconstructed based on sequence analyses of six maternally inherited chloroplast and two low-copy nuclear genes (LFY and ACO). We found that the genus Cypripedium with a wide distribution in the northern temperate and subtropical zones diverged first, followed by Selenipedium endemic to South America, and finally conduplicate-leaved genera in the Tropics. Mexipedium and Phragmipedium from the neotropics are most closely related, and form a clade sister to Paphiopedilum from tropical Asia. According to molecular clock estimates, the genus Selenipedium originated in Palaeocene, while the most recent common ancestor of conduplicate-leaved slipper orchids could be dated back to the Eocene. Ancestral area reconstruction indicates that vicariance is responsible for the disjunct distribution of conduplicate slipper orchids in palaeotropical and neotropical regions. Our study sheds some light on mechanisms underlying generic and species diversification in the orchid family and tropical disjunctions of herbaceous plant groups. In addition, we suggest that the biogeographical study should sample both regional endemics and their widespread relatives.

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.

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.

Characterization of the chloroplast genome sequence of oil palm (Elaeis guineensis Jacq.)

Oil palm (Elaeis guineensis Jacq.) is an economically important crop, which is grown for oil production. To better understand the molecular basis of oil palm chloroplasts, we characterized the complete chloroplast (cp) genome sequence obtained from 454 pyrosequencing. The oil palm cp genome is 156,973 bp in length consisting of a large single-copy region of 85,192 bp flanked on each side by inverted repeats of 27,071 bp with a small single-copy region of 17,639 bp joining the repeats. The genome contains 112 unique genes: 79 protein-coding genes, 4 ribosomal RNA genes and 29 tRNA genes. By aligning the cp genome sequence with oil palm cDNA sequences, we observed 18 non-silent and 10 silent RNA editing events among 19 cp protein-coding genes. Creation of an initiation codon by RNA editing in rpl2 has been reported in several monocots and was also found in the oil palm cp genome. Fifty common chloroplast protein-coding genes from 33 plant taxa were used to construct ML and MP phylogenetic trees. Their topologies are similar and strongly support for the position of E. guineensis as the sister of closely related species Phoenix dactylifera in Arecaceae (palm families) of monocot subtrees.

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.

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.

Sequencing and analysis of plastid genome in mycoheterotrophic orchid Neottia nidus-avis

Plastids are the semiautonomous organelles that possess their own genome inherited from the cyanobacterial ancestor. The primary function of plastids is photosynthesis so the structure and evolution of plastid genomes are extensively studied in photosynthetic plants. In contrast, little is known about the plastomes of nonphotosynthetic species. In higher plants, plastid genome sequences are available for only three strictly nonphotosynthetic species, the liverwort Aneura mirabilis and two flowering plants, Epifagus virginiana and Rhizanthella gardneri. We report here the complete sequence of a plastid genome of nonphotosynthetic mycoheterotrophic orchid Neottia nidus-avis, determined using 454 pyrosequencing technology. It was found to be reduced in both genome size and gene content; this reduction is however not as drastic as in the other nonphotosynthetic orchid, R. gardneri. Neottia plastome lacks all genes encoding photosynthetic proteins, RNA polymerase subunits but retains most genes of translational apparatus. Those genes that are retained have an increased rate of both synonymous and nonsynonymous substitutions but do not exhibit relaxation of purifying selection either in Neottia or in Rhizanthella.

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.

Chloroplast genome variation in upland and lowland switchgrass

Switchgrass (Panicum virgatum L.) exists at multiple ploidies and two phenotypically distinct ecotypes. To facilitate interploidal comparisons and to understand the extent of sequence variation within existing breeding pools, two complete switchgrass chloroplast genomes were sequenced from individuals representative of the upland and lowland ecotypes. The results demonstrated a very high degree of conservation in gene content and order with other sequenced plastid genomes. The lowland ecotype reference sequence (Kanlow Lin1) was 139,677 base pairs while the upland sequence (Summer Lin2) was 139,619 base pairs. Alignments between the lowland reference sequence and short-read sequence data from existing sequence datasets identified as either upland or lowland confirmed known polymorphisms and indicated the presence of other differences. Insertions and deletions principally occurred near stretches of homopolymer simple sequence repeats in intergenic regions while most Single Nucleotide Polymorphisms (SNPs) occurred in intergenic regions and introns within the single copy portions of the genome. The polymorphism rate between upland and lowland switchgrass ecotypes was found to be similar to rates reported between chloroplast genomes of indica and japonica subspecies of rice which were believed to have diverged 0.2–0.4 million years ago.

Recent loss of plastid-encoded ndh genes within Erodium (Geraniaceae)

Plastid genomes in the flowering plant family Geraniaceae are known to be highly rearranged based on complete sequences representing the four major genera Erodium, Geranium, Monsonia, and Pelargonium. In this paper we report on the genome sequence of a second species of Erodium, E. carvifolium, representing the second major clade (clade II) in the phylogeny of this genus. Comparison of this genome sequence to the previously published sequence of E. texanum from clade I demonstrates that the plastid genomes of these two species encode the same number of proteins but differ greatly in their relative degree of rearrangement; 14 kb of additional sequence in E. texanum contains complex repeats associated with rearrangement endpoints, whereas the plastid genome of E. carvifolium is streamlined at 116 kb and displays no unique alterations in gene order. Furthermore, these species from both major clades of Erodium contain intact NADH dehydrogenase (ndh) genes, but the 11 ndh genes are represented as pseudogenes in a small clade of 13 species. It is unclear whether plastid-encoded ndh genes have been lost entirely or functionally transferred to the nucleus. This is the third report of the absence of functional ndh genes, and the current study describes the most recent loss of these genes among photosynthetic seed plants and the second such loss among angiosperms. The other ndh losses from Pinaceae/Gnetales and Orchidaceae are much more ancient. Comparative biochemistry between Erodium species with and without plastid-encoded ndh genes may elucidate changes in photosynthetic function and the role of the Ndh complex.

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.

The chloroplast genome of the green alga Schizomeris leibleinii (Chlorophyceae) provides evidence for bidirectional DNA replication from a single origin in the chaetophorales

In the Chlorophyceae, the chloroplast genome is extraordinarily fluid in architecture and displays unique features relative to other groups of green algae. For the Chaetophorales, 1 of the 5 major lineages of the Chlorophyceae, it has been shown that the distinctive architecture of the 223,902-bp genome of Stigeoclonium helveticum is consistent with bidirectional DNA replication from a single origin. Here, we report the 182,759-bp chloroplast genome sequence of Schizomeris leibleinii, a member of the earliest diverging lineage of the Chaetophorales. Like its Stigeoclonium homolog, the Schizomeris genome lacks a large inverted repeat encoding the rRNA operon and displays a striking bias in coding regions that is associated with a bias in base composition along each strand. Our results support the notion that these two chaetophoralean genomes replicate bidirectionally from a putative origin located in the vicinity of the small subunit ribosomal RNA gene. Their shared structural characteristics were most probably inherited from the common ancestor of all chaetophoralean algae. Short dispersed repeats account for most of the 41-kb size variation between the Schizomeris and Stigeoclonium genomes, and there is no indication that homologous recombination between these repeated elements led to the observed gene rearrangements. A comparison of the extent of variation sustained by the Stigeoclonium and Schizomeris chloroplast DNAs (cpDNAs) with that observed for the cpDNAs of the chlamydomonadalean Chlamydomonas and Volvox suggests that gene rearrangements as well as changes in the abundance of intergenic and intron sequences occurred at a slower pace in the Chaetophorales than in the Chlamydomonadales.

A global picture of tRNA genes in plant genomes

Although transfer RNA (tRNA) has a fundamental role in cell life, little is known about tRNA gene organization and expression on a genome-wide scale in eukaryotes, particularly plants. Here, we analyse the content and distribution of tRNA genes in five flowering plants and one green alga. The tRNA gene content is homogenous in plants, and is mostly correlated with genome size. The number of tRNA pseudogenes and organellar-like tRNA genes present in nuclear genomes varies greatly from one plant species to another. These pseudogenes or organellar-like genes appear to be generated or inserted randomly during evolution. Interestingly, we identified a new family of tRNA-related short interspersed nuclear elements (SINEs) in the Populus trichocarpa nuclear genome. In higher plants, intron-containing tRNA genes are rare, and correspond to genes coding for tRNA(Tyr) and tRNA(Mete) . By contrast, in green algae, more than half of the tRNA genes contain an intron. This suggests divergent means of intron acquisition and the splicing process between green algae and land plants. Numerous tRNAs are co-transcribed in Chlamydomonas, but they are mostly transcribed as a single unit in flowering plants. The only exceptions are tRNA(Gly) -snoRNA and tRNA(Mete) -snoRNA cotranscripts in dicots and monocots, respectively. The internal or external motifs required for efficient transcription of tRNA genes by RNA polymerase III are well conserved among angiosperms. A brief analysis of the mitochondrial and plastidial tRNA gene populations is also provided.

Comparative chloroplast genomes of pinaceae: insights into the mechanism of diversified genomic organizations

Pinaceae, the largest family of conifers, has diversified organizations of chloroplast genomes (cpDNAs) with the two typical inverted repeats (IRs) highly reduced. To unravel the mechanism of this genomic diversification, we examined the cpDNA organizations from 53 species of the ten Pinaceous genera, including those of Larix decidua (122,474 bp), Picea morrisonicola (124,168 bp), and Pseudotsuga wilsoniana (122,513 bp), which were firstly elucidated. The results uncovered four distinct cpDNA forms (A−C and P) that are due to rearrangements of two ∼20 and ∼21 kb specific fragments. The C form was documented for the first time and the A form might be the most ancestral one. In addition, only the individuals of Ps. macrocarpa and Ps. wilsoniana were detected to have isomeric cpDNA forms. Three types (types 1−3) of Pinaceae-specific repeats situated nearby the rearranged fragments were found to be syntenic. We hypothesize that type 1 (949 ± 343 bp) and type 3 (608 ± 73 bp) repeats are substrates for homologous recombination (HR), whereas type 2 repeats are likely inactive for HR because of their relatively short sizes (151 ± 30 bp). Conversions among the four distinct forms may be achieved by HR and mediated by type 1 or 3 repeats, thus resulting in increased diversity of cpDNA organizations. We propose that in the Pinaceae cpDNAs, the reduced IRs have lost HR activity, then decreasing the diversity of cpDNA organizations, but the specific repeats that the evolution endowed Pinaceae complement the reduced IRs and increase the diversity of cpDNA organizations.

Characterization of the complete chloroplast genome of Hevea brasiliensis reveals genome rearrangement, RNA editing sites and phylogenetic relationships

Rubber tree (Hevea brasiliensis) is an economical plant and widely grown for natural rubber production. However, genomic research of rubber tree has lagged behind other species in the Euphorbiaceae family. We report the complete chloroplast genome sequence of rubber tree as being 161,191 bp in length including a pair of inverted repeats of 26,810 bp separated by a small single copy region of 18,362 bp and a large single copy region of 89,209 bp. The chloroplast genome contains 112 unique genes, 16 of which are duplicated in the inverted repeat. Of the 112 unique genes, 78 are predicted protein-coding genes, 4 are ribosomal RNA genes and 30 are tRNA genes. Relative to other plant chloroplast genomes, we observed a unique rearrangement in the rubber tree chloroplast genome: a 30-kb inversion between the trnE(UUC)-trnS(GCU) and the trnT(GGU)-trnR(UCU). A comparison between the rubber tree chloroplast genes and cDNA sequences revealed 51 RNA editing sites in which most (48 sites) were located in 26 protein coding genes and the other 3 sites were in introns. Phylogenetic analysis based on chloroplast genes demonstrated a close relationship between Hevea and Manihot in Euphorbiaceae and provided a strong support for a monophyletic group of the eurosid I.

Plant plastid engineering

Genetic material in plants is distributed into nucleus, plastids and mitochondria. Plastid has a central role of carrying out photosynthesis in plant cells. Plastid transformation is becoming more popular and an alternative to nuclear gene transformation because of various advantages like high protein levels, the feasibility of expressing multiple proteins from polycistronic mRNAs, and gene containment through the lack of pollen transmission. Recently, much progress in plastid engineering has been made. In addition to model plant tobacco, many transplastomic crop plants have been generated which possess higher resistance to biotic and abiotic stresses and molecular pharming. In this mini review, we will discuss the features of the plastid DNA and advantages of plastid transformation. We will also present some examples of transplastomic plants developed so far through plastid engineering, and the various applications of plastid transformation.

The complete chloroplast genome sequence of date palm (Phoenix dactylifera L.)

BACKGROUND

Date palm (Phoenix dactylifera L.), a member of Arecaceae family, is one of the three major economically important woody palms--the two other palms being oil palm and coconut tree--and its fruit is a staple food among Middle East and North African nations, as well as many other tropical and subtropical regions. Here we report a complete sequence of the data palm chloroplast (cp) genome based on pyrosequencing.

METHODOLOGY/PRINCIPAL FINDINGS

After extracting 369,022 cp sequencing reads from our whole-genome-shotgun data, we put together an assembly and validated it with intensive PCR-based verification, coupled with PCR product sequencing. The date palm cp genome is 158,462 bp in length and has a typical quadripartite structure of the large (LSC, 86,198 bp) and small single-copy (SSC, 17,712 bp) regions separated by a pair of inverted repeats (IRs, 27,276 bp). Similar to what has been found among most angiosperms, the date palm cp genome harbors 112 unique genes and 19 duplicated fragments in the IR regions. The junctions between LSC/IRs and SSC/IRs show different features of sequence expansion in evolution. We identified 78 SNPs as major intravarietal polymorphisms within the population of a specific cp genome, most of which were located in genes with vital functions. Based on RNA-sequencing data, we also found 18 polycistronic transcription units and three highly expression-biased genes--atpF, trnA-UGC, and rrn23.

CONCLUSIONS

Unlike most monocots, date palm has a typical cp genome similar to that of tobacco--with little rearrangement and gene loss or gain. High-throughput sequencing technology facilitates the identification of intravarietal variations in cp genomes among different cultivars. Moreover, transcriptomic analysis of cp genes provides clues for uncovering regulatory mechanisms of transcription and translation in chloroplasts.