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

Chloroplast genomes of Simarouba Aubl., molecular evolution and comparative analyses within Sapindales

Simarouba, a neotropical genus in the family Simaroubaceae, currently lacks comprehensive genomic data in existing databases. This study aims to fill this gap by providing genomic resources for three Simarouba species, S. amara, S. versicolor, and S. glauca. It also aims to perform comparative molecular evolutionary analyses in relation to other species within the order Sapindales. The analysis of these three Simarouba species revealed the presence of the typical quadripartite structure expected in plastomes. However, some pseudogenization events were identified in the psbC, infA, rpl22, and ycf1 genes. In particular, the CDS of the psbC gene in S. amara was reduced from 1422 bp to 584 bp due to a premature stop codon. Nucleotide diversity data pointed to gene and intergenic regions as promising candidates for species and family discrimination within the group, specifically matK, ycf1, ndhF, rpl32, petA-psbJ, and trnS-trnG. Selection signal analyses showed strong evidence for positive selection on the rpl23 gene. Phylogenetic analyses indicated that S. versicolor and S. glauca have a closer phylogenetic relationship than S. amara. We provide chloroplast genomes of three Simaruba species and use them to elucidate plastome evolution, highlight the presence of pseudogenization, and identify potential DNA barcode regions.

Structure, gene composition, divergence time and phylogeny analysis of the woody desert species Neltuma alba, Neltuma chilensis and Strombocarpa strombulifera

Neltuma alba (Algarrobo blanco), Neltuma chilensis (Algarrobo Chileno) and Strombocarpa strombulifera (Fortuna) are some of the few drought resistant trees and shrubs found in small highly fragmented populations, throughout the Atacama Desert. We reconstructed their plastid genomes using de novo assembly of paired-end reads from total genomic DNA. We found that the complete plastid genomes of N. alba and N. chilensis are larger in size compared to species of the Strombocarpa genus. The Strombocarpa species presented slightly more GC content than the Neltuma species. Therefore, we assume that Strombocarpa species have been exposed to stronger natural selection than Neltuma species. We observed high variation values in the number of cpSSRs (chloroplast simple sequence repeats) and repeated elements among Neltuma and Strombocarpa species. The p-distance results showed a low evolutionary divergence within the genus Neltuma, whereas a high evolutionary divergence was observed between Strombocarpa species. The molecular divergence time found in Neltuma and Strombocarpa show that these genera diverged in the late Oligocene. With this study we provide valuable information about tree species that provide important ecosystem services in hostile environments which can be used to determine these species in the geographically isolated communities, and keep the highly fragmented populations genetically healthy.

Comparative complete chloroplast genome of Geum japonicum: evolution and phylogenetic analysis

Geum japonicum (Rosaceae) has been widely used in China as a traditional herbal medicine due to its high economic and medicinal value. However, the appearance of Geum species is relatively similar, making identification difficult by conventional phenotypic methods, and the studies of genomics and species evolution are lacking. To better distinguish the medicinal varieties and fill this gap, we carried out relevant research on the chloroplast genome of G. japonicum. Results show a typical quadripartite structure of the chloroplast genome of G. japonicum with a length of 156,042 bp. There are totally 131 unique genes in the genome, including 87 protein-coding genes, 36 tRNA genes, and 8 rRNA genes, and there were also 87 SSRs identified and mostly mononucleotide Adenine-Thymine. We next compared the plastid genomes among four Geum species and obtained 14 hypervariable regions, including ndhF, psbE, trnG-UCC, ccsA, trnQ-UUG, rps16, psbK, trnL-UAA, ycf1, ndhD, atpA, petN, rps14, and trnK-UUU. Phylogenetic analysis revealed that G. japonicum is most closely related to Geum aleppicum, and possibly has some evolutionary relatedness with an ancient relic plant Taihangia rupestris. This research enriched the genome resources and provided fundamental insights for evolutionary studies and the phylogeny of Geum.

Complete Chloroplast Genomes of Pterodon emarginatus Vogel and Pterodon pubescens Benth: Comparative and Phylogenetic Analyses

Background

The species Pterodon emarginatus and P. pubescens, popularly known as white sucupira or faveira, are native to the Cerrado biome and have the potential for medicinal use and reforestation. They are sister species with evolutionary proximity.

Objective

Considering that the chloroplast genome exhibits a conserved structure and genes, the analysis of its sequences can contribute to the understanding of evolutionary, phylogenetic, and diversity issues.

Methods

The chloroplast genomes of P. emarginatus and P. pubescens were sequenced on the Illumina MiSeq platform. The genomes were assembled based on the de novo strategy. We performed the annotation of the genes and the repetitive regions of the genomes. The nucleotide diversity and phylogenetic relationships were analyzed using the gene sequences of these species and others of the Leguminosae family, whose genomes are available in databases.

Results

The complete chloroplast genome of P. emarginatus is 159,877 bp, and that of P. pubescens is 159,873 bp. The genomes of both species have circular and quadripartite structures. A total of 127 genes were predicted in both species, including 110 single-copy genes and 17 duplicated genes in the inverted regions. 141 microsatellite regions were identified in P. emarginatus and 140 in P. pubescens. The nucleotide diversity estimates of the gene regions in twenty-one species of the Leguminosae family were 0.062 in LSC, 0.086 in SSC, and 0.036 in IR. The phylogenetic analysis demonstrated the proximity between the genera Pterodon and Dipteryx, both from the clade Dipterygeae. Ten pairs of primers with potential for the development of molecular markers were designed.

