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

The complete chloroplast genome sequence of the medicinal plant Morinda officinalis (Rubiaceae), an endemic to China

The complete chloroplast genome of Morinda officinalis, an endangered and important Chinese medicine with great economic value, has been sequenced in this article. The genome size is 153 398 bp in length, with 38.05% GC content. A pair of inverted repeats (IRs, 51 834 bp) are separated by a large single copy region (LSC, 83 996 bp) and a small single copy region (SSC, 17 566 bp). The chloroplast genome contains 103 unique genes, 80 protein-coding genes, 19 tRNA genes, and 4 rRNA genes. In these genes, 8 genes contained 1 intron, and 2 genes comprised of 2 introns.

Characterization of the complete chloroplast genome of Juniperus cedrus (Cupressaceae)

The complete chloroplast genome of Juniperus cedrus is a circular molecule of 127 126 bp in length with 115 single copy genes and two duplicated genes (trnI-CAU and trnQ-UUG). The genome contains 82 protein-coding genes, four ribosomal RNA genes and 33 transfer RNA genes. In these genes, eight genes (rpl16, rpl2, rpoC1, petD, petB, ndhB, ndhA and atpF) harbor a single intron and two genes (rps12 and ycf3) harbor two introns. Like other sequenced chloroplast genomes of conifers, this genome does not contain canonical inverted repeats (IRs), and the overall GC content of J. cedrus chloroplast DNA is 35%. The phylogenetic analysis revealed that J. cedrus is more closely related to J. scopulorum and J. bermudiana.

The Chloroplast Genome of Elaeagnus macrophylla and trnH Duplication Event in Elaeagnaceae

Elaeagnaceae, which harbor nitrogen-fixing actinomycetes, is a plant family of the Rosales and sister to Rhamnaceae, Barbeyaceae and Dirachmaceae. The results of previous molecular studies have not strongly supported the families of Elaeagnaceae, Rhamnaceae, Barbeyaceae and Dirachmaceae. However, chloroplast genome studies provide valuable phylogenetic information; therefore, we determined the chloroplast genome of Elaeaganus macrophylla and compared it to that of Rosales such as IR junction and infA gene. The chloroplast genome of Elaeagnus macrophylla is 152,224 bp in length and the infA gene of E. macrophylla was psuedogenation. Phylogenetic analyses based on 79 genes in 30 species revealed that Elaeagnus was closely related to Morus. Comparison of the IR junction in six other rosids revealed that the trnH gene contained the LSC region, whereas E. macrophylla contained a trnH gene duplication in the IR region. Comparison of the LSC/IRb (JLB) and the IRa/LSC (JLA) regions of Elaeagnaceae (Elaeagnus and Shephedia) and Rhamnaceae (Rhamnus) showed that trnH gene duplication only occurred in the Elaeagnaceae. The complete chloroplast genome of Elaeagnus macrophylla provides unique characteristics in rosids. The infA gene has been lost or transferred to the nucleus in rosids, while E. macrophylla lost the infA gene. Evaluation of the chloroplast genome of Elaeagnus revealed trnH gene duplication for the first time in rosids. The availability of Elaeagnus cp genomes provides valuable information describing the relationship of Elaeagnaceae, Barbeyaceae and Dirachmaceae, IR junction that will be valuable to future systematics studies.

Development of chloroplast genomic resources for Cynara

In this study, new chloroplast (cp) resources were developed for the genus Cynara, using whole cp genomes from 20 genotypes, by means of high-throughput sequencing technologies. Our target species included seven globe artichokes, two cultivated cardoons, eight wild artichokes, and three other wild Cynara species (C. baetica, C. cornigera and C. syriaca). One complete cp genome was isolated using short reads from a whole-genome sequencing project, while the others were obtained by means of long-range PCR, for which primer pairs are provided here. A de novo assembly strategy combined with a reference-based assembly allowed us to reconstruct each cp genome. Comparative analyses among the newly sequenced genotypes and two additional Cynara cp genomes ('Brindisino' artichoke and C. humilis) retrieved from public databases revealed 126 parsimony informative characters and 258 singletons in Cynara, for a total of 384 variable characters. Thirty-nine SSR loci and 34 other INDEL events were detected. After data analysis, 37 primer pairs for SSR amplification were designed, and these molecular markers were subsequently validated in our Cynara genotypes. Phylogenetic analysis based on all cp variable characters provided the best resolution when compared to what was observed using only parsimony informative characters, or only short 'variable' cp regions. The evaluation of the molecular resources obtained from this study led us to support the 'super-barcode' theory and consider the total cp sequence of Cynara as a reliable and valuable molecular marker for exploring species diversity and examining variation below the species level.

Plastome Evolution in Hemiparasitic Mistletoes

Santalales is an order of plants consisting almost entirely of parasites. Some, such as Osyris, are facultative root parasites whereas others, such as Viscum, are obligate stem parasitic mistletoes. Here, we report the complete plastome sequences of one species of Osyris and three species of Viscum, and we investigate the evolutionary aspects of structural changes and changes in gene content in relation to parasitism. Compared with typical angiosperms plastomes, the four Santalales plastomes are all reduced in size (10–22% compared with Vitis), and they have experienced rearrangements, mostly but not exclusively in the border areas of the inverted repeats. Additionally, a number of protein-coding genes (matK, infA, ccsA, rpl33, and all 11 ndh genes) as well as two transfer RNA genes (trnG-UCC and trnV-UAC) have been pseudogenized or completely lost. Most of the remaining plastid genes have a significantly changed selection pattern compared with other dicots, and the relaxed selection of photosynthesis genes is noteworthy. Although gene loss obviously reduces plastome size, intergenic regions were also shortened. As plastome modifications are generally most prominent in Viscum, they are most likely correlated with the increased nutritional dependence on the host compared with Osyris.

The complete chloroplast genome sequence of Euonymus japonicus (Celastraceae)

The complete chloroplast (cp) genome sequence of the Euonymus japonicus, the first sequenced of the genus Euonymus, was reported in this study. The total length was 157 637 bp, containing a pair of 26 678 bp inverted repeat region (IR), which were separated by small single copy (SSC) region and large single copy (LSC) region of 18 340 bp and 85 941 bp, respectively. This genome contains 107 unique genes, including 74 coding genes, four rRNA genes, and 29 tRNA genes. Seventeen genes contain intron of E. japonicus, of which three genes (clpP, ycf3, and rps12) include two introns. The maximum likelihood (ML) phylogenetic analysis revealed that E. japonicus was closely related to Manihot and Populus.

Short-range inversions: rethinking organelle genome stability: template switching events during DNA replication destabilize organelle genomes

In the organelles of plants and mammals, recent evidence suggests that genomic instability stems in large part from template switching events taking place during DNA replication. Although more than one mechanism may be responsible for this, some similarities exist between the different proposed models. These can be separated into two main categories, depending on whether they involve a single-strand-switching or a reciprocal-strand-switching event. Single-strand-switching events lead to intermediates containing Y junctions, whereas reciprocal-strand-switching creates Holliday junctions. Common features in all the described models include replication stress, fork stalling and the presence of inverted repeats, but no single element appears to be required in all cases. We review the field, and examine the ideas that several mechanisms may take place in any given genome, and that the presence of palindromes or inverted repeats in certain regions may favor specific rearrangements.

The Highly Reduced Plastome of Mycoheterotrophic Sciaphila (Triuridaceae) Is Colinear with Its Green Relatives and Is under Strong Purifying Selection

