DNA Barcoding and Phylogenomic Analysis of the Genus Fritillaria in China Based on Complete Chloroplast Genomes

Wasabi Gone Wild? Origin and Characterization of the Complete Plastomes of Ulleung Island Wasabi (Eutrema japonicum; Brassicaceae) and Other Cultivars in Korea

Korean wasabi occurs naturally on the young oceanic, volcanic Ulleung Island off the east coast of the Korean Peninsula. Although the Ulleung Island wasabi is reported as Eutrema japonicum and has been suggested to be morphologically identical to cultivars in Korea, very little is known about its taxonomic identity and relationship with other cultivars. In this study, we sequenced the complete chloroplast DNA sequences of three naturally occurring Ulleung Island wasabi plants and six cultivars ('Daewang', 'Daruma', 'Micado', 'Orochi', 'Green Thumb', and 'Shogun') from continental Korea and determined the taxonomic identity of Korean wasabi on Ulleung Island. The size and organization of the complete chloroplast genomes of the nine accessions were nearly identical to those of previously reported wasabi cultivars. In addition, phylogenetic analysis based on the complete plastomes suggested that Ulleung Island wasabi most likely comprises various wasabi cultivars with three chlorotypes ('Shogun', 'Green Thumb', and a unique Chusan type). Based on the complete plastomes, we identified eight chlorotypes for the major wasabi cultivars and the Ulleung Island wasabi. Two major groups (1-'Mazuma' and 'Daruma', and 2-'Fujidaruma'/'Shimane No. 3'/Ulleung Island wasabi/five cultivars in Korea) were also identified based on mother line genealogical history. Furthermore, different types of variations (mutations, insertions/deletions (indels), mononucleotide repeats, and inversions) in plastomes were identified to distinguish different cultivar lines and five highly divergent hotspots. The nine newly obtained complete plastomes are valuable organelle genomic resources for species identification and infraspecific phylogeographic studies on wild and cultivated wasabi.

The newly assembled chloroplast genome of Aeluropus littoralis: molecular feature characterization and phylogenetic analysis with related species

Aeluropus littoralis, a halophyte grass, is widely distributed from the Mediterranean to the Indian subcontinent through the Mongolian Gobi. This model halophyte has garnered increasing attention owing to its use as forage and its high tolerance to environmental stressors. The chloroplast genomes of many plants have been extensively examined for molecular, phylogenetic and transplastomic applications. However, no published research on the A. littoralis chloroplast (cp) genome was discovered. Here, the entire chloroplast genome of A. littoralis was assembled implementing accurate long-read sequences. The entire chloroplast genome, with an estimated length of 135,532 bp (GC content: 38.2%), has a quadripartite architecture and includes a pair of inverted repeat (IR) regions, IRa and IRb (21,012 bp each), separated by a large and a small single-copy regions (80,823 and 12,685 bp, respectively). The features of A. littoralis consist of 133 genes that synthesize 87 peptides, 38 transfer RNAs, and 8 ribosomal RNAs. Of these genes, 86 were unique, whereas 19 were duplicated in IR regions. Additionally, a total of forty-six simple sequence repeats, categorized into 32-mono, four-di, two-tri, and eight-tetranucleotides, were discovered. Furthermore, ten sets of repeats greater than 20 bp were located primarily in the LSC region. Evolutionary analysis based on chloroplast sequence data revealed that A. littoralis with A. lagopoides and A. sinensis belong to the Aeluropodinae subtribe, which is a sister to the Eleusininae in the tribe Cynodonteae and the subfamily Chloridoideae. This subfamily belongs to the PACMAD clade, which contains the majority of the C4 photosynthetic plants in the Poaceae. The newly constructed A. littoralis cp genome offers valuable knowledge for DNA barcoding, phylogenetic, transplastomic research, and other biological studies.

Authentication of Allium ulleungense, A. microdictyon and A. ochotense based on super-barcoding of plastid genome and 45S nrDNA

