Plastid genome evolution in tribe Desmodieae (Fabaceae: Papilionoideae)

Campylotropis xinfeniae (Fabaceae, Papilionoideae), a new species from Yunnan, China, based on morphological and molecular evidence

Campylotropis xinfeniae, a new species from the dry-hot valley of the Jinsha River in the Yunnan province, China, is described and illustrated. It is morphologically similar to C. wilsonii and C. brevifolia in having glabrescent old branches, absent stipels, 3-foliolate leaves, and adaxially puberulent leaflets, while it differs from the latter two in having often paniculate inflorescences, obviously white standard, not incurved sickle keel, larger narrowly oblique legumes, and longer legume beak. The complete chloroplast genome of this new species is 149,073 bp in length and exhibits a typical quadripartite structure. Phylogenetic analyses based on the complete chloroplast genome also supported C. xinfeniae as a new species located at the basal distinct clade of the genus Campylotropis, clearly separated from the remaining members of the genus and its allied genera. A conservation assessment of data deficient (DD) is recommended for the new species without extensive exploring of similar habitats according to the IUCN Categories and Criteria.

Chloroplast genome sequence of Nothapodytes nimmoniana (Icacinaceae), a vulnerable reservoir of camptothecin from Western Ghats

Nothapodytes nimmoniana belongs to family Icacinaceae and is a major source of compound Camptothecin. The global demand for Camptothecin has caused large-scale exploitation of N. nimmoniana from its wild habitat in Western Ghats of India, thereby making it vulnerable. The species is known to exhibit genetic diversity among the populations in Western Ghats. In this study, we report plastome sequence of N. nimmoniana, first for the genus. For the study, the species was collected from Western Ghats of Maharashtra. The plastome of N. nimmoniana was 150,726 bp in length and exhibited typical quadripartite structure with 83,771 bp LSC, 18,513 bp SSC and 24,221 IR region. The plastome was characterized by presence of 124 unique genes, 87 protein coding genes, 29 tRNA genes and four rRNA genes. Further, the plastome was compared with the available basal lamiid plastomes for gene order and composition. N. nimmoniana plastome exhibited SSC region in an inverted configuration. Phylogenomic study placed N. nimmoniana sister to Mappia mexicana. The SSR markers identified in this study, might help to distinguish genetically diverse populations, prioritizing the populations which need immediate conservation effects as well as for checking adulteration.

Comparative Analysis of Chloroplast Genome of Desmodium stryacifolium with Closely Related Legume Genome from the Phaseoloid Clade

Desmodium styracifolium is a medicinal plant from the Desmodieae tribe, also known as Grona styracifolia. Its role in the treatment of urolithiasis, urinary infections, and cholelithiasis has previously been widely documented. The complete chloroplast genome sequence of D. Styracifolium is 149,155 bp in length with a GC content of 35.2%. It is composed of a large single copy (LSC) of 82,476 bp and a small single copy (SSC) of 18,439 bp, which are separated by a pair of inverted repeats (IR) of 24,120 bp each and has 128 genes. We performed a comparative analysis of the D. styracifolium cpDNA with the genome of previously investigated members of the Sesamoidea tribe and on the outgroup from its Phaseolinae sister tribe. The size of all seven cpDNAs ranged from 148,814 bp to 151,217 bp in length due to the contraction and expansion of the IR/SC boundaries. The gene orientation of the SSC region in D. styracifolium was inverted in comparison with the other six studied species. Furthermore, the sequence divergence of the IR regions was significantly lower than that of the LSC and the SSC, and five highly divergent regions, trnL-UAA-trnT-UGU, psaJ-ycf4, psbE-petL, rpl36-rps8, and rpl32-trnL-UGA, were identified that could be used as valuable molecular markers in future taxonomic studies and phylogenetic constructions.

Comparative analysis of the chloroplast genomes of eight Piper species and insights into the utilization of structural variation in phylogenetic analysis

