Plastome phylogenomic study of Gentianeae (Gentianaceae): widespread gene tree discordance and its association with evolutionary rate heterogeneity of plastid genes

Phylogenomics and adaptive evolution of hydrophytic umbellifers (tribe Oenantheae, Apioideae) revealed from chloroplast genomes

BACKGROUND

Tribe Oenantheae consists mainly of aquatic species within the Apioideae. The unique morphology and habitat distinguish this group from other Apioideae groups. However, the genomic information of these group species has not been widely developed, and the molecular mechanisms of adaptive evolution remain unclear.

RESULTS

We provide comparative analyses on 30 chloroplast genomes of this tribe representing five genera to explore the molecular variation response to plant adaptations. The Oenantheae chloroplast genomes presented typical quadripartite structures, with sizes ranging from 153,024 bp to 155,006 bp. Gene content and order were highly conserved with no significant expansion or contraction observed. Seven regions (rps16 intron-trnK, rpoB-trnC, trnE-trnT-psbD, petA-psbJ, ndhF-rpl32-trnL, ycf1a-rps15, and ycf1a gene) were identified as remarkable candidate DNA markers for future studies on species identification, biogeography, and phylogeny of tribe Oenantheae. Our study elucidated the relationships among the genera of tribe Oenantheae and subdivided the genera of Sium and Oenanthe. However, relationships among the Oenanthe I clade remain to be further clarified. Eight positively selected genes (accD, rbcL, rps8, ycf1a, ycf1b, ycf2, ndhF, and ndhK) were persuasively detected under site models tests, and these genes might have played roles in Oenantheae species adaptation to the aquatic environments.

CONCLUSIONS

Our results provide sufficient molecular markers for the subsequent molecular studies of the tribe Oenantheae, and promote the understanding of the adaptation of the Oenantheae species to aquatic environments.

Phylogenomic insights into Neotropical Magnolia relationships

Despite extensive research into the phylogenetic relationships of the genus Magnolia, Neotropical taxa have been neglected. This is partly because their numbers have recently doubled and now account for almost half of the global richness. Therefore, by sampling one-third of all Neotropical taxa their relationships were studied using morphological, nuclear, and plastome data. Two major clades were identified: Clade I, comprising Magnolia sect. Talauma, Magnolia sect. Splendentes from the Neotropics, and the Asian Magnolia sect. Gwillimia; while Clade II included the Neotropical clades Magnolia sect. Macrophylla and Magnolia sect. Magnolia, along with the remaining non-Neotropical sections. Within Clade I, Magnolia sect. Talauma was geographically divided into a northern subclade grouping Mexican and Central American taxa, and a southern subclade comprising South American and Caribbean taxa. Magnolia sect. Splendentes was also dichotomously divided, corresponding to the former Magnolia sect. Talauma subsect. Cubenses and Magnolia sect. Talauma subsect. Dugandiodendron. In Clade II, the relationships within Magnolia sect. Macrophylla and Magnolia sect. Magnolia were unclear, suggesting a species complex in all Magnolia sect. Macrophylla taxa. In total, 25 morphological traits were assessed, and ancestral state reconstructions were carried out. Only the joined clustering of mature follicles was a synapomorphy for the southern subclade of Magnolia sect. Talauma. In conclusion, this highlights the need to re-assess the taxonomic delimitation of certain groups, to update the infrageneric classification of Neotropical clades and to explore morphological traits to support them.

Phylogenomics analysis of Scutellaria (Lamiaceae) of the world

BACKGROUND

Scutellaria, a sub-cosmopolitan genus, stands as one of the Lamiaceae family's largest genera, encompassing approximately 500 species found in both temperate and tropical montane regions. Recognized for its significant medicinal properties, this genus has garnered attention as a research focus, showcasing anti-cancer, anti-inflammatory, antioxidant, and hepatoprotective qualities. Additionally, it finds application in agriculture and horticulture. Comprehending Scutellaria's taxonomy is pivotal for its effective utilization and conservation. However, the current taxonomic frameworks, primarily based on morphological characteristics, are inadequate. Despite several phylogenetic studies, the species relationships and delimitations remain ambiguous, leaving the genus without a stable and reliable classification system.

RESULTS

This study analyzed 234 complete chloroplast genomes, comprising 220 new and 14 previously published sequences across 206 species, subspecies, and varieties worldwide. Phylogenetic analysis was conducted using six data matrices through Maximum Likelihood and Bayesian Inference, resulting in a robustly supported phylogenetic framework for Scutellaria. We propose three subgenera, recommending the elevation of Section Anaspis to subgeneric rank and the merging of Sections Lupulinaria and Apeltanthus. The circumscription of Subgenus Apeltanthus and Section Perilomia needs to be reconsidered. Comparative analysis of chloroplast genomes highlighted the IR/SC boundary feature as a significant taxonomic indicator. We identified a total of 758 SSRs, 558 longer repetitive sequences, and ten highly variable regions, including trnK-rps16, trnC-petN, petN-psbM, accD-psaI, petA-psbJ, rpl32-trnL, ccsA-ndhD, rps15-ycf1, ndhF, and ycf1. These findings serve as valuable references for future research on species identification, phylogeny, and population genetics.

CONCLUSIONS

The phylogeny of Scutellaria, based on the most comprehensive sample collection to date and complete chloroplast genome analysis, has significantly enhanced our understanding of its infrageneric relationships. The extensive examination of chloroplast genome characteristics establishes a solid foundation for the future development and utilization of Scutellaria, an important medicinal plant globally.

Polyphyly of Boehmeria (Urticaceae) congruent with plastome structural variation