Allium ulleungense (AU) and A. microdictyon (AM) are valuable medicinal and edible vegetables, referred to as mountain garlic in Korea. The identification of AU, AM and a neighboring species A. ochotense (AO) is difficult because of their morphological similarities. We collected samples from three species and 46 cultivated collections to understand the genetic diversity of these valuable Allium species. Among them, we sequenced six collections, including three species and three cultivating collections to obtain data from the plastid genome (plastome) and nuclear 45S ribosomal DNA (nrDNA) for super-barcoding. The AM and AO showed around 60 single nucleotide polymorphisms (SNPs) and 39 Insertion/Deletion (InDels) in the plastome but no variations in the nrDNA sequences. Conversely, the AU and AM showed more than 170 SNPs and 80 InDels in the plastomes, and 20 SNPs and 1 InDel were found in the 45S nrDNA sequences. Among the three cultivating collections, one TB collection was determined to be the AU type in both plastome and nrDNA sequences. However, the other two collections, JB and SA, showed the AM type plastome but were heterozygous in the 45S nrDNA sequences, indicating both AU and AM types (putative AM x AU hybrid). Ten molecular markers were developed based on sequence variations to identify these three species and assess their genetic diversity. A total of 49 collections were genotyped using the ten developed markers and classified into five groups: 14 AU, 22 AM, 1 AO, 3 putative AM x AU hybrids, and 9 putative AU x AM hybrid collections. Super-barcoding with plastomes and nrDNAs revealed the genetic diversity of the three Allium species and putative hybrids between species. The newly developed markers will facilitate species and hybrid identification, thereby benefiting marker-assisted molecular breeding of Allium species.

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.

Comparative analyses of Linderniaceae plastomes, with implications for its phylogeny and evolution

Introduction

The recently established Linderniaceae, separated from the traditionally defined Scrophulariaceae, is a taxonomically complicated family. Although previous phylogenetic studies based on a few short DNA markers have made great contributions to the taxonomy of Linderniaceae, limited sampling and low resolution of the phylogenetic tree have failed to resolve controversies between some generic circumscriptions. The plastid genome exhibits a powerful ability to solve phylogenetic relationships ranging from shallow to deep taxonomic levels. To date, no plastid phylogenomic studies have been carried out in Linderniaceae.

Methods

In this study, we newly sequenced 26 plastid genomes of Linderniaceae, including eight genera and 25 species, to explore the phylogenetic relationships and genome evolution of the family through plastid phylogenomic and comparative genomic analyses.

Results

The plastid genome size of Linderniaceae ranged from 152,386 bp to 154,402 bp, exhibiting a typical quartile structure. All plastomes encoded 114 unique genes, comprising 80 protein-coding genes, 30 tRNA genes, and four rRNA genes. The inverted repeat regions were more conserved compared with the single-copy regions. A total of 1803 microsatellites and 1909 long sequence repeats were identified, and five hypervariable regions (petN-psbM, rps16-trnQ, rpl32-trnL, rpl32, and ycf1) were screened out. Most protein-coding genes were relatively conserved, with only the ycf2 gene found under positive selection in a few species. Phylogenomic analyses confirmed that Linderniaceae was a distinctive lineage and revealed that the presently circumscribed Vandellia and Torenia were non-monophyletic.

Discussion

Comparative analyses showed the Linderniaceae plastomes were highly conservative in terms of structure, gene order, and gene content. Combining morphological and molecular evidence, we supported the newly established Yamazakia separating from Vandellia and the monotypic Picria as a separate genus. These findings provide further evidence to recognize the phylogenetic relationships among Linderniaceae and new insights into the evolution of the plastid genomes.

Phylogenomics and plastome evolution of Indigofera (Fabaceae)

Introduction

Indigofera L. is the third largest genus in Fabaceae and includes economically important species that are used for indigo dye-producing, medicinal, ornamental, and soil and water conservation. The genus is taxonomically difficult due to the high level of overlap in morphological characters of interspecies, fewer reliability states for classification, and extensive adaptive evolution. Previous characteristic-based taxonomy and nuclear ITS-based phylogenies have contributed to our understanding of Indigofera taxonomy and evolution. However, the lack of chloroplast genomic resources limits our comprehensive understanding of the phylogenetic relationships and evolutionary processes of Indigofera.

Methods

Here, we newly assembled 18 chloroplast genomes of Indigofera. We performed a series of analyses of genome structure, nucleotide diversity, phylogenetic analysis, species pairwise Ka/Ks ratios, and positive selection analysis by combining with allied species in Papilionoideae.

Results and discussion

