Exploring the History of Chloroplast Capture in Arabis Using Whole Chloroplast Genome Sequencing

Phylogenetic incongruence in an Asiatic species complex of the genus Caryodaphnopsis (Lauraceae)

BACKGROUND

Caryodaphnopsis, a group of tropical trees (ca. 20 spp.) in the family Lauraceae, has an amphi-Pacific disjunct distribution: ten species are distributed in Southeast Asia, while eight species are restricted to tropical rainforests in South America. Previously, phylogenetic analyses using two nuclear markers resolved the relationships among the five species from Latin America. However, the phylogenetic relationships between the species in Asia remain poorly known.

RESULTS

Here, we first determined the complete mitochondrial genome (mitogenome), plastome, and the nuclear ribosomal cistron (nrDNA) sequences of C. henryi with lengths of 1,168,029 bp, 154,938 bp, and 6495 bp, respectively. We found 2233 repeats and 368 potential SSRs in the mitogenome of C. henryi and 50 homologous DNA fragments between its mitogenome and plastome. Gene synteny analysis revealed a mass of rearrangements in the mitogenomes of Magnolia biondii, Hernandia nymphaeifolia, and C. henryi and only six conserved clustered genes among them. In order to reconstruct relationships for the ten Caryodaphnopsis species in Asia, we created three datasets: one for the mitogenome (coding genes and ten intergenic regions), another for the plastome (whole genome), and the other for the nuclear ribosomal cistron. All of the 22 Caryodaphnopsis individuals were divided into four, five, and six different clades in the phylogenies based on mitogenome, plastome, and nrDNA datasets, respectively.

CONCLUSIONS

The study showed phylogenetic conflicts within and between nuclear and organellar genome data of Caryodaphnopsis species. The sympatric Caryodaphnopsis species in Hekou and Malipo SW China may be related to the incomplete lineage sorting, chloroplast capture, and/or hybridization, which mixed the species as a complex in their evolutionary history.

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.

Structural divergence and phylogenetic relationships of Ajania (Asteraceae) from plastomes and ETS

BACKGROUND

Ajania Poljakov, an Asteraceae family member, grows mostly in Asia's arid and semi-desert areas and is a significant commercial and decorative plant. Nevertheless, the genus' classification has been disputed, and the evolutionary connections within the genus have not been thoroughly defined. Hence, we sequenced and analyzed Ajania's plastid genomes and combined them with ETS data to assess their phylogenetic relationships.

RESULTS

We obtained a total of six new Ajania plastid genomes and nine ETS sequences. The whole plastome lengths of the six species sampled ranged from 151,002 bp to 151,115 bp, showing conserved structures. Combined with publicly available data from GenBank, we constructed six datasets to reconstruct the phylogenetic relationships, detecting nucleoplasmic clashes. Our results reveal the affinities of Artemisia, Chrysanthemum and Stilpnolepis to Ajania and validate the early taxonomy reclassification. Some of the plastid genes with low phylogenetic information and gene trees with topological differences may have contributed to the ambiguous phylogenetic results of Ajania. There is extensive evolutionary rate heterogeneity in plastid genes. The psbH and ycf2 genes, which are involved in photosynthesis and ATP transport, are under selective pressure. Plastomes from Ajania species diverged, and structural aspects of plastomes may indicate some of the real evolutionary connections. We suggest the ycf1 gene as a viable plastid DNA barcode because it has significant nucleotide diversity and better reflects evolutionary connections.

CONCLUSION

Our findings validate the early Ajania taxonomy reclassification and show evolutionary rate heterogeneity, genetic variety, and phylogenetic heterogeneity of plastid genes. This research might provide new insights into the taxonomy and evolution of Ajania, as well as provide useful information for germplasm innovation and genetic enhancement in horticultural species.

Genomic Analysis of Plastid-Nuclear Interactions and Differential Evolution Rates in Coevolved Genes across Juglandaceae Species

