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

Assembly of chloroplast genomes with long- and short-read data: a comparison of approaches using Eucalyptus pauciflora as a test case

BACKGROUND

Chloroplasts are organelles that conduct photosynthesis in plant and algal cells. The information chloroplast genome contained is widely used in agriculture and studies of evolution and ecology. Correctly assembling chloroplast genomes can be challenging because the chloroplast genome contains a pair of long inverted repeats (10-30 kb). Typically, it is simply assumed that the gross structure of the chloroplast genome matches the most commonly observed structure of two single-copy regions separated by a pair of inverted repeats. The advent of long-read sequencing technologies should remove the need to make this assumption by providing sufficient information to completely span the inverted repeat regions. Yet, long-reads tend to have higher error rates than short-reads, and relatively little is known about the best way to combine long- and short-reads to obtain the most accurate chloroplast genome assemblies. Using Eucalyptus pauciflora, the snow gum, as a test case, we evaluated the effect of multiple parameters, such as different coverage of long-(Oxford nanopore) and short-(Illumina) reads, different long-read lengths, different assembly pipelines, with a view to determining the most accurate and efficient approach to chloroplast genome assembly.

RESULTS

Hybrid assemblies combining at least 20x coverage of both long-reads and short-reads generated a single contig spanning the entire chloroplast genome with few or no detectable errors. Short-read-only assemblies generated three contigs (the long single copy, short single copy and inverted repeat regions) of the chloroplast genome. These contigs contained few single-base errors but tended to exclude several bases at the beginning or end of each contig. Long-read-only assemblies tended to create multiple contigs with a much higher single-base error rate. The chloroplast genome of Eucalyptus pauciflora is 159,942 bp, contains 131 genes of known function.

CONCLUSIONS

Our results suggest that very accurate assemblies of chloroplast genomes can be achieved using a combination of at least 20x coverage of long- and short-reads respectively, provided that the long-reads contain at least ~5x coverage of reads longer than the inverted repeat region. We show that further increases in coverage give little or no improvement in accuracy, and that hybrid assemblies are more accurate than long-read-only or short-read-only assemblies.

Complete chloroplast genome sequence of Betula platyphylla: gene organization, RNA editing, and comparative and phylogenetic analyses

BACKGROUND

Betula platyphylla is a common tree species in northern China that has high economic and medicinal value. Our laboratory has been devoted to genome research on B. platyphylla for approximately 10 years. As primary organelle genomes, the complete genome sequences of chloroplasts are important to study the divergence of species, RNA editing and phylogeny. In this study, we sequenced and analyzed the complete chloroplast (cp) genome sequence of B. platyphylla.

RESULTS

The complete cp genome of B. platyphylla was 160,518 bp in length, which included a pair of inverted repeats (IRs) of 26,056 bp that separated a large single copy (LSC) region of 89,397 bp and a small single copy (SSC) region of 19,009 bp. The annotation contained a total of 129 genes, including 84 protein-coding genes, 37 tRNA genes and 8 rRNA genes. There were 3 genes using alternative initiation codons. Comparative genomics showed that the sequence of the Fagales species cp genome was relatively conserved, but there were still some high variation regions that could be used as molecular markers. The IR expansion event of B. platyphylla resulted in larger cp genomes and rps19 pseudogene formation. The simple sequence repeat (SSR) analysis showed that there were 105 SSRs in the cp genome of B. platyphylla. RNA editing sites recognition indicated that at least 80 RNA editing events occurred in the cp genome. Most of the substitutions were C to U, while a small proportion of them were not. In particular, three editing loci on the rRNA were converted to more than two other bases that had never been reported. For synonymous conversion, most of them increased the relative synonymous codon usage (RSCU) value of the codons. The phylogenetic analysis suggested that B. platyphylla had a closer evolutionary relationship with B. pendula than B. nana.

CONCLUSIONS

In this study, we not only obtained and annotated the complete cp genome sequence of B. platyphylla, but we also identified new RNA editing sites and predicted the phylogenetic relationships among Fagales species. These findings will facilitate genomic, genetic engineering and phylogenetic studies of this important species.

Characterization of chloroplast genomes of Alnus rubra and Betula cordifolia, and their use in phylogenetic analyses in Betulaceae

OBJECTIVE

Betulaceae is a relatively small birch family that comprises about 160 deciduous trees and shrubs. Chloroplast (cp) genome sequencing of Alnus rubra and Betula cordifolia was carried out to elucidate their molecular features and phylogenetic relationship among species in Betulaceae family.

METHODS

Chloroplast genome sequencing was carried out using next generation sequencing method. Molecular and genomic features of the two cp genomes were characterized with other cp genomes in Betulaceae. Also, molecular phylogenetic analysis was performed using the whole cp genome sequences.

RESULTS

The average cp genome length was 160,136 bp among the Betulaceae species. Base compositions of the cp genomes were skewed toward a high AT ratio, with an average of 63.4%. We identified 117 different genes 83 with protein coding, 4 with ribosomal RNA, and 30 with tRNA. Eighteen genes contained introns which were conserved among the cp genomes of all Betulaceae. We mined 82 SSRs from the cp genomes of A. rubra, A. cordifolia, and A. nana. The SSRs were variable in motif repeat numbers and presence/absence among the cp genomes.

CONCLUSION

Chloroplast genome-wide sequence comparison from 11 Betulaceae species and one cp genome of evergreen oak revealed that the patterns of sequence variations were congruent with two subfamily classification Betuloideae (Alnus and Betula) and Corylaceae (Corylus, Ostrya, and Carpinus). Subsequent phylogenetic analysis also supports the sub-classifications of these species.

Complete chloroplast genome sequence of Dryopteris fragrans (L.) Schott and the repeat structures against the thermal environment

Dryopteris fragrans (L.) Schott is a fern growing on the surface of hot rocks and lava. It is exposed to sunlight directly and bears local hot environment. We sequenced the complete nucleotide sequence of its chloroplast (cp) genome. The cp genome was 151,978 bp in length, consisting of a large single-copy region (85,332 bp), a small single-copy region (31,947 bp) and a pair of inverted repeats (17,314 bp). The cp genome contained 112 genes and 345 RNA editing sites in protein-coding genes. Simple sequence repeats (SSRs) and long repeat structure pairs (30–55 bp) were identified. The number and percent of repeat structures are extremely high in ferns. Thermal denaturation experiments showed its cp genome to have numerous, dispersed and high GC percent repeat structures, which conferred the strongest thermal stability. This repeat-heavy genome may provide the molecular basis of how D. fragrans cp survives its hot environment.

Regulated chloroplast transcription termination

Transcription termination by the RNA polymerase (RNAP) is a fundamental step of gene expression that involves the release of the nascent transcript and dissociation of the RNAP from the DNA template. However, the functional importance of termination extends beyond the mere definition of the gene borders. Chloroplasts originate from cyanobacteria and possess their own gene expression system. Plastids have a unique hybrid transcription system consisting of two different types of RNAPs of dissimilar phylogenetic origin together with several additional nuclear encoded components. Although the basic components involved in chloroplast transcription have been identified, little attention has been paid to the chloroplast transcription termination. Recent identification and functional characterization of novel factors in regulating transcription termination in Arabidopsis chloroplasts via genetic and biochemical approaches have provided insights into the mechanisms and significance of transcription termination in chloroplast gene expression. This review provides an overview of the current knowledge of the transcription termination in chloroplasts.

