Complete chloroplast genome sequences of Mongolia medicine Artemisia frigida and phylogenetic relationships with other plants

Exploring Structural Plastome Evolution in Asterales: Insights from Off-Target Hybrid Enrichment Data on the Small Single-Copy Region

The massive increase in the amount of plastid genome data have allowed researchers to address a variety of evolutionary questions within a wide range of plant groups. While plastome structure is generally conserved, some angiosperm lineages exhibit structural changes. Such is the case of the megadiverse order Asterales, where rearrangements in plastome structure have been documented. This study investigates the possibility of recovering plastid loci from off-target reads obtained through hybrid enrichment techniques. Our sampling includes 63 species from the eleven currently recognized families in Asterales derived from previously published studies. We assembled and annotated complete and partial plastomes using custom pipelines and estimate phylogenomic relationships. We retrieved plastid information from 60 of the 63 sampled species including a complete plastome from Tithonia tubaeformis (Asteraceae), circular partial (with gaps) plastomes from seven species, and non-circular partial plastomes from other 52 species. We focused on the small single-copy region because it could be recovered for over 29 species. Within the small single-copy region, we assessed intron losses and presence of putative pseudogenes. Comparative genomics revealed a relocated fragment of ~ 6500 bp in two Campanulaceae lineages (i. e. subfamily Lobelioideae and Pseudonemacladus oppositifolium), involving the genes rbcL, atpB, atpE, trnM-CAU, and trnV-UAC. Obtained phylogenetic hypotheses were congruent across the applied methods and consistent with previously published results. Our study demonstrates the feasibility of recovering plastid information, both complete and partial, from off-target hybrid enrichment data and provides insights on the structural plastome changes that have occurred throughout the evolution of the order Asterales.

Kazakhstan tulips: comparative analysis of complete chloroplast genomes of four local and endangered species of the genus Tulipa L

Species of Tulipa are important ornamental plants used for horticultural purposes in various countries, across Asia, Europe, and North Africa. The present study is the first report on typical features of the complete chloroplast genome sequence of four local and endangered species including T. alberti, T. kaufmanniana, T. greigii, and T. dubia from Kazakhstan using Illumina sequencing technology. The comparative analyses revealed that the complete genomes of four species were highly conserved in terms of total genome size (152. 006 bp - 152. 382 bp), including a pair of inverted repeat regions (26. 330 bp - 26. 371 bp), separated by a large single copy region (82.169 bp - 82,378 bp) and a small copy region (17.172 bp -17.260 bp). Total GC content (36.58-36.62 %), gene number (131), and intron length (540 bp - 2620 bp) of 28 genes. The complete genomes of four species showed nucleotide diversity (π =0,003257). The total number of SSR loci was 159 in T. alberti, 158 in T. kaufmanniana, 174 in T. greigii, and 163 in T. dubia. The result indicated that ten CDS genes, namely rpoC2, cemA, rbcL, rpl36, psbH, rps3, rpl22, ndhF, ycf1, and matK, with effective polymorphic simple sequence repeats (SSRs), high sequence variability (SV) ranging from 2.581 to 6.102, and high nucleotide diversity (Pi) of these loci ranging from 0,004 to 0,010. For all intergenic regions longer than 150 bp, twenty one most variable regions were found with high sequence variability (SV) ranging from 4,848 to 11,862 and high nucleotide diversity (Pi) ranging from 0,01599 to 0,01839. Relative synonymous codon usage (RSCU) analysis was used to identify overrepresented and underrepresented codons for each amino acid. Based on the phylogenic analysis, the sequences clustered into four major groups, reflecting distinct evolutionary lineages corresponding to the subgenera Eriostemons, Tulipa, and Orithyia. Notably, T. greigii was distinctively grouped with species from Orithyia and Eriostemons rather than with other Tulipa species, suggesting a unique evolutionary history potentially shaped by geographical isolation or specific ecological pressures. The complete chloroplast genome of the four Tulipa species provides fundamental information for future research studies, even for designing the high number of available molecular markers.

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.

The first complete chloroplast genome for the species-rich genus Anthemis (Asteraceae)

Anthemis L. 1753 is a species-rich genus of Asteraceae. It includes species of economic value as ornamental, edible and medicinal plants. We report here the complete chloroplast genome of Anthemis cretica subsp. calabrica (Arcang.) R. Fern. 1975, the first chloroplast genome of any Anthemis species. For the scope, we have used long reads obtained with Oxford Nanopore technology. The assembled plastome is 149,509 bp long and is subdivided in a large single-copy region (LSC, 82,317 bp), a small single-copy region (SSC, 18,212 bp) and two inverted repeat regions (IR: 24,490 bp). It comprehends 136 predicted genes, of which 92 protein-coding genes, 35 tRNAs, and 8 rRNAs. The overall GC content is 37.54%. The position of A. cretica subsp. calabrica in the phylogenetic tree is congruent with those of previous analyses based on a few chloroplast or nuclear regions.

Comparative analysis of complete chloroplast genomes of Synotis species (Asteraceae, Senecioneae) for identification and phylogenetic analysis

BACKGROUND

The Synotis (C. B. Clarke) C. Jeffrey & Y. L. Chen is an ecologically important genus of the tribe Senecioneae, family Asteraceae. Because most species of the genus bear similar morphology, traditional morphological identification methods are very difficult to discriminate them. Therefore, it is essential to develop a reliable and effective identification method for Synotis species. In this study, the complete chloroplast (cp.) genomes of four Synotis species, S. cavaleriei (H.Lév.) C. Jeffrey & Y.L. Chen, S. duclouxii (Dunn) C. Jeffrey & Y.L. Chen, S. nagensium (C.B. Clarke) C. Jeffrey & Y.L. Chen and S. erythropappa (Bureau & Franch.) C. Jeffrey & Y. L. Chen had been sequenced using next-generation sequencing technology and reported here.

RESULTS

These four cp. genomes exhibited a typical quadripartite structure and contained the large single-copy regions (LSC, 83,288 to 83,399 bp), the small single-copy regions (SSC, 18,262 to 18,287 bp), and the inverted repeat regions (IR, 24,837 to 24,842 bp). Each of the four cp. genomes encoded 134 genes, including 87 protein-coding genes, 37 tRNA genes, 8 rRNA genes, and 2 pseudogenes (ycf1 and rps19). The highly variable regions (trnC-GCA-petN, ccsA-psaC, trnE-UUC-rpoB, ycf1, ccsA and petN) may be used as potential molecular barcodes. The complete cp. genomes sequence of Synotis could be used as the potentially effective super-barcode to accurately identify Synotis species. Phylogenetic analysis demonstrated that the four Synotis species were clustered into a monophyletic group, and they were closed to the Senecio, Crassocephalum and Dendrosenecio in tribe Senecioneae.

CONCLUSIONS

This study will be useful for further species identification, evolution, genetic diversity and phylogenetic studies within this genus Synotis and the tribe Senecioneae.

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.

The pan-plastome of Hemerocallis citrina reveals new insights into the genetic diversity and cultivation history of an economically important food plant

BACKGROUND

Hemerocallis citrina Baroni (Huang hua cai in Chinese) is a perennial herbaceous plant grown for its flower buds that are eaten fresh or dried and is known as the vegetarian three treasures. The nuclear genome of H. citrina has been reported, but the intraspecific variation of the plastome (plastid genome) has not yet been studied. Therefore, the panplastome of this species collected from diverse locations is reported here for the first time.