Conclusion

The genetic information obtained on the chloroplast genomes of P. emarginatus and P. pubescens presented here reinforces the similarity and evolutionary proximity between these species, with a similarity percentage of 99.8%.

Assembly, annotation and analysis of the chloroplast genome of the Algarrobo tree Neltuma pallida (subfamily: Caesalpinioideae)

BACKGROUND

Neltuma pallida is a tree that grows in arid soils in northwestern Peru. As a predominant species of the Equatorial Dry Forest ecoregion, it holds significant economic and ecological value for both people and environment. Despite this, the species is severely threatened and there is a lack of genetic and genomic research, hindering the proposal of evidence-based conservation strategies.

RESULTS

In this work, we conducted the assembly, annotation, analysis and comparison of the chloroplast genome of a N. pallida specimen with those of related species. The assembled chloroplast genome has a length of 162,381 bp with a typical quadripartite structure (LSC-IRA-SSC-IRB). The calculated GC content was 35.97%. However, this is variable between regions, with a higher GC content observed in the IRs. A total of 132 genes were annotated, of which 19 were duplicates and 22 contained at least one intron in their sequence. A substantial number of repetitive sequences of different types were identified in the assembled genome, predominantly tandem repeats (> 300). In particular, 142 microsatellites (SSR) markers were identified. The phylogenetic reconstruction showed that N. pallida grouped with the other Neltuma species and with Prosopis cineraria. The analysis of sequence divergence between the chloroplast genome sequences of N. pallida, N. juliflora, P. farcta and Strombocarpa tamarugo revealed a high degree of similarity.

CONCLUSIONS

The N. pallida chloroplast genome was found to be similar to those of closely related species. With a size of 162,831 bp, it had the classical chloroplast quadripartite structure and GC content of 35.97%. Most of the 132 identified genes were protein-coding genes. Additionally, over 800 repetitive sequences were identified, including 142 SSR markers. In the phylogenetic analysis, N. pallida grouped with other Neltuma spp. and P. cineraria. Furthermore, N. pallida chloroplast was highly conserved when compared with genomes of closely related species. These findings can be of great potential for further diversity studies and genetic improvement of N. pallida.

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.

The complete chloroplast genome sequences of six Hylotelephium species: Comparative genomic analysis and phylogenetic relationships

To evaluate the phylogenetic relationships between Hylotelephium and Orostachys, and to provide important information for further studies, we analyzed the complete chloroplast genomes of six Hylotelephium species and compared the sequences to those of published chloroplast genomes of congeneric species and species of the closely related genus, Orostachys. The total chloroplast genome length of nineteen species, including the six Hylotelephium species analyzed in this study and the thirteen Hylotelephium and Orostachys species analyzed in previous studies, ranged from 150,369 bp (O. minuta) to 151,739 bp (H. spectabile). Their overall GC contents were almost identical (37.7-37.8%). The chloroplast genomes of the nineteen species contained 113 unique genes comprising 79 protein-coding genes (PCGs), 30 transfer RNA genes (tRNAs), and four ribosomal RNA genes (rRNAs). Among the annotated genes, fourteen genes contained one intron, and two genes contained two introns. The chloroplast genomes of the nineteen Hylotelephium and Orostachys species had identical structures. Additionally, the large single copy (LSC), inverted repeat (IR), and small single copy (SSC) junction regions were conserved in the Hylotelephium and Orostachys species. The nucleotide diversity between the Hylotelephium chloroplast genomes was extremely low in all regions, and only one region showed a high Pi value (>0.03). In all nineteen chloroplast genomes, six regions had a high Pi value (>0.03). The phylogenetic analysis showed that the genus delimitation could not be clearly observed even in this study because Hylotelephium formed a paraphyly with subsect. Orostachys of the genus Orostachys. Additionally, the data supported the taxonomic position of Sedum taqeutii, which was treated as a synonym for H. viridescens in previous studies, as an independent taxon.

Chloroplast genome assembly of Serjania erecta Raldk: comparative analysis reveals gene number variation and selection in protein-coding plastid genes of Sapindaceae

Serjania erecta Raldk is an essential genetic resource due to its anti-inflammatory, gastric protection, and anti-Alzheimer properties. However, the genetic and evolutionary aspects of the species remain poorly known. Here, we sequenced and assembled the complete chloroplast genome of S. erecta and used it in a comparative analysis within the Sapindaceae family. S. erecta has a chloroplast genome (cpDNA) of 159,297 bp, divided into a Large Single Copy region (LSC) of 84,556 bp and a Small Single Copy region (SSC) of 18,057 bp that are surrounded by two Inverted Repeat regions (IRa and IRb) of 28,342 bp. Among the 12 species used in the comparative analysis, S. erecta has the fewest long and microsatellite repeats. The genome structure of Sapindaceae species is relatively conserved; the number of genes varies from 128 to 132 genes, and this variation is associated with three main factors: (1) Expansion and retraction events in the size of the IRs, resulting in variations in the number of rpl22, rps19, and rps3 genes; (2) Pseudogenization of the rps2 gene; and (3) Loss or duplication of genes encoding tRNAs, associated with the duplication of trnH-GUG in X. sorbifolium and the absence of trnT-CGU in the Dodonaeoideae subfamily. We identified 10 and 11 mutational hotspots for Sapindaceae and Sapindoideae, respectively, and identified six highly diverse regions (tRNA-Lys - rps16, ndhC - tRNA-Val, petA - psbJ, ndhF, rpl32 - ccsA, and ycf1) are found in both groups, which show potential for the development of DNA barcode markers for molecular taxonomic identification of Serjania. We identified that the psaI gene evolves under neutrality in Sapindaceae, while all other chloroplast genes are under strong negative selection. However, local positive selection exists in the ndhF, rpoC2, ycf1, and ycf2 genes. The genes ndhF and ycf1 also present high nucleotide diversity and local positive selection, demonstrating significant potential as markers. Our findings include providing the first chloroplast genome of a member of the Paullinieae tribe. Furthermore, we identified patterns in variations in the number of genes and selection in genes possibly associated with the family's evolutionary history.

