The complete chloroplast genome sequence of the medicinal plant Sophora tonkinensis

Comparative plastome analyses and evolutionary relationships of 25 East Asian species within the medicinal plant genus Scrophularia (Scrophulariaceae)

Backgroud

Scrophularia L., a genus of the Scrophulariaceae, is a group of important medicinal plants used for eliminating heat and detoxifying. East Asia has an abundance of potentially medicinal Scrophularia species, and it serves as a secondary diversity center of the genus. However, the genomic resources available for germplasm identification and pharmaceutical exploration of East Asian Scrophularia are insufficient, hindering its commercial and industrial development. Additionally, the interspecific relationships of most East Asian Scrophularia species remain unclear.

Methods

In this study, we sequenced the leaves of 25 East Asian species of the genus Scrophularia, assembled and annotated the complete chloroplast genomes, and subsequently performed comparative and phylogenetic analyses on these genomes.

Results and discussion

The conserved plastome length of these 25 species ranged from 151,582 bp to 153,239 bp, containing a total of 132 coding genes, including 18 duplicated genes and 114 unique genes. Through genome alignment of these 25 species, 38-53 repeated sequences and 7 shared SSRs were identified, along with regions with high nucleotide polymorphism (Pi), which could potentially serve as molecular markers for species identification. The genome structure, gene content, and arrangement showed conservation, while variations were observed in the IR boundary regions and IGS. Phylogenetic inferences based on whole plastomes or on coding sequences (CDS) only yielded congruent results. We categorized the 25 East Asian Scrophularia species into six distinct clades and further explored their interspecies relationships using morphological characteristics, such as flower color, the relative position of stamens and corolla, and plant height. This could lay a genetic basis for future resource development of Scrophularia in East Asia.

Deciphering the complete chloroplast genome sequence of Meconopsis torquata Prain: Insights into genome structure, comparative analysis and phylogenetic relationship

In the present study, we have characterized the complete chloroplast (Cp) genome of Meconopsis torquata Prain (family Papaveraceae), revealing the plastome size of 153,290 bp, and a GC content of 38.72 %. The cp genome features the typical circular quadripartite structure found in flowering plants, including a pair of inverted repeat regions (25,816 bp), isolated by a small single-copy region (17,740 bp) and a large single-copy (83,918 bp). Genome annotation revealed 132 genes: 87 protein-coding genes, 37 tRNAs and eight rRNAs. This comparative study demonstrated that the genome structure, gene number and GC ratio are consistent with several other cp genomes of Meconopsis and Papaver genera. A total of 120 SSRs were detected in the plastome, the majority (111) of which were mononucleotide repeats. Among the longer repeats, palindromic sequences were most common, followed by forward, reverse, and complement repeats. The whole genome alignment revealed the conserved nature of the inverted repeat region over single-copy zones. Nucleotide diversity unveiled hypervariable sites (ycf1, rps16, accD, atpB and psbD) in both the small and large single-copy regions, which could be useful for designing molecular markers for taxonomic identification. Phylogenetic analysis revealed a close alliance of M. torquata with other Meconopsis species, such as M. pinnatifolia and M. paniculata, with strong bootstrap support. Molecular dating suggests that M. torquata originated during the Tortonian age of the Miocene epoch of the Cenozoic era. These findings provide valuable insights for biological research, especially in understanding the genetic and evolutionary divergence within the Papaveraceae family.

Characterization and Comparative Analysis of Complete Chloroplast Genomes of Four Bromus (Poaceae, Bromeae) Species

Bromus (Poaceae Bromeae) is a forage grass with high adaptability and ecological and economic value. Here, we sequenced Bromus ciliatus, Bromus benekenii, Bromus riparius, and Bromus rubens chloroplast genomes and compared them with four previously described species. The genome sizes of Bromus species ranged from 136,934 bp (Bromus vulgaris) to 137,189 bp (Bromus ciliates, Bromus biebersteinii), with a typical quadripartite structure. The studied species had 129 genes, consisting of 83 protein-coding, 38 tRNA-coding, and 8 rRNA-coding genes. The highest GC content was found in the inverted repeat (IR) region (43.85-44.15%), followed by the large single-copy (LSC) region (36.25-36.65%) and the small single-copy (SSC) region (32.21-32.46%). There were 33 high-frequency codons, with those ending in A/U accounting for 90.91%. A total of 350 simple sequence repeats (SSRs) were identified, with single-nucleotide repeats being the most common (61.43%). A total of 228 forward and 141 palindromic repeats were identified. No reverse or complementary repeats were detected. The sequence identities of all sequences were very similar, especially with respect to the protein-coding and inverted repeat regions. Seven highly variable regions were detected, which could be used for molecular marker development. The constructed phylogenetic tree indicates that Bromus is a monophyletic taxon closely related to Triticum. This comparative analysis of the chloroplast genome of Bromus provides a scientific basis for species identification and phylogenetic studies.