RESULTS

In this study, 65 H. citrina samples were resequenced, de novo assembled, and aligned with the published plastome of H. citrina to resolve the H. citrina panplastome. The sizes of the 65 newly assembled complete plastomes of H. citrina ranged from 156,048 bp to 156,263 bp, and the total GC content ranged from 37.31 to 37.34%. The structure of the complete plastomes showed a typical tetrameric structure, including a large single copy (LSC), a small single copy (SSC), and a pair of inverted repeat regions (IRA and IRB). Many nucleotide variants were identified between plastomes, among which the variants in the intergenic spacer region were the most abundant, with the highest number of variants concentrated in the LSC region. Based on the phylogenetic tree constructed using the ML method, population structure analysis, and principal component analysis (PCA), the panplastome data were subdivided into five genetic clusters. The C5 genetic cluster was mostly represented by samples from Qidong, Hunan Province, while samples from Shanxi and Shaanxi Provinces were classified into the C4 genetic cluster. The greatest genetic diversity was found in the C1 genetic cluster, and the greatest genetic distance between any two clusters was found between the C4 and C5 clusters.

CONCLUSION

The resolution of the panplastome and the analysis of the population structure of H. citrina plastomes provide important data for future breeding projects and germplasm preservation.

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

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

Complete chloroplast genome sequence of Rhododendronmariesii and comparative genomics of related species in the family Ericaeae

Rhododendronmariesii Hemsley et Wilson, 1907, a typical member of the family Ericaeae, possesses valuable medicinal and horticultural properties. In this research, the complete chloroplast (cp) genome of R.mariesii was sequenced and assembled, which proved to be a typical quadripartite structure with the length of 203,480 bp. In particular, the lengths of the large single copy region (LSC), small single copy region (SSC), and inverted repeat regions (IR) were 113,715 bp, 7,953 bp, and 40,918 bp, respectively. Among the 151 unique genes, 98 were protein-coding genes, 8 were tRNA genes, and 45 were rRNA genes. The structural characteristics of the R.mariesiicp genome was similar to other angiosperms. Leucine was the most representative amino acid, while cysteine was the lowest representative. Totally, 30 codons showed obvious codon usage bias, and most were A/U-ending codons. Six highly variable regions were observed, such as trnK-pafI and atpE-rpoB, which could serve as potential markers for future barcoding and phylogenetic research of R.mariesii species. Coding regions were more conserved than non-coding regions. Expansion and contraction in the IR region might be the main length variation in R.mariesii and related Ericaeae species. Maximum-likelihood (ML) phylogenetic analysis revealed that R.mariesii was relatively closed to the R.simsii Planchon, 1853 and R.pulchrum Sweet,1831. This research will supply rich genetic resource for R.mariesii and related species of the Ericaeae.

Phylogeny of Leontopodium (Asteraceae) in China-with a reference to plastid genome and nuclear ribosomal DNA

The infrageneric taxonomy system, species delimitation, and interspecies systematic relationships of Leontopodium remain controversial and complex. However, only a few studies have focused on the molecular phylogeny of this genus. In this study, the characteristics of 43 chloroplast genomes of Leontopodium and its closely related genera were analyzed. Phylogenetic relationships were inferred based on chloroplast genomes and nuclear ribosomal DNA (nrDNA). Finally, together with the morphological characteristics, the relationships within Leontopodium were identified and discussed. The results showed that the chloroplast genomes of Filago, Gamochaeta, and Leontopodium were well-conserved in terms of gene number, gene order, and GC content. The most remarkable differences among the three genera were the length of the complete chloroplast genome, large single-copy region, small single-copy region, and inverted repeat region. In addition, the chloroplast genome structure of Leontopodium exhibited high consistency and was obviously different from that of Filago and Gamochaeta in some regions, such as matk, trnK (UUU)-rps16, petN-psbM, and trnE (UUC)-rpoB. All the phylogenetic trees indicated that Leontopodium was monophyletic. Except for the subgeneric level, our molecular phylogenetic results were inconsistent with the previous taxonomic system, which was based on morphological characteristics. Nevertheless, we found that the characteristics of the leaf base, stem types, and carpopodium base were phylogenetically correlated and may have potential value in the taxonomic study of Leontopodium. In the phylogenetic trees inferred using complete chloroplast genomes, the subgen. Leontopodium was divided into two clades (Clades 1 and 2), with most species in Clade 1 having herbaceous stems, amplexicaul, or sheathed leaves, and constricted carpopodium; most species in Clade 2 had woody stems, not amplexicaul and sheathed leaves, and not constricted carpopodium.

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.

The pan-plastome of tartary buckwheat (fagopyrum tataricum): key insights into genetic diversity and the history of lineage divergence

BACKGROUND

Tartary buckwheat (Fagopyrum tataricum) is an important food and medicine crop plant, which has been cultivated for 4000 years. A nuclear genome has been generated for this species, while an intraspecific pan-plastome has yet to be produced. As such a detailed understanding of the maternal genealogy of Tartary buckwheat has not been thoroughly investigated.

RESULTS

In this study, we de novo assembled 513 complete plastomes of Fagopyrum and compared with 8 complete plastomes of Fagopyrum downloaded from the NCBI database to construct a pan-plastome for F. tartaricum and resolve genomic variation. The complete plastomes of the 513 newly assembled Fagopyrum plastome sizes ranged from 159,253 bp to 159,576 bp with total GC contents ranged from 37.76 to 37.97%. These plastomes all maintained the typical quadripartite structure, consisting of a pair of inverted repeat regions (IRA and IRB) separated by a large single copy region (LSC) and a small single copy region (SSC). Although the structure and gene content of the Fagopyrum plastomes are conserved, numerous nucleotide variations were detected from which population structure could be resolved. The nucleotide variants were most abundant in the non-coding regions of the genome and of those the intergenic regions had the most. Mutational hotspots were primarily found in the LSC regions. The complete 521 Fagopyrum plastomes were divided into five genetic clusters, among which 509 Tartary buckwheat plastomes were divided into three genetic clusters (Ft-I/Ft-II/Ft-III). The genetic diversity in the Tartary buckwheat genetic clusters was the greatest in Ft-III, and the genetic distance between Ft-I and Ft-II was the largest. Based on the results of population structure and genetic diversity analysis, Ft-III was further subdivided into three subgroups Ft-IIIa, Ft-IIIb, and Ft-IIIc. Divergence time estimation indicated that the genera Fagopyrum and Rheum (rhubarb) shared a common ancestor about 48 million years ago (mya) and that intraspecies divergence in Tartary buckwheat began around 0.42 mya.

CONCLUSIONS

The resolution of pan-plastome diversity in Tartary buckwheat provides an important resource for future projects such as marker-assisted breeding and germplasm preservation.

Genetic relationship of Pleione based on the chloroplast genome