The interaction between the nuclear and chloroplast genomes in plants is crucial for preserving essential cellular functions in the face of varying rates of mutation, levels of selection, and modes of transmission. Despite this, identifying nuclear genes that coevolve with chloroplast genomes at a genome-wide level has remained a challenge. In this study, we conducted an evolutionary rate covariation analysis to identify candidate nuclear genes coevolving with chloroplast genomes in Juglandaceae. Our analysis was based on 4,894 orthologous nuclear genes and 76 genes across seven chloroplast partitions in nine Juglandaceae species. Our results indicated that 1,369 (27.97%) of the nuclear genes demonstrated signatures of coevolution, with the Ycf1/2 partition yielding the largest number of hits (765) and the ClpP1 partition yielding the fewest (13). These hits were found to be significantly enriched in biological processes related to leaf development, photoperiodism, and response to abiotic stress. Among the seven partitions, AccD, ClpP1, MatK, and RNA polymerase partitions and their respective hits exhibited a narrow range, characterized by dN/dS values below 1. In contrast, the Ribosomal, Photosynthesis, Ycf1/2 partitions and their corresponding hits, displayed a broader range of dN/dS values, with certain values exceeding 1. Our findings highlight the differences in the number of candidate nuclear genes coevolving with the seven chloroplast partitions in Juglandaceae species and the correlation between the evolution rates of these genes and their corresponding chloroplast partitions.

The evolutionary history of Cardamine bulbifera shows a successful rapid postglacial Eurasian range expansion in the absence of sexual reproduction

BACKGROUND AND AIMS

Sexual reproduction is known to drive plant diversification and adaptation. Here we investigate the evolutionary history and spatiotemporal origin of a dodecaploid (2n = 12x = 96) Eurasian deciduous woodland species, Cardamine bulbifera, which reproduces and spreads via vegetative bulb-like structures only. The species has been among the most successful range-expanding understorey woodland plants in Europe, which raises the question of the genetic architecture of its gene pool, since its hexaploid (2n = 6x = 48) but putatively outcrossing closest relative, C. quinquefolia, displays a smaller distribution range in Eastern Europe towards the Caucasus region. Cardamine bulbifera belongs to a small monophyletic clade of four species comprising also C. abchasica (2n = 2x = 16) and C. bipinnata (unknown ploidy) from the Caucasus region.

METHODS

We sequenced the genomes of the two polyploids and their two putative ancestors using Illumina short-read sequencing technology (×7-8 coverage). Covering the entire distribution range, genomic data were generated for 67 samples of the two polyploids (51 samples of C. bulbifera, 16 samples of C. quinquefolia) and 6 samples of the putative diploid taxa (4 samples of C. abchasica, 2 samples of C. bipinnata) to unravel the evolutionary origin of the polyploid taxa using phylogenetic reconstructions of biparentally and maternally inherited genetic sequence data. Ploidy levels of C. bulbifera and C. quinquefolia were analysed by comparative chromosome painting. We used genetic assignment analysis (STRUCTURE) and approximate Bayesian computation (ABC) modelling to test whether C. bulbifera represents genetically differentiated lineages and addressed the hypothesis of its hybrid origin. Comparative ecological modelling was applied to unravel possible niche differentiation among the two polyploid species.

KEY RESULTS

Cardamine bulbifera was shown to be a non-hybridogenous, auto-dodecaploid taxon of early Pleistocene origin, but with a history of past gene flow with its hexaploid sister species C. quinquefolia, likely during the last glacial maximum in shared refuge areas in Eastern Europe towards Western Turkey and the Crimean Peninsula region. The diploid Caucasian endemic C. abchasica is considered an ancestral species, which also provides evidence for the origin of the species complex in the Caucasus region. Cardamine bulbifera successfully expanded its distribution range postglacially towards Central and Western Europe accompanied by a transition to exclusively vegetative propagation.

CONCLUSIONS

A transition to vegetative propagation in C. bulbifera is hypothesized as the major innovation to rapidly expand its distribution range following postglacially progressing woodland vegetation throughout Europe. Preceding and introgressive gene flow from its sister species C. quinquefolia in the joint refuge area is documented. This transition and ecological differentiation may have been triggered by preceding introgressive gene flow from its sister species in the joint East European refuge areas.

Sequence Characteristics and Phylogenetic Analysis of the Artemisia argyi Chloroplast Genome