Two new cytochalasins from the endophytic fungus Xylaria sp. GDGJ-77B

Two new open-chain cytochalasins, xylarchalasins A and B (1 and 2), together with six known analogues (3-8), were isolated from the endophytic fungus Xylaria sp. GDGJ-77B from the Chinese medicinal plant Sophora tonkinensis. Their structures were elucidated on the basis of comprehensive spectroscopic analysis. Compound 2 displayed moderate antibacterial activities against Bacillus subtilis and Escherichia coli with MIC values of 25 and 12.5 μg/mL, respectively.

Cycloaspeptide H, a cyclopentapeptide from the endophytic fungus Penicillium virgatum

One new cyclopentapeptide, cycloaspeptide H (1), featuring a serine residue, along with seven known compounds (2-8), was isolated from the endophytic fungus Penicillium virgatum GDGJ-227. The planar structure of 1 was elucidated by a comprehensive analysis of NMR and MS spectroscopic spectra, and the absolute configuration was determined by single-crystal X-ray diffraction (Cu Kα) analysis. Compounds 7 and 8 displayed antibacterial activities against Bacillus subtilis and Enterobacter aerogenes with MIC values ranging from 12.5 to 50 μg/mL.

Applications of some advanced sequencing, analytical, and computational approaches in medicinal plant research: a review

The potential active chemicals found in medicinal plants, which have long been employed as natural medicines, are abundant. Exploring the genes responsible for producing these compounds has given new insights into medicinal plant research. Previously, the authentication of medicinal plants was done via DNA marker sequencing. With the advancement of sequencing technology, several new techniques like next-generation sequencing, single molecule sequencing, and fourth-generation sequencing have emerged. These techniques enshrined the role of molecular approaches for medicinal plants because all the genes involved in the biosynthesis of medicinal compound(s) could be identified through RNA-seq analysis. In several research insights, transcriptome data have also been used for the identification of biosynthesis pathways. miRNAs in several medicinal plants and their role in the biosynthesis pathway as well as regulation of the disease-causing genes were also identified. In several research articles, an in silico study was also found to be effective in identifying the inhibitory effect of medicinal plant-based compounds against virus' gene(s). The use of advanced analytical methods like spectroscopy and chromatography in metabolite proofing of secondary metabolites has also been reported in several recent research findings. Furthermore, advancement in molecular and analytic methods will give new insight into studying the traditionally important medicinal plants that are still unexplored.

Characteristics and phylogenetic analysis of the complete chloroplast genome of Primulina hedyotidea

Primulina hedyotidea (Woon Young Chun) Yin Zheng Wang 2011 is an important medicinal plant that has a long history of medicinal use in China. In this experiment, the whole chloroplast genome of P. hedyotidea was determined by next-generation sequencing technology. The total base length of P. hedyotidea was 153,297 bp, the GC content was 37.62%, the inverted repeat (IR) region length was 25,494 bp, the large single copy (LSC) region was 84,158 bp and the small single copy (SSC) region was 18,151 bp. In addition, the genome consisted of 80 protein-coding genes, 4 rRNA genes, and 28 tRNA genes, for a total of 112 genes. A phylogenetic tree was constructed to explore the evolutionary relationship between P. hedyotidea and other species. The findings of phylogenetic tree analysis show that Primulina huaijiensis and P. hedyotidea have a close relationship, and this study can help with species identification and phylogenetic analysis within Primulina and Gesneriaceae species.

The known, unknown, and the intriguing about members of a critically endangered traditional medicinal plant genus Aconitum