Chloroplast genome characterization of Uncaria guianensis and Uncaria tomentosa and evolutive dynamics of the Cinchonoideae subfamily

Uncaria species are used in traditional medicine and are considered of high therapeutic value and economic importance. This work describes the assembly and annotation of the chloroplast genomes of U. guianensis and U. tomentosa, as well as a comparative analysis. The genomes were sequenced on MiSeq Illumina, assembled with NovoPlasty, and annotated using CHLOROBOX GeSeq. Addictionaly, comparative analysis were performed with six species from NCBI databases and primers were designed in Primer3 for hypervariable regions based on the consensus sequence of 16 species of the Rubiaceae family and validated on an in-silico PCR in OpenPrimeR. The genome size of U. guianensis and U. tomentosa was 155,505 bp and 156,390 bp, respectively. Both Species have 131 genes and GC content of 37.50%. The regions rpl32-ccsA, ycf1, and ndhF-ccsA showed the three highest values of nucleotide diversity within the species of the Rubiaceae family and within the Uncaria genus, these regions were trnH-psbA, psbM-trnY, and rps16-psbK. Our results indicates that the primer of the region ndhA had an amplification success for all species tested and can be promising for usage in the Rubiaceae family. The phylogenetic analysis recovered a congruent topology to APG IV. The gene content and the chloroplast genome structure of the analyzed species are conserved and most of the genes are under negative selection. We provide the cpDNA of Neotropical Uncaria species, an important genomic resource for evolutionary studies of the group.

The complete chloroplast genome sequence of Eugenia klotzschiana O. Berg unveils the evolutionary dynamics in plastomes of Myrteae DC. tribe (Myrtaceae)

Myrteae is the most diversified tribe in the Myrtaceae family and has great ecological and economic importance. Here, we performed the assembly and annotation of the chloroplast genome of Eugenia klotzschiana O. Berg and used this in a comparative analysis with other 13 species from the Myrteae tribe. The E. klotzschiana plastome exhibited a length of 158,977 bp and a very conserved structure and gene composition when compared with other Myrteae genomes. We identified 34 large repetitive sequences and 94 SSR repeats in E. klotzschiana plastome. The trnT-trnL, rpl32-trnL, ndhF-rpl32, psbE-petL, and ycf1 regions were identified as mutational hotspots. A negative selection signal was detected in 74 protein-coding genes while neutral selection was detected in two genes (rps12 and psaI). Furthermore, 222 RNA editing sites were identified in the E. klotzschiana plastome. We also obtained a plastome-based Myrtales phylogenetic tree, including E. klotzschiana for the first time in a molecular phylogeny, recovering its sister relationship for all other Eugenia species. Our results illuminate how evolution shaped the chloroplast genome structure and composition in the Myrteae tribe, especially in the E. klotzschiana plastome.

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.

Plastid Phylogenomic Analyses Reveal the Taxonomic Position of Peucedanum franchetii

Peucedanum franchetii is a famous folk medicinal plant in China. However, the taxonomy of the P. franchetii has not been sufficiently resolved. Due to similar morphological features between P. franchetii and Ligusticopsis members, the World Flora Online (WFO) Plant List suggested that this species transformed into the genus Ligusticopsis and merged with Ligusticopsis likiangensis. However, both species are obviously diverse in leaf shape, bracts, and bracteoles. To check the taxonomic position of P. franchetii, we newly sequenced and assembled the plastome of P. franchetii and compared it with nine other plastomes of the genus Ligusticopsis. Ten plastomes were highly conserved and similar in gene order, codon bias, RNA editing sites, IR borders, and SSRs. Nevertheless, 10 mutation hotspot regions (infA, rps8, matK, ndhF, rps15, psbA-trnH, rps2-rpoC2, psbA-trnK, ycf2-trnL, and ccsA-ndhD) were still detected. In addition, both phylogenetic analyses based on plastome data and ITS sequences robustly supported that P. franchetii was not clustered with members of Peucedanum but nested in Ligusticopsis. P. franchetii was sister to L. likiangensis in the ITS topology but clustered with L. capillacea in the plastome tree. These findings implied that P. franchetii should be transferred to genus Ligusticopsis and not merged with L. likiangensis, but as an independent species, which was further verified by morphological evidences. Therefore, transferring P. franchetii under the genus Ligusticopsis as an independent species was reasonable, and a new combination was presented.

The complete chloroplast genome sequences of eight Orostachys species: Comparative analysis and assessment of phylogenetic relationships

