Cp Genome Research Articles

Characterization and phylogenetic analysis of the complete chloroplast genome of Emilia sonchifolia (L.) DC.

Emilia sonchifolia is a herb with antioxidant, anti-inflammatory, antitumor, and wound healing properties. The complete chloroplast genome (cp genome) of the genus Emilia was sequenced for the first time. The cp genome of E. sonchifolia is 151,474 bp in length. It contained a large single-copy (LSC) region (84,004 bp), and small single-copy (SSC) region (17,980 bp), and two inverted repeats (IRs, 24,745 bp). Phylogenetic analysis of 24 species was conducted. E. sonchifolia was found to be closely related to Pericallis hybrida and Dendrosenecio spp. The sequenced cp genome would be useful to understand the phylogeny and genomic studies of the genus Emilia.

Structural characterization of four Rhododendron spp. chloroplast genomes and comparative analyses with other azaleas

Azalea is a general designation of Rhododendron in the Ericaceae family. Rhododendron not only has high ornamental value but also has application value in ecological protection, medicine, and scientific research. In this study, we used Illumina and PacBio sequencing to assemble and annotate the entire chloroplast genomes (cp genomes) of four Rhododendron species. The chloroplast genomes of R. concinnum, R. henanense subsp. lingbaoense, R. micranthum, and R. simsii were assembled into 207,236, 208,015, 207,233, and 206,912 bp, respectively. All chloroplast genomes contain eight rRNA genes, with either 88 or 89 protein-coding genes. The four cp genomes were compared and analyzed by bioinformatics, and the phylogenetic analysis based on chloroplast genomes of 26 species of Ericaceae, Actinidiaceae, and Primulaceae under Ericales was conducted. A comparison of the linear structure of cp genomes of four Rhododendron showed that there were substantial sequence similarities in coding regions, but high differences in non-coding regions. A phylogenetic analysis, based on chloroplast whole genome sequences, showed that all Rhododendron species are in the clade Ericaceae. This study provides valuable genetic information for the study of population genetics and evolutionary relationships in Rhododendron and other azalea species.

The complete chloroplast genome sequence of Vincetoxicum mongolicum (Apocynaceae), a perennial medicinal herb.

Vincetoxicum mongolicum Maxim. (1876), is a perennial medicinal herb, widely distributed in the Loess Plateau of China. Here, we sequenced, assembled, and annotated the complete chloroplast (cp) genome of V. mongolicum, and compared the highly variable gene regions and phylogenetic positions between V. mongolicum and other related species. Results showed that the complete cp genome of V. mongolicum was 160,157 bp in length, containing a large single copy (LSC) region of 91,263 bp, a pair of inverted repeats (IR) region of 23,892 bp, and a small single copy (SSC) region of 21,110 bp. The GC content accounts for 37.8%, and we annotated 131 single genes, which include 86 protein-coding genes, 8 rRNA genes, and 37 tRNA genes. By comparing and analyzing the variable region of the cp gene of V. mongolicum and other Vincetoxicum, we found that the variable sequences of rpoC1-rpoB, ycf4-cemA, ndhF, ndhF-rpl32, and rpl32-ccsA fragments were highly significant, which could be targeted as the DNA barcodes for evidence of V. mongolicum and its relatives in Apocynaceae. Maximum-likelihood (ML) phylogenetic tree analysis elucidated that V. mongolicum was sister to V. pycnostelma with strong support. Our results provide useful information for future phylogenetic studies and plastid super-barcodes of the family Apocynaceae.

Extensive characterization of 28 complete chloroplast genomes of Hydrangea species: A perspective view of their organization and phylogenetic and evolutionary relationships

The tribe Hydrangeeae displays a unique, distinctive disjunct distribution encompassing East Asia, North America and Hawaii. Despite its complex trait variations and polyphyletic nature, comprehensive phylogenomic and biogeographical studies on this tribe have been lacking. To address this gap, we sequenced and characterized 28 plastomes of Hydrangeeae. Our study highlights the highly conserved nature of Hydrangeaceae chloroplast (cp) genomes in terms of gene content and arrangement. Notably, synapomorphic characteristics of tandem repeats in the conserved domain of accD were observed in the Macrophyllae, Chinenses, and Dichroa sections within the Hydrangeeae tribe. Additionally, we found lower expression of accD in these sections using structure prediction and quantitative real-time PCR analysis. Phylogenomic analyses revealed the subdivision of the Hydrangeeae tribe into two clades with robust support values. Consistent with polyphyletic relationships, sect. Broussaisia was identified as the basal group in the tribe Hydrangeeae. Our study also provides insights into the phylogenetic relationships of Hydrangea petiolaris in the Jeju and Ulleung Island populations, suggesting the need for further studies with more samples and molecular data. Divergence time estimation and biogeographical analyses suggested that the common ancestors of the tribe Hydrangeeae likely originated from North America and East Asia during the Paleocene period via the Bering Land Bridge, potentially facilitating migration within the tribe between these regions. In conclusion, this study enhances our understanding of the evolutionary history and biogeography of the tribe Hydrangeeae, shedding light on the dispersal patterns and origins of this intriguing plant group with its unique disjunct distribution.

