The Complete Chloroplast Genome of Wild Rice (Oryza minuta) and Its Comparison to Related Species

Initial Characterization of the Chloroplast Genome of Vicia sepium, an Important Wild Resource Plant, and Related Inferences About Its Evolution

Lack of complete genomic information concerning Vicia sepium (Fabaceae: Fabeae) precludes investigations of evolution and populational diversity of this perennial high-protein forage plant suitable for cultivation in extreme conditions. Here, we present the complete and annotated chloroplast genome of this important wild resource plant. V. sepium chloroplast genome includes 76 protein-coding genes, 29 tRNA genes, 4 rRNA genes, and 1 pseudogene. Its 124,095 bp sequence has a loss of one inverted repeat (IR). The GC content of the whole genome, the protein-coding, intron, tRNA, rRNA, and intergenic spacer regions was 35.0%, 36.7%, 34.6%, 52.3%, 54.2%, and 29.2%, respectively. Comparative analyses with plastids from related genera belonging to Fabeae demonstrated that the greatest variation in the V. sepium genome length occurred in protein-coding regions. In these regions, some genes and introns were lost or gained; for example, ycf4, clpP intron, and rpl16 intron deletions and rpl20 and ORF292 insertions were observed. Twelve highly divergent regions, 66 simple sequence repeats (SSRs) and 27 repeat sequences were also found in these regions. Detailed evolutionary rate analysis of protein-coding genes showed that Vicia species exhibit additional interesting characteristics including positive selection of ccsA, clpP, rpl32, rpl33, rpoC1, rps15, rps2, rps4, and rps7, and the evolutionary rates of atpA, accD, and rps2 in Vicia are significantly accelerated. These genes are important candidate genes for understanding the evolutionary strategies of Vicia and other genera in Fabeae. The phylogenetic analysis showed that Vicia and Lens are included in the same clade and that Vicia is paraphyletic. These results provide evidence regarding the evolutionary history of the chloroplast genome.

Comparative Chloroplast Genomics of Endangered Euphorbia Species: Insights into Hotspot Divergence, Repetitive Sequence Variation, and Phylogeny

Euphorbia is one of the largest genera in the Euphorbiaceae family, comprising 2000 species possessing commercial, medicinal, and ornamental importance. However, there are very little data available on their molecular phylogeny and genomics, and uncertainties still exist at a taxonomic level. Herein, we sequence the complete chloroplast (cp) genomes of two species, E. larica and E. smithii, of the genus Euphorbia through next-generation sequencing and perform a comparative analysis with nine related genomes in the family. The results revealed that the cp genomes had similar quadripartite structure, gene content, and genome organization with previously reported genomes from the same family. The size of cp genomes ranged from 162,172 to 162,358 bp with 132 and 133 genes, 8 rRNAs, 39 tRNA in E. smithii and E. larica, respectively. The numbers of protein-coding genes were 85 and 86, with each containing 19 introns. The four-junction regions were studied and results reveal that rps19 was present at J_LB (large single copy region and inverted repeat b junction) in E. larica where its complete presence was located in the IRb (inverted repeat b) region in E. smithii. The sequence comparison revealed that highly divergent regions in rpoC1, rpocB, ycf3, clpP, petD, ycf1, and ndhF of the cp genomes might provide better understanding of phylogenetic inferences in the Euphorbiaceae and order Malpighiales. Phylogenetic analyses of this study illustrate sister clades of E. smithii with E. tricullii and these species form a monophyletic clade with E. larica. The current study might help us to understand the genome architecture, genetic diversity among populations, and evolutionary depiction in the genera.

Comparison among three methods for obtaining chloroplast genome sequences from the conifer Pinus massoniana

