Phylogenetic relationships, selective pressure and molecular markers development of six species in subfamily Polygonoideae based on complete chloroplast genomes.

The subfamily Polygonoideae encompasses a diverse array of medicinal and horticultural plants that hold significant economic value. However, due to the lack of a robust taxonomy based on phylogenetic relationships, the classification within this family is perplexing, and there is also a scarcity of reports on the chloroplast genomes of many plants falling under this classification. In this study, we conducted a comprehensive analysis by sequencing and characterizing the complete chloroplast genomes of six Polygonoideae plants, namely Pteroxygonum denticulatum, Pleuropterus multiflorus, Pleuropterus ciliinervis, Fallopia aubertii, Fallopia dentatoalata, and Fallopia convolvulus. Our findings revealed that these six plants possess chloroplast genomes with a typical quadripartite structure, averaging 162,931 bp in length. Comparative chloroplast analysis, codon usage analysis, and repetitive sequence analysis demonstrated a high level of conservation within the chloroplast genomes of these plants. Furthermore, phylogenetic analysis unveiled a distinct clade occupied by P. denticulatum, while P. ciliinrvis displayed a closer relationship to the three plants belonging to the Fallopia genus. Selective pressure analysis based on maximum likelihood trees showed that a total of 14 protein-coding genes exhibited positive selection, with psbB and ycf1 having the highest number of positive amino acid sites. Additionally, we identified four molecular markers, namely petN-psbM, psal-ycf4, ycf3-trnS-GGA, and trnL-UAG-ccsA, which exhibit high variability and can be utilized for the identification of these six plants.

Effective dispersal patterns in prairie plant species across human-modified landscapes.

Effective dispersal among plant populations is dependent on vector behaviour, landscape features and availability of adequate habitats. To capture landscape feature effects on dispersal, studies must be conducted at scales reflecting single-generation dispersal events (mesoscale). Many studies are conducted at large scales where genetic differentiation is due to dispersal occurring over multiple generations, making it difficult to interpret the effects of specific landscape features on vector behaviour. Genetic structure at the mesoscale may be determined by ecological and evolutionary processes, such as the consequences of vector behaviour on patterns of gene flow. We used chloroplast haplotypes and nuclear genome SNP surveys to identify landscape features influencing seed and pollen dispersal at a mesoscale within the Rogue River Valley in southern Oregon. We evaluated biotic and abiotic vector behaviour by contrasting two annual species with differing dispersal mechanisms; Achyrachaena mollis (Asteraceae) is a self-pollinating and anemochoric species, and Plectritis congesta (Caprifoliaceae) is biotically pollinated with barochoric seeds. Using landscape genetics methods, we identified features of the study region that conduct or restrict dispersal. We found chloroplast haplotypes were indicative of historic patterns of gene flow prior to human modification of landscapes. Seed dispersal of A. mollis was best supported by models of isolation by distance, while seed-driven gene flow of P. congesta was determined by the distribution of preserved natural spaces and quality habitat. Nuclear genetic structure was driven by both pollen and seed dispersal, and both species responded to contemporary landscape changes, such as urban and agricultural conversion, and habitat availability.

Genome-wide DNA methylation dynamics following recent polyploidy in the allotetraploid Tragopogon miscellus (Asteraceae).

Polyploidy is an important evolutionary force, yet epigenetic mechanisms, such as DNA methylation, that regulate genome-wide expression of duplicated genes remain largely unknown. Here, we use Tragopogon (Asteraceae) as a model system to discover patterns and temporal dynamics of DNA methylation in recently formed polyploids. The naturally occurring allotetraploid Tragopogon miscellus formed in the last 95-100 yr from parental diploids Tragopogon dubius and T. pratensis. We profiled the DNA methylomes of these three species using whole-genome bisulfite sequencing. Genome-wide methylation levels in T. miscellus were intermediate between its diploid parents. However, nonadditive CG and CHG methylation occurred in transposable elements (TEs), with variation among TE types. Most differentially methylated regions (DMRs) showed parental legacy, but some novel DMRs were detected in the polyploid. Differentially methylated genes (DMGs) were also identified and characterized. This study provides the first assessment of both overall and locus-specific patterns of DNA methylation in a recent natural allopolyploid and shows that novel methylation variants can be generated rapidly after polyploid formation. Together, these results demonstrate that mechanisms to regulate duplicate gene expression may arise soon after allopolyploid formation and that these mechanisms vary among genes.

