Comparative analyses of chloroplast genomes among three Firmiana species: Identification of mutational hotspots and phylogenetic relationship with other species of Malvaceae

Chloroplast genomics unravels taxonomic status and evolutionary relationships in Euphorbia fauriei and E. pekinensis

Euphorbia fauriei, a perennial plant endemic to South Korea, exhibits both morphological characteristics and intricate genetic identities akin to E. pekinensis. This study aimed to provide clarity on the taxonomic status of E. faurirei and E. pekinensis through a comprehensive chloroplast genome analysis. Additionally, we sequenced the Acalypha australis chloroplast genome as an outgroup for the construction of a phylogenetic tree with other Euphorbia species. The three chloroplast genomes, ranging from 162,834 bp to 168,832 bp, displayed typical quadripartite structures. The Euphorbia genomes contained 111 unique genes, whereas the A. australis genome contained two additional protein-coding genes (PCGs), rpl32 and rps16. Comparative analysis unveiled the loss of rpl32 and rps16 as synapomorphic characteristics in Euphorbia, whereas the loss of infA occurred across Euphorbiaceae. High collinearity and sequence similarity were observed among Euphorbia species, accompanied by significant inversion regions in the E. pekinensis chloroplast genomes from China and Japan, indicating regional genetic variability. Nucleotide substitution analysis revealed that the ndh group exhibited the highest KA/KS values (0.224), whereas the atp, psa, psb and rps groups had the lowest. Phylogenomic analysis utilizing whole genomes, PCGs, and intron and intergenic regions consistently demonstrated that E. pekinensis from South Korea clusters closely with E. fauriei. These findings challenge the current taxonomic distinction between E. pekinensis and E. fauriei in Korea, suggesting that while they exhibit distinct characteristics, E. fauriei should be considered a closely related subspecies rather than the same species as E. pekinensis. This study emphasizes the need for population studies to clarify the taxonomic relationships between E. pekinensis and E. fauriei.

Chloroplast genome evolution of Berberis (Berberidaceae): Implications for phylogeny and metabarcoding

Berberidis Radix (Sankezhen), a typical multi-origin Chinese medicinal material, originates from the dried roots of plants of the Berberis genus and is used to treat various ailments. These species have similar morphologies, potentially leading to misidentifications that can impact medicine efficacy. Therefore, developing suitable molecular markers to identify medicinal species is imperative. Furthermore, discrepancies exist in the taxonomy of the Berberis genus. In the present study, we de novo assembled the chloroplast genomes of six Berberis species (Berberis woomungensis C. Y. Wu, Berberis pruinosa Franch., Berberis thunbergii DC., Berberis chinensis Poir., Berberis wilsoniae Hemsl., and Berberis sp.) that commonly constitute Berberidis Radix and compared them with previously reported genomes. Our comparative analysis revealed similarities in genome structure, relative synonymous codon usage, amino acid frequency, repeats, and substitutions. Higher synonymous substitutions, indicative of predominant purifying selection on protein-coding genes, were observed compared to non-synonymous substitutions. However, positive selection was identified in six genes across 29 Berberis species-accD, matK, ndhD, rbcL, ycf1, and ycf2-highlighting their potential roles in adaptive responses to specific environmental conditions within the genus. Inverted repeats expansion and contraction affected the rate of mutations and were associated with the phylogenetic classification of Berberis. Our phylogenetic analysis supported the division of the Berberis complex into four genera, which corroborates previous studies involving extensive sampling. We identified the ndhD-ccsA region as the most polymorphic region and applied this region to Chinese patent medicines containing Berberidis Radix through metabarcoding. The metabarcoding analysis confirmed that five Berberis species commonly constitute Berberidis Radix in Chinese patent medicines. In conclusion, this study provides insight into the molecular evolution of the chloroplast genome and the phylogeny of the Berberis genus. In addition, metabarcoding provides insight into the species composition of Berberidis Radix in Chinese patent medicines.

Complete chloroplast genome sequence of Artemisia littoricola (Asteraceae) from Dokdo Island Korea: genome structure, phylogenetic analysis, and biogeography study

The Asteraceae family, particularly the Artemisia genus, presents taxonomic challenges due to limited morphological characteristics and frequent natural hybridization. Molecular tools, such as chloroplast genome analysis, offer solutions for accurate species identification. In this study, we sequenced and annotated the chloroplast genome of Artemisia littoricola sourced from Dokdo Island, employing comparative analyses across six diverse Artemisia species. Our findings reveal conserved genome structures with variations in repeat sequences and junction boundaries. Notably, the chloroplast genome of A. littoricola spans 150,985 bp, consistent with other Artemisia species, and comprises 131 genes, including 86 protein-coding, 37 tRNA, and 8 rRNA genes. Among these genes, 16 possess a single intron, while clpP and ycf3 exhibit two introns each. Furthermore, 18 genes display duplicated copies within the IR regions. Moreover, the genome possesses 42 Simple Sequence Repeats (SSRs), predominantly abundant in A/T content and located within intergenic spacer regions. The analysis of codon usage revealed that the codons for leucine were the most frequent, with a preference for ending with A/U. While the chloroplast genome exhibited conservation overall, non-coding regions showed lower conservation compared to coding regions, with the Inverted Repeat (IR) region displaying higher conservation than single-copy regions. Phylogenetic analyses position A. littoricola within subgenus Dracunculus, indicating a close relationship with A. scoparia and A. desertorum. Additionally, biogeographic reconstructions suggest ancestral origins in East Asia, emphasizing Mongolia, China (North East and North Central and South Central China), and Korea. This study underscores the importance of chloroplast genomics in understanding Artemisia diversity and evolution, offering valuable insights into taxonomy, evolutionary patterns, and biogeographic history. These findings not only enhance our understanding of Artemisia's intricate biology but also contribute to conservation efforts and facilitate the development of molecular markers for further research and applications in medicine and agriculture.

Comparative chloroplast genome analyses of cultivated and wild Capsicum species shed light on evolution and phylogeny

The chloroplast (cp.) genome, also known as plastome, plays crucial roles in plant survival, adaptation, and evolution. The stable genetic structure of cp. genomes provides an ideal system for investigating species evolution. We sequenced three complete cp. genome sequences of Capsicum species and analyzed them using sequences of various Capsicum species retrieved from the NCBI database. The cp. genome of Capsicum species maintains a well-preserved quadripartite structure consisting of two inverted repeats (IRs) flanked by a large single copy (LSC) region and a small single copy (SSC) region. The sizes of cp. genome sequences ranged from 156,583 bp (C. lycianthoides) to 157,390 bp (C.pubescens). A total of 127-132 unique genes, including 83-87 protein-coding, 36-37 tRNA, and eight rRNA genes, were predicted. Comparison of cp. genomes of 10 Capsicum species revealed high sequence similarity in genome-wide organization and gene arrangements. Fragments of trnT-UGU/trnL-UAA, ccsA, ndhD, rps12, and ycf1 were identified as variable regions, and nucleotide variability of LSC and SSC was higher than that of IR. Phylogenetic speciation analysis showed that the major domesticated C. annuum species were the most extensively divergent species and closely related to C. tovarii and C. frutescens. Analysis of divergent times suggested that a substantial range of speciation events started occurring ~ 25.79 million years ago (Mya). Overall, comparative analysis of cp. genomes of Capsicum species not only offers new insights into their genetic variation and phylogenetic relationships, but also lays a foundation for evolutionary history, genetic diversity, conservation, and biological breeding of Capsicum species.

The chloroplast genome inheritance pattern of the Deli-Nigerian prospection material (NPM) × Yangambi population of Elaeis guineensis Jacq

Background

The chloroplast genome has the potential to be genetically engineered to enhance the agronomic value of major crops. As a crop plant with major economic value, it is important to understand every aspect of the genetic inheritance pattern among Elaeis guineensis individuals to ensure the traceability of agronomic traits.

Methods

Two parental E. guineensis individuals and 23 of their F1 progenies were collected and sequenced using the next-generation sequencing (NGS) technique on the Illumina platform. Chloroplast genomes were assembled de novo from the cleaned raw reads and aligned to check for variations. The sequences were compared and analyzed with programming language scripting and relevant bioinformatic softwares. Simple sequence repeat (SSR) loci were determined from the chloroplast genome.

Results