The enigmatic monocot family Triuridaceae provides a potentially useful model system for studying the effects of an ancient loss of photosynthesis on the plant plastid genome, as all of its members are mycoheterotrophic and achlorophyllous. However, few studies have placed the family in a comparative context, and its phylogenetic placement is only partly resolved. It was also unclear whether any taxa in this family have retained a plastid genome. Here, we used genome survey sequencing to retrieve plastid genome data for Sciaphila densiflora (Triuridaceae) and ten autotrophic relatives in the orders Dioscoreales and Pandanales. We recovered a highly reduced plastome for Sciaphila that is nearly colinear with Carludovica palmata, a photosynthetic relative that belongs to its sister group in Pandanales, Cyclanthaceae-Pandanaceae. This phylogenetic placement is well supported and robust to a broad range of analytical assumptions in maximum-likelihood inference, and is congruent with recent findings based on nuclear and mitochondrial evidence. The 28 genes retained in the S. densiflora plastid genome are involved in translation and other nonphotosynthetic functions, and we demonstrate that nearly all of the 18 protein-coding genes are under strong purifying selection. Our study confirms the utility of whole plastid genome data in phylogenetic studies of highly modified heterotrophic plants, even when they have substantially elevated rates of substitution.

The chloroplast genome of the hexaploid Spartina maritima (Poaceae, Chloridoideae): Comparative analyses and molecular dating

The history of many plant lineages is complicated by reticulate evolution with cases of hybridization often followed by genome duplication (allopolyploidy). In such a context, the inference of phylogenetic relationships and biogeographic scenarios based on molecular data is easier using haploid markers like chloroplast genome sequences. Hybridization and polyploidization occurred recurrently in the genus Spartina (Poaceae, Chloridoideae), as illustrated by the recent formation of the invasive allododecaploid S. anglica during the 19th century in Europe. Until now, only a few plastid markers were available to explore the history of this genus and their low variability limited the resolution of species relationships. We sequenced the complete chloroplast genome (plastome) of S. maritima, the native European parent of S. anglica, and compared it to the plastomes of other Poaceae. Our analysis revealed the presence of fast-evolving regions of potential taxonomic, phylogeographic and phylogenetic utility at various levels within the Poaceae family. Using secondary calibrations, we show that the tetraploid and hexaploid lineages of Spartina diverged 6-10 my ago, and that the two parents of the invasive allopolyploid S. anglica separated 2-4 my ago via long distance dispersal of the ancestor of S. maritima over the Atlantic Ocean. Finally, we discuss the meaning of divergence times between chloroplast genomes in the context of reticulate evolution.

The Highly Reduced Plastome of Mycoheterotrophic Sciaphila (Triuridaceae) Is Colinear with Its Green Relatives and Is under Strong Purifying Selection

The enigmatic monocot family Triuridaceae provides a potentially useful model system for studying the effects of an ancient loss of photosynthesis on the plant plastid genome, as all of its members are mycoheterotrophic and achlorophyllous. However, few studies have placed the family in a comparative context, and its phylogenetic placement is only partly resolved. It was also unclear whether any taxa in this family have retained a plastid genome. Here, we used genome survey sequencing to retrieve plastid genome data for Sciaphila densiflora (Triuridaceae) and ten autotrophic relatives in the orders Dioscoreales and Pandanales. We recovered a highly reduced plastome for Sciaphila that is nearly colinear with Carludovica palmata, a photosynthetic relative that belongs to its sister group in Pandanales, Cyclanthaceae-Pandanaceae. This phylogenetic placement is well supported and robust to a broad range of analytical assumptions in maximum-likelihood inference, and is congruent with recent findings based on nuclear and mitochondrial evidence. The 28 genes retained in the S. densiflora plastid genome are involved in translation and other nonphotosynthetic functions, and we demonstrate that nearly all of the 18 protein-coding genes are under strong purifying selection. Our study confirms the utility of whole plastid genome data in phylogenetic studies of highly modified heterotrophic plants, even when they have substantially elevated rates of substitution.

YCF1: A Green TIC?

A pivotal step in the transformation of an endosymbiotic cyanobacterium to a plastid some 1.5 billion years ago was the evolution of a protein import apparatus, the TOC/TIC machinery, in the common ancestor of Archaeplastida. Recently, a putative new TIC member was identified in Arabidopsis thaliana: TIC214. This finding is remarkable for a number of reasons: (1) TIC214 is encoded by ycf1, so it would be the first plastid-encoded protein of this apparatus; (2) ycf1 is unique to the green lineage (Chloroplastida) but entirely lacking in glaucophytes (Glaucophyta) and the red lineage (Rhodophyta) of the Archaeplastida; (3) ycf1 has been shown to be one of the few indispensable plastid genes (aside from the ribosomal machinery), yet it is missing in the grasses; and (4) 30 years of previous TOC/TIC research missed it. These observations prompted us to survey the evolution of ycf1. We found that ycf1 is not only lacking in grasses and some parasitic plants, but also for instance in cranberry (Ericaceae). The encoded YCF proteins are highly variable, both in sequence length and in the predicted number of N-terminal transmembrane domains. The evolution of the TOC/TIC machinery in the green lineage experienced specific modifications, but our analysis does not support YCF1 to be a general green TIC. It remains to be explained how the apparent complete loss of YCF1 can be tolerated by some embryophytes and whether what is observed for YCF1 function in a member of the Brassicaceae is also true for, e.g., algal and noncanonical YCF1 homologs.

YCF1: A Green TIC: Response to the de Vries et al. Commentary

This response to a recent Commentary article by de Vries et al. highlights critical errors in the annotation and identification of Ycf1 homologs in the sequenced chloroplast genomes. Contrary to what is reported by de Vries et al., the majority of chloroplast genomes sequenced to date appear to have retained a typical Ycf1 sequence (i.e., including the N-terminal 6TM domain and a variable hydrophilic C-terminal domain) as my group previously reported. Our evidence continues to support the model that Ycf1 forms an essential component of a "green TIC" that is largely conserved among the Chlorophyta and land plants. Since the establishment of this green TIC with Tic20 as the core component, some cases of loss of Ycf1 during the evolution of the green lineages might be regarded as modifications or alterations of the complex. Here, I discuss our working model that the presence of an alternative "nonphotosynthetic-type" or "ancestral-type" TIC might explain other (or specific) cases of the lack of Ycf1, not only in early lineages, including Glaucophyta and Rhodophyta, but also in the grasses.

Complete chloroplast genome sequence of Adenophora remotiflora (Campanulaceae)

The complete chloroplast genome of Adenophora remotiflora was determined in this study. The total genome size was 171 724 bp in length, containing a pair of inverted repeats of 27 437 bp, which were separated by large single copy and small single copy of 105 555 bp and 11 295 bp, respectively. The overall GC contents of the plastid genome were 38.8%. One hundred and twelve unique genes were annotated, including 78 protein coding genes, 30 tRNA genes and 4 rRNA genes. In these genes, 15 genes contained one intron and 2 genes had two introns. A phylogenetic tree showed that Adenophora remotiflora was closely related to Hanabusaya asiatica.

Evidence for horizontal transfer of mitochondrial DNA to the plastid genome in a bamboo genus

In flowering plants, three genomes (nuclear, mitochondrial, and plastid) coexist and intracellular horizontal transfer of DNA is prevalent, especially from the plastid to the mitochondrion genome. However, the plastid genomes are generally conserved in evolution and have long been considered immune to foreign DNA. Recently, the opposite direction of DNA transfer from the mitochondrial to the plastid genome has been reported in two eudicot lineages. Here we sequenced 6 plastid genomes of bamboos, three of which are neotropical woody species and three are herbaceous ones. Several unusual features were found, including the duplication of trnT-GGU and loss of one copy of rps19 due to contraction of inverted repeats (IRs). The most intriguing was the ~2.7 kb insertion in the plastid IR regions in the three herbaceous bamboos. Furthermore, the insertion was documented to be horizontally transferred from the mitochondrial to the plastid genome. Our study provided evidence of the mitochondrial-to-plastid DNA transfer in the monocots, demonstrating again that this rare event does occur in other angiosperm lineages. However, the mechanism underlying the transfer remains obscure, and more studies in other plants may elucidate it in the future.