Boehmeria is a taxonomically challenging group within the nettle family (Urticaceae). The polyphyly of the genus has been proposed by previous studies with respect to five genera (Debregeasia, Cypholophus, Sarcochlamys, Archiboehmeria, and Astrothalamus). Extensive homoplasy of morphological characters has made generic delimitation problematic. Previous studies in other plant groups suggest that plastome structural variations have the potential to provide characters useful in reconstructing evolutionary relationships. We aimed to test this across Boehmeria and its allied genera by mapping plastome structural variations onto a resolved strongly supported phylogeny. In doing so, we expanded the sampling of the plastome to include Cypholophus, Sarcochlamys, Archiboehmeria, and Astrothalamus for the first time. The results of our phylogenomic analyses provide strong support for Sarcochlamys as being more closely related to Leucosyke puya than to Boehmeria and for the clustering of Boehmeria s.l. into four subclades. The sizes of the plastomes in Boehmeria s.l. ranged from 142,627 bp to 170,958 bp. The plastomes recovered a typical quadripartite structure comprising 127~146 genes. We observe several obvious structural variations across the taxa such as gene loss and multiple gene duplication, inverted repeat (IR) contraction and wide expansions, and inversions. Moreover, we recover a trend for these variations that the early clades were relatively conserved in evolution, whereas the later diverging clades were variable. We propose that the structural variations documented may be linked to the adaptation of Boehmeria s.l. to a wide range of habitats, from moist broadleaf forests in Asia to xeric shrublands and deserts in Africa. This study confirms that variation in plastome gene loss/duplication, IR contraction/expansion, and inversions can provide evidence useful for the reconstruction of evolutionary relationships.

Comparative and phylogenetic analyses of Loranthaceae plastomes provide insights into the evolutionary trajectories of plastome degradation in hemiparasitic plants

BACKGROUND

The lifestyle transition from autotrophy to heterotrophy often leads to extensive degradation of plastomes in parasitic plants, while the evolutionary trajectories of plastome degradation associated with parasitism in hemiparasitic plants remain poorly understood. In this study, phylogeny-oriented comparative analyses were conducted to investigate whether obligate Loranthaceae stem-parasites experienced higher degrees of plastome degradation than closely related facultative root-parasites and to explore the potential evolutionary events that triggered the 'domino effect' in plastome degradation of hemiparasitic plants.

RESULTS

Through phylogeny-oriented comparative analyses, the results indicate that Loranthaceae hemiparasites have undergone varying degrees of plastome degradation as they evolved towards a heterotrophic lifestyle. Compared to closely related facultative root-parasites, all obligate stem-parasites exhibited an elevated degree plastome degradation, characterized by increased downsizing, gene loss, and pseudogenization, thereby providing empirical evidence supporting the theoretical expectation that evolution from facultative parasitism to obligate parasitism may result in a higher degree of plastome degradation in hemiparasites. Along with infra-familial divergence in Loranthaceae, several lineage-specific gene loss/pseudogenization events occurred at deep nodes, whereas further independent gene loss/pseudogenization events were observed in shallow branches.

CONCLUSIONS

The findings suggest that in addition to the increasing levels of nutritional reliance on host plants, cladogenesis can be considered as another pivotal evolutionary event triggering the 'domino effect' in plastome degradation of hemiparasitic plants. These findings provide new insights into the evolutionary trajectory of plastome degradation in hemiparasitic plants.

The complete chloroplast genome and phylogentic results support the species position of Swertia banzragczii and Swertia marginata (Gentianaceae) in Mongolia

BACKGROUND

Swertia banzragczii and S. marginata are important medicinal species in Mongolia. However, their taxonomic positions and genetic backgrounds remain unknown. In this study, we explored the complete chloroplast genomes and DNA barcoding of these species and compared them with those of closely related species within the subgenus to determine their taxonomic positions and phylogenetic relationships.

RESULT

The chloroplast genomes of S. banzragczii and S. marginata encoded 114 genes, including 80 protein-coding genes, 30 tRNA genes, and 4 rRNA genes. Among them, 16 genes contained a single intron, and 2 genes had two introns. Closely related species had a conserved genome structure and gene content. Only differences in genome length were noticed, which were caused by the expansion and contraction of the inverted repeat (IR) region and loss of exons in some genes. The trnH-GUG-psbA and trnD-GUC-trnY-GUA intergenic regions had high genetic diversity within Swertia plastomes. Overall, S. banzragczii and S. marginata are true species and belong to the subgenus Swertia.

CONCLUSIONS

These results provide valuable genetic and morphological information on rare and subendemic Swertia species in Mongolia, which can be used for further advanced studies on the Swertia genus.

Species delimitation of tea plants (Camellia sect. Thea) based on super-barcodes

BACKGROUND

The era of high throughput sequencing offers new paths to identifying species boundaries that are complementary to traditional morphology-based delimitations. De novo species delimitation using traditional or DNA super-barcodes serve as efficient approaches to recognizing putative species (molecular operational taxonomic units, MOTUs). Tea plants (Camellia sect. Thea) form a group of morphologically similar species with significant economic value, providing the raw material for tea, which is the most popular nonalcoholic caffeine-containing beverage in the world. Taxonomic challenges have arisen from vague species boundaries in this group.

RESULTS

Based on the most comprehensive sampling of C. sect. Thea by far (165 individuals of 39 morphospecies), we applied three de novo species delimitation methods (ASAP, PTP, and mPTP) using plastome data to provide an independent evaluation of morphology-based species boundaries in tea plants. Comparing MOTU partitions with morphospecies, we particularly tested the congruence of MOTUs resulting from different methods. We recognized 28 consensus MOTUs within C. sect. Thea, while tentatively suggesting that 11 morphospecies be discarded. Ten of the 28 consensus MOTUs were uncovered as morphospecies complexes in need of further study integrating other evidence. Our results also showed a strong imbalance among the analyzed MOTUs in terms of the number of molecular diagnostic characters.

CONCLUSION

This study serves as a solid step forward for recognizing the underlying species boundaries of tea plants, providing a needed evidence-based framework for the utilization and conservation of this economically important plant group.

Comparative Analysis of Six Complete Plastomes of Tripterospermum spp

