Complete chloroplast genome sequence of MD-2 pineapple and its comparative analysis among nine other plants from the subclass Commelinidae

Plastomes of Garcinia mangostana L. and Comparative Analysis with Other Garcinia Species

The two varieties of mangosteen (Garcinia mangostana L.) cultivated in Malaysia are known as Manggis and Mesta. The latter is preferred for its flavor, texture, and seedlessness. Here, we report a complete plastome (156,580 bp) of the Mesta variety that was obtained through a hybrid assembly approach using PacBio and Illumina sequencing reads. It encompasses a large single-copy (LSC) region (85,383 bp) and a small single-copy (SSC) region (17,137 bp) that are separated by 27,230 bp of inverted repeat (IR) regions at both ends. The plastome comprises 128 genes, namely, 83 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. The plastome of the Manggis variety (156,582 bp) obtained from reference-guided assembly of Illumina reads was found to be nearly identical to Mesta except for two indels and the presence of a single-nucleotide polymorphism (SNP). Comparative analyses with other publicly available Garcinia plastomes, including G. anomala, G. gummi-gutta, G. mangostana var. Thailand, G. oblongifolia, G. paucinervis, and G. pedunculata, found that the gene content, gene order, and gene orientation were highly conserved among the Garcinia species. Phylogenomic analysis divided the six Garcinia plastomes into three groups, with the Mesta and Manggis varieties clustered closer to G. anomala, G. gummi-gutta, and G. oblongifolia, while the Thailand variety clustered with G. pedunculata in another group. These findings serve as future references for the identification of species or varieties and facilitate phylogenomic analysis of lineages from the Garcinia genus to better understand their evolutionary history.

The population genomic analyses of chloroplast genomes shed new insights on the complicated ploidy and evolutionary history in Fragaria

The genus Fragaria consists of a rich diversity of ploidy levels with diploid (2x), tetraploid (4x), pentaploid (5x), hexaploidy (6x), octoploid (8x) and decaploid (10x) species. Only a few studies have explored the origin of diploid and octoploid strawberry, and little is known about the roles of tetraploidy and hexaploidy during the evolution of octoploid strawberry. The chloroplast genome is usually a stable circular genome and is widely used in investigating the evolution and matrilineal identification. Here, we assembled the chloroplast genomes of F. x ananassa cv. 'Benihoppe' (8x) using Illumina and HiFi data seperately. The genome alignment results showed that more InDels were located in the chloroplast genomes based on the PacBio HiFi data than Illumina data. We obtain highly accurate chloroplast genomes assembled through GetOrganelle using Illumina reads. We assembled 200 chloroplast genomes including 198 Fragaria (21 species) and 2 Potentilla samples. Analyses of sequence variation, phylogenetic and PCA analyses showed that Fragaria was divided into five groups. F. iinumae, F. nilgerrensis and all octoploid accessions formed Group A, C and E separately. Species native to western China were clustered into Group B. Group D consisted of F. virdis, F. orientalis, F. moschata, and F. vesca. STRUCTURE and haplotype network confirmed that the diploid F. vesca subsp. bracteata was the last maternal donator of octoploid strawberry. The dN/dS ratio estimated for the protein-coding genes revealed that genes involved in ATP synthase and photosystem function were under positive selection. These findings demonstrate the phylogeny of totally 21 Fragaria species and the origin of octoploid species. F. vesca was the last female donator of octoploid, which confirms the hypothesis that the hexaploid species F. moschata may be an evolutionary intermediate between the diploids and wild octoploid species.

The complete chloroplast genome of Hibiscus syriacus using long-read sequencing: Comparative analysis to examine the evolution of the tribe Hibisceae

Hibiscus syriacus, a member of the tribe Hibisceae, is considered an important ornamental and medicinal plant in east Asian countries. Here, we sequenced and assembled the complete chloroplast genome of H. syriacus var. Baekdansim using the PacBio long-read sequencing platform. A quadripartite structure with 161,026 base pairs was obtained, consisting of a pair of inverted repeats (IRA and IRB) with 25,745 base pairs, separated by a large single-copy region of 89,705 base pairs and a short single-copy region of 19,831 base pairs. This chloroplast genome had 79 protein-coding genes, 30 transfer RNA genes, 4 ribosomal RNA genes, and 109 simple sequence repeat regions. Among them, ndhD and rpoC1, containing traces of RNA-editing events associated with adaptive evolution, were identified by analysis of putative RNA-editing sites. Codon usage analysis revealed a preference for A/U-terminated codons. Furthermore, the codon usage pattern had a clustering tendency similar to that of the phylogenetic analysis of the tribe Hibisceae. This study provides clues for understanding the relationships and refining the taxonomy of the tribe Hibisceae.

Plastid Phylogenomic Analyses Reveal the Taxonomic Position of Peucedanum franchetii

Peucedanum franchetii is a famous folk medicinal plant in China. However, the taxonomy of the P. franchetii has not been sufficiently resolved. Due to similar morphological features between P. franchetii and Ligusticopsis members, the World Flora Online (WFO) Plant List suggested that this species transformed into the genus Ligusticopsis and merged with Ligusticopsis likiangensis. However, both species are obviously diverse in leaf shape, bracts, and bracteoles. To check the taxonomic position of P. franchetii, we newly sequenced and assembled the plastome of P. franchetii and compared it with nine other plastomes of the genus Ligusticopsis. Ten plastomes were highly conserved and similar in gene order, codon bias, RNA editing sites, IR borders, and SSRs. Nevertheless, 10 mutation hotspot regions (infA, rps8, matK, ndhF, rps15, psbA-trnH, rps2-rpoC2, psbA-trnK, ycf2-trnL, and ccsA-ndhD) were still detected. In addition, both phylogenetic analyses based on plastome data and ITS sequences robustly supported that P. franchetii was not clustered with members of Peucedanum but nested in Ligusticopsis. P. franchetii was sister to L. likiangensis in the ITS topology but clustered with L. capillacea in the plastome tree. These findings implied that P. franchetii should be transferred to genus Ligusticopsis and not merged with L. likiangensis, but as an independent species, which was further verified by morphological evidences. Therefore, transferring P. franchetii under the genus Ligusticopsis as an independent species was reasonable, and a new combination was presented.

Dynamic evolution of the plastome in the Elm family (Ulmaceae)

This study compared the plastomes of Ulmaceae allowing analyses of the dynamic evolution, including genome structure, codon usage bias, repeat sequences, molecular mutation rates, and phylogenetic inferences. Ulmaceae is a small family in the order Rosales. This family consists of seven genera, including Ulmus, Zelkova, Planera, Hemiptelea, Phyllostylon, Ampelocera, and Holoptelea. Ulmaceae is an interesting lineage from plant biogeographic, systematic, evolutionary, and paleobotanic perspectives. It is also a good model to investigate the evolution of the plastomes in woody plants. In this study, we sequenced and assembled the complete plastomes of the six Ulmaceae genera to compare genomic structures and reveal the molecular evolutionary patterns. The size of the quadripartite plastomes ranged from 158,290 bp to 161,886 bp. The genomes contained 131 genes, including 87 coding genes, 36 tRNA, and 8 rRNA. The gene number, gene content, and genomic structure were highly consistent among the Ulmaceae genera. Nine variable regions including ndhA intron, ndhF-rpl32, ycf1, psbK-trnS, rps16-trnQ, trnT-trnL, trnT-psbD, trnS-trnG, and rpl32-trnL, were identified in Ulmaceae plastomes according to the nucleotide diversity values. Condon usage was biased among the genes and showed consistent trends in the seven genera. Molecular evolution analyses revealed that most of the genes and all gene groups were under widespread purifying selection. Twelve genes (ccsA, matK, psbH, psbK, rbcL, rpl22, rpl32, rpoA, rps12, rps15, rps16, and ycf2) were under positive selection. Phylogenetic analyses supported that Ulmaceae should be divided into two main clades, such as the temperate clade, including Ulmus, Zelkova, Planera, and Hemiptelea and the tropical clade, including Phyllostylon, Ampelocera and Holoptelea. This study reports the structure and evolutionary characteristics of the Elm family. These new genomic data will benefit assessments of genomic evolution and provide information to elucidate the phylogenetic relationships among Ulmaceae species.

A reappraisal of the phylogeny and historical biogeography of Sparganium (Typhaceae) using complete chloroplast genomes

BACKGROUND

Sparganium (Typhaceae) is a widespread temperate genus of ecologically important aquatic plants. Previous reconstructions of the phylogenetic relationships among Sparganium species are incompletely resolved partly because they were based on molecular markers comprising < 7,000 bp. Here, we sequenced and assembled the complete chloroplast genomes from 19 Sparganium samples representing 15 putative species and three putative subspecies in order to explore chloroplast genome evolution in this genus, clarify taxonomic lineages, estimate the divergence times of Sparganium species, and reconstruct aspects of the biogeographic history of the genus.

RESULTS