High-Quality Assembly and Analysis of the Complete Mitogenomes of German Chamomile (Matricaria recutita) and Roman Chamomile (Chamaemelum nobile).

German chamomile (Matricaria chamomilla L.) and Roman chamomile (Chamaemelum nobile) are the two well-known chamomile species from the Asteraceae family. Owing to their essential oils and higher medicinal value, these have been cultivated widely across Europe, Northwest Asia, North America, and Africa. Regarding medicinal applications, German chamomile is the most commonly utilized variety and is frequently recognized as the "star among medicinal species". The insufficient availability of genomic resources may negatively impact the progression of chamomile industrialization. Chamomile's mitochondrial genome is lacking in extensive empirical research. In this study, we achieved the successful sequencing and assembly of the complete mitochondrial genome of M. chamomilla and C. nobile for the first time. An analysis was conducted on codon usage, sequence repeats within the mitochondrial genome of M. chamomilla and C. nobile. The phylogenetic analysis revealed a consistent positioning of M. chamomilla and C. nobile branches within both mitochondrial and plastid-sequence-based phylogenetic trees. Furthermore, the phylogenetic analysis also showed a close relationship between M. chamomilla and C. nobile within the clade comprising species from the Asteraceae family. The results of our analyses provide valuable resources for evolutionary research and molecular barcoding in chamomile.

Phylogenomics and morphology of Celmisiinae (Asteraceae: Astereae): Taxonomic and evolutionary implications.

The tribe Astereae (Asteraceae) includes 36 subtribes and 252 genera, and is distributed worldwide in temperate and tropical regions. One of the subtribes, Celmisiinae Saldivia, has been recently circumscribed to include six genera and ca. 160 species, and is restricted to eastern Australia, New Zealand, and New Guinea. The species show an impressive range of growth habit, from small herbs and ericoid subshrubs to medium-sized trees. They live in a wide range of habitats and are often dominant in subalpine and alpine vegetation. Despite the well-supported circumscription of Celmisiinae, uncertainties have remained about their internal relationships and classification at genus and species levels. This study exploited recent advances in high-throughput sequencing to build a robust multi-gene phylogeny for the subtribe Celmisiinae. The target enrichment Angiosperms353 bait set and the hybpiper-nf and paragone-nf pipelines were used to retrieve, infer, and assemble orthologous loci from 75 taxa representing all the main putative clades within the subtribe. Because of the diploidised ploidy level in Celmisiinae, as well as missing data in the assemblies, uncertainty remains surrounding the inference of orthology detection. However, based on a variety of gene-family sets, coalescent and concatenation-based phylogenetic reconstructions recovered similar topologies. Paralogy and missing data in the gene-families caused some problems, but the estimated phylogenies were well-supported and well-resolved. The phylogenomic evidence supported Celmisiinae and three main clades: the Pleurophyllum clade (Pleurophyllum, Macrolearia and Damnamenia), mostly in the New Zealand Subantarctic Islands, Celmisia of mainland New Zealand and Australia, and Shawia (including 'Olearia pro parte' and Pachystegia) of New Zealand, Australia and New Guinea. The results presented here add to the accumulating support for the Angiosperms353 bait set as an efficient method for documenting plant diversity.

Overexpression of tocopherol biosynthesis genes in guayule (Parthenium argentatum) reduces rubber, resin and argentatins content in stem and leaf tissues.