Pleione (Orchidaceae) is a very famous horticultural plant with a high international reputation for its unique flower shape and abundant color. The small difference in morphological characteristics among Pleione species caused by weak reproductive isolation and easy hybridization makes the taxonomic status of individual species very confusing. Chloroplast (cp) genome information is of great significance for the study of plant phylogeny and taxonomy. In this study, the cp genomes of Pleione were sequenced and the complete cp structure and sequence characteristics of 19 species were compared and analyzed. The cp genome of Pleione species exhibited a conserved tetrad structure and each species encoded 135 protein-coding genes, 38 tRNA and 8 RNA genes. The cp genome sizes of 19 Pleione were 157964-159269 bp and the length of LSC, SSC and IR were 85953-87172 bp, 18499-18712 bp, 26459-26756 bp, respectively. Palindromic and forward repeats accounted for a high proportion and the SSRs were mainly mononucleotide repeats in Pleione. Analysis of chloroplast sequence differences indicated that the differences in coding regions were smaller than those in non-coding regions, and the variation in LSC and SSC regions was greater than that in IR regions. Phylogenetic analysis showed that all Pleione species inferred from the cp genome were clustered together and received high support. However, the genetic relationship of Pleione plants is different from the current update system of this genus. Therefore, the demarcation of Pleione interspecific relationships still needs further investigation due to the lack of sufficient evidence. The cp genome serves as valuable information for the identification of Pleione species and the study of phylogenetic relationships.

Insights into Synonymous Codon Usage Bias in Hepatitis C Virus and Its Adaptation to Hosts

Hepatitis C virus (HCV) is enveloped RNA virus, encoding for a polyprotein that is processed by cellular proteases. The virus is responsible for liver cirrhosis, allograft rejection, and human hepatocellular carcinoma. Based on studies including compositional analysis, odds ratio analysis, parity analysis, skew analysis, relative synonymous codon usage, codon bias, and protein properties, it was evident that codon usage bias in HCV is dependent upon the nucleotide composition. Codon context analysis revealed CTC-CTG as a preferred codon pair. While CGA and CGT codons were rare, none of the codons were rare in HCV-like viruses envisaged in the present study. Many of the preferred codon pairs were valine amino acid-initiated, which possibly infers viral infectivity; hence the role of selection forces appears to act on the HCV genome, which was further validated by neutrality analysis where selection accounted for 87.28%, while mutation accounted for 12.72% force shaping codon usage. Furthermore, codon usage was correlated with the length of the genome. HCV viruses prefer valine-initiated codon pairs, while HCV-like viruses prefer alanine-initiated codon pairs. The HCV host range is very narrow and is confined to only humans and chimpanzees. Based on indices including codon usage correlation analysis, similarity index, and relative codon deoptimization index, it is evident in the study that the chimpanzee is the primary host of the virus. The present study helped elucidate the preferred host for HCV. The information presented in the study paved the way for generating an attenuated vaccine candidate through viral recoding, with finely tuned nucleotide composition and a perfect balance of preferred and rare codons.

The complete chloroplast genome of an endangered plant Artemisia borotalensis (Asteraceae) and phylogenetic analysis

Artemisia borotalensis Poljakov is an endemic and endangered herb in China. In this study, we sequenced and analyzed the complete chloroplast genome of this species. Sequencing revealed the genome to be 151,179 bp in length, containing a large single copy region (82,862 bp), a small single copy region (18,377 bp), and a pair of inverted repeat regions (24,970 bp each). Our analyses demonstrated that it contained 133 genes, including 87 protein-coding genes, 37 transfer RNA genes, eight ribosomal RNA genes, and one pseudogene (ycf1). Furthermore, we found the genome to have an overall GC content of 37.4%. A phylogenetic analysis indicated that A. borotalensis and A. maritima clustered together as sister group to A. annua and A. fukudo clade.

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.

Chloroplast phylogenomics and the taxonomy of Saxifraga section Ciliatae (Saxifragaceae)

Comprising ca. 200 species, Saxifraga sect. Ciliatae is the most species-rich section of Saxifraga s.str., whose center of diversity is in the Tibeto-Himalayan region. The infra-sectional classification of sect. Ciliatae is still in debate due to the high level of species richness, as well as remarkable variations of habitat, morphology, physiology and life cycles. Subdivisions of this section proposed in various taxonomic systems have not been adequately tested in previous phylogenetic studies, partly due to low taxonomic sampling density, but also to the use of few DNA markers. In order to achieve a more robust infra-sectional classification of sect. Ciliatae, complete chloroplast genomes of 94 taxa from this section were analyzed, of which 93 were newly sequenced, assembled and annotated. The length of the 94 plastomes of sect. Ciliatae taxa range from 143,479 to 159,938 bp, encoding 75 to 79 unique protein-coding genes (PCGs). Analyses of the 94 plastomes revealed high conservation in structural organization, gene arrangement, and gene content. Gene loss and changes of IR boundaries were detected but in extremely low frequency. The molecular phylogenetic tree from concatenated PCGs and complete chloroplast genome sequences exhibited high resolution and support values and confirms that sect. Ciliatae is monophyletic. Three well-supported clades were revealed within the section that agree relatively well with the subsectional taxonomy of Gornall (1987), but some minor modifications should be made. Firstly, the monotypic subsection Cinerascentes should be abandoned and its constituent species, S. cinerascens, assigned to subsect. Gemmiparae. Secondly, subsections Rosulares and Serpyllifoliae should be merged and become subsect. Rosulares. Section Ciliatae thus comprises: subsect. Hirculoideae Engl. & Irmsch.; subsect. Rosulares Gornall; subsect. Gemmiparae Engl. & Irmsch.; subsect. Flagellares (C. B. Clarke) Engl. & Irmsch. and subsect. Hemisphaericae (Engl. & Irmsch.) Gornall.

The Chloroplast Genome of Wild Saposhnikovia divaricata: Genomic Features, Comparative Analysis, and Phylogenetic Relationships

Saposhnikovia divaricata, a well-known Chinese medicinal herb, is the sole species under the genus Saposhnikovia of the Apiaceae subfamily Apioideae Drude. However, information regarding its genetic diversity and evolution is still limited. In this study, the first complete chloroplast genome (cpDNA) of wild S. divaricata was generated using de novo sequencing technology. Similar to the characteristics of Ledebouriella seseloides, the 147,834 bp-long S. divaricata cpDNA contained a large single copy, a small single copy, and two inverted repeat regions. A total of 85 protein-coding, 8 ribosomal RNA, and 36 transfer RNA genes were identified. Compared with five other species, the non-coding regions in the S. divaricata cpDNA exhibited greater variation than the coding regions. Several repeat sequences were also discovered, namely, 33 forward, 14 reverse, 3 complement, and 49 microsatellite repeats. Furthermore, phylogenetic analysis using 47 cpDNA sequences of Apioideae members revealed that L. seseloides and S. divaricata clustered together with a 100% bootstrap value, thereby supporting the validity of renaming L. seseloides to S. divaricata at the genomic level. Notably, S. divaricata was most closely related to Libanotis buchtormensis, which contradicts previous reports. Therefore, these findings provide a valuable foundation for future studies on the genetic diversity and evolution of S. divaricata.

Genetic and chemical markers for authentication of three Artemisia species: A. capillaris, A. gmelinii, and A. fukudo

The genus Artemisia is an important source of medicines in both traditional and modern pharmaceutics, particularly in East Asia. Despite the great benefits of herbal medicine, quality assessment methods for these medicinal herbs are lacking. The young leaves from Artemisia species are generally used, and most of the species have similar morphology, which often leads to adulteration and misuse. This study assembled five complete chloroplast genomes of three Artemisia species, two accessions of A. gmelinii and A. capillaris, and one A. fukudo. Through comparative analysis, we revealed genomic variations and phylogenetic relationships between these species and developed seven InDel-based barcode markers which discriminated the tested species from each other. Additionally, we analyzed specialized metabolites from the species using LC-MS and suggested chemical markers for the identification and authentication of these herbs. We expect that this integrated and complementary authentication method would aid in reducing the misuse of Artemisia species.