The Tripterospermum, comprising 34 species, is a genus of Gentianaceae. Members of Tripterospermum are mostly perennial, entwined herbs with high medicinal value and rich in iridoids, xanthones, flavonoids, and triterpenes. However, our inadequate understanding of the differences in the plastid genome sequences of Tripterospermum species has severely hindered the study of their evolution and phylogeny. Therefore, we first analyzed the 86 Gentianae plastid genomes to explore the phylogenetic relationships within the Gentianae subfamily where Tripterospermum is located. Then, we analyzed six plastid genomes of Tripterospermum, including two newly sequenced plastid genomes and four previously published plastid genomes, to explore the plastid genomes' evolution and phylogenetic relationships in the genus Tripterospermum. The Tripterospermum plastomes have a quadripartite structure and are between 150,929 and 151,350 bp in size. The plastomes of Tripterospermum encoding 134 genes were detected, including 86 protein-coding genes (CDS), 37 transfer RNA (tRNA) genes, eight ribosomal RNA (rRNA) genes, and three pseudogenes (infA, rps19, and ycf1). The result of the comparison shows that the Tripterospermum plastomes are very conserved, with the total plastome GC content ranging from 37.70% to 37.79%. In repeat sequence analysis, the number of single nucleotide repeats (A/T) varies among the six Tripterospermum species, and the identified main long repeat types are forward and palindromic repeats. The degree of conservation is higher at the SC/IR boundary. The regions with the highest divergence in the CDS and the intergenic region (IGS) are psaI and rrn4.5-rrn5, respectively. The average pi of the CDS and the IGS are only 0.071% and 0.232%, respectively, indicating that the Tripterospermum plastomes are highly conserved. Phylogenetic analysis indicated that Gentianinae is divided into two clades, with Tripterospermum as a sister to Sinogeniana. Phylogenetic trees based on CDS and CDS + IGS combined matrices have strong support in Tripterospermum. These findings contribute to the elucidation of the plastid genome evolution of Tripterospermum and provide a foundation for further exploration and resource utilization within this genus.

A new perspective on codon usage, selective pressure, and phylogenetic implications of the plastomes in the Telephium clade (Crassulaceae)

The Telephium clade of the Crassulaceae family contains many medicinal, ornamental, and ecologically restorative plants. However, the phylogenetic relationships within the clade remain debated, and comprehensive analyses of codon usage and selection pressure in Telephium plastomes are limited. In this study, we assembled and annotated four plastomes and performed extensive analyses. The plastomes exhibited a typical quadripartite structure and high conservation. The lengths ranged from 151,357 bp to 151,641 bp with 134 genes identified. The GC content was the highest within IR, followed by LSC, and lowest in the SSC region. Meanwhile, a unique inversion was observed within the LSC region of Meterostachys sikokianus. Polymorphisms analysis revealed minimum nucleotide diversity in the IR regions, with over ten highly polymorphic regions identified. Phylogenetically, two subclades formed within the monophyletic Telephium clade, with Umbilicus as the sister group to the remaining Hylotelephium subclade members. Notably, no significant positive selection was found among the 79 plastid genes, which showed varying evolutionary patterns. However, 19 genes contained codons under positive selection. The specific functions of these sites require further investigation. Synonymous codon usage was biased and conserved across the tested plastomes, shaped by natural selection, mutations and other factors of varying influence. We also identified 34 taxon-specific codon aversion motifs from 49 plastid genes. Our plastomic analyses elucidate phylogenetic relationships and evolutionary patterns in this medicinal clade, providing a foundation for further research on these ecologically and pharmaceutically important plants.

Evolution of the Araliaceae family involved rapid diversification of the Asian Palmate group and Hydrocotyle specific mutational pressure

The Araliaceae contain many valuable species in medicinal and industrial aspects. We performed intensive phylogenomics using the plastid genome (plastome) and 45S nuclear ribosomal DNA sequences. A total of 66 plastome sequences were used, 13 of which were newly assembled in this study, 12 from new sequences, and one from existing data. While Araliaceae plastomes showed conserved genome structure, phylogenetic reconstructions based on four different plastome datasets revealed phylogenetic discordance within the Asian Palmate group. The divergence time estimation revealed that splits in two Araliaceae subfamilies and the clades exhibiting phylogenetic discordances in the Asian Palmate group occurred at two climatic optima, suggesting that global warming events triggered species divergence, particularly the rapid diversification of the Asian Palmate group during the Middle Miocene. Nucleotide substitution analyses indicated that the Hydrocotyloideae plastomes have undergone accelerated AT-biased mutations (C-to-T transitions) compared with the Aralioideae plastomes, and the acceleration may occur in their mitochondrial and nuclear genomes as well. This implies that members of the genus Hydrocotyle, the only aquatic plants in the Araliaceae, have experienced a distinct evolutionary history from the other species. We also discussed the intercontinental disjunction in the genus Panax and proposed a hypothesis to complement the previously proposed hypothesis. Our results provide the evolutionary trajectory of Araliaceae and advance our current understanding of the evolution of Araliaceae species.

Insights into the phylogenetic relationships and species boundaries of the Myricaria squamosa complex (Tamaricaceae) based on the complete chloroplast genome

Myricaria plants are widely distributed in Eurasia and are helpful for windbreak and embankment protection. Current molecular evidence has led to controversy regarding species boundaries within the Myricaria genus and interspecific phylogenetic relationships between three specific species-M. bracteata, M. paniculata and M. squamosa-which have remained unresolved. This study treated these three unresolved taxa as a species complex, named the M. squamosa complex. The genome skimming approach was used to determine 35 complete plastome sequences and nuclear ribosomal DNA sequences for the said complex and other closely related species, followed by de novo assembly. Comparative analyses were conducted across Myricaria to identify the genome size, gene content, repeat type and number, SSR (simple sequence repeat) abundance, and codon usage bias of chloroplast genomes. Tree-based species delimitation results indicated that M. bracteata, M. paniculata and M. squamosa could not be distinguished and formed two monophyletic lineages (P1 and P2) that were clustered together. Compared to plastome-based species delimitation, the standard nuclear DNA barcode had the lowest species resolution, and the standard chloroplast DNA barcode and group-specific barcodes delimitated a maximum of four out of the five species. Plastid phylogenomics analyses indicated that the monophyletic M. squamosa complex is comprised of two evolutionarily significant units: one in the western Tarim Basin and the other in the eastern Qinghai-Tibet Plateau. This finding contradicts previous species discrimination and promotes the urgent need for taxonomic revision of the threatened genus Myricaria. Dense sampling and plastid genomes will be essential in this effort. The super-barcodes and specific barcode candidates outlined in this study will aid in further studies of evolutionary history.

More than a spiny morphology: plastome variation in the prickly pear cacti (Opuntieae)

BACKGROUND

