Detecting useful genetic markers and reconstructing the phylogeny of an important medicinal resource plant, Artemisia selengensis, based on chloroplast genomics

The chloroplast genome of Chrozophora sabulosa Kar. & Kir. and its exploration in the evolutionary position uncertainty of genus Chrozophora

Chrozophora sabulosa Kar. & Kir. is a biennial herbaceous plant that belongs to the Euphorbiaceae family and has medicinal properties. This research aimed to identify the genetic characteristics and phylogenetic position of the Chrozophora genus within the Euphorbiaceae family. The evolutionary position of the Chrozophora genus was previously unknown due to insufficient research. Therefore, to determine the evolutionary link between C. sabulosa and other related species, we conducted a study using the NGS Illumina platform to sequence the C. sabulosa chloroplast (cp.) genome. The study results showed that the genome was 156,488 bp in length. It had a quadripartite structure consisting of two inverted repeats (IRb and IRa) of 24,649-bp, separated by an 87,696-bp LSC region and a 19,494-bp SSC region. The CP genome contained 113 unique genes, including four rRNA genes, 30 tRNA genes, and 79 CDS genes. In the second copy of the inverted repeat, there were 18 duplicated genes. The C. sabulosa lacks the petD, petB, rpl2, and rps16 intron. The analysis of simple sequence repeats (SSRs) revealed 93 SSR loci of 22 types and 78 oligonucleotide repeats of four kinds. The phylogenetic investigation showed that the Chrozophora genus evolved paraphyletically from other members of the Euphorbiaceae family. To support the phylogenetic findings, we selected species from the Euphorbiaceae and Phyllanthaceae families to compare with C. sabulosa for Ks and Ka substitution rates, InDels investigation, IR contraction and expansion, and SNPs analysis. The results of these comparative studies align with the phylogenetic findings. We identified six highly polymorphic regions shared by both families, which could be used as molecular identifiers for the Chrozophora genus (rpl33-rps18, rps18-rpl20, rps15-ycf1, ndhG-ndhI, psaI-ycf4, petA-psbJ). The cp. genome sequence of C. sabulosa reveals the evolution of plastid sequences in Chrozophora species. This is the first time the cp. genome of a Chrozophora genus has been sequenced, serving as a foundation for future sequencing of other species within the Chrozophoreae tribe and facilitating in-depth taxonomic research. The results of this research will also aid in identifying new Chrozophora species.

Comparative and phylogenetic analysis of the complete chloroplast genomes of 10 Artemisia selengensis resources based on high-throughput sequencing

BACKGROUND

Artemisia selengensis, classified within the genus Artemisia of the Asteraceae family, is a perennial herb recognized for its dual utility in culinary and medicinal domains. There are few studies on the chloroplast genome of A. selengensis, and the phylogeographic classification is vague, which makes phylogenetic analysis and evolutionary studies very difficult.

RESULTS

The chloroplast genomes of 10 A. selengensis in this study were highly conserved in terms of gene content, gene order, and gene intron number. The genome lengths ranged from 151,148 to 151,257 bp and were typical of a quadripartite structure with a total GC content of approximately 37.5%. The chloroplast genomes of all species encode 133 genes, including 88 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. Due to the contraction and expansion of the inverted repeats (IR), the overlap of ycf1 and ndhF genes occurred at the inverted repeats B (IRB) and short single copy sequence (SSC) boundaries. According to a codon use study, the frequent base in the chloroplast genome of A. selengensis' third codon position was A/T. The number of SSR repeats was 42-44, most of which were single nucleotide A/T repeats. Sequence alignment analysis of the chloroplast genome showed that variable regions were mainly distributed in single copy regions, nucleotide diversity values of 0 to 0.009 were calculated by sliding window analysis, 8 mutation hotspot regions were detected, and coding regions were more conserved than non-coding regions. Analysis of non-synonymous substitution (Ka) and synonymous substitution (Ks) revealed that accD, rps12, petB, and atpF genes were affected by positive selection and no genes were affected by neutral selection. Based on the findings of the phylogenetic analysis, Artemisia selengensis was sister to the genus Artemisia Chrysanthemum and formed a monophyletic group with other Artemisia genera.

CONCLUSIONS

In this research, the present study systematically compared the chloroplast genomic features of A. selengensis and provided important information for the study of the chloroplast genome of A. selengensis and the evolutionary relationships among Asteraceae species.