The chloroplast genome assembly resulted in 156,983 bp, 156,988 bp, 156,982 bp, and 156,984 bp. The gene content and arrangements were consistent with the reference genome published in the GenBank database. Seventy-eight SSRs were detected in the chloroplast genome, with most located in the intergenic spacer region.The chloroplast genomes of 17 F1 progenies were exact copies of the maternal parent, while six individuals showed a single variation in the sequence. Despite the significant variation displayed by the male parent, all the nucleotide variations were synonymous. This study show highly conserve gene content and sequence in Elaeis guineensis chloroplast genomes. Maternal inheritance of chloroplast genome among F1 progenies are robust with a low possibility of mutations over generations. The findings in this study can enlighten inheritance pattern of Elaeis guineensis chloroplast genome especially among crops' scientists who consider using chloroplast genome for agronomic trait modifications.

Characterization of Firmiana danxiaensis plastomes and comparative analysis of Firmiana: insight into its phylogeny and evolution

BACKGROUND

Firmiana danxiaensis is a critically endangered and ecologically important tree currently only found in four locations in Danxia or Karst habitats in northern Guangdong Province, China. The specialized habitat preference makes it an ideal model species for study of adaptive evolution. Meanwhile, the phylogenetic relationships of F. danxiaensis in four locations under two landforms are unclear. Therefore, we sequenced its complete chloroplast (cp.) genomes and conducted comprehensive interspecific and intrageneric plastome studies.

RESULTS

The F. danxiaensis plastomes in four locations showed a typical quadripartite and circular structure that ranged from 160,832 to 161,206 bp in size, with 112 unique genes encoded. Comparative genomics showed that the plastomes of F. danxiaensis were relatively conserved with high similarity of genome organization, gene number, GC content and SSRs. While the genomes revealed higher biased codon preferences in Karst habitat than those in Danxia habitats. Eighteen and 11 divergent hotpots were identified at interspecific and intrageneric levels for species identification and further phylogenetic studies. Seven genes (clpP, accD, ccsA, ndhH, rpl20, rpoC2, and rps4) were under positive selection and may be related to adaptation. Phylogenetic analysis revealed that F. danxiaensis is sister to F. major and F. simplex. However, the interspecific relationships are not consistent with the habitat types.

CONCLUSIONS

The characteristics and interspecific relationship of F. danxiaensis plastomes provide new insights into further integration of geographical factors, environmental factors, and genetic variations on the genomic study of F. danxiaensis. Together, our study will contribute to the study of species identification, population genetics, and conservation biology of F. danxiaensis.

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

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

The complete chloroplast genome of Diplodiscus trichospermus and phylogenetic position of Brownlowioideae within Malvaceae

BACKGROUND

Malvaceae is an economically important plant family of 4,225 species in nine subfamilies. Phylogenetic relationships among the nine subfamilies have always been controversial, especially for Brownlowioideae, whose phylogenetic position remains largely unknown due to the lack of samples in previous analysis datasets. To greatly clarify the phylogenetic relationship of Malvaceae, we newly sequenced and assembled the plastome of Diplodiscus trichospermus taxonomically located in Brownlowioideae, and downloaded the allied genomes from public database to build a dataset covering all subfamily members of Malvaceae.

RESULTS

The annotation results showed that the plastome of Diplodiscus trichospermus has a typical quadripartite structure, comprising 112 unique genes, namely 78 protein-coding genes, 30 tRNA genes and 4 rRNA genes. The total length was 158,570 bp with 37.2% GC content. Based on the maximum likelihood method and Bayesian inference, a robust phylogenetic backbone of Malvaceae was reconstructed. The topology showed that Malvaceae was divided distinctly into two major branches which were previously recognized as Byttneriina and Malvadendrina. In the Malvadendrina clade, Malvoideae and Bombacoideae formed, as always, a close sister clade named as Malvatheca. Subfamily Helicteroideae occupied the most basal position and was followed by Sterculioideae which was sister to the alliance of Malvatheca, Brownlowioideae, Dombeyoideae, and Tilioideae. Brownlowioideae together with the clade comprising Dombeyoideae and Tilioideae formed a sister clade to Malvatheca. In addition, one specific conservation SSR and three specific palindrome sequences were observed in Brownlowioideae.

CONCLUSIONS

In this study, the phylogenetic framework of subfamilies in Malvaceae has been resolved clearly based on plastomes, which may contribute to a better understanding of the classification and plastome evolution for Malvaceae.

Characterization and Comparison of the Two Mitochondrial Genomes in the Genus Rana

The mitochondrial genome (mitogenome) possesses several invaluable attributes, including limited recombination, maternal inheritance, a fast evolutionary rate, compact size, and relatively conserved gene arrangement, all of which make it particularly useful for applications in phylogenetic reconstruction, population genetics, and evolutionary research. In this study, we aimed to determine the complete mitogenomes of two morphologically similar Rana species (Rana hanluica and Rana longicrus) using next-generation sequencing. The entire circular mitogenome was successfully identified, with a length of 19,395 bp for R. hanluica and 17,833 bp for R. longicrus. The mitogenomes of both species contained 37 genes, including 13 protein-coding genes (PCGs), two ribosomal RNA genes, 22 transfer RNA genes, and one control region; mitogenome size varied predominantly with the length of the control region. The two synonymous codon usages in 13 PCGs showed that T and A were used more frequently than G and C. The ratios of non-synonymous to synonymous substitutions of all 13 PCGs were <1 in the Rana species, indicating that the PCGs were under purifying selection. Finally, phylogenetic relationship analyses suggested that R. hanluica and R. longicrus were classified in the R. japonica group. Our study provides valuable reference material for the taxonomy of the genus Rana.

Complete chloroplast genomes of three wild perennial Hordeum species from Central Asia: genome structure, mutation hotspot, phylogenetic relationships, and comparative analysis

Hordeum L. is widely distributed in mountain or plateau of subtropical and warm temperate regions around the world. Three wild perennial Hordeum species, including H. bogdanii, H. brevisubulatum, and H. violaceum, have been used as forage and for grassland ecological restoration in high-altitude areas in recent years. To date, the degree of interspecies sequence variation in the three Hordeum species within existing gene pools is still not well-defined. Herein, we sequenced and assembled chloroplast (cp) genomes of the three species. The results revealed that the cp genome of H. bogdanii showed certain sequence variations compared with the cp genomes of the other two species (H. brevisubulatum and H. violaceum), and the latter two were characterized by a higher relative affinity. Parity rule 2 plot (PR2) analysis illuminated that most genes of all ten Hordeum species were concentrated in nucleotide T and G. Numerous single nucleotide polymorphism (SNP) and insertion/deletion (In/Del) events were detected in the three Hordeum species. A series of hotspots regions (tRNA-GGU ~ tRNA-GCA, tRNA-UGU ~ ndhJ, psbE ~ rps18, ndhF ~ tRNA-UAG, etc.) were identified by mVISTA procedures, and the five highly polymorphic genes (tRNA-UGC, tRNA-UAA, tRNA-UUU, tRNA-UAC, and ndhA) were proved by the nucleotide diversity (Pi). Although the distribution and existence of cp simple sequence repeats (cpSSRs) were predicted in the three Hordeum cp genomes, no rearrangement was found between them. A similar phenomenon has been found in the cp genome of the other seven Hordeum species, which has been published so far. In addition, evolutionary relationships were reappraised based on the currently reported cp genome of Hordeum L. This study offers a framework for gaining a better understanding of the evolutionary history of Hordeum species through the re-examination of their cp genomes, and by identifying highly polymorphic genes and hotspot regions that could provide important insights into the genetic diversity and differentiation of these species.

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

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

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

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

Comparative plastome genomics, taxonomic delimitation and evolutionary divergences of Tetraena hamiensis var. qatarensis and Tetraena simplex (Zygophyllaceae)

The Zygophyllum and Tetraena genera are intriguingly important ecologically and medicinally. Based on morphological characteristics, T. hamiensis var. qatarensis, and T. simplex were transferred from Zygophyllum to Tetraena with the least genomic datasets available. Hence, we sequenced the T. hamiensis and T. simplex and performed in-depth comparative genomics, phylogenetic analysis, and estimated time divergences. The complete plastomes ranged between 106,720 and 106,446 bp-typically smaller than angiosperms plastomes. The plastome circular genomes are divided into large single-copy regions (~ 80,964 bp), small single-copy regions (~ 17,416 bp), and two inverted repeats regions (~ 4170 bp) in both Tetraena species. An unusual shrinkage of IR regions 16-24 kb was identified. This resulted in the loss of 16 genes, including 11 ndh genes which encode the NADH dehydrogenase subunits, and a significant size reduction of Tetraena plastomes compared to other angiosperms. The inter-species variations and similarities were identified using genome-wide comparisons. Phylogenetic trees generated by analyzing the whole plastomes, protein-coding genes, matK, rbcL, and cssA genes exhibited identical topologies, indicating that both species are sisters to the genus Tetraena and may not belong to Zygophyllum. Similarly, based on the entire plastome and proteins coding genes datasets, the time divergence of Zygophyllum and Tetraena was 36.6 Ma and 34.4 Ma, respectively. Tetraena stem ages were 31.7 and 18.2 Ma based on full plastome and protein-coding genes. The current study presents the plastome as a distinguishing and identification feature among the closely related Tetraena and Zygophyllum species. It can be potentially used as a universal super-barcode for identifying plants.