Identification of Medicinal Bidens Plants for Quality Control Based on Organelle Genomes

Bidens plants are annuals or perennials of Asteraceae and usually used as medicinal materials in China. They are difficult to identify by using traditional identification methods because they have similar morphologies and chemical components. Universal DNA barcodes also cannot identify Bidens species effectively. This situation seriously hinders the development of medicinal Bidens plants. Therefore, developing an accurate and effective method for identifying medicinal Bidens plants is urgently needed. The present study aims to use phylogenomic approaches based on organelle genomes to address the confusing relationships of medicinal Bidens plants. Illumina sequencing was used to sequence 12 chloroplast and eight mitochondrial genomes of five species and one variety of Bidens. The complete organelle genomes were assembled, annotated and analysed. Phylogenetic trees were constructed on the basis of the organelle genomes and highly variable regions. The organelle genomes of these Bidens species had a conserved gene content and codon usage. The 12 chloroplast genomes of the Bidens species were 150,489 bp to 151,635 bp in length. The lengths of the eight mitochondrial genomes varied from each other. Bioinformatics analysis revealed the presence of 50–71 simple sequence repeats and 46–181 long repeats in the organelle genomes. By combining the results of mVISTA and nucleotide diversity analyses, seven candidate highly variable regions in the chloroplast genomes were screened for species identification and relationship studies. Comparison with the complete mitochondrial genomes and common protein-coding genes shared by each organelle genome revealed that the complete chloroplast genomes had the highest discriminatory power for Bidens species and thus could be used as a super barcode to authenticate Bidens species accurately. In addition, the screened highly variable region trnS-GGA-rps4 could be also used as a potential specific barcode to identify Bidens species.

Evolutionary patterns of nucleotide substitution rates in plastid genomes of Quercus

Molecular evolution, including nucleotide substitutions, plays an important role in understanding the dynamics and mechanisms of species evolution. Here, we sequenced whole plastid genomes (plastomes) of Quercus fabri, Quercus semecarpifolia, Quercus engleriana, and Quercus phellos and compared them with 14 other Quercus plastomes to explore their evolutionary relationships using 67 shared protein-coding sequences. While many previously identified evolutionary relationships were found, our findings do not support previous research which retrieve Quercus subg. Cerris sect. Ilex as a monophyletic group, with sect. Ilex found to be polyphyletic and composed of three strongly supported lineages inserted between sections Cerris and Cyclobalanposis. Compared with gymnosperms, Quercus plastomes showed higher evolutionary rates (Dn/Ds = 0.3793). Most protein-coding genes experienced relaxed purifying selection, and the high Dn value (0.1927) indicated that gene functions adjusted to environmental changes effectively. Our findings suggest that gene interval regions play an important role in Quercus evolution. We detected greater variation in the intergenic regions (trnH-psbA, trnK_UUU-rps16, trnfM_CAU-rps14, trnS_GCU-trnG_GCC, and atpF-atpH), intron losses (petB and petD), and pseudogene loss and degradation (ycf15). Additionally, the loss of some genes suggested the existence of gene exchanges between plastid and nuclear genomes, which affects the evolutionary rate of the former. However, the connective mechanism between these two genomes is still unclear.

Complete chloroplast genome of the medicinal plant Evolvulus alsinoides: comparative analysis, identification of mutational hotspots and evolutionary dynamics with species of Solanales

Evolvulus alsinoides, belonging to the family Convolvulaceae, is an important medicinal plant widely used as a nootropic in the Indian traditional medicine system. In the genus Evolvulus, no research on the chloroplast genome has been published. Hence, the present study focuses on annotation, characterization, identification of mutational hotspots, and phylogenetic analysis in the complete chloroplast genome (cp) of E. alsinoides. Genome comparison and evolutionary dynamics were performed with the species of Solanales. The cp genome has 114 genes (80 protein-coding genes, 30 transfer RNA, and 4 ribosomal RNA genes) that were unique with total genome size of 157,015 bp. The cp genome possesses 69 RNA editing sites and 44 simple sequence repeats (SSRs). Predicted SSRs were randomly selected and validated experimentally. Six divergent hotspots such as trnQ-UUG, trnF-GAA, psaI, clpP, ndhF, and ycf1 were discovered from the cp genome. These microsatellites and divergent hot spot sequences of the Taxa 'Evolvulus' could be employed as molecular markers for species identification and genetic divergence investigations. The LSC area was found to be more conserved than the SSC and IR region in genome comparison. The IR contraction and expansion studies show that nine genes rpl2, rpl23, ycf1, ycf2, ycf1, ndhF, ndhA, matK, and psbK were present in the IR-LSC and IR-SSC boundaries of the cp genome. Fifty-four protein-coding genes in the cp genome were under negative selection pressure, indicating that they were well conserved and were undergoing purifying selection. The phylogenetic analysis reveals that E. alsinoides is closely related to the genus Cressa with some divergence from the genus Ipomoea. This is the first time the chloroplast genome of the genus Evolvulus has been published. The findings of the present study and chloroplast genome data could be a valuable resource for future studies in population genetics, genetic diversity, and evolutionary relationship of the family Convolvulaceae.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01051-w.

Towards the Well-Tempered Chloroplast DNA Sequences

With the development of next-generation sequencing technology and bioinformatics tools, the process of assembling DNA sequences has become cheaper and easier, especially in the case of much shorter organelle genomes. The number of available DNA sequences of complete chloroplast genomes in public genetic databases is constantly increasing and the data are widely used in plant phylogenetic and biotechnological research. In this work, we investigated possible inconsistencies in the stored form of publicly available chloroplast genome sequence data. The impact of these inconsistencies on the results of the phylogenetic analysis was investigated and the bioinformatic solution to identify and correct inconsistencies was implemented. The whole procedure was demonstrated using five plant families (Apiaceae, Asteraceae, Campanulaceae, Lamiaceae and Rosaceae) as examples.

Identification and phylogenetic analysis of five Crataegus species (Rosaceae) based on complete chloroplast genomes

The chloroplast genomes of the five Crataegus species were shown to have a conserved genome structure. Complete chloroplast genome sequences were more suitable than highly variable regions for the identification and phylogenetic analysis of Crataegus species. Hawthorn, which is commonly used as a traditional Chinese medicine, is one of the most popular sour fruits and has high economic value. Crataegus pinnatifida var. pinnatifida and C. pinnatifida var. major are frequently adulterated with other Crataegus species on the herbal medicine market. However, most Crataegus plants are difficult to identify using traditional morphological methods. Here, we compared five Crataegus chloroplast (CP) genomes comprising two newly sequenced (i.e., C. pinnatifida var. pinnatifida and C. pinnatifida var. major) and three previously published CP genomes. The CP genomes of the five Crataegus species had a conserved genome structure, gene content and codon usage. The total length of the CP genomes was 159,654-159,865 bp. A total of 129-130 genes, including 84-85 protein-coding genes, 37 tRNA genes and 8 rRNA genes, were annotated. Bioinformatics analysis revealed 96-103 simple sequence repeats (SSRs) and 48-70 long repeats in the five CP genomes. Combining the results of mVISTA and nucleotide diversity, five highly variable regions were screened for species identification and relationship studies. Maximum likelihood trees were constructed on the basis of complete CP genome sequences and highly variable regions. The results showed that the former had higher discriminatory power for Crataegus species, indicating that the complete CP genome could be used as a super-barcode to accurately authenticate the five Crataegus species.