The complete chloroplast genome of Primula calliantha subsp. bryophila, an ornamental alpine plant from China

Primula calliantha subsp. bryophila (Balf. f. et Farrer) W.W. Smith and Forrest (1928) is a perennial alpine species with ornamental value. It is distributed in northwestern Yunnan and adjacent eastern Tibet of China, and northern Myanmar. Here, we sequenced and assembled complete plastid genome of P. calliantha subsp. bryophila, which is a circular molecule of 152,045 bp in length, including a large single-copy region (83,966 bp), a small single-copy region (17,663 bp), and a pair of inverted repeats (25,208 bp). The chloroplast genome contained 113 genes, including 79 protein-coding genes, four rRNA genes, and 30 tRNA genes. The phylogenetic tree based on chloroplast genomes showed the relative relationship of P. calliantha subsp. bryophila and P. calliantha, which further supports P. calliantha subsp. bryophila as a subspecies of P. calliantha in taxonomy. The complete chloroplast (cp) genome of P. calliantha subsp. bryophila provides valuable data for further phylogenetic studies of Primulaceae.

The complete chloroplast genome sequence of Rhododendron fortunei: Structural comparative and phylogenetic analysis in the ericaceae family

Rhododendron fortunei (Ericaceae) possesses valuable horticultural and medicinal values. However, the genomic information on R. fortunei is very limited. In this study, the complete chloroplast genome (cp) of R. fortunei was assembled and annotated, SSR loci were characterised, comparative genomic analysis was carried out, and phylogenetic research was also performed. The results showed that the R. fortunei cp genome was of a typical quadripartite structure (200,997 bp). The lengths of the large single copy region (LSC), the inverted repeat regions (IR), and the small single copy region (SSC) were 109,151 bp, 2,604 bp, and 44,619 bp, respectively. A total of 147 unique genes were identified, including 99 protein-coding genes, 42 tRNA genes, and 6 rRNA genes, respectively. Leucine (11.51%) and cysteine (1.15%) were the highest and lowest representative amino acids, respectively. The total of 30 codons with obvious codon usage bias were all A/U-ending codons. Among the 77 simple sequence repeats, the majority were mononucleotide A/T repeats located in the intergenic spacer region. Five gene regions showed high levels of nucleotide diversity (Pi > 0.03). The comparative genome analysis revealed 7 hotspot intergenic regions (trnI-rpoB, trnTrpl16, rpoA-psbJ, rps7-rrn16, ndhI-rps16, rps16-rps19, and rrn16-trnI), showing great potential as molecular makers for species authentication. Expansion and contraction were detected in the IR region of the R. fortunei cp genome. In the phylogenetic tree, R. fortunei was closely related to R. platypodum. This research will be beneficial for evolutionary and genetic diversity studies of R. fortunei and related species among the Ericaceae family.

Chloroplast genome skimming of a potential agroforestry species Melia dubia. Cav and its comparative phylogenetic analysis with major Meliaceae members.

The online version contains supplementary material available at 10.1007/s13205-022-03447-1.

Complete chloroplast genome sequence of <i>Clematis calcicola</i> (Ranunculaceae), a species endemic to Korea

The complete chloroplast genome (cp genome) sequence of <i>Clematis calcicola</i> J. S. Kim (Ranunculaceae) is 159,655 bp in length. It consists of large (79,451 bp) and small (18,126 bp) single-copy regions and a pair of identical inverted repeats (31,039 bp). The genome contains 92 protein-coding genes, 36 transfer RNA genes, eight ribosomal RNA genes, and two pseudogenes. A phylogenetic analysis based on the cp genome of 19 taxa showed high similarity between our cp genome and data published for <i>C. calcicola</i>, which is recognized as a species endemic to the Korean Peninsula. The complete cp genome sequence of <i>C. calcicola</i> reported here provides important information for future phylogenetic and evolutionary studies of Ranunculaceae.