Identification of Laportea bulbifera using the complete chloroplast genome as a potentially effective super-barcode

Laportea bulbifera, a Miao medicine grown in karst areas, has exerted a unique curative effect on skin itching in the elderly, with an annual sales of > 100 million Yuan. Owing to the shortage of resources and large morphological variations in L. bulbifera, it is difficult to identify the species correctly using only traditional methods, which seriously affects the safety of drug usage for patients. This study obtained the complete high-quality L. bulbifera chloroplast (cp) genome, using second- and third-generation high-throughput sequencing. The cp genome was 149,911 bp in length, with a typical quadripartite structure. A total of 127 genes were annotated, including 83 protein-coding genes, 36 tRNA genes, and 8 rRNA genes. There was an inverted small single copy (SSC) structure in the L. bulbifera cp genome, one large-scale rearrangement of ~ 39 kb excised in the SSC and IR regions. The complete cp genome sequence is used as a potentially effective super-barcode and the highly variable regions (ycf1, matK, and ndhD) can be used as potentially specific barcodes to accurately distinguish L. bulbifera from counterfeits and closely related species. This study is important for the identification of L. bulbifera and lays a theoretical foundation for elucidating the phylogenetic relationship of the species.

The Complete Chloroplast Genome Sequence of Machilus chuanchienensis (Lauraceae): Genome Structure and Phylogenetic Analysis

Machilus chuanchienensis is an ecological tree distributed in southwestern China. It has a significant valuation with making Hawk tea using its leaves, an ethnic traditional tea-like beverage with a long history in Chinese tea culture. The whole chloroplast (cp) genome is an ideal model for the phylogenetic study of Lauraceae because of its simple structure and highly conserved features. There have been numerous reports of complete cp genome sequences in Lauraceae, but little is known about M. chuanchienensis. Here, the next-generation sequencing (NGS) was used to sequence the M. chuanchienensis cp genome. Then, a comprehensive comparative genome analysis was performed. The results revealed that the M. chuanchienensis's cp genome measured 152,748 base pairs (bp) with a GC content of 39.15% and coded 126 genes annotated, including comprising eight ribosomal RNA (rRNA), 36 transporter RNA (tRNA), and 82 protein-coding genes. In addition, the cp genome presented a typical quadripartite structure comprising a large single-copy (LSC; 93,811) region, a small single-copy (SSC; 18,803) region, and the inverted repeats (IRs; 20,067) region and contained 92 simple sequence repeat (SSR) locus in total. Phylogenetic relationships of 37 species indicated that M. chuanchienensis was a sister to M. balansae, M. melanophylla, and M. minutiflora. Further research on this crucial species may benefit significantly from these findings.

Dynamic evolution of the plastome in the Elm family (Ulmaceae)

This study compared the plastomes of Ulmaceae allowing analyses of the dynamic evolution, including genome structure, codon usage bias, repeat sequences, molecular mutation rates, and phylogenetic inferences. Ulmaceae is a small family in the order Rosales. This family consists of seven genera, including Ulmus, Zelkova, Planera, Hemiptelea, Phyllostylon, Ampelocera, and Holoptelea. Ulmaceae is an interesting lineage from plant biogeographic, systematic, evolutionary, and paleobotanic perspectives. It is also a good model to investigate the evolution of the plastomes in woody plants. In this study, we sequenced and assembled the complete plastomes of the six Ulmaceae genera to compare genomic structures and reveal the molecular evolutionary patterns. The size of the quadripartite plastomes ranged from 158,290 bp to 161,886 bp. The genomes contained 131 genes, including 87 coding genes, 36 tRNA, and 8 rRNA. The gene number, gene content, and genomic structure were highly consistent among the Ulmaceae genera. Nine variable regions including ndhA intron, ndhF-rpl32, ycf1, psbK-trnS, rps16-trnQ, trnT-trnL, trnT-psbD, trnS-trnG, and rpl32-trnL, were identified in Ulmaceae plastomes according to the nucleotide diversity values. Condon usage was biased among the genes and showed consistent trends in the seven genera. Molecular evolution analyses revealed that most of the genes and all gene groups were under widespread purifying selection. Twelve genes (ccsA, matK, psbH, psbK, rbcL, rpl22, rpl32, rpoA, rps12, rps15, rps16, and ycf2) were under positive selection. Phylogenetic analyses supported that Ulmaceae should be divided into two main clades, such as the temperate clade, including Ulmus, Zelkova, Planera, and Hemiptelea and the tropical clade, including Phyllostylon, Ampelocera and Holoptelea. This study reports the structure and evolutionary characteristics of the Elm family. These new genomic data will benefit assessments of genomic evolution and provide information to elucidate the phylogenetic relationships among Ulmaceae species.

Phylogenomics, plastome structure and species identification in Mahonia (Berberidaceae)

BACKGROUND

Elucidating the phylogenetic relationships within species-rich genera is essential but challenging, especially when lineages are assumed to have been going through radiation events. Mahonia Nutt. (Berberidaceae) is a genus with cosmopolitan distribution, comprising approximately 100 species, two of which are known as Caulis Mahoniae (M. bealei and M. fortunei) with crucial pharmacological significance in Chinese herbal medicine. Mahonia is a taxonomically challenging genus, and intrageneric phylogenetic relationships still need to be explored using genome data. Universal DNA barcodes and floral morphological attributes have limited discriminatory power in Mahonia.

RESULTS

We sequenced 17 representative plastomes and integrated three published plastome data together to conduct comparative and phylogenetic analyses. We found that Mahonia and Berberis share a large IR expansion (~ 12 kb), which is recognized as a typical character of Berberideae. Repeated sequences are revealed in the species of Mahonia, which are valuable for further population genetic studies. Using a comparative plastome analysis, we determined eight hypervariable regions whose discriminative power is comparable to that of the whole plastid genomes. The incongruence of the ITS and the plastome tree topologies may be ascribed to ancestral hybridization events and/or to incomplete lineage sorting. In addition, we suggest that leaf epidermal characters could help to distinguish closely related species in Mahonia.

CONCLUSIONS

We propose an integrative approach combining special barcodes and micromorphological traits to circumscribe Mahonia species. The results cast a new light on the development of an integrative method for accurate species circumscription and provide abundant genetic resources for further research on Mahonia.

Chloroplast Genomes of Genus Tilia: Comparative Genomics and Molecular Evolution

Tilia is a complex genus in the family Malvaceae that has high ecological and economical values. Owing to the lack of sufficient distinguishable morphological and molecular characteristics, interspecific relationships in this genus are not clear. Chloroplast (cp) genomes are small, meanwhile most angiosperms usually undergo matrilineal inheritance. Consequently, they can be used in molecular systematics and phylogenetic analyses. Here, we sequenced and assembled cp genomes from T. endochrysea, T. tomentosa, T. miqueliana, T. americana and T. cordata, and compared them with those of seven previously reported Tilia species. Similar gene contents, gene orders and GC contents existed among the 12 cp genomes, which ranged from 162,564 to 162,855 bp and encoded 113 unique genes. Abundant simple sequence repeats (119–127) and dispersed repeats (97–135) were detected in Tilia cp genomes. In total, 11 hypervariable regions were identified that could be suitable for species identification and phylogenetic studies. A phylogenetic analysis of Malvaceae based on 5 hypervariable genes (matK + ndhF + rpoB + rpoC2+ycf1) revealed that all eight subfamilies were monophyletic groups. Additionally, the genus Tilia was divided into three groups on the basis of all 521 molecular variation loci. The current study provides valuable insights into the genomic evolution of the genus Tilia.

Comparative Chloroplast Genome Analysis of Wax Gourd (Benincasa hispida) with Three Benincaseae Species, Revealing Evolutionary Dynamic Patterns and Phylogenetic Implications