Artemisia argyi Levl. et Van is an important Asteraceae species with a high medicinal value. There are abundant A. argyi germplasm resources in Asia, especially in China, but the evolutionary relationships of these varieties and the systematic localization of A. argyi in the family Asteraceae are still unclear. In this study, the chloroplast (cp) genomes of 72 A. argyi varieties were systematically analyzed. The 72 varieties originated from 47 regions in China at different longitudes, latitudes and altitudes, and included both wild and cultivated varieties. The A. argyi cp genome was found to be ∼151 kb in size and to contain 114 genes, including 82 protein-coding, 28 tRNA, and 4 rRNA genes. The number of short sequence repeats (SSRs) in A. argyi cp genomes ranged from 35 to 42, and most of them were mononucleotide A/T repeats. A total of 196 polymorphic sites were detected in the cp genomes of the 72 varieties. Phylogenetic analysis demonstrated that the genetic relationship between A. argyi varieties had a weak relationship with their geographical distribution. Furthermore, inverted repeat (IR) boundaries of 10 Artemisia species were found to be significantly different. A sequence divergence analysis of Asteraceae cp genomes showed that the variable regions were mostly located in single-copy (SC) regions and that the coding regions were more conserved than the non-coding regions. A phylogenetic tree was constructed using 43 protein-coding genes common to 67 Asteraceae species. The resulting tree was consistent with the traditional classification system; Artemisia species were clustered into one group, and A. argyi was shown to be closely related to Artemisia lactiflora and Artemisia montana. In summary, this study systematically analyzed the cp genome characteristics of A. argyi and compared cp genomes of Asteraceae species. The results provide valuable information for the definitive identification of A. argyi varieties and for the understanding of the evolutionary relationships between Asteraceae species.

Complete chloroplast genome features of the model heavy metal hyperaccumulator Arabis paniculata Franch and its phylogenetic relationships with other Brassicaceae species

Arabis paniculata Franch (Brassicaceae) has been widely used for the phytoremediation of heavy mental, owing to its hyper tolerance of extreme Pb, Zn, and Cd concentrations. However, studies on its genome or plastid genome are scarce. In the present study, we obtained the complete chloroplast (cp) genome of A. paniculata via de novo assembly through the integration of Illumina reads and PacBio subreads. The cp genome presents a typical quadripartite cycle with a length of 153,541 bp, and contains 111 unigenes, with 79 protein-coding genes, 28 tRNAs and 4 rRNAs. Codon usage analysis showed that the codons for leucine were the most frequent codons and preferentially ended with A/U. Synonymous (Ks) and non-synonymous (Ka) substitution rate analysis indicated that the unigenes, ndhF and rpoC2, related to "NADH-dehydrogenase" and "RNA polymerase" respectively, underwent the lowest purifying selection pressure. Phylogenetic analysis demonstrated that Arabis flagellosa and A. hirsuta are more similar to each other than to A. paniculata, and Arabis is the closest relative of Draba among all Brassicaceae genera. These findings provide valuable information for the optimal exploitation of this model species as a heavy-metal hyperaccumulator.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01151-1.

Approximate Bayesian computation untangles signatures of contemporary and historical hybridization between two endangered species

Contemporary gene flow, when resumed after a period of isolation, can have crucial consequences for endangered species, as it can both increase the supply of adaptive alleles and erode local adaptation. Determining the history of gene flow and thus the importance of contemporary hybridization, however, is notoriously difficult. Here, we focus on two endangered plant species, Arabis nemorensis and A. sagittata, which hybridize naturally in a sympatric population located on the banks of the Rhine. Using reduced genome sequencing, we determined the phylogeography of the two taxa but report only a unique sympatric population. Molecular variation in chloroplast DNA indicated that A. sagittata is the principal receiver of gene flow. Applying classical D-statistics and its derivatives to whole-genome data of 35 accessions, we detect gene flow not only in the sympatric population but also among allopatric populations. Using an Approximate Bayesian computation approach, we identify the model that best describes the history of gene flow between these taxa. This model shows that low levels of gene flow have persisted long after speciation. Around 10 000 years ago, gene flow stopped and a period of complete isolation began. Eventually, a hotspot of contemporary hybridization was formed in the unique sympatric population. Occasional sympatry may have helped protect these lineages from extinction in spite of their extremely low diversity.

Variation and Evolution of the Whole Chloroplast Genomes of Fragaria spp. (Rosaceae)