Plastid genomes (plastomes) have long been recognized as highly conserved in their overall structure, size, gene arrangement and content among land plants. However, recent studies have shown that some lineages present unusual variations in some of these features. Members of the cactus family are one of these lineages, with distinct plastome structures reported across disparate lineages, including gene losses, inversions, boundary movements or loss of the canonical inverted repeat (IR) region. However, only a small fraction of cactus diversity has been analysed so far.

METHODS

Here, we investigated plastome features of the tribe Opuntieae, the remarkable prickly pear cacti, which represent one of the most diverse and important lineages of Cactaceae. We assembled de novo the plastome of 43 species, representing a comprehensive sampling of the tribe, including all seven genera, and analysed their evolution in a phylogenetic comparative framework. Phylogenomic analyses with different datasets (full plastome sequences and genes only) were performed, followed by congruence analyses to assess signals underlying contentious nodes.

KEY RESULTS

Plastomes varied considerably in length, from 121 to 162 kbp, with striking differences in the content and size of the IR region (contraction and expansion events), including a lack of the canonical IR in some lineages and the pseudogenization or loss of some genes. Overall, nine different types of plastomes were reported, deviating in the presence of the IR region or the genes contained in the IR. Overall, plastome sequences resolved phylogenetic relationships within major clades of Opuntieae with high bootstrap values but presented some contentious nodes depending on the dataset analysed (e.g. whole plastome vs. genes only). Congruence analyses revealed that most plastidial regions lack phylogenetic resolution, while few markers are supporting the most likely topology. Likewise, alternative topologies are driven by a handful of plastome markers, suggesting recalcitrant nodes in the phylogeny.

CONCLUSIONS

Our study reveals a dynamic nature of plastome evolution across closely related lineages, shedding light on peculiar features of plastomes. Variation of plastome types across Opuntieae is remarkable in size, structure and content and can be important for the recognition of species in some major clades. Unravelling connections between the causes of plastome variation and the consequences for species biology, physiology, ecology, diversification and adaptation is a promising and ambitious endeavour in cactus research. Although plastome data resolved major phylogenetic relationships, the generation of nuclear genomic data is necessary to confront these hypotheses and assess the recalcitrant nodes further.

On Pattern-Cladistic Analyses Based on Complete Plastid Genome Sequences

The fundamental Hennigian principle, grouping solely on synapomorphy, is seldom used in modern phylogenetics. In the submitted paper, we apply this principle in reanalyzing five datasets comprising 197 complete plastid genomes (plastomes). We focused on the latter because plastome-based DNA sequence data gained dramatic popularity in molecular systematics during the last decade. We show that pattern-cladistic analyses based on complete plastid genome sequences can successfully resolve affinities between plant taxa, simultaneously simplifying both the genomic and analytical frameworks of phylogenetic studies. We developed "Matrix to Newick" (M2N), a program to represent the standard molecular alignment of plastid genomes in the form of trees or relationships directly. Thus, massive plastome-based DNA sequence data can be successfully represented in a relational form rather than as a standard molecular alignment. Application of methods of median supertree construction (the Average Consensus method has been used as an example in this study) or Maximum Parsimony analysis to relational representations of plastome sequence data may help systematist to avoid the complicated assumption-based frameworks of Maximum Likelihood or Bayesian phylogenetics that are most used today in massive plastid sequence data analyses. We also found that significant amounts of pure genomic information that typically accommodate the majority of current plastid phylogenomic studies can be effectively dropped by systematists if they focus on the pattern-cladistics or relational analyses of plastome-based molecular data. The proposed pattern-cladistic approach is a powerful and straightforward heuristic alternative to modern plastome-based phylogenetics.

The complete chloroplast genome sequence of Swertia japonica (Schult.) Makino (Gentianaceae)

Swertia japonica (Schult.) Makino is a traditional medicinal plant in Japan for which the chloroplast genome has not been previously reported. The complete chloroplast genome of S. japonica was determined using a high-throughput sequencing technique. The total length of the S. japonica chloroplast genome was 153,208 bp, and comprised a large single-copy region of 83,319 bp, and a small single-copy region of 18,375 bp, separated by a pair of 25,757 bp inverted repeat regions. A phylogenetic analysis, based on the obtained chloroplast genome, indicated that S. japonica is closely related to S. diluta, S. franchetiana, S. kouitchensis, S. mussotii, and S. punicea. The presented chloroplast genome will be useful for further taxonomic, pharmacological and evolutionary studies of Swertia.

Characteristics of the complete chloroplast genome of Swertia divaricata Harry Sm. (Gentianaceae) and its phylogenetic inference

Swertia divaricata Harry Sm., 1965, (Gentianaceae) is a perennial herb endemic to Northwest Yunnan, China, belonging to the species-rich genus Swertia. It possesses unique morphological features but its systematic position remains uncertain. To determine its phylogenetic placement, the complete plastid genome of S. divaricata was assembled utilizing high-throughput sequencing data. The genome is circular, spanning 152,073 bp, and comprises a large single-copy (LSC) region of 82,470 bp, a small single-copy (SSC) region of 18,153 bp, and two inverted repeats (IR) regions, each 25,725 bp. A total of 130 genes were annotated, including 85 protein-coding genes, 37 tRNA genes, and eight rRNA genes. The plastome of S. divaricata exhibits a structure and gene composition highly similar to those of other Swertia plastomes. Phylogenetic analysis indicated that S. divaricata is closely related to S. erythrosticta, sister to a subclade comprising species from sections Swertia and Apterae. The plastome sequence described herein constitutes a valuable contribution to phylogenetic and evolutionary research on Swertia.

Phylogeny of Leontopodium (Asteraceae) in China-with a reference to plastid genome and nuclear ribosomal DNA