The Complete Chloroplast Genomes of Two Physalis Species, Physalis macrophysa and P. ixocarpa: Comparative Genomics, Evolutionary Dynamics and Phylogenetic Relationships

Physalis macrophysa (PMA) and Physalis ixocarpa (PIX) have significant economic, medicinal, and ornamental value and are often used in the fields of edible fruits, medicinal herbs, and ornamental plants. In the present study, we obtained the complete chloroplast (cp) genome sequences of PMA and PIX, compared to the cp genome sequences of 10 Physalis species, and constructed the phylogenetic tree among the tribe Physaleae. The results showed that the cp genomes of PMA and PIX consisted of a large single copy (LSC) region (87,115 bp and 87,103 bp, respectively), a small single copy (SSC) region (18,412 bp and 18,420 bp, respectively), and a pair of same-length inverted-repeat (IRa and IRb) regions (25,604 bp and 25,674 bp, respectively). The two species contained 132 genes, including 87 encoding proteins, eight encoding ribosomal RNAs (rRNAs), and 37 encoding transfer RNAs (tRNAs), which indicated that the two species have strong similarities with respect to genome structure and gene content. PMA and PIX contained repeat sequences (35 and 40, respectively) and simple-sequence repeats (SSRs) (61 and 60, respectively). Nine regions with considerable nucleotide divergence were found, most of which were located in the LSC and SSC regions. The gene selective pressure analysis indicated that eight genes were affected by positive selection, the Ka/Ks values of which were greater than one. Our phylogenetic results indicated that PMA and PIX had the closest genetic relationship and are closely adjacent to Physalis philadelphica (PPH) in the subtribe Physalinae. Our analysis of the cp genomes in both Physalis species will be beneficial for further research into species identification, phylogeny, evolution, and the potential for germplasm resource exploitation in Physalis.

Phylogenetic Analysis of Wild Pomegranate (Punica granatum L.) Based on Its Complete Chloroplast Genome from Tibet, China

Pomegranates (Punica granatum L.) are one of the most fashionable fruits and can be cultivated in both tropical and subtropical regions of the world. It is a shiny and attractive pome due to its cultivation. It belongs to the Lythraceae family. In this study, we analyzed the chloroplast genome of wild pomegranate based on whole genome shotgun sequences. In wild pomegranate, we found 158,645 bp in size, 132 genes containing 84 protein genes, 37 tRNA genes, 8 rRNA genes, and 36.92% of GC content, one infA and two duplicated ycf15 pseudogenes. Moreover, 21 chloroplast genes contained intros that are detected in a large single copy (LSC), small single copy (SSC), and two inverted repeats (IRA and IRB) regions, 17 of which were involved in single introns, while four genes (ycf3, rps12, clpP and rsp12) located in LSC, IRA, and IRB region. In total, 26,272 codons are found in protein-coding genes (PCGs); relative synonymous codon usage (RSCU) analysis revealed that the most abundant amino acid is leucine containing 2773 codons (10.55%), less abundant is methionine amino acid containing 1 codon (0.0032) in the PCGs. Furthermore, a total of 233 cpSSRs were identified in the wild pomegranate cp genome, and their distribution was analyzed in three regions, namely IR, LSC, and SSC. However, 155 cpSSR were found in the LSC (66.5%), followed by 40 cpSSR in the SSC (17.2%) and 38 cpSSR in the IR (16.3%) regions. Phylogenetic validation revealed that wild pomegranate is close to the pemphis acidula species. We believe that the cp genome allocates significant information promising for breeding research of wild pomegranate to Lythraceae.

The Complete Chloroplast Genome Sequence of Eupatorium fortunei: Genome Organization and Comparison with Related Species.