Natural rubber produced in stems of the guayule plant (Parthenium argentatum) is susceptible to post-harvest degradation from microbial or thermo-oxidative processes, especially once stems are chipped. As a result, the time from harvest to extraction must be minimized to recover high quality rubber, especially in warm summer months. Tocopherols are natural antioxidants produced in plants through the shikimate and methyl-erythtiol-4-phosphate (MEP) pathways. We hypothesized that increased in vivo guayule tocopherol content might protect rubber from post-harvest degradation, and/or allow reduced use of chemical antioxidants during the extraction process. With the objective of enhancing tocopherol content in guayule, we overexpressed four Arabidopsis thaliana tocopherol pathway genes in AZ-2 guayule via Agrobacterium-mediated transformation. Tocopherol content was increased in leaf and stem tissues of most transgenic lines, and some improvement in thermo-oxidative stability was observed. Overexpression of the four tocopherol biosynthesis enzymes, however, altered other isoprenoid pathways resulting in reduced rubber, resin and argentatins content in guayule stems. The latter molecules are mainly synthesized from precursors derived from the mevalonate (MVA) pathway. Our results suggest the existence of crosstalk between the MEP and MVA pathways in guayule and the possibility that carbon metabolism through the MEP pathway impacts rubber biosynthesis.

Development and evolution of the Asteraceae capitulum.

Asteraceae represent one of the largest and most diverse families of plants. The evolutionary success of this family has largely been contributed to their unique inflorescences, capitula that mimic solitary flowers but are typically aggregates of multiple florets. Here, we summarize the recent molecular and genetic level studies that have promoted our understanding of the development and evolution of capitula. We focus on new results on patterning of the enlarged meristem resulting in the iconic phyllotactic arrangement of florets in Fibonacci numbers of spirals. We also summarize the current understanding of the genetic networks regulating the characteristic reproductive traits in the family such as floral dimorphism and differentiation of highly specialized floral organs. So far, developmental studies in Asteraceae are still limited to a very narrow selection of model species. Along with the recent advancements in genomics and phylogenomics, Asteraceae and its relatives provide an outstanding model clade for extended evo-devo studies to exploit the morphological diversity and the underlying molecular networks and to translate this knowledge to the breeding of the key crops in the family.

The complete plastome sequences of invasive weed Parthenium hysterophorus: genome organization, evolutionary significance, structural features, and comparative analysis.

Parthenium hysterophorus, a globally widespread weed, poses a significant threat to agricultural ecosystems due to its invasive nature. We investigated the chloroplast genome of P. hysterophorus in this study. Our analysis revealed that the chloroplast genome of P. hysterophorus spans a length of 151,881 base pairs (bp). It exhibits typical quadripartite structure commonly found in chloroplast genomes, including inverted repeat regions (IR) of 25,085 bp, a small single copy (SSC) region of 18,052 bp, and a large single copy (LSC) region of 83,588 bp. A total of 129 unique genes were identified in P. hysterophorus chloroplast genomes, including 85 protein-coding genes, 36 tRNAs, and eight rRNAs genes. Comparative analysis of the P. hysterophorus plastome with those of related species from the tribe Heliantheae revealed both conserved structures and intriguing variations. While many structural elements were shared among the species, we identified a rearrangement in the large single-copy region of P. hysterophorus. Moreover, our study highlighted notable gene divergence in several specific genes, namely matK, ndhF, clpP, rps16, ndhA, rps3, and ndhD. Phylogenetic analysis based on the 72 shared genes placed P. hysterophorus in a distinct clade alongside another species, P. argentatum. Additionally, the estimated divergence time between the Parthenium genus and Helianthus (sunflowers) was approximately 15.1 million years ago (Mya). These findings provide valuable insights into the evolutionary history and genetic relationships of P. hysterophorus, shedding light on its divergence and adaptation over time.

Plastid genomes provide insights into the phylogeny and chloroplast evolution of the paper daisy tribe Gnaphalieae (Asteraceae).