The chloroplast genome (CPG) is a powerful tool for phylogenetic studies. Many CPGs have been determined using NGS. However, the large nuclear-genome and difficult CPG-DNA separation in conifers limit their application in related research. In this study, three methods (PCR + Sanger, PCR + HiSeq, cpDNA+HiSeq) for obtaining the CPGs of Pinus massoniana were compared for sequence accuracy, time and cost. PCR + Sanger obtained the most accurate CPGs with advantages in cost (3.08$/kb) and time (2-3 days); PCR + HiSeq generated some DNA fragments with low depth, and the SNPs false-positive-rate (0.44) and sequencing error-rate (0.0265) of this method were higher than those of the cpDNA+HiSeq. Moreover, the cost (~6.17$/kb) and time (4-5 weeks) would significantly increase when HiSeq sequencing were outsourced to sequencing service company. Thus, for the study of intraspecific and interspecies variation in CPGs, CPG sequences can be obtained by comprehensive methods to bridge the method shortcomings. Scuh as sequence accuracy, cost and time.

Analysis of complete chloroplast genomes of Curcuma and the contribution to phylogeny and adaptive evolution

Curcuma is an important member of Zingiberaceae. Many species of this genus are widely used in traditional medicine and have important cultural value in East Asia. Among them, C. longa is considered to be the main source of curcumin and has a very wide range of uses. The rapid development of molecular phylogeny has deepened our understanding of taxonomy and evolution of Curcuma. However, little is known about the chloroplast genome phylogeny and the genetic bases of adaptative evolution. In this work, we sequenced the complete chloroplast genome of 4 Curcuma species. Curcuma chloroplast genomes showed highly conserved structures and the length ranged from 159,423 bp to 152,723 bp. A total of 133 genes were observed. Multiple repeats and simple sequence repeats (SSRs) were detected. By comparing with related species, 7 highly variable regions were identified as potential specific DNA barcodes for species identification. Phylogenetic analysis of complete plastome sequences and specific data sets revealed discordance with expected genus boundary. Chloroplast phylogenetic relationships were better predicted by geography than by morphological and nuclear DNA, indicating a substantial existence of introgression. 9 genes were proved to have high posteriori probability in positive selection analysis, and 4 of them (psbA, psbD, PetA and rbcL) closely related to photosynthesis, implying that chloroplast genes may had undergone positive selection pressure in evolution. These results are of great significance for us to understand the genetic basis, phylogeny and adaptive evolution of Curcuma chloroplast.

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.

Characterization of the complete chloroplast genome of Saussurea involuerata (Compositae), an endangered species endemic to China

Saussurea involuerata (S. involuerata), a kind of composite and perennial herbage plant, is a rare and precious alpine medicinal plants with high medicinal value, it can regulate human’s body and has health care function, therefore attracts attention of many scholars. It is currently on the International Union for Conservation of Nature (IUCN) red list of threatened species. In this study, we first assembled the complete chloroplast (cp) genome of S. involucrata by Illumina paired-end reads data. The whole genome was 152,490 bp, consisting of a pair of inverted repeats of 69,041 bp, large and small single-copy regions of 69,045bp and 18,638 bp in length, respectively. The cp genome contained 134 genes, including 90 protein-coding genes, 36 tRNA genes, and 8 rRNA genes. The overall GC content of the whole genome was 37.7%. A neighbor-joining phylogenetic analysis demonstrated a close relationship between S. involuerata and Arctium lappa.

Characterization of the complete chloroplast genome sequence of Dalbergia species and its phylogenetic implications

The pantropical plant genus Dalbergia comprises approximately 250 species, most of which have a high economic and ecological value. However, these species are among the most threatened due to illegal logging and the timber trade. To enforce protective legislation and ensure effective conservation of Dalbergia species, the identity of wood being traded must be accurately validated. For the rapid and accurate identification of Dalbergia species and assessment of phylogenetic relationships, it would be highly desirable to develop more effective DNA barcodes for these species. In this study, we sequenced and compared the chloroplast genomes of nine species of Dalbergia. We found that these chloroplast genomes were conserved with respect to genome size, structure, and gene content and showed low sequence divergence. We identified eight mutation hotspots, namely, six intergenic spacer regions (trnL-trnT, atpA-trnG, rps16-accD, petG-psaJ, ndhF-trnL, and ndhG-ndhI) and two coding regions (ycf1a and ycf1b), as candidate DNA barcodes for Dalbergia. Phylogenetic analyses based on whole chloroplast genome data provided the best resolution of Dalbergia, and phylogenetic analysis of the Fabaceae showed that Dalbergia was sister to Arachis. Based on comparison of chloroplast genomes, we identified a set of highly variable markers that can be developed as specific DNA barcodes.