Humanity will always be indebted to plants. In the ongoing scientific era, the 'Herbal Revolution' has helped discover several valuable medicinal plants and associated novel secondary metabolites from the diverse unexplored ecosystems, treating several diseases via phytotherapy. The Aconitum genus comprises several economically-important poisonous mountainous medicinal plant species whose unique biodiversity is on the verge of extinction due to illegal human intervention triggered habitat loss, over-harvesting, and unrestricted trading. Owing to its vast diversity of diterpene alkaloids, most species are extensively used to treat several ailments in rural parts of the world. Irrespective of this, many unexplored and intriguing prospects exist to understand and utilize this critical plant for human benefit. This systematic review tries to fill this gap by compiling information from the sporadically available literature known for ~300 Aconitum spp. regarding its nomenclature and classification, endangerment, plant morphology, ploidy, secondary metabolites, drug pharmacokinetics, conservation, and omics-based computational studies. We also depicted the disparity in the studied model organisms for this diverse genus. The absence of genomic/metagenomic data is becoming a limiting factor in understanding its plant physiology, metabolic pathways, and plant-microbes interactions, and therefore must be promoted. Additionally, government support and public participation are crucial in establishing conservation protocols to save this plant from endangerment.

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.

Complete chloroplast genome and phylogenetic analysis of Anemone shikokiana

BACKGROUND

Anemone shikokiana (Makino) Makino, disjunctly distributed in Shandong Peninsula of China and Shikoku Island of Japan, is a rare and endangered species. To provide genetic information and understand its phylogeny, we conducted research on the chloroplast (cp) genome of A. shikokiana.

METHODS AND RESULTS

The complete cp genome sequence of A. shikokiana was constructed in this study. The results showed that the cp genome of A. shikokiana has a typical quadripartite cyclic with a total length of 159,286 bp. In total, 111 unique genes were identified, including 78 protein-coding genes, 29 tRNA-coding genes and 4 rRNA-coding genes. A total of 37 long repeat sequences and 67 microsatellites were found in this cp genome. The cp genome of A. shikokiana was compared with eleven other Anemone cp genomes available from the Genbank database. We found some variations among the different genomes, especially in the LSC and SSC regions, and identified some regions as potential molecular markers such as ycf1, ndhE, ndhD, ndhF-trnL, ndhA and ndhF. The results of phylogenetic analysis suggested that A. narcissiflora was the closest relative of A. shikokiana.

CONCLUSIONS

The results filled the gap of cp genome sequence information of A. shikokiana, laying the foundation to explore the evolutionary relationships of A. shikokiana in future studies. It provided a valuable genetic resource for the molecular identification and phylogenetic study of Anemone.

A comparison of 25 complete chloroplast genomes between sister mangrove species Kandelia obovata and Kandelia candel geographically separated by the South China Sea

In 2003, Kandelia obovata was identified as a new mangrove species differentiated from Kandelia candel. However, little is known about their chloroplast (cp) genome differences and their possible ecological significance. In this study, 25 whole cp genomes, with seven samples of K. candel from Malaysia, Thailand, and Bangladesh and 18 samples of K. obovata from China, were sequenced for comparison. The cp genomes of both species encoded 128 genes, namely 83 protein-coding genes, 37 tRNA genes, and eight rRNA genes, but the cp genome size of K. obovata was ~2 kb larger than that of K. candle due to the presence of more and longer repeat sequences. Of these, tandem repeats and simple sequence repeats exhibited great differences. Principal component analysis based on indels, and phylogenetic tree analyses constructed with homologous protein genes from the single-copy genes, as well as 38 homologous pair genes among 13 mangrove species, gave strong support to the separation of the two species within the Kandelia genus. Homologous genes ndhD and atpA showed intraspecific consistency and interspecific differences. Molecular dynamics simulations of their corresponding proteins, NAD(P)H dehydrogenase chain 4 (NDH-D) and ATP synthase subunit alpha (ATP-A), predicted them to be significantly different in the functions of photosynthetic electron transport and ATP generation in the two species. These results suggest that the energy requirement was a pivotal factor in their adaptation to differential environments geographically separated by the South China Sea. Our results also provide clues for future research on their physiological and molecular adaptation mechanisms to light and temperature.

Comparative plastomes and phylogenetic analysis of seven Korean endemic Saussurea (Asteraceae)

BACKGROUND

Saussurea is one of the most species-rich genera in the Cardueae, Asteraceae. There are approximately 40 Saussurea species distributed in Korea, with nearly 40% of them endemics. Infrageneric relationships remain uncertain due to insufficient resolutions and low statistical support. In this study, we sequenced the plastid genomes of five Korean endemic Saussurea (S. albifolia, S. calcicola, S. diamantica, S. grandicapitula, and S. seoulensis), and comparative analyses including two other endemics (S. chabyoungsanica and S. polylepis) were conducted.