Comparative genome analysis revealed gene inversions, boundary expansions and contractions, and gene loss in the Stemona sessilifolia (Miq.) Miq. chloroplast genome

Stemona sessilifolia (Miq.) Miq., commonly known as Baibu, is one of the most popular herbal medicines in Asia. In the Chinese Pharmacopoeia, Baibu has multiple authentic sources and there are many similar herbs sold as Baibu in herbal medicine markets. The existence of counterfeits of Baibu brings challenges to its identification. To assist in its accurate identification, we sequenced and analyzed the complete chloroplast genome of S. sessilifolia using next-generation sequencing technology. The genome was found to be 154,037 bp in length, possessing a typical quadripartite structure consisting of a pair of inverted repeats (IRs: 27,090 bp) separated by a large single copy (LSC: 81,949 bp) and a small single copy (SSC: 17,908 bp). A total of 112 unique genes were identified, including 80 protein-coding, 28 transfer RNA and four ribosomal RNA genes. In addition, 45 tandem, 27 forward, 23 palindromic and 104 simple sequence repeats were detected in the genome by repeated analysis. Compared with its counterfeits (Asparagus officinalis and Carludovica palmata) we found that IR expansion and SSC contraction events of S. sessilifolia resulted in two copies of the rpl22 gene in the IR regions and a partial duplication of the ndhF gene in the SSC region. An approximately 3-kb-long inversion was also identified in the LSC region, leading to the petA and cemA genes being presented in the complementary strand of the chloroplast DNA molecule. Comparative analysis revealed some highly variable regions, including trnF-GAA_ndhJ, atpB_rbcL, rps15_ycf1, trnG-UCC_trnR-UCU, ndhF_rpl32, accD_psaI, rps2_rpoC2, trnS-GCU_trnG-UCC, trnT-UGU_trnL-UAA and rps16_trnQ-UUG. Finally, gene loss events were investigated in the context of phylogenetic relationships. In summary, the complete plastome of S. sessilifolia will provide valuable information for the distinction between Baibu and its counterfeits and assist in elucidating the evolution of S. sessilifolia.

Characteristics of the completed chloroplast genome sequence of Xanthium spinosum: comparative analyses, identification of mutational hotspots and phylogenetic implications

BACKGROUND

The invasive species Xanthium spinosum has been used as a traditional Chinese medicine for many years. Unfortunately, no extensive molecular studies of this plant have been conducted.

RESULTS

Here, the complete chloroplast (cp) genome sequence of X. spinosum was assembled and analyzed. The cp genome of X. spinosum was 152,422 base pairs (bp) in length, with a quadripartite circular structure. The cp genome contained 115 unique genes, including 80 PCGs, 31 tRNA genes, and 4 rRNA genes. Comparative analyses revealed that X. spinosum contains a large number of repeats (999 repeats) and 701 SSRs in its cp genome. Fourteen divergences (Π > 0.03) were found in the intergenic spacer regions. Phylogenetic analyses revealed that Parthenium is a sister clade to both Xanthium and Ambrosia and an early-diverging lineage of subtribe Ambrosiinae, although this finding was supported with a very weak bootstrap value.

CONCLUSION

The identified hotspot regions could be used as molecular markers for resolving phylogenetic relationships and species identification in the genus Xanthium.

The Complete Chloroplast Genome of the Vulnerable Oreocharis esquirolii (Gesneriaceae): Structural Features, Comparative and Phylogenetic Analysis

Oreocharis esquirolii, a member of Gesneriaceae, is known as Thamnocharis esquirolii, which has been regarded a synonym of the former. The species is endemic to Guizhou, southwestern China, and is evaluated as vulnerable (VU) under the International Union for Conservation of Nature (IUCN) criteria. Until now, the sequence and genome information of O. esquirolii remains unknown. In this study, we assembled and characterized the complete chloroplast (cp) genome of O. esquirolii using Illumina sequencing data for the first time. The total length of the cp genome was 154,069 bp with a typical quadripartite structure consisting of a pair of inverted repeats (IRs) of 25,392 bp separated by a large single copy region (LSC) of 85,156 bp and a small single copy region (SSC) of18,129 bp. The genome comprised 114 unique genes with 80 protein-coding genes, 30 tRNA genes, and four rRNA genes. Thirty-one repeat sequences and 74 simple sequence repeats (SSRs) were identified. Genome alignment across five plastid genomes of Gesneriaceae indicated a high sequence similarity. Four highly variable sites (rps16-trnQ, trnS-trnG, ndhF-rpl32, and ycf 1) were identified. Phylogenetic analysis indicated that O. esquirolii grouped together with O. mileensis, supporting resurrection of the name Oreocharis esquirolii from Thamnocharisesquirolii. The complete cp genome sequence will contribute to further studies in molecular identification, genetic diversity, and phylogeny.

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.

Plastome Evolution in Dolomiaea (Asteraceae, Cardueae) Using Phylogenomic and Comparative Analyses

Dolomiaea is a medicinally important genus of Asteraceae endemic to alpine habitats of the Qinghai-Tibet Plateau (QTP) and adjacent areas. Despite significant medicinal value, genomic resources of Dolomiaea are still lacking, impeding our understanding of its evolutionary history. Here, we sequenced and annotated plastomes of four Dolomiaea species. All analyzed plastomes share the gene content and structure of most Asteraceae plastomes, indicating the conservation of plastome evolutionary history of Dolomiaea. Eight highly divergent regions (rps16-trnQ, trnC-petN, trnE-rpoB, trnT-trnL-trnF, psbE-petL, ndhF-rpl32-trnL, rps15-ycf1, and ycf1), along with a total of 51-61 simple sequence repeats (SSRs) were identified as valuable molecular markers for further species delimitation and population genetic studies. Phylogenetic analyses confirmed the evolutionary position of Dolomiaea as a clade within the subtribe Saussureinae, while revealing the discordance between the molecular phylogeny and morphological treatment. Our analysis also revealed that the plastid genes, rpoC2 and ycf1, which are rarely used in Asteraceae phylogenetic inference, exhibit great phylogenetic informativeness and promise in further phylogenetic studies of tribe Cardueae. Analysis for signatures of selection identified four genes that contain sites undergoing positive selection (atpA, ndhF, rbcL, and ycf4). These genes may play important roles in the adaptation of Dolomiaea to alpine environments. Our study constitutes the first investigation on the sequence and structural variation, phylogenetic utility and positive selection of plastomes of Dolomiaea, which will facilitate further studies of its taxonomy, evolution and conservation.

Can we use it? On the utility of de novo and reference-based assembly of Nanopore data for plant plastome sequencing