Punctuated plastome reduction and host-parasite horizontal gene transfer in the holoparasitic plant genus Aphyllon

Foundational studies of chloroplast genome (plastome) evolution in parasitic plants have focused on broad trends across large clades, particularly among the Orobanchaceae, a species-rich and ecologically diverse family of root parasites. However, the extent to which such patterns and processes of plastome evolution, such as stepwise gene loss following the complete loss of photosynthesis (shift to holoparasitism), are detectable at shallow evolutionary time scale is largely unknown. We used genome skimming to assemble eight chloroplast genomes representing complete taxonomic sampling of Aphyllon sect. Aphyllon, a small clade within the Orobanchaceae that evolved approximately 6 Ma, long after the origin of holoparasitism. We show substantial plastome reduction occurred in the stem lineage, but subsequent change in plastome size, gene content, and structure has been relatively minimal, albeit detectable. This lends additional fine-grained support to existing models of stepwise plastome reduction in holoparasitic plants. Additionally, we report phylogenetic evidence based on an rbcL gene tree and assembled 60+ kb fragments of the Aphyllon epigalium mitochondrial genome indicating host-to-parasite horizontal gene transfers (hpHGT) of several genes originating from the plastome of an ancient Galium host into the mitochondrial genome of a recent common ancestor of A. epigalium Ecologically, this evidence of hpHGT suggests that the host-parasite associations between Galium and A. epigalium have been stable at least since its subspecies diverged hundreds of thousands of years ago.

Insights into population structure of East African sweetpotato cultivars from hybrid assembly of chloroplast genomes

Background: The chloroplast (cp) genome is an important resource for studying plant diversity and phylogeny. Assembly of the cp genomes from next-generation sequencing data is complicated by the presence of two large inverted repeats contained in the cp DNA. Methods: We constructed a complete circular cp genome assembly for the hexaploid sweetpotato using extremely low coverage (<1×) Oxford Nanopore whole-genome sequencing (WGS) data coupled with Illumina sequencing data for polishing. Results: The sweetpotato cp genome of 161,274 bp contains 152 genes, of which there are 96 protein coding genes, 8 rRNA genes and 48 tRNA genes. Using the cp genome assembly as a reference, we constructed complete cp genome assemblies for a further 17 sweetpotato cultivars from East Africa and an I. triloba line using Illumina WGS data. Analysis of the sweetpotato cp genomes demonstrated the presence of two distinct subpopulations in East Africa. Phylogenetic analysis of the cp genomes of the species from the Convolvulaceae Ipomoea section Batatas revealed that the most closely related diploid wild species of the hexaploid sweetpotato is I. trifida. Conclusions: Nanopore long reads are helpful in construction of cp genome assemblies, especially in solving the two long inverted repeats. We are generally able to extract cp sequences from WGS data of sufficiently high coverage for assembly of cp genomes. The cp genomes can be used to investigate the population structure and the phylogenetic relationship for the sweetpotato.

First reported chloroplast genome sequence of Punica granatum (cultivar Helow) from Jabal Al-Akhdar, Oman: phylogenetic comparative assortment with Lagerstroemia

Pomegranate (Punica granatum L.) is one of the oldest known edible fruits. It has grown in popularity and is a profitable fruit crop due to its attractive features including a bright red appearance and its biological activities. Scientific exploration of the genetics and evolution of these beneficial traits has been hampered by limited genomic information. In this study, we sequenced the complete chloroplast (cp) genome of the native P. granatum (cultivar Helow) cultivated in the mountains of Jabal Al-Akhdar, Oman. The results revealed a P. granatum cp genome length of 158,630 bp, characterized by a relatively conserved structure containing 2 inverted repeat regions of 25,466 bp, an 18,686 bp small single copy regions, and an 89,015 bp large single copy region. The 86 protein-coding genes included 37 transfer RNA genes and 8 ribosomal RNA genes. Comparison of the P. granatum whole cp genome with seven Lagerstroemia species revealed an overall high degree of sequence similarity with divergence among intergenic spacers. The location, distribution, and divergence of repeat sequences and shared genes of the Punica and Lagerstroemia species were highly similar. Analyses of nucleotide substitution, insertion/deletions, and highly variable regions in these cp genomes identified potential plastid markers for taxonomic and phylogenetic studies in Myrtales. A phylogenetic study of the cp genomes and 76 shared coding regions generated similar cladograms. The complete cp genome of P. granatum will aid in taxonomical studies of the family Lythraceae.

Commonalities and differences of chloroplast translation in a green alga and land plants

Chloroplast gene expression is a fascinating and highly regulated process, which was mainly studied on specific genes in a few model organisms including the unicellular green alga Chlamydomonas (Chlamydomonas reinhardtii) and the embryophyte (land) plants tobacco (Nicotiana tabacum) and Arabidopsis (Arabidopsis thaliana). However, a direct plastid genome-wide interspecies comparison of chloroplast gene expression that includes translation was missing. We adapted a targeted chloroplast ribosome profiling approach to quantitatively compare RNA abundance and translation output between Chlamydomonas, tobacco and Arabidopsis. The re-analysis of established chloroplast mutants confirmed the capability of the approach by detecting known as well as previously undetected translation defects (including the potential photosystem II assembly-dependent regulation of PsbH). Systematic comparison of the algal and land plant wild-type gene expression showed that, for most genes, the steady-state translation output is highly conserved among the three species, while the levels of transcript accumulation are more distinct. Whereas in Chlamydomonas transcript accumulation and translation output are closely balanced, this correlation is less obvious in embryophytes, indicating more pronounced translational regulation. Altogether, this suggests that green algae and land plants evolved different strategies to achieve conserved levels of protein synthesis.

Plant phylogenomics based on genome-partitioning strategies: Progress and prospects

The rapid expansion of next-generation sequencing (NGS) has generated a powerful array of approaches to address fundamental questions in biology. Several genome-partitioning strategies to sequence selected subsets of the genome have emerged in the fields of phylogenomics and evolutionary genomics. In this review, we summarize the applications, advantages and limitations of four NGS-based genome-partitioning approaches in plant phylogenomics: genome skimming, transcriptome sequencing (RNA-seq), restriction site associated DNA sequencing (RAD-Seq), and targeted capture (Hyb-seq). Of these four genome-partitioning approaches, targeted capture (especially Hyb-seq) shows the greatest promise for plant phylogenetics over the next few years. This review will aid researchers in their selection of appropriate genome-partitioning approaches to address questions of evolutionary scale, where we anticipate continued development and expansion of whole-genome sequencing strategies in the fields of plant phylogenomics and evolutionary biology research.

First chloroplast genomics study of Phoenix dactylifera (var. Naghal and Khanezi): A comparative analysis

