Plastid phylogenomic insights into the evolution of the Caprifoliaceae s.l. (Dipsacales)

Bioactive Compounds from Medicinal Plants as Potential Adjuvants in the Treatment of Mild Acne Vulgaris

In recent years, there has been a growing interest in the use of medicinal plants and phytochemicals as potential treatments for acne vulgaris. This condition, characterized by chronic inflammation, predominantly affects adolescents and young adults. Conventional treatment typically targets the key factors contributing to its development: the proliferation of Cutibacterium acnes and the associated inflammation. However, these treatments often involve the use of potent drugs. As a result, the exploration of herbal medicine as a complementary approach has emerged as a promising strategy. By harnessing the therapeutic properties of medicinal plants and phytochemicals, it may be possible to address acne vulgaris while minimizing the reliance on strong drugs. This approach not only offers potential benefits for individuals seeking alternative treatments but also underscores the importance of natural remedies of plant origin in dermatological care. The primary aim of this study was to assess the antimicrobial, antioxidant, and anti-inflammatory properties of plants and their phytochemical constituents in the management of mild acne vulgaris. A comprehensive search of scientific databases was conducted from 2018 to September 2023. The findings of this review suggest that medicinal plants and their phytochemical components hold promise as treatments for mild acne vulgaris. However, it is crucial to note that further research employing high-quality evidence and standardized methodologies is essential to substantiate their efficacy and safety profiles.

Phylogenomics and plastomics offer new evolutionary perspectives on Kalanchoideae (Crassulaceae)

BACKGROUND AND AIMS

Kalanchoideae is one of three subfamilies within Crassulaceae and contains four genera. Despite previous efforts, the phylogeny of Kalanchoideae remains inadequately resolved with persistent issues including low support, unstructured topologies and polytomies. This study aimed to address two central objectives: (1) resolving the pending phylogenetic questions within Kalanchoideae by using organelle-scale 'barcodes' (plastomes) and nuclear data; and (2) investigating interspecific diversity patterns among Kalanchoideae plastomes.

METHODS

To explore the plastome evolution in Kalanchoideae, we newly sequenced 38 plastomes representing all four constituent genera (Adromischus, Cotyledon, Kalanchoe and Tylecodon). We performed comparative analyses of plastomic features, including GC and gene contents, gene distributions at the IR (inverted repeat) boundaries, nucleotide divergence, plastomic tRNA (pttRNA) structures and codon aversions. Additionally, phylogenetic inferences were inferred using both the plastomic dataset (79 genes) and nuclear dataset (1054 genes).

KEY RESULTS

Significant heterogeneities were observed in plastome lengths among Kalanchoideae, strongly correlated with LSC (large single copy) lengths. Informative diversities existed in the gene content at SSC/IRa (small single copy/inverted repeat a), with unique patterns individually identified in Adromischus leucophyllus and one major Kalanchoe clade. The ycf1 gene was assessed as a shared hypervariable region among all four genera, containing nine lineage-specific indels. Three pttRNAs exhibited unique structures specific to Kalanchoideae and the genera Adromischus and Kalanchoe. Moreover, 24 coding sequences revealed a total of 41 lineage-specific unused codons across all four constituent genera. The phyloplastomic inferences clearly depicted internal branching patterns in Kalanchoideae. Most notably, by both plastid- and nuclear-based phylogenies, our research offers the first evidence that Kalanchoe section Eukalanchoe is not monophyletic.

CONCLUSIONS

This study conducted comprehensive analyses on 38 newly reported Kalanchoideae plastomes. Importantly, our results not only reconstructed well-resolved phylogenies within Kalanchoideae, but also identified highly informative unique markers at the subfamily, genus and species levels. These findings significantly enhance our understanding of the evolutionary history of Kalanchoideae.

A systematic approach for authentication of medicinal Patrinia species using an integration of morphological, chemical and molecular methods

Four common Patrinia species, including P. heterophylla, P. monandra, P. scabiosifolia and P. villosa, have been documented as herbal medicines with various clinical applications, such as anti-cancer, anti-diarrhea and sedative. However, the authentication of medicinal Patrinia species poses a problem, particularly with the processed herbal materials. This study aimed to systematically authenticate the four medicinal Patrinia species in the market using morphological and chemical characterization, as well as DNA markers. We found the species identity authenticated by traditional morphologies were in good agreement with both chemical and molecular results. The four species showed species-specific patterns in chromatographic profiles with distinct chemical markers. We also revealed the power of complete chloroplast genomes in species authentication. The sequences of targeted loci, namely atpB, petA, rpl2-rpl23 and psaI-ycf4, contained informative nucleotides for the species differentiation. Our results also facilitate authentication of medicinal Patrinia species using new DNA barcoding markers. To the best of our knowledge, this is the first report on the application of morphology, chemical fingerprinting, complete chloroplast genomes and species-specific Insertion-Deletions (InDels) in differentiating Patrinia species. This study reported on the power of a systematic, multidisciplinary approach in authenticating medicinal Patrinia species.

Plastid genome data provide new insights into the dynamic evolution of the tribe Ampelopsideae (Vitaceae)

BACKGROUND

Ampelopsideae J. Wen & Z.L. Nie is a small-sized tribe of Vitaceae Juss., including ca. 47 species from four genera showing a disjunct distribution worldwide across all the continents except Antarctica. There are numerous species from the tribe that are commonly used as medicinal plants with immune-modulating, antimicrobial, and anti-hypertensive properties. The tribe is usually recognized into three clades, i.e., Ampelopsis Michx., Nekemias Raf., and the Southern Hemisphere clade. However, the relationships of the three clades differ greatly between the nuclear and the plastid topologies. There has been limited exploration of the chloroplast phylogenetic relationships within Ampelopsideae, and studies on the chloroplast genome structure of this tribe are only available for a few individuals. In this study, we aimed to investigate the evolutionary characteristics of plastid genomes of the tribe, including their genome structure and evolutionary insights.

RESULTS

We sequenced, assembled, and annotated plastid genomes of 36 species from the tribe and related taxa in the family. Three main clades were recognized within Ampelopsideae, corresponding to Ampelopsis, Nekemias, and the Southern Hemisphere lineage, respectively, and all with 100% bootstrap supports. The genome sequences and content of the tribe are highly conserved. However, comparative analyses suggested that the plastomes of Nekemias demonstrate a contraction in the large single copy region and an expansion in the inverted repeat region, and possess a high number of forward and palindromic repeat sequences distinct from both Ampelopsis and the Southern Hemisphere taxa.