The 19 chloroplast genomes shared a conserved genome structure, gene content, and gene order. Our phylogenomic analysis presented a well-resolved phylogeny with robust support for most clades. Non-monophyly was revealed in three species: S. erectum, S. eurycarpum, and S. stoloniferum. Divergence time estimates suggest that the two subgenera of Sparganium split from each other ca. 30.67 Ma in the middle Oligocene. The subgenus Xanthosparganium diversified in the late Oligocene and Miocene, while the subgenus Sparganium diversified in the late Pliocene and Pleistocene. Ancestral area reconstruction suggested that the two subgenera may have originated in East Eurasia and North America.

CONCLUSION

The non-monophyletic nature of three putative species underscores the necessity of taxonomic revision for Sparganium: S. stoloniferum subsp. choui may be more appropriately identified as S. choui, and subspecies of S. erectum may be in fact distinct species. The estimated diversification times of the two subgenera correspond to their species and nucleotide diversities. The likely ancestral area for most of subgenus Xanthosparganium was East Eurasia and North America from where it dispersed into West Eurasia and Australia. Most of subgenus Sparganium likely originated in North America and then dispersed into Eurasia. Our study demonstrates some of the ways in which complete chloroplast genome sequences can provide new insights into the evolution, phylogeny, and biogeography of the genus Sparganium.

Characterization of the Evolutionary Pressure on Anisodus tanguticus Maxim. with Complete Chloroplast Genome Sequence

Anisodus tanguticus Maxim. (Solanaceae), a traditional endangered Tibetan herb, is endemic to the Qinghai-Tibet Plateau. Here, we report the de novo assembled chloroplast (cp) genome sequences of A. tanguticus (155,765 bp). The cp contains a pair of inverted repeated (IRa and IRb) regions of 25,881 bp that are separated by a large single copy (LSC) region (86,516 bp) and a small single copy SSC (17,487 bp) region. A total of 132 functional genes were annotated in the cp genome, including 87 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. Moreover, 199 simple sequence repeats (SSR) and 65 repeat structures were detected. Comparative plastome analyses revealed a conserved gene order and high similarity of protein-coding sequences. The A. tanguticus cp genome exhibits contraction and expansion, which differs from Przewalskia tangutica and other related Solanaceae species. We identified 30 highly polymorphic regions, mostly belonging to intergenic spacer regions (IGS), which may be suitable for the development of robust and cost-effective markers for inferring the phylogeny of the genus Anisodus and family Solanaceae. Analysis of the Ka/Ks ratios of the Hyoscyameae tribe revealed significant positive selection exerted on the cemA, rpoC2, and clpP genes, which suggests that protein metabolism may be an important strategy for A. tanguticus and other species in Hyoscyameae in adapting to the adverse environment on the Qinghai-Tibetan Plateau. Phylogenetic analysis revealed that A. tanguticus clustered closer with Hyoscyamus niger than P. tangutica. Our results provide reliable genetic information for future exploration of the taxonomy and phylogenetic evolution of the Hyoscyameae tribe and related species.

Newly reported chloroplast genome of Sinosenecio albonervius Y. Liu & Q. E. Yang and comparative analyses with other Sinosenecio species

Background

Sinosenecio B. Nordenstam (Asteraceae) currently comprises 44 species. To investigate the interspecific relationship, several chloroplast markers, including ndhC-trnV, rpl32-trnL, matK, and rbcL, are used to analyze the phylogeny of Sinosenecio. However, the chloroplast genomes of this genus have not been thoroughly investigated. We sequenced and assembled the Sinosenecio albonervius chloroplast genome for the first time. A detailed comparative analysis was performed in this study using the previously reported chloroplast genomes of three Sinosenecio species.

Results

The results showed that the chloroplast genomes of four Sinosenecio species exhibit a typical quadripartite structure. There are equal numbers of total genes, protein-coding genes and RNA genes among the annotated genomes. Per genome, 49–56 simple sequence repeats and 99 repeat sequences were identified. Thirty codons were identified as RSCU values greater than 1 in the chloroplast genome of S. albonervius based on 54 protein-coding genes, indicating that they showed biased usage. Among 18 protein-coding genes, 46 potential RNA editing sites were discovered. By comparing these chloroplast genomes' structures, inverted repeat regions and coding regions were more conserved than single-copy and non-coding regions. The junctions among inverted repeat and single-copy regions showed slight difference. Several hot spots of genomic divergence were detected, which can be used as new DNA barcodes for species identification. Phylogenetic analysis of the whole chloroplast genome showed that the four Sinosenecio species have close interspecific relationships.

Conclusions

The complete chloroplast genome of Sinosenecio albonervius was revealed in this study, which included a comparison of Sinosenecio chloroplast genome structure, variation, and phylogenetic analysis for related species. These will help future research on Sinosenecio taxonomy, identification, origin, and evolution to some extent.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08872-3.

New plastome structural rearrangements discovered in core Tillandsioideae (Bromeliaceae) support recently adopted taxonomy

Full plastome sequences for land plants have become readily accessible thanks to the development of Next Generation Sequencing (NGS) techniques and powerful bioinformatic tools. Despite this vast amount of genomic data, some lineages remain understudied. Full plastome sequences from the highly diverse (>1,500 spp.) subfamily Tillandsioideae (Bromeliaceae, Poales) have been published for only three (i.e., Guzmania, Tillandsia, and Vriesea) out of 22 currently recognized genera. Here, we focus on core Tillandsioideae, a clade within subfamily Tillandsioideae, and explore the contribution of individual plastid markers and data categories to inform deep divergences of a plastome phylogeny. We generated 37 high quality plastome assemblies and performed a comparative analysis in terms of plastome structure, size, gene content and order, GC content, as well as number and type of repeat motifs. Using the obtained phylogenetic context, we reconstructed the evolution of these plastome attributes and assessed if significant shifts on the evolutionary traits' rates have occurred in the evolution of the core Tillandsioideae. Our results agree with previously published phylogenetic hypotheses based on plastid data, providing stronger statistical support for some recalcitrant nodes. However, phylogenetic discordance with previously published nuclear marker-based hypotheses was found. Several plastid markers that have been consistently used to address phylogenetic relationships within Tillandsioideae were highly informative for the retrieved plastome phylogeny and further loci are here identified as promising additional markers for future studies. New lineage-specific plastome rearrangements were found to support recently adopted taxonomic groups, including large inversions, as well as expansions and contractions of the inverted repeats. Evolutionary trait rate shifts associated with changes in size and GC content of the plastome regions were found across the phylogeny of core Tillandsioideae.

Comparative Analysis and Phylogenetic Relationships of Ceriops Species (Rhizophoraceae) and Avicennia lanata (Acanthaceae): Insight into the Chloroplast Genome Evolution between Middle and Seaward Zones of Mangrove Forests

Ceriops and Avicennia are true mangroves in the middle and seaward zones of mangrove forests, respectively. The chloroplast genomes of Ceriops decandra, Ceriops zippeliana, and Ceriops tagal were assembled into lengths of 166,650, 166,083 and 164,432 bp, respectively, whereas Avicennia lanata was 148,264 bp in length. The gene content and gene order are highly conserved among these species. The chloroplast genome contains 125 genes in A. lanata and 129 genes in Ceriops species. Three duplicate genes (rpl2, rpl23, and trnM-CAU) were found in the IR regions of the three Ceriops species, resulting in expansion of the IR regions. The rpl32 gene was lost in C. zippeliana, whereas the infA gene was present in A. lanata. Short repeats (<40 bp) and a lower number of SSRs were found in A. lanata but not in Ceriops species. The phylogenetic analysis supports that all Ceriops species are clustered in Rhizophoraceae and A. lanata is in Acanthaceae. In a search for genes under selective pressures of coastal environments, the rps7 gene was under positive selection compared with non-mangrove species. Finally, two specific primer sets were developed for species identification of the three Ceriops species. Thus, this finding provides insightful genetic information for evolutionary relationships and molecular markers in Ceriops and Avicennia species.

Characterization of the complete chloroplast genome of Brassica oleracea var. italica and phylogenetic relationships in Brassicaceae

Broccoli (Brassica oleracea var. italica) is an important B. oleracea cultivar, with high economic and agronomic value. However, comparative genome analyses are still needed to clarify variation among cultivars and phylogenetic relationships within the family Brassicaceae. Herein, the complete chloroplast (cp) genome of broccoli was generated by Illumina sequencing platform to provide basic information for genetic studies and to establish phylogenetic relationships within Brassicaceae. The whole genome was 153,364 bp, including two inverted repeat (IR) regions of 26,197 bp each, separated by a small single copy (SSC) region of 17,834 bp and a large single copy (LSC) region of 83,136 bp. The total GC content of the entire chloroplast genome accounts for 36%, while the GC content in each region of SSC,LSC, and IR accounts for 29.1%, 34.15% and 42.35%, respectively. The genome harbored 133 genes, including 88 protein-coding genes, 37 tRNAs, and 8 rRNAs, with 17 duplicates in IRs. The most abundant amino acid was leucine and the least abundant was cysteine. Codon usage analyses revealed a bias for A/T-ending codons. A total of 35 repeat sequences and 92 simple sequence repeats were detected, and the SC-IR boundary regions were variable between the seven cp genomes. A phylogenetic analysis suggested that broccoli is closely related to Brassica oleracea var. italica MH388764.1, Brassica oleracea var. italica MH388765.1, and Brassica oleracea NC_0441167.1. Our results are expected to be useful for further species identification, population genetics analyses, and biological research on broccoli.