Date palm (Phoenix dactylifera L.) is one of the oldest fruit crops in the arid regions of the Middle East. However, little information is available regarding its plastid genomes. In this study, we sequenced the chloroplast (cp) genomes of two economically important but genomically unexplored date palm cultivars of Phoenix dactylifera (var. Naghal and Khanezi). The data assembly and genome annotation revealed a typical quadripartite structure similar to Arecaceae, and the genome sizes of Naghal and Khanezi were 158,210 bp and 158,211 bp, respectively. Structurally, both cp genomes were comprised of four regions: a pair of inverted repeats (27,273 bp for Khanezi and for Naghal 27,272 bp), a large single-copy region (86,090 bp and 86,092 bp) and a small single-copy region (17,575 bp and 17,574 bp). Both genomes had 138 representative genes, whereas 227 and 229 randomly distributed microsatellites were also observed in Khanezi and Naghal, respectively. Phylogenetic analysis based on the whole cp genomes and 68 shared genes showed identical phylogenetic trees of Khanezi and Naghal forming clades with Khalas and Aseel cultivars, respectively. The current study showed detailed comparative cp genome analysis, which could be essential for broader population genetics and molecular studies of these four date palm cultivars.

Understanding plastome evolution in Hemiparasitic Santalales: Complete chloroplast genomes of three species, Dendrotrophe varians, Helixanthera parasitica, and Macrosolen cochinchinensis

Santalales is a large order, with over 2200 species, most of which are root or aerial (stem) hemiparasites. In this study, we report the newly assembled chloroplast genome of Dendrotrophe varians (140,666 bp) in the family Amphorogynaceae and the cp genomes of Helixanthera parasitica (124,881 bp) and Macrosolen cochinchinensis (122,986 bp), both in the family Loranthaceae. We compared the cp genomes of 11 Santalales including eight currently available cp genomes. Santalales cp genomes are slightly or not reduced in size (119-147 kb), similar to other hemiparasitic species, when compared with typical angiosperm cp genomes (120-170 kb). In a phylogeny examining gene content, the NADH dehydrogenase gene group is the only one among eight functional gene groups that lost complete functionally in all examined Santalales. This supports the idea that the functional loss of ndh genes is the initial stage in the evolution of the plastome of parasitic plants, but the loss has occurred independently multiple times in angiosperms, while they are not found in some parasites. This suggests that the functional loss of ndh genes is not essential for the transition from autotroph to parasite. We additionally examined the correlation between gene content and type of parasitism (obligate/facultative and stem/root parasites) of all hemiparasitic species in which cp genomes have been reported to date. Correlation was not found in any types of parasitism.

Next generation sequencing reveals packaging of host RNAs by brome mosaic virus

Although RNA viruses evolved the mechanisms of specific encapsidation, miss-packaging of cellular RNAs has been reported in such RNA virus systems as flock house virus or cucumber necrosis virus. To find out if brome mosaic virus (BMV), a tripartite RNA virus, can package cellular RNAs, BMV was propagated in barley and in Nicotiana benthamiana hosts, purified by cesium chloride (CsCl) gradient ultracentrifugation followed by nuclease treatment to remove any contaminating cellular (host) RNAs. The extracted virion RNA was then sequenced by using next-generation sequencing (NGS RNA-Seq) with the Illumina protocol. Bioinformatic analysis revealed the content of host RNAs ranging from 0.07% for BMV extracted from barley to 0.10% for the virus extracted from N. benthamiana. The viruses from two sources appeared to co-encapsidate different patterns of host-RNAs, including ribosomal RNAs (rRNAs), messenger RNAs (mRNAs) but also mitochondrial and plastid RNAs and, interestingly, transposable elements, both transposons and retrotransposons. Our data reveal that BMV virions can carry host RNAs, having a potential to mediate horizontal gene transfer (HGT) in plants.

Structure-Function Analysis of Chloroplast Proteins via Random Mutagenesis Using Error-Prone PCR

Site-directed mutagenesis of chloroplast genes was developed three decades ago and has greatly advanced the field of photosynthesis research. Here, we describe a new approach for generating random chloroplast gene mutants that combines error-prone polymerase chain reaction of a gene of interest with chloroplast complementation of the knockout Chlamydomonas reinhardtii mutant. As a proof of concept, we targeted a 300-bp sequence of the petD gene that encodes subunit IV of the thylakoid membrane-bound cytochrome b6f complex. By sequencing chloroplast transformants, we revealed 149 mutations in the 300-bp target petD sequence that resulted in 92 amino acid substitutions in the 100-residue target subunit IV sequence. Our results show that this method is suited to the study of highly hydrophobic, multisubunit, and chloroplast-encoded proteins containing cofactors such as hemes, iron-sulfur clusters, and chlorophyll pigments. Moreover, we show that mutant screening and sequencing can be used to study photosynthetic mechanisms or to probe the mutational robustness of chloroplast-encoded proteins, and we propose that this method is a valuable tool for the directed evolution of enzymes in the chloroplast.

The complete plastome sequence of Rubus takesimensis endemic to Ulleung Island, Korea: Insights into molecular evolution of anagenetically derived species in Rubus (Rosaceae)

Previous phylogenetic studies have suggested that Rubus takesimensis (Rosaceae), which is endemic to Ulleung Island, Korea, is closely related to R. crataegifolius, which is broadly distributed across East Asia. A recent phylogeographic study also suggested the possible polyphyletic origins of R. takesimensis from multiple source populations of its continental progenitor R. crataegifolius in China, Japan, Korea, and the Russian Far East. However, even though the progenitor-derivative relationship between R. crataegifolius and R. takesimensis has been established, little is known about the chloroplast genome (i.e., plastome) evolution of anagenetically derived species on oceanic islands and their continental progenitor species. In the present study, we characterized the complete plastome of R. takesimensis and compared it to those of R. crataegifolius and four other Rubus species. The R. takesimensis plastome was 155,760 base pairs (bp) long, a total of 46 bp longer than the plastome of R. crataegifolius (28 from LSC and 18 from SSC). No structural or content rearrangements were found between the species pairs. Four highly variable intergenic regions (rpl32/trnL, rps4/trnT, trnT/trnL, and psbZ/trnG) were identified between R. takesimensis and R. crataegifolius. Compared to the plastomes of other congeneric species (R. corchorifolius, R. fockeanus, and R. niveus), six highly variable intergenic regions (ndhC/psaC, rps16/trnQ, trnK/rps16, trnL/trnF, trnM/atpE, and trnQ/psbK) were also identified. A total of 116 simple sequence repeats (SSRs), including 48 mononucleotide, 64 dinucleotide, and four trinucleotide repeat motifs were characterized in R. takesimensis. The plastome resources generated by the present study will help to elucidate plastome evolution within the genus and to resolve phylogenetic relationships within highly complex and reticulated lineages. Phylogenetic analysis supported both the monophyly of Rubus and the sister relationship between R. crataegifolius and R. takesimensis.

The chloroplast genome sequence of bittersweet (Solanum dulcamara): Plastid genome structure evolution in Solanaceae

Bittersweet (Solanum dulcamara) is a native Old World member of the nightshade family. This European diploid species can be found from marshlands to high mountainous regions and it is a common weed that serves as an alternative host and source of resistance genes against plant pathogens such as late blight (Phytophthora infestans). We sequenced the complete chloroplast genome of bittersweet, which is 155,580 bp in length and it is characterized by a typical quadripartite structure composed of a large (85,901 bp) and small (18,449 bp) single-copy region interspersed by two identical inverted repeats (25,615 bp). It consists of 112 unique genes from which 81 are protein-coding, 27 tRNA and four rRNA genes. All bittersweet plastid genes including non-functional ones and even intergenic spacer regions are transcribed in primary plastid transcripts covering 95.22% of the genome. These are later substantially edited in a post-transcriptional phase to activate gene functions. By comparing the bittersweet plastid genome with all available Solanaceae sequences we found that gene content and synteny are highly conserved across the family. During genome comparison we have identified several annotation errors, which we have corrected in a manual curation process then we have identified the major plastid genome structural changes in Solanaceae. Interpreted in a phylogenetic context they seem to provide additional support for larger clades. The plastid genome sequence of bittersweet could help to benchmark Solanaceae plastid genome annotations and could be used as a reference for further studies. Such reliable annotations are important for gene diversity calculations, synteny map constructions and assigning partitions for phylogenetic analysis with de novo sequenced plastomes of Solanaceae.