CONCLUSIONS

Our results highlighted plastome variations in genome length, expansion or contraction of the inverted repeat region, codon usage bias, and repeat sequences, are corresponding to the three lineages of the tribe, which probably faced with different environmental selection pressures and evolutionary history. This study provides valuable insights into understanding the evolutionary patterns of plastid genomes within the Ampelopsideae of Vitaceae.

Plastome phylogenomics and morphological traits analyses provide new insights into the phylogenetic position, species delimitation and speciation of Triplostegia (Caprifoliaceae)

BACKGROUND

The genus Triplostegia contains two recognized species, T. glandulifera and T. grandiflora, but its phylogenetic position and species delimitation remain controversial. In this study, we assembled plastid genomes and nuclear ribosomal DNA (nrDNA) cistrons sampled from 22 wild Triplostegia individuals, each from a separate population, and examined these with 11 recently published Triplostegia plastomes. Morphological traits were measured from herbarium specimens and wild material, and ecological niche models were constructed.

RESULTS

Triplostegia is a monophyletic genus within the subfamily Dipsacoideae comprising three monophyletic species, T. glandulifera, T. grandiflora, and an unrecognized species Triplostegia sp. A, which occupies much higher altitude than the other two. The new species had previously been misidentified as T. glandulifera, but differs in taproot, leaf, and other characters. Triplotegia is an old genus, with stem age 39.96 Ma, and within it T. glandulifera diverged 7.94 Ma. Triplostegia grandiflora and sp. A diverged 1.05 Ma, perhaps in response to Quaternary climate fluctuations. Niche overlap between Triplostegia species was positively correlated with their phylogenetic relatedness.

CONCLUSIONS

Our results provide new insights into the species delimitation of Triplostegia, and indicate that a taxonomic revision of Triplostegia is needed. We also identified that either rpoB-trnC or ycf1 could serve as a DNA barcode for Triplostegia.

Complete chloroplast genome assembly and phylogenetic analysis of blackcurrant (Ribes nigrum), red and white currant (Ribes rubrum), and gooseberry (Ribes uva-crispa) provide new insights into the phylogeny of Grossulariaceae

Background

Blackcurrant (Ribes nigrum), red currant (R. rubrum), white currant (R. rubrum), and gooseberry (R. uva-crispa) belong to Grossulariaceae and are popular small-berry crops worldwide. The lack of genomic data has severely limited their systematic classification and molecular breeding.

Methods

The complete chloroplast (cp) genomes of these four taxa were assembled for the first time using MGI-DNBSEQ reads, and their genome structures, repeat elements and protein-coding genes were annotated. By genomic comparison of the present four and previous released five Ribes cp genomes, the genomic variations were identified. By phylogenetic analysis based on maximum-likelihood and Bayesian methods, the phylogeny of Grossulariaceae and the infrageneric relationships of the Ribes were revealed.

Results

The four cp genomes have lengths ranging from 157,450 to 157,802 bp and 131 shared genes. A total of 3,322 SNPs and 485 Indels were identified from the nine released Ribes cp genomes. Red currant and white currant have 100% identical cp genomes partially supporting the hypothesis that white currant (R. rubrum) is a fruit color variant of red currant (R. rubrum). The most polymorphic genic and intergenic region is ycf1 and trnT-psbD, respectively. The phylogenetic analysis demonstrated the monophyly of Grossulariaceae in Saxifragales and the paraphyletic relationship between Saxifragaceae and Grossulariaceae. Notably, the Grossularia subgenus is well nested within the Ribes subgenus and shows a paraphyletic relationship with the co-ancestor of Calobotrya and Coreosma sections, which challenges the dichotomous subclassification of the Ribes genus based on morphology (subgenus Ribes and subgenus Grossularia). These data, results, and insights lay a foundation for the phylogenetic research and breeding of Ribes species.

Complete chloroplast genome of Lonicera ligustrina Wall. (Caprifoliaceae) and its phylogenetic implications

Lonicera ligustrina is a folk medicinal herb in China and India with highly potential medicinal value. Here, we reported the complete chloroplast (cp) genome of L. ligustrina (GenBank accession number: ON968694). The cp genome was 155,330 bp long, with a large single-copy region (LSC) of 88,855 bp and a small single-copy region (SSC) of 18,647 bp separated by a pair of inverted repeats (IRs) of 23,914 bp. We also reconstructed the phylogeny of Lonicera using maximum likelihood (ML) method, including our data and previously reported cp genomes of related taxa. The current study indicated that L. ligustrina is sister with the Nintooa clade of subgen. Lonicera.

Target sequence capture data shed light on the deeper evolutionary relationships of subgenus Chamaecerasus in Lonicera (Caprifoliaceae)

The genus Lonicera L. is widely distributed in the north temperate zone and is well-known for its high species richness and morphological diversity. Previous studies have suggested that many sections of Lonicera are not monophyletic and phylogenetic relationships within the genus are still poorly resolved. In this study, we sampled 37 accessions of Lonicera, covering four sections of subgenus Chamaecerasus plus six outgroup taxa, to recover the main clades of Lonicera based on sequences of nuclear loci generated by target enrichment and cpDNA from genome skimming. We found extensive cytonuclear discordance across the subgenus. Both nuclear and plastid phylogenetic analyses supported subgenus Chamaecerasus sister to subgenus Lonicera. Within subgenus Chamaecerasus, sections Isika and Nintooa were each polyphyletic. Based on the nuclear and chloroplast phylogenies, we propose to merge Lonicera korolkowii into section Coeloxylosteum and Lonicera caerulea into section Nintooa. In addition, Lonicera is estimated to have originated in the mid Oligocene (26.45 Ma). The stem age of section Nintooa was estimated to be 17.09 Ma (95% HPD: 13.30-24.45). The stem age of subgenus Lonicera was estimated to be 16.35 Ma (95% HPD: 14.12-23.66). Ancestral area reconstruction analyses indicate that subgenus Chamaecerasus originated in East Asia and Central Asia. In addition, sections Coeloxylosteum and Nintooa originated in East Asia, with subsequent dispersals into other areas. The aridification of the Asian interior likely promoted the rapid radiation of sections Coeloxylosteum and Nintooa within this region. Moreover, our biogeographic analysis fully supports the Bering and the North Atlantic Land Bridge hypotheses for the intercontinental migrations in the Northern Hemisphere. Overall, this study provides new insights into the taxonomically complex lineages of subgenus Chamaecerasus and the process of speciation.

