Complete chloroplast genome of Sophora alopecuroides (Papilionoideae): molecular structures, comparative genome analysis and phylogenetic analysis

The complete chloroplast genome of Oxytropis ochrocephala Bunge 1874 (Fabaceae) and its phylogenetic analysis

This study presents the first-ever complete chloroplast (cp) genome sequence of Oxytropis ochrocephala Bunge 1874, a member of the Fabaceae family. The cp genome spans 126,996 base pairs and includes 109 genes, comprising 76 protein-coding genes, 29 tRNA genes, and four rRNA genes. Notably, the genome lacks an inverted repeat (IR) region. Additionally, we constructed phylogenetic trees for 34 species within Trib. Galegeae, employing both maximum likelihood (ML) and Bayesian inference (BI) methods. These analyses robustly support the monophyly of the Oxytropis species, evidenced by high bootstrap values (BP = 100) and posterior probabilities (PP = 1). Within this clade, O. ochrocephala exhibits a sister relationship with other Oxytropis species. Our findings provide valuable insights into the genetic makeup and evolutionary relationships of O. ochrocephala within the Galegeae tribe.

Characteristics of chloroplast and mitochondrial genomes and intracellular gene transfer in the Korean endemic shrub, Sophora koreensis Nakai (Fabaceae)

Sophora koreensis Nakai, an endemic species distributed only in the Korean Peninsula, is of great geographical, economic, and taxonomic importance. Although its complete chloroplast (cp) genome sequence has been reported, its mitochondrial (mt) genome sequence has not yet been studied. Therefore, in this study, we aimed to investigate its mt genome sequence and compare it with those reported for other Fabaceae species. Total genomic DNA was extracted from fresh S. koreensis leaves collected from natural habitats in Gangwon-do Province, South Korea. This was followed by polymerase chain reaction (PCR) amplification of cpDNA insertions in the mt genome and the detection of microsatellites and dispersed repeats in the cp and mt genomes. Finally, the cp and mt genomes of S. koreensis were compared with those reported for other Fabaceae species. The cp sequence of S. koreensis showed identical gene orders and contents as those previously reported. Only six substitutions and one deletion were detected with 99 % homology. Conversely, the complete mt genome sequence, which was 517,845 bp in length and encoded 61 genes, including 43 protein-coding, 15 transfer RNAs, and 3 ribosomal RNA genes, was considerably different from that of S. japonica in terms of gene order and composition. Further, the mt genome of S. koreensis included ca. 7 and 3 kb insertions, representing an intracellular gene transfer (IGT) event, and the regions with these insertions were determined to be originally present in the cp genome. This IGT event was also confirmed via PCR amplification. IGT events can be induced via biological gene expression control or the use of repetitive sequences, and they provide important insights into the evolutionary lineage of S. koreensis. However, further studies are needed to clarify the gene transfer mechanisms between the two organelles.

Sequencing and Phylogenetic Analysis of the Chloroplast Genome of Three Apricot Species

The production and quality of apricots in China is currently limited by the availability of germplasm resource characterizations, including identification at the species and cultivar level. To help address this issue, the complete chloroplast genomes of Prunus armeniaca L., P. sibirica L. and kernel consumption apricot were sequenced, characterized, and phylogenetically analyzed. The three chloroplast (cp) genomes ranged from 157,951 to 158,224 bp, and 131 genes were identified, including 86 protein-coding genes, 37 rRNAs, and 8 tRNAs. The GC content ranged from 36.70% to 36.75%. Of the 170 repetitive sequences detected, 42 were shared by all three species, and 53-57 simple sequence repeats were detected with AT base preferences. Comparative genomic analysis revealed high similarity in overall structure and gene content as well as seven variation hotspot regions, including psbA-trnK-UUU, rpoC1-rpoB, rpl32-trnL-UAG, trnK-rps16, ndhG-ndhI, ccsA-ndhD, and ndhF-trnL. Phylogenetic analysis showed that the three apricot species clustered into one group, and the genetic relationship between P. armeniaca and kernel consumption apricot was the closest. The results of this study provide a theoretical basis for further research on the genetic diversity of apricots and the development and utilization of molecular markers for the genetic engineering and breeding of apricots.