Comparative Genomics, Phylogenetics, Biogeography, and Effects of Climate Change on Toddalia asiatica (L.) Lam. (Rutaceae) from Africa and Asia

In the present study, two samples of Toddalia asiatica species, both collected from Kenya, were sequenced and comparison of their genome structures carried out with T. asiatica species from China, available in the NCBI database. The genome size of both species from Africa was 158, 508 base pairs, which was slightly larger, compared to the reference genome of T. asiatica from Asia (158, 434 bp). The number of genes was 113 for both species from Africa, consisting of 79 protein-coding genes, 30 transfer RNA (tRNA) genes, and 4 ribosomal RNA (rRNA) genes. Toddalia asiatica from Asia had 115 genes with 81 protein-coding genes, 30 transfer RNA (tRNA) genes, and 4 ribosomal RNA (rRNA) genes. Both species compared displayed high similarity in gene arrangement. The gene number, orientation, and order were highly conserved. The IR/SC boundary structures were the same in all chloroplast genomes. A comparison of pairwise sequences indicated that the three regions (trnH-psbA, rpoB, and ycf1) were more divergent and can be useful in developing effective genetic markers. Phylogenetic analyses of the complete cp genomes and 79 protein-coding genes indicated that the Toddalia species collected from Africa were sister to T. asiatica collected from Asia. Both species formed a sister clade to the Southwest Pacific and East Asian species of Zanthoxylum. These results supported the previous studies of merging the genus Toddalia with Zanthoxylum and taxonomic change of Toddalia asiatica to Zanthoxylum asiaticum, which should also apply for the African species of Toddalia. Biogeographic results demonstrated that the two samples of Toddalia species from Africa diverged from T. asiatica from Asia (3.422 Mya, 95% HPD). These results supported an Asian origin of Toddalia species and later dispersal to Africa and Madagascar. The maxent model analysis showed that Asia would have an expansion of favorable areas for Toddalia species in the future. In Africa, there will be contraction and expansion of the favorable areas for the species. The availability of these cp genomes will provide valuable genetic resources for further population genetics and biogeographic studies of these species. However, more T. asiatica species collected from a wide geographical range are required.

Comparative Genomics, Phylogenetics, Biogeography, and Effects of Climate Change on Toddalia asiatica (L.) Lam. (Rutaceae) from Africa and Asia

In the present study, two samples of Toddalia asiatica species, both collected from Kenya, were sequenced and comparison of their genome structures carried out with T. asiatica species from China, available in the NCBI database. The genome size of both species from Africa was 158, 508 base pairs, which was slightly larger, compared to the reference genome of T. asiatica from Asia (158, 434 bp). The number of genes was 113 for both species from Africa, consisting of 79 protein-coding genes, 30 transfer RNA (tRNA) genes, and 4 ribosomal RNA (rRNA) genes. Toddalia asiatica from Asia had 115 genes with 81 protein-coding genes, 30 transfer RNA (tRNA) genes, and 4 ribosomal RNA (rRNA) genes. Both species compared displayed high similarity in gene arrangement. The gene number, orientation, and order were highly conserved. The IR/SC boundary structures were the same in all chloroplast genomes. A comparison of pairwise sequences indicated that the three regions (trnH-psbA, rpoB, and ycf1) were more divergent and can be useful in developing effective genetic markers. Phylogenetic analyses of the complete cp genomes and 79 protein-coding genes indicated that the Toddalia species collected from Africa were sister to T. asiatica collected from Asia. Both species formed a sister clade to the Southwest Pacific and East Asian species of Zanthoxylum. These results supported the previous studies of merging the genus Toddalia with Zanthoxylum and taxonomic change of Toddalia asiatica to Zanthoxylum asiaticum, which should also apply for the African species of Toddalia. Biogeographic results demonstrated that the two samples of Toddalia species from Africa diverged from T. asiatica from Asia (3.422 Mya, 95% HPD). These results supported an Asian origin of Toddalia species and later dispersal to Africa and Madagascar. The maxent model analysis showed that Asia would have an expansion of favorable areas for Toddalia species in the future. In Africa, there will be contraction and expansion of the favorable areas for the species. The availability of these cp genomes will provide valuable genetic resources for further population genetics and biogeographic studies of these species. However, more T. asiatica species collected from a wide geographical range are required.

Complete Chloroplast Genomes of Three Salix Species: Genome Structures and Phylogenetic Analysis

High genetic diversity and low differentiation present challenges in taxonomy and systematics of Salix. Chloroplast (cp) genome sequencing is efficient for providing new genomic information and elucidating phylogenetic relationships. Salix spathulifolia Seemen, S. cupularis Rehder, and S. annulifera C.Marquand & Airy Shaw are three shrubby willows spread in high-altitude regions in western China. In this study, the integrated circular cp genomes were sequenced and analyzed, and a phylogeny of Salix was constructed on the basis of the cp genomes. The results of chloroplast assembly and annotation information were used to characterize genome feature and interspecific variation. The phylogenetic position of the three willows was evaluated using phylogenetic analysis. Full-length cp genomes were 155,566–155,680 bp with a typical double-stranded circular quadripartite structure, containing one large single-copy region (LSC, 84,431–4552 bp), one small single-copy region (SSC: 16,206–16,221 bp), and two inverted repeats (IR: 27,453–27,461 bp). The cp genomes encoded 130 genes, including 8 rRNA genes, 37 tRNA genes, and 85 protein-coding genes. The guanine-cytosine (GC) content of the overall genome was 36.7%. Comparison among the three willows’ cp genomes revealed high similarity. Phylogenetic analysis indicated that S. spathulifolia was a basal taxon of clade I, while S. annulifera formed a monophyletic group with S. rorida Laksch.; S. cupularis was sister to S. suchowensis W.C. Cheng and S. psammophila Z. Wang & Chang Y. Yang. The complete chloroplast genomes of the three willows provides an additional sequence-based resource for studying the phylogeny and evolutionary history of Salicaceae.

Complete Chloroplast Genome Sequence of Erigeron breviscapus and Characterization of Chloroplast Regulatory Elements

Erigeron breviscapus is a famous medicinal plant. However, the limited chloroplast genome information of E. breviscapus, especially for the chloroplast DNA sequence resources, has hindered the study of E. breviscapus chloroplast genome transformation. Here, the complete chloroplast (cp) genome of E. breviscapus was reported. This genome was 152,164bp in length, included 37.2% GC content and was structurally arranged into two 24,699bp inverted repeats (IRs) and two single-copy areas. The sizes of the large single-copy region and the small single-copy region were 84,657 and 18,109bp, respectively. The E. breviscapus cp genome consisted of 127 coding genes, including 83 protein coding genes, 36 transfer RNA (tRNA) genes, and eight ribosomal RNA (rRNA) genes. For those genes, 95 genes were single copy genes and 16 genes were duplicated in two inverted regions with seven tRNAs, four rRNAs, and five protein coding genes. Then, genomic DNA of E. breviscapus was used as a template, and the endogenous 5' and 3' flanking sequences of the trnI gene and trnA gene were selected as homologous recombinant fragments in vector construction and cloned through PCR. The endogenous 5' flanking sequences of the psbA gene and rrn16S gene, the endogenous 3' flanking sequences of the psbA gene, rbcL gene, and rps16 gene and one sequence element from the psbN-psbH chloroplast operon were cloned, and certain chloroplast regulatory elements were identified. Two homologous recombination fragments and all of these elements were constructed into the cloning vector pBluescript SK (+) to yield a series of chloroplast expression vectors, which harbored the reporter gene EGFP and the selectable marker aadA gene. After identification, the chloroplast expression vectors were transformed into Escherichia coli and the function of predicted regulatory elements was confirmed by a spectinomycin resistance test and fluorescence intensity measurement. The results indicated that aadA gene and EGFP gene were efficiently expressed under the regulation of predicted regulatory elements and the chloroplast expression vector had been successfully constructed, thereby providing a solid foundation for establishing subsequent E. breviscapus chloroplast transformation system and genetic improvement of E. breviscapus.

Complete chloroplast genome sequences of five Bruguiera species (Rhizophoraceae): comparative analysis and phylogenetic relationships