We analyzed the complete chloroplast genomes of eight Orostachys species and compared the sequences to those of published chloroplast genomes of the congeneric and closely related genera, Meterostachys and Hylotelephium. The total chloroplast genome length of thirteen species, including the eight species analyzed in this study and the five species analyzed in previous studies, ranged from 149,860 (M. sikokianus) to 151,707 bp (H. verticillatum). The overall GC contents of the genomes were almost identical (37.6 to 37.8%). The thirteen chloroplast genomes each contained 113 unique genes comprising 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. Among the annotated genes, sixteen genes contained one or two introns. Although the genome structures of all Orostachys and Hylotelephium species were identical, Meterostachys differed in structure due to a relatively large gene block (trnS-GCU-trnS-GGA) inversion. The nucleotide diversity among the subsect. Orostachys chloroplast genomes was extremely low in all regions, and among the subsect. Appendiculatae, genus Orostachys, and all thirteen chloroplast genomes showed high values of Pi (>0.03) in one, five, or three regions. The phylogenetic analysis showed that Orostachys formed polyphyly, and subsect. Orostachys and Appendiculatae were clustered with Hylotelephium and Meterostachys, respectively, supporting the conclusion that each subsection should be considered as an independent genus. Furthermore, the data supported the taxonomic position of O. margaritifolia and O. iwarenge f. magnus, which were treated as synonyms for O. iwarenge in a previous study, as independent taxa. Our results suggested that O. ramosa and O. japonica f. polycephala were individual variations of O. malacophylla and O. japonica, respectively. The exact taxonomic position of O. latielliptica and the phylogenetic relationship among the three species, O. chongsunensis, O. malacophylla and O. ramosa, should be a topic of future study.

Comprehensive Comparative Analysis and Development of Molecular Markers for Dianthus Species Based on Complete Chloroplast Genome Sequences

Dianthus spp. is a genus with high economic and ornamental value in the Caryophyllaceae, which include the famous fresh-cut carnation and the traditional Chinese herbal medicine, D. superbus. Despite the Dianthus species being seen everywhere in our daily lives, its genome information and phylogenetic relationships remain elusive. Thus, we performed the assembly and annotation of chloroplast genomes for 12 individuals from seven Dianthus species. On this basis, we carried out the first comprehensive and systematic analysis of the chloroplast genome sequence characteristics and the phylogenetic evolution of Dianthus. The chloroplast genome of 12 Dianthus individuals ranged from 149,192 bp to 149,800 bp, containing 124 to 126 functional genes. Sequence repetition analysis showed the number of simple sequence repeats (SSRs) ranged from 75 to 80, tandem repeats ranged from 23 to 41, and pair-dispersed repeats ranged from 28 to 43. Next, we calculated the synonymous nucleotide substitution rates (Ks) of all 76 protein coding genes to obtain the evolution rate of these coding genes in Dianthus species; rpl22 showed the highest Ks (0.0471), which suggested that it evolved the swiftest. By reconstructing the phylogenetic relationships within Dianthus and other species of Caryophyllales, 16 Dianthus individuals (12 individuals reported in this study and four individuals downloaded from NCBI) were divided into two strongly supported sister clades (Clade A and Clade B). The Clade A contained five species, namely D. caryophyllus, D. barbatus, D. gratianopolitanus, and two cultivars (‘HY’ and ‘WC’). The Clade B included four species, in which D. superbus was a sister branch with D. chinensis, D. longicalyx, and F1 ‘87M’ (the hybrid offspring F1 from D. chinensis and ‘HY’). Further, based on sequence divergence analysis and hypervariable region analysis, we selected several regions that had more divergent sequences, to develop DNA markers. Additionally, we found that one DNA marker can be used to differentiate Clade A and Clade B in Dianthus. Taken together, our results provide useful information for our understanding of Dianthus classification and chloroplast genome evolution.

Complete chloroplast genomes and comparative analyses of Hippeastrum ‘milady’, Hippeastrum albertii and Hippeastrum reticulatum (Amaryllidaceae)

Hippeastrum is a genus of ornamental plants with large, brightly colored flowers. Due to the very high seed-setting rate of the hybridization of Hippeastrum, the large population of hybrid progeny and the existence of superparent inheritance, it is difficult to trace the origin of the varieties collected from the market during breeding. In this study, we analyzed the chloroplast genomes of Hippeastrum ‘Milady’, H. alberti, and H. reticulatum using the Illumina NovaSeq sequencing platform and generated full-length sequences of 158,067, 158,067, and 158,522 bp, respectively. All three genomes had the typical tetrad structure. The large single copy, small single copy, and inverted repeat regions of H. reticulatum were observed to be respectively 277, 138, and 20 bp longer than the corresponding regions of H. ‘Milady’ and H. alberti. The results of comparative analysis of simple sequence repeats (SSRs), Ka/Ks ratios, codon preferences, and complete sequences of chloroplasts of these three taxa and 14 other plant species were as follows. First, the chloroplast genomes of H. ‘Milady’, H. alberti, and H. reticulatum contain 209, 209, and 211 SSR sites, respectively, most of which (123, 123, and 122, respectively) are single nucleotide repeats. Second, leucine, arginine, and serine are the most frequently used amino acids in the three chloroplast genomes. Third, H. ‘Milady’, H. alberti, and H. reticulatum are more closely related to Lycoris and Narcissus than to Allium and Agapanthus. Our results will provide information on the study of origins or relatedness of native species, and the identification of cultivars.

The complete chloroplast genome of critically endangered Chimonobambusa hirtinoda (Poaceae: Chimonobambusa) and phylogenetic analysis

