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

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

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

Comparative genomics revealed new insights into the plastome evolution of Ludwigia (Onagraceae, Myrtales)

The primrose-willow (Ludwigia L.), a well-defined genus of the Onagraceae family, comprises 87 species widely distributed worldwide. In this study, we sequenced and characterized the complete chloroplast (cp) genomes of three species in the genus, including Ludwigia adscendens, Ludwigia hyssopifolia, and Ludwigia prostrata. Three Ludwigia cp genomes ranged from 158,354 to 159,592 bp in size, and each contained 113 genes, including 79 unique protein-coding genes (PCGs), four rRNA genes, and 30 tRNA genes. A comparison of the Ludwigia cp genomes revealed that they were highly conserved in gene composition, gene orientation, and GC content. Moreover, we compared the structure of cp genomes and reconstructed phylogenetic relationships with related species in the Onagraceae family. Regarding contraction/expansion of inverted repeat (IR) region, two kinds of expansion IR region structures were found in Oenothera, Chamaenerion, and Epilobium genera, with primitive IR structures in Ludwigia and Circeae genera. The regions clpP, ycf2, and ycf1 genes possessed highly divergent nucleotides among all available cp genomes of the Onagraceae family. The phylogenetic reconstruction using 79 PCGs from 39 Onagraceae cp genomes inferred that Ludwigia (including L. adscendens, L. hyssopifolia, L. prostrata, and Ludwigia octovalvis) clade was monophyletic and well-supported by the bootstrap and posterior probability values. This study provides the reference cp genomes of three Ludwigia species, which can be used for species identification and phylogenetic reconstruction of Ludwigia and Onagraceae taxa.

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

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

Complete chloroplast genome structure of four Ulmus species and Hemiptelea davidii and comparative analysis within Ulmaceae species

In this study, the chloroplast (cp) genomes of Hemiptelea davidii, Ulmus parvifolia, Ulmus lamellosa, Ulmus castaneifolia, and Ulmus pumila ‘zhonghuajinye’ were spliced, assembled and annotated using the Illumina HiSeq PE150 sequencing platform, and then compared to the cp genomes of other Ulmus and Ulmaceae species. The results indicated that the cp genomes of the five sequenced species showed a typical tetrad structure with full lengths ranging from 159,113 to 160,388 bp. The large single copy (LSC), inverted repeat (IR), and small single copy (SSC) lengths were in the range of 87,736–88,466 bp, 26,317–26,622 bp and 18,485–19,024 bp, respectively. A total of 130–131 genes were annotated, including 85–86 protein-coding genes, 37 tRNA genes and eight rRNA genes. The GC contents of the five species were similar, ranging from 35.30 to 35.62%. Besides, the GC content was different in different region and the GC content in IR region was the highest. A total of 64-133 single sequence repeat (SSR) loci were identified among all 21 Ulmaceae species. The (A)_n and (T)_n types of mononucleotide were highest in number, and the lengths were primarily distributed in 10–12 bp, with a clear AT preference. A branch-site model and a Bayes Empirical Bayes analysis indicated that the rps15 and rbcL had the positive selection sites. Besides, the analysis of mVISTA and sliding windows got a lot of hotspots such as trnH/psbA, rps16/trnQ, trnS/trnG, trnG/trnR and rpl32/trnL, which could be utilized as potential markers for the species identification and phylogeny reconstruction within Ulmus in the further studies. Moreover, the evolutionary tree of Ulmaceae species based on common protein genes, whole cp genome sequences and common genes in IR region of the 23 Ulmaceae species were constructed using the ML method. The results showed that these Ulmaceae species were divided into two branches, one that included Ulmus, Zelkova and Hemiptelea, among which Hemiptelea was the first to differentiate and one that included Celtis, Trema, Pteroceltis, Gironniera and Aphananthe. Besides, these variations found in this study could be used for the classification, identification and phylogenetic study of Ulmus species. Our study provided important genetic information to support further investigations into the phylogenetic development and adaptive evolution of Ulmus and Ulmaceae species.