The chloroplast genomes of Indigofera exhibited highly conserved structures and ranged in size from 157,918 to 160,040 bp, containing 83 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Thirteen highly variable regions were identified, of which trnK-rbcL, ndhF-trnL, and ycf1 were considered as candidate DNA barcodes for species identification of Indigofera. Phylogenetic analysis using maximum likelihood (ML) and Bayesian inference (BI) methods based on complete chloroplast genome and protein-coding genes (PCGs) generated a well-resolved phylogeny of Indigofera and allied species. Indigofera monophyly was strongly supported, and four monophyletic lineages (i.e., the Pantropical, East Asian, Tethyan, and Palaeotropical clades) were resolved within the genus. The species pairwise Ka/Ks ratios showed values lower than 1, and 13 genes with significant posterior probabilities for codon sites were identified in the positive selection analysis using the branch-site model, eight of which were associated with photosynthesis. Positive selection of accD suggested that Indigofera species have experienced adaptive evolution to selection pressures imposed by their herbivores and pathogens. Our study provided insight into the structural variation of chloroplast genomes, phylogenetic relationships, and adaptive evolution in Indigofera. These results will facilitate future studies on species identification, interspecific and intraspecific delimitation, adaptive evolution, and the phylogenetic relationships of the genus Indigofera.

Comparative plastome genomics, taxonomic delimitation and evolutionary divergences of Tetraena hamiensis var. qatarensis and Tetraena simplex (Zygophyllaceae)

The Zygophyllum and Tetraena genera are intriguingly important ecologically and medicinally. Based on morphological characteristics, T. hamiensis var. qatarensis, and T. simplex were transferred from Zygophyllum to Tetraena with the least genomic datasets available. Hence, we sequenced the T. hamiensis and T. simplex and performed in-depth comparative genomics, phylogenetic analysis, and estimated time divergences. The complete plastomes ranged between 106,720 and 106,446 bp-typically smaller than angiosperms plastomes. The plastome circular genomes are divided into large single-copy regions (~ 80,964 bp), small single-copy regions (~ 17,416 bp), and two inverted repeats regions (~ 4170 bp) in both Tetraena species. An unusual shrinkage of IR regions 16-24 kb was identified. This resulted in the loss of 16 genes, including 11 ndh genes which encode the NADH dehydrogenase subunits, and a significant size reduction of Tetraena plastomes compared to other angiosperms. The inter-species variations and similarities were identified using genome-wide comparisons. Phylogenetic trees generated by analyzing the whole plastomes, protein-coding genes, matK, rbcL, and cssA genes exhibited identical topologies, indicating that both species are sisters to the genus Tetraena and may not belong to Zygophyllum. Similarly, based on the entire plastome and proteins coding genes datasets, the time divergence of Zygophyllum and Tetraena was 36.6 Ma and 34.4 Ma, respectively. Tetraena stem ages were 31.7 and 18.2 Ma based on full plastome and protein-coding genes. The current study presents the plastome as a distinguishing and identification feature among the closely related Tetraena and Zygophyllum species. It can be potentially used as a universal super-barcode for identifying plants.

Phylogenomics, phylogeography and germplasms authentication of the Rheum palmatum complex based on complete chloroplast genomes

As a traditional Chinese medicine, rhubarb is used to treat several diseases such as severe acute pancreatitis, sepsis and chronic renal failure. However, few studies focused on the authentication of germplasm for the Rheum palmatum complex, and no studies have been conducted to elucidate the evolutionary history of the R. palmatum complex using plastome datasets. Hence, we aim to develop the potential molecular markers to identify the elite germplasms of rhubarb and explore the divergence and biogeographic history of the R. palmatum complex based on the newly sequenced chloroplast genome datasets. Chloroplast genomes of thirty-five the R. palmatum complex germplasms were sequenced, and the length ranged from 160,858 to 161,204 bp. The structure, gene content and gene order were highly conserved across all genomes. Eight InDels and sixty-one SNPs loci could be used to authenticate the high-quality germplasms of rhubarb in specific areas. Phylogenetic analysis revealed that all rhubarb germplasms were clustered in the same clade with high bootstrap support values and Bayesian posterior probabilities. According to the molecular dating result, the intraspecific divergence of the complex occurred in the Quaternary, which might be affected by climatic fluctuation. The biogeography reconstruction indicated that the ancestor of the R. palmatum complex might originate from the Himalaya-Hengduan Mountains or/and Bashan-Qinling Mountains, and then spread to surrounding areas. Several useful molecular markers were developed to identify rhubarb germplasms, and our study will provide further understanding on speciation, divergence and biogeography of the R. palmatum complex.