Benincasa hispida (wax gourd) is an important Cucurbitaceae crop, with enormous economic and medicinal importance. Here, we report the de novo assembly and annotation of the complete chloroplast genome of wax gourd with 156,758 bp in total. The quadripartite structure of the chloroplast genome comprises a large single-copy (LSC) region with 86,538 bp and a small single-copy (SSC) region with 18,060 bp, separated by a pair of inverted repeats (IRa and IRb) with 26,080 bp each. Comparison analyses among B. hispida and three other species from Benincaseae presented a significant conversion regarding nucleotide content, genome structure, codon usage, synonymous and non-synonymous substitutions, putative RNA editing sites, microsatellites, and oligonucleotide repeats. The LSC and SSC regions were found to be much more varied than the IR regions through a divergent analysis of the species within Benincaseae. Notable IR contractions and expansions were observed, suggesting a difference in genome size, gene duplication and deletion, and the presence of pseudogenes. Intronic gene sequences, such as trnR-UCU–atpA and atpH–atpI, were observed as highly divergent regions. Two types of phylogenetic analysis based on the complete cp genome and 72 genes suggested sister relationships between B. hispida with the Citrullus, Lagenaria, and Cucumis. Variations and consistency with previous studies regarding phylogenetic relationships are discussed. The cp genome of B. hispida provides valuable genetic information for the detection of molecular markers, research on taxonomic discrepancies, and the inference of the phylogenetic relationships of Cucurbitaceae.

The complete chloroplast genome of Nymphaea thermarum (Nymphaeaceae) from Rwanda, Africa

Nymphaea thermarum is classified in the Nymphaeaceae, and is the smallest water lily in the world. It has been extinct its native environment and needs urgent protection. Here, we report and characterize the complete chloroplast genome of N. thermarum. The total length of the chloroplast genome is 159,849 bp and the GC content is 39.2% (A: 30.1%, C: 20.0%, G: 19.2%, T: 30.8%). The chloroplast genome consists of 8 rRNA, 37 tRNA, and 85 protein-coding genes. Phylogenetic analysis of N. thermarum fully resolved this taxon in a clade with Nymphaea capensis. The chloroplast genome of N. thermarum provides scientific guidance for its conservation genetics and also contributes genome resources for the phylogenetic relationship of Nymphaea.

Comparative Analysis of Chloroplast Genome and New Insights Into Phylogenetic Relationships of Polygonatum and Tribe Polygonateae

Members of Polygonatum are perennial herbs that have been widely used in traditional Chinese medicine to invigorate Qi, moisten the lung, and benefit the kidney and spleen among patients. However, the phylogenetic relationships and intrageneric taxonomy within Polygonatum have long been controversial because of the complexity of their morphological variations and lack of high-resolution molecular markers. The chloroplast (cp) genome is an optimal model for deciphering phylogenetic relationships in related families. In the present study, the complete cp genome of 26 species of Trib. Polygonateae were de novo assembled and characterized; all species exhibited a conserved quadripartite structure, that is, two inverted repeats (IR) containing most of the ribosomal RNA genes, and two unique regions, large single sequence (LSC) and small single sequence (SSC). A total of 8 highly variable regions (rps16-trnQ-UUG, trnS-GCU-trnG-UCC, rpl32-trnL-UAG, matK-rps16, petA-psbJ, trnT-UGU-trnL-UAA, accD-psaI, and trnC-GCA-petN) that might be useful as potential molecular markers for identifying Polygonatum species were identified. The molecular clock analysis results showed that the divergence time of Polygonatum might occur at ∼14.71 Ma, and the verticillate leaf might be the ancestral state of this genus. Moreover, phylogenetic analysis based on 88 cp genomes strongly supported the monophyly of Polygonatum. The phylogenetic analysis also suggested that Heteropolygonatum may be the sister group of the Polygonatum, but the Disporopsis, Maianthemum, and Disporum may have diverged earlier. This study provides valuable information for further species identification, evolution, and phylogenetic research of Polygonatum.

Comparative Chloroplast Genome Analyses of the Winter-Blooming Eastern Asian Endemic Genus Chimonanthus (Calycanthaceae) With Implications For Its Phylogeny and Diversification

Chimonanthus of Calycanthaceae is a small endemic genus in China, with unusual winter-blooming sweet flowers widely cultivated for ornamentals and medicinal uses. The evolution of Chimonanthus plastomes and its phylogenetic relationships remain unresolved due to limited availability of genetic resources. Here, we report fully assembled and annotated chloroplast genomes of five Chimonanthus species. The chloroplast genomes of the genus (size range 153,010 – 153,299 bp) reveal high similarities in gene content, gene order, GC content, codon usage, amino acid frequency, simple sequence repeats, oligonucleotide repeats, synonymous and non-synonymous substitutions, and transition and transversion substitutions. Signatures of positive selection are detected in atpF and rpoB genes in C. campanulatus. The correlations among substitutions, InDels, and oligonucleotide repeats reveal weak to strong correlations in distantly related species at the intergeneric levels, and very weak to weak correlations among closely related Chimonanthus species. Chloroplast genomes are used to reconstruct a well-resolved phylogenetic tree, which supports the monophyly of Chimonanthus. Within Chimonanthus, C. praecox and C. campanulatus form one clade, while C. grammatus, C. salicifolius, C. zhejiangensis, and C. nitens constitute another clade. Chimonanthus nitens appears paraphyletic and is closely related to C. salicifolius and C. zhejiangensis, suggesting the need to reevaluate the species delimitation of C. nitens. Chimonanthus and Calycanthus diverged in mid-Oligocene; the radiation of extant Chimonanthus species was dated to the mid-Miocene, while C. grammatus diverged from other Chimonanthus species in the late Miocene. C. salicifolius, C. nitens(a), and C. zhejiangensis are inferred to have diverged in the Pleistocene of the Quaternary period, suggesting recent speciation of a relict lineage in the subtropical forest regions in eastern China. This study provides important insights into the chloroplast genome features and evolutionary history of Chimonanthus and family Calycanthaceae.

Chloroplast Phylogenomic Analyses Reveal a Maternal Hybridization Event Leading to the Formation of Cultivated Peanuts

Peanuts (Arachis hypogaea L.) offer numerous healthy benefits, and the production of peanuts has a prominent role in global food security. As a result, it is in the interest of society to improve the productivity and quality of peanuts with transgenic means. However, the lack of a robust phylogeny of cultivated and wild peanut species has limited the utilization of genetic resources in peanut molecular breeding. In this study, a total of 33 complete peanut plastomes were sequenced, analyzed and used for phylogenetic analyses. Our results suggest that sect. Arachis can be subdivided into two lineages. All the cultivated species are contained in Lineage I with AABB and AA are the two predominant genome types present, while species in Lineage II possess diverse genome types, including BB, KK, GG, etc. Phylogenetic studies also indicate that all allotetraploid cultivated peanut species have been derived from a possible maternal hybridization event with one of the diploid Arachis duranensis accessions being a potential AA sub-genome ancestor. In addition, Arachis monticola, a tetraploid wild species, is placed in the same group with all the cultivated peanuts, and it may represent a transitional species, which has been through the recent hybridization event. This research could facilitate a better understanding of the taxonomic status of various Arachis species/accessions and the evolutionary relationship among them, and assists in the correct and efficient use of germplasm resources in breeding efforts to improve peanuts for the benefit of human beings.

Fossil fruits of Firmiana and Tilia from the middle Miocene of South Korea and the efficacy of the Bering land bridge for the migration of mesothermal plants

Determining whether the high-latitude Bering land bridge (BLB) was ecologically suitable for the migration of mesothermal plants is significant for Holarctic phytogeographic inferences. Paleobotanical studies provide a critical source of data on the latitudinal positions of different plant lineages at different times, permitting assessment of the efficacy of the BLB for migration. Here we report exceptionally preserved fossils of Firmiana and Tilia endochrysea from the middle Miocene of South Korea. This represents a new reliable record of Firmiana and the first discovery of the T. endochrysea lineage in the fossil record of Asia. The occurrence of these fossils in South Korea indicates that the two lineages had a distribution that extended much farther north during the middle Miocene, but they were still geographically remote from the BLB. In light of the broader fossil record of Asia, our study shows that, in the middle Miocene, some mesothermal plants apparently inhabited the territory adjacent to the BLB and thus they were possibly capable of utilizing the BLB as a migratory corridor. Some other mesothermal plants, such as Firmiana and the T. endochrysea lineages, however, are restricted to more southern regions relative to the BLB based on current fossil evidence. These lineages may have been ecologically unable to traverse the BLB, which raises questions about the efficacy of the BLB as a universal exchange route for mesothermal plants between Asia and North America during the middle Miocene.

