Complete chloroplast genomes from apomictic Taraxacum (Asteraceae): Identity and variation between three microspecies

Comparative Plastome Analysis Between Endangered Mangrove Species Acanthus ebracteatus and Acanthus Relatives Provides Insights into Its Origin and Adaptive Evolution

Acanthus ebracteatus is a typical true mangrove species with great ecological and medicinal values. However, it has become endangered in China. Moreover, because of the similar morphology and distribution, it is commonly confused with the congeneric mangrove species, A. ilicifolius, which poses challenges to the protection and proper medicinal utilization of A. ebracteatus. Plastomes provide a solution for molecular identification and adaptive evolution investigation of plants. In this study, we dissected the complete plastome for A. ebracteatus and performed comparative analysis to A. ilicifolius and three non-mangrove relatives (A. montanus, A. leucostachyus and A. mollis). Both plastome sequences and structure are highly conserved between the two mangrove species, while less similar between mangrove and non-mangrove species. Phylogenetic analysis suggested that the mangrove species were divergent from the non-mangrove groups at approximately 15.15 million years ago (Mya), where early to middle Miocene global warming and high sea level might act as one of the main forces driving the mangrove lineage entering into intertidal environments. Furthermore, 12 single nucleotide polymorphisms (SNPs) and 10 insertions/deletions (indels) were detected between the plastomes of A. ebracteatus and A. ilicifolius. PCR validation further demonstrated the effectiveness of the plastid marker in distinguishing the two sibling mangrove species. Taken together, our study broadens the understanding of the origin and evolution of Acanthus mangrove plants, and provided valuable information on the correct identification and protection of endangered mangrove species A. ebracteatus.

Comparison of organelle genomes between endangered mangrove plant Dolichandrone spathacea to terrestrial relative provides insights into its origin and adaptative evolution

Dolichandrone spathacea is a mangrove associate with high medicinal and ecological values. However, due to the dual-pressure of climate change and human activities, D. spathacea has become endangered in China. Moreover, misidentification between D. spathacea and its terrestrial relative D. cauda-felina poses further challenges to field protection and proper medicinal usage of D. spathacea. Thus, to address these problems, we sequenced and assembled mitochondrial (mt) and chloroplast (cp) genomes for both D. spathacea and D. cauda-felina. Comparative analysis revealed apparently different size and scaffold number between the two mt genomes, but a high similarity between the cp genomes. Eight regions with high sequence divergence were identified between the two cp genomes, which might be used for developing candidate DNA markers for distinguishing the two species. The splitting between D. spathacea and D. cauda-felina was inferred to occur at ~6.8 - 7.7 million years ago (Mya), which may be driven by the environment fluctuations in late Miocene. In the cp genome, 12 genes related to the expression of photosynthesis-associated proteins were detected with signatures of positive selection, which may contribute to the origin and evolutionary adaptation of Dolichandrone mangrove species. These new findings do not only enrich organelle genomic resources of Dolichandrone species, but also provide important genetic clues for improving the conservation and proper usage of endangered mangrove associate D. spathacea.

The complete chloroplast genome of Taraxacum albidum (Asteraceae), a Japanese endemic dandelion

Taraxacum albidum, a perennial herb of the Asteraceae family, exhibits both tetraploid and pentaploid in Japan. This study sequenced and characterized the complete chloroplast genome of T. albidum, revealing a 151,451 bp sequence with a typical quadripartite structure, comprising one large single-copy (LSC) region of 84,052 bp, one small single-copy (SSC) region of 18,541 bp, and two inverted repeat (IR) regions, IRa and IRb, each 24,429 bp in length. The chloroplast genome, excluding duplicates, contained 113 unique genes, including 79 protein-coding genes, 30 transfer RNA genes, and four ribosomal RNA genes. The GC content of this genome was 37.7%. Phylogenetic analysis revealed that T. albidum is most closely related to T. mongolicum, with the chloroplast genome sequences being nearly identical, differing by only one nucleotide. These findings suggest that the maternal lineage of T. albidum likely originates from T. mongolicum or its closely related species.

