Highly variable chloroplast markers for evaluating plant phylogeny at low taxonomic levels and for DNA barcoding

Characterization, comparison, and phylogenetic analyses of chloroplast genomes of Euphorbia species

The genus Euphorbia (Euphorbiaceae) has near-cosmopolitan distribution and serves as a significant resource for both ornamental and medicinal purposes. Despite its economic importance, Euphorbia's taxonomy has long been challenged by the intricate nature of morphological traits exhibiting high levels of convergence. While molecular markers are essential for phylogenetic studies, their availability for Euphorbia has been limited. To address this gap, we conducted comparative analyses focusing on the chloroplast (CP) genomes of nine Euphorbia species, incorporating three newly sequenced and annotated accessions. In addition, phylogenetic informativeness and nucleotide diversity were computed to identify candidate markers for phylogenetic analyses among closely related taxa in the genus. Our investigation revealed relatively conserved sizes and structures of CP genomes across the studied species, with notable interspecific variations observed primarily in non-coding regions and IR/SC borders. By leveraging phylogenetic informativeness and nucleotide diversity, we identified rpoB gene as the optimal candidate for species delimitation and shallow-level phylogenetic inference within the genus. Through this comprehensive analysis of CP genomes across multiple taxa, our study sheds light on the evolutionary dynamics and taxonomic intricacies of Euphorbia, offering valuable insights into its CP genome evolution and taxonomy.

Phylogenomic analysis and dynamic evolution of chloroplast genomes of Clematis nannophylla

Clematis nannophylla is a perennial shrub of Clematis with ecological, ornamental, and medicinal value, distributed in the arid and semi-arid areas of northwest China. This study successfully determined the chloroplast (cp) genome of C. nannophylla, reconstructing a phylogenetic tree of Clematis. This cp genome is 159,801 bp in length and has a typical tetrad structure, including a large single-copy, a small single-copy, and a pair of reverse repeats (IRa and IRb). It contains 133 unique genes, including 89 protein-coding, 36 tRNA, and 8 rRNA genes. Additionally, 66 simple repeat sequences, 50 dispersed repeats, and 24 tandem repeats were found; many of the dispersed and tandem repeats were between 20-30 bp and 10-20 bp, respectively, and the abundant repeats were located in the large single copy region. The cp genome was relatively conserved, especially in the IR region, where no inversion or rearrangement was observed, further revealing that the coding regions were more conserved than the noncoding regions. Phylogenetic analysis showed that C. nannophylla is more closely related to C. fruticosa and C. songorica. Our analysis provides reference data for molecular marker development, phylogenetic analysis, population studies, and cp genome processes to better utilise C. nannophylla.

Complete chloroplast genomes of eight Delphinium taxa (Ranunculaceae) endemic to Xinjiang, China: insights into genome structure, comparative analysis, and phylogenetic relationships

BACKGROUND

Delphinium L. represents a taxonomically intricate genus of significant phylogenetic and economic importance in Ranunculaceae. Despite the existence of few chloroplast genome datasets, a comprehensive understanding of genome structures and selective pressures within the genus remains unknown. Furthermore, several taxa in this genus are exclusively found in Xinjiang, China, a region renowned for its distribution and diversity of Chinese and Central Asian Delphinium species. Therefore, investigating the features of chloroplast genomes in this area will provide valuable insights into the evolutionary processes and phylogenetic relationships of the genus.

RESULTS

In this study, the eight newly completed chloroplast genomes are examined, ranging in length from 153,979 bp to 154,284 bp. Alongside these, analysing six previously reported taxa re-annotated in Delphinium, 111 unique genes are identified across all samples. Genome structure, distributions of simple sequence repeats and short dispersed repeats, as well as gene content are similar among these Delphinium taxa. Nine hypervariable intergenic spacers and protein coding regions, including ndhF-trnL(TAG), rpl16-intron, rpl33, rps15, rps18, trnK(TTT)-trnQ(TTG), trnP(TGG)-psaJ, trnT(GGT)-psbD and ycf1, are identified among 13 perennial Delphinium. Selective pressure and codon usage bias of all the plastid genes are performed within 14 Delphinium taxa. Phylogenetic analysis based on 14 Delphinium plastomes, alongside two Aconitum (Ranunculaceae) species serving as outgroup taxa, reveals the monophyletic nature of Delphinium. Our findings further discern Delphinium into two distinct clades: perennial species (clade I) and annual species (clade II). In addition, compared with the nrDNA ITS topology, cytological data and morphological characters, D. mollifolium and D. maackianum showed potential involvement in hybridization or polyploidization processes. Excluding these two species, the perennial Delphinium (clade I) exhibits a stronger consistency with the morphology-based system that utilized seed morphology.

CONCLUSION

This study represents the first comprehensive analysis of plastomic variations among Delphinium taxa, based on the examination of 14 complete plastomes. The chloroplast genome structure of Delphinium is similar to other angiosperms and possesses the typical quadripartite structure with the conserved genome arrangement and gene features. In addition, the variation of non-coding regions is larger than coding regions of the chloroplast genome. Through DNA sequence divergence across Delphinium plastomes and subsequent phylogenomic analyses ndhF-trnL(TAG) and ycf1 are identified as promising molecular markers. These highly variable loci held significant potential for future phylogenetic and phylogeographic studies on Delphinium. Our phylogenomic analyses based on the whole plastomes, concatenation of 132 unique intergenic spacer regions, concatenation of 77 unique protein-coding genes and nrDNA ITS, all support the monophyly of Delphinium and perennial taxa clusters together into one clade within this genus. These findings provide crucial data for systematic, phylogenomic and evolutionary research in the genus for future studies.

DNA barcoding combined with high-resolution melting analysis to discriminate rhubarb species and its traditional Chinese patent medicines

Introduction: Rhubarb is a frequently used and beneficial traditional Chinese medicine. Wild resources of these plants are constantly being depleted, meaning that rhubarb products have been subjected to an unparalleled level of adulteration. Consequentially, reliable technology is urgently required to verify the authenticity of rhubarb raw materials and commercial botanical drugs. Methods: In this study, the barcode-DNA high-resolution melting (Bar-HRM) method was applied to characterize 63 rhubarb samples (five Polygonaceae species: Rheum tanguticum, Rh. palmatum, Rh. officinale, Rumex japonicus and Ru. sp.) and distinguish the rhubarb contents of 24 traditional Chinese patent medicine (TCPM) samples. Three markers, namely ITS2, rbcL and psbA-trnH, were tested to assess the candidate DNA barcodes for their effectiveness in distinguishing rhubarb from its adulterants. A segment from ITS2 was selected as the most suitable mini-barcode to identify the botanical drug rhubarb in TCPMs. Then, rhubarbs and TCPM samples were subjected to HRM analysis based on the ITS2 barcode. Results: Among the tested barcoding loci, ITS2 displayed abundant sites of variation and was effective in identifying Polygonaceae species and their botanical origins. HRM analysis based on the ITS2 mini-barcode region successfully distinguished the authenticity of five Polygonaceae species and eight batches of TCPMs. Of the 18 TCPM samples, 66.7 % (12 samples) were identified as containing Rh. tanguticum or Rh. officinale. However, 33.3 % were shown to consist of adulterants. Conclusions: These results demonstrated that DNA barcoding combined with HRM is a specific, suitable and powerful approach for identifying rhubarb species and TCPMs, which is crucial to guaranteeing the security of medicinal plants being traded internationally.

Phylogeny, Genetic Diversity and Population Structure of Fritillaria cirrhosa and Its Relatives Based on Chloroplast Genome Data

Fritillaria cirrhosa and its relatives have been utilized in traditional Chinese medicine for many years and are under priority protection in China. Despite their medicinal and protective value, research on their phylogeny, genetic diversity, and divergence remains limited. Here, we investigate the chloroplast genome variation architecture of 46 samples of F. cirrhosa and its relatives collected from various regions, encompassing the majority of wild populations across diverse geographical areas. The results indicate abundant variations in 46 accessions including 1659 single-nucleotide polymorphisms and 440 indels. Six variable markers (psbJ, ndhD, ycf1, ndhG, trnT-trnL, and rpl32-trnL) were identified. Phylogenetic and network analysis, population structure analysis, and principal component analysis showed that the 46 accessions formed five clades with significant divergence, which were related to their geographical distribution. The regions spanning from the southern Hengduan Mountains to the Qinghai-Tibet Plateau exhibited the highest levels of genetic diversity. F. cirrhosa and its relatives may have suffered a genetic bottleneck and have a relatively low genetic diversity level. Moreover, geographical barriers and discrete patches may have accelerated population divergence. The study offers novel perspectives on the phylogeny, genetic diversity, and population structure of F. cirrhosa and its relatives, information that can inform conservation and utilization strategies in the future.