Comparison of plastid genomes and ITS of two sister species in Gentiana and a discussion on potential threats for the endangered species from hybridization

BACKGROUND

Gentiana rigescens Franchet is an endangered medicinal herb from the family Gentianaceae with medicinal values. Gentiana cephalantha Franchet is a sister species to G. rigescens possessing similar morphology and wider distribution. To explore the phylogeny of the two species and reveal potential hybridization, we adopted next-generation sequencing technology to acquire their complete chloroplast genomes from sympatric and allopatric distributions, as along with Sanger sequencing to produce the nrDNA ITS sequences.

RESULTS

The plastid genomes were highly similar between G. rigescens and G. cephalantha. The lengths of the genomes ranged from 146,795 to 147,001 bp in G. rigescens and from 146,856 to 147,016 bp in G. cephalantha. All genomes consisted of 116 genes, including 78 protein-coding genes, 30 tRNA genes, four rRNA genes and four pseudogenes. The total length of the ITS sequence was 626 bp, including six informative sites. Heterozygotes occurred intensively in individuals from sympatric distribution. Phylogenetic analysis was performed based on chloroplast genomes, coding sequences (CDS), hypervariable sequences (HVR), and nrDNA ITS. Analysis based on all the datasets showed that G. rigescens and G. cephalantha formed a monophyly. The two species were well separated in phylogenetic trees using ITS, except for potential hybrids, but were mixed based on plastid genomes. This study supports that G. rigescens and G. cephalantha are closely related, but independent species. However, hybridization was confirmed to occur frequently between G. rigescens and G. cephalantha in sympatric distribution owing to the lack of stable reproductive barriers. Asymmetric introgression, along with hybridization and backcrossing, may probably lead to genetic swamping and even extinction of G. rigescens.

CONCLUSION

G. rigescens and G. cephalantha are recently diverged species which might not have undergone stable post-zygotic isolation. Though plastid genome shows obvious advantage in exploring phylogenetic relationships of some complicated genera, the intrinsic phylogeny was not revealed because of matrilineal inheritance here; nuclear genomes or regions are hence crucial for uncovering the truth. As an endangered species, G. rigescens faces serious threats from both natural hybridization and human activities; therefore, a balance between conservation and utilization of the species is extremely critical in formulating conservation strategies.

The Curious Case of Fritillaria sonnikovae (Liliaceae) in South Siberia: New Insights into Its Origin and Phylogeny

Fritillaria Tourn. ex L. is a genus of Liliaceae including a little more than 150 species occurring in the temperate Holarctic. Fritillaria sonnikovae Shaulo & Erst is the most recently described Siberian species in the genus. In the affinity section of the F. sonnikovae diagnosis, only F. dagana Turcz. and F. roylei Hook. are mentioned. Our study is an original attempt to shed light on the F. sonnikovae origin and its evolutionary relationships with other Fritillaria using nuclear (ITS) and plastid (matK + rps16 + trnH-psbA) DNA markers. Our results showed that F. sonnikovae together with F. dagana and F. maximowiczii Freyn belongs to the North Asian lineage of the Liliorhiza subgenus and produced no evidence supporting relationship between F. sonnikovae and F. roylei. Monophyly of Fritillaria sonnikovae was not reliably confirmed in our study since its close affinity with F. maximowiczii was demonstrated by phylogenetic analysis and morphology. Fritillaria dagana was shown to be a sister to the F. maximowiczii + F. sonnikovae group. Most authors of the present study suggest considering F. sonnikovae as a synonym for F. maximowiczii. In this view, F. sonnikovae may be considered a narrow endemic and one of the light-perianth morphs of F. maximowiczii, which has emerged in the Western Sayan and remained there as a tertiary relict.

The plastid genome of twenty-two species from Ferula, Talassia, and Soranthus: comparative analysis, phylogenetic implications, and adaptive evolution

BACKGROUND