The BAHD Gene Family in Cacao (Theobroma cacao, Malvaceae): Genome-Wide Identification and Expression Analysis

The benzyl alcohol O-acetyl transferase, anthocyanin O-hydroxycinnamoyl transferase, N-hydroxycinnamoyl anthranilate benzoyl transferase, and deacetylvindoline 4-O-acetyltransferase (BAHD) enzymes play a critical role in regulating plant metabolites and affecting cell stability. In the present study, members of the BAHD gene family were recognized in the genome of Theobroma cacao and characterized using various bioinformatics tools. We found 27 non-redundant putative tcBAHD genes in cacao for the first time. Our findings indicate that tcBAHD genes are diverse based on sequence structure, physiochemical properties, and function. When analyzed with BAHDs of Gossypium raimondii and Corchorus capsularis clustered into four main groups. According to phylogenetic analysis, BAHD genes probably evolved drastically after their divergence. The divergence time of duplication events with purifying selection pressure was predicted to range from 1.82 to 15.50 MYA. Pocket analysis revealed that serine amino acid is more common in the binding site than other residuals, reflecting its key role in regulating the activity of tcBAHDs. Furthermore, cis-acting elements related to the responsiveness of stress and hormone, particularly ABA and MeJA, were frequently observed in the promoter region of tcBAHD genes. RNA-seq analysis further illustrated that tcBAHD13 and tcBAHD26 are involved in response to Phytophthora megakarya fungi. In conclusion, it is likely that evolutionary processes, such as duplication events, have caused high diversity in the structure and function of tcBAHD genes.

Pseudogenization of the chloroplast threonine (trnT-GGU) gene in the sunflower family (Asteraceae)

The chloroplast genome evolves through the course of evolution. Various types of mutational events are found within the chloroplast genome, including insertions-deletions (InDels), substitutions, inversions, gene rearrangement, and pseudogenization of genes. The pseudogenization of the chloroplast threonine (trnT-GGU) gene was previously reported in Cryptomeria japonica (Cupressaceae), Pelargonium × hortorum (Geraniaceae), and Anaphalis sinica and Leontopodium leiolepis of the tribe Gnaphalieae (Asteroideae, Asteraceae). Here, we performed a broad analysis of the trnT-GGU gene among the species of 13 subfamilies of Asteraceae and found this gene as a pseudogene in core Asteraceae (Gymnarrhenoideae, Cichorioideae, Corymbioideae, and Asteroideae), which was linked to an insertion event within the 5' acceptor stem and is not associated with ecological factors such as habit, habitat, and geographical distribution of the species. The pseudogenization of trnT-GGU was not predicted in codon usage, indicating that the superwobbling phenomenon occurs in core Asteraceae in which a single transfer RNA (trnT-UGU) decodes all four codons of threonine. To the best of our knowledge, this is the first evidence of a complete clade of a plant species using the superwobbling phenomenon for translation.

On the Species Delimitation of the Maddenia Group of Prunus (Rosaceae): Evidence From Plastome and Nuclear Sequences and Morphology

The recognition, identification, and differentiation of closely related plant species present significant and notorious challenges to taxonomists. The Maddenia group of Prunus, which comprises four to seven species, is an example of a group in which species delimitation and phylogenetic reconstruction have been difficult, due to the lack of clear morphological distinctions, limited sampling, and low informativeness of molecular evidence. Thus, the precise number of species in the group and the relationships among them remain unclear. Here, we used genome skimming to generate the DNA sequence data for 22 samples, including 17 Maddenia individuals and five outgroups in Amygdaloideae of Rosaceae, from which we assembled the plastome and 446 single-copy nuclear (SCN) genes for each sample. The phylogenetic relationships of the Maddenia group were then reconstructed using both concatenated and coalescent-based methods. We also identified eight highly variable regions and detected simple sequence repeats (SSRs) and repeat sequences in the Maddenia species plastomes. The phylogenetic analysis based on the complete plastomes strongly supported three main subclades in the Maddenia group of Prunus, while five subclades were recognized based on the nuclear tree. The phylogenetic network analysis detected six hybridization events. Integrating the nuclear and morphological evidence, we proposed to recognize five species within the Maddenia group, i.e., Prunus fujianensis, P. himalayana, P. gongshanensis, P. hypoleuca, and P. hypoxantha. Within this group, the first three species are well-supported, while the gene flow occurring throughout the Maddenia group seems to be especially frequent between P. hypoleuca and P. hypoxantha, eroding the barrier between them. The phylogenetic trees based on eight concatenated hypervariable regions had a similar topology with the complete plastomes, showing their potential as molecular markers and effective barcodes for further phylogeographic studies on Maddenia.

The complete chloroplast genome of Wurfbainia neoaurantiaca (Zingiberaceae: Zingiberea) from Yunnan Province, China

Wurfbainia neoaurantiaca is a medicinal plant endemic to Yunnan Province, China. In this study, its complete chloroplast genome was assembled and characterized. The total genome size of W. neoaurantiaca was 158,484 bp in length, consisting of a large single-copy region (LSC), a small single-copy region (SSC) and two inverted repeat regions (IRs) with 88,605 bp, 15,285 bp and 29,822 bp, respectively. Its GC content was 36.08%. The chloroplast genome encoded 113 unique genes, including 79 protein-coding, 30 tRNA, and four rRNA genes. The result of the phylogenetic analysis indicated that W. neoaurantiaca was related to W. villosa var. xanthioides and supported de Boer’s classification that W. compacta, W. longiligularis, W. neoaurantiaca, W. villosa, W. villosa var. xanthioides and Amomum krervanh belonged to the Wurfbainia Clade.

Comparative plastome analysis of Blumea, with implications for genome evolution and phylogeny of Asteroideae

The genus Blumea (Asteroideae, Asteraceae) comprises about 100 species, including herbs, shrubs, and small trees. Previous studies have been unable to resolve taxonomic issues and the phylogeny of the genus Blumea due to the low polymorphism of molecular markers. Therefore, suitable polymorphic regions need to be identified. Here, we de novo assembled plastomes of the three Blumea species B. oxyodonta, B. tenella, and B. balsamifera and compared them with 26 other species of Asteroideae after correction of annotations. These species have quadripartite plastomes with similar gene content, genome organization, and inverted repeat contraction and expansion comprising 113 genes, including 80 protein-coding, 29 transfer RNA, and 4 ribosomal RNA genes. The comparative analysis of codon usage, amino acid frequency, microsatellite repeats, oligonucleotide repeats, and transition and transversion substitutions has revealed high resemblance among the newly assembled species of Blumea. We identified 10 highly polymorphic regions with nucleotide diversity above 0.02, including rps16-trnQ, ycf1, ndhF-rpl32, petN-psbM, and rpl32-trnL, and they may be suitable for the development of robust, authentic, and cost-effective markers for barcoding and inference of the phylogeny of the genus Blumea. Among these highly polymorphic regions, five regions also co-occurred with oligonucleotide repeats and support use of repeats as a proxy for the identification of polymorphic loci. The phylogenetic analysis revealed a close relationship between Blumea and Pluchea within the tribe Inuleae. At tribe level, our phylogeny supports a sister relationship between Astereae and Anthemideae rooted as Gnaphalieae, Calenduleae, and Senecioneae. These results are contradictory to recent studies which reported a sister relationship between "Senecioneae and Anthemideae" and "Astereae and Gnaphalieae" or a sister relationship between Astereae and Gnaphalieae rooted as Calenduleae, Anthemideae, and then Senecioneae using nuclear genome sequences. The conflicting phylogenetic signals observed at the tribal level between plastidt and nuclear genome data require further investigation.

Chloroplast genomes of two Mediterranean Bupleurum species and the phylogenetic relationship inferred from combined analysis with East Asian species