Chloroplast genomes, as an essential source of phylogenetic information, are increasingly utilized in the evolutionary study of angiosperms. Gnaphalieae is a medium-sized tribe of the sunflower family of Asteraceae, with about 2,100 species in 178 genera distributed in temperate habitats worldwide. There has been considerable progress in our understanding of their phylogenetic evolution using both nuclear and chloroplast sequences, but no focus on chloroplast genomic data. In this study, we performed sequencing, assembly, and annotation of 16 representative chloroplast genomes from all the major lineages of Gnaphalieae. Our results showed that the plastomes exhibited a typical circular tetrad structure with similar genomic structure gene content. But there were differences in genome size, SSRs, and codon usage within the tribe. Phylogenetic analysis revealed Relhania clade is the earliest diverged lineages with the Lasiopogon clade and the Gnaphalium s.s. clade diverged subsequently. The core group includes FLAG clade sister to the HAP and Australasian group. Compared with the outgroup species, chloroplast genome size of the FLAG clade is much reduced whereas those of Australasian, HAP, Gnaphalium s.s., Lasiopogon and Relhania clades are relatively expanded. Insertions and deletions in the intergenic regions associated with repetitive sequence variations are supposed to be the main factor leading to length variations in the chloroplast genomes of Gnaphalieae. The comparative analyses of chloroplast genomes would provide useful implications into understanding the taxonomic and evolutionary history of Gnaphalieae.

Investigation of regulatory divergence between homoeologs in the recently formed allopolyploids, Tragopogon mirus and T. miscellus (Asteraceae).

Polyploidy is an important evolutionary process throughout eukaryotes, particularly in flowering plants. Duplicated gene pairs (homoeologs) in allopolyploids provide additional genetic resources for changes in molecular, biochemical, and physiological mechanisms that result in evolutionary novelty. Therefore, understanding how divergent genomes and their regulatory networks reconcile is vital for unraveling the role of polyploidy in plant evolution. Here, we compared the leaf transcriptomes of recently formed natural allotetraploids (Tragopogon mirus and T. miscellus) and their diploid parents (T. porrifolius X T. dubius and T. pratensis X T. dubius, respectively). Analysis of 35 400 expressed loci showed a significantly higher level of transcriptomic additivity compared to old polyploids; only 22% were non-additively expressed in the polyploids, with 5.9% exhibiting transgressive expression (lower or higher expression in the polyploids than in the diploid parents). Among approximately 7400 common orthologous regions (COREs), most loci in both allopolyploids exhibited expression patterns that were vertically inherited from their diploid parents. However, 18% and 20.3% of the loci showed novel expression bias patterns in T. mirus and T. miscellus, respectively. The expression changes of 1500 COREs were explained by cis-regulatory divergence (the condition in which the two parental subgenomes do not interact) between the diploid parents, whereas only about 423 and 461 of the gene expression changes represent trans-effects (the two parental subgenomes interact) in T. mirus and T. miscellus, respectively. The low degree of both non-additivity and trans-effects on gene expression may present the ongoing evolutionary processes of the newly formed Tragopogon polyploids (~80-90 years).

Comparative Physiological and Transcriptome Analysis of Crossostephium chinense Reveals Its Molecular Mechanisms of Salt Tolerance.

Crossostephium chinense is a wild species with strong salt tolerance that has great potential to improve the salt tolerance of cultivated chrysanthemums. Conversely, the unique salt-tolerant molecular mechanisms of Cr. chinense are still unclear. This study performed a comparative physiological and transcriptome analysis of Cr. chinense, Chrysanthemum lavandulifolium, and three hybrids to investigate the salt-tolerant molecular mechanisms of Cr. chinense. The physiological results showed that Cr. chinense maintained higher superoxide dismutase (SOD) activity, alleviating oxidative damage to the membrane. KEGG enrichment analysis showed that plant hormone signaling transduction and the MAPK signaling pathway were mostly enriched in Cr. chinense and hybrids under salt stress. Further weighted gene co-expression network analysis (WGCNA) of DEGs suggested that abscisic acid (ABA) signaling transduction may play a significant role in the salt-tolerant mechanisms of Cr. chinense and hybrids. The tissue-specific expression patterns of the candidate genes related to ABA signaling transduction and the MAPK signaling pathway indicate that genes related to ABA signaling transduction demonstrated significant expression levels under salt stress. This study offers important insights into exploring the underlying salt-tolerant mechanisms of Cr. chinense mediated by ABA signaling transduction and broadens our understanding of the breeding strategies for developing salt-tolerant cultivars utilizing salt-tolerant chrysanthemum germplasms.