Complete chloroplast genome of Lens lamottei reveals intraspecies variation among with Lens culinaris

Lens lamottei is a member of the Fabaceae family and the second gene pool of the genus Lens. The environmental factors that drove the divergence among wild and cultivated species have been studied extensively. Recent research has focused on genomic signatures associated with various phenotypes with the acceleration of next-generation techniques in molecular profiling. Therefore, in this study, we provide the complete sequence of the chloroplast genome sequence in the wild Lens species L. lamottei with a deep coverage of 713 × next-generation sequencing (NGS) data for the first time. Compared to the cultivated species, Lens culinaris, we identified synonymous, and nonsynonymous changes in the protein-coding regions of the genes ndhB, ndhF, ndhH, petA, rpoA, rpoC2, rps3, and ycf2 in L. lamottei. Phylogenetic analysis of chloroplast genomes of various plants under Leguminosae revealed that L. lamottei and L. culinaris are closest to one another than to other species. The complete chloroplast genome of L. lamottei also allowed us to reanalyze previously published transcriptomic data, which showed high levels of gene expression for ATP-synthase, rubisco, and photosystem genes. Overall, this study provides a deeper insight into the diversity of Lens species and the agricultural importance of these plants through their chloroplast genomes.

Analysis of Complete Chloroplast Genome: Structure, Phylogenetic Relationships of Galega orientalis and Evolutionary Inference of Galegeae

Galega orientalis, a leguminous herb in the Fabaceae family, is an ecologically and economically important species widely cultivated for its strong stress resistance and high protein content. However, genomic information of Galega orientalis has not been reported, which limiting its evolutionary analysis. The small genome size makes chloroplast relatively easy to obtain genomic sequence for phylogenetic studies and molecular marker development. Here, the chloroplast genome of Galega orientalis was sequenced and annotated. The results showed that the chloroplast genome of G. orientalis is 125,280 bp in length with GC content of 34.11%. A total of 107 genes were identified, including 74 protein-coding genes, 29 tRNAs and four rRNAs. One inverted repeat (IR) region was lost in the chloroplast genome of G. orientalis. In addition, five genes (rpl22, ycf2, rps16, trnE-UUC and pbf1) were lost compared with the chloroplast genome of its related species G. officinalis. A total of 84 long repeats and 68 simple sequence repeats were detected, which could be used as potential markers in the genetic studies of G. orientalis and related species. We found that the Ka/Ks values of three genes petL, rpl20, and ycf4 were higher than one in the pairwise comparation of G. officinalis and other three Galegeae species (Calophaca sinica, Caragana jubata, Caragana korshinskii), which indicated those three genes were under positive selection. A comparative genomic analysis of 15 Galegeae species showed that most conserved non-coding sequence regions and two genic regions (ycf1 and clpP) were highly divergent, which could be used as DNA barcodes for rapid and accurate species identification. Phylogenetic trees constructed based on the ycf1 and clpP genes confirmed the evolutionary relationships among Galegeae species. In addition, among the 15 Galegeae species analyzed, Galega orientalis had a unique 30-bp intron in the ycf1 gene and Tibetia liangshanensis lacked two introns in the clpP gene, which is contrary to existing conclusion that only Glycyrrhiza species in the IR lacking clade (IRLC) lack two introns. In conclusion, for the first time, the complete chloroplast genome of G. orientalis was determined and annotated, which could provide insights into the unsolved evolutionary relationships within the genus Galegeae.

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

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

The complete chloroplast genome of Flemingia macrophylla (Willd.) Prain (Fabaceae) from Guangxi, China

Flemingia macrophylla (Willd.) Prain is an ethnomedicinal plant with high nutritional and medicinal values. In this study, we report the complete chloroplast genome of F. macrophylla. The chloroplast genome has a typical quadripartite structure with a genome size of 152,988 bp, including a large single-copy (LSC) of 83,634 bp, a small single-copy (SSC) of 17,774 bp and two inverted repeats (IRs) of 25,790 bp. The genome contains 129 genes, including 84 protein-coding, 37 tRNA and 8 rRNA genes. The overall GC content is 35.1%. Phylogenetic analysis showed that F. macrophylla grouped with a clade containing the genera of Fagelia, Dolichos, Eriosema, Dunbaria and Cajanus in Fabaceae. This study provides essential data and insight for understanding the phylogenetic placement of Flemingia.