With more than 2000 species, Piper is regarded as having high medicinal, cosmetic, and edible value. There also remain some taxonomic and evolutionary uncertainties about the genus. This study performed chloroplast genome sequencing of eight poorly studied Piper species and a comparative analysis with black pepper (Piper nigrum). All examined species were highly similar in gene content, with 79 protein-coding genes, 24 tRNAs, and four rRNAs. They also harbored significant structural differences: The number of SSRs ranged from 63 to 87, over 10,000 SNPs were detected, and over 1,000 indels were found. The spatial distribution of structural differences was uneven, with the IR and LSC being relatively more conserved and the SSC region highly variable. Such structural variations of the chloroplast genome can help in evaluating the phylogenetic relationships between species, deciding some hard-to-distinguish evolutionary relationships, or eliminating improper markers. The SSC region may be evolving at high speed, and some species showed a high degree of sequence variation in the SSC region, which seriously affected marker sequence detection. Conversely, CDS sequences tended to lack variation, and some CDSs can serve as ideal markers for phylogenetic reconstruction. All told, this study provides an effective strategy for selecting chloroplast markers, analyzing difficult-to-distinguish phylogenetic relationships and avoiding the taxonomic errors caused by high degree of sequence variations.

Complete Chloroplast Genomes Provide Insights Into Evolution and Phylogeny of Campylotropis (Fabaceae)

The genus Campylotropis Bunge (Desmodieae, Papilionoideae) comprises about 37 species distributed in temperate and tropical Asia. Despite the great potential in soil conservation, horticulture, and medicine usage, little is known about the evolutionary history and phylogenetic relationships of Campylotropis due to insufficient genetic resources. Here, we sequenced and assembled 21 complete chloroplast genomes of Campylotropis species. In combination with the previously published chloroplast genomes of C. macrocarpa and closely related species, we conducted comparative genomics and phylogenomic analysis on these data. Comparative analysis of the genome size, structure, expansion and contraction of inverted repeat (IR) boundaries, number of genes, GC content, and pattern of simple sequence repeats (SSRs) revealed high similarities among the Campylotropis chloroplast genomes. The activities of long sequence repeats contributed to the variation in genome size and gene content in Campylotropis chloroplast genomes. The Campylotropis chloroplast genomes showed moderate sequence variation, and 13 highly variable regions were identified for species identification and further phylogenetic studies. We also reported one more case of matK pseudogene in the legume family. The phylogenetic analysis confirmed the monophyly of Campylotropis and the sister relationship between Lespedeza and Kummerowia, the latter two genera were then sister to Campylotropis. The intrageneric relationships of Campylotropis based on genomic scale data were firstly reported in this study. The two positively selected genes (atpF and rps19) and eight fast-evolving genes identified in this study may help us to understand the adaptation of Campylotropis species. Overall, this study enhances our understanding of the chloroplast genome evolution and phylogenetic relationships of Campylotropis.

Highly Resolved Papilionoid Legume Phylogeny Based on Plastid Phylogenomics

Comprising 501 genera and around 14,000 species, Papilionoideae is not only the largest subfamily of Fabaceae (Leguminosae; legumes), but also one of the most extraordinarily diverse clades among angiosperms. Papilionoids are a major source of food and forage, are ecologically successful in all major biomes, and display dramatic variation in both floral architecture and plastid genome (plastome) structure. Plastid DNA-based phylogenetic analyses have greatly improved our understanding of relationships among the major groups of Papilionoideae, yet the backbone of the subfamily phylogeny remains unresolved. In this study, we sequenced and assembled 39 new plastomes that are covering key genera representing the morphological diversity in the subfamily. From 244 total taxa, we produced eight datasets for maximum likelihood (ML) analyses based on entire plastomes and/or concatenated sequences of 77 protein-coding sequences (CDS) and two datasets for multispecies coalescent (MSC) analyses based on individual gene trees. We additionally produced a combined nucleotide dataset comprising CDS plus matK gene sequences only, in which most papilionoid genera were sampled. A ML tree based on the entire plastome maximally supported all of the deep and most recent divergences of papilionoids (223 out of 236 nodes). The Swartzieae, ADA (Angylocalyceae, Dipterygeae, and Amburaneae), Cladrastis, Andira, and Exostyleae clades formed a grade to the remainder of the Papilionoideae, concordant with nine ML and two MSC trees. Phylogenetic relationships among the remaining five papilionoid lineages (Vataireoid, Dermatophyllum, Genistoid s.l., Dalbergioid s.l., and Baphieae + Non-Protein Amino Acid Accumulating or NPAAA clade) remained uncertain, because of insufficient support and/or conflicting relationships among trees. Our study fully resolved most of the deep nodes of Papilionoideae, however, some relationships require further exploration. More genome-scale data and rigorous analyses are needed to disentangle phylogenetic relationships among the five remaining lineages.