Bruguiera is a genus of true mangroves that are mostly distributed in the Indo-West Pacific region. However, the number of published whole chloroplast genome sequences of Bruguiera species are limited. Here, the complete chloroplast sequences of five Bruguiera species were sequenced and assembled using Illumina data. The chloroplast genomes of B. gymnorhiza, B. hainesii, B. cylindrica, B. parviflora and B. sexangula were assembled into 161,195, 164,295, 164,297, 163,228 and 164,170 bp, respectively. All chloroplast genomes contain 37 tRNA and eight rRNA genes, with either 84 or 85 protein-coding genes. A comparative analysis of these genomes revealed high similarity in gene structure, gene order and boundary position of the LSC, SSC and two IR regions. Interestingly, B. gymnorhiza lost a rpl32 gene in the SSC region. In addition, a ndhF gene in B. parviflora straddles both the SSC and IRB boundary regions. These genes reveal differences in chloroplast evolution among Bruguiera species. Repeats and SSRs in the chloroplast genome sequences were found to be highly conserved between B. cylindrica and B. hainesii as well as B. gymnorhiza and B. sexangula indicating close genetic relationships based on maternal inheritance. Notably, B. hainesii, which is considered a hybrid between B. gymnorhiza and B. cylindrica, appears to have inherited the chloroplast from B. cylindrica. Investigating the effects of selection events on shared protein-coding genes showed a positive selection in rps7 and rpl36 genes in all species compared to land-plant species. A phylogenetic analysis, based on 59 conserved chloroplast protein-coding genes, showed strong support that all Bruguiera species are in the clade Rhizophoraceae. This study provides valuable genetic information for the study of evolutionary relationships and population genetics in Bruguiera and other mangrove species.

Plastome Characterization and Phylogenomic Analysis Yield New Insights into the Evolutionary Relationships among the Species of the Subgenus Bryocles (Hosta; Asparagaceae) in East Asia

The genus Hosta, which has a native distribution in temperate East Asia and a number of species ranging from 23 to 40, represents a taxonomically important and ornamentally popular plant. Despite its taxonomic and horticultural importance, the genus Hosta has remained taxonomically challenging owing to insufficient diagnostic features, continuous morphological variation, and the process of hybridization and introgression, making species circumscription and phylogenetic inference difficult. In this study, we sequenced 11 accessions of Hosta plastomes, including members of three geographically defined subgenera, Hosta, Bryocles, and Giboshi, determined the characteristics of plastomes, and inferred their phylogenetic relationships. We found highly conserved plastomes among the three subgenera, identified several mutation hotspots that can be used as barcodes, and revealed the patterns of codon usage bias and RNA editing sites. Five positively selected plastome genes (rbcL, rpoB, rpoC2, rpl16, and rpl20) were identified. Phylogenetic analysis suggested (1) the earliest divergence of subg. Hosta, (2) non-monophyly of subg. Bryocles and its two sections (Lamellatae and Stoloniferae), (3) a sister relationship between H. sieboldiana (subg. Giboshi) and H. ventricosa (subg. Bryocles), and (4) reciprocally monophyletic and divergent lineages of H. capitata in Korea and Japan, requiring further studies of their taxonomic distinction.

Unprecedented Intraindividual Structural Heteroplasmy in Eleocharis (Cyperaceae, Poales) Plastomes

Plastid genomes (plastomes) of land plants have a conserved quadripartite structure in a gene-dense unit genome consisting of a large inverted repeat that separates two single copy regions. Recently, alternative plastome structures were suggested in Geraniaceae and in some conifers and Medicago the coexistence of inversion isomers has been noted. In this study, plastome sequences of two Cyperaceae, Eleocharis dulcis (water chestnut) and Eleocharis cellulosa (gulf coast spikerush), were completed. Unlike the conserved plastomes in basal groups of Poales, these Eleocharis plastomes have remarkably divergent features, including large plastome sizes, high rates of sequence rearrangements, low GC content and gene density, gene duplications and losses, and increased repetitive DNA sequences. A novel finding among these features was the unprecedented level of heteroplasmy with the presence of multiple plastome structural types within a single individual. Illumina paired-end assemblies combined with PacBio single-molecule real-time sequencing, long-range polymerase chain reaction, and Sanger sequencing data identified at least four different plastome structural types in both Eleocharis species. PacBio long read data suggested that one of the four E. dulcis plastome types predominates.

Complete Chloroplast Genomes of Acanthochlamys bracteata (China) and Xerophyta (Africa) (Velloziaceae): Comparative Genomics and Phylogenomic Placement

Acanthochlamys P.C. Kao is a Chinese endemic monotypic genus, whereas Xerophyta Juss. is a genus endemic to Africa mainland, Arabian Peninsula and Madagascar with ca.70 species. In this recent study, the complete chloroplast genome of Acanthochlamys bracteata was sequenced and its genome structure compared with two African Xerophyta species (Xerophyta spekei and Xerophyta viscosa) present in the NCBI database. The genomes showed a quadripartite structure with their sizes ranging from 153,843 bp to 155,498 bp, having large single-copy (LSC) and small single-copy (SSC) regions divided by a pair of inverted repeats (IR regions). The total number of genes found in A. bracteata, X. spekei and X. viscosa cp genomes are 129, 130, and 132, respectively. About 50, 29, 28 palindromic, forward and reverse repeats and 90, 59, 53 simple sequence repeats (SSRs) were found in the A. bracteata, X. spekei, and X. viscosa cp genome, respectively. Nucleotide diversity analysis in all species was 0.03501, Ka/Ks ratio average score was calculated to be 0.26, and intergeneric K2P value within the Order Pandanales was averaged to be 0.0831. Genomic characterization was undertaken by comparing the genomes of the three species of Velloziaceae and it revealed that the coding regions were more conserved than the non-coding regions. However, key variations were noted mostly at the junctions of IRs/SSC regions. Phylogenetic analysis suggests that A. bracteata species has a closer genetic relationship to the genus Xerophyta. The present study reveals the complete chloroplast genome of A. bracteata and gives a genomic comparative analysis with the African species of Xerophyta. Thus, can be useful in developing DNA markers for use in the study of genetic variabilities and evolutionary studies in Velloziaceae.

Chloroplast genome assembly of Handroanthus impetiginosus: comparative analysis and molecular evolution in Bignoniaceae

Bignoniaceae species have conserved chloroplast structure, with hotspots of nucleotide diversity. Several genes are under positive selection, and can be targets for evolutionary studies. Bignoniaceae is one of the most species-rich family of woody plants in Neotropical seasonally dry forests. Here we report the assembly of Handroanthus impetiginosus chloroplast genome and evolutionary comparative analyses of ten Bignoniaceae species representing the genera for which whole-genome chloroplast sequences were available. The chloroplast genome of H. impetiginosus is 159,462 bp in size and has a similar structure compared to the other nine species. The total number of genes was slightly variable amongst the Bignoniaceae, ranging from 124 in H. impetiginosus to 144 in Anemopaegma acutifolium. The inverted repeat (IR) size was variable, ranging from 24,657 bp (Tecomaria capensis) to 40,481 bp (A. acutifolium), due to the contraction and retraction at its boundaries. However, gene boundaries were very similar among the ten species. We found 98 forward and palindromic dispersed repeats, and 85 simple sequence repeats (SSRs). In general, chloroplast sequences were highly conserved, with few nucleotide diversity hotspots in the genes accD, clpP, rpoA, ycf1, ycf2. The phylogenetic analysis based on 77 coding genes was highly consistent with Angiosperm Phylogeny Group (APG) IV. Our results also indicate that most genes are under negative selection or neutral evolution. We found no evidence of branch-site selection, implying that H. impetiginosus is not evolving faster than the other species analyzed, notwithstanding we found site positive selection signal in several genes. These genes can provide targets for evolutionary studies in Bignoniaceae and Lamiales species.

Complete Chloroplast Genome of Argania spinosa: Structural Organization and Phylogenetic Relationships in Sapotaceae

Argania spinosa (Sapotaceae), an important endemic Moroccan oil tree, is a primary source of argan oil, which has numerous dietary and medicinal proprieties. The plant species occupies the mid-western part of Morocco and provides great environmental and socioeconomic benefits. The complete chloroplast (cp) genome of A. spinosa was sequenced, assembled, and analyzed in comparison with those of two Sapotaceae members. The A. spinosa cp genome is 158,848 bp long, with an average GC content of 36.8%. The cp genome exhibits a typical quadripartite and circular structure consisting of a pair of inverted regions (IR) of 25,945 bp in length separating small single-copy (SSC) and large single-copy (LSC) regions of 18,591 and 88,367 bp, respectively. The annotation of A. spinosa cp genome predicted 130 genes, including 85 protein-coding genes (CDS), 8 ribosomal RNA (rRNA) genes, and 37 transfer RNA (tRNA) genes. A total of 44 long repeats and 88 simple sequence repeats (SSR) divided into mononucleotides (76), dinucleotides (7), trinucleotides (3), tetranucleotides (1), and hexanucleotides (1) were identified in the A. spinosa cp genome. Phylogenetic analyses using the maximum likelihood (ML) method were performed based on 69 protein-coding genes from 11 species of Ericales. The results confirmed the close position of A. spinosa to the Sideroxylon genus, supporting the revisiting of its taxonomic status. The complete chloroplast genome sequence will be valuable for further studies on the conservation and breeding of this medicinally and culinary important species and also contribute to clarifying the phylogenetic position of the species within Sapotaceae.