A phylogenomic analysis of Lonicera and its bearing on the evolution of organ fusion

PREMISE

The ~140 species of Lonicera are characterized by variously fused leaves, bracteoles, and ovaries, making it a model system for studying the evolution and development of organ fusion. However, previous phylogenetic analyses, based mainly on chloroplast DNA markers, have yielded uncertain and conflicting results. A well-supported phylogeny of Lonicera will allow us to trace the evolutionary history of organ fusion.

METHODS

We inferred the phylogeny of Lonicera using restriction site-associated DNA sequencing (RADSeq), sampling all major clades and 18 of the 23 subsections. This provided the basis for inferring the evolution of five fusion-related traits.

RESULTS

RADSeq data yielded a well-resolved and well-supported phylogeny. The two traditionally recognized subgenera (Periclymenum and Chamaecerasus), three of the four sections (Isoxylosteum, Coeloxylosteum, and Nintooa), and half of the subsections sampled were recovered as monophyletic. However, the large and heterogeneous section Isika was strongly supported as paraphyletic. Nintooa, a clade of ~22 mostly vine-forming species, including L. japonica, was recovered in a novel position, raising the possibility of cytonuclear discordance. We document the parallel evolution of fused leaves, bracteoles, and ovaries, with rare reversals. Most strikingly, complete cupules, in which four fused bracteoles completely enclose two unfused ovaries, arose at least three times. Surprisingly, these appear to have evolved directly from ancestors with free bracteoles instead of partial cupules.

CONCLUSIONS

We provide the most comprehensive and well-supported phylogeny of Lonicera to date. Our inference of multiple evolutionary shifts in organ fusion provides a solid foundation for in depth developmental and functional analyses.

Complete chloroplast genome and phylogenetic analysis of Lonicera caerulea var. edulis (Caprifoliaceae)

Lonicera caerulea var. edulis, known as "blue honeysuckle" or "Haskap," is a deciduous shrub that belongs to the Caprifoliaceae family. Characterized by the high cold hardiness and high quality of fruit, it has become a novel cash crop in cold regions worldwide. The lack of available chloroplast (cp) genome information limits studies of its molecular breeding and phylogeny. Here, the complete cp genome of Lonicera caerulea var. edulis was assembled and characterized for the first time. It was 155,142 bp in length with 38.43% of GC content, including 23,841 bp inverted repeat regions (IRs), an 88,737 bp large single-copy region (LSC), and an 18,723 bp small single-copy region (SSC). A total of 132 genes, including 85 protein-coding genes, 8 rRNA genes, and 39 tRNA genes were annotated. Phylogenetic analysis indicated that L. caerulea var. edulis was closely related to L. tangutica. These data and results provide a valuable resource for the development of breeding tools and genetic diversity studies for L. caerulea.

Phylogenomics and phylogeography of Menispermum (Menispermaceae)

INTRODUCTION

Phylogenomics have been widely used to resolve ambiguous and controversial evolutionary relationships among plant species and genera, and the identification of unique indels in plastomes may even help to understand the evolution of some plant families. Menispermum L. (Menispermaceae) consists of three species, M. dauricum DC., M. canadense L., and M. mexicanum Rose, which are disjuncly distributed among East Asia, Eastern North America and Mexico. Taxonomists continue to debate whether M. mexicanum is a distinct species, a variety of M. dauricum, or simply a synonym of M. canadense. To date, no molecular systematics studies have included this doubtful species in phylogenetic analyses.

METHODS

In this study, we examined phylogenomics and phylogeography of Menispermum across its entire range using 29 whole plastomes of Menispermaceae and 18 ITS1&ITS2 sequences of Menispermeae. We reconstructed interspecific relationships of Menispermum and explored plastome evolution in Menispermaceae, revealing several genomic hotspot regions for the family.

RESULTS AND DISCUSSION

Phylogenetic and network analyses based on whole plastome and ITS1&ITS2 sequences show that Menispermum clusters into two clades with high support values, Clade A (M. dauricum) and Clade B (M. canadense + M. mexicanum). However, M. mexicanum is nested within M. canadense and, as a result, we support that M. mexicanum is a synonym of M. canadense. We also identified important molecular variations in the plastomes of Menispermaceae. Several indels and consequently premature terminations of genes occur in Menispermaceae. A total of 54 regions were identified as the most highly variable plastome regions, with nucleotide diversity (Pi) values > 0.05, including two coding genes (matK, ycf1), four introns (trnK intron, rpl16 intron, rps16 intron, ndhA intron), and 48 intergenic spacer (IGS) regions. Of these, four informative hotspot regions (trnH-psbA, ndhF-rpl32, trnK-rps16, and trnP-psaJ) should be especially useful for future studies of phylogeny, phylogeography and conservation genetics of Menispermaceae.

Population Genomics Study and Implications for the Conservation of Zabelia tyaihyonii Based on Genotyping-By-Sequencing

Zabelia tyaihyonii (Nakai) Hisauti and H. Hara is a perennial shrub endemic to Republic of Korea that grows naturally in only a very limited region of the dolomite areas of Gangwon-do and Chungcheongbuk-do Provinces in the Republic of Korea. Given its geographical characteristics, it is more vulnerable than more widely distributed species. Despite the need for comprehensive information to support conservation, population genetic information for this species is very scarce. In this study, we analyzed the genetic diversity and population structure of 94 individuals from six populations of Z. tyaihyonii using a genotyping-by-sequencing (GBS) approach to provide important information for proper conservation and management. Our results, based on 3088 single nucleotide polymorphisms (SNPs), showed a mean expected heterozygosity (He) of 0.233, no sign of within-population inbreeding (G_IS that was close to or even below zero in all populations), and a high level of genetic differentiation (F_ST = 0.170). Analysis of molecular variance (AMOVA) indicated that the principal molecular variance existed within populations (84.5%) rather than among populations (17.0%). We suggested that six management units were proposed for conservation considering Bayesian structure analysis and phylogenetic analysis, and given the various current situations faced by Z. tyaihyonii, it is believed that not only the in situ conservation but also the ex situ conservation should be considered.

Plastome phylogenomics provide new perspective into the phylogeny and evolution of Betulaceae (Fagales)

BACKGROUND