The chloroplast genome harbors plenty of valuable information for phylogenetic research. Illumina short-read data is generally used for de novo assembly of whole plastomes. PacBio or Oxford Nanopore long reads are additionally employed in hybrid approaches to enable assembly across the highly similar inverted repeats of a chloroplast genome. Unlike for PacBio, plastome assemblies based solely on Nanopore reads are rarely found, due to their high error rate and non-random error profile. However, the actual quality decline connected to their use has rarely been quantified. Furthermore, no study has employed reference-based assembly using Nanopore reads, which is common with Illumina data. Using Leucanthemum Mill. as an example, we compared the sequence quality of seven chloroplast genome assemblies of the same species, using combinations of two sequencing platforms and three analysis pipelines. In addition, we assessed the factors which might influence Nanopore assembly quality during sequence generation and bioinformatic processing. The consensus sequence derived from de novo assembly of Nanopore data had a sequence identity of 99.59% compared to Illumina short-read de novo assembly. Most of the errors detected were indels (81.5%), and a large majority of them is part of homopolymer regions. The quality of reference-based assembly is heavily dependent upon the choice of a close-enough reference. When using a reference with 0.83% sequence divergence from the studied species, mapping of Nanopore reads results in a consensus comparable to that from Nanopore de novo assembly, and of only slightly inferior quality compared to a reference-based assembly with Illumina data. For optimal de novo assembly of Nanopore data, appropriate filtering of contaminants and chimeric sequences, as well as employing moderate read coverage, is essential. Based on these results, we conclude that Nanopore long reads are a suitable alternative to Illumina short reads in plastome phylogenomics. Few errors remain in the finalized assembly, which can be easily masked in phylogenetic analyses without loss in analytical accuracy. The easily applicable and cost-effective technology might warrant more attention by researchers dealing with plant chloroplast genomes.

Long-Reads Reveal That the Chloroplast Genome Exists in Two Distinct Versions in Most Plants

The chloroplast genome usually has a quadripartite structure consisting of a large single copy region and a small single copy region separated by two long inverted repeats. It has been known for some time that a single cell may contain at least two structural haplotypes of this structure, which differ in the relative orientation of the single copy regions. However, the methods required to detect and measure the abundance of the structural haplotypes are labor-intensive, and this phenomenon remains understudied. Here, we develop a new method, Cp-hap, to detect all possible structural haplotypes of chloroplast genomes of quadripartite structure using long-read sequencing data. We use this method to conduct a systematic analysis and quantification of chloroplast structural haplotypes in 61 land plant species across 19 orders of Angiosperms, Gymnosperms, and Pteridophytes. Our results show that there are two chloroplast structural haplotypes which occur with equal frequency in most land plant individuals. Nevertheless, species whose chloroplast genomes lack inverted repeats or have short inverted repeats have just a single structural haplotype. We also show that the relative abundance of the two structural haplotypes remains constant across multiple samples from a single individual plant, suggesting that the process which maintains equal frequency of the two haplotypes operates rapidly, consistent with the hypothesis that flip-flop recombination mediates chloroplast structural heteroplasmy. Our results suggest that previous claims of differences in chloroplast genome structure between species may need to be revisited.

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.

Complete chloroplast genome sequence of the medicinal plant Arctium lappa

Arctium lappa, commonly called burdock, has a long medicinal and edible history. It has recently gained increasing attention because of its economic value. In this study, we obtained the complete chloroplast genome of A. lappa by Illumina Hiseq. The complete chloroplast genome of A. lappa is a typical circular structure with 152 708 bp in length. The GC content in the whole chloroplast genome of A. lappa is 37.7%. A total of 37 tRNA genes, 8 rRNA genes, and 87 protein-coding genes were successfully annotated. And the chloroplast genome contains 113 unique genes, 19 of which are duplicated in the inverted repeat. The distribution of 39 simple sequence repeats was analysed, and most of them are in the large single-copy (LSC) sequence. An inversion comprising 16 genes was found in the LSC region, which is 26 283 bp long. We performed multiple sequence alignments using 72 common protein-coding genes of 29 species and constructed a Maximum Parsimony (MP) tree. The MP phylogenetic result shows that A. lappa grouped together with Carthamus tinctorius, Centaurea diffusa, and Saussurea involucrata. The chloroplast genome of A. lappa is a valuable resource for further studies in Asteraceae.

Characterization of the complete chloroplast genome of Malus baccata var. xiaojinensis

Malus baccata var. xiaojinensis belonging to Rosaceae is not only an ornamental tree but also an apple genotype that is highly tolerant to Fe deficiency. In this study, we reported the complete chloroplast (cp) genome of M. baccata var. xiaojinensis using ILLUMINA sequencing. The whole cp genome is 160,067 bp in length, containing a pair of inverted repeats (IRs) of 26,358 bp, a large single copy (LSC) region of 88,157 bp and a small single copy (SSC) region of 19,194 bp. And, the overall GC content of the cp genome was 36.56%. A total of 110 unique genes were found in the cp genome. 17 genes were duplicated in the IRs, including six protein-coding genes, seven tRNA genes, and four rRNA genes. Fourteen genes contained one intron, whereas three genes (rps12, clpP, and ycf3) contained two introns. The phylogenetic analysis demonstrated a close relationship between M. baccata var. xiaojinensis and Malus hupehensis.

Plastome phylogenomics of Saussurea (Asteraceae: Cardueae)

BACKGROUND

Saussurea DC. is one of the largest and most morphologically heterogeneous genera in Asteraceae. The relationships within Saussurea have been poorly resolved, probably due an early, rapid radiation. To examine plastome evolution and resolve backbone relationships within Saussurea, we sequenced the complete plastomes of 17 species representing all four subgenera.

RESULTS

All Saussurea plastomes shared the gene content and structure of most Asteraceae plastomes. Molecular evolutionary analysis showed most of the plastid protein-coding genes have been under purifying selection. Phylogenomic analyses of 20 Saussurea plastomes that alternatively included nucleotide or amino acid sequences of all protein-coding genes, vs. the nucleotide sequence of the entire plastome, supported the monophyly of Saussurea and identified three clades within it. Three of the four traditional subgenera were recovered as paraphyletic. Seven plastome regions were identified as containing the highest nucleotide variability.

CONCLUSIONS

Our analyses reveal both the structural conservatism and power of the plastome for resolving relationships in congeneric taxa. It is very likely that differences in topology among data sets is due primarily to differences in numbers of parsimony-informative characters. Our study demonstrates that the current taxonomy of Saussurea is likely based at least partly on convergent morphological character states. Greater taxon sampling will be necessary to explore character evolution and biogeography in the genus. Our results here provide helpful insight into which loci will provide the most phylogenetic signal in Saussurea and Cardueae.

Molecular phylogenomics of the tribe Shoreeae (Dipterocarpaceae) using whole plastid genomes

BACKGROUND AND AIMS

Phylogenetic relationships within tribe Shoreeae, containing the main elements of tropical forests in Southeast Asia, present a long-standing problem in the systematics of Dipterocarpaceae. Sequencing whole plastomes using next-generation sequencing- (NGS) based genome skimming is increasingly employed for investigating phylogenetic relationships of plants. Here, the usefulness of complete plastid genome sequences in resolving phylogenetic relationships within Shoreeae is evaluated.

METHODS

A pipeline to obtain alignments of whole plastid genome sequences across individuals with different amounts of available data is presented. In total, 48 individuals, representing 37 species and four genera of the ecologically and economically important tribe Shoreeae sensu Ashton, were investigated. Phylogenetic trees were reconstructed using maximum parsimony, maximum likelihood and Bayesian inference.

KEY RESULTS

Here, the first fully sequenced plastid genomes for the tribe Shoreeae are presented. Their size, GC content and gene order are comparable with those of other members of Malvales. Phylogenomic analyses demonstrate that whole plastid genomes are useful for inferring phylogenetic relationships among genera and groups of Shorea (Shoreeae) but fail to provide well-supported phylogenetic relationships among some of the most closely related species. Discordance in placement of Parashorea was observed between phylogenetic trees obtained from plastome analyses and those obtained from nuclear single nucleotide polymorphism (SNP) data sets identified in restriction-site associated sequencing (RADseq).