The infrageneric taxonomy system, species delimitation, and interspecies systematic relationships of Leontopodium remain controversial and complex. However, only a few studies have focused on the molecular phylogeny of this genus. In this study, the characteristics of 43 chloroplast genomes of Leontopodium and its closely related genera were analyzed. Phylogenetic relationships were inferred based on chloroplast genomes and nuclear ribosomal DNA (nrDNA). Finally, together with the morphological characteristics, the relationships within Leontopodium were identified and discussed. The results showed that the chloroplast genomes of Filago, Gamochaeta, and Leontopodium were well-conserved in terms of gene number, gene order, and GC content. The most remarkable differences among the three genera were the length of the complete chloroplast genome, large single-copy region, small single-copy region, and inverted repeat region. In addition, the chloroplast genome structure of Leontopodium exhibited high consistency and was obviously different from that of Filago and Gamochaeta in some regions, such as matk, trnK (UUU)-rps16, petN-psbM, and trnE (UUC)-rpoB. All the phylogenetic trees indicated that Leontopodium was monophyletic. Except for the subgeneric level, our molecular phylogenetic results were inconsistent with the previous taxonomic system, which was based on morphological characteristics. Nevertheless, we found that the characteristics of the leaf base, stem types, and carpopodium base were phylogenetically correlated and may have potential value in the taxonomic study of Leontopodium. In the phylogenetic trees inferred using complete chloroplast genomes, the subgen. Leontopodium was divided into two clades (Clades 1 and 2), with most species in Clade 1 having herbaceous stems, amplexicaul, or sheathed leaves, and constricted carpopodium; most species in Clade 2 had woody stems, not amplexicaul and sheathed leaves, and not constricted carpopodium.

Phylogeny of Trigonotis in China-with a special reference to its nutlet morphology and plastid genome

The genus Trigonotis comprises nearly 60 species mainly distributed in East and Southeast Asia. China has the largest number of Trigonotis species in the world, with a total of 44 species, of which 38 are endemic. Nutlet morphology is useful for the taxonomic delimitation of Trigonotis. However, there are still controversial circumscriptions of nutlet shape in some species. In previous studies, interspecies phylogenetic relationships were inferred using few DNA markers and very few taxa, which possibly led to erroneous or incomplete conclusions. In this study, the nutlet morphology of 39 Trigonotis taxa and the characteristics of 34 complete chloroplast genomes (29 taxa) were investigated and analyzed. Then, the phylogenetic relationships were discussed within this genus based on complete chloroplast genomes. To the best of our knowledge, this study is the first comprehensive analysis of nutlet morphology and complete chloroplast genome of Trigonotis. Based on nutlet morphology, Trigonotis can be divided into two groups: Group 1, hemispherical or oblique tetrahedron with carpopodiums, and Group 2, inverted tetrahedron without carpopodiums. The chloroplast genome of Trigonotis exhibited a typical quadripartite structure, including 84-86 protein-coding, 37 transfer RNA, and 8 ribosomal RNA genes, with a total length of 147,247-148,986 bp. Genes in the junctions were well conserved in Trigonotis, similar to those in other Boraginaceae s.str. species. Furthermore, Trigonotis chloroplast genomes showed relatively high diversity, with more conserved genic regions than intergenic regions; in addition, we detected 14 hot spots (Pi > 0.005) in non-coding regions. Phylogenetic analyses based on chloroplast genome data identified highly resolved relationships between Trigonotis species. Specifically, Trigonotis was divided into two clades with strong support: one clade included species with hemispherical or oblique tetrahedron nutlets with carpopodiums and bracts, whereas the other clade included species with inverted tetrahedron nutlets without carpopodiums or bracts. Our results may inform future taxonomic, phylogenetic, and evolutionary studies on Boraginaceae.

A comprehensive metabolic map reveals major quality regulations in red-flesh kiwifruit (Actinidia chinensis)

Kiwifruit (Actinidia chinensis) is one of the popular fruits world-wide, and its quality is mainly determined by key metabolites (sugars, flavonoids, and vitamins). Previous works on kiwifruit are mostly done via a single omics approach or involve only limited metabolites. Consequently, the dynamic metabolomes during kiwifruit development and ripening and the underlying regulatory mechanisms are poorly understood. In this study, using high-resolution metabolomic and transcriptomic analyses, we investigated kiwifruit metabolic landscapes at 11 different developmental and ripening stages and revealed a parallel classification of 515 metabolites and their co-expressed genes into 10 distinct metabolic vs gene modules (MM vs GM). Through integrative bioinformatics coupled with functional genomic assays, we constructed a global map and uncovered essential transcriptomic and transcriptional regulatory networks for all major metabolic changes that occurred throughout the kiwifruit growth cycle. Apart from known MM vs GM for metabolites such as soluble sugars, we identified novel transcription factors that regulate the accumulation of procyanidins, vitamin C, and other important metabolites. Our findings thus shed light on the kiwifruit metabolic regulatory network and provide a valuable resource for the designed improvement of kiwifruit quality.

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.

Untying the Gordian knot of plastid phylogenomic conflict: A case from ferns

Phylogenomic studies based on plastid genome have resolved recalcitrant relationships among various plants, yet the phylogeny of Dennstaedtiaceae at the level of family and genera remains unresolved due to conflicting plastid genes, limited molecular data and incomplete taxon sampling of previous studies. The present study generated 30 new plastid genomes of Dennstaedtiaceae (9 genera, 29 species), which were combined with 42 publicly available plastid genomes (including 24 families, 27 genera, 42 species) to explore the evolution of Dennstaedtiaceae. In order to minimize the impact of systematic errors on the resolution of phylogenetic inference, we applied six strategies to generate 30 datasets based on CDS, intergenic spacers, and whole plastome, and two tree inference methods (maximum-likelihood, ML; and multispecies coalescent, MSC) to comprehensively analyze the plastome-scale data. Besides, the phylogenetic signal among all loci was quantified for controversial nodes using ML framework, and different topologies hypotheses among all datasets were tested. The species trees based on different datasets and methods revealed obvious conflicts at the base of the polypody ferns. The topology of the "CDS-codon-align-rm3" (CDS with the removal of the third codon) matrix was selected as the primary reference or summary tree. The final phylogenetic tree supported Dennstaedtiaceae as the sister group to eupolypods, and Dennstaedtioideae was divided into four clades with full support. This robust reconstructed phylogenetic backbone establishes a framework for future studies on Dennstaedtiaceae classification, evolution and diversification. The present study suggests considering plastid phylogenomic conflict when using plastid genomes. From our results, reducing saturated genes or sites can effectively mitigate tree conflicts for distantly related taxa. Moreover, phylogenetic trees based on amino acid sequences can be used as a comparison to verify the confidence of nucleotide-based trees.