The First Complete Chloroplast Genome Sequences in Actinidiaceae: Genome Structure and Comparative Analysis

Actinidia chinensis is an important economic plant belonging to the basal lineage of the asterids. Availability of a complete Actinidia chloroplast genome sequence is crucial to understanding phylogenetic relationships among major lineages of angiosperms and facilitates kiwifruit genetic improvement. We report here the complete nucleotide sequences of the chloroplast genomes for Actinidia chinensis and A. chinensis var deliciosa obtained through de novo assembly of Illumina paired-end reads produced by total DNA sequencing. The total genome size ranges from 155,446 to 157,557 bp, with an inverted repeat (IR) of 24,013 to 24,391 bp, a large single copy region (LSC) of 87,984 to 88,337 bp and a small single copy region (SSC) of 20,332 to 20,336 bp. The genome encodes 113 different genes, including 79 unique protein-coding genes, 30 tRNA genes and 4 ribosomal RNA genes, with 16 duplicated in the inverted repeats, and a tRNA gene (trnfM-CAU) duplicated once in the LSC region. Comparisons of IR boundaries among four asterid species showed that IR/LSC borders were extended into the 5’ portion of the psbA gene and IR contraction occurred in Actinidia. The clap gene has been lost from the chloroplast genome in Actinidia, and may have been transferred to the nucleus during chloroplast evolution. Twenty-seven polymorphic simple sequence repeat (SSR) loci were identified in the Actinidia chloroplast genome. Maximum parsimony analyses of a 72-gene, 16 taxa angiosperm dataset strongly support the placement of Actinidiaceae in Ericales within the basal asterids.

The chloroplast genome sequence of an important medicinal plant Dioscorea nipponica

Dioscorea nipponica is an important medicinal plant belonging to Dioscoreaceae, a family which is vital for the evolution of monocotyledon. In the present study, the nucleotide sequence of the D. nipponica chloroplast genome was determined. It was an AT-rich (63.3%) chloroplast genome with 152,946 bp in length, containing a pair of 23,113 bp inverted repeats, which were separated by a large and a small single copy region of 83,557 bp and 23,064 bp in length, respectively. It encodes 120 unique genes, including 89 protein-coding genes, 27 tRNA genes and 4 rRNA genes. The predicted gene-coding regions covered 58.7% of the genome sequences. Ten genes contained one intron, while two genes had two introns. Phylogenetic analyses showed the present chloroplast genome can be used as a potential supper barcode to distinguish D. nipponica from its closely related species. Furthermore, the chloroplast genome provides a molecular base for the next investigation on this important medicinal species.

Transcriptome and genome size analysis of the Venus flytrap

The insectivorous Venus flytrap (Dionaea muscipula) is renowned from Darwin’s studies of plant carnivory and the origins of species. To provide tools to analyze the evolution and functional genomics of D. muscipula, we sequenced a normalized cDNA library synthesized from mRNA isolated from D. muscipula flowers and traps. Using the Oases transcriptome assembler 79,165,657 quality trimmed reads were assembled into 80,806 cDNA contigs, with an average length of 679 bp and an N50 length of 1,051 bp. A total of 17,047 unique proteins were identified, and assigned to Gene Ontology (GO) and classified into functional categories. A total of 15,547 full-length cDNA sequences were identified, from which open reading frames were detected in 10,941. Comparative GO analyses revealed that D. muscipula is highly represented in molecular functions related to catalytic, antioxidant, and electron carrier activities. Also, using a single copy sequence PCR-based method, we estimated that the genome size of D. muscipula is approx. 3 Gb. Our genome size estimate and transcriptome analyses will contribute to future research on this fascinating, monotypic species and its heterotrophic adaptations.

Mutation rates in plastid genomes: they are lower than you might think

Within plastid-bearing species, the mutation rate of the plastid genome is often assumed to be greater than that of the mitochondrial genome. This assumption is based on early, pioneering studies of land plant molecular evolution, which uncovered higher rates of synonymous substitution in plastid versus mitochondrial DNAs. However, much of the plastid-containing eukaryotic diversity falls outside of land plants, and the patterns of plastid DNA evolution for embryophytes do not necessarily reflect those of other groups. Recent analyses of plastid and mitochondrial substitution rates in diverse lineages have uncovered very different trends than those recorded for land plants. Here, I explore these new data and argue that for many protists the plastid mutation rate is lower than that of the mitochondrion, including groups with primary or secondary plastids as well as nonphotosynthetic algae. These findings have far-reaching implications for how we view plastid genomes and how their sequences are used for evolutionary analyses, and might ultimately reflect a general tendency toward more efficient DNA repair mechanisms in plastids than in mitochondria.

Ancient nuclear plastid DNA in the yew family (taxaceae)

Plastid-to-nucleus DNA transfer provides a rich genetic resource to the complexity of plant nuclear genome architecture. To date, the evolutionary route of nuclear plastid DNA (nupt) remain unknown in conifers. We have sequenced the complete plastomes of two yews, Amentotaxus formosana and Taxus mairei (Taxaceae of coniferales). Our comparative genomic analyses recovered an evolutionary scenario for plastomic reorganization from ancestral to extant plastomes in the three sampled Taxaceae genera, Amentotaxus, Cephalotaxus, and Taxus. Specific primers were designed to amplify nonsyntenic regions between ancestral and extant plastomes, and 12.6 kb of nupts were identified based on phylogenetic analyses. These nupts have significantly accumulated GC-to-AT mutations, reflecting a nuclear mutational environment shaped by spontaneous deamination of 5-methylcytosin. The ancestral initial codon of rps8 is retained in the T. nupts, but its corresponding extant codon is mutated and requires C-to-U RNA-editing. These findings suggest that nupts can help recover scenarios of the nucleotide mutation process. We show that the Taxaceae nupts we retrieved may have been retained because the Cretaceous and they carry information of both ancestral genomic organization and nucleotide composition, which offer clues for understanding the plastome evolution in conifers.

Mitochondrial and plastid genome architecture: Reoccurring themes, but significant differences at the extremes

Mitochondrial and plastid genomes show a wide array of architectures, varying immensely in size, structure, and content. Some organelle DNAs have even developed elaborate eccentricities, such as scrambled coding regions, nonstandard genetic codes, and convoluted modes of posttranscriptional modification and editing. Here, we compare and contrast the breadth of genomic complexity between mitochondrial and plastid chromosomes. Both organelle genomes have independently evolved many of the same features and taken on similar genomic embellishments, often within the same species or lineage. This trend is most likely because the nuclear-encoded proteins mediating these processes eventually leak from one organelle into the other, leading to a high likelihood of processes appearing in both compartments in parallel. However, the complexity and intensity of genomic embellishments are consistently more pronounced for mitochondria than for plastids, even when they are found in both compartments. We explore the evolutionary forces responsible for these patterns and argue that organelle DNA repair processes, mutation rates, and population genetic landscapes are all important factors leading to the observed convergence and divergence in organelle genome architecture.

ycf1, the most promising plastid DNA barcode of land plants