Chloroplast Genomes and Phylogenetic Analysis of Three Carthamus (Asteraceae) Species.

The genus Carthamus Linnaeus, which belongs to the tribe Cardueae in the Asteraceae family, originated in the Mediterranean region and consists of approximately 20 species worldwide. Understanding the phylogeny of the Carthamus is crucial for the cultivation of C. tinctorius. Although chloroplast genomes are widely used for species identification and evolutionary studies, there have been limited investigations on the chloroplast genomes of Carthamus species. In this study, we assembled the chloroplast genomes of C. persicus, C. tinctorius × C. persicus, and C. lanatus and combined them with the five chloroplast genomes of C. tinctorius for comparative genomic analysis. The sizes of the chloroplast genomes of C. lanatus, C. persicus, and C. tinctorius × C. persicus were 152,602 bp, 153,177 bp, and 153,177 bp, respectively. Comparative analysis showed that the chloroplast genome structures of the four Carthamus species were highly conserved. Additionally, the phylogenomic analysis demonstrated that the plastid genome and angiosperms353 dataset significantly improved the phylogenetic support of Carthamus species. This analysis supported Carthamus as a monophyletic taxon and its internal division into the sect. Carthamus and sect. Atractylis. The Carthamus was closely related to Carduncellus, Femeniasia, Phonus, and Centaurea. In conclusion, this study not only expands our understanding of the cp genomes of Carthamus species but also provides support for more comprehensive phylogenetic studies of Carthamus.

Nuclear and Mitochondrial Genome Assemblies of the Beetle, Zygogramma bicolorata, a Globally Important Biocontrol Agent of Invasive Weed Parthenium hysterophorus.

Implementing a genetic-based approach to achieve the full potential of classical biocontrol programs has been advocated for decades. The availability of genome-level information brings the opportunity to scrutinize biocontrol traits for their efficacy and evolvability. However, implementation of this advocacy remains limited to few instances. Biocontrol of a globally noxious weed, Parthenium hysterophorus, by the leaf-feeding beetle, Zygogramma bicolorata, has been in place for more than four decades now, with varying levels of success. As the first step in providing genetic-based improvement to this biocontrol program, we describe the nuclear and mitochondrial assemblies of Z. bicolorata. We assembled the genome from the long-read sequence data, error corrected with high-throughput short reads and checked for contaminants and sequence duplication to produce a 936 Mb nuclear genome. With 96.5% Benchmarking Universal Single-Copy Orthologs completeness and the long terminal repeat assembly index 12.91, we present a reference-quality assembly that appeared to be repeat rich at 62.7% genome-wide and consists of 29,437 protein-coding regions. We detected signature of nuclear insertion of mitochondrial fragments in 80 nuclear positions comprising 13 kb out of 17.9 kb mitochondria genome sequence. This genome, along with its annotations, provides a valuable resource to gain further insights into the biocontrol traits of Z. bicolorata for improving the control of the invasive weed P. hysterophorus.

Development and characterization of expressed sequence tag-simple sequence repeat markers for Anaphalis margaritacea var. yedoensis (Asteraceae).

RNA-sequencing was used to develop 16 microsatellite markers for the pearly everlasting, Anaphalis margaritacea var. yedoensis (Franch. et Sav.) Ohwi (Asteraceae), which inhabits gravel bars throughout the Japanese archipelago. The mean number of alleles for these 16 markers in two populations in the Hokkaido and Shizuoka Prefectures, was 3.5 and 4.0, respectively, while the mean expected heterozygosity was 0.525 and 0.560, respectively, with a significant genetic differentiation between the two populations. All markers could also be amplified in two conspecific taxa, A. margaritacea var. margaritacea and var. angustifolia, whereas 11 of the 16 markers were amplifiable in two congeneric species, A. sinica and A. alpicola. These newly developed microsatellite markers will support understanding of population genetics and mating systems in A. margaritacea var. yedoensis, and several will potentially be of use in similar studies in other Anaphalis species.