Phylogenetic analysis based on single-copy orthologous proteins in highly variable chloroplast genomes of Corydalis

Corydalis is one of the few lineages that have been reported to have extensive large-scale chloroplast genome (cp-genome) rearrangements. In this study, novel cp-genome rearrangements of Corydalis pinnata, C. mucronate, and C. sheareri are described. C. pinnata is a narrow endemic species only distributed at Qingcheng Mountain in southwest China. Two independent relocations of the same four genes (trnM-CAU-rbcL) were found relocated from the typically posterior part of the large single-copy region to the front of it. A uniform inversion of an 11-14-kb segment (ndhB-trnR-ACG) was found in the inverted repeat region; and extensive losses of accD, clpP, and trnV-UAC genes were detected in all cp-genomes of all three species of Corydalis. In addition, a phylogenetic tree was reconstructed based on 31 single-copy orthologous proteins in 27 cp-genomes. This study provides insights into the evolution of cp-genomes throughout the genus Corydalis and also provides a reference for further studies on the taxonomy, identification, phylogeny, and genetic transformation of other lineages with extensive rearrangements in cp-genomes.

Comparative Analyses of Six Complete Chloroplast Genomes from the Genus Cupressus and Juniperus (Cupressaceae)

Cupressaceae is a conifer family distributed around the world. Cupressus and Juniperus are the main genera of the Cupressaceae family and have important medicinal value. This leads to confusion between Cupressus and Juniperus due to similar morphologies. Here, the complete cp genomes of two Cupressus (C. duclouxiana and C. funebri) and four Juniperus (J. chinensis, J. gaussenii J. pingii and J. procumbens) were sequenced. The results revealed that the length of the cp genomes ranged from 126,996 bp to 129,959 bp, with 119 genes comprising 82 protein-coding genes, 33 transfer RNAs and 4 ribosomal RNAs. All chloroplast genomes of Cupressus and Juniperus lost whole IR regions, which is consistent with gymnosperm cp genome studies. In addition, the number of SSRs per species ranged from 54 to 73 and was dominated by mononucleotide repeats. In the six cp genomes of Cupressus and Juniperus, five highly divergent regions, including accD, accD-rpl2, ycf1, ycf2 and rrn23-rrn4.5, can be used as DNA barcodes of interspecific relationships and potential genetic markers. We compared the gene selection pressures (C. chengiana as reference species), and 6 genes underwent positive selection, the majority of which were related to photosynthesis. Phylogenetic results showed that the monophyly of Cupressus and Juniperus supported most bootstrap support. Cupressus funebris and J. chinensis were resolved to be early diverging species within Cupressus and Juniperus, and the two genera were sister groups to each other. This research revealed a new understanding of the structural pluralism and phylogenetic relationships of Cupressaceae cp genomes. These results will facilitate comprehension of the complexity and diversity of conifer cp genomes. SIGNIFICANCE:: Phylogenetic relationships among Cupressus, Juniperus, and their closest relatives are controversial, and generic delimitations have been in flux for the past decade. To address relationships and attempt to produce a more robust classification, we sequenced 6 new plastid genomes (plastomes) from the two variously described genera in this complex (Cupressus and Juniperus) and compared them with additional plastomes from diverse members of Cupressaceae. Our study corroborated the accD of Cupressophytes have a tendency to expand in size and strongly supported a sister relationship between Cupressus and Juniperus. The disparity in these results could be traced to the facts that the chloroplast genome is uniparentally inherited, also the usage of the whole chloroplast genome for this research is of a better advantage compared to usage of selected genes or portion of the plastome. The complete CP genomic data will provide useful information for studying genetic diversity and species identification, which is important for the overarching goal of biodiversity conservation.

New Methods for Testing/Determining the Environmental Exposure to Glyphosate in Sunflower (Helianthus annuus L.) Plants