Strong plastid degradation is consistent within section Chondrophyllae, the most speciose lineage of Gentiana

Recovering phylogenetic relationships in lineages experiencing intense diversification has always been a persistent challenge in evolutionary studies, including in Gentiana section Chondrophyllae sensu lato (s.l.). Indeed, this subcosmopolitan taxon encompasses more than 180 mostly annual species distributed around the world. We sequenced and assembled 22 new plastomes representing 21 species in section Chondrophyllae s.l. In addition to previously released plastome data, our study includes all main lineages within the section. We reconstructed their phylogenetic relationships based on protein-coding genes and recombinant DNA (rDNA) cistron sequences, and then investigated plastome structural evolution as well as divergence time. Despite an admittedly humble species cover overall, we recovered a well-supported phylogenetic tree based on plastome data, and found significant discordance between phylogenetic relationships and taxonomic treatments. Our results show that G. capitata and G. leucomelaena diverged early within the section, which is then further divided into two clades. The divergence time estimation showed that section Chondrophyllae s.l. evolved in the second half of the Oligocene. We found that section Chondrophyllae s.l. had the smallest average plastome size (128 KB) in tribe Gentianeae (Gentianaceae), with frequent gene and sequence losses such as the ndh complex and its flanking regions. In addition, we detected both expansion and contraction of the inverted repeat (IR) regions. Our study suggests that plastome degradation parallels the diversification of this group, and illustrates the strong discordance between phylogenetic relationships and taxonomic treatments, which now need to be carefully revised.

Comparative Chloroplast Genomics of Seven Endangered Cypripedium Species and Phylogenetic Relationships of Orchidaceae

The species in the genus Cypripedium (Orchidaceae) are considered endangered, mainly distributed in the temperate regions of the Northern Hemisphere, with high ornamental and economic value. Despite previous extensive studies based on both morphology and molecular data, species and sections relationships within Cypripedium remain controversial. Here, we employed two newly generated Cypripedium chloroplast genomes with five other published genomes to elucidate their genomic characteristics. The two genomes were 162,773-207,142 bp in length and contained 128-130 genes, including 82-84 protein-coding genes, 38 tRNA genes, and 8 rRNA genes. We identified 2,192 simple sequence repeats, 786 large repeat sequences, and 7,929 variable loci. The increase of repeat sequences (simple sequence repeats and large repeat sequences) causes a significant amplification in the chloroplast genome size of Cypripedium. The expansion of the IR region led to the pseudogenization or loss of genes in the SSC region. In addition, we identified 12 highly polymorphic loci (Pi > 0.09) suitable for inferring the phylogeny of Cypripedium species. Based on data sets of whole chloroplast genomes (IRa excluded) and protein-coding sequences, a well-supported phylogenetic tree was reconstructed, strongly supporting the five subfamilies of Orchidaceae and the genus Cypripedium as monophyletic taxa. Our findings also supported that C. palangshanense belonged to sect. Palangshanensia rather than sect. Retinervia. This study also enriched the genomic resources of Cypripedium, which may help to promote the conservation efforts of these endangered species.

Plastid phylogenomic analyses of the Selaginella sanguinolenta group (Selaginellaceae) reveal conflict signatures resulting from sequence types, outlier genes, and pervasive RNA editing

Different from the generally conserved plastomes (plastid genomes) of most land plants, the Selaginellaceae plastomes exhibit dynamic structure, high GC content and high substitution rates. Previous plastome analyses identified strong conflict on several clades in Selaginella, however the factors causing the conflictions and the impact on the phylogenetic inference have not been sufficiently investigated. Here, we dissect the distribution of phylogenetic signals and conflicts in Selaginella sanguinolenta group, the plastome of which is DR (direct repeats) structure and with genome-wide RNA editing. We analyzed the data sets including 22 plastomes representing all species of the S. sanguinolenta group, covering the entire geographical distribution from the Himalayas to Siberia and the Russian Far East regions. We recovered four different topologies by applying multispecies coalescent (ASTRAL) and concatenation methods (IQ-TREE and RAxML) on four data sets of PC (protein-coding genes), NC (non-coding sequences), PCN (the concatenated PC and NC), and RC (predicted RNA editing sites "C" were corrected by "T"), respectively. Six monophyletic clades, S. nummularifolia clade, S. rossii clade, S. sajanensis clade, S. sanguinolenta I clade, S. sanguinolenta II clade, and S. sanguinolenta III clade, were consistently resolved and supported by the characteristics of GC content, RNA editing frequency, and gene content. However, the relationships among these clades varied across the four topologies. To explore the underlying causes of the uncertainty, we compared the phylogenetic signals of the four topologies. We identified that the sequence types (coding versus non-coding), outlier genes (genes with extremely high |ΔGLS| values), and C-to-U RNA editing frequency in the protein-coding genes were responsible for the unstable phylogenomic relationship. We further revealed a significant positive correlation between the |ΔGLS| values and the variation coefficient of the RNA editing number. Our results demonstrated that the coalescent method performed better than the concatenation method in overcoming the problems caused by outlier genes and extreme RNA editing events. Our study particularly focused on the importance of exploring the plastid phylogenomic conflicts and suggested conducting concatenated analyses cautiously when adopting organelle genome data.

Complete chloroplast genome of Gentianopsis barbata and comparative analysis with related species from Gentianaceae

Gentianopsis barbata is an essential medicinal plant in China with high ornamental and medicinal values. Unfortunately, the study of the chloroplast genome of this plant still has a gap. This study sequenced and characterized the complete chloroplast genome of G. barbata. The complete chloroplast genome of G. barbata is a typical circular structure with 151,123 bp. It consists of a large single-copy region (82,690 bp) and a small single-copy region (17,887 bp) separated by a pair of inverted repeats (25,273 bp), which covers 78 protein-coding genes, 30 tRNAs, and 4 rRNAs. Repeat analysis showed the highest frequency of palindrome. Thirty-seven simple sequence repeats were identified, most of which were single nucleotides. The bayesian inference tree, maximum likelihood tree, and neighbor joining tree suggested that G. barbata is grouped with Gentianopsis grandis and Gentianopsis paludosa. The divergence time analysis showed that G. barbata diverged at 1.243 Mya. Comparative chloroplast analysis can reveal interspecific diversity, and regions with high variation can be used to develop molecular markers applicable to various research areas. Our results provide new insight into plastome evolution and valuable resource for further studies on G. barbata.