Species identification is vital for protecting species diversity and selecting high-quality germplasm resources. Wild Fragaria spp. comprise rich and excellent germplasm resources; however, the variation and evolution of the whole chloroplast (cp) genomes in the genus Fragaria have been ignored. In the present study, 27 complete chloroplast genomes of 11 wild Fragaria species were sequenced using the Illumina platform. Then, the variation among complete cp genomes of Fragaria was analyzed, and phylogenetic relationships were reconstructed from those genome sequences. There was an overall high similarity of sequences, with some divergence. According to analysis with mVISTA, non-coding regions were more variable than coding regions. Inverted repeats (IRs) were observed to contract or expand to different degrees, which resulted in different sizes of cp genomes. Additionally, five variable loci, trnS-trnG, trnR-atpA, trnC-petN, rbcL-accD, and psbE-petL, were identified that could be used to develop DNA barcoding for identification of Fragaria species. Phylogenetic analyses based on the whole cp genomes supported clustering all species into two groups (A and B). Group A species were mainly distributed in western China, while group B contained several species from Europe and Americas. These results support allopolyploid origins of the octoploid species F. chiloensis and F. virginiana and the tetraploid species F. moupinensis and F. tibetica. The complete cp genomes of these Fragaria spp. provide valuable information for selecting high-quality Fragaria germplasm resources in the future.

Plastid Phylogenomic Data Offers Novel Insights Into the Taxonomic Status of the Trichosanthes kirilowii Complex (Cucurbitaceae) in South Korea

Trichosanthes is a genus in Cucurbitaceae comprising 90-100 species. Trichosanthes species are valuable as herbaceous medicinal ingredients. The fruits, seeds, and roots of species such as T. kirilowii and T. rosthornii are used in Korean traditional herbal medicines. T. rosthornii is only found in China, whereas in South Korea two varieties, T. kirilowii var. kirilowii and T. kirilowii var. japonica, are distributed. T. kirilowii var. kirilowii and T. kirilowii var. japonica have different fruit and leaf shapes but are recognized as belonging to the same species. Furthermore, although its members have herbal medicine applications, genomic information of the genus is still limited. The broad goals of this study were (i) to evaluate the taxonomy of Trichosanthes using plastid phylogenomic data and (ii) provide molecular markers specific for T. kirilowii var. kirilowii and T. kirilowii var. japonica, as these have differences in their pharmacological effectiveness and thus should not be confused and adulterated. Comparison of five Trichosanthes plastid genomes revealed locally divergent regions, mainly within intergenic spacer regions (trnT-UGU-trnL-UAA: marker name Tri, rrn4.5-rrn5: TRr, trnE-UUC-trnT-GGU: TRtt). Using these three markers as DNA-barcodes for important herbal medicine species in Trichosanthes, the identity of Trichosanthes material in commercial medicinal products in South Korea could be successfully determined. Phylogenetic analysis of the five Trichosanthes species revealed that the species are clustered within tribe Sicyoeae. T. kirilowii var. kirilowii and T. rosthornii formed a clade with T. kirilowii var. japonica as their sister group. As T. kirilowii in its current circumscription is paraphyletic and as the two varieties can be readily distinguished morphologically (e.g., in leaf shape), T. kirilowii var. japonica should be treated (again) as an independent species, T. japonica.

Molecular Phylogenetics and Micromorphology of Australasian Stipeae (Poaceae, Subfamily Pooideae), and the Interrelation of Whole-Genome Duplication and Evolutionary Radiations in This Grass Tribe

The mainly Australian grass genus Austrostipa (tribe Stipeae) comprising approximately 64 species represents a remarkable example of an evolutionary radiation. To investigate aspects of diversification, macro- and micromorphological variation in this genus, we conducted molecular phylogenetic and scanning electron microscopy (SEM) analyses including representatives from most of Austrostipa's currently accepted subgenera. Because of its taxonomic significance in Stipeae, we studied the lemma epidermal pattern (LEP) in 34 representatives of Austrostipa. Plastid DNA variation within Austrostipa was low and only few lineages were resolved. Nuclear ITS and Acc1 yielded comparable groupings of taxa and resolved subgenera Arbuscula, Petaurista, and Bambusina in a common clade and as monophyletic. In most of the Austrostipa species studied, the LEP was relatively uniform (typical maize-like), but six species had a modified cellular structure. The species representing subgenera Lobatae, Petaurista, Bambusina as well as A. muelleri from subg. Tuberculatae were well-separated from all the other species included in the analysis. We suggest recognizing nine subgenera in Austrostipa (with number of species): Arbuscula (4), Aulax (2), Austrostipa (36), Bambusina (2), Falcatae (10), Lobatae (5), Longiaristatae (2), Petaurista (2) and the new subgenus Paucispiculatae (1) encompassing A. muelleri. Two paralogous sequence copies of Acc1, forming two distinct clades, were found in polyploid Austrostipa and Anemanthele. We found analogous patterns for our samples of Stipa s.str. with their Acc1 clades strongly separated from those of Austrostipa and Anemanthele. This underlines a previous hypothesis of Tzvelev (1977) that most extant Stipeae are of hybrid origin. We also prepared an up-to-date survey and reviewed the chromosome number variation for our molecularly studied taxa and the whole tribe Stipeae. The chromosome base number patterns as well as dysploidy and whole-genome duplication events were interpreted in a phylogenetic framework. The rather coherent picture of chromosome number variation underlines the enormous phylogenetic and evolutionary significance of this frequently ignored character.

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.