Complete Chloroplast Genomes of Vachellia nilotica and Senegalia senegal: Comparative Genomics and Phylogenomic Placement in a New Generic System

Vachellia and Senegalia are the most important genera in the subfamily Mimosoideae (Fabaceae). Recently, species from both genera were separated from the long-characterized Acacia due to their macro-morphological characteristics. However, this morpho-taxonomic differentiation struggles to discriminate some species, for example, Vachellia nilotica and Senegalia senegal. Therefore, sequencing the chloroplast (cp) genomes of these species and determining their phylogenetic placement via conserved genes may help to validate the taxonomy. Hence, we sequenced the cp genomes of V. nilotica and S. senegal, and the results showed that the sizes of the genomes are 165.3 and 162.7 kb, respectively. The cp genomes of both species comprised large single-copy regions (93,849~91,791 bp) and pairs of inverted repeats (IR; 26,093~26,008 bp). The total numbers of genes found in the V. nilotica and S. senegal cp genomes were 135 and 132, respectively. Approximately 123:130 repeats and 290:281 simple sequence repeats were found in the S. senegal and V. nilotica cp genomes, respectively. Genomic characterization was undertaken by comparing these genomes with those of 17 species belonging to related genera in Fabaceae. A phylogenetic analysis of the whole genome dataset and 56 shared genes was undertaken by generating cladograms with the same topologies and placing both species in a new generic system. These results support the likelihood of identifying segregate genera from Acacia with phylogenomic disposition of both V. nilotica and S. senegal in the subfamily Mimosoideae. The current study is the first to obtain complete genomic information on both species and may help to elucidate the genome architecture of these species and evaluate the genetic diversity among species.

Complete chloroplast genome sequence of the mangrove species Kandelia obovata and comparative analyses with related species

As one of the most cold and salt-tolerant mangrove species, Kandelia obovata is widely distributed in China. Here, we report the complete chloroplast genome sequence K. obovata (Rhizophoraceae) obtained via next-generation sequencing, compare the general features of the sampled plastomes of this species to those of other sequenced mangrove species, and perform a phylogenetic analysis based on the protein-coding genes of these plastomes. The complete chloroplast genome of K. obovata is 160,325 bp in size and has a 35.22% GC content. The genome has a typical circular quadripartite structure, with a pair of inverted repeat (IR) regions 26,670 bp in length separating a large single-copy (LSC) region (91,156 bp) and a small single-cope (SSC) region (15,829 bp). The chloroplast genome of K. obovata contains 128 unique genes, including 80 protein-coding genes, 38 tRNA genes, 8 rRNA genes and 2 pseudogenes (ycf1 in the IRA region and rpl22 in the IRB region). In addition, a simple sequence repeat (SSR) analysis found 108 SSR loci in the chloroplast genome of K. obovata, most of which are A/T rich. IR expansion and contraction regions were compared between K. obovata and five related species: two from Malpighiales and three mangrove species from different orders. The mVISTA results indicated that the genome structure, gene order and gene content are highly conserved among the analyzed species. The phylogenetic analysis using 54 common protein-coding genes from the chloroplast genome showed that the plant most closely related to K. obovata is Ceriops tagal of Rhizophoraceae. The results of this study provide useful molecular information about the evolution and molecular biology of these mangrove trees.

Complete chloroplast genome of Castanopsis sclerophylla (Lindl.) Schott: Genome structure and comparative and phylogenetic analysis