Screening and identification of salt-tolerance genes in Sophora alopecuroides and functional verification of SaAQP

Overexpression of SaAQP can improve the salt tolerance of transgenic soybean hairy roots and A. thaliana. Salt stress severely affects crop yield and food security. There is a need to improve the salt tolerance of crops, but the discovery and utilization of salt-tolerance genes remains limited. Owing to its strong stress tolerance, Sophora alopecuroides is ideal for the identification of salt-tolerance genes. Therefore, we aimed to screen and identify the salt-tolerance genes in S. alopecuroides. With a yeast expression library of seedlings, salt-tolerant genes were screened using a salt-containing medium to simulate salt stress. By combining salt-treatment screening and transcriptome sequencing, 11 candidate genes related to salt tolerance were identified, including genes for peroxidase, inositol methyltransferase, aquaporin, cysteine synthase, pectinesterase, and WRKY. The expression dynamics of candidate genes were analyzed after salt treatment of S. alopecuroides, and salt tolerance was verified in yeast BY4743. The candidate genes participated in the salt-stress response in S. alopecuroides, and their overexpression significantly improved the salt tolerance of yeast. Salt tolerance mediated by SaAQP was further verified in soybean hairy roots and Arabidopsis thaliana, and it was found that SaAQP might enhance the salt tolerance of A. thaliana by participating in a reactive oxygen species scavenging mechanism. This result provides new genetic resources in plant breeding for salt resistance.

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.

The complete chloroplast genome of Euphorbia hirta (Euphorbiaceae), a commonly used medicinal plant in China

Plants in the genus Euphorbia have been widely used as herbal medicine, and for ornamental horticulture and biofuel production. In this study, we report the complete chloroplast genome of Euphorbia hirta which is known as the ‘asthma-plant’ due to its medicinal use. The chloroplast genome of this species is 164,340 bp in length, including a pair of inverted repeat regions (IRs) (27,354 bp) that are divided by a large single-copy region (LSC) (91,373 bp) and a small single-copy region (SSC) (18,259 bp). The chloroplast genome of E. hirta contains 111 unique genes (77 protein-coding, 30 tRNA, and four rRNA), 19 of which are duplicated in the IR regions. The overall GC content is 35.4%. Phylogenetic analysis fully resolved E. hirta groups with other species of Euphorbia. The complete chloroplast genome of E. hirta provides useful information that can be used to distinguish related species and reconstruct evolutionary relationships.

Chloroplast genome sequence of Chongming lima bean (Phaseolus lunatus L.) and comparative analyses with other legume chloroplast genomes

Background

Lima bean (Phaseolus lunatus L.) is a member of subfamily Phaseolinae belonging to the family Leguminosae and an important source of plant proteins for the human diet. As we all know, lima beans have important economic value and great diversity. However, our knowledge of the chloroplast genome level of lima beans is limited.

Results

The chloroplast genome of lima bean was obtained by Illumina sequencing technology for the first time. The Cp genome with a length of 150,902 bp, including a pair of inverted repeats (IRA and IRB 26543 bp each), a large single-copy (LSC 80218 bp) and a small single-copy region (SSC 17598 bp). In total, 124 unique genes including 82 protein-coding genes, 34 tRNA genes, and 8 rRNA genes were identified in the P. lunatus Cp genome. A total of 61 long repeats and 290 SSRs were detected in the lima bean Cp genome. It has a typical 50 kb inversion of the Leguminosae family and an 70 kb inversion to subtribe Phaseolinae. rpl16, accD, petB, rsp16, clpP, ndhA, ndhF and ycf1 genes in coding regions was found significant variation, the intergenic regions of trnk-rbcL, rbcL-atpB, ndhJ-rps4, psbD-rpoB, atpI-atpA, atpA-accD, accD-psbJ, psbE-psbB, rsp11-rsp19, ndhF-ccsA was found in a high degree of divergence. A phylogenetic analysis showed that P. lunatus appears to be more closely related to P. vulgaris, V.unguiculata and V. radiata.