Comparing and phylogenetic analysis chloroplast genome of three Achyranthes species

In this study, the chloroplast genome sequencing of the Achyranthes longifolia, Achyranthes bidentata and Achyranthes aspera were performed by Next-generation sequencing technology. The results revealed that there were a length of 151,520 bp (A. longifolia), 151,284 bp (A. bidentata), 151,486 bp (A. aspera), respectively. These chloroplast genome have a highly conserved structure with a pair of inverted repeat (IR) regions (25,150 bp; 25,145 bp; 25,150 bp), a large single copy (LSC) regions (83,732 bp; 83,933 bp; 83,966 bp) and a small single copy (SSC) regions (17,252 bp; 17,263 bp; 17,254 bp) in A. bidentate, A. aspera and A. longifolia. There were 127 genes were annotated, which including 8 rRNA genes, 37 tRNA genes and 82 functional genes. The phylogenetic analysis strongly revealed that Achyranthes is monophyletic, and A. bidentata was the closest relationship with A. aspera and A. longifolia. A. bidentata and A. longifolia were clustered together, the three Achyranthes species had the same origin, then the gunes of Achyranthes is the closest relative to Alternanthera, and that forms a group with Alternanthera philoxeroides. The research laid a foundation and provided relevant basis for the identification of germplasm resources in the future.

The complete chloroplast genome sequence of yellow mustard (Sinapis alba L.) and its phylogenetic relationship to other Brassicaceae species

As a member of the large Brassicaceae family, yellow mustard (Sinapis alba L.) has been used as an important gene pool for the genetic improvement of cash crops in Brassicaceae. Understanding the phylogenetic relationship between Sinapis alba (S. alba) and other Brassicaceae crops can provide guidance on the introgression of its favorable alleles into related species. The chloroplast (cp) genome is an ideal model for assessing genome evolution and the phylogenetic relationships of complex angiosperm families. Herein, we de novo assembled the complete cp genome of S. alba by integrating the PacBio and Illumina sequencing platforms. A 153,760 bp quadripartite cycle without any gap was obtained, including a pair of inverted repeats (IRa and IRb) of 26,221 bp, separated by a large single copy (LSC) region of 83,506 bp and a small single copy (SSC) region of 17,821 bp. A total of 78 protein-coding genes, 30 tRNA genes, and four rRNA genes were identified in this cp genome, as were 89 simple sequence repeat (SSR) loci of 18 types. The codon usage analysis revealed a preferential use of the Leu codon with the A/U ending. The phylogenetic analysis using 82 Brassicaceae species demonstrated that S. alba had a close relationship with important Brassica and Raphanus species; moreover, it likely originated from a separate evolutionary pathway compared with the congeneric Sinapis arvensis. The synonymous (Ks) and non-synonymous (Ks) substitution rate analysis showed that genes encoding "Subunits of cytochrome b/f complex" were under the lowest purifying selection pressure, whereas those associated with "Maturase", "Subunit of acetyl-CoA", and "Subunits of NADH-dehydrogenase" underwent relatively higher purifying selection pressures. Our results provide valuable information for fully utilizing the S. alba cp genome as a potential genetic resource for the genetic improvement of Brassica and Raphanus species.

Complete chloroplast genomes of two Siraitia Merrill species: Comparative analysis, positive selection and novel molecular marker development