An Overview to the Index to Chromosome Numbers in Asteraceae Database: Revisiting Base Chromosome Numbers, Polyploidy, Descending Dysploidy, and Hybridization.

A brief overview to the Index to Chromosome Numbers in Asteraceae database is provided. The database contains karyological information on Asteraceae and has been repeatedly improved and updated and is now hosted at the National Bioscience Database center. Also, we take the opportunity to revisit the evolution of base chromosome numbers in Asteraceae, emphasizing the phenomena of polyploidy, descending dysploidy, and hybridization, common in the family. Chromosome numbers for species included in one of the most recent phylogenetic treatments of the Asteraceae were obtained from the Index to Chromosome Numbers in Asteraceae database were mapped on to the modified phylogeny diagram, and base chromosome numbers were determined for each branch of the phylogeny. Results for tribal base numbers were the same as those hypothesized in our previous work with additional base numbers added for tribes not previously recognized but supported by newer phylogenetic methods. The Asteraceae show an ancestral base chromosome number of x = 9 and originated in the Antarctica (Gondowanaland) in Cretaceous (80 Mys ago). The x = 9 number has been retained through successive South American lineages of the Barnadesieeae, Gochnatieae, Stiffieae, Wunderlichieae, Astereae, and Senecioneae following northward migration. Northward migration to Africa was accompanied with x = 10 becoming the dominant base chromosome number as the family evolved multiple additional tribes. Northward migration to Australasia with x = 9 was in Astereae and the families Goodeneaseae, Menyanthaceae, and Stylydiaceae. The evolution of the North American Heliantheae alliance began with the appearance of x2 = 19 which persisted in multiple additional new tribes. Frequent dysploidy decreases, polyploidy and hybridization occurred throughout the history of the family.

Comparative genomics reveals a unique nitrogen-carbon balance system in Asteraceae.

The Asteraceae (daisy family) is one of the largest families of plants. The genetic basis for its high biodiversity and excellent adaptability has not been elucidated. Here, we compare the genomes of 29 terrestrial plant species, including two de novo chromosome-scale genome assemblies for stem lettuce, a member of Asteraceae, and Scaevola taccada, a member of Goodeniaceae that is one of the closest outgroups of Asteraceae. We show that Asteraceae originated ~80 million years ago and experienced repeated paleopolyploidization. PII, the universal regulator of nitrogen-carbon (N-C) assimilation present in almost all domains of life, has conspicuously lost across Asteraceae. Meanwhile, Asteraceae has stepwise upgraded the N-C balance system via paleopolyploidization and tandem duplications of key metabolic genes, resulting in enhanced nitrogen uptake and fatty acid biosynthesis. In addition to suggesting a molecular basis for their ecological success, the unique N-C balance system reported for Asteraceae offers a potential crop improvement strategy.

Comparative transcriptomics and gene expression divergence associated with homoploid hybrid speciation in Argyranthemum.

Ecological isolation is increasingly thought to play an important role in speciation, especially for the origin and reproductive isolation of homoploid hybrid species. However, the extent to which divergent and/or transgressive gene expression changes are involved in speciation is not well studied. In this study, we employ comparative transcriptomics to investigate gene expression changes associated with the origin and evolution of two homoploid hybrid plant species, Argyranthemum sundingii and A. lemsii (Asteraceae). As there is no standard methodology for comparative transcriptomics, we examined five different pipelines for data assembly and analysing gene expression across the four species (two hybrid and two parental). We note biases and problems with all pipelines, and the approach used affected the biological interpretation of the data. Using the approach that we found to be optimal, we identify transcripts showing DE between the parental taxa and between the homoploid hybrid species and their parents; in several cases, putative functions of these DE transcripts have a plausible role in ecological adaptation and could be the cause or consequence of ecological speciation. Although independently derived, the homoploid hybrid species have converged on similar expression phenotypes, likely due to adaptation to similar habitats.