Comparative Analysis of Complete Chloroplast Genomes of Rubus in China: Hypervariable Regions and Phylogenetic Relationships

With more than 200 species of native Rubus, China is considered a center of diversity for this genus. Due to a paucity of molecular markers, the phylogenetic relationships for this genus are poorly understood. In this study, we sequenced and assembled the plastomes of 22 out of 204 Chinese Rubus species (including varieties) from three of the eight sections reported in China, i.e., the sections Chamaebatus, Idaeobatus, and Malachobatus. Plastomes were annotated and comparatively analyzed with the inclusion of two published plastomes. The plastomes of all 24 Rubus species were composed of a large single-copy region (LSC), a small single-copy region (SSC), and a pair of inverted repeat regions (IRs), and ranged in length from 155,464 to 156,506 bp. We identified 112 unique genes, including 79 protein-coding genes, 29 transfer RNAs, and four ribosomal RNAs. With highly consistent gene order, these Rubus plastomes showed strong collinearity, and no significant changes in IR boundaries were noted. Nine divergent hotspots were identified based on nucleotide polymorphism analysis: trnH-psbA, trnK-rps16, rps16-trnQ-psbK, petN-psbM, trnT-trnL, petA-psbJ, rpl16 intron, ndhF-trnL, and ycf1. Based on whole plastome sequences, we obtained a clearer phylogenetic understanding of these Rubus species. All sampled Rubus species formed a monophyletic group; however, sections Idaeobatus and Malachobatus were polyphyletic. These data and analyses demonstrate the phylogenetic utility of plastomes for systematic research within Rubus.

The complete chloroplast genome of Illicium simonsii Maxim. (Illiciaceae), a species with important medicinal properties

Illicium simonsii Maxim (1888) is a medicinal species of the genus Illicium in the Illiciaceae family. It is commonly used to cure gastro-frigid vomiting, cystic hernia, gas pains in the chest, and scabies as folk medicine. To utilize its resources efficiently, the complete chloroplast genome of I. simonsii was sequenced, assembled, and annotated by using high-throughput sequencing data. The complete chloroplast genome was 143,038 bp in length, with a large single-copy region (LSC) of 101,094 bp, a short single-copy region (SSC) of 20,070 bp, and a pair of inverted repeats (IRs) of 21,874 bp. A total of 113 genes were annotated, including 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. The phylogenetic tree exhibited that I. simonsii and Illicium burmanicum form a sister group, and were nested in the monophyletic clade of the Illicium genus.

The Phylogeography of Deciduous Tree Ulmus macrocarpa (Ulmaceae) in Northern China

Disentangling how climate oscillations and geographical events significantly influence plants' genetic architecture and demographic history is a central topic in phytogeography. The deciduous ancient tree species Ulmus macrocarpa is primarily distributed throughout Northern China and has timber and horticultural value. In the current study, we studied the phylogenic architecture and demographical history of U. macrocarpa using chloroplast DNA with ecological niche modeling. The results indicated that the populations' genetic differentiation coefficient (NST) value was significantly greater than the haplotype frequency (GST) (p < 0.05), suggesting that U. macrocarpa had a clear phylogeographical structure. Phylogenetic inference showed that the putative chloroplast haplotypes could be divided into three groups, in which the group Ⅰ was considered to be ancestral. Despite significant genetic differentiation among these groups, gene flow was detected. The common ancestor of all haplotypes was inferred to originate in the middle-late Miocene, followed by the haplotype overwhelming diversification that occurred in the Quaternary. Combined with demography pattern and ecological niche modeling, we speculated that the surrounding areas of Shanxi and Inner Mongolia were potential refugia for U. macrocarpa during the glacial period in Northern China. Our results illuminated the demography pattern of U. macrocarpa and provided clues and references for further population genetics investigations of precious tree species distributed in Northern China.

Phylogenetic relationships, selective pressure and molecular markers development of six species in subfamily Polygonoideae based on complete chloroplast genomes

The subfamily Polygonoideae encompasses a diverse array of medicinal and horticultural plants that hold significant economic value. However, due to the lack of a robust taxonomy based on phylogenetic relationships, the classification within this family is perplexing, and there is also a scarcity of reports on the chloroplast genomes of many plants falling under this classification. In this study, we conducted a comprehensive analysis by sequencing and characterizing the complete chloroplast genomes of six Polygonoideae plants, namely Pteroxygonum denticulatum, Pleuropterus multiflorus, Pleuropterus ciliinervis, Fallopia aubertii, Fallopia dentatoalata, and Fallopia convolvulus. Our findings revealed that these six plants possess chloroplast genomes with a typical quadripartite structure, averaging 162,931 bp in length. Comparative chloroplast analysis, codon usage analysis, and repetitive sequence analysis demonstrated a high level of conservation within the chloroplast genomes of these plants. Furthermore, phylogenetic analysis unveiled a distinct clade occupied by P. denticulatum, while P. ciliinrvis displayed a closer relationship to the three plants belonging to the Fallopia genus. Selective pressure analysis based on maximum likelihood trees showed that a total of 14 protein-coding genes exhibited positive selection, with psbB and ycf1 having the highest number of positive amino acid sites. Additionally, we identified four molecular markers, namely petN-psbM, psal-ycf4, ycf3-trnS-GGA, and trnL-UAG-ccsA, which exhibit high variability and can be utilized for the identification of these six plants.

The first complete chloroplast genome of Thalictrum fargesii: insights into phylogeny and species identification

Introduction

Thalictrum fargesii is a medicinal plant belonging to the genus Thalictrum of the Ranunculaceae family and has been used in herbal medicine in the Himalayan regions of China and India. This species is taxonomically challenging because of its morphological similarities to other species within the genus. Thus, herbal drugs from this species are frequently adulterated, substituted, or mixed with other species, thereby endangering consumer safety.

Methods

The present study aimed to sequence and assemble the entire chloroplast (cp) genome of T. fargesii using the Illumina HiSeq 2500 platform to better understand the genomic architecture, gene composition, and phylogenetic relationships within the Thalictrum.

Results and discussion

The cp genome was 155,929 bp long and contained large single-copy (85,395 bp) and small single-copy (17,576 bp) regions that were segregated by a pair of inverted repeat regions (26,479 bp) to form a quadripartite structure. The cp genome contains 133 genes, including 88 protein-coding genes (PCGs), 37 tRNA genes, and 8 rRNA genes. Additionally, this genome contains 64 codons that encode 20 amino acids, the most preferred of which are alanine and leucine. We identified 68 SSRs, 27 long repeats, and 242 high-confidence C-to-U RNA-editing sites in the cp genome. Moreover, we discovered seven divergent hotspot regions in the cp genome of T. fargesii, among which ndhD-psaC and rpl16-rps3 may be useful for developing molecular markers for identifying ethnodrug species and their contaminants. A comparative study with eight other species in the genus revealed that pafI and rps19 had highly variable sites in the cp genome of T. fargesii. Additionally, two special features, (i) the shortest length of the ycf1 gene at the IRA-SSC boundary and (ii) the distance between the rps19 fragment and trnH at the IRA-LSC junction, distinguish the cp genome of T. fargesii from those of other species within the genus. Furthermore, phylogenetic analysis revealed that T. fargesii was closely related to T. tenue and T. petaloidium.

Conclusion

Considering all these lines of evidence, our findings offer crucial molecular and evolutionary information that could play a significant role in further species identification, evolution, and phylogenetic studies on T. fargesii.

Complete Chloroplast Genome of Alternanthera sessilis and Comparative Analysis with Its Congeneric Invasive Weed Alternanthera philoxeroides

Alternanthera sessilis is considered the closest relative to the invasive weed Alternanthera philoxeroides in China, making it an important native species for studying the invasive mechanisms and adaptations of A. philoxeroides. Chloroplasts play a crucial role in a plant's environmental adaptation, with their genomes being pivotal in the evolution and adaptation of both invasive and related species. However, the chloroplast genome of A. sessilis has remained unknown until now. In this study, we sequenced and assembled the complete chloroplast genome of A. sessilis using high-throughput sequencing. The A. sessilis chloroplast genome is 151,935 base pairs long, comprising two inverted repeat regions, a large single copy region, and a small single copy region. This chloroplast genome contains 128 genes, including 8 rRNA-coding genes, 37 tRNA-coding genes, 4 pseudogenes, and 83 protein-coding genes. When compared to the chloroplast genome of the invasive weed A. philoxeroides and other Amaranthaceae species, we observed significant variations in the ccsA, ycf1, and ycf2 regions in the A. sessilis chloroplast genome. Moreover, two genes, ccsA and accD, were found to be undergoing rapid evolution due to positive selection pressure. The phylogenetic trees were constructed for the Amaranthaceae family, estimating the time of independent species formation between A. philoxeroides and A. sessilis to be approximately 3.5186-8.8242 million years ago. These findings provide a foundation for understanding the population variation within invasive species among the Alternanthera genus.