Divergence of 10 satellite repeats in Artemisia (Asteraceae: Anthemideae) based on sequential fluorescence in situ hybridization analysis: evidence for species identification and evolution

Artemisia is a large genus encompassing about 400 diverse species, many of which have considerable medicinal and ecological value. However, complex morphological information and variation in ploidy level and nuclear DNA content have presented challenges for evolution studies of this genus. Consequently, taxonomic inconsistencies within the genus persist, hindering the utilization of such large plant resources. Researchers have utilized satellite DNAs to aid in chromosome identification, species classification, and evolutionary studies due to their significant sequence and copy number variation between species and close relatives. In the present study, the RepeatExplorer2 pipeline was utilized to identify 10 satellite DNAs from three species (Artemisia annua, Artemisia vulgaris, Artemisia viridisquama), and fluorescence in situ hybridization confirmed their distribution on chromosomes in 24 species, including 19 Artemisia species with 5 outgroup species from Ajania and Chrysanthemum. Signals of satellite DNAs exhibited substantial differences between species. We obtained one genus-specific satellite from the sequences. Additionally, molecular cytogenetic maps were constructed for Artemisia vulgaris, Artemisia leucophylla, and Artemisia viridisquama. One species (Artemisia verbenacea) showed a FISH distribution pattern suggestive of an allotriploid origin. Heteromorphic FISH signals between homologous chromosomes in Artemisia plants were observed at a high level. Additionally, the relative relationships between species were discussed by comparing ideograms. The results of the present study provide new insights into the accurate identification and taxonomy of the Artemisia genus using molecular cytological methods.

Complete chloroplast genomes of eight Artemisia species: Comparative analysis, molecular identification, and phylogenetic analysis

Artemisia L. is the largest genus in the Asteraceae, and well known for its high medicinal value. The morphological features of Artemisia species are similar, making taxonomic identification and evolutionary research difficult. We sequenced chloroplast genomes of eight Artemisia species, all of which are common adulterants of A. argyi. We used novel genetic data and compared these data to the published A. argyi chloroplast genome in to develop molecular markers for species identification and reconstructing phylogenetic relationships between Artemisia species. The eight chloroplast sequences were highly similar in gene order, content, and structure, encoding a total of 114 genes (82 protein-coding genes, 28 tRNAs, and four rRNAs). All species harboured similar repeat sequences and simple sequence repeats (SSRs), ranging from 47 to 49 and 38 to 40 repeats, respectively. In addition, we identified five hypervariable regions (rpl32-trnL, rps16-trnQ, petN-psbM, trnE-rpoB, and atpA-trnR) and ten variable coding genes (ycf1, psbG, rpl36, psaC, psaI, accD, psbT, ndhD, ndhE, and psbH), which can be used to develop chloroplast molecular markers. Finally, phylogenetic reconstructions based on six datasets produced similar topologies, revealing A. argyi is closely related to species often found as adulterants, as expected. Our research provides valuable new information on the evolution and phylogenetic relationships between Artemisia chloroplast genomes and identifies valuable molecular makers to distinguish it from closely related species.

The complete chloroplast genome sequences of four Liparis species (Orchidaceae) and phylogenetic implications

Liparis Richard (Malaxideae, Epidendroideae) is a large and diverse genus of the family Orchidaceae, the taxonomy of which is complicated and controversial. In this study, we sequenced, assembled and analyzed four complete chloroplast genomes of Liparis species including L. kumokiri, L. makinoana, L. pauliana, and L. viridiflora, and evaluated their phylogenetic relationships with related species for the first time. These four chloroplast genomes (size range 153,095 to 158,239 bp) possess typical quadripartite structures that consist of a large single copy (LSC, 83,533-86,752 bp), a small single copy (SSC, 17,938-18,156 bp) and a pair of inverted repeats (IRs, 26,421-26,933 bp). The genomes contain 133 genes, including 87 protein coding genes, 38 tRNAs and 8 rRNA genes. The genome arrangements, gene contents, gene order, long repeats and simple sequence repeats were similar with small differences observed among these four chloroplast genomes. Five highly variable regions including ycf1, ndhA, ndhF, trnQ and trnK were identified from the comparative analysis with other nine related Liparis species, which had the potential to be used as DNA markers for species identification and phylogenetic studies of Liparis species. Phylogenetic analysis based on the complete chloroplast genome sequences strongly supported the polyphyly of Liparis and its further division into three branches. These results provided valuable information to illustrate the complicated taxonomy, phylogeny and evolution process of the Liparis genus.