Siraitia grosvenorii fruit, known as Luo-Han-Guo, has been used as a traditional Chinese medicine for many years, and mogrosides are its primary active ingredients. Unfortunately, Siraitia siamensis, its wild relative, might be misused due to its indistinguishable appearance, not only threatening the reliability of the medication but also partly exacerbating wild resource scarcity. Therefore, high-resolution genetic markers must be developed to discriminate between these species. Here, the complete chloroplast genomes of S. grosvenorii and S. siamensis were assembled and analyzed for the first time; they were 158,757 and 159,190 bp in length, respectively, and possessed conserved quadripartite circular structures. Both contained 134 annotated genes, including 8 rRNA, 37 tRNA and 89 protein-coding genes. Twenty divergences (Pi > 0.03) were found in the intergenic regions. Nine protein-coding genes, accD, atpA, atpE, atpF, clpP, ndhF, psbH, rbcL, and rpoC2, underwent selection within Cucurbitaceae. Phylogenetic relationship analysis indicated that these two species originated from the same ancestor. Finally, four pairs of molecular markers were developed to distinguish the two species. The results of this study will be beneficial for taxonomic research, identification and conservation of Siraitia Merrill wild resources in the future.

Next-Generation Genome Sequencing of Sedum plumbizincicola Sheds Light on the Structural Evolution of Plastid rRNA Operon and Phylogenetic Implications within Saxifragales

The genus Sedum, with about 470 recognized species, is classified in the family Crassulaceae of the order Saxifragales. Phylogenetic relationships within the Saxifragales are still unresolved and controversial. In this study, the plastome of S. plumbizincicola was firstly presented, with a focus on the structural analysis of rrn operon and phylogenetic implications within the order Saxifragaceae. The assembled complete plastome of S. plumbizincicola is 149,397 bp in size, with a typical circular, double-stranded, and quadripartite structure of angiosperms. It contains 133 genes, including 85 protein-coding genes (PCGs), 36 tRNA genes, 8 rRNA genes, and four pseudogenes (one ycf1, one rps19, and two ycf15). The predicted secondary structure of S. plumbizincicola 16S rRNA includes three main domains organized in 74 helices. Further, our results confirm that 4.5S rRNA of higher plants is associated with fragmentation of 23S rRNA progenitor. Notably, we also found the sequence of putative rrn5 promoter has some evolutionary implications within the order Saxifragales. Moreover, our phylogenetic analyses suggested that S. plumbizincicola had a closer relationship with S. sarmentosum than S. oryzifolium, and supported the taxonomic revision of Phedimus. Our findings of the present study will be useful for further investigation of the evolution of plastid rRNA operon and phylogenetic relationships within Saxifragales.

Sequencing and Analysis of Chrysanthemum carinatum Schousb and Kalimeris indica. The Complete Chloroplast Genomes Reveal Two Inversions and rbcL as Barcoding of the Vegetable

Chrysanthemum carinatum Schousb and Kalimeris indica are widely distributed edible vegetables and the sources of the Chinese medicine Asteraceae. The complete chloroplast (cp) genome of Asteraceae usually occurs in the inversions of two regions. Hence, the cp genome sequences and structures of Asteraceae species are crucial for the cp genome genetic diversity and evolutionary studies. Hence, in this paper, we have sequenced and analyzed for the first time the cp genome size of C. carinatum Schousb and K. indica, which are 149,752 bp and 152,885 bp, with a pair of inverted repeats (IRs) (24,523 bp and 25,003) separated by a large single copy (LSC) region (82,290 bp and 84,610) and a small single copy (SSC) region (18,416 bp and 18,269), respectively. In total, 79 protein-coding genes, 30 distinct transfer RNA (tRNA) genes, four distinct rRNA genes and two pseudogenes were found not only in C. carinatum Schousb but also in the K. indica cp genome. Fifty-two (52) and fifty-nine (59) repeats, and seventy (70) and ninety (90) simple sequence repeats (SSRs) were found in the C. carinatum Schousb and K. indica cp genomes, respectively. Codon usage analysis showed that leucine, isoleucine, and serine are the most frequent amino acids and that the UAA stop codon was the significantly favorite stop codon in both cp genomes. The two inversions, the LSC region ranging from trnC-GCA to trnG-UCC and the whole SSC region were found in both of them. The complete cp genome comparison with other Asteraceae species showed that the coding area is more conservative than the non-coding area. The phylogenetic analysis revealed that the rbcL gene is a good barcoding marker for identifying different vegetables. These results give an insight into the identification, the barcoding, and the understanding of the evolutionary model of the Asteraceae cp genome.