Twelve newly assembled jasmine chloroplast genomes: unveiling genomic diversity, phylogenetic relationships and evolutionary patterns among Oleaceae and Jasminum species

BACKGROUND

Jasmine (Jasminum), renowned for its ornamental value and captivating fragrance, has given rise to numerous species and accessions. However, limited knowledge exists regarding the evolutionary relationships among various Jasminum species.

RESULTS

In the present study, we sequenced seven distinct Jasminum species, resulting in the assembly of twelve high-quality complete chloroplast (cp) genomes. Our findings revealed that the size of the 12 cp genomes ranged from 159 to 165 kb and encoded 134-135 genes, including 86-88 protein-coding genes, 38-40 tRNA genes, and 8 rRNA genes. J. nudiflorum exhibited a larger genome size compared to other species, mainly attributed to the elevated number of forward repeats (FRs). Despite the typically conservative nature of chloroplasts, variations in the presence or absence of accD have been observed within J. sambac. The calculation of nucleotide diversity (Pi) values for 19 cp genomes indicated that potential mutation hotspots were more likely to be located in LSC regions than in other regions, particularly in genes ycf2, rbcL, atpE, ndhK, and ndhC (Pi > 0.2). Ka/Ks values revealed strong selection pressure on the genes rps2, atpA, rpoA, rpoC1, and rpl33 when comparing J. sambac with the three most closely related species (J. auriculatum, J. multiflorum, and J. dichotomum). Additionally, SNP identification, along with the results of Structure, PCA, and phylogenetic tree analyses, divided the Jasminum cp genomes into six groups. Notably, J. polyanthum showed gene flow signals from both the G5 group (J. nudiflorum) and the G3 group (J. tortuosum and J. fluminense). Phylogenetic tree analysis reflected that most species from the same genus clustered together with robust support in Oleaceae, strongly supporting the monophyletic nature of cp genomes within the genus Jasminum.

CONCLUSION

Overall, this study provides comprehensive insights into the genomic composition, variation, and phylogenetic relationships among various Jasminum species. These findings enhance our understanding of the genetic diversity and evolutionary history of Jasminum.

Chloroplast genome structure analysis of Equisetum unveils phylogenetic relationships to ferns and mutational hotspot region

Equisetum is one of the oldest extant group vascular plants and is considered to be the key to understanding vascular plant evolution. Equisetum is distributed almost all over the world and has a high degree of adaptability to different environments. Despite the fossil record of horsetails (Equisetum, Equisetaceae) dating back to the Carboniferous, the phylogenetic relationship of this genus is not well, and the chloroplast evolution in Equisetum remains poorly understood. In order to fill this gap, we sequenced, assembled, and annotated the chloroplast genomes of 12 species of Equisetum, and compared them to 13 previously published vascular plants chloroplast genomes to deeply examine the plastome evolutionary dynamics of Equisetum. The chloroplast genomes have a highly conserved quadripartite structure across the genus, but these chloroplast genomes have a lower GC content than other ferns. The size of Equisetum plastomes ranges from 130,773 bp to 133,684 bp and they encode 130 genes. Contraction/expansion of IR regions and the number of simple sequences repeat regions underlie large genomic variations in size among them. Comparative analysis revealed we also identified 13 divergence hotspot regions. Additionally, the genes accD and ycf1 can be used as potential DNA barcodes for the identification and phylogeny of the genus Equisetum. Twelve photosynthesis-related genes were specifically selected in Equisetum. Comparative genomic analyses implied divergent evolutionary patterns between Equisetum and other ferns. Phylogenomic analyses and molecular dating revealed a relatively distant phylogenetic relationship between Equisetum and other ferns, supporting the division of pteridophyte into Lycophytes, Equisetaceae and ferns. The results show that the chloroplast genome can be used to solve phylogenetic problems within or between Equisetum species, and also provide genomic resources for the study of Equisetum systematics and evolution.

Chloroplast genome of Calamus tetradactylus revealed rattan phylogeny

BACKGROUND

Calamus tetradactylus, a species primarily distributed in Vietnam, Laos, and southern China, is highly valued for its utilization as a small-diameter rattan material. While its physical and mechanical properties have been extensively studied, the genomic characteristics of C. tetradactylus remain largely unexplored.

RESULTS

To gain a better understanding of its chloroplast genomic features and evolutionary relationships, we conducted sequencing and assembly of the chloroplast genome of C. tetradactylus. The complete chloroplast genome exhibited the typical highly conserved quartile structure, with specific variable regions identified in the single-copy region (like psbF-psbE, π = 0.10327, ndhF-rpl32, π = 0.10195), as well as genes such as trnT-GGU (π = 0.05764) and ycf1 (π = 0.03345) and others. We propose that these regions and genes hold potential as markers for species identification. Furthermore, phylogenetic analysis revealed that C. tetradactylus formed a distinct clade within the phylogenetic tree, alongside other Calamus species, and C. tetradactylus was most closely related to C. walkeri, providing support for the monophyly of the genus.

CONCLUSION

The analysis of the chloroplast genome conducted in this study provides valuable insights that can contribute to the improvement of rattan breeding programs and facilitate sustainable development in the future.

Development of an allele-specific PCR (AS-PCR) method for identifying high-methyl eugenol-containing purple Tulsi (Ocimum tenuiflorum L.) in market samples

BACKGROUND

Ocimum tenuiflorum L. is a highly traded medicinal with several therapeutic values. Green Tulsi and purple Tulsi are two subtypes in O. tenuiflorum and both have the same medicinal properties. Recent reports have revealed that purple Tulsi contains higher quantities of methyl eugenol (ME), which is moderately toxic and potentially carcinogenic. Therefore, we developed an allele-specific PCR (AS-PCR) method to distinguish the green and purple Tulsi.

METHODS AND RESULT

Using the green Tulsi as a reference, 12 single nucleotide polymorphisms (SNPs) and 10 insertions/deletions (InDels) were identified in the chloroplast genome of the purple Tulsi. The C > T SNP at the 1,26,029 position in the ycf1 gene was selected for the development of the AS-PCR method. The primers were designed to amplify 521 bp and 291 bp fragments specific to green and purple Tulsi, respectively. This AS-PCR method was validated in 10 accessions from each subtype and subsequently verified using Sanger sequencing. Subsequently, 30 Tulsi powder samples collected from the market were subjected to molecular identification by AS-PCR. The results showed that 80% of the samples were purple Tulsi, and only 3.5% were green Tulsi. About 10% of the samples were a mixture of both green and purple Tulsi. Two samples (6.5%) did not contain O. tenuiflorum and were identified as O. gratissimum.

CONCLUSION

The market samples of Tulsi were predominantly derived from purple Tulsi. The AS-PCR method will be helpful for quality control and market surveillance of Tulsi herbal powders.

Three complete chloroplast genomes from two north American Rhus species and phylogenomics of Anacardiaceae

BACKGROUND

The suamc genus Rhus (sensu stricto) includes two subgenera, Lobadium (ca. 25 spp.) and Rhus (ca. 10 spp.). Their members, R. glabra and R. typhina (Rosanae: Sapindales: Anacardiaceae), are two economic important species. Chloroplast genome information is of great significance for the study of plant phylogeny and taxonomy.

RESULTS

The three complete chloroplast genomes from two Rhus glabra and one R. typhina accessions were obtained with a total of each about 159k bp in length including a large single-copy region (LSC, about 88k bp), a small single-copy regions (SSC, about 19k bp) and a pair of inverted repeats regions (IRa/IRb, about 26k bp), to form a canonical quadripartite structure. Each genome contained 88 protein-coding genes, 37 transfer RNA genes, eight ribosomal RNA genes and two pseudogenes. The overall GC content of the three genomes all were same (37.8%), and RSCU values showed that they all had the same codon prefers, i.e., to use codon ended with A/U (93%) except termination codon. Three variable hotspots, i.e., ycf4-cemA, ndhF-rpl32-trnL and ccsA-ndhD, and a total of 152-156 simple sequence repeats (SSR) were identified. The nonsynonymous (Ka)/synonymous (Ks) ratio was calculated, and cemA and ycf2 genes are important indicators of gene evolution. The phylogenetic analyses of the family Anacardiaceae showed that the eight genera were grouped into three clusters, and supported the monophyly of the subfamilies and all the genera. The accessions of five Rhus species formed four clusters, while, one individual of R. typhina grouped with the R. glabra accessions instead of clustering into the two other individuals of R. typhina in the subgenus Rhus, which showed a paraphyletic relationship.

CONCLUSIONS