Characterization of the Chloroplast Genome of Argyranthemum frutescens and a Comparison with Other Species in Anthemideae

Argyranthemum frutescens, which belongs to the Anthemideae (Asteraceae), is widely cultivated as an ornamental plant. In this study, the complete chloroplast genome of A. frutescens was obtained based on the sequences generated by Illumina HiSeq. The chloroplast genome of A. frutescens was 149,626 base pairs (bp) in length, containing a pair of inverted repeats (IR, 24,510 bp) regions separated by a small single-copy (SSC, 18,352 bp) sequence and a large single-copy (LSC, 82,254 bp) sequence. The genome contained 132 genes, consisting of 85 coding DNA sequences, 37 tRNA genes, and 8 rRNA genes, with nineteen genes duplicated in the IR region. A comparison chloroplast genome analysis among ten species from the tribe of Anthemideae revealed that the chloroplast genome size varied, but the genome structure, gene content, and oligonucleotide repeats were highly conserved. Highly divergent regions, e.g., ycf1, trnK-psbK, petN-psbM intronic, were detected. Phylogenetic analysis supported Argyranthemum as a separate genus. The findings of this study will be helpful in the exploration of the phylogenetic relationships of the tribe of Anthemideae and contribute to the breeding improvement of A. frutescens.

Plastome structure of 8 Calanthe s.l. species (Orchidaceae): comparative genomics, phylogenetic analysis

Background

Calanthe (Epidendroideae, Orchidaceae) is a pantropical genus distributed in Asia and Africa. Its species are of great importance in terms of economic, ornamental and medicinal values. However, due to limited and confusing delimitation characters, the taxonomy of the Calanthe alliance (Calanthe, Cephalantheropsis, and Phaius) has not been sufficiently resolved. Additionally, the limited genomic information has shown incongruences in its systematics and phylogeny. In this study, we used illumina platform sequencing, performed a de novo assembly, and did a comparative analysis of 8 Calanthe group species' plastomes: 6 Calanthe and 2 Phaius species. Phylogenetic analyses were used to reconstruct the relationships of the species as well as with other species of the family Orchidaceae.

Results

The complete plastomes of the Calanthe group species have a quadripartite structure with varied sizes ranging between 150,105bp-158,714bp, including a large single-copy region (LSC; 83,364bp- 87,450bp), a small single-copy region (SSC; 16,297bp -18,586bp), and a pair of inverted repeat regions (IRs; 25,222bp - 26,430bp). The overall GC content of these plastomes ranged between 36.6-36.9%. These plastomes encoded 131-134 differential genes, which included 85-88 protein-coding genes, 37-38 tRNA genes, and 8 rRNA genes. Comparative analysis showed no significant variations in terms of their sequences, gene content, gene order, sequence repeats and the GC content hence highly conserved. However, some genes were lost in C. delavayi (P. delavayi), including ndhC, ndhF, and ndhK genes. Compared to the coding regions, the non-coding regions had more sequence repeats hence important for species DNA barcoding. Phylogenetic analysis revealed a paraphyletic relationship in the Calanthe group, and confirmed the position of Phaius delavayi in the genus Calanthe as opposed to its previous placement in Phaius.

Conclusion

This study provides a report on the complete plastomes of 6 Calanthe and 2 Phaius species and elucidates the structural characteristics of the plastomes. It also highlights the power of plastome data to resolve phylogenetic relationships and clarifies taxonomic disputes among closely related species to improve our understanding of their systematics and evolution. Furthermore, it also provides valuable genetic resources and a basis for studying evolutionary relationships and population genetics among orchid species.

Plastome structure of 8 Calanthe s.l. species (Orchidaceae): comparative genomics, phylogenetic analysis

BACKGROUND

Calanthe (Epidendroideae, Orchidaceae) is a pantropical genus distributed in Asia and Africa. Its species are of great importance in terms of economic, ornamental and medicinal values. However, due to limited and confusing delimitation characters, the taxonomy of the Calanthe alliance (Calanthe, Cephalantheropsis, and Phaius) has not been sufficiently resolved. Additionally, the limited genomic information has shown incongruences in its systematics and phylogeny. In this study, we used illumina platform sequencing, performed a de novo assembly, and did a comparative analysis of 8 Calanthe group species' plastomes: 6 Calanthe and 2 Phaius species. Phylogenetic analyses were used to reconstruct the relationships of the species as well as with other species of the family Orchidaceae.