Complete chloroplast genomes of four Physalis species (Solanaceae): lights into genome structure, comparative analysis, and phylogenetic relationships

BACKGROUND

Physalis L. is a genus of herbaceous plants of the family Solanaceae, which has important medicinal, edible, and ornamental values. The morphological characteristics of Physalis species are similar, and it is difficult to rapidly and accurately distinguish them based only on morphological characteristics. At present, the species classification and phylogeny of Physalis are still controversial. In this study, the complete chloroplast (cp) genomes of four Physalis species (Physalis angulata, P. alkekengi var. franchetii, P. minima and P. pubescens) were sequenced, and the first comprehensive cp genome analysis of Physalis was performed, which included the previously published cp genome sequence of Physalis peruviana.

RESULTS

The Physalis cp genomes exhibited typical quadripartite and circular structures, and were relatively conserved in their structure and gene synteny. However, the Physalis cp genomes showed obvious variations at four regional boundaries, especially those of the inverted repeat and the large single-copy regions. The cp genomes' lengths ranged from 156,578 bp to 157,007 bp. A total of 114 different genes, 80 protein-coding genes, 30 tRNA genes, and 4 rRNA genes, were observed in four new sequenced Physalis cp genomes. Differences in repeat sequences and simple sequence repeats were detected among the Physalis cp genomes. Phylogenetic relationships among 36 species of 11 genera of Solanaceae based on their cp genomes placed Physalis in the middle and upper part of the phylogenetic tree, with a monophyletic evolution having a 100% bootstrap value.

CONCLUSION

Our results enrich the data on the cp genomes of the genus Physalis. The availability of these cp genomes will provide abundant information for further species identification, increase the taxonomic and phylogenetic resolution of Physalis, and assist in the investigation and utilization of Physalis plants.

The Chloroplast Genome of Carya illinoinensis: Genome Structure, Adaptive Evolution, and Phylogenetic Analysis

Research Highlights: For the first time, the complete chloroplast (cp) genome of Carya illinoinensis cv. ‘Pawnee’ was de novo assembled. Comprehensive analysis the cp genome of C. illinoinensis revealed potential cpDNA markers for intraspecies identification, genes involved in adaptation, and its phylogenetic position. Background and Objectives: C. illinoinensis is an economically important nut tree in the family Juglandaceae. Cp-derived markers are helpful for genetic research, but they still need to be developed in C. illinoinensis. Additionally, the adaptation and phylogenetic relationships of C. illinoinensis have not been revealed based on the complete cp genome. Materials and Methods: Chloroplast genomic DNA of C. illinoinensis cv. ‘Pawnee’ was extracted and subjected to Illumina sequencing. Results: The cp genome is 160,819 bp in size, exhibiting a typical quadripartite structure with a large single copy (LSC) of 90,022 bp, a small single copy (SSC) of 18,791 bp, and a pair of inverted repeats (IRA and IRB) regions of 26,003 bp each. The genome was predicted to encode 112 unique genes, including 79 protein-coding genes, 29 tRNAs, and four rRNAs, with 19 duplicates in the IR regions. In total, 213 SSRs and 44 long repeats were identified in the cp genome. A comparison of two different C. illinoinensis genotypes, ‘Pawnee’ and 87MX3-2.11, obtained 143 SNPs and 74 indels. The highly variable regions such as atpF, clpP, and ndhA genes, and matK-rps16, trnS-trnG, and trnT-psbD intergenic spacers might be helpful for future intraspecific identification. Positive selection was acting on the ccsA and rps12 cp genes based on the Ka/Ks ratios. Phylogenetic analysis indicated that C. illinoinensis forms a sister clade to Asian Carya species, represented by C. kweichowensis and Annamocarya sinensis. Conclusions: The genome information in our study will have significance for further research on the intraspecies identification and genetic improvement of C. illinoinensis.

Comparative Analysis of Actaea Chloroplast Genomes and Molecular Marker Development for the Identification of Authentic Cimicifugae Rhizoma

Actaea (Ranunculaceae; syn. Cimicifuga) is a controversial and complex genus. Dried rhizomes of Actaea species are used as Korean traditional herbal medicine. Although Actaea species are valuable, given their taxonomic classification and medicinal properties, sequence information of Actaea species is limited. In this study, we determined the complete chloroplast (cp) genome sequences of three Actaea species, including A. simplex, A. dahurica, and A. biternata. The cp genomes of these species varied in length from 159,523 to 159,789 bp and contained 112 unique functional genes, including 78 protein-coding genes, 30 transfer RNA genes, and 4 ribosomal RNA genes. Gene order, orientation, and content were well conserved in the three cp genomes. Comparative sequence analysis revealed the presence of hotspots, including ndhC-trnV-UAC, in Actaea cp genomes. High-resolution phylogenetic relationships were established among Actaea species based on cp genome sequences. Actaea species were clustered into each Actaea section, consistent with the Angiosperm Phylogeny Group (APG) IV system of classification. We also developed a novel indel marker, based on copy number variation of tandem repeats, to facilitate the authentication of the herbal medicine Cimicifugae Rhizoma. The availability Actaea cp genomes will provide abundant information for the taxonomic and phylogenetic analyses of Actaea species, and the Actaea (ACT) indel marker will be useful for the authentication of the herbal medicine.

The complete chloroplast genome sequence of yellow mustard (Sinapis alba L.) and its phylogenetic relationship to other Brassicaceae species

As a member of the large Brassicaceae family, yellow mustard (Sinapis alba L.) has been used as an important gene pool for the genetic improvement of cash crops in Brassicaceae. Understanding the phylogenetic relationship between Sinapis alba (S. alba) and other Brassicaceae crops can provide guidance on the introgression of its favorable alleles into related species. The chloroplast (cp) genome is an ideal model for assessing genome evolution and the phylogenetic relationships of complex angiosperm families. Herein, we de novo assembled the complete cp genome of S. alba by integrating the PacBio and Illumina sequencing platforms. A 153,760 bp quadripartite cycle without any gap was obtained, including a pair of inverted repeats (IRa and IRb) of 26,221 bp, separated by a large single copy (LSC) region of 83,506 bp and a small single copy (SSC) region of 17,821 bp. A total of 78 protein-coding genes, 30 tRNA genes, and four rRNA genes were identified in this cp genome, as were 89 simple sequence repeat (SSR) loci of 18 types. The codon usage analysis revealed a preferential use of the Leu codon with the A/U ending. The phylogenetic analysis using 82 Brassicaceae species demonstrated that S. alba had a close relationship with important Brassica and Raphanus species; moreover, it likely originated from a separate evolutionary pathway compared with the congeneric Sinapis arvensis. The synonymous (Ks) and non-synonymous (Ks) substitution rate analysis showed that genes encoding "Subunits of cytochrome b/f complex" were under the lowest purifying selection pressure, whereas those associated with "Maturase", "Subunit of acetyl-CoA", and "Subunits of NADH-dehydrogenase" underwent relatively higher purifying selection pressures. Our results provide valuable information for fully utilizing the S. alba cp genome as a potential genetic resource for the genetic improvement of Brassica and Raphanus species.

Methods and Tools for Plant Organelle Genome Sequencing, Assembly, and Downstream Analysis

Organelles play an important role in a eukaryotic cell. Among them, the two organelles, chloroplast and mitochondria, are responsible for the critical function of photosynthesis and aerobic respiration. Organellar genomes are also very important for plant systematic studies. Here we have described the methods for isolation of the mitochondrial and plastid DNA and its subsequent sequencing with the help of NGS technology. We have also discussed in detail the various tools available for assembly, annotation, and visualization of the organelle genome sequence.

The complete chloroplast genome of Stryphnodendron adstringens (Leguminosae - Caesalpinioideae): comparative analysis with related Mimosoid species

Stryphnodendron adstringens is a medicinal plant belonging to the Leguminosae family, and it is commonly found in the southeastern savannas, endemic to the Cerrado biome. The goal of this study was to assemble and annotate the chloroplast genome of S. adstringens and to compare it with previously known genomes of the mimosoid clade within Leguminosae. The chloroplast genome was reconstructed using de novo and referenced-based assembly of paired-end reads generated by shotgun sequencing of total genomic DNA. The size of the S. adstringens chloroplast genome was 162,169 bp. This genome included a large single-copy (LSC) region of 91,045 bp, a small single-copy (SSC) region of 19,014 bp and a pair of inverted repeats (IRa and IRb) of 26,055 bp each. The S. adstringens chloroplast genome contains a total of 111 functional genes, including 77 protein-coding genes, 30 transfer RNA genes, and 4 ribosomal RNA genes. A total of 137 SSRs and 42 repeat structures were identified in S. adstringens chloroplast genome, with the highest proportion in the LSC region. A comparison of the S. adstringens chloroplast genome with those from other mimosoid species indicated that gene content and synteny are highly conserved in the clade. The phylogenetic reconstruction using 73 conserved coding-protein genes from 19 Leguminosae species was supported to be paraphyletic. Furthermore, the noncoding and coding regions with high nucleotide diversity may supply valuable markers for molecular evolutionary and phylogenetic studies at different taxonomic levels in this group.