Identification of evolutionary relationships and DNA markers in the medicinally important genus Fritillaria based on chloroplast genomics

The genus Fritillaria has attracted great attention because of its medicinal and ornamental values. At least three reasons, including the accurate discrimination between various Fritillaria species, protection and sustainable development of rare Fritillaria resources as well as understanding of relationship of some perplexing species, have prompted phylogenetic analyses and development of molecular markers for Fritillaria species. Here we determined the complete chloroplast (CP) genomes for F. unibracteata, F. przewalskii, F. delavayi, and F. sinica through Illumina sequencing, followed by de novo assembly. The lengths of the genomes ranged from 151,076 in F. unibracteata to 152,043 in F. przewalskii. Those CP genomes displayed a typical quadripartite structure, all including a pair of inverted repeats (26,078 to 26,355 bp) separated by the large single-copy (81,383 to 81,804 bp) and small single-copy (17,537 to 17,569 bp) regions. Fritillaria przewalskii, F. delavayi, and F. sinica equivalently encoded 133 unique genes consisting of 38 transfer RNA genes, eight ribosomal RNA genes, and 87 protein coding genes, whereas F. unibracteata contained 132 unique genes due to absence of the rps16 gene. Subsequently, comparative analysis of the complete CP genomes revealed that ycf1, trnL, trnF, ndhD, trnN-trnR, trnE-trnT, trnN, psbM-trnD, atpI, and rps19 to be useful molecular markers in taxonomic studies owning to their interspecies variations. Based on the comprehensive CP genome data collected from 53 species in Fritillaria and Lilium genera, a phylogenomic study was carried out with three Cardiocrinum species and five Amana species as outgroups. The results of the phylogenetic analysis showed that Fritillaria was a sister to Lilium, and the interspecies relationships within subgenus Fritillaria were well resolved. Furthermore, phylogenetic analysis based on the CP genome was proved to be a promising method in selecting potential novel medicinal resources to substitute current medicinal species that are on the verge of extinction.

Complete Plastome of Three Korean Asarum (Aristolochiaceae): Confirmation Tripartite Structure within Korean Asarum and Comparative Analyses

The genus Asarum (Aristolochiaceae) is a well-known resource of medicinal and ornamental plants. However, the taxonomy of Korean Asarum is ambiguous due to their considerable morphological variations. Previously, a unique plastome structure has been reported from this genus. Therefore, we investigated the structural change in the plastomes within three Korean Asarum species and inferred their phylogenetic relationships. The plastome sizes of Asarum species assembled here range from 190,168 to 193,356 bp, which are longer than a typical plastome size (160 kb). This is due to the incorporation and duplication of the small single copy into the inverted repeat, which resulted in a unique tripartite structure. We first verified this unique structure using the Illumina Miseq and Oxford Nanopore MinION platforms. We also investigated the phylogeny of 26 Aristolochiaceae species based on 79 plastid protein-coding genes, which supports the monophyly of Korean Asarum species. Although the 79 plastid protein-coding gene data set showed some limitations in supporting the previous classification, it exhibits its effectiveness in delineating some sections and species. Thus, it can serve as an effective tool for resolving species-level phylogeny in Aristolochiaceae. Last, we evaluated variable sites and simple sequence repeats in the plastome as potential molecular markers for species delimitation.