The chloroplast genomes of Mediterranean Bupleurum species are reported for the first time. Phylogenetic analysis supports the species as a basal clade of Bupleurum with divergence time at 35.40 Ma. Bupleurum is one of the most species-rich genus with high medicinal value in Apiaceae. Although infrageneric classifications of Bupleurum have been the subject of numerous studies, it still remains controversial. Chloroplast genome information will prove essential in advancing our understanding on phylogenetic study. Here we report cp genomes of two woody Bupleurum species (Bupleurum gibraltaricum and B. fruticosum) endemic to Mediterranean. The complete cp genomes of the two species were 157,303 and 157,391 bp in size, respectively. They encoded 114 unique genes including 30 tRNA genes, 4 rRNA genes and 80 protein coding genes. Genome structure, distributions of SDRs and SSRs, gene content exhibited similarities among Bupleurum species. High variable hotspots were detected in eight intergenic spacers and four genes. Most of genes were under purifying selection with two exceptions: atpF and clpP. The phylogenetic analysis based on 80 coding genes revealed that the genus was divided into 2 distinct clades corresponding to the 2 subgenera (subg. Penninervia, subg. Bupleurum) with divergence time at the end of collision of India with Eurasia. Most species diversified mainly during the later period of uplift of Qinghai-Tibetan Plateau. The cp genomes of the two Bupleurum species can be significant complementary to insights into the cp genome characteristics of this genus. The comparative chloroplast genomes and phylogenetic analysis advances our understanding of the evolution of cp genomes and phylogeny in Bupleurum.

The GASA Gene Family in Theobroma cacao: Genome wide Identification and Expression Analysis.. [DOI: 10.1101/2021.01.27.425041]

The gibberellic acid-stimulated Arabidopsis (GASA/GAST) gene family is widely distributed in plants. The role of the GASA gene family has been reported previously in various physiological and biological processes, such as cell division, root and seed development, stem growth, and fruit ripening. These genes also provide resistance to abiotic and biotic stresses including antimicrobial, antiviral, and antifungal. Here, we report 17 tcGASA genes in Theobroma cacao L. distributed on six chromosomes. The gene structure, promoter-region sequences, protein structure, and biochemical properties, expression, and phylogenetics of all tcGASAs were analyzed. Phylogenetic analyses divided tcGASA proteins into five groups. The nine segmentally duplicating genes form four pairs and cluster together in phylogenetic tree. Purifying selection pressure was recorded on tcGASA, including duplicated genes. Several stress/hormone-responsive cis-regulatory elements were also recognized in the promoter region of tcGASAs. Differential expression analyses revealed that most of the tcGASA genes showed elevated expression in the seeds (cacao food), implying their role in seed development. The black rod disease of genus Phytophthora caused up to 20–25% loss (700,000 metric tons) in world cacao production. The role of tcGASA genes in conferring fungal resistance was also explored based on RNAseq data against Phytophthora megakarya. The differential expression of tcGASA genes was recorded between the tolerant and susceptible cultivars of cacao plants, which were inoculated with the fungus for 24h and 72h. This differential expression indicating possible role of tcGASA genes to fungal resistant in cacao. Our findings provide new insight into the function, evolution, and regulatory system of the GASA family genes in T. cacao and provide new target genes for development of fungi-resistant cacao varieties in breeding programs.

Mutational Dynamics of Aroid Chloroplast Genomes II

The co-occurrence among single nucleotide polymorphisms (SNPs), insertions-deletions (InDels), and oligonucleotide repeats has been reported in prokaryote, eukaryote, and chloroplast genomes. Correlations among SNPs, InDels, and repeats have been investigated in the plant family Araceae previously using pair-wise sequence alignments of the chloroplast genomes of two morphotypes of one species, Colocasia esculenta belonging to subfamily Aroideae (crown group), and four species from the subfamily Lemnoideae, a basal group. The family Araceae is a large family comprising 3,645 species in 144 genera, grouped into eight subfamilies. In the current study, we performed 34 comparisons using 27 species from 7 subfamilies of Araceae to determine correlation coefficients among the mutational events at the family, subfamily, and genus levels. We express strength of the correlations as: negligible or very weak (0.10–0.19), weak (0.20–0.29), moderate (0.30–0.39), strong (0.40–0.69), very strong (0.70–0.99), and perfect (1.00). We observed strong/very strong correlations in most comparisons, whereas a few comparisons showed moderate correlations. The average correlation coefficient was recorded as 0.66 between “SNPs and InDels,” 0.50 between “InDels and repeats,” and 0.42 between “SNPs and repeats.” In qualitative analyses, 95–100% of the repeats at family and sub-family level, while 36–86% of the repeats at genus level comparisons co-occurred with SNPs in the same bins. Our findings show that such correlations among mutational events exist throughout Araceae and support the hypothesis of distribution of oligonucleotide repeats as a proxy for mutational hotspots.

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 and Evolutionary Analyses on the Complete Plastomes of Five Kalanchoe Horticultural Plants

Many species of the genus Kalanchoe are important horticultural plants. They have evolved the Crassulacean acid metabolism (CAM) photosynthetic pathway to allow them to be better adapted to dry environments. Despite their importance, it is still debating whether Kalanchoe is monophyletic, and understanding the past diversification of this genus requires a tremendous amount of effort and work being devoted to the studies of morphological and molecular characters of this genus. However, molecular information, plastic sequence data, in particular, reported on Kalanchoe species is scarce, and this has posed a great challenge in trying to interpret the evolutionary history of this genus. In this study, plastomes of the five Kalanchoe species, including Kalanchoe daigremontiana, Kalanchoe delagoensis, Kalanchoe fedtschenkoi, Kalanchoe longiflora, and Kalanchoe pinnata, were sequenced and analyzed. The results indicate that the five plastomes are comparable in size, guanine-cytosine (GC) contents and the number of genes, which also demonstrate an insignificant difference in comparison with other species from the family Crassulaceae. About 224 simple sequence repeats (SSRs) and 144 long repeats were identified in the five plastomes, and most of these are distributed in the inverted repeat regions. In addition, highly divergent regions containing either single nucleotide polymorphism (SNP) or insertion or deletion (InDel) mutations are discovered, which could be potentially used for establishing phylogenetic relationships among members of the Kalanchoe genus in future studies. Furthermore, phylogenetic analyses suggest that Bryophyllum should be placed into one single genus as Kalanchoe. Further genomic analyses also reveal that several genes are undergone positive selection. Among them, 11 genes are involved in important cellular processes, such as cell survival, electron transfer, and may have played indispensable roles in the adaptive evolution of Kalanchoe to dry environments.

Chloroplast genome evolution in the Dracunculus clade (Aroideae, Araceae)

Chloroplast (cp) genomes are considered important for the study of lineage-specific molecular evolution, population genetics, and phylogenetics. Our aim here was to elucidate the molecular evolution in cp genomes of species in the Dracunculus clade (Aroideae, Araceae). We report de novo assembled cp genomes for eight species from eight genera and also retrieved cp genomes of four species from the National Center for Biotechnology Information (NCBI). The cp genomes varied in size from 162,424 bp to 176,835 bp. Large Single Copy (LSC) region ranged in size from 87,141 bp to 95,475 bp; Small Single Copy (SSC) from 14,338 bp to 23,981 bp; and Inverted Repeats (IRa and IRb) from 25,131 bp to 32,708 bp. The expansion in inverted repeats led to duplication of ycf1 genes in four species. The genera showed high similarity in gene content and yielded 113 unique genes (79 protein-coding, 4 rRNA, and 30 tRNA genes). Codon usage, amino acid frequency, RNA editing sites, microsatellites repeats, transition and transversion substitutions, and synonymous and non-synonymous substitutions were also similar across the clade. A previous study reported deletion of ycf1, accD, psbE, trnL-CAA, and trnG-GCC genes in four Amorphophallus species. Our study supports conservative structure of cp genomes in the Dracunculus clade including Amorphophallus species and does not support gene deletion mentioned above. We also report suitable polymorphic loci based on comparative analyses of Dracunculus clade species, which could be useful for phylogenetic inference. Overall, the current study broad our knowledge about the molecular evolution of chloroplast genome in aroids.

Complete chloroplast genomes of three important species, Abelmoschus moschatus, A. manihot and A. sagittifolius: Genome structures, mutational hotspots, comparative and phylogenetic analysis in Malvaceae