Chimonobambusa hirtinoda, a threatened species, is only naturally distributed in Doupeng Mountain, Duyun, Guizhou, China. Next-generation sequencing (NGS) is used to obtain the complete chloroplast (cp) genome sequence of C. hirtinoda. The sequence was assembled and analyzed for phylogenetic and evolutionary studies. Additionally, we compared the cp genome of C. hirtinoda with previously published Chimonobambusa species. The cp genome of C. hirtinoda has a total length of 139, 561 bp and 38.90% GC content. This genome included a large single -copy (LSC) region of 83, 166 bp, a small single-copy (SSC) region of 20, 811 bp and a pair of inverted repeats of 21,792 bp each. We discovered 130 genes in the cp genome, including 85 protein-coding genes, 37 tRNA, and 8 rRNA genes. A total of 48 simple sequence repeats (SSRs) were detected. The A/U preference of the third nucleotide in the cp genome of C. hirtinoda was obtained by measuring the codon usage frequency of amino acids. Furthermore, phylogenetic analysis using complete cp sequences and matK gene revealed a genetic relationship within the Chimonobambusa genus. This study reported the chloroplast genome of the C. hirtinoda.

The chloroplast genome sequence and phylogenetic analysis of Apocynum venetum L

Apocynum venetum L. (Apocynaceae) is valuable for its medicinal compounds and fiber content. Native A. venetum populations are threatened and require protection. Wild A. venetum resources are limited relative to market demand and a poor understanding of the composition of A. venetum at the molecular level. The chloroplast genome contains genetic markers for phylogenetic analysis, genetic diversity evaluation, and molecular identification. In this study, the entire genome of the A. venetum chloroplast was sequenced and analyzed. The A. venetum cp genome is 150,878 bp, with a pair of inverted repeat regions (IRA and IRB). Each inverted repeat region is 25,810 bp, which consist of large (LSC, 81,951 bp) and small (SSC, 17,307 bp) single copy areas. The genome-wide GC content was 38.35%, LSC made up 36.49%, SSC made up 32.41%, and IR made up 43.3%. The A. venetum chloroplast genome encodes 131 genes, including 86 protein-coding genes, eight ribosomal RNA genes, and 37 transfer RNA genes. This study identified the unique characteristics of the A. venetum chloroplast genome, which will help formulate effective conservation and management strategies as well as molecular identification approaches for this important medicinal plant.

The complete chloroplast genome of Onobrychis gaubae (Fabaceae-Papilionoideae): comparative analysis with related IR-lacking clade species

BACKGROUND

Plastome (Plastid genome) sequences provide valuable markers for surveying evolutionary relationships and population genetics of plant species. Papilionoideae (papilionoids) has different nucleotide and structural variations in plastomes, which makes it an ideal model for genome evolution studies. Therefore, by sequencing the complete chloroplast genome of Onobrychis gaubae in this study, the characteristics and evolutionary patterns of plastome variations in IR-loss clade were compared.

RESULTS

In the present study, the complete plastid genome of O. gaubae, endemic to Iran, was sequenced using Illumina paired-end sequencing and was compared with previously known genomes of the IRLC species of legumes. The O. gaubae plastid genome was 122,688 bp in length and included a large single-copy (LSC) region of 81,486 bp, a small single-copy (SSC) region of 13,805 bp and one copy of the inverted repeat (IR_b) of 29,100 bp. The genome encoded 110 genes, including 76 protein-coding genes, 30 transfer RNA (tRNA) genes and four ribosome RNA (rRNA) genes and possessed 83 simple sequence repeats (SSRs) and 50 repeated structures with the highest proportion in the LSC. Comparative analysis of the chloroplast genomes across IRLC revealed three hotspot genes (ycf1, ycf2, clpP) which could be used as DNA barcode regions. Moreover, seven hypervariable regions [trnL(UAA)-trnT(UGU), trnT(GGU)-trnE(UUC), ycf1, ycf2, ycf4, accD and clpP] were identified within Onobrychis, which could be used to distinguish the Onobrychis species. Phylogenetic analyses revealed that O. gaubae is closely related to Hedysarum. The complete O. gaubae genome is a valuable resource for investigating evolution of Onobrychis species and can be used to identify related species.

CONCLUSIONS

Our results reveal that the plastomes of the IRLC are dynamic molecules and show multiple gene losses and inversions. The identified hypervariable regions could be used as molecular markers for resolving phylogenetic relationships and species identification and also provide new insights into plastome evolution across IRLC.

Complete chloroplast genomes of Rubus species (Rosaceae) and comparative analysis within the genus

BACKGROUND

Rubus is the largest genus of the family Rosaceae and is valued as medicinal, edible, and ornamental plants. Here, we sequenced and assembled eight chloroplast (cp) genomes of Rubus from the Dabie Mountains in Central China. Fifty-one Rubus species were comparatively analyzed for the cp genomes including the eight newly discovered genomes and forty-three previously reported in GenBank database (NCBI).

RESULTS