A DNA barcode is a DNA fragment used to identify species. For land plants, DNA fragments of plastid genome could be the primary consideration. Unfortunately, most of the plastid candidate barcodes lack species-level resolution. The identification of DNA barcodes of high resolution at species level is critical to the success of DNA barcoding in plants. We searched the available plastid genomes for the most variable regions and tested the best candidates using both a large number of tree species and seven well-sampled plant groups. Two regions of the plastid gene ycf1, ycf1a and ycf1b, were the most variable loci that were better than existing plastid candidate barcodes and can serve as a barcode of land plants. Primers were designed for the amplification of these regions, and the PCR success of these primers ranged from 82.80% to 98.17%. Of 420 tree species, 357 species could be distinguished using ycf1b, which was slightly better than the combination of matK and rbcL. For the well-sampled representative plant groups, ycf1b generally performed better than any of the matK, rbcL and trnH-psbA. We concluded that ycf1a or ycf1b is the most variable plastid genome region and can serve as a core barcode of land plants.

Chloroplast RNA polymerases: Role in chloroplast biogenesis

Plastid genes are transcribed by two types of RNA polymerase in angiosperms: the bacterial type plastid-encoded RNA polymerase (PEP) and one (RPOTp in monocots) or two (RPOTp and RPOTmp in dicots) nuclear-encoded RNA polymerase(s) (NEP). PEP is a bacterial-type multisubunit enzyme composed of core subunits (coded for by the plastid rpoA, rpoB, rpoC1 and rpoC2 genes) and additional protein factors (sigma factors and polymerase associated protein, PAPs) encoded in the nuclear genome. Sigma factors are required by PEP for promoter recognition. Six different sigma factors are used by PEP in Arabidopsis plastids. NEP activity is represented by phage-type RNA polymerases. Only one NEP subunit has been identified, which bears the catalytic activity. NEP and PEP use different promoters. Many plastid genes have both PEP and NEP promoters. PEP dominates in the transcription of photosynthesis genes. Intriguingly, rpoB belongs to the few genes transcribed exclusively by NEP. Both NEP and PEP are active in non-green plastids and in chloroplasts at all stages of development. The transcriptional activity of NEP and PEP is affected by endogenous and exogenous factors. This article is part of a Special Issue entitled: Chloroplast Biogenesis.

Are differences in genomic data sets due to true biological variants or errors in genome assembly: an example from two chloroplast genomes

DNA sequencing has been revolutionized by the development of high-throughput sequencing technologies. Plummeting costs and the massive throughput capacities of second and third generation sequencing platforms have transformed many fields of biological research. Concurrently, new data processing pipelines made rapid de novo genome assemblies possible. However, high quality data are critically important for all investigations in the genomic era. We used chloroplast genomes of one Oryza species (O. australiensis) to compare differences in sequence quality: one genome (GU592209) was obtained through Illumina sequencing and reference-guided assembly and the other genome (KJ830774) was obtained via target enrichment libraries and shotgun sequencing. Based on the whole genome alignment, GU592209 was more similar to the reference genome (O. sativa: AY522330) with 99.2% sequence identity (SI value) compared with the 98.8% SI values in the KJ830774 genome; whereas the opposite result was obtained when the SI values in coding and noncoding regions of GU592209 and KJ830774 were compared. Additionally, the junctions of two single copies and repeat copies in the chloroplast genome exhibited differences. Phylogenetic analyses were conducted using these sequences, and the different data sets yielded dissimilar topologies: phylogenetic replacements of the two individuals were remarkably different based on whole genome sequencing or SNP data and insertions and deletions (indels) data. Thus, we concluded that the genomic composition of GU592209 was heterogeneous in coding and non-coding regions. These findings should impel biologists to carefully consider the quality of sequencing and assembly when working with next-generation data.

Complete plastome sequence of Thalictrum coreanum (Ranunculaceae) and transfer of the rpl32 gene to the nucleus in the ancestor of the subfamily Thalictroideae

BACKGROUND

Plastids originated from cyanobacteria and the majority of the ancestral genes were lost or functionally transferred to the nucleus after endosymbiosis. Comparative genomic investigations have shown that gene transfer from plastids to the nucleus is an ongoing evolutionary process but molecular evidence for recent functional gene transfers among seed plants have only been documented for the four genes accD, infA, rpl22, and rpl32.

RESULTS

The complete plastid genome of Thalictrum coreanum, the first from the subfamily Thalictroideae (Ranunculaceae), was sequenced and revealed the losses of two genes, infA and rpl32. The functional transfer of these two genes to the nucleus in Thalictrum was verified by examination of nuclear transcriptomes. A survey of the phylogenetic distribution of the rpl32 loss was performed using 17 species of Thalictrum and representatives of related genera in the subfamily Thalictroideae. The plastid-encoded rpl32 gene is likely nonfunctional in members of the subfamily Thalictroideae (Aquilegia, Enemion, Isopyrum, Leptopyrum, Paraquilegia, and Semiaquilegia) including 17 Thalictrum species due to the presence of indels that disrupt the reading frame. A nuclear-encoded rpl32 with high sequence identity was identified in both Thalictrum and Aquilegia. The phylogenetic distribution of this gene loss/transfer and the high level of sequence similarity in transit peptides suggest a single transfer of the plastid-encoded rpl32 to the nucleus in the ancestor of the subfamily Thalictroideae approximately 20-32 Mya.

CONCLUSIONS

The genome sequence of Thalictrum coreanum provides valuable information for improving the understanding of the evolution of plastid genomes within Ranunculaceae and across angiosperms. Thalictrum is unusual among the three sequenced Ranunculaceae plastid genomes in the loss of two genes infA and rpl32, which have been functionally transferred to the nucleus. In the case of rpl32 this represents the third documented independent transfer from the plastid to the nucleus with the other two transfers occurring in the unrelated angiosperm families Rhizophoraceae and Salicaceae. Furthermore, the transfer of rpl32 provides additional molecular evidence for the monophyly of the subfamily Thalictroideae.

Complete chloroplast genome of Macadamia integrifolia confirms the position of the Gondwanan early-diverging eudicot family Proteaceae

BACKGROUND

Sequence data from the chloroplast genome have played a central role in elucidating the evolutionary history of flowering plants, Angiospermae. In the past decade, the number of complete chloroplast genomes has burgeoned, leading to well-supported angiosperm phylogenies. However, some relationships, particulary among early-diverging lineages, remain unresolved. The diverse Southern Hemisphere plant family Proteaceae arose on the ancient supercontinent Gondwana early in angiosperm history and is a model group for adaptive radiation in response to changing climatic conditions. Genomic resources for the family are limited, and until now it is one of the few early-diverging 'basal eudicot' lineages not represented in chloroplast phylogenomic analyses.

RESULTS

The chloroplast genome of the Australian nut crop tree Macadamia integrifolia was assembled de novo from Illumina paired-end sequence reads. Three contigs, corresponding to a collapsed inverted repeat, a large and a small single copy region were identified, and used for genome reconstruction. The complete genome is 159,714 bp in length and was assembled at deep coverage (3.29 million reads; ~2000 x). Phylogenetic analyses based on 83-gene and inverted repeat region alignments, the largest sequence-rich datasets to include the basal eudicot family Proteaceae, provide strong support for a Proteales clade that includes Macadamia, Platanus and Nelumbo. Genome structure and content followed the ancestral angiosperm pattern and were highly conserved in the Proteales, whilst size differences were largely explained by the relative contraction of the single copy regions and expansion of the inverted repeats in Macadamia.

CONCLUSIONS

The Macadamia chloroplast genome presented here is the first in the Proteaceae, and confirms the placement of this family with the morphologically divergent Plantanaceae (plane tree family) and Nelumbonaceae (sacred lotus family) in the basal eudicot order Proteales. It provides a high-quality reference genome for future evolutionary studies and will be of benefit for taxon-rich phylogenomic analyses aimed at resolving relationships among early-diverging angiosperms, and more broadly across the plant tree of life.

RHON1 mediates a Rho-like activity for transcription termination in plastids of Arabidopsis thaliana