Betulaceae is a relatively small but morphologically diverse family, with many species having important economic and ecological values. Although plastome structure of Betulaceae has been reported sporadically, a comprehensive exploration for plastome evolution is still lacking. Besides, previous phylogenies had been constructed based on limited gene fragments, generating unrobust phylogenetic framework and hindering further studies on divergence ages, biogeography and character evolution. Here, 109 plastomes (sixteen newly assembled and 93 previously published) were subject to comparative genomic and phylogenomic analyses to reconstruct a robust phylogeny and trace the diversification history of Betulaceae.

RESULTS

All Betulaceae plastomes were highly conserved in genome size, gene order, and structure, although specific variations such as gene loss and IR boundary shifts were revealed. Ten divergent hotspots, including five coding regions (P_i > 0.02) and five noncoding regions (P_i > 0.035), were identified as candidate DNA barcodes for phylogenetic analysis and species delimitation. Phylogenomic analyses yielded high-resolution topology that supported reciprocal monophyly between Betula and Alnus within Betuloideae, and successive divergence of Corylus, Ostryopsis, and Carpinus-Ostrya within Coryloideae. Incomplete lineage sorting and hybridization may be responsible for the mutual paraphyly between Ostrya and Carpinus. Betulaceae ancestors originated from East Asia during the upper Cretaceous; dispersals and subsequent vicariance accompanied by historical environment changes contributed to its diversification and intercontinental disjunction. Ancestral state reconstruction indicated the acquisition of many taxonomic characters was actually the results of parallel or reversal evolution.

CONCLUSIONS

Our research represents the most comprehensive taxon-sampled and plastome-level phylogenetic inference for Betulaceae to date. The results clearly document global patterns of plastome structural evolution, and established a well-supported phylogeny of Betulaceae. The robust phylogenetic framework not only provides new insights into the intergeneric relationships, but also contributes to a perspective on the diversification history and evolution of the family.

Complete chloroplast genome sequences of Acanthocalyx alba and Acanthocalyx nepalensis subsp. delavayi (Caprifoliaceae)

In this study, the chloroplast genomes of Acanthocalyx alba (Hand.-Mazz., 1925) and Acanthocalyx nepalensis subsp. delavayi (Franchet, 1885) were sequenced, and their total lengths were 148,720 bp and 149,253 bp, respectively. The A. alba genome contained two inverted repeat regions (IRs) of 21,849 bp, a large single-copy region (LSC) of 89,084 bp, and a small single-copy region (SSC) of 15,938 bp, whereas A. nepalensis subsp. delavayi contained two IRs of 21,736 bp, one LSC of 89,034 bp, and one SSC of 16,747 bp. The chloroplast genomes of both A. alba and A. nepalensis subsp. delavayi contained 109 genes, including 72 mRNA, 33 tRNA, and four rRNA genes. Phylogenetic analysis suggested that A. alba is in a clade with A. nepalensis subsp. delavayi. This study provides useful data for further phylogenetic studies of A. alba and A. nepalensis subsp. delavayi.

Complete plastome sequence of Lonicera gynochlamydea Hemsl. (Caprifoliaceae)

We report and characterize the complete plastome of Lonicera gynochlamydea Hemsl. L. gynochlamydea is a shrub, belonging to the family Caprifoliaceae. Our results show that the length of the complete plastome is 154,643 bp, including 131 genes consisting of 84 protein-coding genes, 39 tRNA genes, and eight rRNA genes. The plastome exhibits the typical quadripartite structure and gene content of angiosperms, composed of two inverted repeats (IRs) regions of 23,846 bp, a large single-copy (LSC) region of 88,298 bp, and a small single-copy (SSC) region of 18,653 bp. The total G/C content in L. gynochlamydea plastome is 38.4%. The complete plastome sequence of L. gynochlamydea will make contributions to the conservation genetics of this species as well as to phylogenetic studies in Caprifoliaceae.

Plastid Phylogenomics and Plastome Evolution of Nandinoideae (Berberidaceae)

Subfamily Nandinoideae Heintze (Berberidaceae), comprising four genera and ca. 19 species, is disjunctively distributed in eastern North America vs. Eurasia (eastern Asia, Central Asia, Middle East, and southeastern Europe), and represents an ideal taxon to explore plastid phylogenomics and plastome evolution in Berberidaceae. Many species of this subfamily have been listed as national or international rare and endangered plants. In this study, we sequenced and assembled 20 complete plastomes, representing three genera and 13 species of Nandinoideae. Together with six plastomes from GenBank, a total of 26 plastomes, representing all four genera and 16 species of Nandinoideae, were used for comparative genomic and phylogenomic analyses. These plastomes showed significant differences in overall size (156,626-161,406 bp), which is mainly due to the expansion in inverted repeat (IR) regions and/or insertion/deletion (indel) events in intergenic spacer (IGS) regions. A 75-bp deletion in the ndhF gene occurred in Leontice and Gymnospermium when compared with Nandina and Caulophyllum. We found a severe truncation at the 5' end of ycf1 in three G. altaicum plastomes, and a premature termination of ropC1 in G. microrrhynchum. Our phylogenomic results support the topology of {Nandina, [Caulophyllum, (Leontice, Gymnospermium)]}. Within the core genus Gymnospermium, we identified G. microrrhynchum from northeastern Asia (Clade A) as the earliest diverging species, followed by G. kiangnanense from eastern China (Clade B), while the rest species clustered into the two sister clades (C and D). Clade C included three species from West Tianshan (G. albertii, G. darwasicum, G. vitellinum). Clade D consisted of G. altaicum from northern Central Asia, plus one species from the Caucasus Mountains (G. smirnovii) and three from southeastern Europe (G. odessanum, G. peloponnesiacum, G. scipetarum). Overall, we identified 21 highly variable plastome regions, including two coding genes (rpl22, ycf1) and 19 intergenic spacer (IGS) regions, all with nucleotide diversity (Pi) values > 0.02. These molecular markers should serve as powerful tools (including DNA barcodes) for future phylogenetic, phylogeographic and conservation genetic studies.

Complete Chloroplast Genome Sequence of Triosteum sinuatum, Insights into Comparative Chloroplast Genomics, Divergence Time Estimation and Phylogenetic Relationships among Dipsacales