The Ferula genus encompasses 180-185 species and is one of the largest genera in Apiaceae, with many of Ferula species possessing important medical value. The previous studies provided more information for Ferula, but its infrageneric relationships are still confusing. In addition, its genetic basis of its adaptive evolution remains poorly understood. Plastid genomes with more variable sites have the potential to reconstruct robust phylogeny in plants and investigate the adaptive evolution of plants. Although chloroplast genomes have been reported within the Ferula genus, few studies have been conducted using chloroplast genomes, especially for endemic species in China.

RESULTS

Comprehensively comparative analyses of 22 newly sequenced and assembled plastomes indicated that these plastomes had highly conserved genome structure, gene number, codon usage, and repeats type and distribution, but varied in plastomes size, GC content, and the SC/IR boundaries. Thirteen mutation hotspot regions were detected and they would serve as the promising DNA barcodes candidates for species identification in Ferula and related genera. Phylogenomic analyses with high supports and resolutions showed that Talassia transiliensis and Soranthus meyeri were nested in the Ferula genus, and thus they should be transferred into the Ferula genus. Our phylogenies also indicated the monophyly of subgenera Sinoferula and subgenera Narthex in Ferula genus. Twelve genes with significant posterior probabilities for codon sites were identified in the positively selective analysis, and their function may relate to the photosystem II, ATP subunit, and NADH dehydrogenase. Most of them might play an important role to help Ferula species adapt to high-temperatures, strong-light, and drought habitats.

CONCLUSION

Plastome data is powerful and efficient to improve the support and resolution of the complicated Ferula phylogeny. Twelve genes with significant posterior probabilities for codon sites were helpful for Ferula to adapt to the harsh environment. Overall, our study supplies a new perspective for comprehending the phylogeny and evolution of Ferula.

Complete chloroplast genomes of 11 Sabia samples: Genomic features, comparative analysis, and phylogenetic relationship

The genus Sabia is a woody climber belonging to the family Sabiaceae, order Proteales. Several species of this genus have been utilized as medicines for treating diseases, such as rheumatic arthritis, traumatism, hepatitis, etc. However, the lack of molecular data has prevented the accurate identification and refinement of taxonomic relationships in this genus. In this study, chloroplast genomes of 11 samples of the genus Sabia were assembled and analyzed. These chloroplast genomes showed a typical quadripartite structure and ranged in length from 160,956 to 162,209 bp. The structure of the genomes was found to be relatively conserved, with 130 genes annotated, including 85 coding genes, 37 tRNA genes, and eight rRNA genes. A total of 78-98 simple sequence repeats and 52-61 interspersed repeats were detected. Sequence alignment revealed 11 highly variable loci in chloroplast genomes. Among these loci, ndhF-ndhD achieved a remarkably higher resolution than the other regions. In addition, phylogenetic analysis indicated that Sect. Pachydiscus and Sect. Sabia of Sabia did not form two separate monophyletic groups. The divergence time calculated based on the Reltime method indicated that the evolutionary branches of Sabia and Meliosma started to form approximately 85.95 million years ago (Mya), and the species within Sabia began to diverge approximately 7.65 Mya. In conclusion, our study provides a basis for comprehensively exploring the phylogenetic relationships of Sabia. It also provides a methodological basis and data support for establishing a standardized and scientific identification system for this genus.

Super-Mitobarcoding in Plant Species Identification? It Can Work! The Case of Leafy Liverworts Belonging to the Genus Calypogeia

Molecular identification of species is especially important where traditional taxonomic methods fail. The genus Calypogeia belongs to one of the tricky taxons. The simple morphology of these species and a tendency towards environmental plasticity make them complicated in identification. The finding of the universal single-locus DNA barcode in plants seems to be 'the Holy Grail'; therefore, researchers are increasingly looking for multiloci DNA barcodes or super-barcoding. Since the mitochondrial genome has low sequence variation in plants, species delimitation is usually based on the chloroplast genome. Unexpectedly, our research shows that super-mitobarcoding can also work! However, our outcomes showed that a single method of molecular species delimitation should be avoided. Moreover, it is recommended to interpret the results of molecular species delimitation alongside other types of evidence, such as ecology, population genetics or comparative morphology. Here, we also presented genetic data supporting the view that C. suecica is not a homogeneous species.

Characterization of the Dicranostigma leptopodum chloroplast genome and comparative analysis within subfamily Papaveroideae

BACKGROUND