Comparing the complete chloroplast genomes of the Rhus species, it was found that most SSRs were A/T rich and located in the intergenic spacer, and the nucleotide divergence exhibited higher levels in the non-coding region than in the coding region. The Ka/Ks ratio of cemA gene was > 1 for species collected in America, while it was < 1 for other species in China, which dedicated that the Rhus species from North America and East Asia have different evolutionary pressure. The phylogenetic analysis of the complete chloroplast genome clarified the Rhus placement and relationship. The results obtained in this study are expected to provide valuable genetic resources to perform species identification, molecular breeding, and intraspecific diversity of the Rhus species.

Morpho-anatomical characterization and DNA barcoding of Artemesia vulgaris L

Artemisia vulgaris L. belongs to Asteraceae, is a herbal plant that has various benefits in the medical field, so that its use in the medical field can be explored optimally, the plant must be thoroughly identified. This study aims to identify A. vulgaris both in terms of descriptive morpho-anatomy and DNA barcoding using BLAST and phylogenetic tree reconstruction. The morpho-anatomical character was observed on root, stem, and leaf. DNA barcoding analysis was carried out through amplification and alignment of the rbcL and matK genes. All studies were conducted on three samples from Taman Husada (Medicinal Plant Garden) Graha Famili Surabaya, Indonesia. The anatomical slide was prepared by the paraffin method. Morphological studies revealed that the leaves of A. vulgaris both on the lower-middle part and on the upper part of the stem have differences, especially in the character of the stipules, petioles, and incisions they have. Meanwhile, from the study of anatomy, A. vulgaris has an anomocytic type of stomata and its distribution is mostly on the ventral part of the leaves. Through the BLAST process and phylogenetic tree reconstruction, the plant sequences being studied are closely related to several species of the genus Artemisia as indicated by a percentage identity above 98% and branch proximity between taxa in the reconstructed phylogenetic tree.

Demographic patterns of two related desert shrubs with overlapping distributions in response to past climate changes

Numerous studies have revealed that past geological events and climatic fluctuations had profoundly affected the genetic structure and demographic patterns of species. However, related species with overlapping ranges may have responded to such environmental changes in different ways. In this study, we compared the genetic structure and population dynamics of two typical desert shrubs with overlapping distributions in northern China, Nitraria tangutorum and Nitraria sphaerocarpa, based on chloroplast DNA (cpDNA) variations and species distribution models. We sequenced two cpDNA fragments (trnH-trnA and atpH-atpI) in 633 individuals sampled from 52 natural populations. Twenty-four chlorotypes, including eight rare chlorotypes, were identified, and a single dominant haplotype (H4) widely occurred in the entire geographical ranges of the two species. There were also a few distinctive chlorotypes fixed in different geographical regions. Population structure analyses suggested that the two species had significantly different levels of total genetic diversity and interpopulation differentiation, which was highly likely correlated with the special habitat preferences of the two species. A clear phylogeographic structure was identified to exist among populations of N. sphaerocarpa, but not exist for N. tangutorum. The neutral tests, together with the distribution of pairwise differences revealed that N. tangutorum experienced a sudden demographic expansion, and its expansion approximately occurred between 21 and 7 Kya before present, while a rapid range expansion was not identified for N. sphaerocarpa. The ecological niche modeling (ENM) analysis indicated that the potential ranges of two species apparently fluctuated during the past and present periods, with obvious contraction in the Last Glacial Maximum (LGM) and recolonization in the present, respectively, comparing to the Last Interglacial (LIG). These findings suggest that the two species extensively occurred in the Northwest of China before the Quaternary, and the current populations of them originated from a few separated glacial refugia following their habitat fragmentation in the Quarternary. Our results provide new insights on the impact of past geological and climatic fluctuations on the population dynamics of desert plants in northwestern China, and further enforce the hypothesis that there were several independent glacial refugia for these species during the Quaternary glaciations.

Complete Chloroplast Genomes of Four Oaks from the Section Cyclobalanopsis Improve the Phylogenetic Analysis and Understanding of Evolutionary Processes in the Genus Quercus

Quercus is a valuable genus ecologically, economically, and culturally. They are keystone species in many ecosystems. Species delimitation and phylogenetic studies of this genus are difficult owing to frequent hybridization. With an increasing number of genetic resources, we will gain a deeper understanding of this genus. In the present study, we collected four Quercus section Cyclobalanopsis species (Q. poilanei, Q. helferiana, Q. camusiae, and Q. semiserrata) distributed in Southeast Asia and sequenced their complete genomes. Following analysis, we compared the results with those of other species in the genus Quercus. These four chloroplast genomes ranged from 160,784 bp (Q. poilanei) to 161,632 bp (Q. camusiae) in length, with an overall guanine and cytosine (GC) content of 36.9%. Their chloroplast genomic organization and order, as well as their GC content, were similar to those of other Quercus species. We identified seven regions with relatively high variability (rps16, ndhk, accD, ycf1, psbZ-trnG-GCC, rbcL-accD, and rpl32-trnL-UAG) which could potentially serve as plastid markers for further taxonomic and phylogenetic studies within Quercus. Our phylogenetic tree supported the idea that the genus Quercus forms two well-differentiated lineages (corresponding to the subgenera Quercus and Cerris). Of the three sections in the subgenus Cerris, the section Ilex was split into two clusters, each nested in the other two sections. Moreover, Q. camusiae and Q. semiserrata detected in this study diverged first in the section Cyclobalanopsis and mixed with Q. engleriana in the section Ilex. In particular, 11 protein coding genes (atpF, ndhA, ndhD, ndhF, ndhK, petB, petD, rbcL, rpl22, ycf1, and ycf3) were subjected to positive selection pressure. Overall, this study enriches the chloroplast genome resources of Quercus, which will facilitate further analyses of phylogenetic relationships in this ecologically important tree genus.

DNA barcoding for the assessment of marine and coastal fish diversity from the Coast of Mozambique

The ichthyological provinces of Mozambique are understudied hotspots of global fish diversity. In this study, we applied DNA barcoding to identify the composition of the fish fauna from the coast of Mozambique. A total of 143 species belonging to 104 genera, 59 families, and 30 orders were identified. The overall K2P distance of the COI sequences within species ranged from 0.00% to 1.51%, while interspecific distances ranged from 3.64% to 24.49%. Moreover, the study revealed 15 threatened species according to the IUCN Red List of Threatened Species, with elasmobranchs being the most represented group. Additionally, the study also uncovered four new species that were not previously recorded in this geographic area, including Boleophthalmus dussumieri, Maculabatis gerrardi, Hippocampus kelloggi, and Lethrinus miniatus. This study represents the first instance of utilizing molecular references to explore the fish fauna along the Mozambican coast. Our results indicate that DNA barcoding is a dependable technique for the identification and delineation of fish species in the waters of Mozambique. The DNA barcoding library established in this research will be an invaluable asset for advancing the understanding of fish diversity and guiding future conservation initiatives.

The complete chloroplast genome sequence of Nepeta bracteata and comparison with congeneric species

Nepeta bracteata (N. bracteata) is an important medicinal plant used by Chinese ethnic minorities. However, the lack of knowledge regarding the chloroplast genome of N. bracteata has imposed current limitations on our study. Here, we used Next-generation sequencing to obtain the chloroplast genome of N. bracteata. The findings suggested that the 151,588 bp cp genome of N. bracteata comprises 130 genes, including 35 tRNA genes and 87 protein-coding genes. And its chloroplast genome exhibits a typical quadripartite structure, the largest single copy (LSC; 82,819 bp) and the smallest single copy (SSC; 17,557 bp) separate a pair of inverted repeats IR regions (IRa and IRb; 25,606 bp) from one another. Interestingly, palindromic repeats are more common, as shown by the examination of repetition. In the interim, 18 SSRs were discovered in the interim, the bulk of which were Adenine-Thymine (A-T) mononucleotides. Meanwhile, we compared it with five other species from the Nepeta genus. Five hypervariable areas were found by the study, including ndhH-rps15, accD-psal, ndhG-ndhl, trnH-GUG-psbA, and rpoC1-rpoB. Furthermore, the phylogenetic study revealed that N. bracteata and Nepeta stewartiana (N. stewartiana) were linked to each other most closely. In summary, our findings enrich the resources available for chloroplast genomes in the Nepeta genus. Moreover, these hypervariable regions have the potential to be developed into molecular markers, enabling the rapid identification of species within the Nepeta genus. Comparative analysis of species within the Nepeta genus can help enhance our study of their phylogenetic relationships, potential medicinal properties and bioprospecting.

An insight on the complete chloroplast genome of Gomphocarpus siniacus and Duvalia velutina, Asclepiadoideae (Apocynaceae)