Triosteum himalayanum, Triosteum pinnatifidum (Triosteum L., Caprifoliaceae, Dipsacales) are widely distributed in China while Triosteum sinuatum mainly occurrs in northeast China. Few reports have been determined on the genus Triosteum. In the present research, we sequenced 2 chloroplast genomes of Triosteum and analyzed 18 chloroplast genomes, trying to explore the sequence variations and phylogeny of genus Triosteum in the order Dipsacales. The chloroplast genomes of the genus Triosteum ranged from 154,579 bp to 157,178 bp, consisting of 132 genes (86 protein-coding genes, 38 transfer RNA genes, and 8 ribosomal RNA genes). Comparative analyses and phylogenetic analysis supported the division of Dipsacales into two clades, Adoxaceae and six other families. Among the six families, a clade of Valerianaceae+Dipsacaceae was recovered as a sister to a clade of Morinaceae+Linnaeaceae. A closer relationship of T. himalayanum and T. pinnatifidum among three species was revealed. Our research supported that Loniceraferdinandi and Triosteum was closely related. Zabelia had a closer relationship with Linnaea borealis and Dipelta than Morinaceae. The divergence between T. sinuatum and two other species in Triosteum was dated to 13.4 mya.

The complete chloroplast genome of Onobrychis gaubae (Fabaceae-Papilionoideae): comparative analysis with related IR-lacking clade species

BACKGROUND

Plastome (Plastid genome) sequences provide valuable markers for surveying evolutionary relationships and population genetics of plant species. Papilionoideae (papilionoids) has different nucleotide and structural variations in plastomes, which makes it an ideal model for genome evolution studies. Therefore, by sequencing the complete chloroplast genome of Onobrychis gaubae in this study, the characteristics and evolutionary patterns of plastome variations in IR-loss clade were compared.

RESULTS

In the present study, the complete plastid genome of O. gaubae, endemic to Iran, was sequenced using Illumina paired-end sequencing and was compared with previously known genomes of the IRLC species of legumes. The O. gaubae plastid genome was 122,688 bp in length and included a large single-copy (LSC) region of 81,486 bp, a small single-copy (SSC) region of 13,805 bp and one copy of the inverted repeat (IR_b) of 29,100 bp. The genome encoded 110 genes, including 76 protein-coding genes, 30 transfer RNA (tRNA) genes and four ribosome RNA (rRNA) genes and possessed 83 simple sequence repeats (SSRs) and 50 repeated structures with the highest proportion in the LSC. Comparative analysis of the chloroplast genomes across IRLC revealed three hotspot genes (ycf1, ycf2, clpP) which could be used as DNA barcode regions. Moreover, seven hypervariable regions [trnL(UAA)-trnT(UGU), trnT(GGU)-trnE(UUC), ycf1, ycf2, ycf4, accD and clpP] were identified within Onobrychis, which could be used to distinguish the Onobrychis species. Phylogenetic analyses revealed that O. gaubae is closely related to Hedysarum. The complete O. gaubae genome is a valuable resource for investigating evolution of Onobrychis species and can be used to identify related species.

CONCLUSIONS

Our results reveal that the plastomes of the IRLC are dynamic molecules and show multiple gene losses and inversions. The identified hypervariable regions could be used as molecular markers for resolving phylogenetic relationships and species identification and also provide new insights into plastome evolution across IRLC.

Valerianapraecipitis (Caprifoliaceae), a species new to science and endemic to Central Chile

The species Valerianapraecipitis (Caprifoliaceae), new to science and endemic to the Ñuble Region, Central Chile, is formally described. Morphological data support its placement in a new species, clearly different from V.philippiana. A detailed description, insights about its habitat and ecology, distribution map and illustration are provided. A table of comparison is also given with the morphological characters discriminating V.praecipitis from V.philippiana. The species is assessed as Endangered (EN) under the IUCN categories.

Phylogenomic analyses of the East Asian endemic Abelia (Caprifoliaceae) shed insights into the temporal and spatial diversification history with widespread hybridization

BACKGROUND AND AIMS

Abelia (Caprifoliaceae) is a small genus with five species, including one artificial hybrid and several natural hybrids. The genus has a discontinuous distribution in Mainland China, Taiwan Island and the Ryukyu Islands, providing a model system to explore the mechanisms of species dispersal in the East Asian flora. However, the current phylogenetic relationships within Abelia remain uncertain.

METHODS

We reconstructed the phylogenetic relationships within Abelia using nuclear loci generated by target enrichment and plastomes from genome skimming. Divergence time estimation, ancestral area reconstruction and ecological niche modelling (ENM) were used to examine the diversification history of Abelia.

KEY RESULTS

We found extensive cytonuclear discordance across the genus. By integrating lines of evidence from molecular phylogenies, divergence times and morphology, we propose to merge Abelia macrotera var. zabelioides into A. uniflora. Network analyses suggested that there have been multiple widespread hybridization events among Abelia species. These hybridization events may have contributed to the speciation mechanism and resulted in the high observed morphological diversity. The diversification of Abelia began in the early Eocene, followed by A. chinensis var. ionandra colonizing Taiwan Island during the Middle Miocene. The ENM results suggested an expansion of climatically suitable areas during the Last Glacial Maximum and range contraction during the Last Interglacial. Disjunction between the Himalayan-Hengduan Mountain region and Taiwan Island is probably the consequence of topographical isolation and postglacial contraction.

CONCLUSIONS

We used genomic data to reconstruct the phylogeny of Abelia and found a clear pattern of reticulate evolution in the group. In addition, our results suggest that shrinkage of postglacial range and the heterogeneity of the terrain have led to the disjunction between Mainland China and Taiwan Island. This study provides important new insights into the speciation process and taxonomy of Abelia.

Characterization of the chloroplast genome of Lonicera ruprechtiana Regel and comparison with other selected species of Caprifoliaceae

Lonicera ruprechtiana Regel is widely used as a greening tree in China and also displays excellent pharmacological activities. The phylogenetic relationship between L. ruprechtiana and other members of Caprifoliaceae remains unclear. In this study, the complete cp genome of L. ruprechtiana was identified using high-throughput Illumina pair-end sequencing data. The circular cp genome was 154,611 bp long and has a large single-copy region of 88,182 bp and a small single-copy region of 18,713 bp, with the two parts separated by two inverted repeat (IR) regions (23,858 bp each). A total of 131 genes were annotated, including 8 ribosomal RNAs, 39 transfer RNAs, and 84 protein-coding genes (PCGs). In addition, 49 repeat sequences and 55 simple sequence repeat loci of 18 types were also detected. Codon usage analysis demonstrated that the Leu codon is preferential for the A/U ending. Maximum-likelihood phylogenetic analysis using 22 Caprifoliaceae species revealed that L. ruprechtiana was closely related to Lonicera insularis. Comparison of IR regions revealed that the cp genome of L. ruprechtiana was largely conserved with that of congeneric species. Moreover, synonymous (Ks) and non-synonymous (Ka) substitution rate analysis showed that most genes were under purifying selection pressure; ycf3, and some genes associated with subunits of NADH dehydrogenase, subunits of the cytochrome b/f complex, and subunits of the photosystem had been subjected to strong purifying selection pressure (Ka/Ks < 0.1). This study provides useful genetic information for future study of L. ruprechtiana evolution.