Glyphosate is still the subject of much debate, as several studies report its effects on the environment. Sunflower (GK Milia CL) was set up as an experimental plant and treated with glyphosate concentrations of 500 ppm and 1000 ppm in two treatments. Glyphosate was found to be absorbed from the soil into the plant organism through the roots, which was also detectable in the leaf and root. Glyphosate was also significantly detected in the plant 5 weeks after treatment and in plants that did not receive glyphosate treatment directly, so it could be taken up through the soil. Based on the morphological results, treatment with higher concentrations (1000 ppm) of glyphosate increased the dried mass and resulted in shorter, thicker roots. Histological results also showed that basal and transporter tissue distortions were observed in the glyphosate-treated plants compared to the control group. Cells were distorted with increasing concentration, vacuoles formed, and the cell wall was weakened in both the leaf-treated and inter-row-treated groups. In the future, it will be worth exploring alternative agricultural technologies that can reduce the risk of glyphosate while increasing economic outcomes. This may make the use of glyphosate more environmentally conscious.

Extensive reorganization of the chloroplast genome of Corydalis platycarpa: A comparative analysis of their organization and evolution with other Corydalis plastomes

Introduction

The chloroplast (cp) is an autonomous plant organelle with an individual genome that encodes essential cellular functions. The genome architecture and gene content of the cp is highly conserved in angiosperms. The plastome of Corydalis belongs to the Papaveraceae family, and the genome is comprised of unusual rearrangements and gene content. Thus far, no extensive comparative studies have been carried out to understand the evolution of Corydalis chloroplast genomes.

Methods

Therefore, the Corydalis platycarpa cp genome was sequenced, and wide-scale comparative studies were conducted using publicly available twenty Corydalis plastomes.

Results

Comparative analyses showed that an extensive genome rearrangement and IR expansion occurred, and these events evolved independently in the Corydalis species. By contrast, the plastomes of its closely related subfamily Papaveroideae and other Ranunculales taxa are highly conserved. On the other hand, the synapomorphy characteristics of both accD and the ndh gene loss events happened in the common ancestor of the Corydalis and sub-clade of the Corydalis lineage, respectively. The Corydalis-sub clade species (ndh lost) are distributed predominantly in the Qinghai-Tibetan plateau (QTP) region. The phylogenetic analysis and divergence time estimation were also employed for the Corydalis species.

Discussion

The divergence time of the ndh gene in the Corydalis sub-clade species (44.31 - 15.71 mya) coincides very well with the uplift of the Qinghai-Tibet Plateau in Oligocene and Miocene periods, and maybe during this period, it has probably triggered the radiation of the Corydalis species.

Conclusion

To the best of the authors' knowledge, this is the first large-scale comparative study of Corydalis plastomes and their evolution. The present study may provide insights into the plastome architecture and the molecular evolution of Corydalis species.

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

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

Comparative Chloroplast Genomics of Corydalis Species (Papaveraceae): Evolutionary Perspectives on Their Unusual Large Scale Rearrangements

The chloroplast genome (plastome) of angiosperms (particularly photosynthetic members) is generally highly conserved, although structural rearrangements have been reported in a few lineages. In this study, we revealed Corydalis to be another unusual lineage with extensive large-scale plastome rearrangements. In the four newly sequenced Corydalis plastomes that represent all the three subgenera of Corydalis, we detected (1) two independent relocations of the same five genes (trnV-UAC-rbcL) from the typically posterior part of the large single-copy (LSC) region to the front, downstream of either the atpH gene in Corydalis saxicola or the trnK-UUU gene in both Corydalis davidii and Corydalis hsiaowutaishanensis; (2) relocation of the rps16 gene from the LSC region to the inverted repeat (IR) region in Corydalis adunca; (3) uniform inversion of an 11-14 kb segment (ndhB-trnR-ACG) in the IR region of all the four Corydalis species (the same below); (4) expansions (>10 kb) of IR into the small single-copy (SSC) region and corresponding contractions of SSC region; and (5) extensive pseudogenizations or losses of 13 genes (accD, clpP, and 11 ndh genes). In addition, we also found that the four Corydalis plastomes exhibited elevated GC content in both gene and intergenic regions and high number of dispersed repeats. Phylogenomic analyses generated a well-supported topology that was consistent with the result of previous studies based on a few DNA markers but contradicted with the morphological character-based taxonomy to some extent. This study provided insights into the evolution of plastomes throughout the three Corydalis subgenera and will be of value for further study on taxonomy, phylogeny, and evolution of Corydalis.