Apomixis and the paradox of sex in plants

BACKGROUND

The predominance of sex in eukaryotes, despite the high costs of meiosis and mating, remains an evolutionary enigma. Many theories have been proposed, none of them being conclusive on its own, and they are, in part, not well applicable to land plants. Sexual reproduction is obligate in embryophytes for the great majority of species.

SCOPE

This review compares the main forms of sexual and asexual reproduction in ferns and angiosperms, based on the generation cycling of sporophyte and gametophyte (leaving vegetative propagation aside). The benefits of sexual reproduction for maintenance of genomic integrity in comparison to asexuality are discussed in the light of developmental, evolutionary, genetic and phylogenetic studies.

CONCLUSIONS

Asexual reproduction represents modifications of the sexual pathway, with various forms of facultative sexuality. For sexual land plants, meiosis provides direct DNA repair mechanisms for oxidative damage in reproductive tissues. The ploidy alternations of meiosis-syngamy cycles and prolonged multicellular stages in the haploid phase in the gametophytes provide a high efficiency of purifying selection against recessive deleterious mutations. Asexual lineages might buffer effects of such mutations via polyploidy and can purge the mutational load via facultative sexuality. The role of organelle-nuclear genome compatibility for maintenance of genome integrity is not well understood. In plants in general, the costs of mating are low because of predominant hermaphroditism. Phylogenetic patterns in the archaeplastid clade suggest that high frequencies of sexuality in land plants are concomitant with a stepwise increase of intrinsic and extrinsic stress factors. Furthermore, expansion of genome size in land plants would increase the potential mutational load. Sexual reproduction appears to be essential for keeping long-term genomic integrity, and only rare combinations of extrinsic and intrinsic factors allow for shifts to asexuality.

Insights into comparative genomics, structural features, and phylogenetic relationship of species from Eurasian Aster and its related genera (Asteraceae: Astereae) based on complete chloroplast genome

Aster L. is an economically and phylogenetically important genus in the tribe Astereae. Here, the complete plastomes of the eight Aster species were assembled and characterized using next-generation sequencing datasets. The results indicated the complete plastomes of Aster had a quadripartite structure. These genomes were 152,045-152,729 bp in length and contained 132-133 genes, including 87 protein-coding genes, 37-38 tRNA genes, and eight rRNA genes. Expansion or contraction of inverted repeat regions and forward, palindromic, complement, and reverse repeats were detected in the eight Aster species. Additionally, our analyses showed the richest type of simple sequence repeats was A/T mononucleotides, and 14 highly variable regions were discovered by analyzing the border regions, sequence divergence, and hotspots. Phylogenetic analyses indicated that 27 species in Astereae were clustered into six clades, i.e., A to D, North American, and outgroup clades, and supported that the genera Heteropappus, Kalimeris, and Heteroplexis are nested within Aster. The results indicated the clades B to D might be considered as genera. Divergence time estimate showed the clades A, B, C, and D diverged at 23.15 Mya, 15.13 Mya, 24.29 Mya, and 21.66 Mya, respectively. These results shed light on the phylogenetic relationships of Aster and provided new information on species identification of Aster and its related genera.

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.

Taxonomy and distribution of Taraxacum sect. Erythrosperma (Asteraceae) in Poland

The dandelions from Taraxacumsect.Erythrosperma are taxonomically well distinguished and ecologically restricted to warm and sunlit habitats of steppes, dry and sandy grasslands, and distributed in temperate regions of Europe and Central Asia, with some being introduced to North America. Despite the long tradition of botanical research, the taxonomy and distribution of dandelions of T.sect.Erythrosperma is still underexplored in central Europe. In this paper, by combining traditional taxonomic studies supported by micromorphological, molecular and flow cytometry analyses as well as potential distribution modelling we shed light on taxonomical and phylogenetical relationships between members of T.sect.Erythrosperma in Poland. We also provide an identification key, species-checklist, detailed descriptions of morphology and occupated habitats as well as distribution maps for 14 Polish erythrosperms (T.bellicum, T.brachyglossum, T.cristatum, T.danubium, T.disseminatum, T.dissimile, T.lacistophyllum, T.parnassicum, T.plumbeum, T.proximum, T.sandomiriense, T.scanicum, T.tenuilobum, T.tortilobum). Finally, conservation assessments performed using the IUCN method and threat categories for all the examined species are proposed.