Deep Insights Into the Plastome Evolution and Phylogenetic Relationships of the Tribe Urticeae (Family Urticaceae)

Urticeae s.l., a tribe of Urticaceae well-known for their stinging trichomes, consists of more than 10 genera and approximately 220 species. Relationships within this tribe remain poorly known due to the limited molecular and taxonomic sampling in previous studies, and chloroplast genome (CP genome/plastome) evolution is still largely unaddressed. To address these concerns, we used genome skimming data-CP genome and nuclear ribosomal DNA (18S-ITS1-5.8S-ITS2-26S); 106 accessions-for the very first time to attempt resolving the recalcitrant relationships and to explore chloroplast structural evolution across the group. Furthermore, we assembled a taxon rich two-locus dataset of trnL-F spacer and ITS sequences across 291 accessions to complement our genome skimming dataset. We found that Urticeae plastomes exhibit the tetrad structure typical of angiosperms, with sizes ranging from 145 to 161 kb and encoding a set of 110-112 unique genes. The studied plastomes have also undergone several structural variations, including inverted repeat (IR) expansions and contractions, inversion of the trnN-GUU gene, losses of the rps19 gene, and the rpl2 intron, and the proliferation of multiple repeat types; 11 hypervariable regions were also identified. Our phylogenomic analyses largely resolved major relationships across tribe Urticeae, supporting the monophyly of the tribe and most of its genera except for Laportea, Urera, and Urtica, which were recovered as polyphyletic with strong support. Our analyses also resolved with strong support several previously contentious branches: (1) Girardinia as a sister to the Dendrocnide-Discocnide-Laportea-Nanocnide-Zhengyia-Urtica-Hesperocnide clade and (2) Poikilospermum as sister to the recently transcribed Urera sensu stricto. Analyses of the taxon-rich, two-locus dataset showed lower support but was largely congruent with results from the CP genome and nuclear ribosomal DNA dataset. Collectively, our study highlights the power of genome skimming data to ameliorate phylogenetic resolution and provides new insights into phylogenetic relationships and chloroplast structural evolution in Urticeae.

Plastome Diversity and Phylogenomic Relationships in Asteraceae

Plastid genomes are in general highly conserved given their slow evolutionary rate, and thus large changes in their structure are unusual. However, when specific rearrangements are present, they are often phylogenetically informative. Asteraceae is a highly diverse family whose evolution is long driven by polyploidy (up to 48x) and hybridization, both processes usually complicating systematic inferences. In this study, we generated one of the most comprehensive plastome-based phylogenies of family Asteraceae, providing information about the structure, genetic diversity and repeat composition of these sequences. By comparing the whole-plastome sequences obtained, we confirmed the double inversion located in the long single-copy region, for most of the species analyzed (with the exception of basal tribes), a well-known feature for Asteraceae plastomes. We also showed that genome size, gene order and gene content are highly conserved along the family. However, species representative of the basal subfamily Barnadesioideae-as well as in the sister family Calyceraceae-lack the pseudogene rps19 located in one inverted repeat. The phylogenomic analysis conducted here, based on 63 protein-coding genes, 30 transfer RNA genes and 21 ribosomal RNA genes from 36 species of Asteraceae, were overall consistent with the general consensus for the family's phylogeny while resolving the position of tribe Senecioneae and revealing some incongruences at tribe level between reconstructions based on nuclear and plastid DNA data.

Lineage-Specific Variation in IR Boundary Shift Events, Inversions, and Substitution Rates among Caprifoliaceae s.l. (Dipsacales) Plastomes

Caprifoliaceae s.l. plastid genomes (plastomes) show that one inversion and two inverted repeat boundary shifts occurred in the common ancestor of this family, after which the plastomes are generally conserved. This study reports plastome sequences of five additional species, Fedia cornucopiae, Valeriana fauriei, and Valerianella locusta from the subfamily Valerianoideae, as well as Dipsacus japonicus and Scabiosa comosa from the subfamily Dipsacoideae. Combined with the published plastomes, these plastomes provide new insights into the structural evolution of plastomes within the family. Moreover, the three plastomes from the subfamily Valerianoideae exhibited accelerated nucleotide substitution rates, particularly at synonymous sites, across the family. The patterns of accD sequence divergence in the family are dynamic with structural changes, including interruption of the conserved domain and increases in nonsynonymous substitution rates. In particular, the Valeriana accD gene harbors a large insertion of amino acid repeat (AAR) motifs, and intraspecific polymorphism with a variable number of AARs in the Valeriana accD gene was detected. We found a correlation between intron losses and increased ratios of nonsynonymous to synonymous substitution rates in the clpP gene with intensified positive selection. In addition, two Dipsacoideae plastomes revealed the loss of the plastid-encoded rps15, and a potential functional gene transfer to the nucleus was confirmed.

An updated phylogenetic and biogeographic analysis based on genome skimming data reveals convergent evolution of shrubby habit in Clematis in the Pliocene and Pleistocene

Convergent evolution, often viewed as the inevitable outcome of natural selection, has received special attention since the time of Darwin. Clematis is well known for its climbing habit, but it has some shrubby species, known as sect. Fruticella s.l. The shrubby Clematis species are distributed in the dry habitats of Central Asia and adjacent areas showing possible convergent evolution. In this study, we assembled the complete plastome and nuclear ribosomal DNA (nrDNA) sequences of 56 Clematis species, representing most sections and covering most of the shrubby species, to reconstruct their evolutionary histories. Using both maximum likelihood and Bayesian methods, the plastome and nrDNA datasets generated similar, but not identical, phylogenetic relationships, which are better resolved than in previous studies. Then, molecular dating, historical range reconstruction, and character optimization analyses were conducted based on this updated phylogenetic framework. All the morphological characters widely used for taxonomy were shown to have evolved multiple times. Molecular dating inferred that Clematis diverged from its sister in the mid Miocene, and all six major clades of Clematis originated during the late Miocene, with a species radiation during the Pliocene to Pleistocene. The results clearly showed that the shrubby habit evolved independently in four lineages of Clematis in Asia. We also revealed that the shrubby lineages have emerged since the very beginning of Pliocene. Asian monsoon variation in the Pliocene and glacial period fluctuation in the Pleistocene may be the driving forces for the origin and diversification of the shrubby Clematis in Central Asia and adjacent dry areas.