Eupatorium fortunei Turcz, a perennial herb of the Asteraceae family, is one of the horticultural and medicinal plants used for curing various diseases and is widely distributed in China and other Asian countries. It possesses antibacterial, antimetastatic, antiangiogenic, and antioxidant properties along with anticancer potential. However, the intrageneric classification and phylogenetic relationships within Eupatorium have long been controversial due to the lack of high-resolution molecular markers, and the complete chloroplast (cp) genome sequencing has not been reported with new evolutionary insights. In the present study, E. fortunei was used as an experimental material, and its genome was sequenced using high-throughput sequencing technology. We assembled the complete cp genome, and a systematic analysis was conducted for E. fortunei, acquiring the correspondence of its NCBI accession number (OK545755). The results showed that the cp genome of E. fortunei is a typical tetrad structure with a total length of 152,401 bp, and the genome encodes 133 genes. Analysis of the complete cp genomes of 20 Eupatorieae shows that the number of simple sequence repeats (SSRs) ranged from 19 to 36 while the number of long sequence repeats was 50 in all cases. Eleven highly divergent regions were identified and are potentially useful for the DNA barcoding of Eupatorieae. Phylogenetic analysis among 22 species based on protein-coding genes strongly supported that E. fortunei is more closely related to Praxelis clematidea and belongs to the same branch. The genome assembly and analysis of the cp genome of E. fortunei will facilitate the identification, taxonomy, and utilization of E. fortunei as well as provide more accurate evidence for the taxonomic identification and localization of Asteraceae plants.

Complete chloroplast genome and phylogenetic analysis of Anemone shikokiana.

Anemone shikokiana (Makino) Makino, disjunctly distributed in Shandong Peninsula of China and Shikoku Island of Japan, is a rare and endangered species. To provide genetic information and understand its phylogeny, we conducted research on the chloroplast (cp) genome of A. shikokiana. The complete cp genome sequence of A. shikokiana was constructed in this study. The results showed that the cp genome of A. shikokiana has a typical quadripartite cyclic with a total length of 159,286bp. In total, 111 unique genes were identified, including 78 protein-coding genes, 29 tRNA-coding genes and 4 rRNA-coding genes. A total of 37 long repeat sequences and 67 microsatellites were found in this cp genome. The cp genome of A. shikokiana was compared with eleven other Anemone cp genomes available from the Genbank database. We found some variations among the different genomes, especially in the LSC and SSC regions, and identified some regions as potential molecular markers such as ycf1, ndhE, ndhD, ndhF-trnL, ndhA and ndhF. The results of phylogenetic analysis suggested that A. narcissiflora was the closest relative of A. shikokiana. The results filled the gap of cp genome sequence information of A. shikokiana, laying the foundation to explore the evolutionary relationships of A. shikokiana in future studies. It provided a valuable genetic resource for the molecular identification and phylogenetic study of Anemone.

Complete Chloroplast Genome Sequence of Endangered Species in the Genus Opisthopappus C. Shih: Characterization, Species Identification, and Phylogenetic Relationships.

Opisthopappus C. Shih is a rare genus of the Asteraceae family native to the Taihang Mountains in China. Due to the narrow distribution area, poor reproduction ability and human harvesting, Opisthopappus is threatened by extinction. However, the limited genetic information within Opisthopappus impede understanding of the conservation efforts and bioprospecting. Therefore, in this study, we reported the complete chloroplast (cp) genome sequences of two Opisthopappus species, including Opisthopappus taihangensis and Opisthopappus longilobus. The cp genomes of O. taihangensis and O. longilobus were 151,117 and 151,123 bp, which contained 88 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. The repeat sequences, codon usage, RNA-editing sites, and comparative analyses revealed a high degree of conservation between the two species. The ycf1 gene was identified as a potential molecular marker. The phylogenetic tree demonstrated that O. longilobus was a separate species and not a synonym or variety of O. taihangensis. The molecular clock showed that two species diverge over a large time span, O. longilobus diverged at 15.24 Mya (Million years ago), whereas O. taihangensis diverged at 5.40 Mya We found that Opisthopappus and Ajania are closely related, which provides new ideas for the development of Opisthopappus. These results provide biological information and an essential basis to understand the evolutionary history of the Opisthopappus species, which will aid in the future the bioprospecting and conservation of endangered species.

The Complete Chloroplast Genome Sequence of Machilus chuanchienensis (Lauraceae): Genome Structure and Phylogenetic Analysis

Machilus chuanchienensis is an ecological tree distributed in southwestern China. It has a significant valuation with making Hawk tea using its leaves, an ethnic traditional tea-like beverage with a long history in Chinese tea culture. The whole chloroplast (cp) genome is an ideal model for the phylogenetic study of Lauraceae because of its simple structure and highly conserved features. There have been numerous reports of complete cp genome sequences in Lauraceae, but little is known about M. chuanchienensis. Here, the next-generation sequencing (NGS) was used to sequence the M. chuanchienensis cp genome. Then, a comprehensive comparative genome analysis was performed. The results revealed that the M. chuanchienensis's cp genome measured 152,748 base pairs (bp) with a GC content of 39.15% and coded 126 genes annotated, including comprising eight ribosomal RNA (rRNA), 36 transporter RNA (tRNA), and 82 protein-coding genes. In addition, the cp genome presented a typical quadripartite structure comprising a large single-copy (LSC; 93,811) region, a small single-copy (SSC; 18,803) region, and the inverted repeats (IRs; 20,067) region and contained 92 simple sequence repeat (SSR) locus in total. Phylogenetic relationships of 37 species indicated that M. chuanchienensis was a sister to M. balansae, M. melanophylla, and M. minutiflora. Further research on this crucial species may benefit significantly from these findings.