Complete Chloroplast Genome Sequence of Sonchus brachyotus Helps to Elucidate Evolutionary Relationships with Related Species of Asteraceae

Sonchus brachyotus DC. possesses both edible and medicinal properties and is widely distributed throughout China. In this study, the complete cp genome of S. brachyotus was sequenced and assembled. The total length of the complete S. brachyotus cp genome was 151,977 bp, including an LSC region of 84,553 bp, SSC region of 18,138 bp, and IR region of 24,643 bp. Sequence analyses revealed that the cp genome encoded 132 genes, including 87 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. The GC content was 37.6%. One hundred mononucleotide microsatellites, 4 dinucleotide microsatellites, 67 trinucleotide microsatellites, 4 tetranucleotide microsatellites, and 1 long repeat were identified. The SSR frequency of the LSC region was significantly greater than that of the IR and SSC regions. In total, 175 SSRs and highly variable regions were recognized as potential cp markers. By analyzing the IR/LSC and IR/SSC boundaries, structural differences between S. brachyotus and 6 other species were detected. According to phylogenetic analyses, S. brachyotus was most closely related to S. arvensis and S. oleraceus. Overall, this study provides complete cp genome resources for S. brachyotus that will be beneficial for identifying potential molecular markers and evolutionary patterns of S. brachyotus and its closely related species.

Comparative analyses of Mikania (Asteraceae: Eupatorieae) plastomes and impact of data partitioning and inference methods on phylogenetic relationships

We assembled new plastomes of 19 species of Mikania and of Ageratina fastigiata, Litothamnus nitidus, and Stevia collina, all belonging to tribe Eupatorieae (Asteraceae). We analyzed the structure and content of the assembled plastomes and used the newly generated sequences to infer phylogenetic relationships and study the effects of different data partitions and inference methods on the topologies. Most phylogenetic studies with plastomes ignore that processes like recombination and biparental inheritance can occur in this organelle, using the whole genome as a single locus. Our study sought to compare this approach with multispecies coalescent methods that assume that different parts of the genome evolve at different rates. We found that the overall gene content, structure, and orientation are very conserved in all plastomes of the studied species. As observed in other Asteraceae, the 22 plastomes assembled here contain two nested inversions in the LSC region. The plastomes show similar length and the same gene content. The two most variable regions within Mikania are rpl32-ndhF and rpl16-rps3, while the three genes with the highest percentage of variable sites are ycf1, rpoA, and psbT. We generated six phylogenetic trees using concatenated maximum likelihood and multispecies coalescent methods and three data partitions: coding and non-coding sequences and both combined. All trees strongly support that the sampled Mikania species form a monophyletic group, which is further subdivided into three clades. The internal relationships within each clade are sensitive to the data partitioning and inference methods employed. The trees resulting from concatenated analysis are more similar among each other than to the correspondent tree generated with the same data partition but a different method. The multispecies coalescent analysis indicate a high level of incongruence between species and gene trees. The lack of resolution and congruence among trees can be explained by the sparse sampling (~ 0.45% of the currently accepted species) and by the low number of informative characters present in the sequences. Our study sheds light into the impact of data partitioning and methods over phylogenetic resolution and brings relevant information for the study of Mikania diversity and evolution, as well as for the Asteraceae family as a whole.

Extremely low nucleotide diversity among thirty-six new chloroplast genome sequences from Aldama (Heliantheae, Asteraceae) and comparative chloroplast genomics analyses with closely related genera