CONCLUSIONS

Phylogenomic analyses of the entire plastid genomes are useful for inferring phylogenetic relationships at lower taxonomic levels, but are not sufficient for detailed phylogenetic reconstructions of closely related species groups in Shoreeae. Discordance in placement of Parashorea was further investigated for evidence of ancient hybridization.

Sequencing the Plastid Genome of Giant Ragweed (Ambrosia trifida, Asteraceae) From a Herbarium Specimen

We report the first plastome sequence of giant ragweed (Ambrosia trifida); with this new genome information, we assessed the phylogeny of Asteraceae and the transcriptional profiling against glyphosate resistance in giant ragweed. Assembly and genic features show a normal angiosperm quadripartite plastome structure with no signatures of deviation in gene directionality. Comparative analysis revealed large inversions across the plastome of giant ragweed and the previously sequenced members of the plant family. Asteraceae plastid genomes contain two inversions of 22.8 and 3.3 kb; the former is located between trnS-GCU and trnG-UCC genes, and the latter between trnE-UUC and trnT-GGU genes. The plastid genome sequences of A. trifida and the related species, Ambrosia artemisiifolia, are identical in gene content and arrangement, but they differ in length. The phylogeny is well-resolved and congruent with previous hypotheses about the phylogenetic relationship of Asteraceae. Transcriptomic analysis revealed divergence in the relative expressions at the exonic and intronic levels, providing hints toward the ecological adaptation of the genus. Giant ragweed shows various levels of glyphosate resistance, with introns displaying higher expression patterns at resistant time points after the assumed herbicide treatment.

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

Artemisia selengenesis is not only a health food, but also a well-known traditional Chinese medicine. Only a fraction of the chloroplast (cp) genome data of Artemisia has been reported and chloroplast genomic materials have been widely used in genomic evolution studies, molecular marker development, and phylogenetic analysis of the genus Artemisia, which makes evolutionary studies, genetic improvement, and phylogenetic identification very difficult. In this study, the complete chloroplast genome of A. selengensis was compared with that of other species within Artemisia and phylogenetic analyses was conducted with other genera in the Asteraceae family. The results showed that A. selengensis is an AT-rich species and has a typical quadripartite structure that is 151,215 bp in length. Comparative genome analyses demonstrated that the available chloroplast genomes of species of Artemisia were well conserved in terms of genomic length, GC contents, and gene organization and order. However, some differences, which may indicate evolutionary events, were found, such as a re-inversion event within the Artemisia genus, an unequal duplicate phenomenon of the ycf1 gene because of the expansion and contraction of the IR region, and the fast-evolving regions. Repeated sequences analysis showed that Artemisia chloroplast genomes presented a highly similar pattern of SSR or LDR distribution. A total of 257 SSRs and 42 LDRs were identified in the A. selengensis chloroplast genome. The phylogenetic analysis showed that A. selengensis was sister to A. gmelinii. The findings of this study will be valuable in further studies to understand the genetic diversity and evolutionary history of Asteraceae.

Comparative Chloroplast Genomics at Low Taxonomic Levels: A Case Study Using Amphilophium (Bignonieae, Bignoniaceae)

Chloroplast (cp) genome organization, gene order, and content have long been considered conserved among land plants. Despite that, the generation of thousands of complete plastomes through next-generation sequencing (NGS) has challenged their conserved nature. In this study, we analyze 11 new complete plastomes of Amphilophium (Bignonieae, Bignoniaceae), a diverse genus of Neotropical lianas, and that of Anemopaegma prostratum. We explored the structure and content of the assembled plastomes and performed comparative analyses within Amphilophium and among other plastomes available for Bignoniaceae. The overall gene content and orientation of plastomes is similar in all species studied. Plastomes are not conserved among Amphilophium, showing significant differences in length (155,262-164,786 bp), number of genes duplicated in the IRs (eight, 18, or 19), and location of the SC/IR boundaries (i.e., LSC/IRa junction between rps19 and rpl2 genes, within petD, or within petB). Length differences reflect expansions of the IRs and contractions of the LSC regions. The plastome of A. prostratum is 168,172 bp, includes 19 duplicated genes, and has the LSC/IRa boundary located within the petB gene. Amphilophium plastomes show high nucleotide diversity, with many hypervariable regions, and 16 genes with signatures of positive selection. Multiple SSRs and repeat regions were identified for Amphilophium and Anemopaegma prostratum. The differences in structure detected within Amphilophium plastomes in terms of LSC/IR and IR/SSC boundaries, number of duplicated genes, and genome sizes are mostly shared between taxa that belong to the same clade. Our results bring new insights into the evolution of plastomes at low taxonomic levels.

Complete Chloroplast Genome Sequence and Phylogenetic Analysis of Aster tataricus

We sequenced and analyzed the complete chloroplast genome of Aster tataricus (family Asteraceae), a Chinese herb used medicinally to relieve coughs and reduce sputum. The A. tataricus chloroplast genome was 152,992 bp in size, and harbored a pair of inverted repeat regions (IRa and IRb, each 24,850 bp) divided into a large single-copy (LSC, 84,698 bp) and a small single-copy (SSC, 18,250 bp) region. Our annotation revealed that the A. tataricus chloroplast genome contained 115 genes, including 81 protein-coding genes, 4 ribosomal RNA genes, and 30 transfer RNA genes. In addition, 70 simple sequence repeats (SSRs) were detected in the A. tataricus chloroplast genome, including mononucleotides (36), dinucleotides (1), trinucleotides (23), tetranucleotides (1), pentanucleotides (8), and hexanucleotides (1). Comparative chloroplast genome analysis of three Aster species indicated that a higher similarity was preserved in the IR regions than in the LSC and SSC regions, and that the differences in the degree of preservation were slighter between A. tataricus and A. altaicus than between A. tataricus and A. spathulifolius. Phylogenetic analysis revealed that A. tataricus was more closely related to A. altaicus than to A. spathulifolius. Our findings offer valuable information for future research on Aster species identification and selective breeding.

Phylogenomic and Comparative Analyses of Complete Plastomes of Croomia and Stemona (Stemonaceae)

The monocot genus Croomia (Stemonaceae) comprises three herbaceous perennial species that exhibit EA (Eastern Asian)–ENA (Eastern North American) disjunct distribution. However, due to the lack of effective genomic resources, its evolutionary history is still weakly resolved. In the present study, we conducted comparative analysis of the complete chloroplast (cp) genomes of three Croomia species and two Stemona species. These five cp genomes proved highly similar in overall size (154,407–155,261 bp), structure, gene order and content. All five cp genomes contained the same 114 unique genes consisting of 80 protein-coding genes, 30 tRNA genes and 4 rRNA genes. Gene content, gene order, AT content and IR/SC boundary structures were almost the same among the five Stemonaceae cp genomes, except that the Stemona cp genome was found to contain an inversion in cemA and petA. The lengths of five genomes varied due to contraction/expansion of the IR/SC borders. A/T mononucleotides were the richest Simple Sequence Repeats (SSRs). A total of 46, 48, 47, 61 and 60 repeats were identified in C. japonica, C. heterosepala, C. pauciflora, S. japonica and S. mairei, respectively. A comparison of pairwise sequence divergence values across all introns and intergenic spacers revealed that the ndhF–rpl32, psbM–trnD and trnS–trnG regions are the fastest-evolving regions. These regions are therefore likely to be the best choices for molecular evolutionary and systematic studies at low taxonomic levels in Stemonaceae. Phylogenetic analyses of the complete cp genomes and 78 protein-coding genes strongly supported the monophyly of Croomia. Two Asian species were identified as sisters that likely diverged in the Early Pleistocene (1.62 Mya, 95% HPD: 1.125–2.251 Mya), whereas the divergence of C. pauciflora dated back to the Late Miocene (4.77 Mya, 95% HPD: 3.626–6.162 Mya). The availability of these cp genomes will provide valuable genetic resources for further population genetics and phylogeographic studies on Croomia.