The Complete Chloroplast Genomes of Punica granatum and a Comparison with Other Species in Lythraceae

Pomegranates (Punica granatum L.) are one of the most popular fruit trees cultivated in arid and semi-arid tropics and subtropics. In this study, we determined and characterized three complete chloroplast (cp) genomes of P. granatum cultivars with different phenotypes using the genome skimming approach. The complete cp genomes of three pomegranate cultivars displayed the typical quadripartite structure of angiosperms, and their length ranged from 156,638 to 156,639 bp. They encoded 113 unique genes and 17 are duplicated in the inverted regions. We analyzed the sequence diversity of pomegranate cp genomes coupled with two previous reports. The results showed that the sequence diversity is extremely low and no informative sites were detected, which suggests that cp genome sequences may be not be suitable for investigating the genetic diversity of pomegranate genotypes. Further, we analyzed the codon usage pattern and identified the potential RNA editing sites. A comparative cp genome analysis with other species within Lythraceae revealed that the gene content and organization are highly conserved. Based on a site-specific model, 11 genes with positively selected sites were detected, and most of them were photosynthesis-related genes and genetic system-related genes. Together with previously released cp genomes of the order Myrtales, we determined the taxonomic position of P. granatum based on the complete chloroplast genomes. Phylogenetic analysis suggested that P. granatum form a single clade with other species from Lythraceae with a high support value. The complete cp genomes provides valuable information for understanding the phylogenetic position of P. gramatum in the order Myrtales.

The complete chloroplast genome sequence of watercress (Nasturtium officinale R. Br.): Genome organization, adaptive evolution and phylogenetic relationships in Cardamineae

Watercress (Nasturtium officinale R. Br.), an aquatic leafy vegetable of the Brassicaceae family, is known as a nutritional powerhouse. Here, we de novo sequenced and assembled the complete chloroplast (cp) genome of watercress based on combined PacBio and Illumina data. The cp genome is 155,106 bp in length, exhibiting a typical quadripartite structure including a pair of inverted repeats (IRA and IRB) of 26,505 bp separated by a large single copy (LSC) region of 84,265 bp and a small single copy (SSC) region of 17,831 bp. The genome contained 113 unique genes, including 79 protein-coding genes, 30 tRNAs and 4 rRNAs, with 20 duplicate in the IRs. Compared with the prior cp genome of watercress deposited in GenBank, 21 single nucleotide polymorphisms (SNPs) and 27 indels were identified, mainly located in noncoding sequences. A total of 49 repeat structures and 71 simple sequence repeats (SSRs) were detected. Codon usage showed a bias for A/T-ending codons in the cp genome of watercress. Moreover, 45 RNA editing sites were predicted in 16 genes, all for C-to-U transitions. A comparative plastome study with Cardamineae species revealed a conserved gene order and high similarity of protein-coding sequences. Analysis of the Ka/Ks ratios of Cardamineae suggested positive selection exerted on the ycf2 gene in watercress, which might reflect specific adaptations of watercress to its particular living environment. Phylogenetic analyses based on complete cp genomes and common protein-coding genes from 56 species showed that the genus Nasturtium was a sister to Cardamine in the Cardamineae tribe. Our study provides valuable resources for future evolution, population genetics and molecular biology studies of watercress.

Comparative chloroplast genomes of Paris Sect. Marmorata: insights into repeat regions and evolutionary implications

Background

Species of Paris Sect. Marmorata are valuable medicinal plants to synthesize steroidal saponins with effective pharmacological therapy. However, the wild resources of the species are threatened by plundering exploitation before the molecular genetics studies uncover the genomes and evolutionary significance. Thus, the availability of complete chloroplast genome sequences of Sect. Marmorata is necessary and crucial to the understanding the plastome evolution of this section and facilitating future population genetics studies. Here, we determined chloroplast genomes of Sect. Marmorata, and conducted the whole chloroplast genome comparison.

Results

This study presented detailed sequences and structural variations of chloroplast genomes of Sect. Marmorata. Over 40 large repeats and approximately 130 simple sequence repeats as well as a group of genomic hotspots were detected. Inverted repeat contraction of this section was inferred via comparing the chloroplast genomes with the one of P. verticillata. Additionally, almost all the plastid protein coding genes were found to prefer ending with A/U. Mutation bias and selection pressure predominately shaped the codon bias of most genes. And most of the genes underwent purifying selection, whereas photosynthetic genes experienced a relatively relaxed purifying selection.

Conclusions

Repeat sequences and hotspot regions can be scanned to detect the intraspecific and interspecific variability, and selected to infer the phylogenetic relationships of Sect. Marmorata and other species in subgenus Daiswa. Mutation and natural selection were the main forces to drive the codon bias pattern of most plastid protein coding genes. Therefore, this study enhances the understanding about evolution of Sect. Marmorata from the chloroplast genome, and provide genomic insights into genetic analyses of Sect. Marmorata.

Electronic supplementary material

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

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.

Comparative analysis of the complete plastomes of garlic Allium sativum and bulb onion Allium cepa

Sequencing and comparative characterization of plant plastid genomes, or plastomes, is an important tool for modern phylogenetic and taxonomic studies, as well as for understanding the plastome evolution. The genusAlliumL. (family Amaryllidaceae) incorporates more than 900 species, includes economically signifi­cant vegetable crops such as garlicA. sativum, onionA. cepa, leekA. porrum, etc. In this work, the plastome of garlicA. sativumhas been completely sequenced. TheA. sativumplastome is 153172 bp in size. It consists of a large unique (LSC, 82035 bp) and small unique (SSC, 18015 bp) copies, separated by inverted repeats (IRa and IRb) of 26561 bp each. In the garlic plastome, 134 genes have been annotated: 82 protein-coding genes, 38 tRNA genes, 8 rRNA genes, and 6 pseudogenes. Comparative analysis ofA. sativumandA. cepaplastomes reveals differences in the sizes of structural elements and spacers at the inverted repeat bound­aries. The total numbers of genes inA. sativumandA. cepaare the same, but the gene composition is dif­ferent: therpl22gene is functional inA. sativum, being a pseudogene inA. cepa; conversely, therps16gene is a pseudogene inA. sativumand a protein-coding gene inA. cepa. In theA. sativumandA. cepaplastomes, 32 SSR sequences have been identified. More than half of them are dinucleotides, and the remaining are tetra-, penta-, and hexanucleotides at the same time, trinucleotides were absent. The compared plastomes differ in the numbers of certain SSRs, and some are present in only one of the species.

Authentication of Herbal Medicines Dipsacus asper and Phlomoides umbrosa Using DNA Barcodes, Chloroplast Genome, and Sequence Characterized Amplified Region (SCAR) Marker

Dried roots of Dipsacus asper (Caprifoliaceae) are used as important traditional herbal medicines in Korea. However, the roots are often used as a mixture or contaminated with Dipsacus japonicus in Korean herbal markets. Furthermore, the dried roots of Phlomoides umbrosa (Lamiaceae) are used indiscriminately with those of D. asper, with the confusing Korean names of Sok-Dan and Han-Sok-Dan for D. asper and P. umbrosa, respectively. Although D. asper and P. umbrosa are important herbal medicines, the molecular marker and genomic information available for these species are limited. In this study, we analysed DNA barcodes to distinguish among D. asper, D. japonicus, and P. umbrosa and sequenced the chloroplast (CP) genomes of D. asper and D. japonicus. The CP genomes of D. asper and D. japonicus were 160,530 and 160,371 bp in length, respectively, and were highly divergent from those of the other Caprifoliaceae species. Phylogenetic analysis revealed a monophyletic group within Caprifoliaceae. We also developed a novel sequence characterised amplified region (SCAR) markers to distinguish among D. asper, D. japonicus, and P. umbrosa. Our results provide important taxonomic, phylogenetic, and evolutionary information on the Dipsacus species. The SCAR markers developed here will be useful for the authentication of herbal medicines.

Study of spontaneous mutations in the transmission of poplar chloroplast genomes from mother to offspring

Background

Chloroplasts have their own genomes, independent from nuclear genomes, that play vital roles in growth, which is a major targeted trait for genetic improvement in Populus. Angiosperm chloroplast genomes are maternally inherited, but the chloroplast’ variation pattern of poplar at the single-base level during the transmission from mother to offspring remains unknown.

Results