Reconstructing Dipsacales phylogeny using Angiosperms353: issues and insights

PREMISE

Phylogenetic relationships within major angiosperm clades are increasingly well resolved, but largely informed by plastid data. Areas of poor resolution persist within the Dipsacales, including placement of Heptacodium and Zabelia, and relationships within the Caprifolieae and Linnaeeae, hindering our interpretation of morphological evolution. Here, we sampled a significant number of nuclear loci using a Hyb-Seq approach and used these data to infer the Dipsacales phylogeny and estimate divergence times.

METHODS

Sampling all major clades within the Dipsacales, we applied the Angiosperms353 probe set to 96 species. Data were filtered based on locus completeness and taxon recovery per locus, and trees were inferred using RAxML and ASTRAL. Plastid loci were assembled from off-target reads, and 10 fossils were used to calibrate dated trees.

RESULTS

Varying numbers of targeted loci and off-target plastomes were recovered from most taxa. Nuclear and plastid data confidently place Heptacodium with Caprifolieae, implying homoplasy in calyx morphology, ovary development, and fruit type. Placement of Zabelia, and relationships within the Caprifolieae and Linnaeeae, remain uncertain. Dipsacales diversification began earlier than suggested by previous angiosperm-wide dating analyses, but many major splitting events date to the Eocene.

CONCLUSIONS

The Angiosperms353 probe set facilitated the assembly of a large, single-copy nuclear dataset for the Dipsacales. Nevertheless, many relationships remain unresolved, and resolution was poor for woody clades with low rates of molecular evolution. We favor expanding the Angiosperms353 probe set to include more variable loci and loci of special interest, such as developmental genes, within particular clades.

Structural Mutations in the Organellar Genomes of Valeriana sambucifolia f. dageletiana (Nakai. ex Maekawa) Hara Show Dynamic Gene Transfer

Valeriana sambucifolia f. dageletiana (Nakai. ex Maekawa) Hara is a broad-leaved valerian endemic to Ulleung Island, a noted hot spot of endemism in Korea. However, despite its widespread pharmacological use, this plant remains comparatively understudied. Plant cells generally contain two types of organellar genomes (the plastome and the mitogenome) that have undergone independent evolution, which accordingly can provide valuable information for elucidating the phylogenetic relationships and evolutionary histories of terrestrial plants. Moreover, the extensive mega-data available for plant genomes, particularly those of plastomes, can enable researchers to gain an in-depth understanding of the transfer of genes between different types of genomes. In this study, we analyzed two organellar genomes (the 155,179 bp plastome and the 1,187,459 bp mitogenome) of V. sambucifolia f. dageletiana and detected extensive changes throughout the plastome sequence, including rapid structural mutations associated with inverted repeat (IR) contraction and genetic variation. We also described features characterizing the first reported mitogenome sequence obtained for a plant in the order Dipsacales and confirmed frequent gene transfer in this mitogenome. We identified eight non-plastome-originated regions (NPRs) distributed within the plastome of this endemic plant, for six of which there were no corresponding sequences in the current nucleotide sequence databases. Indeed, one of these unidentified NPRs unexpectedly showed certain similarities to sequences from bony fish. Although this is ostensibly difficult to explain, we suggest that this surprising association may conceivably reflect the occurrence of gene transfer from a bony fish to the plastome of an ancestor of V. sambucifolia f. dageletiana mediated by either fungi or bacteria.

Complete plastome sequence of Xylosma longifolia Clos. (Salicaceae)

Xylosma longifolia is a tree species within Salicaceae and is distributed in Guizhou, Yunnan, Fujian, Guangxi, Guangdong, and Hainan provinces of China as well as in Vietnam, Laos, and India. There are no studies utilizing the complete plastome of Xylosma longifolia in the current literature. Therefore, this report provides a reference for the plastid gene sequence of Xylosma longifolia, and it contributes to the phylogenetic placement and species identification. In this report, we described the complete plastome sequence of Xylosma longifolia. The complete plastome length of Xylosma longifolia is 156,938 bp and has the typical quadripartite structure and gene content of angiosperms, including two inverted repeat (IR) regions of 27,514 bp, a large single-copy (LSC) region of 85,221 bp and a small single-copy (SSC) region of 16,689 bp. The plastome contains 130 genes, including 86 protein coding genes, 36 tRNA genes, eight rRNA genes (5S rRNA, 4.5S rRNA, 16S rRNA, and 23S rRNA). The GC content of the plastome is 36.8%. The complete plastome sequence will be a valuable resource for studies involving the phylogenetic inference of Salicaceae.

Phylogenomic and evolutionary dynamics of inverted repeats across Angelica plastomes

BACKGROUND

Angelica L. (family Apiaceae) is an economically important genus comprising ca. One hundred ten species. Angelica species are found on all continents of the Northern Hemisphere, and East Asia hosts the highest number of species. Morphological characters such as fruit anatomy, leaf morphology and subterranean structures of Angelica species show extreme diversity. Consequently, the taxonomic classification of Angelica species is complex and remains controversial, as the classifications proposed by previous studies based on morphological data and molecular data are highly discordant. In addition, the phylogenetic relationships of major clades in the Angelica group, particularly in the Angelica s. s. clade, remain unclear. Chloroplast (cp) genome sequences have been widely used in phylogenetic studies and for evaluating genetic diversity.

RESULTS

In this study, we sequenced and assembled 28 complete cp genomes from 22 species, two varieties and two cultivars of Angelica. Combined with 36 available cp genomes in GenBank from representative clades of the subfamily Apioideae, the characteristics and evolutionary patterns of Angelica cp genomes were studied, and the phylogenetic relationships of Angelica species were resolved. The Angelica cp genomes had the typical quadripartite structure including a pair of inverted repeats (IRs: 5836-34,706 bp) separated by a large single-copy region (LSC: 76,657-103,161 bp) and a small single-copy region (SSC: 17,433-21,794 bp). Extensive expansion and contraction of the IR region were observed among cp genomes of Angelica species, and the pattern of the diversification of cp genomes showed high consistency with the phylogenetic placement of Angelica species. Species of Angelica were grouped into two major clades, with most species grouped in the Angelica group and A. omeiensis and A. sinensis grouped in the Sinodielsia with Ligusticum tenuissimum.