RESULTS

The plastomes of Korean endemics were highly conserved in gene content, order, and numbers. Exceptionally, S. diamantica had mitochondrial DNA sequences including two tRNAs in SSC region. There were no significant differences of the type and numbers of SSRs among the seven Korean endemics except in S. seoulensis. Nine mutation hotspots with high nucleotide diversity value (Pi > 0.0033) were identified, and phylogenetic analysis suggested that those Korean endemic species most likely evolved several times from diverse lineages within the genus. Moreover, molecular dating estimated that the Korean endemic species diverged since the late Miocene.

CONCLUSIONS

This study provides insight into understanding the plastome evolution and evolutionary relationships of highly complex species of Saussurea in Korean peninsula.

Comparative Analysis of the Chloroplast Genome of Cardamine hupingshanensis and Phylogenetic Study of Cardamine

Cardamine hupingshanensis (K. M. Liu, L. B. Chen, H. F. Bai and L. H. Liu) is a perennial herbal species endemic to China with narrow distribution. It is known as an important plant for investigating the metabolism of selenium in plants because of its ability to accumulate selenium. However, the phylogenetic position of this particular species in Cardamine remains unclear. In this study, we reported the chloroplast genome (cp genome) for the species C. hupingshanensis and analyzed its position within Cardamine. The cp genome of C. hupingshanensis is 155,226 bp in length and exhibits a typical quadripartite structure: one large single copy region (LSC, 84,287 bp), one small single copy region (17,943 bp) and a pair of inverted repeat regions (IRs, 26,498 bp). Guanine-Cytosine (GC) content makes up 36.3% of the total content. The cp genome contains 111 unique genes, including 78 protein-coding genes, 29 tRNA genes and 4 rRNA genes. A total of 115 simple sequences repeats (SSRs) and 49 long repeats were identified in the genome. Comparative analyses among 17 Cardamine species identified the five most variable regions (trnH-GUG-psbA, ndhK-ndhC, trnW-CCA-trnP-UGG, rps11-rpl36 and rpl32-trnL-UAG), which could be used as molecular markers for the classification and phylogenetic analyses of various Cardamine species. Phylogenetic analyses based on 79 protein coding genes revealed that the species C. hupingshanensis is more closely related to the species C. circaeoides. This relationship is supported by their shared morphological characteristics.

Cold Resistance of Euonymus japonicus Beihaidao Leaves and Its Chloroplast Genome Structure and Comparison with Celastraceae Species

Euonymus japonicus Beihaidao is one of the most economically important ornamental species of the Euonymus genus. There are approximately 97 genera and 1194 species of plants worldwide in this family (Celastraceae). Using E. japonicus Beihaidao, we conducted a preliminary study of the cold resistance of this species, evaluated its performance during winter, assembled and annotated its chloroplast genome, and performed a series of analyses to investigate its gene structure GC content, sequence alignment, and nucleic acid diversity. Our objectives were to understand the evolutionary relationships of the genus and to identify positive selection genes that may be related to adaptations to environmental change. The results indicated that E. japonicus Beihaidao leaves have certain cold resistance and can maintain their viability during wintering. Moreover, the chloroplast genome of E. japonicus Beihaidao is a typical double-linked ring tetrad structure, which is similar to that of the other four Euonymus species, E. hamiltonianus, E. phellomanus, E. schensianus, and E. szechuanensis, in terms of gene structure, gene species, gene number, and GC content. Compared to other Celastraceae species, the variation in the chloroplast genome sequence was lower, and the gene structure was more stable. The phylogenetic relationships of 37 species inferred that members of the Euonymus genus do not form a clade and that E. japonicus Beihaidao is closely related to E. japonicus and E. fortunei. A total of 11 functional positive selected genes were identified, which may have played an important role in the process of Celastraceae species adapting to environmental changes. Our study provides important genetic information to support further investigations into the phylogenetic development and adaptive evolution of Celastraceae species.

The complete chloroplast genome of Verbascum thapsus L. (Scrophulariaceae) and its phylogenetic affinities