Aldama (Heliantheae, Asteraceae) is a diverse genus in the sunflower family. To date, nearly 200 Asteraceae chloroplast genomes have been sequenced, but the plastomes of Aldama remain undescribed. Plastomes in Asteraceae usually show little sequence divergence, consequently, our hypothesis is that species of Aldama will be overall conserved. In this study, we newly sequenced 36 plastomes of Aldama and of five species belonging to other Heliantheae genera selected as outgroups (i.e., Dimerostemma asperatum, Helianthus tuberosus, Iostephane heterophylla, Pappobolus lanatus var. lanatus, and Tithonia diversifolia). We analyzed the structure and gene content of the assembled plastomes and performed comparative analyses within Aldama and with other closely related genera. As expected, Aldama plastomes are very conserved, with the overall gene content and orientation being similar in all studied species. The length of the plastome is also consistent and the junction between regions usually contain the same genes and have similar lengths. A large ∼20 kb and a small ∼3 kb inversion were detected in the Large Single Copy (LSC) regions of all assembled plastomes, similarly to other Asteraceae species. The nucleotide diversity is very low, with only 1,509 variable sites in 127,466 bp (i.e., 1.18% of the sites in the alignment of 36 Aldama plastomes, with one of the IRs removed, is variable). Only one gene, rbcL, shows signatures of positive selection. The plastomes of the selected outgroups feature a similar gene content and structure compared to Aldama and also present the two inversions in the LSC region. Deletions of different lengths were observed in the gene ycf2. Multiple SSRs were identified for the sequenced Aldama and outgroups. The phylogenetic analysis shows that Aldama is not monophyletic due to the position of the Mexican species A. dentata. All Brazilian species form a strongly supported clade. Our results bring new understandings into the evolution and diversity of plastomes at the species level.

Comparative analysis and implications of the chloroplast genomes of three thistles (Carduus L., Asteraceae)

BACKGROUND

Carduus, commonly known as plumeless thistles, is a genus in the Asteraceae family that exhibits both medicinal value and invasive tendencies. However, the genomic data of Carduus (i.e., complete chloroplast genomes) have not been sequenced.

METHODS

We sequenced and assembled the chloroplast genome (cpDNA) sequences of three Carduus species using the Illumina Miseq sequencing system and Geneious Prime. Phylogenetic relationships between Carduus and related taxa were reconstructed using Maximum Likelihood and Bayesian Inference analyses. In addition, we used a single nucleotide polymorphism (SNP) in the protein coding region of the matK gene to develop molecular markers to distinguish C. crispus from C. acanthoides and C. tenuiflorus.

RESULTS

The cpDNA sequences of C. crispus, C. acanthoides, and C. tenuiflorus ranged from 152,342 bp to 152,617 bp in length. Comparative genomic analysis revealed high conservation in terms of gene content (including 80 protein-coding, 30 tRNA, and four rRNA genes) and gene order within the three focal species and members of subfamily Carduoideae. Despite their high similarity, the three species differed with respect to the number and content of repeats in the chloroplast genome. Additionally, eight hotspot regions, including psbI-trnS_GCU, trnE_UUC-rpoB, trnR_UCU-trnG_UCC, psbC-trnS_UGA, trnT_UGU-trnL_UAA, psbT-psbN, petD-rpoA, and rpl16-rps3, were identified in the study species. Phylogenetic analyses inferred from 78 protein-coding and non-coding regions indicated that Carduus is polyphyletic, suggesting the need for additional studies to reconstruct relationships between thistles and related taxa. Based on a SNP in matK, we successfully developed a molecular marker and protocol for distinguishing C. crispus from the other two focal species. Our study provides preliminary chloroplast genome data for further studies on plastid genome evolution, phylogeny, and development of species-level markers in Carduus.

Genomic Analyses of Phenotypic Differences Between Native and Invasive Populations of Diffuse Knapweed (Centaurea diffusa)