Comparative Mitogenomic Analysis Reveals Gene and Intron Dynamics in Rubiaceae and Intra-Specific Diversification in Damnacanthus indicus

The dynamic evolution of mitochondrial gene and intron content has been reported across the angiosperms. However, a reference mitochondrial genome (mitogenome) is not available in Rubiaceae. The phylogenetic utility of mitogenome data at a species level is rarely assessed. Here, we assembled mitogenomes of six Damnacanthus indicus (Rubiaceae, Rubioideae) representing two varieties (var. indicus and var. microphyllus). The gene and intron content of D. indicus was compared with mitogenomes from representative angiosperm species and mitochondrial contigs from the other Rubiaceae species. Mitogenome structural rearrangement and sequence divergence in D. indicus were analyzed in six individuals. The size of the mitogenome in D. indicus varied from 417,661 to 419,435 bp. Comparing the number of intact mitochondrial protein-coding genes in other Gentianales taxa (38), D. indicus included 32 genes representing several losses. The intron analysis revealed a shift from cis to trans splicing of a nad1 intron (nad1i728) in D. indicus and it is a shared character with the other four Rubioideae taxa. Two distinct mitogenome structures (type A and B) were identified. Two-step direct repeat-mediated recombination was proposed to explain structural changes between type A and B mitogenomes. The five individuals from two varieties in D. indicus diverged well in the whole mitogenome-level comparison with one exception. Collectively, our study elucidated the mitogenome evolution in Rubiaceae along with D. indicus and showed the reliable phylogenetic utility of the whole mitogenome data at a species-level evolution.

Insights into phylogenetic relationships and genome evolution of subfamily Commelinoideae (Commelinaceae Mirb.) inferred from complete chloroplast genomes

Background

Commelinaceae (Commelinales) comprise 41 genera and are widely distributed in both the Old and New Worlds, except in Europe. The relationships among genera in this family have been suggested in several morphological and molecular studies. However, it is difficult to explain their relationships due to high morphological variations and low support values. Currently, many researchers have been using complete chloroplast genome data for inferring the evolution of land plants. In this study, we completed 15 new plastid genome sequences of subfamily Commelinoideae using the Mi-seq platform. We utilized genome data to reveal the structural variations and reconstruct the problematic positions of genera for the first time.

Results

All examined species of Commelinoideae have three pseudogenes (accD, rpoA, and ycf15), and the former two might be a synapomorphy within Commelinales. Only four species in tribe Commelineae presented IR expansion, which affected duplication of the rpl22 gene. We identified inversions that range from approximately 3 to 15 kb in four taxa (Amischotolype, Belosynapsis, Murdannia, and Streptolirion). The phylogenetic analysis using 77 chloroplast protein-coding genes with maximum parsimony, maximum likelihood, and Bayesian inference suggests that Palisota is most closely related to tribe Commelineae, supported by high support values. This result differs significantly from the current classification of Commelinaceae. Also, we resolved the unclear position of Streptoliriinae and the monophyly of Dichorisandrinae. Among the ten CDS (ndhH, rpoC2, ndhA, rps3, ndhG, ndhD, ccsA, ndhF, matK, and ycf1), which have high nucleotide diversity values (Pi > 0.045) and over 500 bp length, four CDS (ndhH, rpoC2, matK, and ycf1) show that they are congruent with the topology derived from 77 chloroplast protein-coding genes.

Conclusions

In this study, we provide detailed information on the 15 complete plastid genomes of Commelinoideae taxa. We identified characteristic pseudogenes and nucleotide diversity, which can be used to infer the family evolutionary history. Also, further research is needed to revise the position of Palisota in the current classification of Commelinaceae.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07541-1.

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.

Comparative Genomics and Phylogenetic Analyses of Christia vespertilionis and Urariopsis brevissima in the Tribe Desmodieae (Fabaceae: Papilionoideae) Based on Complete Chloroplast Genomes