Chloroplast genome analyses and genomic resource development for epilithic sister genera Oresitrophe and Mukdenia (Saxifragaceae), using genome skimming data

Background

Epilithic sister genera Oresitrophe and Mukdenia (Saxifragaceae) have an epilithic habitat (rocky slopes) and a parapatric distribution in East Asia, which makes them an ideal model for a more comprehensive understanding of the demographic and divergence history and the influence of climate changes in East Asia. However, the genetic background and resources for these two genera are scarce.

Results

The complete chloroplast (cp) genomes of two Oresitrophe rupifraga and one Mukdenia rossii individuals were reconstructed and comparative analyses were conducted to examine the evolutionary pattern of chloroplast genomes in Saxifragaceae. The cp genomes ranged from 156,738 bp to 156,960 bp in length and had a typical quadripartite structure with a conserved genome arrangement. Comparative analysis revealed the intron of rpl2 has been lost in Heuchera parviflora, Tiarella polyphylla, M. rossii and O. rupifraga but presents in the reference genome of Penthorum chinense. Seven cp hotspot regions (trnH-psbA, trnR-atpA, atpI-rps2, rps2-rpoC2, petN-psbM, rps4-trnT and rpl33-rps18) were identified between Oresitrophe and Mukdenia, while four hotspots (trnQ-psbK, trnR-atpA, trnS-psbZ and rpl33-rps18) were identified within Oresitrophe. In addition, 24 polymorphic cpSSR loci were found between Oresitrophe and Mukdenia. Most importantly, we successfully developed 126 intergeneric polymorphic gSSR markers between Oresitrophe and Mukdenia, as well as 452 intrageneric ones within Oresitrophe. Twelve randomly selected intergeneric gSSRs have shown that these two genera exhibit a significant genetic structure.

Conclusions

In this study, we conducted genome skimming for Oresitrophe rupifraga and Mukdenia rossii. Using these data, we were able to not only assemble their complete chloroplast genomes, but also develop abundant genetic resources (cp hotspots, cpSSRs, polymorphic gSSRs). The genomic patterns and genetic resources presented here will contribute to further studies on population genetics, phylogeny and conservation biology in Saxifragaceae.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4633-x) contains supplementary material, which is available to authorized users.

Complete plastome sequencing from Toona (Meliaceae) and phylogenomic analyses within Sapindales

PREMISE OF THE STUDY

Toona (Meliaceae, Sapindales) is a small genus of five species of trees native from southern and eastern Asia to New Guinea and Australia. Complete plastomes were sequenced for three Toona species to provide a basis for future plastome genetic studies in threatened species of Toona. In addition, plastome structural evolution and phylogenetic relationships across Sapindales were explored with a larger data set of 29 Sapindales plastomes (including members of six out of nine families).

METHODS

The plastomes were determined using the Illumina sequencing platform; the phylogenetic analyses were conducted using maximum likelihood by RAxML.

RESULTS

The lengths of three Toona plastomes range from 159,185 to 158,196 bp. A total of 113 unique genes were found in each plastome. Across Sapindales, plastome gene structure and content were largely conserved, with the exception of the contraction of the inverted repeat region to exclude ycf1 in some species of Rutaceae and Sapindaceae, and the movement of trnI-GAU and trnA-UGC to a position outside the inverted repeat region in some Rutaceae species.

DISCUSSION

The three Toona plastomes possess the typical structure of angiosperm plastomes. Phylogenomic analysis of Sapindales recovered a mostly strongly supported phylogeny of Sapindales, including most of the backbone relationships, with some improvements compared to previous targeted-gene analyses.

Complete chloroplast genome sequence and comparative analysis of loblolly pine (Pinus taeda L.) with related species

Pinaceae, the largest family of conifers, has a diversified organization of chloroplast (cp) genomes with two typical highly reduced inverted repeats (IRs). In the current study, we determined the complete sequence of the cp genome of an economically and ecologically important conifer tree, the loblolly pine (Pinus taeda L.), using Illumina paired-end sequencing and compared the sequence with those of other pine species. The results revealed a genome size of 121,531 base pairs (bp) containing a pair of 830-bp IR regions, distinguished by a small single copy (42,258 bp) and large single copy (77,614 bp) region. The chloroplast genome of P. taeda encodes 120 genes, comprising 81 protein-coding genes, four ribosomal RNA genes, and 35 tRNA genes, with 151 randomly distributed microsatellites. Approximately 6 palindromic, 34 forward, and 22 tandem repeats were found in the P. taeda cp genome. Whole cp genome comparison with those of other Pinus species exhibited an overall high degree of sequence similarity, with some divergence in intergenic spacers. Higher and lower numbers of indels and single-nucleotide polymorphism substitutions were observed relative to P. contorta and P. monophylla, respectively. Phylogenomic analyses based on the complete genome sequence revealed that 60 shared genes generated trees with the same topologies, and P. taeda was closely related to P. contorta in the subgenus Pinus. Thus, the complete P. taeda genome provided valuable resources for population and evolutionary studies of gymnosperms and can be used to identify related species.

Application of Chloroplast Phylogenomics to Resolve Species Relationships Within the Plant Genus Amaranthus

Amaranthus species are an emerging and promising nutritious traditional vegetable food source. Morphological plasticity and poorly resolved dendrograms have led to the need for well resolved species phylogenies. We hypothesized that whole chloroplast phylogenomics would result in more reliable differentiation between closely related amaranth species. The aims of the study were therefore: to construct a fully assembled, annotated chloroplast genome sequence of Amaranthus tricolor; to characterize Amaranthus accessions phylogenetically by comparing barcoding genes (matK, rbcL, ITS) with whole chloroplast sequencing; and to use whole chloroplast phylogenomics to resolve deeper phylogenetic relationships. We generated a complete A. tricolor chloroplast sequence of 150,027 bp. The three barcoding genes revealed poor inter- and intra-species resolution with low bootstrap support. Whole chloroplast phylogenomics of 59 Amaranthus accessions increased the number of parsimoniously informative sites from 92 to 481 compared to the barcoding genes, allowing improved separation of amaranth species. Our results support previous findings that two geographically independent domestication events of Amaranthus hybridus likely gave rise to several species within the Hybridus complex, namely Amaranthus dubius, Amaranthus quitensis, Amaranthus caudatus, Amaranthus cruentus and Amaranthus hypochondriacus. Poor resolution of species within the Hybridus complex supports the recent and ongoing domestication within the complex, and highlights the limitation of chloroplast data for resolving recent evolution. The weedy Amaranthus retroflexus and Amaranthus powellii was found to share a common ancestor with the Hybridus complex. Leafy amaranth, Amaranthus tricolor, Amaranthus blitum, Amaranthus viridis and Amaranthus graecizans formed a stable sister lineage to the aforementioned species across the phylogenetic trees. This study demonstrates the power of next-generation sequencing data and reference-based assemblies to resolve phylogenies, and also facilitated the identification of unknown Amaranthus accessions from a local genebank. The informative phylogeny of the Amaranthus genus will aid in selecting accessions for breeding advanced genotypes to satisfy global food demand.