Abstract We studied the complete chloroplast genome of Gomphocarpus siniacus and Duvalia velutina from Asclepiadoideae subfamily; due to their medicinal importance and distribution worldwide their interest became high. In this study we analyzed the complete chloroplast genomes of G. siniacus and D. velutina using Illumina sequencing technology. The sequences were compared with the other species from Apocynaceae family. The complete genome of G. siniacus is 162,570 bp while D. velutina has154, 478 bp in length. Both genomes consist of 119 genes; encode 31 tRNA genes, and eight rRNA genes. Comparative studies of the two genomes showed variations in SSR markers in which G. siniacus possesses 223 while D. velutina has 186. This could be used for barcoding in order to aid in easy identification of the species. Phylogenetic analysis on the other hand reaffirms the tribal position of G. siniacus in Asclepiadeae and D. velutina in Ceropegieae. These findings could be used in subsequent research studies of angiosperms identification, genetic engineering, herb genomics and phylogenomic studies of Apocynaceae family.

Evaluation of chloroplast DNA barcoding markers to individualize Papaver somniferum for forensic intelligence purposes

The opium poppy, Papaver somniferum L., is a forensically important plant due to the medicinal and illegal uses for the milky latex stored in the pods. This latex contains the alkaloids morphine, codeine, and thebaine that are used for their analgesic properties and/or for synthesizing other opioids. However, these compounds are highly addictive and have caused a national opioid epidemic. Two other Papaver species, P. setigerum DC. and P. bracteatum Lindl., are also of forensic interest because they pose both forensic and legal issues. They are largely uncontrolled under the Controlled Substances Act, making these species a common defense strategy. Current morphological and chemical identification methods have been moderately successful but have drawbacks. There is also a lack of sequencing data available. Therefore, exploiting the genome using chloroplast DNA barcoding markers could help to accurately identify these species of interest when plant material is taken. This study screened and assessed the genetic variation both between species and within populations of P. somniferum in nine cpDNA barcode regions (ndhF-rpl32, petA-psbJ, rpl32-trnL, rps16-trnQ, trnE-trnT, trnH-psbA, trnL-trnF, rpl16 intron, and psbE-petL). Published reference genomes from the NCBI GenBank database were aligned and compared for an initial in silico screening. Additionally, ten P. somniferum seed samples from various vendors were sequenced and compared across samples and to published reference data at the various barcode regions of interest. This study showed that the regions trnH-psbA and petA-psbJ have promise for utility in individualization for both inter- and intra-species individualization of P. somniferum.

Rhodophyta DNA Barcoding: Ribulose-1, 5-Bisphosphate Carboxylase Gene and Novel Universal Primers

Red algae (Rhodophyta) are a heterogeneous group of marine algal species that have served as a source of high-value molecules, including antioxidants and scaffolds, for novel drug development. However, it is challenging to identify Rhodophytes through morphological features alone, and in most instances, that has been the prevailing approach to identification. Consequently, this study undertook the identification of red algae species in Kenton-on-Sea, South Africa, as a baseline for future research on red algae biodiversity and conservation. The identification was achieved by designing, analysing, and using a set of universal primers through DNA barcoding of the rbcL gene. The PCR products of the rbcL gene were sequenced, and 96% of the amplicons were successfully sequenced from this set and matched with sequences on BOLD, which led to these species being molecularly described. Amongst these species are medicinally essential species, such as Laurencia natalensis and Hypnea spinella, and potential cryptic species. This calls for further investigation into the biodiversity of the studied region. Meanwhile, the availability of these primers will ease the identification process of red algae species from other coastal regions.

Whole mitochondrial and chloroplast genome sequencing of Tunisian date palm cultivars: diversity and evolutionary relationships

BACKGROUND

Date palm (Phoenix dactylifera L.) is the most widespread crop in arid and semi-arid regions and has great traditional and socioeconomic importance, with its fruit well-known for its high nutritional and health value. However, the genetic variation of date palm cultivars is often neglected. The advent of high-throughput sequencing has made possible the resequencing of whole organelle (mitochondria and chloroplast) genomes to explore the genetic diversity and phylogenetic relationships of cultivated plants with unprecedented detail.

RESULTS

Whole organelle genomes of 171 Tunisian accessions (135 females and 36 males) were sequenced. Targeted bioinformatics pipelines were used to identify date palm haplotypes and genome variants, aiming to provide variant annotation and investigate patterns of evolutionary relationship. Our results revealed the existence of unique haplotypes, identified by 45 chloroplastic and 156 mitochondrial SNPs. Estimation of the effect of these SNPs on genes functions was predicted in silico.

CONCLUSIONS

The results of this study have important implications, in the light of ongoing environmental changes, for the conservation and sustainable use of the genetic resources of date palm cultivars in Tunisia, where monoculture threatens biodiversity leading to genetic erosion. These data will be useful for breeding and genetic improvement programs of the date palm through selective cross-breeding.

Comparative plastome genomics and phylogenetic relationships of the genus Trollius

Trollius, a genus in the Ranunculaceae family, has significant medicinal and ornamental value. It is widely distributed in China with 16 different species accepted. However, due to the lack of enough samples and information sites, the molecular phylogenetic relationships of Trollius have been unresolved till now. Here we sequenced, assembled and annotated the plastomes of 16 Trollius species to investigate their genomic characteristics, inverted repeat (IR) boundaries, sequence repeats, and hypervariable loci. In addition, the phylogenetic relationships of this genus was reconstructed based on the whole plastomes and the protein-coding sequences data-sets. The plastomes of Trollius ranged between 159,597 bp and 160,202 bp in length, and contained 113 unique genes, including 79 protein coding, 30 tRNA, and 4 rRNA. The IR boundaries were relatively conserved within the genus Trollius. 959 simple sequence repeats and 657 long sequence repeats were detected in the Trollius plastomes. We identified 12 highly polymorphic loci (Pi > 0.0115) that can be used as plastid markers in molecular identification and phylogenetic investigation of the genus. Besides, Trollius was a monophyletic group with the earliest divergence clade being Trollius lilacinus Bunge, and the remaining species were divided into two strongly-supported clades. The phylogeny in our study supported the traditional classification systems based on the color of sepal, but not the previous classification system based on the types and relative lengths of the nectaries, and distribution. The genomic resources provided in our study can be used in the taxonomy of the genus Trollius, promoting the development and utilization of this genus.

Complete plastid genome structure of 13 Asian Justicia (Acanthaceae) species: comparative genomics and phylogenetic analyses

BACKGROUND

Justicia L. is the largest genus in Acanthaceae Juss. and widely distributed in tropical and subtropical regions of the world. Previous phylogenetic studies have proposed a general phylogenetic framework for Justicia based on several molecular markers. However, their studies were mainly focused on resolution of phylogenetic issues of Justicia in Africa, Australia and South America due to limited sampling from Asia. Additionally, although Justicia plants are of high medical and ornamental values, little research on its genetics was reported. Therefore, to improve the understanding of its genomic structure and relationships among Asian Justicia plants, we sequenced complete chloroplast (cp.) genomes of 12 Asian plants and combined with the previously published cp. genome of Justicia leptostachya Hemsl. for further comparative genomics and phylogenetic analyses.

RESULTS

All the cp. genomes exhibit a typical quadripartite structure without genomic rearrangement and gene loss. Their sizes range from 148,374 to 151,739 bp, including a large single copy (LSC, 81,434-83,676 bp), a small single copy (SSC, 16,833-17,507 bp) and two inverted repeats (IR, 24,947-25,549 bp). GC contents range from 38.1 to 38.4%. All the plastomes contain 114 genes, including 80 protein-coding genes, 30 tRNAs and 4 rRNAs. IR variation and repetitive sequences analyses both indicated that Justicia grossa C. B. Clarke is different from other Justicia species because its lengths of ndhF and ycf1 in IRs are shorter than others and it is richest in SSRs and dispersed repeats. The ycf1 gene was identified as the candidate DNA barcode for the genus Justicia. Our phylogenetic results showed that Justicia is a polyphyletic group, which is consistent with previous studies. Among them, J. grossa belongs to subtribe Tetramerinae of tribe Justicieae while the other Justicia members belong to subtribe Justiciinae. Therefore, based on morphological and molecular evidence, J. grossa should be undoubtedly recognized as a new genus. Interestingly, the evolutionary history of Justicia was discovered to be congruent with the morphology evolution.

CONCLUSION

Our study not only elucidates basic features of Justicia whole plastomes, but also sheds light on interspecific relationships of Asian Justicia plants for the first time.

Comparison of chloroplast genomes and phylogenetic analysis of four species in Quercus section Cyclobalanopsis