Taxonomic and phylogenetic relationships of Christia, Urariopsis, Uraria and related genera within the tribe Desmodieae (Fabaceae: Papilionoideae) have long been controversial. Here, we report the complete chloroplast (cp) genomes of Christia vespertilionis and Urariopsis brevissima and perform comparative and phylogenetic analyses with Uraria lagopodioides and other relatives in the Desmodieae. The cp genomes of C. vespertilionis and U. brevissima are 149,656 and 149,930 bp long, with 128 unique genes (83 protein-coding genes, 37 tRNA genes and 8 rRNA genes), respectively. Comparative analyses revealed 95-129 simple sequence repeats (SSRs) and eleven highly variable regions (trnK-rbcL, rbcL-atpB, ndhJ-trnF, trnL-trnT, psbD-rpoB, accD-cemA, petA-psbL, psbE-petL, rps11-rps19, ndhF-ccsA, and rps15-ycf1) among six Desmodieae species. Phylogenetic analyses clearly resolved two subtribes (Desmodiinae and Lespedezinae) of Desmodieae as monophyletic, and the newly reported C. vespertilionis and U. brevissima clustered in subtribe Desmodiinae. A sister relationship of C. vespertilionis to U. lagopodioides was supported. Evidence was presented to support the treatment of Urariopsis as a distinct genus rather than in synonymy with Uraria. The results provide valuable information for further studies on species delimitation, phylogenetics, population genetics, and the evolutionary process of speciation in the Desmodieae.

Insights Into Chloroplast Genome Evolution Across Opuntioideae (Cactaceae) Reveals Robust Yet Sometimes Conflicting Phylogenetic Topologies

Chloroplast genomes (plastomes) are frequently treated as highly conserved among land plants. However, many lineages of vascular plants have experienced extensive structural rearrangements, including inversions and modifications to the size and content of genes. Cacti are one of these lineages, containing the smallest plastome known for an obligately photosynthetic angiosperm, including the loss of one copy of the inverted repeat (∼25 kb) and the ndh gene suite, but only a few cacti from the subfamily Cactoideae have been sufficiently characterized. Here, we investigated the variation of plastome sequences across the second-major lineage of the Cactaceae, the subfamily Opuntioideae, to address (1) how variable is the content and arrangement of chloroplast genome sequences across the subfamily, and (2) how phylogenetically informative are the plastome sequences for resolving major relationships among the clades of Opuntioideae. Our de novo assembly of the Opuntia quimilo plastome recovered an organelle of 150,347 bp in length with both copies of the inverted repeat and the presence of all the ndh gene suite. An expansion of the large single copy unit and a reduction of the small single copy unit was observed, including translocations and inversion of genes, as well as the putative pseudogenization of some loci. Comparative analyses among all clades within Opuntioideae suggested that plastome structure and content vary across taxa of this subfamily, with putative independent losses of the ndh gene suite and pseudogenization of genes across disparate lineages, further demonstrating the dynamic nature of plastomes in Cactaceae. Our plastome dataset was robust in resolving three tribes with high support within Opuntioideae: Cylindropuntieae, Tephrocacteae and Opuntieae. However, conflicting topologies were recovered among major clades when exploring different assemblies of markers. A plastome-wide survey for highly informative phylogenetic markers revealed previously unused regions for future use in Sanger-based studies, presenting a valuable dataset with primers designed for continued evolutionary studies across Cactaceae. These results bring new insights into the evolution of plastomes in cacti, suggesting that further analyses should be carried out to address how ecological drivers, physiological constraints and morphological traits of cacti may be related with the common rearrangements in plastomes that have been reported across the family.

Extensive survey of the ycf4 plastid gene throughout the IRLC legumes: Robust evidence of its locus and lineage specific accelerated rate of evolution, pseudogenization and gene loss in the tribe Fabeae