The complete chloroplast genome sequence of Epipremnum aureum and its comparative analysis among eight Araceae species

Epipremnum aureum is an important foliage plant in the Araceae family. In this study, we have sequenced the complete chloroplast genome of E. aureum by using Illumina Hiseq sequencing platforms. This genome is a double-stranded circular DNA sequence of 164,831 bp that contains 35.8% GC. The two inverted repeats (IRa and IRb; 26,606 bp) are spaced by a small single-copy region (22,868 bp) and a large single-copy region (88,751 bp). The chloroplast genome has 131 (113 unique) functional genes, including 86 (79 unique) protein-coding genes, 37 (30 unique) tRNA genes, and eight (four unique) rRNA genes. Tandem repeats comprise the majority of the 43 long repetitive sequences. In addition, 111 simple sequence repeats are present, with mononucleotides being the most common type and di- and tetranucleotides being infrequent events. Positive selection pressure on rps12 in the E. aureum chloroplast has been demonstrated via synonymous and nonsynonymous substitution rates and selection pressure sites analyses. Ycf15 and infA are pseudogenes in this species. We constructed a Maximum Likelihood phylogenetic tree based on the complete chloroplast genomes of 38 species from 13 families. Those results strongly indicated that E. aureum is positioned as the sister of Colocasia esculenta within the Araceae family. This work may provide information for further study of the molecular phylogenetic relationships within Araceae, as well as molecular markers and breeding novel varieties by chloroplast genetic-transformation of E. aureum in particular.

Electron Transfer through the Acceptor Side of Photosystem I: Interaction with Exogenous Acceptors and Molecular Oxygen

This review considers the state-of-the-art on mechanisms and alternative pathways of electron transfer in photosynthetic electron transport chains of chloroplasts and cyanobacteria. The mechanisms of electron transport control between photosystems (PS) I and II and the Calvin-Benson cycle are considered. The redistribution of electron fluxes between the noncyclic, cyclic, and pseudocyclic pathways plays an important role in the regulation of photosynthesis. Mathematical modeling of light-induced electron transport processes is considered. Particular attention is given to the electron transfer reactions on the acceptor side of PS I and to interactions of PS I with exogenous acceptors, including molecular oxygen. A kinetic model of PS I and its interaction with exogenous electron acceptors has been developed. This model is based on experimental kinetics of charge recombination in isolated PS I. Kinetic and thermodynamic parameters of the electron transfer reactions in PS I are scrutinized. The free energies of electron transfer between quinone acceptors A_1A/A_1B in the symmetric redox cofactor branches of PS I and iron-sulfur clusters F_X, F_A, and F_B have been estimated. The second-order rate constants of electron transfer from PS I to external acceptors have been determined. The data suggest that byproduct formation of superoxide radical in PS I due to the reduction of molecular oxygen in the A₁ site (Mehler reaction) can exceed 0.3% of the total electron flux in PS I.

Experimental reconstruction of double-stranded break repair-mediated plastid DNA insertion into the tobacco nucleus

The mitochondria and plastids of eukaryotic cells evolved from endosymbiotic prokaryotes. DNA from the endosymbionts has bombarded nuclei since the ancestral prokaryotes were engulfed by a precursor of the nucleated eukaryotic host. An experimental confirmation regarding the molecular mechanisms responsible for organelle DNA incorporation into nuclei has not been performed until the present analysis. Here we introduced double-stranded DNA breaks into the nuclear genome of tobacco through inducible expression of I-SceI, and showed experimentally that tobacco chloroplast DNAs insert into nuclear genomes through double-stranded DNA break repair. Microhomology-mediated linking of disparate segments of chloroplast DNA occurs frequently during healing of induced nuclear double-stranded breaks (DSB) but the resulting nuclear integrants are often immediately unstable. Non-Mendelian inheritance of a selectable marker (neo), used to identify plastid DNA transfer, was observed in the progeny of about 50% of lines emerging from the screen. The instability of these de novo nuclear insertions of plastid DNA (nupts) was shown to be associated with deletion not only of the nupt itself but also of flanking nuclear DNA within one generation of transfer. This deletion of pre-existing nuclear DNA suggests that the genetic impact of organellar DNA transfer to the nucleus is potentially far greater than previously thought.

Independent degradation in genes of the plastid ndh gene family in species of the orchid genus Cymbidium (Orchidaceae; Epidendroideae)

In this paper, we compare ndh genes in the plastid genome of many Cymbidium species and three closely related taxa in Orchidaceae looking for evidence of ndh gene degradation. Among the 11 ndh genes, there were frequently large deletions in directly repeated or AT-rich regions. Variation in these degraded ndh genes occurs between individual plants, apparently at population levels in these Cymbidium species. It is likely that ndh gene transfers from the plastome to mitochondrial genome (chondriome) occurred independently in Orchidaceae and that ndh genes in the chondriome were also relatively recently transferred between distantly related species in Orchidaceae. Four variants of the ycf1-rpl32 region, which normally includes the ndhF genes in the plastome, were identified, and some Cymbidium species contained at least two copies of that region in their organellar genomes. The four ycf1-rpl32 variants seem to have a clear pattern of close relationships. Patterns of ndh degradation between closely related taxa and translocation of ndh genes to the chondriome in Cymbidium suggest that there have been multiple bidirectional intracellular gene transfers between two organellar genomes, which have produced different levels of ndh gene degradation among even closely related species.

Chloroplast genomes of Lilium lancifolium, L. amabile, L. callosum, and L. philadelphicum: Molecular characterization and their use in phylogenetic analysis in the genus Lilium and other allied genera in the order Liliales

Chloroplast (cp) genomes of Lilium amabile, L. callosum, L. lancifolium, and L. philadelphicum were fully sequenced. Using these four novel cp genome sequences and five other previously sequenced cp genomes, features of the cp genomes were characterized in detail among species in the genus Lilium and other related genera in the order Liliales. The lengths and nucleotide composition showed little variation. No structural variation was found among the cp genomes in Liliales. Gene contents were conserved among four newly sequenced cp genome in Lilium species, the only differences being in two pseudogenes. We identified 112 genes in 13 functional categories, 18 of which carried introns that were conserved among the species in Liliales. There were 16-21 SSR loci (>12 bp, >3 repeats) in the cp genomes in Lilium and the genomic locations of these loci were highly variable among the species. Average mutations were 15 SNPs per 1kb and 5 indels per 1kb, respectively, in the cp genomes of the newly sequenced four Lilium species. Phylogenetic classifications revealed some discrepancies between trees based on the cp genomes and previous classifications based on the morphology and geographic distributions.

The complete plastome of tropical fruit Garcinia mangostana (Clusiaceae)