Although transcription termination is essential to generate functional RNAs, its underlying molecular mechanisms are still poorly understood in plastids of vascular plants. Here, we show that the RNA binding protein RHON1 participates in transcriptional termination of rbcL (encoding large subunit of ribulose-1,5-bisphosphate carboxylase/oxygenase) in Arabidopsis thaliana. Inactivation of RHON1 leads to enhanced rbcL read-through transcription and to aberrant accD (encoding β-subunit of the acetyl-CoA carboxylase) transcriptional initiation, which may result from inefficient transcription termination of rbcL. RHON1 can bind to the mRNA as well as to single-stranded DNA of rbcL, displays an RNA-dependent ATPase activity, and terminates transcription of rbcL in vitro. These results suggest that RHON1 terminates rbcL transcription using an ATP-driven mechanism similar to that of Rho of Escherichia coli. This RHON1-dependent transcription termination occurs in Arabidopsis but not in rice (Oryza sativa) and appears to reflect a fundamental difference between plastomes of dicotyledonous and monocotyledonous plants. Our results point to the importance and significance of plastid transcription termination and provide insights into its machinery in an evolutionary context.

Evolutionary origin of highly repetitive plastid genomes within the clover genus (Trifolium)

Background

Some clover species, particularly Trifolium subterraneum, have previously been reported to have highly unusual plastomes, relative to closely related legumes, enlarged with many duplications, gene losses and the presence of DNA unique to Trifolium, which may represent horizontal transfer. In order to pinpoint the evolutionary origin of this phenomenon within the genus Trifolium, we sequenced and assembled the plastomes of eight additional Trifolium species widely sampled from across the genus.

Results

The Trifolium plastomes fell into two groups: those of Trifolium boissieri, T. strictum and T. glanduliferum (representing subgenus Chronosemium and subg. Trifolium section Paramesus) were tractable, assembled readily and were not unusual in the general context of Fabeae plastomes. The other Trifolium species (“core Trifolium”) proved refractory to assembly mainly because of numerous short duplications. These species form a single clade, which we call the “refractory clade” (comprising subg, Trifolium sections Lupinaster, Trifolium, Trichocephalum, Vesicastrum and Trifoliastrum). The characteristics of the refractory clade are the presence of numerous short duplications and 7-15% longer genomes than the tractable species. Molecular dating estimates that the origin of the most recent common ancestor (MRCA) of the refractory clade is approximately 13.1 million years ago (MYA). This is considerably younger than the estimated MRCA ages of Trifolium (c. 18.6 MYA) and Trifolium subg. Trifolium (16.1 MYA).

Conclusions

We conclude that the unusual repetitive plastome type previously characterized in Trifolium subterraneum had a single origin within Trifolium and is characteristic of most (but not all) species of subgenus Trifolium. It appears that an ancestral plastome within Trifolium underwent an evolutionary change resulting in plastomes that either actively promoted, were permissive to, or were unable to control, duplications within the genome. The precise mechanism of this important change in the mode and tempo of plastome evolution deserves further investigation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0228-6) contains supplementary material, which is available to authorized users.

ChloroSSRdb: a repository of perfect and imperfect chloroplastic simple sequence repeats (cpSSRs) of green plants

Simple sequence repeats (SSRs) are regions in DNA sequence that contain repeating motifs of length 1-6 nucleotides. These repeats are ubiquitously present and are found in both coding and non-coding regions of genome. A total of 534 complete chloroplast genome sequences (as on 18 September 2014) of Viridiplantae are available at NCBI organelle genome resource. It provides opportunity to mine these genomes for the detection of SSRs and store them in the form of a database. In an attempt to properly manage and retrieve chloroplastic SSRs, we designed ChloroSSRdb which is a relational database developed using SQL server 2008 and accessed through ASP.NET. It provides information of all the three types (perfect, imperfect and compound) of SSRs. At present, ChloroSSRdb contains 124 430 mined SSRs, with majority lying in non-coding region. Out of these, PCR primers were designed for 118 249 SSRs. Tetranucleotide repeats (47 079) were found to be the most frequent repeat type, whereas hexanucleotide repeats (6414) being the least abundant. Additionally, in each species statistical analyses were performed to calculate relative frequency, correlation coefficient and chi-square statistics of perfect and imperfect SSRs. In accordance with the growing interest in SSR studies, ChloroSSRdb will prove to be a useful resource in developing genetic markers, phylogenetic analysis, genetic mapping, etc. Moreover, it will serve as a ready reference for mined SSRs in available chloroplast genomes of green plants. Database URL: www.compubio.in/chlorossrdb/

Complete chloroplast genome sequence of poisonous and medicinal plant Datura stramonium: organizations and implications for genetic engineering

Datura stramonium is a widely used poisonous plant with great medicinal and economic value. Its chloroplast (cp) genome is 155,871 bp in length with a typical quadripartite structure of the large (LSC, 86,302 bp) and small (SSC, 18,367 bp) single-copy regions, separated by a pair of inverted repeats (IRs, 25,601 bp). The genome contains 113 unique genes, including 80 protein-coding genes, 29 tRNAs and four rRNAs. A total of 11 forward, 9 palindromic and 13 tandem repeats were detected in the D. stramonium cp genome. Most simple sequence repeats (SSR) are AT-rich and are less abundant in coding regions than in non-coding regions. Both SSRs and GC content were unevenly distributed in the entire cp genome. All preferred synonymous codons were found to use A/T ending codons. The difference in GC contents of entire genomes and of the three-codon positions suggests that the D. stramonium cp genome might possess different genomic organization, in part due to different mutational pressures. The five most divergent coding regions and four non-coding regions (trnH-psbA, rps4-trnS, ndhD-ccsA, and ndhI-ndhG) were identified using whole plastome alignment, which can be used to develop molecular markers for phylogenetics and barcoding studies within the Solanaceae. Phylogenetic analysis based on 68 protein-coding genes supported Datura as a sister to Solanum. This study provides valuable information for phylogenetic and cp genetic engineering studies of this poisonous and medicinal plant.

chlB requirement for chlorophyll biosynthesis under short photoperiod in Marchantia polymorpha L

Chlorophylls (Chls) play pivotal roles in energy absorption and transduction and also in charge separation in reaction centers in all photosynthetic organisms. In Chl biosynthesis steps, only a step for the enzymatic reduction of protochlorophyllide (Pchlide) to chlorophyllide (Chlide) is mediated by both nuclear- and chloroplast-encoded genes in land plants. Many plants encode the genes for light-dependent Pchlide reductase (LPOR) and light-independent Pchlide reductase (DPOR) in the nucleus and chloroplast genome, respectively. During the diversification of land plants, the reduction step of Pchlide to Chlide has become solely dependent on LPOR, and the genes for DPOR have been lost from chloroplast genome. It remains unclear why DPOR persists in some land plants, how they were eliminated from chloroplast genomes during the diversification of land plants, and under what environmental conditions DPOR was required. We demonstrate that DPOR is functional in liverwort (Marchantia polymorpha L.) and plays an important role in Chl biosynthesis. Having established a plastid transformation system in liverwort, we disrupted chlB, which encodes a subunit of DPOR in the M. polymorpha chloroplast genome. Morphological and Chl content analysis of a chlB mutant grown under different photoperiods revealed that DPOR is particularly required for Chl biosynthesis under short-day conditions. Our findings suggest that an environmental condition in the form of photoperiod is an important factor that determines the loss or retention of chloroplast-encoded genes mediating Pchlide reduction to Chlide.

Whole plastome sequences from five ginger species facilitate marker development and define limits to barcode methodology