RESULTS

The complete plastomes of the Calanthe group species have a quadripartite structure with varied sizes ranging between 150,105bp-158,714bp, including a large single-copy region (LSC; 83,364bp- 87,450bp), a small single-copy region (SSC; 16,297bp -18,586bp), and a pair of inverted repeat regions (IRs; 25,222bp - 26,430bp). The overall GC content of these plastomes ranged between 36.6-36.9%. These plastomes encoded 131-134 differential genes, which included 85-88 protein-coding genes, 37-38 tRNA genes, and 8 rRNA genes. Comparative analysis showed no significant variations in terms of their sequences, gene content, gene order, sequence repeats and the GC content hence highly conserved. However, some genes were lost in C. delavayi (P. delavayi), including ndhC, ndhF, and ndhK genes. Compared to the coding regions, the non-coding regions had more sequence repeats hence important for species DNA barcoding. Phylogenetic analysis revealed a paraphyletic relationship in the Calanthe group, and confirmed the position of Phaius delavayi in the genus Calanthe as opposed to its previous placement in Phaius.

CONCLUSION

This study provides a report on the complete plastomes of 6 Calanthe and 2 Phaius species and elucidates the structural characteristics of the plastomes. It also highlights the power of plastome data to resolve phylogenetic relationships and clarifies taxonomic disputes among closely related species to improve our understanding of their systematics and evolution. Furthermore, it also provides valuable genetic resources and a basis for studying evolutionary relationships and population genetics among orchid species.

Comparative Analysis of Chloroplast Genome Structure and Phylogenetic Relationships Among Six Taxa Within the Genus Catalpa (Bignoniaceae)

Species within the Genus Catalpa are mostly semievergreen or deciduous trees with opposite or whorled leaves. C. bungei, C. fargesii f. duclouxii and C. fargesii are sources of traditional precious wood in China, known as the “kings of wood”. Due to a lack of phenotypic and molecular studies and insufficient sequence information, intraspecific morphological differences, common DNA barcodes and partial sequence fragments cannot clearly reveal the phylogenetic or intraspecific relationships within Catalpa. Therefore, we sequenced the complete chloroplast genomes of six taxa of the genus Catalpa and analyzed their basic structure and evolutionary relationships. The chloroplast genome of Catalpa shows a typical tetrad structure with a total length ranging from 157,765 bp (C. fargesii) to 158,355 bp (C. ovata). The length of the large single-copy (LSC) region ranges from 84,599 bp (C. fargesii) to 85,004 bp (C. ovata), that of the small single-copy (SSC) region ranges from 12,662 bp (C. fargesii) to 12,675 bp (C. ovata), and that of the inverted repeat (IR) regions ranges from 30,252 bp (C. fargesii) to 30,338 bp (C. ovata). The GC content of the six chloroplast genomes were 38.1%. In total, 113 unique genes were detected, and there were 19 genes in IR regions. The 113 genes included 79 protein-coding genes, 30 tRNA genes and four rRNA genes. Five hypervariable regions (trnH-psbA, rps2-rpoC2, rpl22, ycf15-trnl-CAA and rps15) were identified by analyzing chloroplast nucleotide polymorphisms, which might be serve as potential DNA barcodes for the species. Comparative analysis showed that single nucleotide polymorphisms (SNPs) and simple sequence repeats (SSRs) were highly diverse in the six species. Codon usage patterns were highly similar among the taxa included in the present study. In addition to the stop codons, all codons showed a preference for ending in A or T. Phylogenetic analysis of the entire chloroplast genome showed that all taxa within the genus Catalpa formed a monophyletic group, clearly reflecting the relationships within the genus. This study provides information on the chloroplast genome sequence, structural variation, codon bias and phylogeny of Catalpa, which will facilitate future research efforts.

Analysis of the Complete Plastomes of 31 Species of Hoya Group: Insights Into Their Comparative Genomics and Phylogenetic Relationships