The eight newly obtained cp genomes had the same quadripartite structure as the other cp genomes in Rubus. The length of the eight plastomes ranged from 155,546 bp to 156,321 bp with similar GC content (37.0 to 37.3%). The results indicated 133-134 genes were annotated for the Rubus plastomes, which contained 88 or 89 protein coding genes (PCGs), 37 transfer RNA genes (tRNAs), and eight ribosomal RNA genes (rRNAs). Among them, 16 (or 18) of the genes were duplicated in the IR region. Structural comparative analysis results showed that the gene content and order were relatively preserved. Nucleotide variability analysis identified nine hotspot regions for genomic divergence and multiple simple sequences repeats (SSRs), which may be used as markers for genetic diversity and phylogenetic analysis. Phylogenetic relationships were highly supported within the family Rosaceae, as evidenced by sub-clade taxa cp genome sequences.

CONCLUSION

Thus, the whole plastome may be used as a super-marker in phylogenetic studies of this genus.

Heterogeneous Genetic Diversity Estimation of a Promising Domestication Medicinal Motherwort Leonurus Cardiaca Based on Chloroplast Genome Resources

Leonurus cardiaca has a long history of use in western herbal medicine and is applied for the treatment of gynaecological conditions, anxiety, and heart diseases. Because of its botanical relationship to the primary Chinese species, L. japonicus, and extensive medical indications that go beyond the traditional indications for the Chinese species, it is a promising medicinal resource. Therefore, the features of genetic diversity and variability in the species have been prioritized. To explore these issues, we sequenced the chloroplast genomes of 22 accessions of L. cardiaca from different geographical locations worldwide using high-throughput sequencing. The results indicate that L. cardiaca has a typical quadripartite structure and range from 1,51,236 bp to 1,51,831 bp in size, forming eight haplotypes. The genomes all contain 114 distinct genes, including 80 protein-coding genes, 30 transfer RNA genes and four ribosomal RNA genes. Comparative analysis showed abundant diversity of single nucleotide polymorphisms (SNPs), indels, simple sequence repeats (SSRs) in 22 accessions. Codon usage showed highly similar results for L. cardiaca species. The phylogenetic and network analysis indicated 22 accessions forming four clades that were partly related to the geographical distribution. In summary, our study highlights the advantage of chloroplast genome with large data sets in intraspecific diversity evaluation and provides a new tool to facilitate medicinal plant conservation and domestication.

The complete chloroplast genome of Artocarpus altilis (Moraceae) and phylogenetic relationships

The chloroplast (cp) is an essential organelle in higher plants. The genes of the plastome are well suited to infer phylogenetic relationships among taxa. In this study, we report the assembly of the cp genome of Artocarpus altilis and its phylogeny among species from Moraceae family. The cp genome of A. altilis was 160,822 base pair (bp) in length, comprising one large single-copy region of 88,692 bp, one small single-copy region of 19,290 bp, and a pair of inverted repeat regions (IRs) of 26,420 bp. A total of 113 different genes were predicted, including 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. The phylogenetic analysis of 19 species belonging to the Moraceae family confirmed the phylogenetic proximity of the genus Artocarpus and Morus and the genetic similarity of A. camansi and A. altilis.

Plastome evolution in the Caesalpinia group (Leguminosae) and its application in phylogenomics and populations genetics

The chloroplast genomes of Caesalpinia group species are structurally conserved, but sequence level variation is useful for both phylogenomic and population genetic analyses. Variation in chloroplast genomes (plastomes) has been an important source of information in plant biology. The Caesalpinia group has been used as a model in studies correlating ecological and genomic variables, yet its intergeneric and infrageneric relationships are not fully solved, despite densely sampled phylogenies including nuclear and plastid loci by Sanger sequencing. Here, we present the de novo assembly and characterization of plastomes from 13 species from the Caesalpinia group belonging to eight genera. A comparative analysis was carried out with 13 other plastomes previously available, totalizing 26 plastomes and representing 15 of the 26 known Caesalpinia group genera. All plastomes showed a conserved quadripartite structure and gene repertoire, except for the loss of four ndh genes in Erythrostemon gilliesii. Thirty polymorphic regions were identified for inter- or intrageneric analyses. The 26 aligned plastomes were used for phylogenetic reconstruction, revealing a well-resolved topology, and dividing the Caesalpinia group into two fully supported clades. Sixteen microsatellite (cpSSR) loci were selected from Cenostigma microphyllum for primer development and at least two were cross-amplified in different Leguminosae subfamilies by in vitro or in silico approaches. Four loci were used to assess the genetic diversity of C. microphyllum in the Brazilian Caatinga. Our results demonstrate the structural conservation of plastomes in the Caesalpinia group, offering insights into its systematics and evolution, and provides new genomic tools for future phylogenetic, population genetics, and phylogeographic studies.

Complete chloroplast genome sequence and phylogenetic analysis of Torenia fournieri

Torenia fournieri belongs to the genus Torenia in the family Linderniaceae. The complete chloroplast genome of T. fournieri was sequenced and analyzed by Illumina sequencing in this study. The full length of the complete chloroplast genome is 153,938 bp, containing a pair of inverted repeat regions of 24,805 bp (IRa and IRb) separated by a large single copy region (LSC) of 85,498 bp and a small single copy region (SSC) of 18,830 bp. The T. fournieri chloroplast genome encodes 131 genes, comprising 87 protein-coding genes, 36 tRNA genes, 8 rRNA genes, without pseudogene. Phylogenetic analysis showed that T. fournieri was closely related to T. benthamiana and T. concolor within the genus Torenia in family Linderniaceae.

Complete Chloroplast Genome of Clethra fargesii Franch., an Original Sympetalous Plant from Central China: Comparative Analysis, Adaptive Evolution, and Phylogenetic Relationships