Plants from the Zingiberaceae family are a key source of spices and herbal medicines. Species identification within this group is critical in the search for known and possibly novel bioactive compounds. To facilitate precise characterization of this group, we have sequenced chloroplast genomes from species representing five major groups within Zingiberaceae. Generally, the structure of these genomes is similar to the basal angiosperm excepting an expansion of 3 kb associated with the inverted repeat A region. Portions of this expansion appear to be shared across the entire Zingiberales order, which includes gingers and bananas. We used whole plastome alignment information to develop DNA barcodes that would maximize the ability to differentiate species within the Zingiberaceae. Our computation pipeline identified regions of high variability that were flanked by highly conserved regions used for primer design. This approach yielded hitherto unexploited regions of variability. These theoretically optimal barcodes were tested on a range of species throughout the family and were found to amplify and differentiate genera and, in some cases, species. Still, though these barcodes were specifically optimized for the Zingiberaceae, our data support the emerging consensus that whole plastome sequences are needed for robust species identification and phylogenetics within this family.

Complete chloroplast genome sequence of Hanabusaya asiatica (Campanulaceae), an endemic genus to Korea

The complete chloroplast genome of Hanabusaya asiatica, an endemic genus to Korea, was determined in this study. The total genome size was 167,287 bp in length, containing a pair of inverted repeats (IRs) of 26,877 bp, which were separated by large single copy (LSC) and small single copy (SSC) of 104,955 and 8578, respectively. The overall GC contents of the plastid genome were 38.8%. One hundred and twelve unique genes were annotated, including 78 protein coding genes, 30 tRNA genes and 4 rRNA genes. In these genes, 17 genes contained 1 intron, and 2 genes comprised of 2 introns.

Investigating the path of plastid genome degradation in an early-transitional clade of heterotrophic orchids, and implications for heterotrophic angiosperms

Parasitic organisms exemplify morphological and genomic reduction. Some heterotrophic, parasitic plants harbor drastically reduced and degraded plastid genomes resulting from relaxed selective pressure on photosynthetic function. However, few studies have addressed the initial stages of plastome degradation in groups containing both photosynthetic and nonphotosynthetic species. Corallorhiza is a genus of leafless, heterotrophic orchids that contains both green, photosynthetic species and nongreen, putatively nonphotosynthetic species, and represents an ideal system in which to assess the beginning of the transition to a "minimal plastome." Complete plastomes were generated for nine taxa of Corallorhiza using Illumina paired-end sequencing of genomic DNA to assess the degree of degradation among taxa, and for comparison with a general model of degradation among angiosperms. Quantification of total chlorophyll suggests that nongreen Corallorhiza still produce chlorophyll, but at 10-fold lower concentrations than green congeners. Complete plastomes and partial nuclear rDNA cistrons yielded a fully resolved tree for Corallorhiza, with at least two independent losses of photosynthesis, evidenced by gene deletions and pseudogenes in Co. striata and nongreen Co. maculata. All Corallorhiza show some evidence of degradation in genes of the NAD(P)H dehydrogenase complex. Among genes with open reading frames, photosynthesis-related genes displayed evidence of neutral evolution in nongreen Corallorhiza, whereas genes of the ATP synthase complex displayed some evidence of positive selection in these same groups, though for reasons unknown. Corallorhiza spans the early stages of a general model of plastome degradation and has added critical insight for understanding the process of plastome evolution in heterotrophic angiosperms.

Complete plastome sequences from Glycine syndetika and six additional perennial wild relatives of soybean

Organelle sequences have a long history of utility in phylogenetic analyses. Chloroplast sequences when combined with nuclear data can help resolve relationships among flowering plant genera, and within genera incongruence can point to reticulate evolution. Plastome sequences are becoming plentiful because they are increasingly easier to obtain. Complete plastome sequences allow us to detect rare rearrangements and test the tempo of sequence evolution. Chloroplast sequences are generally considered a nuisance to be kept to a minimum in bacterial artificial chromosome libraries. Here, we sequenced two bacterial artificial chromosomes per species to generate complete plastome sequences from seven species. The plastome sequences from Glycine syndetika and six other perennial Glycine species are similar in arrangement and gene content to the previously published soybean plastome. Repetitive sequences were detected in high frequencies as in soybean, but further analysis showed that repeat sequence numbers are inflated. Previous chloroplast-based phylogenetic trees for perennial Glycine were incongruent with nuclear gene-based phylogenetic trees. We tested whether the hypothesis of introgression was supported by the complete plastomes. Alignment of complete plastome sequences and Bayesian analysis allowed us to date putative hybridization events supporting the hypothesis of introgression and chloroplast "capture."

Evolutionary and biotechnology implications of plastid genome variation in the inverted-repeat-lacking clade of legumes

Land plant plastid genomes (plastomes) provide a tractable model for evolutionary study in that they are relatively compact and gene dense. Among the groups that display an appropriate level of variation for structural features, the inverted-repeat-lacking clade (IRLC) of papilionoid legumes presents the potential to advance general understanding of the mechanisms of genomic evolution. Here, are presented six complete plastome sequences from economically important species of the IRLC, a lineage previously represented by only five completed plastomes. A number of characters are compared across the IRLC including gene retention and divergence, synteny, repeat structure and functional gene transfer to the nucleus. The loss of clpP intron 2 was identified in one newly sequenced member of IRLC, Glycyrrhiza glabra. Using deeply sequenced nuclear transcriptomes from two species helped clarify the nature of the functional transfer of accD to the nucleus in Trifolium, which likely occurred in the lineage leading to subgenus Trifolium. Legumes are second only to cereal crops in agricultural importance based on area harvested and total production. Genetic improvement via plastid transformation of IRLC crop species is an appealing proposition. Comparative analyses of intergenic spacer regions emphasize the need for complete genome sequences for developing transformation vectors for plastid genetic engineering of legume crops.

The mitochondrial and chloroplast genomes of the haptophyte Chrysochromulina tobin contain unique repeat structures and gene profiles

BACKGROUND

Haptophytes are widely and abundantly distributed in both marine and freshwater ecosystems. Few genomic analyses of representatives within this taxon have been reported, despite their early evolutionary origins and their prominent role in global carbon fixation.

RESULTS

The complete mitochondrial and chloroplast genome sequences of the haptophyte Chrysochromulina tobin (Prymnesiales) provide insight into the architecture and gene content of haptophyte organellar genomes. The mitochondrial genome (~34 kb) encodes 21 protein coding genes and contains a complex, 9 kb tandem repeat region. Similar to other haptophytes and rhodophytes, but not cryptophytes or stramenopiles, the mitochondrial genome has lost the nad7, nad9 and nad11 genes. The ~105 kb chloroplast genome encodes 112 protein coding genes, including ycf39 which has strong structural homology to NADP-binding nitrate transcriptional regulators; a divergent 'CheY-like' two-component response regulator (ycf55) and Tic/Toc (ycf60 and ycf80) membrane transporters. Notably, a zinc finger domain has been identified in the rpl36 ribosomal protein gene of all chloroplasts sequenced to date with the exception of haptophytes and cryptophytes--algae that have gained (via lateral gene transfer) an alternative rpl36 lacking the zinc finger motif. The two C. tobin chloroplast ribosomal RNA operon spacer regions differ in tRNA content. Additionally, each ribosomal operon contains multiple single nucleotide polymorphisms (SNPs)--a pattern observed in rhodophytes and cryptophytes, but few stramenopiles. Analysis of small (<200 bp) chloroplast encoded tandem and inverted repeats in C. tobin and 78 other algal chloroplast genomes show that repeat type, size and location are correlated with gene identity and taxonomic clade.

CONCLUSION

The Chrysochromulina tobin organellar genomes provide new insight into organellar function and evolution. These are the first organellar genomes to be determined for the prymnesiales, a taxon that is present in both oceanic and freshwater systems and represents major primary photosynthetic producers and contributors to global ecosystem stability.