Genome-Wide Datasets of Chicories (Cichorium intybus L.) for Marker-Assisted Crop Breeding Applications: A Systematic Review and Meta-Analysis.

Cichorium intybus L. is the most economically important species of its genus and among the most important of the Asteraceae family. In chicory, many linkage maps have been produced, several sets of mapped and unmapped markers have been developed, and dozens of genes linked to traits of agronomic interest have been investigated. This treasure trove of information, properly cataloged and organized, is of pivotal importance for the development of superior commercial products with valuable agronomic potential in terms of yield and quality, including reduced bitter taste and increased inulin production, as well as resistance or tolerance to pathogens and resilience to environmental stresses. For this reason, a systematic review was conducted based on the scientific literature published in chicory during 1980-2023. Based on the results obtained from the meta-analysis, we created two consensus maps capable of supporting marker-assisted breeding (MAB) and marker-assisted selection (MAS) programs. By taking advantage of the recently released genome of C. intybus, we built a 639 molecular marker-based consensus map collecting all the available mapped and unmapped SNP and SSR loci available for this species. In the following section, after summarizing and discussing all the genes investigated in chicory and related to traits of interest such as reproductive barriers, sesquiterpene lactone biosynthesis, inulin metabolism and stress response, we produced a second map encompassing 64 loci that could be useful for MAS purposes. With the advent of omics technologies, molecular data chaos (namely, the situation where the amount of molecular data is so complex and unmanageable that their use becomes challenging) is becoming far from a negligible issue. In this review, we have therefore tried to contribute by standardizing and organizing the molecular data produced thus far in chicory to facilitate the work of breeders.

Phylogeography of the temperate grassland plant Tephroseris kirilowii (Asteraceae) inferred from multiplexed inter-simple sequence repeat genotyping by sequencing (MIG-seq) data.

A group of temperate grassland plant species termed the "Mansen elements" occurs in Japan and is widely distributed in the grasslands of continental East Asia. It has been hypothesized that these species are continental grassland relicts in Japan that stretch back to a colder age, but their migration history has not been elucidated. To assess the migration history of the Mansen elements, we performed phylogeographic analyses of Tephroseris kirilowii, a member of this group, using single-nucleotide polymorphisms (SNPs) obtained from multiplexed inter-simple sequence repeat genotyping by sequencing (MIG-seq). It was estimated that the Japanese populations of T. kirilowii were divided from those of continental East Asia at 25.2 thousand years ago (ka) with 95% highest probability density interval (HPD) of 15.3-40.0 ka and that Japanese clades first diverged at 20.2 ka with 95% HPD of 10.4-30.1 ka. As the climatically suitable range during the last glacial maximum (LGM) estimated using ecological niche modeling (ENM) was limited in Japan and there was a slight genetic differentiation among Japanese populations, a post-glacial expansion of T. kirilowii in the Japanese Archipelago was indicated.

Multiple gene co-options underlie the rapid evolution of sexually deceptive flowers in Gorteria diffusa.

Gene co-option, the redeployment of an existing gene in an unrelated developmental context, is an important mechanism underlying the evolution of morphological novelty. In most cases described to date, novel traits emerged by co-option of a single gene or genetic network. Here, we show that the integration of multiple co-opted genetic elements facilitated the rapid evolution of complex petal spots that mimic female bee-fly pollinators in the sexually deceptive South African daisy Gorteria diffusa. First, co-option of iron homeostasis genes altered petal spot pigmentation, producing a color similar to that of female pollinators. Second, co-option of the root hair gene GdEXPA7 enabled the formation of enlarged papillate petal epidermal cells, eliciting copulation responses from male flies. Third, co-option of the miR156-GdSPL1 transcription factor module altered petal spot placement, resulting in better mimicry of female flies resting on the flower. The three genetic elements were likely co-opted sequentially, and strength of sexual deception in different G. diffusa floral forms strongly correlates with the presence of the three corresponding morphological alterations. Our findings suggest that gene co-options can combine in a modular fashion, enabling rapid evolution of novel complex traits.