Here, we constructed high-quality and almost complete chloroplast genomes for three poplar clones, ‘NL895’ and its parents, ‘I69’ and ‘I45’, from the short-read datasets using multi-pass sequencing (15–16 times per clone) and ultra-high coverage (at least 8500× per clone), with the four-step strategy of Simulation–Assembly–Merging–Correction. Each of the three resulting chloroplast assemblies contained contigs covering > 99% of Populus trichocarpa chloroplast DNA as a reference. A total of 401 variant loci were identified by a hybrid strategy of genome comparison-based and mapping-based single nucleotide polymorphism calling. The genotypes of 94 variant loci were different among the three poplar clones. However, only 1 of the 94 loci was a missense mutation, which was located in the exon region of rpoC1 encoding the β’ subunit of plastid-encoded RNA polymerase. The genotype of the loci in NL895 and its female parent (I69) was different from that of its male parent (I45).

Conclusions

This research provides resources for further chloroplast genomic studies of a F₁ full-sibling family derived from a cross between I69 and I45, and will improve the application of chloroplast genomic information in modern Populus breeding programs.

Electronic supplementary material

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

The Complete Chloroplast Genome of Heimia myrtifolia and Comparative Analysis within Myrtales

Heimia myrtifolia is an important medicinal plant with several pharmacologically active alkaloids and is also used as an ornamental landscape plant. The purpose of this study is to complete and characterize the chloroplast (cp) genome of H. myrtifolia and compare genomic features to other Myrtales species' cp genomes. The analysis showed that H. myrtifolia has a total length of 159,219 bp with a typical quadripartite structure containing two identical inverted repeats (IRs) of 25,643 bp isolated by one large single copy (LSC) of 88,571 bp and one small single copy (SSC) of 18,822 bp. The H. myrtifolia cp genome contains 129 genes with eight ribosomal RNAs, 30 transfer RNAs, and 78 protein coding genes, in which 17 genes are duplicated in two IR regions. The genome organization including gene type and number and guanine-cytosine (GC) content is analyzed among the 12 cp genomes in this study. Approximately 255 simple sequence repeats (SSRs) and 16 forward, two reverses, and two palindromic repeats were identified in the H. myrtifolia cp genome. By comparing the whole H. myrtifolia cp genome with 11 other Myrtales species, the results showed that the sequence similarity was high between coding regions while sequence divergence was high between intergenic regions. By employing the full cp genomes for phylogenetic analysis, structural and sequence differences were characterized between H. myrtifolia and 11 Myrtales species illustrating what patterns are common in the evolution of cp genomes within the Myrtales. The first entire cp genome in the genus Heimia provides a valuable resource for further studies in these medicinally and ornamentally important taxa.

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.

The Complete Chloroplast Genome of Catha edulis: A Comparative Analysis of Genome Features with Related Species

Qat (Catha edulis, Celastraceae) is a woody evergreen species with great economic and cultural importance. It is cultivated for its stimulant alkaloids cathine and cathinone in East Africa and southwest Arabia. However, genome information, especially DNA sequence resources, for C. edulis are limited, hindering studies regarding interspecific and intraspecific relationships. Herein, the complete chloroplast (cp) genome of Catha edulis is reported. This genome is 157,960 bp in length with 37% GC content and is structurally arranged into two 26,577 bp inverted repeats and two single-copy areas. The size of the small single-copy and the large single-copy regions were 18,491 bp and 86,315 bp, respectively. The C. edulis cp genome consists of 129 coding genes including 37 transfer RNA (tRNA) genes, 8 ribosomal RNA (rRNA) genes, and 84 protein coding genes. For those genes, 112 are single copy genes and 17 genes are duplicated in two inverted regions with seven tRNAs, four rRNAs, and six protein coding genes. The phylogenetic relationships resolved from the cp genome of qat and 32 other species confirms the monophyly of Celastraceae. The cp genomes of C. edulis, Euonymus japonicus and seven Celastraceae species lack the rps16 intron, which indicates an intron loss took place among an ancestor of this family. The cp genome of C. edulis provides a highly valuable genetic resource for further phylogenomic research, barcoding and cp transformation in Celastraceae.

Chloroplast genomes of Byrsonima species (Malpighiaceae): comparative analysis and screening of high divergence sequences

Byrsonima is the third largest genus (about 200 species) in the Malpighiaceae family, and one of the most common in Brazilian savannas. However, there is no molecular phylogeny available for the genus and taxonomic uncertainties at the generic and family level still remain. Herein, we sequenced the complete chloroplast genome of B. coccolobifolia and B. crassifolia, the first ones described for Malpighiaceae, and performed comparative analyses with sequences previously published for other families in the order Malpighiales. The chloroplast genomes assembled had a similar structure, gene content and organization, even when compared with species from other families. Chloroplast genomes ranged between 160,212 bp in B. crassifolia and 160,329 bp in B. coccolobifolia, both containing 115 genes (four ribosomal RNA genes, 28 tRNA genes and 83 protein-coding genes). We also identified sequences with high divergence that might be informative for phylogenetic inferences in the Malpighiales order, Malpighiaceae family and within the genus Byrsonima. The phylogenetic reconstruction of Malpighiales with these regions highlighted their utility for phylogenetic studies. The comparative analyses among species in Malpighiales provided insights into the chloroplast genome evolution in this order, including the presence/absence of three genes (infA, rpl32 and rps16) and two pseudogenes (ycf1 and rps19).

The Complete Chloroplast Genome Sequences of Aconitum pseudolaeve and Aconitum longecassidatum, and Development of Molecular Markers for Distinguishing Species in the Aconitum Subgenus Lycoctonum

Aconitum pseudolaeve Nakai and Aconitum longecassidatum Nakai, which belong to the Aconitum subgenus Lycoctonum, are distributed in East Asia and Korea. Aconitum species are used in herbal medicine and contain highly toxic components, including aconitine. A. pseudolaeve, an endemic species of Korea, is a commercially valuable material that has been used in the manufacture of cosmetics and perfumes. Although Aconitum species are important plant resources, they have not been extensively studied, and genomic information is limited. Within the subgenus Lycoctonum, which includes A. pseudolaeve and A. longecassidatum, a complete chloroplast (CP) genome is available for only one species, Aconitum barbatum Patrin ex Pers. Therefore, we sequenced the complete CP genomes of two Aconitum species, A. pseudolaeve and A. longecassidatum, which are 155,628 and 155,524 bp in length, respectively. Both genomes have a quadripartite structure consisting of a pair of inverted repeated regions (51,854 and 52,108 bp, respectively) separated by large single-copy (86,683 and 86,466 bp) and small single-copy (17,091 and 16,950 bp) regions similar to those in other Aconitum CP genomes. Both CP genomes consist of 112 unique genes, 78 protein-coding genes, 4 ribosomal RNA (rRNA) genes, and 30 transfer RNA (tRNA) genes. We identified 268 and 277 simple sequence repeats (SSRs) in A. pseudolaeve and A. longecassidatum, respectively. We also identified potential 36 species-specific SSRs, 53 indels, and 62 single-nucleotide polymorphisms (SNPs) between the two CP genomes. Furthermore, a comparison of the three Aconitum CP genomes from the subgenus Lycoctonum revealed highly divergent regions, including trnK-trnQ, ycf1-ndhF, and ycf4-cemA. Based on this finding, we developed indel markers using indel sequences in trnK-trnQ and ycf1-ndhF. A. pseudolaeve, A. longecassidatum, and A. barbatum could be clearly distinguished using the novel indel markers AcoTT (Aconitum trnK-trnQ) and AcoYN (Aconitum ycf1-ndhF). These two new complete CP genomes provide useful genomic information for species identification and evolutionary studies of the Aconitum subgenus Lycoctonum.

The complete chloroplast genome sequence of the CAM epiphyte Spanish moss (Tillandsia usneoides, Bromeliaceae) and its comparative analysis

Spanish moss (Tillandsia usneoides) is an epiphytic bromeliad widely distributed throughout tropical and warm temperate America. This plant is highly adapted to extreme environmental conditions. Striking features of this species include specialized trichomes (scales) covering the surface of its shoots aiding the absorption of water and nutrients directly from the atmosphere and a specific photosynthesis using crassulacean acid metabolism (CAM). Here we report the plastid genome of Spanish moss and present the comparison of genome organization and sequence evolution within Poales. The plastome of Spanish moss has a quadripartite structure consisting of a large single copy (LSC, 87,439 bp), two inverted regions (IRa and IRb, 26,803 bp) and short single copy (SSC, 18,612 bp) region. The plastid genome had 37.2% GC content and 134 genes with 88 being unique protein-coding genes and 20 of these are duplicated in the IR, similar to other reported bromeliads. Our study shows that early diverging lineages of Poales do not have high substitution rates as compared to grasses, and plastid genomes of bromeliads show structural features considered to be ancestral in graminids. These include the loss of the introns in the clpP and rpoC1 genes and the complete loss or partial degradation of accD and ycf genes in the Graminid clade. Further structural rearrangements appeared in the graminids lacking in Spanish moss, which include a 28-kb inversion between the trnG-UCC-rps14 region and 6-kb in the trnG-UCC-psbD, followed by a third <1kb inversion in the trnT sequence.