The identification in Quercus L. species was considered to be difficult all the time. The fundamental phylogenies of Quercus have already been discussed by morphological and molecular means. However, the morphological characteristics of some Quercus groups may not be consistent with the molecular results (such as the group Helferiana), which may lead to blurring of species relationships and prevent further evolutionary researches. To understand the interspecific relationships and phylogenetic positions, we sequenced and assembled the CPGs (160,715 bp-160842 bp) of four Quercus section Cyclobalanopsis species by Illumina pair-end sequencing. The genomic structure, GC content, and IR/SC boundaries exhibited significant conservatism. Six highly variable hotspots were detected in comparison analysis, among which rpoC1, clpP and ycf1 could be used as molecular markers. Besides, two genes (petA, ycf2) were detected to be under positive selection pressure. The phylogenetic analysis showed: Trigonobalanus genus and Fagus genus located at the base of the phylogeny tree; The Quercus genus species were distincted to two clades, including five sections. All Compound Trichome Base species clustered into a single branch, which was in accordance with the results of the morphological studies. But neither of group Gilva nor group Helferiana had formed a monophyly. Six Compound Trichome Base species gathered together in pairs to form three branch respectively (Quercus kerrii and Quercus chungii; Quercus austrocochinchinensis with Quercus gilva; Quercus helferiana and Quercus rex). Due to a low support rate (0.338) in the phylogeny tree, the interspecies relationship between the two branches differentiated by this node remained unclear. We believe that Q. helferiana and Q. kerrii can exist as independent species due to their distance in the phylogeny tree. Our study provided genetic information in Quercus genus, which could be applied to further studies in taxonomy and phylogenetics.

Comparative analysis of chloroplast genome and new insights into phylogenetic relationships of Ajuga and common adulterants

Introduction

The potential contamination of herbal medicinal products poses a significant concern for consumer health. Given the limited availability of genetic information concerning Ajuga species, it becomes imperative to incorporate supplementary molecular markers to enhance and ensure accurate species identification.

Methods

In this study, the chloroplast (cp) genomes of seven species of the genus Ajuag were sequenced, de novo assembled and characterized.

Results

exhibiting lengths ranging from 150,342 bp to 150,472 bp, encompassing 86 - 88 protein-coding genes (PCGs), 35 - 37 transfer RNA, and eight ribosomal RNA. The repetitive sequences, codon uses, and cp genomes of seven species were highly conserved, and PCGs were the reliable molecular markers for investigating the phylogenetic relationship within the Ajuga genus. Moreover, four mutation hotspot regions (accD-psaI, atpH-atpI, ndhC-trnV(UAC), and ndhF-rpl23) were identified within cp genomes of Ajuga, which could help distinguish A. bracteosa and its contaminants. Based on cp genomes and PCGs, the phylogenetic tree preliminary confirmed the position of Ajuga within the Lamiaceae family. It strongly supported a sister relationship between Subsect. Genevense and Subsect. Biflorae, suggesting the merger of Subsect. Biflorae and Subsect. Genevenses into one group rather than maintaining separate categorizations. Additionally, molecular clock analysis estimated the divergence time of Ajuga to be around 7.78 million years ago.

Discussion

The species authentication, phylogeny, and evolution analyses of the Ajuga species may benefit from the above findings.

Comparative Analyses of Chloroplast Genome Provide Effective Molecular Markers for Species and Cultivar Identification in Bougainvillea

Bougainvillea is popular in ornamental horticulture for its colorful bracts and excellent adaptability, but the complex genetic relationship among this genus is fuzzy due to limited genomic data. To reveal more genomic resources of Bougainvillea, we sequenced and assembled the complete chloroplast (cp) genome sequences of Bougainvillea spectabilis 'Splendens'. The cp genome size was 154,869 bp in length, containing 86 protein-coding genes, 38 tRNAs, and eight rRNAs. Cp genome comparison across 12 Bougainvillea species (B. spectabilis, B. glabra, B. peruviana, B. arborea, B. praecox, B. stipitata, B. campanulata, B. berberidifolia, B. infesta, B. modesta, B. spinosa, and B. pachyphylla) revealed five mutational hotspots. Phylogenetic analysis suggested that B. spectabilis published previously and B. glabra clustered into one subclade as two distinct groups, sister to the subclade of B. spectabilis 'Splendens'. We considered the phylogeny relationships between B. spectabilis and B. glabra to be controversial. Based on two hypervariable regions and three common plastid regions, we developed five molecular markers for species identification in Bougainvillea and applied them to classify 53 ornamental Bougainvillea cultivars. This study provides a valuable genetic resource for Bougainvillea breeding and offers effective molecular markers to distinguish the representative ornamental species of Bougainvillea.

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.

Comparative analysis of the medicinal plant Polygonatum kingianum (Asparagaceae) with related verticillate leaf types of the Polygonatum species based on chloroplast genomes

Background

Polygonatum kingianum has been widely used as a traditional Chinese medicine as well as a healthy food. Because of its highly variable morphology, this medicinal plant is often difficult to distinguish from other related verticillate leaf types of the Polygonatum species. The contaminants in P. kingianum products not only decrease the products' quality but also threaten consumer safety, seriously inhibiting the industrial application of P. kingianum.

Methods

Nine complete chloroplast (cp) genomes of six verticillate leaf types of the Polygonatum species were de novo assembled and systematically analyzed.

Results

The total lengths of newly sequenced cp genomes ranged from 155,437 to 155,977 bp, including 86/87 protein-coding, 38 tRNA, and 8 rRNA genes, which all exhibited well-conserved genomic structures and gene orders. The differences in the IR/SC (inverted repeats/single-copy) boundary regions and simple sequence repeats were detected among the verticillate leaf types of the Polygonatum cp genomes. Comparative cp genomes analyses revealed that a higher similarity was conserved in the IR regions than in the SC regions. In addition, 11 divergent hotspot regions were selected, providing potential molecular markers for the identification of the Polygonatum species with verticillate leaf types. Phylogenetic analysis indicated that, as a super barcode, plastids realized a fast and efficient identification that clearly characterized the relationships within the verticillate leaf types of the Polygonatum species. In brief, our results not only enrich the data on the cp genomes of the genus Polygonatum but also provide references for the P. kingianum germplasm resource protection, herbal cultivation, and drug production.

Conclusion

This study not only accurately identifies P. kingianum species, but also provides valuable information for the development of molecular markers and phylogenetic analyses of the Polygonatum species with verticillate leaf types.

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.

Complete Chloroplast Genome of Corethrodendron fruticosum (Papilionoideae: Fabaceae): Comparative and Phylogenetic Analysis

Corethrodendron fruticosum is an endemic forage grasses in China with high ecological value. In this study, the complete chloroplast genome of C. fruticosum was sequenced using Illumina paired-end sequencing. The C. fruticosum chloroplast genome was 123,100 bp and comprised 105 genes, including 74 protein-coding genes, 4 rRNA-coding genes, and 27 tRNA-coding genes. The genome had a GC content of 34.53%, with 50 repetitive sequences and 63 simple repeat repetitive sequences that did not contain reverse repeats. The simple repeats included 45 single-nucleotide repeats, which accounted for the highest proportion and primarily comprised A/T repeats. A comparative analysis of C. fruticosum, C. multijugum, and four Hedysarum species revealed that the six genomes were highly conserved, with differentials primarily located in the conserved non-coding regions. Moreover, the accD and clpP genes in the coding regions exhibited high nucleotide variability. Accordingly, these genes may serve as molecular markers for the classification and phylogenetic analysis of Corethrodendron species. Phylogenetic analysis further revealed that C. fruticosum and C. multijugum appeared in different clades than the four Hedysarum species. The newly sequenced chloroplast genome provides further insights into the phylogenetic position of C. fruticosum, which is useful for the classification and identification of Corethrodendron.

Comparative chloroplast genome analysis of Sambucus L. (Viburnaceae): inference for phylogenetic relationships among the closely related Sambucus adnata Wall. ex DC Sambucus javanica Blume

Sambucus L. is found in the family Viburnaceae (syn. Adoxaceae) and encompasses approximately 29 accepted species. The complex morphology of these species has caused continued confusion concerning their nomenclature, classification, and identification. Despite previous attempts to resolve taxonomic complexities in the Sambucus genus, there are still unclear phylogenetic relationships among several species. In this study, the newly obtained plastome of Sambucus williamsii Hance. as well as the populations of Sambucus canadensis L., Sambucus javanica Blume, and Sambucus adnata Wall. ex DC were sequenced, and their sizes, structural similarity, gene order, gene number, and guanine-cytosine (GC) contents were analyzed. The phylogenetic analyses were conducted using the whole chloroplast genomes and protein-coding genes (PCGs). The findings revealed that the chloroplast genomes of Sambucus species exhibited typical quadripartite double-stranded DNA molecules. Their lengths ranged from 158,012 base pairs (bp) (S. javanica) to 158,716 bp (S. canadensis L). Each genome comprised a pair of inverted repeats (IRs), which separated the large single-copy (LSC) and small single-copy (SSC) regions. In addition, the plastomes contained 132 genes, encompassing 87 protein-coding, 37 tRNA, and four rRNA genes. In the simple sequence repeat (SSR) analysis, A/T mononucleotides had the highest proportion, with the most repetitive sequences observed in S. williamsii. The comparative genome analyses showed high similarities in structure, order, and gene contents. The hypervariable regions in the studied chloroplast genomes were trnT-GGU, trnF-GAA, psaJ, trnL-UAG, ndhF, and ndhE, which may be used as candidate barcodes for species discrimination in Sambucus genus. Phylogenetic analyses supported the monophyly of Sambucus and revealed the separation of S. javanica and S. adnata populations. Sambucus chinensis Lindl. was nested within S. javanica in the same clade, collaborating their conspecific treatment. These outcomes indicate that the chloroplast genome of Sambucus plants is a valuable genetic resource for resolving taxonomic discrepancies at the lower taxonomic levels and can be applied in molecular evolutionary studies.