CONCLUSIONS

Our results further demonstrate the power of plastid phylogenomics in enhancing the phylogenetic reconstructions of complex genera and provide new insights into plastome evolution across Angelica L.

Comparative and Evolutionary Analyses on the Complete Plastomes of Five Kalanchoe Horticultural Plants

Many species of the genus Kalanchoe are important horticultural plants. They have evolved the Crassulacean acid metabolism (CAM) photosynthetic pathway to allow them to be better adapted to dry environments. Despite their importance, it is still debating whether Kalanchoe is monophyletic, and understanding the past diversification of this genus requires a tremendous amount of effort and work being devoted to the studies of morphological and molecular characters of this genus. However, molecular information, plastic sequence data, in particular, reported on Kalanchoe species is scarce, and this has posed a great challenge in trying to interpret the evolutionary history of this genus. In this study, plastomes of the five Kalanchoe species, including Kalanchoe daigremontiana, Kalanchoe delagoensis, Kalanchoe fedtschenkoi, Kalanchoe longiflora, and Kalanchoe pinnata, were sequenced and analyzed. The results indicate that the five plastomes are comparable in size, guanine-cytosine (GC) contents and the number of genes, which also demonstrate an insignificant difference in comparison with other species from the family Crassulaceae. About 224 simple sequence repeats (SSRs) and 144 long repeats were identified in the five plastomes, and most of these are distributed in the inverted repeat regions. In addition, highly divergent regions containing either single nucleotide polymorphism (SNP) or insertion or deletion (InDel) mutations are discovered, which could be potentially used for establishing phylogenetic relationships among members of the Kalanchoe genus in future studies. Furthermore, phylogenetic analyses suggest that Bryophyllum should be placed into one single genus as Kalanchoe. Further genomic analyses also reveal that several genes are undergone positive selection. Among them, 11 genes are involved in important cellular processes, such as cell survival, electron transfer, and may have played indispensable roles in the adaptive evolution of Kalanchoe to dry environments.

Phylogenetic and Comparative Analyses of Complete Chloroplast Genomes of Chinese Viburnum and Sambucus (Adoxaceae)

Phylogenetic analyses of complete chloroplast genome sequences have yielded significant improvements in our understanding of relationships in the woody flowering genus Viburnum (Adoxaceae, Dipsacales); however, these relationships were evaluated focusing only on Viburnum species within Central and South America and Southeast Asia. By contrast, despite being a hotspot of Viburnum diversity, phylogenetic relationships of Viburnum species in China are less well known. Here, we characterized the complete chloroplast (cp) genomes of 21 Viburnum species endemic to China, as well as three Sambucus species. These 24 plastomes were highly conserved in genomic structure, gene order and content, also when compared with other Adoxaceae. The identified repeat sequences, simple sequence repeats (SSRs) and highly variable plastid regions will provide potentially valuable genetic resources for further population genetics and phylogeographic studies on Viburnum and Sambucus. Consistent with previous combined phylogenetic analyses of 113 Viburnum species, our phylogenomic analyses based on the complete cp genome sequence dataset confirmed the sister relationship between Viburnum and the Sambucus-Adoxa-Tetradoxa-Sinadoxa group, the monophyly of four recognized sections in Flora of China (i.e., Viburnum sect. Tinus, Viburnum sect. Solenotinus, Viburnum sect. Viburnum and Viburnum sect. Pseudotinus) and the nonmonophyly of Viburnum sect. Odontotinus and Viburnum sect. Megalotinus. Additionally, our study confirmed the sister relationships between the clade Valvatotinus and Viburnum sect. Pseudotinus, as well as between Viburnum sect. Opulus and the Odontotinus-Megalotinus group. Overall, our results clearly document the power of the complete cp genomes in improving phylogenetic resolution, and will contribute to a better understanding of plastome evolution in Chinese Adoxaceae.

Characteristic of the complete chloroplast genome of Symphoricarpos sinensis, an important ornamental plant endemic to China

Symphoricarpos sinensis is an important ornamental plant in China. Due to the lack of efficient molecular markers, although recent studies supported that Symphoricarpos was phylogenetically closely related to Lonicera, Triosteum, Heptacodium, and Leycesteria, relationships among these taxa remained uncertain. To examine the phylogenetic position of Symphoricarpos within Caprifoliaceae, we characterized the complete chloroplast (cp) genome of S. sinensis. The results showed that S. sinensis shared a typical quadripartite structure as with most angiosperms. It was 155,738 bp in length, including two inverted repeat (IR) regions (24,010 bp), a small single-copy (SSC) region (18,942 bp) and a large single-copy (LSC) region (88,776 bp). Phylogenetic analysis demonstrated that Caprifoliaceae and Adoxaceae formed two distinct monophyletic clades, and S. sinensis was closely related to lonicera and Triosteum, albeit with low support. The whole cp genome of S. sinensis will be useful resources for future studies on phylogeny and conservation in Symphoricarpos.

Long-reads reveal that Rhododendron delavayi plastid genome contains extensive repeat sequences, and recombination exists among plastid genomes of photosynthetic Ericaceae

BACKGROUND

Rhododendron delavayi Franch. var. delavayi is a wild ornamental plant species in Guizhou Province, China. The lack of its plastid genome information seriously hinders the further application and conservation of the valuable resource.

METHODS

The complete plastid genome of R. delavayi was assembled from long sequence reads. The genome was then characterized, and compared with those of other photosynthetic Ericaceae species.

RESULTS

The plastid genome of R. delavayi has a typical quadripartite structure, and a length of 202,169 bp. It contains a large number of repeat sequences and shows preference for codon usage. The comparative analysis revealed the irregular recombination of gene sets, including rearrangement and inversion, in the large single copy region. The extreme expansion of the inverted repeat region shortened the small single copy, and expanded the full length of the genome. In addition, consistent with traditional taxonomy, R. delavayi with nine other species of the same family were clustered into Ericaceae based on the homologous protein-coding sequences of the plastid genomes. Thus, the long-read assembly of the plastid genome of R. delavayi would provide basic information for the further study of the evolution, genetic diversity, and conservation of R. delavayi and its relatives.