Castanopsis sclerophylla (Lindl.) Schott is an important species of evergreen broad-leaved tree in subtropical areas and has high ecological and economic value. However, there are few studies on its chloroplast genome. In this study, the complete chloroplast genome sequence of C. sclerophylla was determined using the Illumina HiSeq 2500 platform. The complete chloroplast genome of C. sclerophylla is 160,497 bp long, including a pair of inverted repeat (IR) regions (25,675 bp) separated by a large single-copy (LSC) region of 90,255 bp and a small single-copy (SSC) region of 18,892 bp. The overall GC content of the chloroplast genome is 36.82%. A total of 131 genes were found; of these, 111 genes are unique and annotated, including 79 protein-coding genes, 27 transfer RNA genes (tRNAs), and four ribosomal RNA genes (rRNAs). Twenty-one genes were found to be duplicated in the IR regions. Comparative analysis indicated that IR contraction might be the reason for the smaller chloroplast genome of C. sclerophylla compared to three congeneric species. Sequence analysis indicated that the LSC and SSC regions are more divergent than IR regions within Castanopsis; furthermore, greater divergence was found in noncoding regions than in coding regions. The maximum likelihood phylogenetic analysis showed that four species of the genus Castanopsis form a monophyletic clade and that C. sclerophylla is closely related to Castanopsis hainanensis with strong bootstrap values. These results not only provide a basic understanding of Castanopsis chloroplast genomes, but also illuminate Castanopsis species evolution within the Fagaceae family. Furthermore, these findings will be valuable for future studies of genetic diversity and enhance our understanding of the phylogenetic evolution of Castanopsis.

Complete Chloroplast Genome of Pinus densiflora Siebold & Zucc. and Comparative Analysis with Five Pine Trees

Pinus densiflora (Korean red pine) is widely distributed in East Asia and considered one of the most important species in Korea. In this study, the complete chloroplast genome of P. densiflora was sequenced by combining the advantages of Oxford Nanopore MinION and Illumina MiSeq. The sequenced genome was then compared with that of a previously published conifer plastome. The chloroplast genome was found to be circular and comprised of a quadripartite structure, including 113 genes encoding 73 proteins, 36 tRNAs and 4 rRNAs. It had short inverted repeat regions and lacked ndh gene family genes, which is consistent with other Pinaceae species. The gene content of P. densiflora was found to be most similar to that of P. sylvestris. The newly attempted sequencing method could be considered an alternative method for obtaining accurate genetic information, and the chloroplast genome sequence of P. densiflora revealed in this study can be used in the phylogenetic analysis of Pinus species.

Characterization of the complete chloroplast genome of Ginkgo biloba L. (Ginkgoaceae), an endangered species endemic to China

Ginkgo biloba L. is the oldest relict plant among the gymnosperms, left after the quaternary glacial movement. There are few living G. biloba and few old trees over a hundred years old. It is currently on the International Union for Conservation of Nature (IUCN) red list of threatened species. In this study, we first assembled the complete chloroplast (cp) genome of G. biloba L. by Illumina paired-end reads data. The whole genome was 156,988 bp, consisting of a pair of inverted repeats of 34,056 bp, large and small single copy regions of 56,819 and 22,763 bp in length, respectively. The cp genome contained 131 genes, including 84 protein-coding genes, 31 trRNA genes and 5 rRNA genes. The overall GC content of the whole genome was 39.6%. A neighbour-joining phylogenetic analysis demonstrated a close relationship between G. biloba L. and Cycas revoluta.

The Complete Chloroplast Genome of Euphrasia regelii, Pseudogenization of ndh Genes and the Phylogenetic Relationships Within Orobanchaceae

Euphrasia (Orobanchaceae) is a genus which is widely distributed in temperate regions of the southern and northern hemisphere. The taxonomy of Euphrasia is still controversial due to the similarity of morphological characters and a lack of genomic resources. Here, we present the first complete chloroplast (cp) genome of this taxonomically challenging genus. The cp genome of Euphrasia regelii consists of 153,026 bp, including a large single-copy region (83,893 bp), a small single-copy region (15,801 bp) and two inverted repeats (26,666 bp). There are 105 unique genes, including 71 protein-coding genes, 30 tRNA and 4 rRNA genes. Although the structure and gene order is comparable to the one in other angiosperm cp genomes, genes encoding the NAD(P)H dehydrogenase complex are widely pseudogenized due to mutations resulting in frameshifts, and stop codon positions. We detected 36 dispersed repeats, 7 tandem repeats and 65 simple sequence repeat loci in the E. regelii plastome. Comparative analyses indicated that the cp genome of E. regelii is more conserved compared to other hemiparasitic taxa in the Pedicularideae and Buchnereae. No structural rearrangements or loss of genes were detected. Our analyses suggested that three genes (clpP, ycf2 and rps14) were under positive selection and other genes under purifying selection. Phylogenetic analysis of monophyletic Orobanchaceae based on 45 plastomes indicated a close relationship between E. regelii and Neobartsia inaequalis. In addition, autotrophic lineages occupied the earliest diverging branches in our phylogeny, suggesting that autotrophy is the ancestral trait in this parasitic family.