Comparative chloroplast genome analysis of Sambucus L. (Viburnaceae): inference for phylogenetic relationships among the closely related Sambucus adnata Wall. ex DC Sambucus javanica Blume

Sambucus L. is found in the family Viburnaceae (syn. Adoxaceae) and encompasses approximately 29 accepted species. The complex morphology of these species has caused continued confusion concerning their nomenclature, classification, and identification. Despite previous attempts to resolve taxonomic complexities in the Sambucus genus, there are still unclear phylogenetic relationships among several species. In this study, the newly obtained plastome of Sambucus williamsii Hance. as well as the populations of Sambucus canadensis L., Sambucus javanica Blume, and Sambucus adnata Wall. ex DC were sequenced, and their sizes, structural similarity, gene order, gene number, and guanine–cytosine (GC) contents were analyzed. The phylogenetic analyses were conducted using the whole chloroplast genomes and protein-coding genes (PCGs). The findings revealed that the chloroplast genomes of Sambucus species exhibited typical quadripartite double-stranded DNA molecules. Their lengths ranged from 158,012 base pairs (bp) (S. javanica) to 158,716 bp (S. canadensis L). Each genome comprised a pair of inverted repeats (IRs), which separated the large single-copy (LSC) and small single-copy (SSC) regions. In addition, the plastomes contained 132 genes, encompassing 87 protein-coding, 37 tRNA, and four rRNA genes. In the simple sequence repeat (SSR) analysis, A/T mononucleotides had the highest proportion, with the most repetitive sequences observed in S. williamsii. The comparative genome analyses showed high similarities in structure, order, and gene contents. The hypervariable regions in the studied chloroplast genomes were trnT-GGU, trnF-GAA, psaJ, trnL-UAG, ndhF, and ndhE, which may be used as candidate barcodes for species discrimination in Sambucus genus. Phylogenetic analyses supported the monophyly of Sambucus and revealed the separation of S. javanica and S. adnata populations. Sambucus chinensis Lindl. was nested within S. javanica in the same clade, collaborating their conspecific treatment. These outcomes indicate that the chloroplast genome of Sambucus plants is a valuable genetic resource for resolving taxonomic discrepancies at the lower taxonomic levels and can be applied in molecular evolutionary studies.

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.

Pueraria montana Population Structure and Genetic Diversity Based on Chloroplast Genome Data

Despite having a generally conserved structure, chloroplast genome data have been helpful for plant population genetics and evolution research. To mine Pueraria montana chloroplast genome variation architecture and phylogeny, we investigated the chloroplast variation architecture of 104 P. montana accessions from across China. P. montana's chloroplast genome showed high diversity levels, with 1674 variations, including 1118 single nucleotide polymorphisms and 556 indels. The intergenic spacers, psbZ-trnS and ccsA-ndhD, are the two mutation hotspot regions in the P. montana chloroplast genome. Phylogenetic analysis based on the chloroplast genome dataset supported four P. montana clades. P. montana variations were conserved among and within clades, which showed high gene flow levels. Most P. montana clades were estimated to have diverged at 3.82-5.17 million years ago. Moreover, the East Asian summer monsoon and South Asian summer monsoon may have accelerated population divergence. Our results show that chloroplast genome sequences were highly variable and can be used as molecular markers to assess genetic variation and relationships in P. montana.

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

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

Intraspecific Chloroplast Genome Variation and Domestication Origins of Major Cultivars of Styphnolobium japonicum

Styphnolobium japonicum is a significant resource of ornamental and medicinal plants. In this study, we employed high-throughput sequencing to assemble nine chloroplast genomes of S. japonicum. We compared and reconstructed the phylogenetic relationships of these genomes, along with three publicly available chloroplast genomes. Our results showed that the length of the 12 S. japonicum chloroplast genomes ranged from 158,613 bp to 158,837 bp, all containing 129 unique functional genes. The genetic diversity within S. japonicum chloroplast genomes was relatively low, with π = 0.00029, Theta-W = 0.00028, and an indel frequency of 0.62 indels/1 kb. Among the four regions, the SSC region exhibited the highest genetic diversity and indel frequency, while the IR region had the lowest. Non-coding regions displayed greater genetic variation compared to coding regions, with a few highly variable regions identified. The phylogenetic tree constructed revealed that the major cultivars of S. japonicum originated from two genetic 'sources. S. japonicum 'JinhuaiJ2' had an independent origin and showed close relatedness to S. japonicum var. violacea, S. japonicum var. japonicum, and S. japonicum f. oligophylla. On the other hand, other major cultivars shared a common genetic origin and were closely related to S. japonicum f. pendula. This study highlights the variability of chloroplast genomes within S. japonicum and provides insights into the genetic origins of major cultivars and their relationships with different varieties and forma.

Characterization and comparative analysis of the complete plastid genomes of four Astragalus species

Astragalus is the largest flowering plant genus. We assembled the plastid genomes of four Astragalus species (Astragalus iranicus, A. macropelmatus, A. mesoleios, A. odoratus) using next-generation sequencing and analyzed their plastomes including genome organization, codon usage, nucleotide diversity, prediction of RNA editing and etc. The total length of the newly sequenced Astragalus plastomes ranged from 121,050 bp to 123,622 bp, with 110 genes comprising 76 protein-coding genes, 30 transfer RNA (tRNA) genes and four ribosome RNA (rRNA) genes. Comparative analysis of the chloroplast genomes of Astragalus revealed several hypervariable regions comprising three non-coding sites (trnQ(UUG)-accD, rps7 -trnV(GAC) and trnR(ACG)-trnN(GUU)) and four protein-coding genes (ycf1, ycf2, accD and clpP), which have potential as molecular markers. Positive selection signatures were found in five genes in Astragalus species including rps11, rps15, accD, clpP and ycf1. The newly sequenced species, A. macropelmatus, has an approximately 13-kb inversion in IR region. Phylogenetic analysis based on 75 protein-coding gene sequences confirmed that Astragalus form a monophyletic clade within the tribe Galegeae and Oxytropis is sister group to the Coluteoid clade. The results of this study may helpful in elucidating the chloroplast genome structure, understanding the evolutionary dynamics at genus Astragalus and IRLC levels and investigating the phylogenetic relationships. Moreover, the newly plastid genomes sequenced have been increased the plastome data resources on Astragalus that can be useful in further phylogenomic studies.

Chloroplast Pan-Genomes and Comparative Transcriptomics Reveal Genetic Variation and Temperature Adaptation in the Cucumber

Although whole genome sequencing, genetic variation mapping, and pan-genome studies have been done on a large group of cucumber nuclear genomes, organelle genome information is largely unclear. As an important component of the organelle genome, the chloroplast genome is highly conserved, which makes it a useful tool for studying plant phylogeny, crop domestication, and species adaptation. Here, we have constructed the first cucumber chloroplast pan-genome based on 121 cucumber germplasms, and investigated the genetic variations of the cucumber chloroplast genome through comparative genomic, phylogenetic, haplotype, and population genetic structure analysis. Meanwhile, we explored the changes in expression of cucumber chloroplast genes under high- and low-temperature stimulation via transcriptome analysis. As a result, a total of 50 complete chloroplast genomes were successfully assembled from 121 cucumber resequencing data, ranging in size from 156,616-157,641 bp. The 50 cucumber chloroplast genomes have typical quadripartite structures, consisting of a large single copy (LSC, 86,339-86,883 bp), a small single copy (SSC, 18,069-18,363 bp), and two inverted repeats (IRs, 25,166-25,797 bp). Comparative genomic, haplotype, and population genetic structure results showed that there is more genetic variation in Indian ecotype cucumbers compared to other cucumber cultivars, which means that many genetic resources remain to be explored in Indian ecotype cucumbers. Phylogenetic analysis showed that the 50 cucumber germplasms could be classified into 3 types: East Asian, Eurasian + Indian, and Xishuangbanna + Indian. The transcriptomic analysis showed that matK were significantly up-regulated under high- and low-temperature stresses, further demonstrating that cucumber chloroplasts respond to temperature adversity by regulating lipid metabolism and ribosome metabolism. Further, accD has higher editing efficiency under high-temperature stress, which may contribute to the heat tolerance. These studies provide useful insight into genetic variation in the chloroplast genome, and established the foundation for exploring the mechanisms of temperature-stimulated chloroplast adaptation.