Identification of Ligularia Herbs Using the Complete Chloroplast Genome as a Super-Barcode

More than 30 Ligularia Cass. (Asteraceae) species have long been used in folk medicine in China. Morphological features and common DNA regions are both not ideal to identify Ligularia species. As some Ligularia species contain pyrrolizidine alkaloids, which are hazardous to human and animal health and are involved in metabolic toxification in the liver, it is important to find a better way to distinguish these species. Here, we report complete chloroplast (CP) genomes of six Ligularia species, L. intermedia, L. jaluensis, L. mongolica, L. hodgsonii, L. veitchiana, and L. fischeri, obtained through high-throughput Illumina sequencing technology. These CP genomes showed typical circular tetramerous structure and their sizes range from 151,118 to 151,253 bp. The GC content of each CP genome is 37.5%. Every CP genome contains 134 genes, including 87 protein-coding genes, 37 tRNA genes, eight rRNA genes, and two pseudogenes (ycf1 and rps19). From the mVISTA, there were no potential coding or non-coding regions to distinguish these six Ligularia species, but the maximum likelihood tree of the six Ligularia species and other related species showed that the whole CP genome can be used as a super-barcode to identify these six Ligularia species. This study provides invaluable data for species identification, allowing for future studies on phylogenetic evolution and safe medical applications of Ligularia.

Comparative chloroplast genomics between the invasive weed Mikania micrantha and its indigenous congener Mikania cordata: Structure variation, identification of highly divergent regions, divergence time estimation, and phylogenetic analysis

Mikania micrantha and Mikania cordata are the only two species in genus Mikania (Asteraceae) in China. They share very similar morphological and life-history characteristics but occupy quite different habitats. Most importantly, they generate totally different ecological consequences. While M. micrantha has become an exotic invasive weed, M. cordata exists as an indigenous species with no harmful effects on native plants or habitats. As a continuous study of our previously reported M. micrantha chloroplast (cp) genome, in this study we have further sequenced the M. cordata cp genome to (1) conduct a comparative genome analysis to gain insights into the mechanism of invasiveness; (2) develop cp markers to examine the population genetic adaptation of M. micrantha; and (3) screen variable genome regions of phylogenetic utility. The M. cordata chloroplast genome is 151,984 bp in length and displays a typical quadripartite structure. The number and distribution of protein coding genes, tRNA genes, and rRNA genes of M. cordata are identical to those of M. micrantha. The main difference lays in that the pseudogenization of ndhF and a 118-bp palindromic repeat only arises in M. cordata. Fourteen highly divergent regions, 235 base substitutions, and 58 indels were identified between the two cp genomes. Phylogenetic inferences revealed a sister relationship between M. micrantha and M. cordata whose divergence was estimated to occur around 1.78 million years ago (MYA). Twelve cpSSR loci were detected to be polymorphic and adopted to survey the genetic adaptation of M. micrantha populations. No cpSSR loci were found to undergo selection. Our results build a foundation to examine the invasive mechanism of Mikania weed.

Different Natural Selection Pressures on the atpF Gene in Evergreen Sclerophyllous and Deciduous Oak Species: Evidence from Comparative Analysis of the Complete Chloroplast Genome of Quercus aquifolioides with Other Oak Species

Quercus is an economically important and phylogenetically complex genus in the family Fagaceae. Due to extensive hybridization and introgression, it is considered to be one of the most challenging plant taxa, both taxonomically and phylogenetically. Quercus aquifolioides is an evergreen sclerophyllous oak species that is endemic to, but widely distributed across, the Hengduanshan Biodiversity Hotspot in the Eastern Himalayas. Here, we compared the fully assembled chloroplast (cp) genome of Q. aquifolioides with those of three closely related species. The analysis revealed a cp genome ranging in size from 160,415 to 161,304 bp and with a typical quadripartite structure, composed of two inverted repeats (IRs) separated by a small single copy (SSC) and a large single copy (LSC) region. The genome organization, gene number, gene order, and GC content of these four Quercus cp genomes are similar to those of many angiosperm cp genomes. We also analyzed the Q. aquifolioides repeats and microsatellites. Investigating the effects of selection events on shared protein-coding genes using the Ka/Ks ratio showed that significant positive selection had acted on the atpF gene of Q. aquifolioides compared to two deciduous oak species, and that there had been significant purifying selection on the atpF gene in the chloroplast of evergreen sclerophyllous oak trees. In addition, site-specific selection analysis identified positively selected sites in 12 genes. Phylogenetic analysis based on shared protein-coding genes from 14 species defined Q. aquifolioides as belonging to sect. Heterobalanus and being closely related to Q. rubra and Q. aliena. Our findings provide valuable genetic information for use in accurately identifying species, resolving taxonomy, and reconstructing the phylogeny of the genus Quercus.

Application of Chloroplast Phylogenomics to Resolve Species Relationships Within the Plant Genus Amaranthus

Amaranthus species are an emerging and promising nutritious traditional vegetable food source. Morphological plasticity and poorly resolved dendrograms have led to the need for well resolved species phylogenies. We hypothesized that whole chloroplast phylogenomics would result in more reliable differentiation between closely related amaranth species. The aims of the study were therefore: to construct a fully assembled, annotated chloroplast genome sequence of Amaranthus tricolor; to characterize Amaranthus accessions phylogenetically by comparing barcoding genes (matK, rbcL, ITS) with whole chloroplast sequencing; and to use whole chloroplast phylogenomics to resolve deeper phylogenetic relationships. We generated a complete A. tricolor chloroplast sequence of 150,027 bp. The three barcoding genes revealed poor inter- and intra-species resolution with low bootstrap support. Whole chloroplast phylogenomics of 59 Amaranthus accessions increased the number of parsimoniously informative sites from 92 to 481 compared to the barcoding genes, allowing improved separation of amaranth species. Our results support previous findings that two geographically independent domestication events of Amaranthus hybridus likely gave rise to several species within the Hybridus complex, namely Amaranthus dubius, Amaranthus quitensis, Amaranthus caudatus, Amaranthus cruentus and Amaranthus hypochondriacus. Poor resolution of species within the Hybridus complex supports the recent and ongoing domestication within the complex, and highlights the limitation of chloroplast data for resolving recent evolution. The weedy Amaranthus retroflexus and Amaranthus powellii was found to share a common ancestor with the Hybridus complex. Leafy amaranth, Amaranthus tricolor, Amaranthus blitum, Amaranthus viridis and Amaranthus graecizans formed a stable sister lineage to the aforementioned species across the phylogenetic trees. This study demonstrates the power of next-generation sequencing data and reference-based assemblies to resolve phylogenies, and also facilitated the identification of unknown Amaranthus accessions from a local genebank. The informative phylogeny of the Amaranthus genus will aid in selecting accessions for breeding advanced genotypes to satisfy global food demand.