Mutational Dynamics of Aroid Chloroplast Genomes II

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

Comparative Chloroplast Genomics of Sophora Species: Evolution and Phylogenetic Relationships in the Early-Diverging Legume Subfamily Papilionoideae (Fabaceae)

The taxonomy and evolutionary history of Sophora L., a genus with high economic and medicinal value, remain uncertain due to the absence of genetic resource (especially in China) and low polymorphism of molecular markers. Our aim was to elucidate the molecular evolution and phylogenetic relationships in chloroplast genomes of Sophora species in the early-diverging legume subfamily Papilionoideae (Fabaceae). We reported nine Sophora chloroplast genome from China using Illumina sequencing. We performed a series of analyses with previously published genomes of Sophora species to investigate their genomic characteristics, identified simple sequence repeats, large repeat sequences, tandem repeats, and highly polymorphic loci. The genomes were 152,953-158,087 bp in length, and contained 111-113 unique genes, including 76-78 protein coding, 31 tRNA, and 4 rRNA. The expansion of inverted repeat boundary of Sophora resulted in rps12 entering into the LSC region and loss of trnT-CGU gene in some species. Also, we found an approximately 23 kb inversion between trnC-GCA and trnF-GAA within the genus. In addition, we identified seven highly polymorphic loci (pi (π) > 0.035) suitable for inferring the phylogeny of Sophora species. Among these, three regions also co-occurred with large repeat sequences and support use of repeats as a proxy for the identification of polymorphic loci. Based on whole chloroplast genome and protein-coding sequences data-set, a well-supported phylogenetic tree of Sophora and related taxa showed that this genus is monophyletic, but sect. Disamaea and sect. Sophora, are incongruent with traditional taxonomic classifications based on fruit morphology. Our finding provides significant genetic resources to support further investigation into the phylogenetic relationship and evolution of the genus Sophora.

Plastid Phylogenomics of Dendroseris (Cichorieae; Asteraceae): Insights Into Structural Organization and Molecular Evolution of an Endemic Lineage From the Juan Fernández Islands

Dendroseris D. Don comprises 11 species endemic to the Juan Fernández islands in Chile. They demonstrate spectacular and unusual growth forms of rosette trees with extremely variable morphology and occupy wide ecological ranges on the islands. These unique plants are now highly threatened with extinction with very small population sizes, typically consisting of 10 or fewer individuals in wild. Despite morphological and ecological divergence among species of Dendroseris, their monophyly has been supported in previous studies, but with little resolution among subgeneric groups. We assembled seven complete plastome sequences from seven species of Dendroseris, including representatives from three subgenera, and carried out comparative phylogenomic analyses. The plastomes are highly conserved in gene content and order, with size ranging from 152,199 to 152,619 bp and containing 130 genes (87 coding genes, 6 rRNA genes, and 37 tRNA genes). Plastid phylogenomic analyses based on both the complete plastome sequences and 81 concatenated coding genes only show Dendroseris nested within Sonchus sensu lato, and also that inter-subgeneric relationships are fully resolved. Subg. Phoenicoseris is resolved as sister to the remaining species of the genus and a sister relationship between the two subgenera Dendroseris and Rea. Ten mutation hotspots from LSC and SSC regions and variable SSRs are identified as potential chloroplast markers for future phylogenetic and phylogeographic studies of Sonchus and related groups.

Ancient Plant Genomics in Archaeology, Herbaria, and the Environment

The ancient DNA revolution of the past 35 years has driven an explosion in the breadth, nuance, and diversity of questions that are approachable using ancient biomolecules, and plant research has been a constant, indispensable facet of these developments. Using archaeological, paleontological, and herbarium plant tissues, researchers have probed plant domestication and dispersal, plant evolution and ecology, paleoenvironmental composition and dynamics, and other topics across related disciplines. Here, we review the development of the ancient DNA discipline and the role of plant research in its progress and refinement. We summarize our understanding of long-term plant DNA preservation and the characteristics of degraded DNA. In addition, we discuss challenges in ancient DNA recovery and analysis and the laboratory and bioinformatic strategies used to mitigate them. Finally, we review recent applications of ancient plant genomic research.