The dynamic history of plastid genomes in the Campanulaceae sensu lato is unique among angiosperms

Why have some plants lost the organizational stability in plastid genomes (plastomes) that evolved in their algal ancestors? During the endosymbiotic transformation of a cyanobacterium into the eukaryotic plastid, most cyanobacterial genes were transferred to the nucleus or otherwise lost from the plastome, and the resulting plastome architecture in land plants confers organizational stability, as evidenced by the conserved gene order among bryophytes and lycophytes, whereas ferns, gymnosperms, and angiosperms share a single, 30-kb inversion. Although some additional gene losses have occurred, gene additions to angiosperm plastomes were previously unknown. Plastomes in the Campanulaceae sensu lato have incorporated dozens of large ORFs (putative protein-coding genes). These insertions apparently caused many of the 125+ large inversions now known in this small eudicot clade. This phylogenetically restricted phenomenon is not biogeographically localized, which indicates that these ORFs came from the nucleus or (less likely) a cryptic endosymbiont.

Transplastomic integration of a cyanobacterial bicarbonate transporter into tobacco chloroplasts

Improving global yields of agricultural crops is a complex challenge with evidence indicating benefits in productivity are achieved by enhancing photosynthetic carbon assimilation. Towards improving rates of CO2 capture within leaf chloroplasts, this study shows the versatility of plastome transformation for expressing the Synechococcus PCC7002 BicA bicarbonate transporter within tobacco plastids. Fractionation of chloroplast membranes from transplastomic tob(BicA) lines showed that ~75% of the BicA localized to the thylakoid membranes and ~25% to the chloroplast envelope. BicA levels were highest in young emerging tob(BicA) leaves (0.12 μmol m(-2), ≈7mg m(-2)) accounting for ~0.1% (w/w) of the leaf protein. In these leaves, the molar amount of BicA was 16-fold lower than the abundant thylakoid photosystem II D1 protein (~1.9 μmol m(-2)) which was comparable to the 9:1 molar ratio of D1:BicA measured in air-grown Synechococcus PCC7002 cells. The BicA produced had no discernible effect on chloroplast ultrastructure, photosynthetic CO2-assimilation rates, carbon isotope discrimination, or growth of the tob(BicA) plants, implying that the bicarbonate transporter had little or no activity. These findings demonstrate the utility of plastome transformation for targeting bicarbonate transporter proteins into the chloroplast membranes without impeding growth or plastid ultrastructure. This study establishes the span of experimental measurements required to verify heterologous bicarbonate transporter function and location in chloroplasts and underscores the need for more detailed understanding of BicA structure and function to identify solutions for enabling its activation and operation in leaf chloroplasts.

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.

Highly rearranged and size-variable chloroplast genomes in conifers II clade (cupressophytes): evolution towards shorter intergenic spacers

Although conifers are of immense ecological and economic value, bioengineering of their chloroplasts remains undeveloped. Understanding the chloroplast genomic organization of conifers can facilitate their bioengineering. Members of the conifer II clade (or cupressophytes) are highly diverse in both morphologic features and chloroplast genomic organization. We compared six cupressophyte chloroplast genomes (cpDNAs) that represent four of the five cupressophyte families, including three genomes that are first reported here (Agathis dammara, Calocedrus formosana and Nageia nagi). The six cupressophyte cpDNAs have lost a pair of large inverted repeats (IRs) and vary greatly in size, organization and tRNA copies. We demonstrate that cupressophyte cpDNAs have evolved towards reduced size, largely due to shrunken intergenic spacers. In cupressophytes, cpDNA rearrangements are capable of extending intergenic spacers, and synonymous mutations are negatively associated with the size and frequency of rearrangements. The variable cpDNA sizes of cupressophytes may have been shaped by mutational burden and genomic rearrangements. On the basis of cpDNA organization, our analyses revealed that in gymnosperms, cpDNA rearrangements are phylogenetically informative, which supports the 'gnepines' clade. In addition, removal of a specific IR influences the minimal rearrangements required for the gnepines and cupressophyte clades, whereby Pinaceae favours the removal of IRB but cupressophytes exclusion of IRA. This result strongly suggests that different IR copies have been lost from conifers I and II. Our data help understand the complexity and evolution of cupressophyte cpDNAs.

A plastid without a genome: evidence from the nonphotosynthetic green algal genus Polytomella

Polytomella spp. are free-living, nonphotosynthetic green algae closely related to the model organism Chlamydomonas reinhardtii. Although colorless, Polytomella spp. have a plastid, but it is still unknown whether they harbor a plastid genome. We took a next generation sequencing approach, along with transcriptome sequencing, to search for a plastid genome and an associated gene expression system in Polytomella spp. Illumina sequencing of total DNA from four Polytomella spp. did not produce any recognizable plastid-derived reads but did generate a large number of mitochondrial DNA sequences. Transcriptomic analysis of Polytomella parva uncovered hundreds of putative nuclear-encoded, plastid-targeted proteins, which support the presence of plastid-based metabolic functions, similar to those observed in the plastids of other nonphotosynthetic algae. Conspicuously absent, however, were any plastid-targeted proteins involved in the expression, replication, or repair of plastid DNA. Based on these findings and earlier findings, we argue that the Polytomella genus represents the first well-supported example, to our knowledge, of a primary plastid-bearing lineage without a plastid genome.

Analyses of charophyte chloroplast genomes help characterize the ancestral chloroplast genome of land plants

Despite the significance of the relationships between embryophytes and their charophyte algal ancestors in deciphering the origin and evolutionary success of land plants, few chloroplast genomes of the charophyte algae have been reconstructed to date. Here, we present new data for three chloroplast genomes of the freshwater charophytes Klebsormidium flaccidum (Klebsormidiophyceae), Mesotaenium endlicherianum (Zygnematophyceae), and Roya anglica (Zygnematophyceae). The chloroplast genome of Klebsormidium has a quadripartite organization with exceptionally large inverted repeat (IR) regions and, uniquely among streptophytes, has lost the rrn5 and rrn4.5 genes from the ribosomal RNA (rRNA) gene cluster operon. The chloroplast genome of Roya differs from other zygnematophycean chloroplasts, including the newly sequenced Mesotaenium, by having a quadripartite structure that is typical of other streptophytes. On the basis of the improbability of the novel gain of IR regions, we infer that the quadripartite structure has likely been lost independently in at least three zygnematophycean lineages, although the absence of the usual rRNA operonic synteny in the IR regions of Roya may indicate their de novo origin. Significantly, all zygnematophycean chloroplast genomes have undergone substantial genomic rearrangement, which may be the result of ancient retroelement activity evidenced by the presence of integrase-like and reverse transcriptase-like elements in the Roya chloroplast genome. Our results corroborate the close phylogenetic relationship between Zygnematophyceae and land plants and identify 89 protein-coding genes and 22 introns present in the chloroplast genome at the time of the evolutionary transition of plants to land, all of which can be found in the chloroplast genomes of extant charophytes.

A pangenomic analysis of the Nannochloropsis organellar genomes reveals novel genetic variations in key metabolic genes

BACKGROUND

Microalgae in the genus Nannochloropsis are photosynthetic marine Eustigmatophytes of significant interest to the bioenergy and aquaculture sectors due to their ability to efficiently accumulate biomass and lipids for utilization in renewable transportation fuels, aquaculture feed, and other useful bioproducts. To better understand the genetic complement that drives the metabolic processes of these organisms, we present the assembly and comparative pangenomic analysis of the chloroplast and mitochondrial genomes from Nannochloropsis salina CCMP1776.

RESULTS