Verbascum thapsus L. has extensive pharmacological effects, including antioxidative and antineoplastic action, memory improvement and neuroprotection. However, its phylogenetic position is not established in Scrophulariaceae. In this study, we reported the complete chloroplast genome sequence of V. thapsus L. for the first time and investigate its phylogenetic relationship in Scrophulariaceae. The assembled chloroplast genome is a circular 153,338 bp sequence, including a large single copy (LSC) region of 84,627 bp, a small single copy (SSC) region of 17,829 bp and a pair of inverted repeats (IRs) of 25,441 bp. The genome contains 135 genes, including 86 protein coding genes, 37 tRNA genes, and eight rRNA genes. The phylogenetic tree showed that V. thapsus is closely associated with V. chinense and V. phoeniceum.

The Complete Chloroplast Genome of Endangered Species Stemona parviflora: Insight into the Phylogenetic Relationship and Conservation Implications

Stemona parviflora is an endangered species, narrowly endemic to Hainan and Southwest Guangdong. The taxonomic classification of S. parviflora remains controversial. Moreover, studying endangered species is helpful for current management and conservation. In this study, the first complete chloroplast genome of S. parviflora was assembled and compared with other Stemona species. The chloroplast genome size of S. parviflora was 154,552 bp, consisting of 87 protein-coding genes, 38 tRNA genes, 8 rRNA genes, and one pseudogene. The ψycf1 gene was lost in the cp genome of S. sessilifolia, but it was detected in four other species of Stemona. The inverted repeats (IR) regions have a relatively lower length variation compared with the large single copy (LSC) and small single copy (SSC) regions. Long repeat sequences and simple sequence repeat (SSR) were detected, and most SSR were distributed in the LSC region. Codon usage bias analyses revealed that the RSCU value of the genus Stemona has almost no difference. As with most angiosperm chloroplast genomes, protein-coding regions were more conservative than the inter-gene spacer. Seven genes (atpI, ccsA, cemA, matK, ndhA, petA, and rpoC1) were detected under positive selection in different Stemona species, which may result from adaptive evolution to different habitats. Phylogenetic analyses show the Stemona cluster in two main groups; S. parviflora were closest to S. tuberosa. A highly suitable region of S. parviflora was simulated by Maxent in this study; it is worth noting that the whole territory of Taiwan has changed to a low fitness area and below in the 2050 s, which may not be suitable for the introduction and cultivation of S. parviflora. In addition, limited by the dispersal capacity of S. parviflora, it is necessary to carry out artificial grafts to expand the survival areas of S. parviflora. Our results provide valuable information on characteristics of the chloroplast genome, phylogenetic relationships, and potential distribution range of the endangered species S. parviflora.

Sequencing and characterization of the chloroplast genome of Aconitum forrestii Stapf provide insights into phylogenetics in Aconitum

Aconitum forrestii Stapf is an essential traditional Chinese medicine, and is beneficial in dispelling wind, removing dampness, warming, and relieving pain. However, its phylogenetic position of Aconitum is not accepted yet. In order to clarify the evolutionary relationship of A. forrestii, complete sequencing of chloroplast genome was carried out using Illumina sequencing technology. In total, the chloroplast genome was about 155,869 base pair (bp) in length and carried a typical tetrad structure that included a large single-copy, a small-single copy and two inverted repeat regions. A total of 132 genes were annotated, that included 85 protein -coding genes, 37 transfer RNA genes, eight ribosomal RNA genes, and two pseudogenes. The phylogenetic tree analysis indicated that Aconitum forrestii is closely related to Aconitum episcopale and Aconitum delavayi.

Comparative analysis of chloroplast genomes reveals phylogenetic relationships and intraspecific variation in the medicinal plant Isodon rubescens