The complete plastome sequence of Garcinia mangostana L. (Clusiaceae) is completed in this study (NCBI acc. no. KX822787). This is a first complete plastome sequence from the Clusiaceae. The complete plastome size is 158,179 bp in length and consists of a large single copy of 86,458 bp and a small single copy of 17,703 bp, separated by two inverted repeats of 27,009 bp. The G. mangostana plastome shows four minor structural modifications including infA gene loss, rpl32 gene loss, ycf3 gene intron loss and a 363 bp inversion between trnV-UAC and atpE gene. The plastome contains 111 genes, of which 77 are protein-coding genes, 30 are tRNA genes and four are rRNA genes. The average A-T content of the plastome is 63.9%. A total of 110 simple sequence loci are identified from the genome. Phylogenetic analysis reveals that G. mangostana is a sister group of Erythroxylum novogranatense (Erythroxylaceae) with 78% bootstrap support.

The complete plastome sequence of the endangered orchid Cymbidium macrorhizon (Orchidaceae)

In this study, we determined the complete chloroplast sequence of Cymbidium macrorhizon Lindl. (Orchidaceae) (NCBI acc. no. KY354040), an endangered plant species protected by the national law of Korea. The gene order and number in the C. macrorhizon plastome were similar to a typical Orchid plastome. The complete plastome was 149,859 bp in length and consisted of a large single copy region of 85,187 bp and a small single copy region of 13,766 bp; these were separated by two inverted repeats of 25,453 bp. The plastome contained 103 genes of which 69 were protein-coding genes, 30 were tRNA genes and four were rRNA genes. Fourteen genes contained one intron and two genes (clpP, and ycf3) had two introns. The AT content of the C. macrorhizon plastome was 60.0% and a total of 62 simple sequence repeat regions were identified in the plastome. Phylogenetic analysis also identified C. lancifolium as a closely related sister to C. macrorhizon, suggesting that the hemisaprophytic nature of C. macrorhizon is derived recently from a common leafy ancestor.

The complete plastome sequence of the endangered orchid Oberonia japonica (Orchidaceae)

In this study, we report the complete chloroplast sequence of Oberonia japonica (Maxim.) Makino (Orchidaceae) (NCBI acc. no. KX871235), which is an endangered plant species currently protected by the National Law of Korea. The gene order and content of the O. japonica plastome are similar to those of a typical orchid plastome. The 11 ndh genes are pseudogenized or lost completely from the plastome of O. japonica. The plastome contains 102 genes, of which 68 are protein-coding genes, 30 and four are transfer RNA (tRNA) and ribosomal RNA (rRNA) genes, respectively. Sixteen genes contain one intron and two genes (clpP and ycf3) have two introns. The complete plastome is 142,996 bp long and consists of one large and small single copy each of 81,669 and 10,969 bp, respectively, separated by two inverted repeats of 25,179 bp. The AT content of the O. japonica plastome is 62.6%. Sixty-five simple sequence repeat (SSR) loci, consisting of 56 mono-SSR, seven di-SSR, and two tri-SSR are scattered along the O. japonica plastome. Some of these plastome SSR loci may be useful for developing genetic markers for the O. japonica populations. Phylogenetic analysis has shown that O. japonica is a member of the Malaxideae. The genus Oberonia forms a monophyletic clade with Dendrobium.

The complete plastome sequence of the endangered orchid Kuhlhasseltia nakaiana (Orchidaceae)

In this study, we report the complete plastome sequence of Kuhlhasseltia nakaiana (F.Maek.) Ormerod (Orchidaceae) (NCBI acc. no. KY354041), an endangered plant species protected by the national law of Korea. The gene order and number in the K. nakaiana plastome were similar to a typical orchid plastome. The complete plastome was 147,614 bp in length and consisted of a large single copy region of 81,617 bp and a small single copy region of 13,673 bp, separated by two inverted repeats of 26,162 bp. The plastome contained 103 genes, of which 69 were protein-coding genes, 30 were tRNA genes, and four were rRNA genes. Fourteen genes contained one intron and two genes (clpP and ycf3) had two introns. The AT content of the plastome was 60.5%. A total of 74 simple sequence repeat regions were identified from the plastome. Phylogenetic analysis determined that K. nakaiana was a member of the tribe Cranichideae and revealed the sister group relationship between K. nakaiana and Ludisia discolor within the tribe Cranichideae.

The complete plastome sequences of Mangifera indica L. (Anacardiaceae)

In this study, we determined the complete plastome sequence of Mangifera indica L. (Anacardiaceae) (NCBI acc. no. KX871231). The complete plastome is 157,780 bp in length, and consists of a large single copy of 86,673 bp and a small single copy of 18,349 bp, separated by two inverted repeats of 25,792 bp. The plastome contains 112 genes, of which 78 are protein-coding genes, 30 are tRNA genes, and four are rRNA genes. Sixteen genes contain one intron and two genes have two introns. The average A-T content of the plastome is 62.1%. The M. indica plastome has approximately 15 kb inversion between trnT-UGU and trnT-GGU. We identify a palindromic repeat of 18 bp (ATTCTTTTTTTTTTTTTT/AAAAAAAAAAAAAAGAAT) near the inversion breakpoints of M. indica plastome. Phylogenetic analysis revealed that M. indica is a sister group of Rhus chinensis with 100% bootstrap support. Anacardiaceae clade is a sister group of Boswellia sacra (Burseraceae) with 100% bootstrap support.

The complete plastome sequence of the endangered orchid Habenaria radiata (Orchidaceae)

In this study, we determined the complete chloroplast sequence of Habenaria radiata (Thunb.) Spreng. (Orchidaceae) (NCBI acc. no. KX871237), an endangered plant species protected by the national law of Korea. The gene order and gene content of the H. radiata plastome are similar to those of typical angiosperm plastomes. The 11 ndh genes, which are usually lost in orchid plastomes, are intact in the H. radiata plastome. The complete plastome is 155,353 bp in length and consists of a large single copy of 84,833 bp and a small single copy of 17,718 bp, separated by two inverted repeats of 26,401 bp. The plastome contains 113 genes, of which 79 are protein-coding genes, 30 are tRNA genes, and four are rRNA genes. Sixteen genes contain one intron and two genes (clpP, ycf3) have two introns. A total of 76 simple sequence repeat (SSR) loci, which consist of 58 mono-SSR, 17 di-SSR, and 1 tri-SSR, are scattered along the H. radiata plastome. Some of these plastome SSR and high sequence divergent regions may be useful for development of genetic markers for the populations of H. radiata and other congeneric species. Phylogenetic analysis identified the sister relationship between H. radiata and H. pantlingiana within the tribe Orchideae.

Celebrating wobble decoding: Half a century and still much is new

A simple post-transcriptional modification of tRNA, deamination of adenosine to inosine at the first, or wobble, position of the anticodon, inspired Francis Crick's Wobble Hypothesis 50 years ago. Many more naturally-occurring modifications have been elucidated and continue to be discovered. The post-transcriptional modifications of tRNA's anticodon domain are the most diverse and chemically complex of any RNA modifications. Their contribution with regards to chemistry, structure and dynamics reveal individual and combined effects on tRNA function in recognition of cognate and wobble codons. As forecast by the Modified Wobble Hypothesis 25 years ago, some individual modifications at tRNA's wobble position have evolved to restrict codon recognition whereas others expand the tRNA's ability to read as many as four synonymous codons. Here, we review tRNA wobble codon recognition using specific examples of simple and complex modification chemistries that alter tRNA function. Understanding natural modifications has inspired evolutionary insights and possible innovation in protein synthesis.