Hoya is a genus in Apocynaceae-Asclepiadoideae, known for its showy wax flowers, making it a popular ornamental plant. However, phylogenetic relationships among most Hoya species are not yet fully resolved. In this study, we sequenced 31 plastomes of Hoya group species using genome skimming data and carried out multiple analyses to understand genome variation to resolve the phylogenetic positions of some newly sequenced Chinese endemic species. We also screened possible hotspots, trnT-trnL-trnF, psba-trnH, and trnG-UCC, ndhF, ycf1, matK, rps16, and accD genes that could be used as molecular markers for DNA barcoding and species identification. Using maximum likelihood (ML) and Bayesian Inference (BI), a species phylogeny was constructed. The newly assembled plastomes genomes showed the quasi-tripartite structure characteristic for Hoya and Dischidia with a reduced small single copy (SSC) and extremely enlarged inverted repeats (IR). The lengths ranged from 175,404 bp in Hoya lacunosa to 179,069 bp in H. ariadna. The large single copy (LSC) regions ranged from 80,795 bp (Hoya liangii) to 92,072 bp (Hoya_sp2_ZCF6006). The massively expanded IR regions were relatively conserved in length, with the small single-copy region reduced to a single gene, ndhF. We identified 235 long dispersed repeats (LDRs) and ten highly divergent hotspots in the 31 Hoya plastomes, which can be used as DNA barcodes for species identification. The phylogeny supports Clemensiella as a distinct genus. Hoya ignorata is resolved as a relative to Clade VI species. This study discloses the advantages of using Plastome genome data to study phylogenetic relationships.

Comparative genomics and phylogenetic relationships of two endemic and endangered species (Handeliodendron bodinieri and Eurycorymbus cavaleriei) of two monotypic genera within Sapindales

BACKGROUND

Handeliodendron Rehder and Eurycorymbus Hand.-Mazz. are the monotypic genera in the Sapindaceae family. The phylogenetic relationship of these endangered species Handeliodendron bodinieri (Lévl.) Rehd. and Eurycorymbus cavaleriei (Lévl.) Rehd. et Hand.-Mazz. with other members of Sapindaceae s.l. is not well resolved. A previous study concluded that the genus Aesculus might be paraphyletic because Handeliodendron was nested within it based on small DNA fragments. Thus, their chloroplast genomic information and comparative genomic analysis with other Sapindaceae species are necessary and crucial to understand the circumscription and plastome evolution of this family.

RESULTS

The chloroplast genome sizes of Handeliodendron bodinieri and Eurycorymbus cavaleriei are 151,271 and 158,690 bp, respectively. Results showed that a total of 114 unique genes were annotated in H. bodinieri and E. cavaleriei, and the ycf1 gene contained abundant SSRs in both genomes. Comparative analysis revealed that gene content, PCGs, and total GC content were remarkably similar or identical within 13 genera from Sapindaceae, and the chloroplast genome size of four genera was generally smaller within the family, including Acer, Dipteronia, Aesculus, and Handeliodendron. IR boundaries of the H. bodinieri showed a significant contraction, whereas it presented a notable expansion in E. cavaleriei cp genome. Ycf1, ndhC-trnV-UAC, and rpl32-trnL-UAG-ccsA were remarkably divergent regions in the Sapindaceae species. Analysis of selection pressure showed that there are a few positively selected genes. Phylogenetic analysis based on different datasets, including whole chloroplast genome sequences, coding sequences, large single-copy, small single-copy, and inverted repeat regions, consistently demonstrated that H. bodinieri was sister to the clade consisting of Aesculus chinensis and A. wangii and strongly support Eurycorymbus cavaleriei as sister to Dodonaea viscosa.

CONCLUSION

This study revealed that the cp genome size of the Hippocastanoideae was generally smaller compared to the other subfamilies within Sapindaceae, and three highly divergent regions could be used as the specific DNA barcodes within Sapindaceae. Phylogenetic results strongly support that the subdivision of four subfamilies within Sapindaceae, and Handeliodendron is not nested within the genus Aesculus.

Exploring the evolutionary characteristics between cultivated tea and its wild relatives using complete chloroplast genomes

Background

Cultivated tea is one of the most important economic and ecological trees distributed worldwide. Cultivated tea suffer from long-term targeted selection of traits and overexploitation of habitats by human beings, which may have changed its genetic structure. The chloroplast is an organelle with a conserved cyclic genomic structure, and it can help us better understand the evolutionary relationship of Camellia plants.

Results