Conclusions

The characteristics of the lima bean Cp genome was identified for the first time, these results will provide useful insights for species identification, evolutionary studies and molecular biology research.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07467-8.

Chemosystematic and evolutionary trends of the genistoid clade sensu stricto (Papilionoideae, Fabaceae)

The Papilionoideae, which comprises 503 genera and approximately 14,000 species, is the largest and most diverse subfamily of the Fabaceae family. In this subfamily, the Crotalarieae, Genisteae, Podalyrieae, Thermopsideae, Sophoreae and Euchresteae tribes are closely related by micro and macromolecular features, thus forming the genistoid clade. This group combines well-known genera, whereas other genera lack phytochemical and chemotaxonomic studies. Thus, this work aimed to characterize the special metabolites in these genera in order to define the chemical profile, the micromolecular markers and the chemical diversity, as well as to evaluate the group evolutionary trends. Flavonoids and alkaloids were identified as chemosystematic markers for the studied tribes due to high occurrence number and structural diversity. Among flavonoids, the flavones and isoflavones predominated. Low protection indexes of flavonoid hydroxyls by O-glycosylation or O-methylation were observed, whereas C-prenylation and C-glycosylation were frequent, mainly at C-6 and C-8 positions. The flavone/flavonol ratio shows the predominance of the flavones. Quinolizidine and piperidine alkaloids were present in most genera. Pyrrolizidine alkaloids were found in a few genera from Thermopsideae, Genisteae and Crotalarieae, which suggests a mechanism of adaptive convergence. Cluster analysis allowed separation of genera for each tribe by chemical similarities. The micromolecular trends of protection of flavonoid hydroxyls and alkaloid oxidation indicate the genistoid clade is through evolutionary transition, which is consistent with its phylogenetic position in the Papilionoideae subfamily.

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.

The complete chloroplast genome sequence of the medicinal plant Sophora tonkinensis

Sophora tonkinensis belongs to genus Sophora of the Fabaceae family. It is mainly distributed in the ridge and peak regions of limestone areas in western China and has high medicinal value and important ecological functions. Wild populations of S. tonkinensis are in danger and need urgent conservation. Furthermore, wild S. tonkinensis resources are very limited relative to the needs of the market, and many adulterants are present on the market. Therefore, a method for authenticating S. tonkinensis and its adulterants at the molecular level is needed. Chloroplast genomes are valuable sources of genetic markers for phylogenetic analyses, genetic diversity evaluation, and plant molecular identification. In this study, we report the complete chloroplast genome of S. tonkinensis. The circular complete chloroplast genome was 154,644 bp in length, containing an 85,810 bp long single-copy (LSC) region, an 18,321 bp short single-copy (SSC) region and two inverted repeat (IR) regions of 50,513 bp. The S. tonkinensis chloroplast genome comprised 129 genes, including 83 protein-coding genes, 38 transfer RNA (tRNA) genes, and 8 ribosomal RNA (rRNA) genes. The structure, gene order and guanine and cytosine (GC) content of the S. tonkinensis chloroplast genome were similar to those of the Sophora alopecuroides and Sophora flavescens chloroplast genomes. A total of 1,760 simple sequence repeats (SSRs) were identified in the chloroplast genome of S. tonkinensis, and most of them (93.1%) were mononucleotides. Moreover, the identified SSRs were mainly distributed in the LSC region, accounting for 60% of the total number of SSRs, while 316 (18%) and 383 (22%) were located in the SSC and IR regions, respectively. Only one complete copy of the rpl2 gene was present at the LSC/IRB boundary, while another copy was absent from the IRA region because of the incomplete structure caused by IR region expansion and contraction. The phylogenetic analysis placed S. tonkinensis in Papilionoideae, sister to S. flavescens, and the genera Sophora and Ammopiptanthus were closely related. The complete genome sequencing and chloroplast genome comparative analysis of S. tonkinensis and its closely related species presented in this paper will help formulate effective conservation and management strategies as well as molecular identification approaches for this important medicinal plant.