Complete Chloroplast Genome Sequences of Kaempferia Galanga and Kaempferia Elegans: Molecular Structures and Comparative Analysis

Kaempferia galanga and Kaempferia elegans, which belong to the genus Kaempferia family Zingiberaceae, are used as valuable herbal medicine and ornamental plants, respectively. The chloroplast genomes have been used for molecular markers, species identification and phylogenetic studies. In this study, the complete chloroplast genome sequences of K. galanga and K. elegans are reported. Results show that the complete chloroplast genome of K. galanga is 163,811 bp long, having a quadripartite structure with large single copy (LSC) of 88,405 bp and a small single copy (SSC) of 15,812 bp separated by inverted repeats (IRs) of 29,797 bp. Similarly, the complete chloroplast genome of K. elegans is 163,555 bp long, having a quadripartite structure in which IRs of 29,773 bp length separates 88,020 bp of LSC and 15,989 bp of SSC. A total of 111 genes in K. galanga and 113 genes in K. elegans comprised 79 protein-coding genes and 4 ribosomal RNA (rRNA) genes, as well as 28 and 30 transfer RNA (tRNA) genes in K. galanga and K. elegans, respectively. The gene order, GC content and orientation of the two Kaempferia chloroplast genomes exhibited high similarity. The location and distribution of simple sequence repeats (SSRs) and long repeat sequences were determined. Eight highly variable regions between the two Kaempferia species were identified and 643 mutation events, including 536 single-nucleotide polymorphisms (SNPs) and 107 insertion/deletions (indels), were accurately located. Sequence divergences of the whole chloroplast genomes were calculated among related Zingiberaceae species. The phylogenetic analysis based on SNPs among eleven species strongly supported that K. galanga and K. elegans formed a cluster within Zingiberaceae. This study identified the unique characteristics of the entire K. galanga and K. elegans chloroplast genomes that contribute to our understanding of the chloroplast DNA evolution within Zingiberaceae species. It provides valuable information for phylogenetic analysis and species identification within genus Kaempferia.

First complete chloroplast genomics and comparative phylogenetic analysis of Commiphora gileadensis and C. foliacea: Myrrh producing trees

Commiphora gileadensis and C. foliacea (family Burseraceae) are pantropical in nature and known for producing fragrant resin (myrrh). Both the tree species are economically and medicinally important however, least genomic understanding is available for this genus. Herein, we report the complete chloroplast genome sequences of C. gileadensis and C. foliacea and comparative analysis with related species (C. wightii and Boswellia sacra). A modified chloroplast DNA extraction method was adopted, followed with next generation sequencing, detailed bioinformatics and PCR analyses. The results revealed that the cp genome sizes of C. gileadensis and C. foliacea, are 160,268 and 160,249 bp, respectively, with classic quadripartite structures that comprises of inverted repeat's pair. Overall, the organization of these cp genomes, GC contents, gene order, and codon usage were comparable to other cp genomes in angiosperm. Approximately, 198 and 175 perfect simple sequence repeats were detected in C. gileadensis and C. foliacea genomes, respectively. Similarly, 30 and 25 palindromic, 15 and 25 forward, and 20 and 25 tandem repeats were determined in both the cp genomes, respectively. Comparison of these complete cp genomes with C. wightii and B. sacra revealed significant sequence resemblance and comparatively highest deviation in intergenic spacers. The phylo-genomic comparison showed that C. gileadensis and C. foliacea form a single clade with previously reported C. wightii and B. sacra from family Burseraceae. Current study reports for the first time the cp genomics of species from Commiphora, which could be helpful in understanding genetic diversity and phylogeny of this myrrh producing species.