We conducted comparative and evolutionary analyses on cultivated tea and wild tea, and we detected the evolutionary characteristics of cultivated tea. The chloroplast genome sizes of cultivated tea were slightly different, ranging from 157,025 to 157,100 bp. In addition, the cultivated species were more conserved than the wild species, in terms of the genome length, gene number, gene arrangement and GC content. However, comparing Camellia sinensis var. sinensis and Camellia sinensis var. assamica with their cultivars, the IR length variation was approximately 20 bp and 30 bp, respectively. The nucleotide diversity of 14 sequences in cultivated tea was higher than that in wild tea. Detailed analysis on the genomic variation and evolution of Camellia sinensis var. sinensis cultivars revealed 67 single nucleotide polymorphisms (SNPs), 46 insertions/deletions (indels), and 16 protein coding genes with nucleotide substitutions, while Camellia sinensis var. assamica cultivars revealed 4 indels. In cultivated tea, the most variable gene was ycf1. The largest number of nucleotide substitutions, five amino acids exhibited site-specific selection, and a 9 bp sequence insertion were found in the Camellia sinensis var. sinensis cultivars. In addition, phylogenetic relationship in the ycf1 tree suggested that the ycf1 gene has diverged in cultivated tea. Because C. sinensis var. sinensis and its cultivated species were not tightly clustered.

Conclusions

The cultivated species were more conserved than the wild species in terms of architecture and linear sequence order. The variation of the chloroplast genome in cultivated tea was mainly manifested in the nucleotide polymorphisms and sequence insertions. These results provided evidence regarding the influence of human activities on tea.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01800-1.

Analysis of the Complete Plastomes of 31 Species of Hoya Group: Insights Into Their Comparative Genomics and Phylogenetic Relationships

Hoya is a genus in Apocynaceae-Asclepiadoideae, known for its showy wax flowers, making it a popular ornamental plant. However, phylogenetic relationships among most Hoya species are not yet fully resolved. In this study, we sequenced 31 plastomes of Hoya group species using genome skimming data and carried out multiple analyses to understand genome variation to resolve the phylogenetic positions of some newly sequenced Chinese endemic species. We also screened possible hotspots, trnT-trnL-trnF, psba-trnH, and trnG-UCC, ndhF, ycf1, matK, rps16, and accD genes that could be used as molecular markers for DNA barcoding and species identification. Using maximum likelihood (ML) and Bayesian Inference (BI), a species phylogeny was constructed. The newly assembled plastomes genomes showed the quasi-tripartite structure characteristic for Hoya and Dischidia with a reduced small single copy (SSC) and extremely enlarged inverted repeats (IR). The lengths ranged from 175,404 bp in Hoya lacunosa to 179,069 bp in H. ariadna. The large single copy (LSC) regions ranged from 80,795 bp (Hoya liangii) to 92,072 bp (Hoya_sp2_ZCF6006). The massively expanded IR regions were relatively conserved in length, with the small single-copy region reduced to a single gene, ndhF. We identified 235 long dispersed repeats (LDRs) and ten highly divergent hotspots in the 31 Hoya plastomes, which can be used as DNA barcodes for species identification. The phylogeny supports Clemensiella as a distinct genus. Hoya ignorata is resolved as a relative to Clade VI species. This study discloses the advantages of using Plastome genome data to study phylogenetic relationships.

The large single-copy (LSC) region functions as a highly effective and efficient molecular marker for accurate authentication of medicinal Dendrobium species

Having great medicinal values, Dendrobium species of "Fengdou" (DSFs) are a taxonomically complex group in Dendrobium genus including many closely related and recently diverged species. Traditionally used DNA markers have been proved to be insufficient in authenticating many species of this group. Here, we investigated 101 complete plastomes from 23 DSFs, comprising 72 newly sequenced and 29 documented, which all exhibited well-conserved genomic organization and gene order. Plastome-wide comparison showed the co-occurrence of single nucleotide polymorphisms (SNPs) and insertions/deletions (indels), which can be explained by both the repeat-associated and indel-associated mutation hypotheses. Moreover, guanine-cytosine (GC) content was found to be negatively correlated with the three divergence variables (SNPs, indels and repeats), indicating that GC content may reflect the level of the local sequence divergence. Our species authentication analyses revealed that the relaxed filtering strategies of sequence alignment had no negative impact on species identification. By assessing the maximum likelihood (ML) trees inferred from different datasets, we found that the complete plastome and large single-copy (LSC) datasets both successfully identified all 23 DSFs with the maximum bootstrap values. However, owing to the high efficiency of LSC in species identification, we recommend using LSC for accurate authentication of DSFs.