Comparative and phylogenetic analysis of the complete chloroplast genome sequences of Allium mongolicum

Allium mongolicum Regel is a wild and sandy vegetable with unique flavours. In this study, a complete chloroplast (cp) genome of A. mongolicum was obtained (Genbank accession number: OM630416), and contained 153,609 base pairs with the GC ratio as 36.8%. 130 genes were annotated including 84 protein-coding genes, 38 tRNA, and 8 rRNA genes. The large single-copy (LSC) region was 82,644 bp, and a small single-copy (SSC) region was 18,049 bp, which were separated by two inverted repeats (IRs, including IRa and IRb) of 26,458 bp. Comparative genome analyses of 55 Allium species suggested that genomic structure of genus Allium was conserved, and LSC and SSC regions were outstanding with high variability. Among them, more divergent loci were in the SSC region covering ycf1-rrn4.5 and ndhF-ccsA. Phylogenetic analysis on cp genomes of 55 Allium determined that all members were clustered into 13 clades, and A. mongolicum had close relationship with A. senescens. Corresponding analyses of four protein-coding genes (ycf1, ndhF, rpl32, and ccsA) in aforementioned divergent loci confirmed that ycf1 was finally chosen as the candidate gene for species identification and evolutionary classification of genus Allium. These data provide valuable genetic resources for future research on Allium.

Characterization and molecular evolution analysis of Periploca forrestii inferred from its complete chloroplast genome sequence

Periploca forrestii, a medicinal plant of the family Apocynaceae, is known as an effective and widely used clinical prescription for the treatment of rheumatoid diseases. In this study, we de novo sequenced and assembled the completement chloroplast (cp) genome of P. forrestii based on combined Oxford Nanopore PromethION and Illumina data. The cp genome was 153 724bp in length and had four subregions. Moreover, an 84 433 bp large single-copy and a 17 731 bp small single-copy were separated by 25 780 bp inverted repeats (IRs). The cp genome included 132 genes with 18 duplicates in the IRs. A total of 45 repeat structures and 183 simple sequence repeats were detected. Codon usage showed a bias toward A/T-ending codons. A comparative study of Apocynaceae revealed that an IR expansion occurred on P. forrestii. The Ka/Ks values of eight species of Apocynaceae suggested that positive selection was exerted on the psaI and ycf2 genes, which might reflect specific adaptions to the P. forrestii particular growth environment. Phylogenetic analysis indicated that Periplocoideae was a sister to Asclepiadoideae, forming a monophyletic group in the family Apocynaceae. This study provided an important P. forrestii genomic resource for future evolutionary studies and the phylogenetic reconstruction of the family Apocynaceae.

Chloroplast gene expression level is negatively correlated with evolutionary rates and selective pressure while positively with codon usage bias in Ophioglossum vulgatum L

BackgroundCharacterization of the key factors determining gene expression level has been of significant interest. Previous studies on the relationship among evolutionary rates, codon usage bias, and expression level mostly focused on either nuclear genes or unicellular/multicellular organisms but few in chloroplast (cp) genes. Ophioglossum vulgatum is a unique fern and has important scientific and medicinal values. In this study, we sequenced its cp genome and transcriptome to estimate the evolutionary rates (dN and dS), selective pressure (dN/dS), gene expression level, codon usage bias, and their correlations.ResultsThe correlation coefficients between dN, dS, and dN/dS, and Transcripts Per Million (TPM) average values were -0.278 (P = 0.027 < 0.05), -0.331 (P = 0.008 < 0.05), and -0.311 (P = 0.013 < 0.05), respectively. The codon adaptation index (CAI) and tRNA adaptation index (tAI) were significantly positively correlated with TPM average values (P < 0.05).ConclusionsOur results indicated that when the gene expression level was higher, the evolutionary rates and selective pressure were lower, but the codon usage bias was stronger. We provided evidence from cp gene data which supported the E-R (E stands for gene expression level and R stands for evolutionary rate) anti-correlation.