Abelmoschus is an economically and phylogenetically valuable genus in the family Malvaceae. Owing to coexistence of wild and cultivated form and interspecific hybridization, this genus is controversial in systematics and taxonomy and requires detailed investigation. Here, we present whole chloroplast genome sequences and annotation of three important species: A. moschatus, A. manihot and A. sagittifolius, and compared with A. esculentus published previously. These chloroplast genome sequences ranged from 163121 bp to 163453 bp in length and contained 132 genes with 87 protein-coding genes, 37 transfer RNA and 8 ribosomal RNA genes. Comparative analyses revealed that amino acid frequency and codon usage had similarity among four species, while the number of repeat sequences in A. esculentus were much lower than other three species. Six categories of simple sequence repeats (SSRs) were detected, but A. moschatus and A. manihot did not contain hexanucleotide SSRs. Single nucleotide polymorphisms (SNPs) of A/T, T/A and C/T were the largest number type, and the ratio of transition to transversion was from 0.37 to 0.55. Abelmoschus species showed relatively independent inverted-repeats (IR) boundary traits with different boundary genes compared with the other related Malvaceae species. The intergenic spacer regions had more polymorphic than protein-coding regions and intronic regions, and thirty mutational hotpots (≥200 bp) were identified in Abelmoschus, such as start-psbA, atpB-rbcL, petD-exon2-rpoA, clpP-intron1 and clpP-exon2.These mutational hotpots could be used as polymorphic markers to resolve taxonomic discrepancies and biogeographical origin in genus Abelmoschus. Moreover, phylogenetic analysis of 33 Malvaceae species indicated that they were well divided into six subfamilies, and genus Abelmoschus was a well-supported clade within genus Hibiscus.

Complete Chloroplast Genomes of Anthurium huixtlense and Pothos scandens (Pothoideae, Araceae): Unique Inverted Repeat Expansion and Contraction Affect Rate of Evolution

The subfamily Pothoideae belongs to the ecologically important plant family Araceae. Here, we report the chloroplast genomes of two species of the subfamily Pothoideae: Anthurium huixtlense (size: 163,116 bp) and Pothos scandens (size: 164,719 bp). The chloroplast genome of P. scandens showed unique contraction and expansion of inverted repeats (IRs), thereby increasing the size of the large single-copy region (LSC: 102,956 bp) and decreasing the size of the small single-copy region (SSC: 6779 bp). This led to duplication of many single-copy genes due to transfer to IR regions from the small single-copy (SSC) region, whereas some duplicate genes became single copy due to transfer to large single-copy regions. The rate of evolution of protein-coding genes was affected by the contraction and expansion of IRs; we found higher mutation rates for genes that exist in single-copy regions as compared to those in IRs. We found a 2.3-fold increase of oligonucleotide repeats in P. scandens when compared with A. huixtlense, whereas amino acid frequency and codon usage revealed similarities. The ratio of transition to transversion mutations was 2.26 in P. scandens and 2.12 in A. huixtlense. Transversion mutations mostly translated in non-synonymous substitutions. The phylogenetic inference of the limited species showed the monophyly of the Araceae subfamilies. Our study provides insight into the molecular evolution of chloroplast genomes in the subfamily Pothoideae and family Araceae.

Plastid genomics of Nicotiana (Solanaceae): insights into molecular evolution, positive selection and the origin of the maternal genome of Aztec tobacco (Nicotiana rustica)

Species of the genus Nicotiana (Solanaceae), commonly referred to as tobacco plants, are often cultivated as non-food crops and garden ornamentals. In addition to the worldwide production of tobacco leaves, they are also used as evolutionary model systems due to their complex development history tangled by polyploidy and hybridization. Here, we assembled the plastid genomes of five tobacco species: N. knightiana, N. rustica, N. paniculata, N. obtusifolia and N. glauca. De novo assembled tobacco plastid genomes had the typical quadripartite structure, consisting of a pair of inverted repeat (IR) regions (25,323-25,369 bp each) separated by a large single-copy (LSC) region (86,510-86,716 bp) and a small single-copy (SSC) region (18,441-18,555 bp). Comparative analyses of Nicotiana plastid genomes with currently available Solanaceae genome sequences showed similar GC and gene content, codon usage, simple sequence and oligonucleotide repeats, RNA editing sites, and substitutions. We identified 20 highly polymorphic regions, mostly belonging to intergenic spacer regions (IGS), which could be suitable for the development of robust and cost-effective markers for inferring the phylogeny of the genus Nicotiana and family Solanaceae. Our comparative plastid genome analysis revealed that the maternal parent of the tetraploid N. rustica was the common ancestor of N. paniculata and N. knightiana, and the later species is more closely related to N. rustica. Relaxed molecular clock analyses estimated the speciation event between N. rustica and N. knightiana appeared 0.56 Ma (HPD 0.65-0.46). Biogeographical analysis supported a south-to-north range expansion and diversification for N. rustica and related species, where N. undulata and N. paniculata evolved in North/Central Peru, while N. rustica developed in Southern Peru and separated from N. knightiana, which adapted to the Southern coastal climatic regimes. We further inspected selective pressure on protein-coding genes among tobacco species to determine if this adaptation process affected the evolution of plastid genes. These analyses indicate that four genes involved in different plastid functions, including DNA replication (rpoA) and photosynthesis (atpB, ndhD and ndhF), came under positive selective pressure as a result of specific environmental conditions. Genetic mutations in these genes might have contributed to better survival and superior adaptations during the evolutionary history of tobacco species.

Comparative analyses of 32 complete plastomes of Tef (Eragrostis tef ) accessions from Ethiopia: phylogenetic relationships and mutational hotspots

Eragrostis tef is an important cereal crop in Ethiopia with excellent storage properties, high–quality food, and the unique ability to thrive in extreme environmental conditions. However, the application of advanced molecular tools for breeding and conservation of these species is extremely limited. Therefore, developing chloroplast genome resources and high-resolution molecular markers are valuable to E. tef population and biogeographic studies. In the current study, we assembled and compared the complete plastomes of 32 E. tef accessions. The size of the plastomes ranged from 134,349 to 134,437 bp with similar GC content (∼38.3%). Genomes annotations revealed 112 individual genes, including 77 protein-coding, 31 tRNA, and 4 rRNA genes. Comparison of E. tef plastomes revealed a low degree of intraspecific sequence variations and no structural differentiations. Furthermore, we found 34 polymorphic sites (13 cpSSRs, 12 InDels, and 9 SNPs) that can be used as valuable DNA barcodes. Among them, the majority (88%) of the polymorphic sites were identified in the noncoding genomic regions. Nonsynonymous (ka) and synonymous (ks) substitution analysis showed that all PCGs were under purifying selection (ka/ks <1). The phylogenetic analyses of the whole plastomes and polymorphic region sequences were able to distinguish the accession from the southern population, indicating its potential to be used as a super-barcode. In conclusion, the newly generated plastomes and polymorphic markers developed here could be a useful genomic resource in molecular breeding, population genetics and the biogeographical study of E. tef.

Comparative Plastomics of Ashwagandha (Withania, Solanaceae) and Identification of Mutational Hotspots for Barcoding Medicinal Plants

Within the family Solanaceae, Withania is a small genus belonging to the Solanoideae subfamily. Here, we report the de novo assembled chloroplast genome sequences of W. coagulans, W. adpressa, and W. riebeckii. The length of these genomes ranged from 154,162 to 154,364 base pairs (bp). These genomes contained a pair of inverted repeats (IRa and IRb) ranging from 25,029 to 25,071 bp that were separated by a large single-copy (LSC) region of 85,635-85,765 bp and a small single-copy (SSC) region of 18,457-18,469 bp. We analyzed the structural organization, gene content and order, guanine-cytosine content, codon usage, RNA-editing sites, microsatellites, oligonucleotide and tandem repeats, and substitutions of Withania plastomes, which revealed high similarities among the species. Comparative analysis among the Withania species also highlighted 10 divergent hotspots that could potentially be used for molecular marker development, phylogenetic analysis, and species identification. Furthermore, our analyses showed that even three mutational hotspots (rps4-trnT, trnM-atpE, and rps15) were sufficient to discriminate the Withania species included in current study.

Comparison of Chloroplast Genomes among Species of Unisexual and Bisexual Clades of the Monocot Family Araceae