Multiple lines of evidence supports the two varieties of Halenia elliptica (Gentianaceae) as two species

Delimiting species requires multiple sources of evidence. Here, we delimited two varieties of Halenia elliptica (Gentianaceae) using several lines of evidence, including morphological traits and mating system in a sympatric population, phylogenetic relationships based on nrITS and cpDNA (rpl16) data, and complete chloroplast genome sequences. Comparative analysis of 21 morphological traits clearly separates the two varieties of H. elliptica. Examination of the flowering process and pollination treatments indicate that H. elliptica var. grandiflora produces seeds via outcrossing, whereas H. elliptica var. elliptica produces seeds via mixed mating. Furthermore, hand-pollinated hybridization of the two varieties produced no seeds. Observations of pollinators showed that when bees began a pollination bout on H. elliptica var. grandiflora they preferred to continue pollinating this variety; however, when they began a pollination bout on H. elliptica var. elliptica, they showed no preference for either variety. Phylogenetic analysis confirmed the monophyly of H. elliptica, which was further divided into two monophyletic clades corresponding to the two varieties. A large number of variants from the chloroplast genomes reflected remarkable genetic dissimilarities between the two varieties of H. elliptica. We recommend that the two varieties of H. elliptica should be revised as two species (H. elliptica and H. grandiflora). Our findings indicate that H. elliptica varieties may have split into two separate species due to a shift in mating system, changes in flowering phenology and/or post-pollination reproductive isolation.

Plastome structure, phylogenomics and evolution of plastid genes in Swertia (Gentianaceae) in the Qing-Tibetan Plateau

BACKGROUND

The genus Swertia is of great medicinal importance and one of the most taxonomically challenging taxa within Gentianaceae, largely due to the morphological similarities of species within this genus and with its closely related genera. Previous molecular studies confirmed its polyphyly but suffered from low phylogenetic resolutions because only limited sequence loci were used. Thus, we conducted the structural, gene evolutionary, and phylogenetic analyses of 11 newly obtained plastomes of Swertia. Our result greatly improved the phylogenetic resolutions in Swertia, shed new light on the plastome evolution and phylogenetic relationships of this genus.

RESULTS

The 11 Swertia plastomes together with the published seven species proved highly similar in overall size, structure, gene order, and content, but revealed some structural variations caused by the expansion and contraction of the IRb region into the LSC region, due to the heterogeneous length of the ψycf1. The gene rps16 was found to be in a state flux with pseudogenes or completely lost. Similar situation was also documented in other genera of Gentianaceae. This might imply loss of the gene in the common ancestor of Gentianaceae. The distribution plot of ENC vs. GC3 showed all these plastomes arranging very close in the Wright line with an expected ENC value (49-52%), suggesting the codon usage of Swertia was mainly constrained by a GC mutation bias. Most of the genes remained under the purifying selection, however, the cemA was identified under positive selection, possibly reflecting an adaptive response to low CO₂ atmospheric conditions during the Late Miocene. Our phylogenomic analyses, based on 74 protein-coding genes (CDS), supported the polyphyly of Swertia with its close allies in the subtribe Swertiinae, presumably due to recent rapid radiation. The topology inferred from our phylogenetic analyses partly supported the current taxonomic treatment. Finally, several highly variable loci were identified, which can be used in future phylogenetic studies and accurate identification of medicinal genuineness of Swertia.

CONCLUSIONS

Our study confirmed the polyphyly of Swertia and demonstrated the power of plastome phylogenomics in improvement of phylogenetic resolution, thus contributing to a better understanding of the evolutionary history of Swertia.

Comparative plastome analysis of Musaceae and new insights into phylogenetic relationships

Background

Musaceae is an economically important family consisting of 70-80 species. Elucidation of the interspecific relationships of this family is essential for a more efficient conservation and utilization of genetic resources for banana improvement. However, the scarcity of herbarium specimens and quality molecular markers have limited our understanding of the phylogenetic relationships in wild species of Musaceae. Aiming at improving the phylogenetic resolution of Musaceae, we analyzed a comprehensive set of 49 plastomes for 48 species/subspecies representing all three genera of this family.

Results

Musaceae plastomes have a relatively well-conserved genomic size and gene content, with a full length ranging from 166,782 bp to 172,514 bp. Variations in the IR borders were found to show phylogenetic signals to a certain extent in Musa. Codon usage bias analysis showed different preferences for the same codon between species and three genera and a common preference for A/T-ending codons. Among the two genes detected under positive selection (dN/dS > 1), ycf2 was indicated under an intensive positive selection. The divergent hotspot analysis allowed the identification of four regions (ndhF-trnL, ndhF, matK-rps16, and accD) as specific DNA barcodes for Musaceae species.

Bayesian and maximum likelihood phylogenetic analyses using full plastome resulted in nearly identical tree topologies with highly supported relationships between species. The monospecies genus Musella is sister to Ensete, and the genus Musa was divided into two large clades, which corresponded well to the basic number of n = x = 11 and n = x =10/9/7, respectively. Four subclades were divided within the genus Musa. A dating analysis covering the whole Zingiberales indicated that the divergence of Musaceae family originated in the Palaeocene (59.19 Ma), and the genus Musa diverged into two clades in the Eocene (50.70 Ma) and then started to diversify from the late Oligocene (29.92 Ma) to the late Miocene. Two lineages (Rhodochlamys and Australimusa) radiated recently in the Pliocene /Pleistocene periods.

Conclusions

The plastome sequences performed well in resolving the phylogenetic relationships of Musaceae and generated new insights into its evolution. Plastome sequences provided valuable resources for population genetics and phylogenetics at lower taxon.

Supplementary Information

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

Towards the plastome evolution and phylogeny of Cycas L. (Cycadaceae): molecular-morphology discordance and gene tree space analysis

BACKGROUND

Plastid genomes (plastomes) present great potential in resolving multiscale phylogenetic relationship but few studies have focused on the influence of genetic characteristics of plastid genes, such as genetic variation and phylogenetic discordance, in resolving the phylogeny within a lineage. Here we examine plastome characteristics of Cycas L., the most diverse genus among extant cycads, and investigate the deep phylogenetic relationships within Cycas by sampling 47 plastomes representing all major clades from six sections.