Complete Chloroplast Genome of Pinus massoniana (Pinaceae): Gene Rearrangements, Loss of ndh Genes, and Short Inverted Repeats Contraction, Expansion

The chloroplast genome (CPG) of Pinus massoniana belonging to the genus Pinus (Pinaceae), which is a primary source of turpentine, was sequenced and analyzed in terms of gene rearrangements, ndh genes loss, and the contraction and expansion of short inverted repeats (IRs). P. massoniana CPG has a typical quadripartite structure that includes large single copy (LSC) (65,563 bp), small single copy (SSC) (53,230 bp) and two IRs (IRa and IRb, 485 bp). The 108 unique genes were identified, including 73 protein-coding genes, 31 tRNAs, and 4 rRNAs. Most of the 81 simple sequence repeats (SSRs) identified in CPG were mononucleotides motifs of A/T types and located in non-coding regions. Comparisons with related species revealed an inversion (21,556 bp) in the LSC region; P. massoniana CPG lacks all 11 intact ndh genes (four ndh genes lost completely; the five remained truncated as pseudogenes; and the other two ndh genes remain as pseudogenes because of short insertions or deletions). A pair of short IRs was found instead of large IRs, and size variations among pine species were observed, which resulted from short insertions or deletions and non-synchronized variations between "IRa" and "IRb". The results of phylogenetic analyses based on whole CPG sequences of 16 conifers indicated that the whole CPG sequences could be used as a powerful tool in phylogenetic analyses.

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.

The Complete Chloroplast Genome Sequences of Fritillaria ussuriensis Maxim. and Fritillaria cirrhosa D. Don, and Comparative Analysis with Other Fritillaria Species

The genus Fritillaria belongs to the widely distributed Liliaceae. The bulbs of Fritillaria, F. ussuriensis and F. cirrhosa are valuable herbaceous medicinal ingredients. However, they are still used indiscriminately in herbal medicine. Identification and molecular phylogenic analysis of Fritillaria species are therefore required. Here, we report the complete chloroplast (CP) genome sequences of F. ussuriensis and F. cirrhosa. The two Fritillaria CP genomes were 151,524 and 151,083 bp in length, respectively, and each included a pair of inverted repeated regions (52,678 and 52,156 bp) that was separated by a large single copy region (81,732 and 81,390 bp), and a small single copy region (17,114 and 17,537 bp). A total of 111 genes in F. ussuriensis and 112 in F. cirrhosa comprised 77 protein-coding regions in F. ussuriensis and 78 in F. cirrhosa, 30 transfer RNA (tRNA) genes, and four ribosomal RNA (rRNA) genes. The gene order, content, and orientation of the two Fritillaria CP genomes exhibited the general structure of flowering plants, and were similar to those of other Fritillaria species. Comparison of the six Fritillaria species’ CP genomes indicated seven highly divergent regions in intergenic spacers and in the matK, rpoC1, rpoC2, ycf1, ycf2, ndhD, and ndhF coding regions. We established the position of the six species through phylogenic analysis. The complete chloroplast genome sequences of the two Fritillaria species and a comparison study are useful genomic information for identifying and for studying the phylogenetic relationship among Fritillaria species within the Liliaceae.

Complete plastid genome of Eriobotrya japonica (Thunb.) Lindl and comparative analysis in Rosaceae

Eriobotrya japonica (Thunb.) Lindl (loquat) is an evergreen Rosaceae fruit tree widely distributed in subtropical regions. Its leaves are considered as traditional Chinese medicine and are of high medical value especially for cough and emesis. Thus, we sequenced the complete plastid genome of E. japonica to better utilize this important species. The complete plastid genome of E. japonica is 159,137 bp in length, which contains a typical quadripartite structure with a pair of inverted repeats (IR, 26,326 bp) separated by large (LSC, 89,202 bp) and small (SSC, 19,283 bp) single-copy regions. The E. japonica plastid genome encodes 112 unique genes which consist of 78 protein-coding genes, 30 tRNA genes and 4 rRNA genes. Gene structure and content of E. japonica plastid genome are quite conserved and show similarity among Rosaceous species. Five large indels are unique to E. japonica in comparison with Pyrus pyrifolia and Prunus persica, which could be utilized as molecular markers. A total of 72 simple sequence repeats (SSRs) were detected and most of them are mononucleotide repeats composed of A or T, indicating a strong A or T bias for base composition. The Ka and Ks ratios of most genes are lower than 1, which suggests that most genes are under purifying selection. The phylogenetic analysis described the evolutionary relationship within Rosaceae and fully supported a close relationship between E. japonica and P. pyrifolia.

Comparative sequence analysis of Cyclospora cayetanensis apicoplast genomes originating from diverse geographical regions

BACKGROUND

Cyclospora cayetanensis is an emerging coccidian parasite that causes endemic and epidemic diarrheal disease called cyclosporiasis, and this infection is associated with consumption of contaminated produce or water in developed and developing regions. Food-borne outbreaks of cyclosporiasis have occurred almost every year in the USA since the 1990s. Investigations of these outbreaks are currently hampered due to lack of molecular epidemiological tools for trace back analysis. The apicoplast of C. cayetanensis, a relict non-photosynthetic plastid with an independent genome, provides an attractive target to discover sequence polymorphisms useful as genetic markers for detection and trace back analysis of the parasite. Distinct differences in the apicoplast genomes of C. cayetanensis could be useful in designing advanced molecular methods for rapid detection and, subtyping and geographical source attribution, which would aid outbreak investigations and surveillance studies.

METHODS

To obtain the genome sequence of the C. cayetanensis apicoplast, we sequenced the C. cayetanensis genomic DNA extracted from clinical stool samples, assembled and annotated a 34,146 bp-long circular sequence, and used this sequence as a reference genome in this study. We compared the genome and the predicted proteome to the data available from other apicomplexan parasites. To initialize the search for genetic markers, we mapped the raw sequence reads from an additional 11 distinct clinical stool samples originating from Nepal, New York, Texas, and Indonesia to the apicoplast reference genome.

RESULTS

We identified several high quality single nucleotide polymorphisms (SNPs) and small insertion/deletions spanning the apicoplast genome supported by extensive sequencing reads data, and a 30 bp sequence repeat at the terminal spacer region in a Nepalese sample. The predicted proteome consists of 29 core apicomplexan peptides found in most of the apicomplexans. Cluster analysis of these C. cayetanensis apicoplast genomes revealed a familiar pattern of tight grouping with Eimeria and Toxoplasma, separated from distant species such as Plasmodium and Babesia.

CONCLUSIONS

SNPs and sequence repeats identified in this study may be useful as genetic markers for identification and differentiation of C. cayetanensis isolates found and could facilitate outbreak investigations.

The draft genome of MD-2 pineapple using hybrid error correction of long reads

The introduction of the elite pineapple variety, MD-2, has caused a significant market shift in the pineapple industry. Better productivity, overall increased in fruit quality and taste, resilience to chilled storage and resistance to internal browning are among the key advantages of the MD-2 as compared with its previous predecessor, the Smooth Cayenne. Here, we present the genome sequence of the MD-2 pineapple (Ananas comosus (L.) Merr.) by using the hybrid sequencing technology from two highly reputable platforms, i.e. the PacBio long sequencing reads and the accurate Illumina short reads. Our draft genome achieved 99.6% genome coverage with 27,017 predicted protein-coding genes while 45.21% of the genome was identified as repetitive elements. Furthermore, differential expression of ripening RNASeq library of pineapple fruits revealed ethylene-related transcripts, believed to be involved in regulating the process of non-climacteric pineapple fruit ripening. The MD-2 pineapple draft genome serves as an example of how a complex heterozygous genome is amenable to whole genome sequencing by using a hybrid technology that is both economical and accurate. The genome will make genomic applications more feasible as a medium to understand complex biological processes specific to pineapple.

The Complete Chloroplast Genome Sequence of the Medicinal Plant Swertia mussotii Using the PacBio RS II Platform

Swertia mussotii is an important medicinal plant that has great economic and medicinal value and is found on the Qinghai Tibetan Plateau. The complete chloroplast (cp) genome of S. mussotii is 153,431 bp in size, with a pair of inverted repeat (IR) regions of 25,761 bp each that separate an large single-copy (LSC) region of 83,567 bp and an a small single-copy (SSC) region of 18,342 bp. The S. mussotii cp genome encodes 84 protein-coding genes, 37 transfer RNA (tRNA) genes, and eight ribosomal RNA (rRNA) genes. The identity, number, and GC content of S. mussotii cp genes were similar to those in the genomes of other Gentianales species. Via analysis of the repeat structure, 11 forward repeats, eight palindromic repeats, and one reverse repeat were detected in the S. mussotii cp genome. There are 45 SSRs in the S. mussotii cp genome, the majority of which are mononucleotides found in all other Gentianales species. An entire cp genome comparison study of S. mussotii and two other species in Gentianaceae was conducted. The complete cp genome sequence provides intragenic information for the cp genetic engineering of this medicinal plant.

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

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

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

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