Phylogenetic relationships and biogeography of Asia Callicarpa (Lamiaceae), with consideration of a long-distance dispersal across the Pacific Ocean -insights into divergence modes of pantropical flora

There are about 140 species of Callicarpa L. 1753 (Lamiaceae), with more species richness in tropical to subtropical Asia and the New World. The genus might provide an insight into the amphi-Pacific disjunction pattern of tropical and subtropical vegetation. This study has greatly improved the phylogenetic underpinning for Callicarpa, derived from more inclusive taxonomic samplings, and employing data on both two-nuclear and eight-chloroplast regions. To address time and patterns of diversification in Callicarpa, we conducted divergence time and biogeographic analyses, and inferred shifts in the distribution areas across the phylogenetic clades. Our phylogenetic results show that Callicarpa is monophyletic with respect to the groups considered, and eight well-supported primary clades were discerned in the combined analyses. Our estimates indicated that the crown group of Callicarpa originates around the Late-Eocene (ca. 36.23 Ma) and diversification within most clades is concentrated in the Miocene and continued to the Pleistocene. In addition, our biogeographic analyses suggested that the probable ancestor of the Callicarpa crown clade originated in East Asia and Southeast Asia. Multiple dispersal and vicariance events contributed to the current distribution of the taxa. Furthermore, this genus expanded eastward out of East and Southeast Asia to the New World by long-distance dispersal, which inspired us to better understand the amphi-Pacific disjunct distribution.

Comparative analysis of dioecious Amaranthus plastomes and phylogenomic implications within Amaranthaceae s.s

BACKGROUND

The genus Amaranthus L. consists of 70-80 species distributed across temperate and tropical regions of the world. Nine species are dioecious and native to North America; two of which are agronomically important weeds of row crops. The genus has been described as taxonomically challenging and relationships among species including the dioecious ones are poorly understood. In this study, we investigated the phylogenetic relationships among the dioecious amaranths and sought to gain insights into plastid tree incongruence. A total of 19 Amaranthus species' complete plastomes were analyzed. Among these, seven dioecious Amaranthus plastomes were newly sequenced and assembled, an additional two were assembled from previously published short reads sequences and 10 other plastomes were obtained from a public repository (GenBank).

RESULTS

Comparative analysis of the dioecious Amaranthus species' plastomes revealed sizes ranged from 150,011 to 150,735 bp and consisted of 112 unique genes (78 protein-coding genes, 30 transfer RNAs and 4 ribosomal RNAs). Maximum likelihood trees, Bayesian inference trees and splits graphs support the monophyly of subgenera Acnida (7 dioecious species) and Amaranthus; however, the relationship of A. australis and A. cannabinus to the other dioecious species in Acnida could not be established, as it appears a chloroplast capture occurred from the lineage leading to the Acnida + Amaranthus clades. Our results also revealed intraplastome conflict at some tree branches that were in some cases alleviated with the use of whole chloroplast genome alignment, indicating non-coding regions contribute valuable phylogenetic signals toward shallow relationship resolution. Furthermore, we report a very low evolutionary distance between A. palmeri and A. watsonii, indicating that these two species are more genetically related than previously reported.

CONCLUSIONS

Our study provides valuable plastome resources as well as a framework for further evolutionary analyses of the entire Amaranthus genus as more species are sequenced.

The First Complete Chloroplast Genome of Cordia monoica: Structure and Comparative Analysis

Cordia monoica is a member of the Boraginaceae family. This plant is widely distributed in tropical regions and has a great deal of medical value as well as economic importance. In the current study, the complete chloroplast (cp) genome of C. monoica was sequenced, assembled, annotated, and reported. This circular chloroplast genome had a size of 148,711 bp, with a quadripartite structure alternating between a pair of repeated inverted regions (26,897-26,901 bp) and a single copy region (77,893 bp). Among the 134 genes encoded by the cp genome, there were 89 protein-coding genes, 37 transfer RNA (tRNA) genes, and 8 ribosomal RNA (rRNA) genes. A total of 1387 tandem repeats were detected, with the hexanucleotides class making up 28 percent of the repeats. Cordia monoica has 26,303 codons in its protein-coding regions, and leucine amino acid was the most frequently encoded amino acid in contrast to cysteine. In addition, 12 of the 89 protein-coding genes were found to be under positive selection. The phyloplastomic taxonomical clustering of the Boraginaceae species provides further evidence that chloroplast genome data are reliable not only at family level but also in deciphering the phylogeny at genus level (e.g., Cordia).

Gene Losses and Homology of the Chloroplast Genomes of Taxillus and Phacellaria Species

Research on the chloroplast genome of parasitic plants is limited. In particular, the homology between the chloroplast genomes of parasitic and hyperparasitic plants has not been reported yet. In this study, three chloroplast genomes of Taxillus (Taxillus chinensis, Taxillus delavayi, and Taxillus thibetensis) and one chloroplast genome of Phacellaria (Phacellaria rigidula) were sequenced and analyzed, among which T. chinensis is the host of P. rigidula. The chloroplast genomes of the four species were 119,941-138,492 bp in length. Compared with the chloroplast genome of the autotrophic plant Nicotiana tabacum, all of the ndh genes, three ribosomal protein genes, three tRNA genes and the infA gene were lost in the three Taxillus species. Meanwhile, in P. rigidula, the trnV-UAC gene and the ycf15 gene were lost, and only one ndh gene (ndhB) existed. The results of homology analysis showed that the homology between P. rigidula and its host T. chinensis was low, indicating that P. rigidula grows on its host T. chinensis but they do not share the chloroplast genome. In addition, horizontal gene transfer was not found between P. rigidula and its host T. chinensis. Several candidate highly variable regions in the chloroplast genomes of Taxillus and Phacellaria species were selected for species identification study. Phylogenetic analysis revealed that the species of Taxillus and Scurrula were closely related and supported that Scurrula and Taxillus should be treated as congeneric, while species in Phacellaria had a close relationship with that in Viscum.

New insights into the phylogenetic relationships among wild onions (Allium, Amaryllidaceae), with special emphasis on the subgenera Anguinum and Rhizirideum, as revealed by plastomes

The genus Allium, with over 900 species, is one of the largest monocotyledonous genera and is widely accepted with 15 recognized subgenera and 72 sections. The robust subgeneric and sectional relationships within Allium have long been not resolved. Based on 76 species of Allium (a total of 84 accessions), we developed a highly resolved plastome phylogenetic framework by integrating 18 newly sequenced species (20 accessions) in this study and assessed their subgeneric and sectional relationships, with special emphasis on the two subgenera Anguinum and Rhizirideum. We retrieved the three major evolutionary lines within Allium and found that the two subgenera Anguinum and Rhizirideum are monophyletic whereas others are highly polyphyletic (e.g., Allium, Cepa, Polyprason, and Melanocrommyum). Within the subgenus Anguinum, two strongly supported sublineages in East Asian and Eurasian-American were found. Allium tricoccum in North America belonged to the Eurasian clade. The distinct taxonomic status of A. ulleungense and its sister taxon were further determined. In subg. Rhizirideum, the Ulleung Island endemic A. dumebuchum shared its most recent common ancestor with the species from Mongolia and the narrow Korean endemic A. minus. Two Ulleung Island endemics were estimated to originate independently during the Pleistocene. In addition, a separate monotypic sectional treatment of the east Asian A. macrostemon (subg. Allium) and sister relationship between A. condensatum and A. chinense was suggested.

Multiple domestication events explain the origin of Gossypium hirsutum landraces in Mexico

Several Mesoamerican crops constitute wild-to-domesticated complexes generated by multiple initial domestication events, and continuous gene flow among crop populations and between these populations and their wild relatives. It has been suggested that the domestication of cotton (Gossypium hirsutum) started in the northwest of the Yucatán Peninsula, from where it spread to other regions inside and outside of Mexico. We tested this hypothesis by assembling chloroplast genomes of 23 wild, landraces, and breeding lines (transgene-introgressed and conventional). The phylogenetic analysis showed that the evolutionary history of cotton in Mexico involves multiple events of introgression and genetic divergence. From this, we conclude that Mexican landraces arose from multiple wild populations. Our results also revealed that their structural and functional chloroplast organizations had been preserved. However, genetic diversity decreases as a consequence of domestication, mainly in transgene-introgressed (TI) individuals (π = 0.00020, 0.00001, 0.00016, 0, and 0, of wild, TI-wild, landraces, TI-landraces, and breeding lines, respectively). We identified homologous regions that differentiate wild from domesticated plants and indicate a relationship among the samples. A decrease in genetic diversity associated with transgene introgression in cotton was identified for the first time, and our outcomes are therefore relevant to both biosecurity and agrobiodiversity conservation.