The chloroplast genome provides insight into the evolution of plant species. We de novo assembled and annotated chloroplast genomes of four genera representing three subfamilies of Araceae: Lasia spinosa (Lasioideae), Stylochaeton bogneri, Zamioculcas zamiifolia (Zamioculcadoideae), and Orontium aquaticum (Orontioideae), and performed comparative genomics using these chloroplast genomes. The sizes of the chloroplast genomes ranged from 163,770 bp to 169,982 bp. These genomes comprise 113 unique genes, including 79 protein-coding, 4 rRNA, and 30 tRNA genes. Among these genes, 17-18 genes are duplicated in the inverted repeat (IR) regions, comprising 6-7 protein-coding (including trans-splicing gene rps12), 4 rRNA, and 7 tRNA genes. The total number of genes ranged between 130 and 131. The infA gene was found to be a pseudogene in all four genomes reported here. These genomes exhibited high similarities in codon usage, amino acid frequency, RNA editing sites, and microsatellites. The oligonucleotide repeats and junctions JSB (IRb/SSC) and JSA (SSC/IRa) were highly variable among the genomes. The patterns of IR contraction and expansion were shown to be homoplasious, and therefore unsuitable for phylogenetic analyses. Signatures of positive selection were seen in three genes in S. bogneri, including ycf2, clpP, and rpl36. This study is a valuable addition to the evolutionary history of chloroplast genome structure in Araceae.

Molecular evolution of chloroplast genomes in Monsteroideae (Araceae)

This study provides broad insight into the chloroplast genomes of the subfamily Monsteroideae. The identified polymorphic regions may be suitable for designing unique and robust molecular markers for phylogenetic inference. Monsteroideae is the third largest subfamily (comprises 369 species) and one of the early diverging lineages of the monocot plant family Araceae. The phylogeny of this important subfamily is not well resolved at the species level due to scarcity of genomic resources and suitable molecular markers. Here, we report annotated chloroplast genome sequences of four Monsteroideae species: Spathiphyllum patulinervum, Stenospermation multiovulatum, Monstera adansonii, and Rhaphidophora amplissima. The quadripartite chloroplast genomes (size range 163,335-164,751 bp) consist of a pair of inverted repeats (25,270-25,931 bp), separating a small single copy region (21,448-22,346 bp) from a large single copy region (89,714-91,841 bp). The genomes contain 114 unique genes, including four rRNA genes, 80 protein-coding genes, and 30 tRNA genes. Gene features, amino acid frequencies, codon usage, GC contents, oligonucleotide repeats, and inverted repeats dynamics exhibit similarities among the four genomes. Higher rate of synonymous substitutions was observed as compared to non-synonymous substitutions in 76 protein-coding genes. Positive selection was observed in seven protein-coding genes, including psbK, ndhK, ndhD, rbcL, accD, rps8, and ycf2. Our included species of Araceae showed the monophyly in Monsteroideae and other subfamilies. We report 30 suitable polymorphic regions. The polymorphic regions identified here might be suitable for designing unique and robust markers for inferring the phylogeny and phylogeography among closely related species within the genus Spathiphyllum and among distantly related species within the subfamily Monsteroideae. The chloroplast genomes presented here are a valuable contribution towards understanding the molecular evolutionary dynamics in the family Araceae.

Evolutionary dynamics of chloroplast genomes in subfamily Aroideae (Araceae)

Aroideae is the largest and most diverse subfamily of the plant family Araceae. Despite its agricultural and horticultural importance, the genomic resources are sparse for this subfamily. Here, we report de novo assembled and fully annotated chloroplast genomes of 13 Aroideae species. The quadripartite chloroplast genomes (size range of 158,177-170,037 bp) are comprised of a large single copy (LSC; 75,594-94,702 bp), a small single copy (SSC; 12,903-23,981 bp) and a pair of inverted repeats (IRs; 25,266-34,840 bp). Notable gene rearrangements and IRs contraction / expansions were found for Anchomanes hookeri and Zantedeschia aethiopica. Codon usage, amino acid frequencies, oligonucleotide repeats, GC contents, and gene features revealed similarities among the 13 species. The number of oligonucleotide repeats was uncorrelated with genome size or phylogenetic position of the species. Phylogenetic analyses corroborated the monophyly of Aroideae but were unable to resolve the positions of Calla and Schismatoglottis.

Comparative Analyses of Euonymus Chloroplast Genomes: Genetic Structure, Screening for Loci With Suitable Polymorphism, Positive Selection Genes, and Phylogenetic Relationships Within Celastrineae

In this study, we assembled and annotated the chloroplast (cp) genome of the Euonymus species Euonymus fortunei, Euonymus phellomanus, and Euonymus maackii, and performed a series of analyses to investigate gene structure, GC content, sequence alignment, and nucleic acid diversity, with the objectives of identifying positive selection genes and understanding evolutionary relationships. The results indicated that the Euonymus cp genome was 156,860-157,611bp in length and exhibited a typical circular tetrad structure. Similar to the majority of angiosperm chloroplast genomes, the results yielded a large single-copy region (LSC) (85,826-86,299bp) and a small single-copy region (SSC) (18,319-18,536bp), separated by a pair of sequences (IRA and IRB; 26,341-26,700bp) with the same encoding but in opposite directions. The chloroplast genome was annotated to 130-131 genes, including 85-86 protein coding genes, 37 tRNA genes, and eight rRNA genes, with GC contents of 37.26-37.31%. The GC content was variable among regions and was highest in the inverted repeat (IR) region. The IR boundary of Euonymus happened expanding resulting that the rps19 entered into IR region and doubled completely. Such fluctuations at the border positions might be helpful in determining evolutionary relationships among Euonymus. The simple-sequence repeats (SSRs) of Euonymus species were composed primarily of single nucleotides (A)n and (T)n, and were mostly 10-12bp in length, with an obvious A/T bias. We identified several loci with suitable polymorphism with the potential use as molecular markers for inferring the phylogeny within the genus Euonymus. Signatures of positive selection were seen in rpoB protein encoding genes. Based on data from the whole chloroplast genome, common single copy genes, and the LSC, SSC, and IR regions, we constructed an evolutionary tree of Euonymus and related species, the results of which were consistent with traditional taxonomic classifications. It showed that E. fortunei sister to the Euonymus japonicus, whereby E. maackii appeared as sister to Euonymus hamiltonianus. Our study provides important genetic information to support further investigations into the phylogenetic development and adaptive evolution of Euonymus species.

Chloroplast Genome Sequence of Artemisia scoparia: Comparative Analyses and Screening of Mutational Hotspots

Artemisia L. is among the most diverse and medicinally important genera of the plant family Asteraceae. Discrepancies arise in the taxonomic classification of Artemisia due to the occurrence of multiple polyploidy events in separate lineages and its complex morphology. The discrepancies could be resolved by increasing the genomic resources. A. scoparia is one of the most medicinally important species in Artemisia. In this paper, we report the complete chloroplast genome sequence of Artemisia scoparia. The genome was 151,060 bp (base pairs), comprising a large single copy (82,834 bp) and small single copy (18,282 bp), separated by a pair of long inverted repeats (IRa and IRb: 24,972 bp each). We identified 114 unique genes, including four ribosomal RNAs, 30 transfer RNAs, and 80 protein-coding genes. We analysed the chloroplast genome features, including oligonucleotide repeats, microsatellites, amino acid frequencies, RNA editing sites, and codon usage. Transversion substitutions were twice as frequent as transition substitutions. Mutational hotspot loci included ccsA-ndhD, trnH-psbA, ndhG-ndhI, rps18-rpl20, and rps15-ycf1. These loci can be used to develop cost-effective and robust molecular markers for resolving the taxonomic discrepancies. The reconstructed phylogenetic tree supported previous findings of Artemisia as a monophyletic genus, sister to the genus Chrysanthemum, whereby A. scoparia appeared as sister to A. capillaris.

Chloroplast genome sequences of Artemisia maritima and Artemisia absinthium: Comparative analyses, mutational hotspots in genus Artemisia and phylogeny in family Asteraceae

Artemisia L. is a complex genus of medicinal importance. Publicly available chloroplast genomes of few Artemisia species are insufficient to resolve taxonomic discrepancies at species level. We report chloroplast genome sequences of two further Artemisia species: A. maritima (151,061 bp) and A. absinthium (151,193 bp). Both genomes possess typical quadripartite structure comprising of a large single copy, a small single copy and a pair of long inverted repeats. The two genomes exhibited high similarities in genome sizes, gene synteny, GC content, synonymous and non-synonymous substitutions, codon usage, amino acids frequencies, RNA editing sites, microsatellites, and oligonucleotide repeats. Transition to transversion ratio was <1. Maximum likelihood tree showed Artemisia a monophyletic genus, sister to genus Chrysanthemum. We also identified 20 highly polymorphic regions including rpoC2-rps2, trnR-UCU-trnG-UCC, rps18-rpl20, and trnL-UAG-rpl32 that could be used to develop authentic and cost-effective markers to resolve taxonomic discrepancies and infer phylogenetic relationships among Artemisia species.