Long-reads reveal that Rhododendron delavayi plastid genome contains extensive repeat sequences, and recombination exists among plastid genomes of photosynthetic Ericaceae

BACKGROUND

Rhododendron delavayi Franch. var. delavayi is a wild ornamental plant species in Guizhou Province, China. The lack of its plastid genome information seriously hinders the further application and conservation of the valuable resource.

METHODS

The complete plastid genome of R. delavayi was assembled from long sequence reads. The genome was then characterized, and compared with those of other photosynthetic Ericaceae species.

RESULTS

The plastid genome of R. delavayi has a typical quadripartite structure, and a length of 202,169 bp. It contains a large number of repeat sequences and shows preference for codon usage. The comparative analysis revealed the irregular recombination of gene sets, including rearrangement and inversion, in the large single copy region. The extreme expansion of the inverted repeat region shortened the small single copy, and expanded the full length of the genome. In addition, consistent with traditional taxonomy, R. delavayi with nine other species of the same family were clustered into Ericaceae based on the homologous protein-coding sequences of the plastid genomes. Thus, the long-read assembly of the plastid genome of R. delavayi would provide basic information for the further study of the evolution, genetic diversity, and conservation of R. delavayi and its relatives.

Development of chloroplast microsatellite markers for giant ragweed (Ambrosia trifida)

PREMISE

Plant invasions are increasing globally, and extensive study of the genetic background of the source and invading populations is needed to understand such biological processes. For this reason, chloroplast microsatellite markers were identified to explore the genetic diversity of the noxious weed Ambrosia trifida (Asteraceae).

METHODS AND RESULTS

The complete chloroplast genome of A. trifida was mined for microsatellite loci, and 15 novel chloroplast primers were identified to assess the genetic diversity of 49 Ambrosia samples. The number of alleles amplified ranged from two to six, with an average of 3.2 alleles per locus. Shannon's information index varied from 0.305 and 1.467, expected heterozygosity ranged from 0.178 to 0.645, and the polymorphism information content value ranged from 0.211 to 0.675 (average 0.428). The cross-species transferability of the 15 microsatellite loci was also evaluated in four related Ambrosia species (A. artemisiifolia, A. maritima, A. psilostachya, and A. tenuifolia).

CONCLUSIONS

The novel chloroplast microsatellite markers developed in the current study demonstrate substantial cross-species transferability and will be helpful in future genetic diversity studies of A. trifida and related species.

Sow Thistle Chloroplast Genomes: Insights into the Plastome Evolution and Relationship of Two Weedy Species, Sonchus asper and Sonchus oleraceus (Asteraceae)

Prickly sow thistle, Sonchus asper (L.) Hill, and common sow thistle, Sonchus oleraceus L., are noxious weeds. Probably originating from the Mediterranean region, they have become widespread species. They share similar morphology and are closely related. However, they differ in their chromosome numbers and the precise relationship between them remains uncertain. Understanding their chloroplast genome structure and evolution is an important initial step toward determining their phylogenetic relationships and analyzing accelerating plant invasion processes on a global scale. We assembled four accessions of chloroplast genomes (two S. asper and two S. oleraceus) by the next generation sequencing approach and conducted comparative genomic analyses. All the chloroplast genomes were highly conserved. Their sizes ranged from 151,808 to 151,849 bp, containing 130 genes including 87 coding genes, 6 rRNA genes, and 37 tRNA genes. Phylogenetic analysis based on the whole chloroplast genome sequences showed that S. asper shares a recent common ancestor with S. oleraceus and suggested its likely involvement in a possible amphidiploid origin of S. oleraceus. In total, 79 simple sequence repeats and highly variable regions were identified as the potential chloroplast markers to determine genetic variation and colonization patterns of Sonchus species.