Dicranostigma leptopodum (Maxim.) Fedde is a perennial herb with bright yellow flowers, well known as "Hongmao Cao" for its medicinal properties, and is an excellent early spring flower used in urban greening. However, its molecular genomic information remains largely unknown. Here, we sequenced and analyzed the chloroplast genome of D. leptopodum to discover its genome structure, organization, and phylogenomic position within the subfamily Papaveroideae.

RESULTS

The chloroplast genome size of D. leptopodum was 162,942 bp, and D. leptopodum exhibited a characteristic circular quadripartite structure, with a large single-copy (LSC) region (87,565 bp), a small single-copy (SSC) region (18,759 bp) and a pair of inverted repeat (IR) regions (28,309 bp). The D. leptopodum chloroplast genome encoded 113 genes, including 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. The dynamics of the genome structures, genes, IR contraction and expansion, long repeats, and single sequence repeats exhibited similarities, with slight differences observed among the eight Papaveroideae species. In addition, seven interspace regions and three coding genes displayed highly variable divergence, signifying their potential to serve as molecular markers for phylogenetic and species identification studies. Molecular evolution analyses indicated that most of the genes were undergoing purifying selection. Phylogenetic analyses revealed that D. leptopodum formed a clade with the tribe Chelidonieae.

CONCLUSIONS

Our study provides detailed information on the D. leptopodum chloroplast genome, expanding the available genomic resources that may be used for future evolution and genetic diversity studies.

Species discrimination in Schima (Theaceae): Next-generation super-barcodes meet evolutionary complexity

Plastid genome and nrDNA arrays, proposed recently as "super barcodes", might provide additional discriminatory power and overcome the limitations of traditional barcoding loci, yet super barcodes need to be tested for their effectiveness in more plant groups. Morphological homoplasy among Schima species makes the genus a model for testing the efficacy of super barcodes. In this study, we generated multiple datasets comprising standard DNA barcodes (matK, rbcL, trnH-psbA, nrITS) and super-barcodes (plastid genome, nrDNA arrays) across 58 individuals from 12 out of 13 species of Schima from China. No samples were correctly assigned to species using standard DNA barcodes and nrDNA arrays, while only 27.27% of species with multiple accessions were distinguished using the plastid genome and its partitioned datasets-the lowest estimated rate of super barcode success in the literature so far. For Schima and other taxa with similarly recently divergence and low levels of genetic variation, incomplete lineage sorting, hybridization, or taxonomic oversplitting are all possible causes of the failure. Taken together, our study suggests that by no means are super barcodes immune to the challenges imposed by evolutionary complexity. We therefore call for developing multi-locus nuclear markers for species discrimination in plant groups.

Characterization of the Complete Chloroplast Genome Sequence of the Socotra Dragon`s Blood Tree (Dracaena cinnabari Balf.)

The Socotra dragon`s blood tree (Dracaena cinnabari Balf.) is endemic to the island of Socotra in Yemen. This iconic species plays an essential role in the survival of associated organisms, acting as an umbrella tree. Overexploitation, overgrazing by livestock, global climate change, and insufficient regeneration mean that the populations of this valuable species are declining in the wild. Although there are many studies on the morphology, anatomy, and physiology of D. cinnabari, no genomic analysis of this endangered species has been performed so far. Therefore, the main aim of this study was to characterize the complete chloroplast sequence genome of D. cinnabari for conservation purposes. The D. cinnabari chloroplast genome is 155,371 bp with a total GC content of 37.5%. It has a quadripartite plastid genome structure composed of one large single-copy region of 83,870 bp, one small single-copy region of 18,471 bp, and two inverted repeat regions of 26,515 bp each. One hundred and thirty-two genes were annotated, 86 of which are protein-coding genes, 38 are transfer RNAs, and eight are ribosomal RNAs. Forty simple sequence repeats have also been identified in this chloroplast genome. Comparative analysis of complete sequences of D. cinnabari chloroplast genomes with other species of the genus Dracaena showed a very high conservativeness of their structure and organization. Phylogenetic inference showed that D. cinnabari is much closer to D. draco, D. cochinchinensis, and D. cambodiana than to D. terniflora, D. angustifolia, D. hokouensis, and D. elliptica. The results obtained in this study provide new and valuable omics data for further phylogenetic studies of the genus Dracaena as well as enable the protection of genetic resources of highly endangered D. cinnabari.