Green Biotechnology: A Brief Update on Plastid Genome Engineering

Plant genetic engineering has become an inevitable tool in the molecular breeding of crops. Significant progress has been made in the generation of novel plastid transformation vectors and optimized transformation protocols. There are several advantages of plastid genome engineering over conventional nuclear transformation. Some of the advantages include multigene engineering by expression of biosynthetic pathway genes as operons, extremely high-level expression of protein accumulation, lack of transgene silencing, etc. Transgene containment owing to maternal inheritance is another important advantage of plastid genome engineering. Chloroplast genome modification usually results in alteration of several thousand plastid genome copies in a cell. Several therapeutic proteins, edible vaccines, antimicrobial peptides, and industrially important enzymes have been successfully expressed in chloroplasts so far. Here, we critically recapitulate the latest developments in plastid genome engineering. Latest advancements in plastid genome sequencing are briefed. In addition, advancement of extending the toolbox for plastid engineering for selected applications in the area of molecular farming and production of industrially important enzyme is briefed.

Draft genome sequence of first monocot-halophytic species Oryza coarctata reveals stress-specific genes

Oryza coarctata (KKLL; 2n = 4x = 48, 665 Mb) also known as Porteresia coarctata is an extreme halophyte species of genus Oryza. Using Illumina and Nanopore reads, we achieved the assembled genome size of 569.9 Mb, accounting 85.69% of the estimated genome size with N50 of 1.85 Mb and 19.89% repetitive region. We also found 230,968 simple sequence repeats (SSRs) and 5,512 non-coding RNAs (ncRNAs). The functional annotation of predicted 33,627 protein-coding genes and 4,916 transcription factors revealed that high salinity adaptation of this species is due to the exclusive or excessive presence of stress-specific genes as compared to rice. We have identified 8 homologs to salt-tolerant SOS1 genes, one of the three main components of salt overly sensitive (SOS) signal pathway. On the other hand, the phylogenetic analysis of the assembled chloroplast (134.75 kb) and mitochondrial genome (491.06 kb) favours the conservative nature of these organelle genomes within Oryza taxon.

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.

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.

Comparative Analysis of the Chloroplast Genomes of the Chinese Endemic Genus Urophysa and Their Contribution to Chloroplast Phylogeny and Adaptive Evolution

Urophysa is a Chinese endemic genus comprising two species, Urophysa rockii and Urophysa henryi. In this study, we sequenced the complete chloroplast (cp) genomes of these two species and of their relative Semiquilegia adoxoides. Illumina sequencing technology was used to compare sequences, elucidate the intra- and interspecies variations, and infer the phylogeny relationship with other Ranunculaceae family species. A typical quadripartite structure was detected, with a genome size from 158,473 to 158,512 bp, consisting of a pair of inverted repeats separated by a small single-copy region and a large single-copy region. We analyzed the nucleotide diversity and repeated sequences components and conducted a positive selection analysis by the codon-based substitution on single-copy coding sequence (CDS). Seven regions were found to possess relatively high nucleotide diversity, and numerous variable repeats and simple sequence repeats (SSR) markers were detected. Six single-copy genes (atpA, rpl20, psaA, atpB, ndhI, and rbcL) resulted to have high posterior probabilities of codon sites in the positive selection analysis, which means that the six genes may be under a great selection pressure. The visualization results of the six genes showed that the amino acid properties across each column of all species are variable in different genera. All these regions with high nucleotide diversity, abundant repeats, and under positive selection will provide potential plastid markers for further taxonomic, phylogenetic, and population genetics studies in Urophysa and its relatives. Phylogenetic analyses based on the 79 single-copy genes, the whole complete genome sequences, and all CDS sequences showed same topologies with high support, and U. rockii was closely clustered with U. henryi within the Urophysa genus, with S. adoxoides as their closest relative. Therefore, the complete cp genomes in Urophysa species provide interesting insights and valuable information that can be used to identify related species and reconstruct their phylogeny.