ycf1-ndhF genes, the most promising plastid genomic barcode, sheds light on phylogeny at low taxonomic levels in Prunus persica

Background

Chloroplast genome sequencing is becoming a valuable process for developing several DNA barcodes. At present, plastid DNA barcode for systematics and evolution in flowering plant rely heavily on the use of non-coding genes. The present study was performed to verify the novelty and suitability of the two hotspot barcode plastid coding gene ycf1 and ndhF, to estimate the rate of molecular evolution in the Prunus genus at low taxonomic levels.

Results

Here, 25 chloroplast genomes of Prunus genus were selected for sequences annotation to search for the highly variable coding DNA barcode regions. Among them, 5 genera were of our own data, including the ornamental, cultivated, and wild haplotype, while 20 genera have been downloaded from the GenBank database. The results indicated that the two hotspot plastid gene ycf1 and ndhF were the most variable regions within the coding genes in Prunus with an average of 3268 to 3416 bp in length, which have been predicted to have the highest nucleotide diversity, with the overall transition/transversion bias (R = 1.06). The ycf1-ndhF structural domains showed a positive trend evident in structure variation among the 25 specimens tested, due to the variant overlap^’s gene annotation and insertion or deletion with a broad trend of the full form of IGS sequence. As a result, the principal component analysis (PCA) and the ML tree data drew an accurate monophyletic annotations cluster in Prunus species, offering unambiguous identification without overlapping groups between peach, almond, and cherry.

Conclusion

To this end, we put forward the domain of the two-locus ycf1-ndhF genes as the most promising coding plastid DNA barcode in P. persica at low taxonomic levels. We believe that the discovering of further variable loci with high evolutionary rates is extremely useful and potential uses as a DNA barcode in P. persica for further phylogeny study and species identification.

A comparative analysis of the complete chloroplast genomes of three Chrysanthemum boreale strains

BACKGROUND

Chrysanthemum boreale Makino (Anthemideae, Asteraceae) is a plant of economic, ornamental and medicinal importance. We characterized and compared the chloroplast genomes of three C. boreale strains. These were collected from different geographic regions of Korea and varied in floral morphology.

METHODS

The chloroplast genomes were obtained by next-generation sequencing techniques, assembled de novo, annotated, and compared with one another. Phylogenetic analysis placed them within the Anthemideae tribe.

RESULTS

The sizes of the complete chloroplast genomes of the C. boreale strains were 151,012 bp (strain 121002), 151,098 bp (strain IT232531) and 151,010 bp (strain IT301358). Each genome contained 80 unique protein-coding genes, 4 rRNA genes and 29 tRNA genes. Comparative analyses revealed a high degree of conservation in the overall sequence, gene content, gene order and GC content among the strains. We identified 298 single nucleotide polymorphisms (SNPs) and 106 insertions/deletions (indels) in the chloroplast genomes. These variations were more abundant in non-coding regions than in coding regions. Long dispersed repeats and simple sequence repeats were present in both coding and noncoding regions, with greater frequency in the latter. Regardless of their location, these repeats can be used for molecular marker development. Phylogenetic analysis revealed the evolutionary relationship of the species in the Anthemideae tribe. The three complete chloroplast genomes will be valuable genetic resources for studying the population genetics and evolutionary relationships of Asteraceae species.

Comparative chloroplast genome analysis of Artemisia (Asteraceae) in East Asia: insights into evolutionary divergence and phylogenomic implications

BACKGROUND

Artemisia in East Asia includes a number of economically important taxa that are widely used for food, medicinal, and ornamental purposes. The identification of taxa, however, has been hampered by insufficient diagnostic morphological characteristics and frequent natural hybridization. Development of novel DNA markers or barcodes with sufficient resolution to resolve taxonomic issues of Artemisia in East Asia is significant challenge.

RESULTS