Invasive species represent excellent opportunities to study the evolutionary potential of traits important to success in novel environments. Although some ecologically important traits have been identified in invasive species, little is typically known about the genetic mechanisms that underlie invasion success in non-model species. Here, we use a genome-wide association (GWAS) approach to identify the genetic basis of trait variation in the non-model, invasive, diffuse knapweed [Centaurea diffusa Lam. (Asteraceae)]. To assist with this analysis, we have assembled the first draft genome reference and fully annotated plastome assembly for this species, and one of the first from this large, weedy, genus, which is of major ecological and economic importance. We collected phenotype data from 372 individuals from four native and four invasive populations of C. diffusa grown in a common environment. Using these individuals, we produced reduced-representation genotype-by-sequencing (GBS) libraries and identified 7,058 SNPs. We identify two SNPs associated with leaf width in these populations, a trait which significantly varies between native and invasive populations. In this rosette forming species, increased leaf width is a major component of increased biomass, a common trait in invasive plants correlated with increased fitness. Finally, we use annotations from Arabidopsis thaliana to identify 98 candidate genes that are near the associated SNPs and highlight several good candidates for leaf width variation.

Assessing the genetic and chemical diversity of Taraxacum species in the Korean Peninsula

The genetic relationship between Taraxacum species, also known as the dandelion, is complicated because of asexual and mixed sexual apomictic reproduction. The usage of Taraxacum species in traditional medicines make their specialized metabolism important, but interspecific chemical difference has rarely been reported for the genus. In this study, we assembled the chloroplast genome and 45S rDNA of six Taraxacum species that occur in Korea (T. campylodes, T. coreanum, T. erythrospermum, T. mongolicum, T. platycarpum, and T. ussuriense), and performed a comparative analysis, which revealed their phylogenetic relationships and possible natural hybridity. We also performed a liquid chromatography-mass spectrometry-based phytochemical analysis to reveal interspecific chemical diversity. The comparative metabolomics analysis revealed that Taraxacum species could be separated into three chemotypes according to their major defensive specialized metabolites, which were the sesquiterpene lactones, the phenolic inositols, and chlorogenic acid derivatives. The CP DNA- and 45S rDNA-based phylogenetic trees showed a tangled relationship, which supports the notion of ongoing hybridization of wild Taraxacum species. The untargeted LC-MS analysis revealed that each Taraxacum plant exhibits species-specific defensive specialized metabolism. Moreover, 45S rDNA-based phylogenetic tree correlated with the hierarchical cluster relied on metabolite compositions. Given the coincidence between these analyses, we represented that 45S rDNA could well reflect overall nuclear genome variation in Taraxacum species.

Comparative chloroplast genome analyses of Avena: insights into evolutionary dynamics and phylogeny

BACKGROUND

Oat (Avena sativa L.) is a recognized health-food, and the contributions of its different candidate A-genome progenitor species remain inconclusive. Here, we report chloroplast genome sequences of eleven Avena species, to examine the plastome evolutionary dynamics and analyze phylogenetic relationships between oat and its congeneric wild related species.

RESULTS

The chloroplast genomes of eleven Avena species (size range of 135,889-135,998 bp) share quadripartite structure, comprising of a large single copy (LSC; 80,014-80,132 bp), a small single copy (SSC; 12,575-12,679 bp) and a pair of inverted repeats (IRs; 21,603-21,614 bp). The plastomes contain 131 genes including 84 protein-coding genes, eight ribosomal RNAs and 39 transfer RNAs. The nucleotide sequence diversities (Pi values) range from 0.0036 (rps19) to 0.0093 (rpl32) for ten most polymorphic genes and from 0.0084 (psbH-petB) to 0.0240 (petG-trnW-CCA) for ten most polymorphic intergenic regions. Gene selective pressure analysis shows that all protein-coding genes have been under purifying selection. The adjacent position relationships between tandem repeats, insertions/deletions and single nucleotide polymorphisms support the evolutionary importance of tandem repeats in causing plastome mutations in Avena. Phylogenomic analyses, based on the complete plastome sequences and the LSC intermolecular recombination sequences, support the monophyly of Avena with two clades in the genus.

CONCLUSIONS

Diversification of Avena plastomes is explained by the presence of highly diverse genes and intergenic regions, LSC intermolecular recombination, and the co-occurrence of tandem repeat and indels or single nucleotide polymorphisms. The study demonstrates that the A-genome diploid-polyploid lineage maintains two subclades derived from different maternal ancestors, with A. longiglumis as the first diverging species in clade I. These genome resources will be helpful in elucidating the chloroplast genome structure, understanding the evolutionary dynamics at genus Avena and family Poaceae levels, and are potentially useful to exploit plastome variation in making hybrids for plant breeding.