The genome organization and gene content of plastome (plastid genome) are highly conserved among most flowering plant species. Plastome variation (in size and gene order) is rare in photosynthetic species but size variation, rearrangements and gene/intron losses is attributed to groups of seed plants. Fabaceae (legume family), in particular the subfamily Papilionoideae and the inverted repeat lacking clade (IRLC), a largest legume lineage, display the most dramatic and structural change which providing an excellent model for understanding of mechanisms of genomic evolution. The IRLC comprises 52 genera and ca 4000 species divided into seven tribes. In present study, we have sampled several representatives from each tribe across the IRLC from various herbaria and field. The ycf4 gene, which plays a role in regulating and assembly of photosystem I, is more variable in the tribe Fabeae than in other tribes. In certain species of Lathyrus, Pisum and Vavilovia, all belonging to Fabeae, the gene is either absent or a pseudogene. Our study suggests that ycf4 gene has undergone positive selection. Furthermore, the rapid evolution of the gene is locus and lineage specific and is not a shared character of the IRLC in legumes.

The Complete Chloroplast Genomes of Two Lespedeza Species: Insights into Codon Usage Bias, RNA Editing Sites, and Phylogenetic Relationships in Desmodieae (Fabaceae: Papilionoideae)

The genus Lespedeza (tribe: Desmodieae) consists of about 40 species that have high medicinal and economic value. However, in this genus, using morphological characters, the species identification is quite complicated, which can be solved by the analysis of the complete chloroplast genomes. As primary organelle genomes, the complete genome sequences of chloroplasts (cp) provide unique molecular information to study the divergence of species, RNA editing, and phylogeny. Therefore, to the best of our knowledge, for the first time, we sequenced the complete cp genomes of two representative Lespedeza species: Lespedeza davurica and Lespedeza cuneata. The cp genomes of both the species were found to be 149,010 bp in length, exhibiting the typical angiosperm chloroplast structure containing four regions. The Lespedeza cp genomes showed similar conserved gene contents, order, and orientations with a total GC content of 35.0%. A total of 128 genes, including 83 protein-coding genes, 37 tRNAs, and eight rRNAs, were identified from each genome. Unique molecular features of the two Lespedeza cp genome sequences were obtained by performing the analysis of repeats, sequence divergence, codon usage, and predicting the RNA editing sites in addition to phylogenetic analysis with other key genera in tribe Desmodieae. Using the two datasets, the phylogenetic relationship of Lespedeza species among Deasmodieae was discovered, suggesting that whole cp genomes provided useful information for phylogenetic studies of these species.

The complete chloroplast genome of Stryphnodendron adstringens (Leguminosae - Caesalpinioideae): comparative analysis with related Mimosoid species

Stryphnodendron adstringens is a medicinal plant belonging to the Leguminosae family, and it is commonly found in the southeastern savannas, endemic to the Cerrado biome. The goal of this study was to assemble and annotate the chloroplast genome of S. adstringens and to compare it with previously known genomes of the mimosoid clade within Leguminosae. The chloroplast genome was reconstructed using de novo and referenced-based assembly of paired-end reads generated by shotgun sequencing of total genomic DNA. The size of the S. adstringens chloroplast genome was 162,169 bp. This genome included a large single-copy (LSC) region of 91,045 bp, a small single-copy (SSC) region of 19,014 bp and a pair of inverted repeats (IRa and IRb) of 26,055 bp each. The S. adstringens chloroplast genome contains a total of 111 functional genes, including 77 protein-coding genes, 30 transfer RNA genes, and 4 ribosomal RNA genes. A total of 137 SSRs and 42 repeat structures were identified in S. adstringens chloroplast genome, with the highest proportion in the LSC region. A comparison of the S. adstringens chloroplast genome with those from other mimosoid species indicated that gene content and synteny are highly conserved in the clade. The phylogenetic reconstruction using 73 conserved coding-protein genes from 19 Leguminosae species was supported to be paraphyletic. Furthermore, the noncoding and coding regions with high nucleotide diversity may supply valuable markers for molecular evolutionary and phylogenetic studies at different taxonomic levels in this group.