RESULTS

All Cycas plastomes shared consistent gene content and structure with only one gene loss detected in Philippine species C. wadei. Three novel plastome regions (psbA-matK, trnN-ndhF, chlL-trnN) were identified as containing the highest nucleotide variability. Molecular evolutionary analysis showed most of the plastid protein-coding genes have been under purifying selection except ndhB. Phylogenomic analyses that alternatively included concatenated and coalescent methods, both identified four clades but with conflicting topologies at shallow nodes. Specifically, we found three species-rich Cycas sections, namely Stangerioides, Indosinenses and Cycas, were not or only weakly supported as monophyly based on plastomic phylogeny. Tree space analyses based on different tree-inference methods both revealed three gene clusters, of which the cluster with moderate genetic properties showed the best congruence with the favored phylogeny.

CONCLUSIONS

Our exploration in plastomic data for Cycas supports the idea that plastid protein-coding genes may exhibit discordance in phylogenetic signals. The incongruence between molecular phylogeny and morphological classification reported here may largely be attributed to the uniparental attribute of plastid, which cannot offer sufficient information to resolve the phylogeny. Contrasting to a previous consensus that genes with longer sequences and a higher proportion of variances are superior for phylogeny reconstruction, our result implies that the most effective phylogenetic signals could come from loci that own moderate variation, GC content, sequence length, and underwent modest selection.

Phylogenetic position of Bupleurum sikangense inferred from the complete chloroplast genome sequence

Bupleurum sikangense is an endemic species to China distributed in Xizang (Tibet), which has high saikosaponin content and potential medicinal value. Morphologically, it extremely resembles B. commelynoideum. In order to get a better understanding of the relationship between B. sikangense and B. commelynoideum, and on the phylogenetic status of the two species in the genus, the complete chloroplast (cp) genomes of them were sequenced. The genome organization, repeat sequences, codon usage, RNA-editing sites, and variation of their cp genomes revealed high similarity between the species. Some highly variable regions like trnK-UUU_rps16, rps16_trnQ-UUG, ndhC_trnV-UAC, petA_psbJ, accD_psaI, and petL_psbE were identified, providing potential molecular markers for differentiating the two species. Phylogenetic analysis indicated that B. commelynoideum has a closer relationship to B. chinese than that to B. sikangense. Overall, this study will not only improve our knowledge about cp genomes of these two species, and but also provide data for further research on species identification, safe medical application, conservation genetics, etc., of Bupleurum plants.

Comparative analysis of chloroplast genomes of cultivars and wild species of sweetpotato (Ipomoea batatas [L.] Lam)

Background

Sweetpotato (Ipomoea batatas [L.] Lam.) is an important food crop. However, the genetic information of the nuclear genome of this species is difficult to determine accurately because of its large genome and complex genetic background. This drawback has limited studies on the origin, evolution, genetic diversity and other relevant studies on sweetpotato.

Results

The chloroplast genomes of 107 sweetpotato cultivars were sequenced, assembled and annotated. The resulting chloroplast genomes were comparatively analysed with the published chloroplast genomes of wild species of sweetpotato. High similarity and certain specificity were found among the chloroplast genomes of Ipomoea spp. Phylogenetic analysis could clearly distinguish wild species from cultivars. Ipomoea trifida and Ipomoea tabascana showed the closest relationship with the cultivars, and different haplotypes of ycf1 could be used to distinguish the cultivars from their wild relatives. The genetic structure was analyzed using variations in the chloroplast genome. Compared with traditional nuclear markers, the chloroplast markers designed based on the InDels on the chloroplast genome showed significant advantages.

Conclusions

Comparative analysis of chloroplast genomes of 107 cultivars and several wild species of sweetpotato was performed to help analyze the evolution, genetic structure and the development of chloroplast DNA markers of sweetpotato.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07544-y.

Comparative chloroplast genomes: insights into the evolution of the chloroplast genome of Camellia sinensis and the phylogeny of Camellia

BACKGROUND

Chloroplast genome resources can provide useful information for the evolution of plant species. Tea plant (Camellia sinensis) is among the most economically valuable member of Camellia. Here, we determined the chloroplast genome of the first natural triploid Chinary type tea ('Wuyi narcissus' cultivar of Camellia sinensis var. sinensis, CWN) and conducted the genome comparison with the diploid Chinary type tea (Camellia sinensis var. sinensis, CSS) and two types of diploid Assamica type teas (Camellia sinensis var. assamica: Chinese Assamica type tea, CSA and Indian Assamica type tea, CIA). Further, the evolutionary mechanism of the chloroplast genome of Camellia sinensis and the relationships of Camellia species based on chloroplast genome were discussed.

RESULTS

Comparative analysis showed the evolutionary dynamics of chloroplast genome of Camellia sinensis were the repeats and insertion-deletions (indels), and distribution of the repeats, indels and substitutions were significantly correlated. Chinese tea and Indian tea had significant differences in the structural characteristic and the codon usage of the chloroplast genome. Analysis of sequence characterized amplified region (SCAR) using sequences of the intergenic spacers (trnE/trnT) showed none of 292 different Camellia sinensis cultivars had similar sequence characteristic to triploid CWN, but the other four Camellia species did. Estimations of the divergence time showed that CIA diverged from the common ancestor of two Assamica type teas about 6.2 Mya (CI: 4.4-8.1 Mya). CSS and CSA diverged to each other about 0.8 Mya (CI: 0.4-1.5 Mya). Moreover, phylogenetic clustering was not exactly consistent with the current taxonomy of Camellia.

CONCLUSIONS

The repeat-induced and indel-induced mutations were two important dynamics contributed to the diversification of the chloroplast genome in Camellia sinensis, which were not mutually exclusive. Chinese tea and Indian tea might have undergone different selection pressures. Chloroplast transfer occurred during the polyploid evolution in Camellia sinensis. In addition, our results supported the three different domestication origins of Chinary type tea, Chinese Assamica type tea and Indian Assamica type tea. And, the current classification of some Camellia species might need to be further discussed.