To establish a molecular basis for taxonomic identification and comparative phylogenomic analysis of Artemisia, we newly determined 19 chloroplast genome (plastome) sequences of 18 Artemisia taxa in East Asia, de novo-assembled and annotated the plastomes of two taxa using publicly available Illumina reads, and compared them with 11 Artemisia plastomes reported previously. The plastomes of Artemisia were 150,858-151,318 base pairs (bp) in length and harbored 87 protein-coding genes, 37 transfer RNAs, and 8 ribosomal RNA genes in conserved order and orientation. Evolutionary analyses of whole plastomes and 80 non-redundant protein-coding genes revealed that the noncoding trnH-psbA spacer was highly variable in size and nucleotide sequence both between and within taxa, whereas the coding sequences of accD and ycf1 were under weak positive selection and relaxed selective constraints, respectively. Phylogenetic analysis of the whole plastomes based on maximum likelihood and Bayesian inference analyses yielded five groups of Artemisia plastomes clustered in the monophyletic subgenus Dracunculus and paraphyletic subgenus Artemisia, suggesting that the whole plastomes can be used as molecular markers to infer the chloroplast haplotypes of Artemisia taxa. Additionally, analysis of accD and ycf1 hotspots enabled the development of novel markers potentially applicable across the family Asteraceae with high discriminatory power.

CONCLUSIONS

The complete sequences of the Artemisia plastomes are sufficiently polymorphic to be used as super-barcodes for this genus. It will facilitate the development of new molecular markers and study of the phylogenomic relationships of Artemisia species in the family Asteraceae.

Comparative Analysis of the Complete Chloroplast Genome of Mainland Aster spathulifolius and Other Aster Species

Aster spathulifolius, a common ornamental and medicinal plant, is widely distributed in Korea and Japan, and is genetically classified into mainland and island types. Here, we sequenced the whole chloroplast genome of mainland A. spathulifolius and compared it with those of the island type and other Aster species. The chloroplast genome of mainland A. spathulifolius is 152,732 bp with a conserved quadripartite structure, has 37.28% guanine-cytosine (GC) content, and contains 114 non-redundant genes. Comparison of the chloroplast genomes between the two A. spathulifolius lines and the other Aster species revealed that their sequences, GC contents, gene contents and orders, and exon-intron structure were well conserved; however, differences were observed in their lengths, repeat sequences, and the contraction and expansion of the inverted repeats. The variations were mostly in the single-copy regions and non-coding regions, which, together with the detected simple sequence repeats, could be used for the development of molecular markers to distinguish between these plants. All Aster species clustered into a monophyletic group, but the chloroplast genome of mainland A. spathulifolius was more similar to the other Aster species than to that of the island A. spathulifolius. The accD and ndhF genes were detected to be under positive selection within the Aster lineage compared to other related taxa. The complete chloroplast genome of mainland A. spathulifolius presented in this study will be helpful for species identification and the analysis of the genetic diversity, evolution, and phylogenetic relationships in the Aster genus and the Asteraceae.

Initial Characterization of the Chloroplast Genome of Vicia sepium, an Important Wild Resource Plant, and Related Inferences About Its Evolution

Lack of complete genomic information concerning Vicia sepium (Fabaceae: Fabeae) precludes investigations of evolution and populational diversity of this perennial high-protein forage plant suitable for cultivation in extreme conditions. Here, we present the complete and annotated chloroplast genome of this important wild resource plant. V. sepium chloroplast genome includes 76 protein-coding genes, 29 tRNA genes, 4 rRNA genes, and 1 pseudogene. Its 124,095 bp sequence has a loss of one inverted repeat (IR). The GC content of the whole genome, the protein-coding, intron, tRNA, rRNA, and intergenic spacer regions was 35.0%, 36.7%, 34.6%, 52.3%, 54.2%, and 29.2%, respectively. Comparative analyses with plastids from related genera belonging to Fabeae demonstrated that the greatest variation in the V. sepium genome length occurred in protein-coding regions. In these regions, some genes and introns were lost or gained; for example, ycf4, clpP intron, and rpl16 intron deletions and rpl20 and ORF292 insertions were observed. Twelve highly divergent regions, 66 simple sequence repeats (SSRs) and 27 repeat sequences were also found in these regions. Detailed evolutionary rate analysis of protein-coding genes showed that Vicia species exhibit additional interesting characteristics including positive selection of ccsA, clpP, rpl32, rpl33, rpoC1, rps15, rps2, rps4, and rps7, and the evolutionary rates of atpA, accD, and rps2 in Vicia are significantly accelerated. These genes are important candidate genes for understanding the evolutionary strategies of Vicia and other genera in Fabeae. The phylogenetic analysis showed that Vicia and Lens are included in the same clade and that Vicia is paraphyletic. These results provide evidence regarding the evolutionary history of the chloroplast genome.

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.