The complete chloroplast genome of Cinnamomum camphora and its comparison with related Lauraceae species

Cinnamomum camphora, a member of the Lauraceae family, is a valuable aromatic and timber tree that is indigenous to the south of China and Japan. All parts of Cinnamomum camphora have secretory cells containing different volatile chemical compounds that are utilized as herbal medicines and essential oils. Here, we reported the complete sequencing of the chloroplast genome of Cinnamomum camphora using illumina technology. The chloroplast genome of Cinnamomum camphora is 152,570 bp in length and characterized by a relatively conserved quadripartite structure containing a large single copy region of 93,705 bp, a small single copy region of 19,093 bp and two inverted repeat (IR) regions of 19,886 bp. Overall, the genome contained 123 coding regions, of which 15 were repeated in the IR regions. An analysis of chloroplast sequence divergence revealed that the small single copy region was highly variable among the different genera in the Lauraceae family. A total of 40 repeat structures and 83 simple sequence repeats were detected in both the coding and non-coding regions. A phylogenetic analysis indicated that Calycanthus is most closely related to Lauraceae, both being members of Laurales, which forms a sister group to Magnoliids. The complete sequence of the chloroplast of Cinnamomum camphora will aid in in-depth taxonomical studies of the Lauraceae family in the future. The genetic sequence information will also have valuable applications for chloroplast genetic engineering.

The complete chloroplast genome sequence of Aconitum coreanum and Aconitum carmichaelii and comparative analysis with other Aconitum species

Aconitum species (belonging to the Ranunculaceae) are well known herbaceous medicinal ingredients and have great economic value in Asian countries. However, there are still limited genomic resources available for Aconitum species. In this study, we sequenced the chloroplast (cp) genomes of two Aconitum species, A. coreanum and A. carmichaelii, using the MiSeq platform. The two Aconitum chloroplast genomes were 155,880 and 157,040 bp in length, respectively, and exhibited LSC and SSC regions separated by a pair of inverted repeat regions. Both cp genomes had 38% GC content and contained 131 unique functional genes including 86 protein-coding genes, eight ribosomal RNA genes, and 37 transfer RNA genes. The gene order, content, and orientation of the two Aconitum cp genomes exhibited the general structure of angiosperms, and were similar to those of other Aconitum species. Comparison of the cp genome structure and gene order with that of other Aconitum species revealed general contraction and expansion of the inverted repeat regions and single copy boundary regions. Divergent regions were also identified. In phylogenetic analysis, Aconitum species positon among the Ranunculaceae was determined with other family cp genomes in the Ranunculales. We obtained a barcoding target sequence in a divergent region, ndhC–trnV, and successfully developed a SCAR (sequence characterized amplified region) marker for discrimination of A. coreanum. Our results provide useful genetic information and a specific barcode for discrimination of Aconitum species.

Strategies for complete plastid genome sequencing

Plastid sequencing is an essential tool in the study of plant evolution. This high-copy organelle is one of the most technically accessible regions of the genome, and its sequence conservation makes it a valuable region for comparative genome evolution, phylogenetic analysis and population studies. Here, we discuss recent innovations and approaches for de novo plastid assembly that harness genomic tools. We focus on technical developments including low-cost sequence library preparation approaches for genome skimming, enrichment via hybrid baits and methylation-sensitive capture, sequence platforms with higher read outputs and longer read lengths, and automated tools for assembly. These developments allow for a much more streamlined assembly than via conventional short-range PCR. Although newer methods make complete plastid sequencing possible for any land plant or green alga, there are still challenges for producing finished plastomes particularly from herbarium material or from structurally divergent plastids such as those of parasitic plants.

Characterization of the complete chloroplast genome of Arabis stellari and comparisons with related species

Arabis stellari var. japonica is an ornamental plant of the Brassicaceae family, and is widely distributed in South Korea. However, no information is available about its molecular biology and no genomic study has been performed on A. stellari. In this paper, the authors report the complete chloroplast genome sequence of A. stellari. The plastome of A. stellari was 153,683 bp in length with 36.4% GC and included a pair of inverted repeats (IRs) of 26,423 bp that separated a large single-copy (LSC) region of 82,807 bp and a small single-copy (SSC) region of 18,030 bp. It was also found to contain 113 unique genes, of which 79 were protein-coding genes, 30 were transfer RNAs, and four were ribosomal RNAs. The gene content and organization of the A. stellari chloroplast genome were similar to those of other Brassicaceae genomes except for the absence of the rps16 protein-coding gene. A total of 991 SSRs were identified in the genome. The chloroplast genome of A. stellari was compared with closely related species of the Brassicaceae family. Comparative analysis showed a minor divergence occurred in the protein-coding matK, ycf1, ccsA, accD and rpl22 genes and that the K_A/K_S nucleotide substitution ratio of the ndhA genes of A. stellari and A. hirsuta was 1.35135. The genes infA and rps16 were absent in the Arabis genus and phylogenetic evolutionary studies revealed that these genes evolved independently. However, phylogenetic analysis showed that the positions of Brassicaceae species are highly conserved. The present study provides A. stellari genomic information that may be found useful in conservation and molecular phylogenetic studies on Brassicaceae.

Complete plastome sequencing of both living species of Circaeasteraceae (Ranunculales) reveals unusual rearrangements and the loss of the ndh gene family

BACKGROUND

Among the 13 families of early-diverging eudicots, only Circaeasteraceae (Ranunculales), which consists of the two monotypic genera Circaeaster and Kingdonia, lacks a published complete plastome sequence. In addition, the phylogenetic position of Circaeasteraceae as sister to Lardizabalaceae has only been weakly or moderately supported in previous studies using smaller data sets. Moreover, previous plastome studies have documented a number of novel structural rearrangements among early-divergent eudicots. Hence it is important to sequence plastomes from Circaeasteraceae to better understand plastome evolution in early-diverging eudicots and to further investigate the phylogenetic position of Circaeasteraceae.

RESULTS

Using an Illumina HiSeq 2000, complete plastomes were sequenced from both living members of Circaeasteraceae: Circaeaster agrestis and Kingdonia uniflora . Plastome structure and gene content were compared between these two plastomes, and with those of other early-diverging eudicot plastomes. Phylogenetic analysis of a 79-gene, 99-taxon data set including exemplars of all families of early-diverging eudicots was conducted to resolve the phylogenetic position of Circaeasteraceae. Both plastomes possess the typical quadripartite structure of land plant plastomes. However, a large ~49 kb inversion and a small ~3.5 kb inversion were found in the large single-copy regions of both plastomes, while Circaeaster possesses a number of other rearrangements, particularly in the Inverted Repeat. In addition, infA was found to be a pseudogene and accD was found to be absent within Circaeaster, whereas all ndh genes, except for ndhE and ndhJ, were found to be either pseudogenized (ΨndhA, ΨndhB, ΨndhD, ΨndhH and ΨndhK) or absent (ndhC, ndhF, ndhI and ndhG) in Kingdonia. Circaeasteraceae was strongly supported as sister to Lardizabalaceae in phylogenetic analyses.

CONCLUSION

The first plastome sequencing of Circaeasteraceae resulted in the discovery of several unusual rearrangements and the loss of ndh genes, and confirms the sister relationship between Circaeasteraceae and Lardizabalaceae. This research provides new insight to characterize plastome structural evolution in early-diverging eudicots and to better understand relationships within Ranunculales .