Isodon rubescens (Hemsley) H. Hara (Lamiaceae) is a traditional Chinese medicine plant that has been used to treat various human diseases and conditions such as inflammation, respiratory and gastrointestinal bacterial infections, and malignant tumors. However, the contents of the main active components of I. rubescens from different origins differ significantly, which greatly affected its quality. Therefore, a molecular method to identify and classify I. rubescens is needed. Here, we report the DNA sequence of the chloroplast genome of I. rubescens collected from Lushan, Henan province. The genome is 152,642 bp in length and has a conserved structure that includes a pair of IR regions (25,726 bp), a LSC region (83,527 bp) and a SSC region (17,663 bp). The chloroplast genome contains 113 unique genes, four rRNA genes, 30 tRNA genes, and 79 protein-coding genes, 23 of which contain introns. The protein-coding genes account for a total of 24,412 codons, and most of them are A/T biased usage. We identified 32 simple sequence repeats (SSRs) and 48 long repeats. Furthermore, we developed valuable chloroplast molecular resources by comparing chloroplast genomes from three Isodon species, and both mVISTA and DnaSP analyses showed that rps16-trnQ, trnS-trnG, and ndhC-trnM are candidate regions that will allow the identification of intraspecific differences within I. rubescens. Also 14 candidate fragments can be used to identify interspecific differences between species in Isodon. A phylogenetic analysis of the complete chloroplast genomes of 24 species in subfamily Nepetoideae was performed using the maximum likelihood method, and shows that I. rubescens clustered closer to I. serra than I. lophanthoides. Interestingly, our analysis showed that I. rubescens (MW018469.1) from Xianyang, Shaanxi Province (IR-X), is closer to I. serra than to the other two I. rubescens accessions. These results strongly indicate that intraspecific diversity is present in I. rubescens. Therefore, our results provide further insight into the phylogenetic relationships and interspecific diversity of species in the genus Isodon.

The Complete Chloroplast Genomes of Primula obconica Provide Insight That Neither Species nor Natural Section Represent Monophyletic Taxa in Primula (Primulaceae)

The genus Primula (Primulaceae) comprises more than 500 species, with 300 species distributed in China. The contradictory results between systematic analyses and morphology-based taxonomy make taxonomy studies difficult. Furthermore, frequent introgression between closely related species of Primula can result in non-monophyletic species. In this study, the complete chloroplast genome of sixteen Primula obconica subsp. obconica individuals were assembled and compared with 84 accessions of 74 species from 21 sections of the 24 sections of the genus in China. The plastome sizes of P. obconica subsp. obconica range from 153,584 bp to 154,028 bp. Genome-wide variations were detected, and 1915 high-quality SNPs and 346 InDels were found. Most SNPs were detected in downstream and upstream gene regions (45.549% and 41.91%). Two cultivated accessions, ZP1 and ZP2, were abundant with SSRs. Moreover, 12 SSRs shared by 9 accessions showed variations that may be used as molecular markers for population genetic studies. The phylogenetic tree showed that P. obconica subsp. obconica cluster into two independent clades. Two subspecies have highly recognizable morphological characteristics, isolated geographical distribution areas, and distinct phylogenetic relationships compared with P. obconica subsp. obconica. We elevate the two subspecies of P. obconica to separate species. Our phylogenetic tree is largely inconsistent with morphology-based taxonomy. Twenty-one sections of Primula were mainly divided into three clades. The monophyly of Sect. Auganthus, Sect. Minutissimae, Sect. Sikkimensis, Sect. Petiolares, and Sect. Ranunculoides are well supported in the phylogenetic tree. The Sect. Obconicolisteri, Sect. Monocarpicae, Sect. Carolinella, Sect. Cortusoides, Sect. Aleuritia, Sect. Denticulata, Sect. Proliferae Pax, and Sect. Crystallophlomis are not a monophyletic group. The possible explanations for non-monophyly may be hybridization, polyploidization, recent introgression, incorrect taxonomy, or chloroplast capture. Multiple genomic data and population genetic studies are therefore needed to reveal the evolutionary history of Primula. Our results provided valuable information for intraspecific variation and phylogenetic relationships within Primula.

Characterization and Comparative Analysis of Complete Chloroplast Genomes of Three Species From the Genus Astragalus (Leguminosae)

Astragalus is the largest genus in Leguminosae. Several molecular studies have investigated the potential adulterants of the species within this genus; nonetheless, the evolutionary relationships among these species remain unclear. Herein, we sequenced and annotated the complete chloroplast genomes of three Astragalus species—Astragalus adsurgens, Astragalus mongholicus var. dahuricus, and Astragalus melilotoides using next-generation sequencing technology and plastid genome annotator (PGA) tool. All species belonged to the inverted repeat lacking clade (IRLC) and had similar sequences concerning gene contents and characteristics. Abundant simple sequence repeat (SSR) loci were detected, with single-nucleotide repeats accounting for the highest proportion of SSRs, most of which were A/T homopolymers. Using Astragalus membranaceus var. membranaceus as reference, the divergence was evident in most non-coding regions of the complete chloroplast genomes of these species. Seven genes (atpB, psbD, rpoB, rpoC1, trnV, rrn16, and rrn23) showed high nucleotide variability (Pi), and could be used as DNA barcodes for Astragalus sp. cemA and rpl33 were found undergoing positive selection by the section patterns in the coded protein. Phylogenetic analysis showed that Astragalus is a monophyletic group closely related to the genus Oxytropis within the tribe Galegeae. The newly sequenced chloroplast genomes provide insight into the unresolved evolutionary relationships within Astragalus spp. and are expected to contribute to species identification.