A Large Intergenic Spacer Leads to the Increase in Genome Size and Sequential Gene Movement around IR/SC Boundaries in the Chloroplast Genome of Adiantum malesianum (Pteridaceae)

Expansion and contraction (ebb and flow events) of inverted repeat (IR) boundaries occur and are generally considered to be major factors affecting chloroplast (cp) genome size changes. Nonetheless, the Adiantum malesianum cp genome does not seem to follow this pattern. We sequenced, assembled and corrected the A. flabellulatum and A. malesianum cp genomes using the Illumina NovaSeq6000 platform, and we performed a comparative genome analysis of six Adiantum species. The results revealed differences in the IR/SC boundaries of A. malesianum caused by a 6876 bp long rpoB-trnD-GUC intergenic spacer (IGS) in the LSC. This IGS may create topological tension towards the LSC/IRb boundary in the cp genome, resulting in a sequential movement of the LSC genes. Consequently, this leads to changes of the IR/SC boundaries and may even destroy the integrity of trnT-UGU, which is located in IRs. This study provides evidence showing that it is the large rpoB-trnD-GUC IGS that leads to A. malesianum cp genome size change, rather than ebb and flow events. Then, the study provides a model to explain how the rpoB-trnD-GUC IGS in LSC affects A. malesianum IR/SC boundaries. Moreover, this study also provides useful data for dissecting the evolution of cp genomes of Adiantum. In future research, we can expand the sample to Pteridaceae to test whether this phenomenon is universal in Pteridaceae.

Extensive reorganization of the chloroplast genome of Corydalis platycarpa: A comparative analysis of their organization and evolution with other Corydalis plastomes.

The chloroplast (cp) is an autonomous plant organelle with an individual genome that encodes essential cellular functions. The genome architecture and gene content of the cp is highly conserved in angiosperms. The plastome of Corydalis belongs to the Papaveraceae family, and the genome is comprised of unusual rearrangements and gene content. Thus far, no extensive comparative studies have been carried out to understand the evolution of Corydalis chloroplast genomes. Therefore, the Corydalis platycarpa cp genome was sequenced, and wide-scale comparative studies were conducted using publicly available twenty Corydalis plastomes. Comparative analyses showed that an extensive genome rearrangement and IR expansion occurred, and these events evolved independently in the Corydalis species. By contrast, the plastomes of its closely related subfamily Papaveroideae and other Ranunculales taxa are highly conserved. On the other hand, the synapomorphy characteristics of both accD and the ndh gene loss events happened in the common ancestor of the Corydalis and sub-clade of the Corydalis lineage, respectively. The Corydalis-sub clade species (ndh lost) are distributed predominantly in the Qinghai-Tibetan plateau (QTP) region. The phylogenetic analysis and divergence time estimation were also employed for the Corydalis species. The divergence time of the ndh gene in the Corydalis sub-clade species (44.31 - 15.71 mya) coincides very well with the uplift of the Qinghai-Tibet Plateau in Oligocene and Miocene periods, and maybe during this period, it has probably triggered the radiation of the Corydalis species. To the best of the authors' knowledge, this is the first large-scale comparative study of Corydalis plastomes and their evolution. The present study may provide insights into the plastome architecture and the molecular evolution of Corydalis species.

Comparative chloroplast genomics provides insights into the genealogical relationships of endangered Tetraena mongolica and the chloroplast genome evolution of related Zygophyllaceae species.

A comprehensive understanding of genetic background for rare species will provide an important theoretical basis for the future species management, monitoring and conservation. Tetraena mongolica is restrictedly distributed in the western Ordos plateau of China and has been listed as a national protected plant. We generated 13 chloroplast (cp) genomes of T. mongolica (size range of 106,062-106,230bp) and conducted a series of comparative analyses of six Zygophyllaceae cp genomes. T. mongolica cp genome exhibited a quadripartite structure with drastically reduced inverted repeats (IRs, 4,315bp) and undergone the loss of a suit of ndh genes and a copy of rRNAs. Furthermore, all the T. mongolica populations were divided into two genetic groups based on complete cp phylogenomics. In addition, notably variable genome size, gene order and structural changes had been observed among the six Zygophyllaceae cp genomes. Overall, our findings provide insights into the cp genome evolution mode and intraspecific relationships of T. mongolica, and provide a molecular basis for scientific conservation of this endangered plant.