Species Identification of Dracaena Using the Complete Chloroplast Genome as a Super-Barcode

The taxonomy and nomenclature of Dracaena plants are much disputed, particularly for several Dracaena species in Asia. However, neither morphological features nor common DNA regions are ideal for identification of Dracaena spp. Meanwhile, although multiple Dracaena spp. are sources of the rare traditional medicine dragon's blood, the Pharmacopoeia of the People's Republic of China has defined Dracaena cochinchinensis as the only source plant. The inaccurate identification of Dracaena spp. will inevitably affect the clinical efficacy of dragon's blood. It is therefore important to find a better method to distinguish these species. Here, we report the complete chloroplast (CP) genomes of six Dracaena spp., D. cochinchinensis, D. cambodiana, D. angustifolia, D. terniflora, D. hokouensis, and D. elliptica, obtained through high-throughput Illumina sequencing. These CP genomes exhibited typical circular tetramerous structure, and their sizes ranged from 155,055 (D. elliptica) to 155,449 bp (D. cochinchinensis). The GC content of each CP genome was 37.5%. Furthermore, each CP genome contained 130 genes, including 84 protein-coding genes, 38 tRNA genes, and 8 rRNA genes. There were no potential coding or non-coding regions to distinguish these six species, but the maximum likelihood tree of the six Dracaena spp. and other related species revealed that the whole CP genome can be used as a super-barcode to identify these Dracaena spp. This study provides not only invaluable data for species identification and safe medical application of Dracaena but also an important reference and foundation for species identification and phylogeny of Liliaceae plants.

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.

Comparative genomics of 11 complete chloroplast genomes of Senecioneae (Asteraceae) species: DNA barcodes and phylogenetics

BACKGROUND

Majority of the species within Senecioneae are classified in Senecio, making it the tribe's largest genus. Certain intergeneric relationships within the tribe are vaguely defined, with the genus Senecio being partly linked to this ambiguity. Infrageneric relationships within Senecio remain largely unknown and consequently, the genus has undergone continuous expansion and contraction over the recent past due to addition and removal of taxa. Dendrosenecio, an endemic genus in Africa, is one of its segregate genera. To heighten the understanding of species divergence and phylogeny within the tribe, the complete chloroplast genomes of the first five Senecio and six Dendrosenecio species were sequenced and analyzed in this study.

RESULTS

The entire length of the complete chloroplast genomes was ~ 150 kb and ~ 151 kb in Dendrosenecio and Senecio respectively. Characterization of the 11 chloroplast genomes revealed a significant degree of similarity particularly in their organization, gene content, repetitive sequence composition and patterns of codon usage. The chloroplast genomes encoded an equal number of unique genes out of which 80 were protein-coding genes, 30 transfer ribonucleic acid, and four ribosomal ribonucleic acid genes. Based on comparative sequence analyses, the level of divergence was lower in Dendrosenecio. A total of 331 and 340 microsatellites were detected in Senecio and Dendrosenecio, respectively. Out of which, 25 and five chloroplast microsatellites (cpSSR) were identified as potentially valuable molecular markers. Also, through whole chloroplast genome comparisons and DNA polymorphism tests, ten divergent hotspots were identified. Potential primers were designed creating genomic tools to further molecular studies within the tribe. Intergeneric relationships within the tribe were firmly resolved using genome-scale dataset in partitioned and unpartitioned schemes. Two main clades, corresponding to two subtribes within the Senecioneae, were formed with the genus Ligularia forming a single clade while the other had Dendrosenecio, Pericallis, Senecio and Jacobaea. A sister relationship was revealed between Dendrosenecio and Pericallis whereas Senecio, and Jacobaea were closely placed in a different clade.

CONCLUSION

Besides emphasizing on the potential of chloroplast genome data in resolving intergeneric relationships within Senecioneae, this study provides genomic resources to facilitate species identification and phylogenetic reconstructions within the respective genera.

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.

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.