Clethra fargesii, an essential ecological and endemic woody plant of the genus Clethra in Clethraceae, is widely distributed in Central China. So far, there have been a paucity of studies on its chloroplast genome. In the present study, we sequenced and assembled the complete chloroplast genome of C. fargesii. We also analyzed the chloroplast genome features and compared them to Clethra delavayi and other closely related species in Ericales. The complete chloroplast genome is 157,486 bp in length, including a large single-copy (LSC) region of 87,034 bp and a small single-copy (SSC) region of 18,492 bp, separated by a pair of inverted repeat (IR) regions of 25,980 bp. The GC content of the whole genome is 37.3%, while those in LSC, SSC, and IR regions are 35.4%, 30.7%, and 43.0%, respectively. The chloroplast genome of C. fargesii encodes 132 genes in total, including 87 protein-coding genes (PCGs), 37 tRNA genes, and eight rRNA genes. A total of 26,407 codons and 73 SSRs were identified in C. fargesii chloroplast genome. Additionally, we postulated and demonstrated that the structure of the chloroplast genome in Clethra species may present evolutionary conservation based on the comparative analysis of genome features and genome alignment among eight Ericales species. The low Pi values revealed evolutionary conservation based on the nucleotide diversity analysis of chloroplast genome in two Clethra species. The low selection pressure was shown by a few positively selected genes by adaptive evolution analysis using 80 coding sequences (CDSs) of the chloroplast genomes of two Clethra species. The phylogenetic tree showed that Clethraceae and Ericaceae are sister clades, which reconfirm the previous hypothesis that Clethra is highly conserved in the chloroplast genome using 75 CDSs of chloroplast genome among 40 species. The genome information and analysis results presented in this study are valuable for further study on the intraspecies identification, biogeographic analysis, and phylogenetic relationship in Clethraceae.

Complete chloroplast genome of Myracrodruon urundeuva and its phylogenetics relationships in Anacardiaceae family

Continuous exploratory use of tree species is threatening the existence of several plants in South America. One of these threatened species is Myracroduron urundeuva, highly exploited due to the high quality and durability of its wood. The chloroplast (cp) has been used for several evolutionary studies as well traceability of timber origin, based on its gene sequences and simple sequence repeats (SSR) variability. Cp genome organization is usually consisting of a large single copy and a small single copy region separated by two inverted repeats regions. We sequenced the complete cp genome from M. urundeuva based on Illumina next-generation sequencing. Our results show that the cp genome is 159,883 bp in size. The 36 SSR identified ranging from mono- to hexanucleotides. Positive selection analysis revealed nine genes related to photosystem, protein synthesis, and DNA replication, and protease are under positive selection. Genome comparison a other Anacardiaceae chloroplast genomes showed great variability in the family. The phylogenetic analysis using complete chloroplast genome sequences of other Anacardiaceae family members showed a close relationship with two other economically important genera, Pistacia and Rhus. These results will help future investigations of timber monitoring and population and evolutionary studies.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-00989-1.

Complete chloroplast genome sequence of Adenophora racemosa (Campanulaceae): Comparative analysis with congeneric species

Adenophora racemosa, belonging to the Campanulaceae, is an important species because it is endemic to Korea. The goal of this study was to assemble and annotate the chloroplast genome of A. racemosa and compare it with published chloroplast genomes of congeneric species. The chloroplast genome was reconstructed using de novo assembly of paired-end reads generated by the Illumina MiSeq platform. The chloroplast genome size of A. racemosa was 169,344 bp. In total, 112 unique genes (78 protein-coding genes, 30 tRNAs, and 4 rRNAs) were identified. A Maximum likelihood (ML) tree based on 76 protein-coding genes divided the five Adenophora species into two clades, showing that A. racemosa is more closely related to Adenophora stricta than to Adenophora divaricata. The gene order and contents of the LSC region of A. racemosa were identical to those of A. divaricata and A. stricta, but the structure of the SSC and IRs was unique due to IR contraction. Nucleotide diversity (Pi) >0.05 was found in eleven regions among the three Adenophora species not included in sect. Remotiflorae and in six regions between two species (A. racemosa and A. stricta).

The complete chloroplast genome sequence of Camellia chuongtsoensis

Camellia chuongtsoensis is an evergreen shrub with a single-petaled flower and golden yellow color. The complete chloroplast genome of C. chuongtsoensis was sequenced and analyzed in this study by Illumina sequencing. The chloroplast genome is 156,504 bp in length with a quadripartite structure containing a large single copy (LSC) region of 86,215 bp, a small single copy (SSC) region of 18,253 bp, and a pair of inverted repeat regions of 26,018 bp (IRa and IRb). The chloroplast genome of C. chuongtsoensis encodes 135 genes, comprising 87 protein-coding genes, 37 tRNA genes, 8 rRNA genes, and 3 pseudogenes.

Complete chloroplast genome sequence of Camellia rhytidophylla, comparative and phylogenetic analysis

Camellia rhytidophylla is an endangered plant with economic value. Using Illumina sequencing, the chloroplast genome of C. rhytidophylla was sequenced and analyzed in this study. The complete chloroplast genome is 157,073 bp in length, which consisted of a pair of inverted repeat regions of 26,055 bp (IRa and IRb) separated by a large single-copy region (LSC) of 86,680 bp and a small single-copy region (SSC) of 18,283 bp. The C. rhytidophylla chloroplast genome encodes 135 genes, including 87 protein-coding genes, 37 tRNA genes, 8 rRNA genes, and 3 pseudogenes. Sequence comparison analysis with the chloroplast sequences of 28 other Camellia plants found that C. rhytidophylla had the closest relationship with C. szechuanensis. This study provides a theoretical basis for the analysis of the distant relationship of Camellia.