Comparative analysis of complete plastid genomes from wild soybean (Glycine soja) and nine other Glycine species

The plastid genomes of different plant species exhibit significant variation, thereby providing valuable markers for exploring evolutionary relationships and population genetics. Glycine soja (wild soybean) is recognized as the wild ancestor of cultivated soybean (G. max), representing a valuable genetic resource for soybean breeding programmes. In the present study, the complete plastid genome of G. soja was sequenced using Illumina paired-end sequencing and then compared it for the first time with previously reported plastid genome sequences from nine other Glycine species. The G. soja plastid genome was 152,224 bp in length and possessed a typical quadripartite structure, consisting of a pair of inverted repeats (IRa/IRb; 25,574 bp) separated by small (178,963 bp) and large (83,181 bp) single-copy regions, with a 51-kb inversion in the large single-copy region. The genome encoded 134 genes, including 87 protein-coding genes, eight ribosomal RNA genes, and 39 transfer RNA genes, and possessed 204 randomly distributed microsatellites, including 15 forward, 25 tandem, and 34 palindromic repeats. Whole-plastid genome comparisons revealed an overall high degree of sequence similarity between G. max and G. gracilis and some divergence in the intergenic spacers of other species. Greater numbers of indels and SNP substitutions were observed compared with G. cyrtoloba. The sequence of the accD gene from G. soja was highly divergent from those of the other species except for G. max and G. gracilis. Phylogenomic analyses of the complete plastid genomes and 76 shared genes yielded an identical topology and indicated that G. soja is closely related to G. max and G. gracilis. The complete G. soja genome sequenced in the present study is a valuable resource for investigating the population and evolutionary genetics of Glycine species and can be used to identify related species.

Phylogenetic studies and comparative chloroplast genome analyses elucidate the basal position of halophyte Nitraria sibirica (Nitrariaceae) in the Sapindales

Nitraria sibirica is a halophyte and belongs to the family Nitrariaceae. The chloroplast genome of Nitraria sibirica (159,466 bp) has a quadripartite structure, which consists of a large single-copy (87,914 bp) region, a small single-copy (18,316 bp) region, and a pair of inverted repeats (26,618 bp). Sequencing analyses indicate that the chloroplast genome contains 113 distinct genes, including 79 peptide-coding genes, 30 transfer RNA genes, and 4 ribosomal RNA genes. We also identified 105 perfect simple sequence repeats, 12 most divergent non-coding regions, and 6 most divergent coding regions when compared to the chloroplast genomes of the Sapindales plants. Phylogenetic analyses using the concatenated amino acid sequences of 58 protein-coding genes from 48 species suggest that the 'basal' position of Nitraria sibirica belongs to the Sapindales clade. We also found that the inverted repeat expansion resulted in a duplication of rps19 in Nitraria sibirica when comparing its chloroplast genome structure with Theobroma cacao, Vitis vinifera, Eucalyptus erythrocorys and Arabidopsis thaliana. The duplication of rps19 gene was consistent with that in Zanthoxylum piperitum, Azadirachta indica, Sapindus mukorossi and Citrus sinensis, all of which belong to the order Sapindales, but different from most Rosids plants. In summary, the analyses of Nitraria sibirica chloroplast genome not only provide insights into comparative genome analysis, but also pave the way for a better understanding of the phylogenetic relationships within the Sapindales.

Alternative electron transport pathways in photosynthesis: a confluence of regulation

Photosynthetic reactions proceed along a linear electron transfer chain linking water oxidation at photosystem II (PSII) to CO₂ reduction in the Calvin-Benson-Bassham cycle. Alternative pathways poise the electron carriers along the chain in response to changing light, temperature and CO₂ inputs, under prolonged hydration stress and during development. We describe recent literature that reports the physiological functions of new molecular players. Such highlights include the flavodiiron proteins and their important role in the green lineage. The parsing of the proton-motive force between ΔpH and Δψ, regulated in many different ways (cyclic electron flow, ATPsynthase conductivity, ion/H⁺ transporters), is comprehensively reported. This review focuses on an integrated description of alternative electron transfer pathways and how they contribute to photosynthetic productivity in the context of plant fitness to the environment.

EVOLUTIONARY SIGNIFICANCE OF THE LOSS OF THE CHLOROPLAST-DNA INVERTED REPEAT IN THE LEGUMINOSAE SUBFAMILY PAPILIONOIDEAE

The distribution of a rare chloroplast-DNA structural mutation, the loss of a large inverted repeat, has been determined for 95 species representing 77 genera and 25 of the 31 tribes in the legume subfamily Papilionoideae. This mutation, which is regarded as a derived feature of singular origin within the subfamily, marks a group comprising six temperate tribes, the Galegeae, Hedysareae, Carmichaelieae, Vicieae, Cicereae, and Trifolieae, an assemblage traditionally considered to be monophyletic. This mutation also occurs in the chloroplast genome of Wisteria, a member of the tropical tribe Millettieae whose other members so far surveyed lack the mutation. These new DNA data, together with traditional evidence, support the hypothesis that Wisteria is an unspecialized member of a lineage that gave rise to the temperate tribes marked by the chloroplast-DNA mutation; the probable paraphylesis of Millettieae is revealed. Two other tribes, Loteae and Coronilleae (traditionally regarded as a derived element of the aforesaid temperate tribes) do not possess this chloroplast-DNA structural mutation and, therefore, presumably represent a distinct temperate lineage. This hypothesis is supported by additional evidence from pollen, inflorescence, and root-nodule morphology that suggests that the Loteae and Coronilleae share a more recent ancestry with tropical tribes such as Phaseoleae and Millettieae than with other temperate tribes.

SIX INDEPENDENT LOSSES OF THE CHLOROPLAST DNA rpl2 INTRON IN DICOTYLEDONS: MOLECULAR AND PHYLOGENETIC IMPLICATIONS

Previous studies have shown that in several angiosperms and the liverwort Marchantia the chloroplast gene rpl2, encoding ribosomal protein L2, is interrupted by an intron, but that in spinach (Spinacia oleracea, Caryophyllales) this intron has been lost. We have determined the distribution of the rpl2 intron for 390 species representing 116 angiosperm families. Filter hybridizations reveal that the intron is absent from the chloroplast genomes of all examined families of the Caryophyllales, suggesting that the intron was lost in the common ancestor of the order. Sequencing of the rpl2 gene in five genera of the Caryophyllales and in Rumex (Polygonales) not only confirms the filter hybridization results, but also shows that for all taxa lacking the intron, the rpl2 gene has undergone a precise deletion of the intron. In all cases, it is the original rpl2 gene that has sustained loss of its intron. This implies that in chloroplast DNA, integration of exogenous genes (e.g., a reverse transcript of a spliced mRNA) occurs mainly by homologous, replacement recombination, rather than by illegitimate recombination elsewhere in the genome. Filter hybridizations also reveal that the rpl2 intron was lost independently in the common ancestors of at least five other lineages of dicotyledons: Saxifragaceae (s.s.), Convolvulaceae (including Cuscuta), Menyanthaceae, two genera of Geraniaceae, and one genus of Droseraceae. The molecular and phylogenetic implications of these independent intron losses are discussed.