Dynamic evolution and phylogenomic analysis of the chloroplast genome in Schisandraceae

Chloroplast genomes of plants are highly conserved in both gene order and gene content, are maternally inherited, and have a lower rate of evolution. Chloroplast genomes are considered to be good models for testing lineage-specific molecular evolution. In this study, we use Schisandraceae as an example to generate insights into the overall evolutionary dynamics in chloroplast genomes and to establish the phylogenetic relationship of Schisandraceae based on chloroplast genome data using phylogenomic analysis. By comparing three Schisandraceae chloroplast genomes, we demonstrate that the gene order, gene content, and length of chloroplast genomes in Schisandraceae are highly conserved but experience dynamic evolution among species. The number of repeat variations were detected, and the Schisandraceae chloroplast genome was revealed as unusual in having a 10 kb contraction of the IR due to the genome size variations compared with other angiosperms. Phylogenomic analysis based on 82 protein-coding genes from 66 plant taxa clearly elucidated that Schisandraceae is a sister to a clade that includes magnoliids, monocots, and eudicots within angiosperms. As to genus relationships within Schisandraceae, Kadsura and Schisandra formed a monophyletic clade which was sister to Illicium.

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.

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.

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.

Complete Chloroplast Genome Sequence of Decaisnea insignis: Genome Organization, Genomic Resources and Comparative Analysis

Decaisnea insignis is a wild resource plant and is used as an ornamental, medicinal, and fruit plant. High-throughput sequencing of chloroplast genomes has provided insight into the overall evolutionary dynamics of chloroplast genomes and has enhanced our understanding of the evolutionary relationships within plant families. In the present study, we sequenced the complete chloroplast genome of D. insignis and used the data to assess its genomic resources. The D. insignis chloroplast genome is 158,683 bp in length and includes a pair of inverted repeats of 26,167 bp that are separated by small and large single copy regions of 19,162 bp and 87,187 bp, respectively. We identified 83 simple sequence repeats and 18 pairs of large repeats. Most simple-sequence repeats were located in the noncoding sections of the large single-copy/small single-copy region and exhibited a high A/T content. The D. insignis chloroplast genome bias was skewed towards A/T on the basis of codon usage. A phylogenetic tree based on 82 protein-coding genes of 33 angiosperms showed that D. insignis was clustered with Akebia in Lardizabalaceae. Overall, the results of this study will contribute to better understanding the evolution, molecular biology and genetic improvement of D. insignis.

Chloroplast genomes of Arabidopsis halleri ssp. gemmifera and Arabidopsis lyrata ssp. petraea: Structures and comparative analysis

We investigated the complete chloroplast (cp) genomes of non-model Arabidopsis halleri ssp. gemmifera and Arabidopsis lyrata ssp. petraea using Illumina paired-end sequencing to understand their genetic organization and structure. Detailed bioinformatics analysis revealed genome sizes of both subspecies ranging between 154.4~154.5 kbp, with a large single-copy region (84,197~84,158 bp), a small single-copy region (17,738~17,813 bp) and pair of inverted repeats (IRa/IRb; 26,264~26,259 bp). Both cp genomes encode 130 genes, including 85 protein-coding genes, eight ribosomal RNA genes and 37 transfer RNA genes. Whole cp genome comparison of A. halleri ssp. gemmifera and A. lyrata ssp. petraea, along with ten other Arabidopsis species, showed an overall high degree of sequence similarity, with divergence among some intergenic spacers. The location and distribution of repeat sequences were determined, and sequence divergences of shared genes were calculated among related species. Comparative phylogenetic analysis of the entire genomic data set and 70 shared genes between both cp genomes confirmed the previous phylogeny and generated phylogenetic trees with the same topologies. The sister species of A. halleri ssp. gemmifera is A. umezawana, whereas the closest relative of A. lyrata spp. petraea is A. arenicola.

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.