Complete chloroplast genome sequence and phylogenetic analysis of Camellia fraterna

Camellia fraterna belongs to the genus Camellia in the family Theaceae. We sequenced and analyzed the complete chloroplast genome of C. fraterna by Illumina sequencing in this study. The full length of the complete chloroplast genome is 156,902 bp, containing a pair of inverted repeat regions of 26,030 bp (IRa and IRb) separated by a large single-copy (LSC) region of 86,583 bp and a small single-copy (SSC) region of 18,259 bp. The C. fraterna chloroplast genome encodes 135 genes, comprising 87 protein-coding genes, 37 tRNA genes, eight rRNA genes, and three pseudogenes. This study will be useful for further study on genetic diversity and molecular breeding.

Simple innovative adaptor to improve genome walking with convenient PCR

BACKGROUND

Various polymerase chain reaction (PCR)-based methods have been applied for the development of genome walking (GW) technique. These methods which could be based on the application of restriction enzymes or primers have various efficiencies to identify the unknown nucleotide sequences. The present study was conducted to design a new innovative double-strand adaptor using MAP30 gene sequence of Momordica charantia plant as a model to improve genome walking with convenient PCR.

RESULTS

The adaptor was designed using multiple restriction sites of Hind III, BamH I, EcoR I, and Bgl II enzymes with no restriction site in a known sequence of the MAP30 gene. In addition, no modification was required to add phosphate, amine, or other groups to the adaptor, since restriction enzyme digestion of double-strand adaptor provided the 5' phosphate group. Here, preparation of the phosphate group in the genomic DNA of the plant digestion with restriction enzymes was performed followed by ligation with digested adaptor containing 5' phosphate group.

CONCLUSION

PCR was done to amplify the unknown sequence using MAP30 gene-specific primer and adaptor primer. Results confirmed the ability of the technique for successful identification of the sequence. Consequently, a newly designed adaptor in the developed technique reduced the time and cost of the method compared to the conventional genome walking; also, cloning and culturing of bacterial steps could be eliminated.

The First Glimpse of Streptocarpus ionanthus (Gesneriaceae) Phylogenomics: Analysis of Five Subspecies' Chloroplast Genomes

Streptocarpus ionanthus (Gesneriaceae) comprise nine herbaceous subspecies, endemic to Kenya and Tanzania. The evolution of Str. ionanthus is perceived as complex due to morphological heterogeneity and unresolved phylogenetic relationships. Our study seeks to understand the molecular variation within Str. ionanthus using a phylogenomic approach. We sequence the chloroplast genomes of five subspecies of Str. ionanthus, compare their structural features and identify divergent regions. The five genomes are identical, with a conserved structure, a narrow size range (170 base pairs (bp)) and 115 unique genes (80 protein-coding, 31 tRNAs and 4 rRNAs). Genome alignment exhibits high synteny while the number of Simple Sequence Repeats (SSRs) are observed to be low (varying from 37 to 41), indicating high similarity. We identify ten divergent regions, including five variable regions (psbM, rps3, atpF-atpH, psbC-psbZ and psaA-ycf3) and five genes with a high number of polymorphic sites (rps16, rpoC2, rpoB, ycf1 and ndhA) which could be investigated further for phylogenetic utility in Str. ionanthus. Phylogenomic analyses here exhibit low polymorphism within Str. ionanthus and poor phylogenetic separation, which might be attributed to recent divergence. The complete chloroplast genome sequence data concerning the five subspecies provides genomic resources which can be expanded for future elucidation of Str. ionanthus phylogenetic relationships.

Complete Chloroplast Genomes of Chlorophytum comosum and Chlorophytum gallabatense: Genome Structures, Comparative and Phylogenetic Analysis

The genus Chlorophytum includes many economically important species well-known for medicinal, ornamental, and horticultural values. However, to date, few molecular genomic resources have been reported for this genus. Therefore, there is limited knowledge of phylogenetic studies, and the available chloroplast (cp) genome of Chlorophytum (C. rhizopendulum) does not provide enough information on this genus. In this study, we present genomic resources for C. comosum and C. gallabatense, which had lengths of 154,248 and 154,154 base pairs (bp), respectively. They had a pair of inverted repeats (IRa and IRb) of 26,114 and 26,254 bp each in size, separating the large single-copy (LSC) region of 84,004 and 83,686 bp from the small single-copy (SSC) region of 18,016 and 17,960 bp in C. comosum and C. gallabatense, respectively. There were 112 distinct genes in each cp genome, which were comprised of 78 protein-coding genes, 30 tRNA genes, and four rRNA genes. The comparative analysis with five other selected species displayed a generally high level of sequence resemblance in structural organization, gene content, and arrangement. Additionally, the phylogenetic analysis confirmed the previous phylogeny and produced a phylogenetic tree with similar topology. It showed that the Chlorophytum species (C. comosum, C. gallabatense and C. rhizopendulum) were clustered together in the same clade with a closer relationship than other plants to the Anthericum ramosum. This research, therefore, presents valuable records for further molecular evolutionary and phylogenetic studies which help to fill the gap in genomic resources and resolve the taxonomic complexes of the genus.