Characterization and phylogenetic analysis of the complete chloroplast genome of Juncus effusus L

Juncus effusus L., a perennial herbaceous species of family Juncaceae, is distributed mainly in warm areas worldwide. We studied the complete chloroplast (cp) genome of J. effusus through next-generation sequencing technology. The whole cp genome contained 170,612 base pairs, with the GC ratio as 35.99%. The 70 genes annotated from the cp genome include 32 protein coding genes, eight rRNA genes and 30 tRNA genes. The genome’s large single-copy region (LSC) was 80,640 bp, with the small single-copy region (SSC) 64,718 bp, and inverted repeat (IR) 12,627 bp. Furthermore, a phylogenetic tree was generated to evaluate evolutionary relationship between J. effusus and relevant species. This study will be beneficial for the further understanding and application of J. effusus.

Comparative chloroplast genomes and phylogenetic analysis of Aquilegia

PREMISE

Aquilegia is an ideal taxon for studying the evolution of adaptive radiation. Current phylogenies of Aquilegia based on different molecular markers are inconsistent, and therefore a clear and accurate phylogeny remains uncertain. Analyzing the chloroplast genome, with its simple structure and low recombination rate, may help solve this problem.

METHODS

Next-generation sequencing data were generated or downloaded for Aquilegia species, enabling their chloroplast genomes to be assembled. The assemblies were used to estimate the genome characteristics and infer the phylogeny of Aquilegia.

RESULTS

In this study, chloroplast genome sequences were assembled for Aquilegia species distributed across Asia, North America, and Europe. Three of the genes analyzed (petG, rpl36, and atpB) were shown to be under positive selection and may be related to adaptation. The phylogenetic tree of Aquilegia showed that its member species formed two clades with high support, North American and European species, with the Asian species being paraphyletic; A. parviflora and A. amurensis clustered with the North American species, while the remaining Asian species were found in the European clade. In addition, A. oxysepala var. kansuensis should be considered as a separate species rather than a variety.

DISCUSSION

The complete chloroplast genomes of these Aquilegia species provide new insights into the reconstruction of the phylogeny of related species and contribute to the further study of this genus.

Comparative Chloroplast Genomics of Sophora Species: Evolution and Phylogenetic Relationships in the Early-Diverging Legume Subfamily Papilionoideae (Fabaceae)

The taxonomy and evolutionary history of Sophora L., a genus with high economic and medicinal value, remain uncertain due to the absence of genetic resource (especially in China) and low polymorphism of molecular markers. Our aim was to elucidate the molecular evolution and phylogenetic relationships in chloroplast genomes of Sophora species in the early-diverging legume subfamily Papilionoideae (Fabaceae). We reported nine Sophora chloroplast genome from China using Illumina sequencing. We performed a series of analyses with previously published genomes of Sophora species to investigate their genomic characteristics, identified simple sequence repeats, large repeat sequences, tandem repeats, and highly polymorphic loci. The genomes were 152,953-158,087 bp in length, and contained 111-113 unique genes, including 76-78 protein coding, 31 tRNA, and 4 rRNA. The expansion of inverted repeat boundary of Sophora resulted in rps12 entering into the LSC region and loss of trnT-CGU gene in some species. Also, we found an approximately 23 kb inversion between trnC-GCA and trnF-GAA within the genus. In addition, we identified seven highly polymorphic loci (pi (π) > 0.035) suitable for inferring the phylogeny of Sophora species. Among these, three regions also co-occurred with large repeat sequences and support use of repeats as a proxy for the identification of polymorphic loci. Based on whole chloroplast genome and protein-coding sequences data-set, a well-supported phylogenetic tree of Sophora and related taxa showed that this genus is monophyletic, but sect. Disamaea and sect. Sophora, are incongruent with traditional taxonomic classifications based on fruit morphology. Our finding provides significant genetic resources to support further investigation into the phylogenetic relationship and evolution of the genus Sophora.