New chromosome counts and genome size estimates for 28 species of Taraxacum sect. Taraxacum

The species-rich and widespread genus Taraxacum F. H. Wiggers, 1780 (Asteraceae subfamily Cichorioideae) is one of the most taxonomically complex plant genera in the world, mainly due to its combination of different sexual and asexual reproduction strategies. Polyploidy is usually confined to apomictic microspecies, varying from 3x to 6x (rarely 10x). In this study, we focused on Taraxacum sect. Taraxacum (= T.sect.Ruderalia; T.officinale group), i.e., the largest group within the genus. We counted chromosome numbers and measured the DNA content for species sampled in Central Europe, mainly in Czechia. The chromosome number of the 28 species (T.aberrans Hagendijk, Soest & Zevenbergen, 1974, T.atroviride Štěpánek & Trávníček, 2008, T.atrox Kirschner & Štěpánek, 1997, T.baeckiiforme Sahlin, 1971, T.chrysophaenum Railonsala, 1957, T.coartatum G.E. Haglund, 1942, T.corynodes G.E. Haglund, 1943, T.crassum H. Øllgaard & Trávníček, 2003, T.deltoidifrons H. Øllgaard, 2003, T.diastematicum Marklund, 1940, T.gesticulans H. Øllgaard, 1978, T.glossodon Sonck & H. Øllgaard, 1999, T.guttigestans H. Øllgaard in Kirschner & Štěpánek, 1992, T.huelphersianum G.E. Haglund, 1935, T.ingens Palmgren, 1910, T.jugiferum H. Øllgaard, 2003, T.laticordatum Marklund, 1938, T.lojoense H. Lindberg, 1944 (= T.debrayi Hagendijk, Soest & Zevenbergen, 1972, T.lippertianum Sahlin, 1979), T.lucidifrons Trávníček, ineditus, T.obtusifrons Marklund, 1938, T.ochrochlorum G.E. Haglund, 1942, T.ohlsenii G.E. Haglund, 1936, T.perdubium Trávníček, ineditus, T.praestabile Railonsala, 1962, T.sepulcrilobum Trávníček, ineditus, T.sertatum Kirschner, H. Øllgaard & Štěpánek, 1997, T.subhuelphersianum M.P. Christiansen, 1971, T.valens Marklund, 1938) is 2n = 3x = 24. The DNA content ranged from 2C = 2.60 pg (T.atrox) to 2C = 2.86 pg (T.perdubium), with an average value of 2C = 2.72 pg. Chromosome numbers are reported for the first time for 26 species (all but T.diastematicum and T.obtusifrons), and genome size estimates for 26 species are now published for the first time.

Species Identification of Conyza bonariensis Assisted by Chloroplast Genome Sequencing

Flaxleaf fleabane (Conyza bonariensis [L.] Cronquist) is one of the most difficult weeds to control worldwide. There are more than 150 Conyza species in the world and eight species in Australia. Correct identification of these species can be problematic due to their morphological similarities especially at seedling stage. Developing a robust genetics - based species identification method to distinguish C. bonariensis from other closely related species is important for early control of weeds. We thus examined the chloroplast (cp) genome of C. bonariensis, aiming to identify novel DNA barcodes from the genome sequences, and use the entire cp genome as a super-barcode for molecular identification. The C. bonariensis chloroplast genome is 152,076 bp in size, encodes 133 genes including 88 protein-coding genes, 37 tRNA genes and 8 ribosomal RNA genes. A total of 151 intergenic regions and 19 simple sequence repeats were identified in the cp genome of C. bonariensis, which provides a useful genetic resource to develop robust markers for the genetic diversity studies of Conyza species. The sequence information was used to design a robust DNA barcode rps16 and trnQ-UUG which successfully separated three predominant Conyza species (C. bonariensis, C. canadensis, and C. sumatrensis). Phylogenetic analyses based on the cp genomes of C. bonariensis, C. canadensis and 18 other Asteraceae species revealed the potential of using entire cp genome as a plant super-barcode to distinguish closely-related weed species.

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.

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.