The chloroplast and mitochondrial genomes of N. salina are 98.4% and 97% identical to their counterparts in Nannochloropsis gaditana. Comparison of the Nannochloropsis pangenome to other algae within and outside of the same phyla revealed regions of significant genetic divergence in key genes that encode proteins needed for regulation of branched chain amino synthesis (acetohydroxyacid synthase), carbon fixation (RuBisCO activase), energy conservation (ATP synthase), protein synthesis and homeostasis (Clp protease, ribosome).

CONCLUSIONS

Many organellar gene modifications in Nannochloropsis are unique and deviate from conserved orthologs found across the tree of life. Implementation of secondary and tertiary structure prediction was crucial to functionally characterize many proteins and therefore should be implemented in automated annotation pipelines. The exceptional similarity of the N. salina and N. gaditana organellar genomes suggests that N. gaditana be reclassified as a strain of N. salina.

Combined analysis of the chloroplast genome and transcriptome of the Antarctic vascular plant Deschampsia antarctica Desv

BACKGROUND

Antarctic hairgrass (Deschampsia antarctica Desv.) is the only natural grass species in the maritime Antarctic. It has been researched as an important ecological marker and as an extremophile plant for studies on stress tolerance. Despite its importance, little genomic information is available for D. antarctica. Here, we report the complete chloroplast genome, transcriptome profiles of the coding/noncoding genes, and the posttranscriptional processing by RNA editing in the chloroplast system.

RESULTS

The complete chloroplast genome of D. antarctica is 135,362 bp in length with a typical quadripartite structure, including the large (LSC: 79,881 bp) and small (SSC: 12,519 bp) single-copy regions, separated by a pair of identical inverted repeats (IR: 21,481 bp). It contains 114 unique genes, including 81 unique protein-coding genes, 29 tRNA genes, and 4 rRNA genes. Sequence divergence analysis with other plastomes from the BEP clade of the grass family suggests a sister relationship between D. antarctica, Festuca arundinacea and Lolium perenne of the Poeae tribe, based on the whole plastome. In addition, we conducted high-resolution mapping of the chloroplast-derived transcripts. Thus, we created an expression profile for 81 protein-coding genes and identified ndhC, psbJ, rps19, psaJ, and psbA as the most highly expressed chloroplast genes. Small RNA-seq analysis identified 27 small noncoding RNAs of chloroplast origin that were preferentially located near the 5'- or 3'-ends of genes. We also found >30 RNA-editing sites in the D. antarctica chloroplast genome, with a dominance of C-to-U conversions.

CONCLUSIONS

We assembled and characterized the complete chloroplast genome sequence of D. antarctica and investigated the features of the plastid transcriptome. These data may contribute to a better understanding of the evolution of D. antarctica within the Poaceae family for use in molecular phylogenetic studies and may also help researchers understand the characteristics of the chloroplast transcriptome.

Predominant and substoichiometric isomers of the plastid genome coexist within Juniperus plants and have shifted multiple times during cupressophyte evolution

Most land plant plastomes contain two copies of a large inverted repeat (IR) that promote high-frequency homologous recombination to generate isomeric genomic forms. Among conifer plastomes, this canonical IR is highly reduced in Pinaceae and completely lost from cupressophytes. However, both lineages have acquired short, novel IRs, some of which also exhibit recombinational activity to generate genomic structural diversity. This diversity has been shown to exist between, and occasionally within, cupressophyte species, but it is not known whether multiple genomic forms coexist within individual plants. To examine the recombinational potential of the novel cupressophyte IRs within individuals and between species, we sequenced the plastomes of four closely related species of Juniperus. The four plastomes have identical gene content and genome organization except for a large 36 kb inversion between approximately 250 bp IR containing trnQ-UUG. Southern blotting showed that different isomeric versions of the plastome predominate among individual junipers, whereas polymerase chain reaction and high-throughput read-pair mapping revealed the substoichiometric presence of the alternative isomeric form within each individual plant. Furthermore, our comparative genomic studies demonstrate that the predominant and substoichiometric arrangements of this IR have changed several times in other cupressophytes as well. These results provide compelling evidence for substoichiometric shifting of plastomic forms during cupressophyte evolution and suggest that substoichiometric shifting activity in plastid genomes may be adaptive.

From algae to angiosperms-inferring the phylogeny of green plants (Viridiplantae) from 360 plastid genomes

BACKGROUND

Next-generation sequencing has provided a wealth of plastid genome sequence data from an increasingly diverse set of green plants (Viridiplantae). Although these data have helped resolve the phylogeny of numerous clades (e.g., green algae, angiosperms, and gymnosperms), their utility for inferring relationships across all green plants is uncertain. Viridiplantae originated 700-1500 million years ago and may comprise as many as 500,000 species. This clade represents a major source of photosynthetic carbon and contains an immense diversity of life forms, including some of the smallest and largest eukaryotes. Here we explore the limits and challenges of inferring a comprehensive green plant phylogeny from available complete or nearly complete plastid genome sequence data.

RESULTS

We assembled protein-coding sequence data for 78 genes from 360 diverse green plant taxa with complete or nearly complete plastid genome sequences available from GenBank. Phylogenetic analyses of the plastid data recovered well-supported backbone relationships and strong support for relationships that were not observed in previous analyses of major subclades within Viridiplantae. However, there also is evidence of systematic error in some analyses. In several instances we obtained strongly supported but conflicting topologies from analyses of nucleotides versus amino acid characters, and the considerable variation in GC content among lineages and within single genomes affected the phylogenetic placement of several taxa.

CONCLUSIONS

Analyses of the plastid sequence data recovered a strongly supported framework of relationships for green plants. This framework includes: i) the placement of Zygnematophyceace as sister to land plants (Embryophyta), ii) a clade of extant gymnosperms (Acrogymnospermae) with cycads + Ginkgo sister to remaining extant gymnosperms and with gnetophytes (Gnetophyta) sister to non-Pinaceae conifers (Gnecup trees), and iii) within the monilophyte clade (Monilophyta), Equisetales + Psilotales are sister to Marattiales + leptosporangiate ferns. Our analyses also highlight the challenges of using plastid genome sequences in deep-level phylogenomic analyses, and we provide suggestions for future analyses that will likely incorporate plastid genome sequence data for thousands of species. We particularly emphasize the importance of exploring the effects of different partitioning and character coding strategies.

Possible loss of the chloroplast genome in the parasitic flowering plant Rafflesia lagascae (Rafflesiaceae)

Rafflesia is a genus of holoparasitic plants endemic to Southeast Asia that has lost the ability to undertake photosynthesis. With short-read sequencing technology, we assembled a draft sequence of the mitochondrial genome of Rafflesia lagascae Blanco, a species endemic to the Philippine island of Luzon, with ∼350× sequencing depth coverage. Using multiple approaches, however, we were only able to identify small fragments of plastid sequences at low coverage depth (<2×) and could not recover any substantial portion of a chloroplast genome. The gene fragments we identified included photosynthesis and energy production genes (atp, ndh, pet, psa, psb, rbcL), ribosomal RNA genes (rrn16, rrn23), ribosomal protein genes (rps7, rps11, rps16), transfer RNA genes, as well as matK, accD, ycf2, and multiple nongenic regions from the inverted repeats. None of the identified plastid gene sequences had intact reading frames. Phylogenetic analysis suggests that ∼33% of these remnant plastid genes may have been horizontally transferred from the host plant genus Tetrastigma with the rest having ambiguous phylogenetic positions (<50% bootstrap support), except for psaB that was strongly allied with the plastid homolog in Nicotiana. Our inability to identify substantial plastid genome sequences from R. lagascae using multiple approaches—despite success in identifying and developing a draft assembly of the much larger mitochondrial genome—suggests that the parasitic plant genus Rafflesia may be the first plant group for which there is no recognizable plastid genome, or if present is found in cryptic form at very low levels.