Intra-leaf modeling of Cannabis leaflet shape produces leaf models that predict genetic and developmental identities.

The iconic, palmately compound leaves of Cannabis have attracted significant attention in the past. However, investigations into the genetic basis of leaf shape or its connections to phytochemical composition have yielded inconclusive results. This is partly due to prominent changes in leaflet number within a single plant during development, which has so far prevented the proper use of common morphometric techniques. Here, we present a new method that overcomes the challenge of nonhomologous landmarks in palmate, pinnate, and lobed leaves, using Cannabis as an example. We model corresponding pseudo-landmarks for each leaflet as angle-radius coordinates and model them as a function of leaflet to create continuous polynomial models, bypassing the problems associated with variable number of leaflets between leaves. We analyze 341 leaves from 24 individuals from nine Cannabis accessions. Using 3591 pseudo-landmarks in modeled leaves, we accurately predict accession identity, leaflet number, and relative node number. Intra-leaf modeling offers a rapid, cost-effective means of identifying Cannabis accessions, making it a valuable tool for future taxonomic studies, cultivar recognition, and possibly chemical content analysis and sex identification, in addition to permitting the morphometric analysis of leaves in any species with variable numbers of leaflets or lobes.

Strategies of diaspore dispersal investment in Compositae: the case of the Andean highlands.

BACKGROUND AND AIMS: Understanding diaspore morphology and how much a species invests on dispersal appendages is key for improving our knowledge of dispersal in fragmented habitats. We investigate diaspore morphological traits in high-Andean Compositae and their main abiotic and biotic drivers and test whether they play a role in species distribution patterns across the naturally fragmented high-Andean grasslands. METHODS: We collected diaspore trait data for 125 Compositae species across 47 tropical high-Andean summits, focusing on achene length and pappus-to-achene length ratio, with the latter as a proxy of dispersal investment. We analysed the role of abiotic (temperature, elevation and latitude) and biotic factors (phylogenetic signal and differences between tribes) on diaspore traits and whether they are related to distribution patterns across the Andes, using phylogenomics, distribution modelling and community ecology analyses. KEY RESULTS: Seventy-five percent of the studied species show small achenes (length <3.3 mm) and 67% have high dispersal investment (pappus length at least two times the achene length). Dispersal investment increases with elevation, possibly to compensate for lower air density, and achene length increases towards the equator, where non-seasonal climate prevails. Diaspore traits show significant phylogenetic signal, and higher dispersal investment is observed in Gnaphalieae, Astereae and Senecioneae, which together represent 72% of our species. High-Andean-restricted species found across the tropical Andes have, on average, the pappus four times longer than the achene, a significantly higher dispersal investment than species present only in the northern Andes or only in the central Andes. CONCLUSIONS: Small achenes and high diaspore dispersal investment dominate among high-Andean Compositae, traits typical of mostly three tribes of African origin; but traits are also correlated with the environmental gradients within the high-Andean grasslands. Our results also suggest that diaspore dispersal investment is likely to shape species distribution patterns in naturally fragmented habitats.

Giant Fern Genomes Show Complex Evolution Patterns: A Comparative Analysis in Two Species of Tmesipteris (Psilotaceae).

Giant genomes are rare across the plant kingdom and their study has focused almost exclusively on angiosperms and gymnosperms. The scarce genetic data that are available for ferns, however, indicate differences in their genome organization and a lower dynamism compared to other plant groups. Tmesipteris is a small genus of mainly epiphytic ferns that occur in Oceania and several Pacific Islands. So far, only two species with giant genomes have been reported in the genus, T. tannensis (1C = 73.19 Gbp) and T. obliqua (1C = 147.29 Gbp). Low-coverage genome skimming sequence data were generated in these two species and analyzed using the RepeatExplorer2 pipeline to identify and quantify the repetitive DNA fraction of these genomes. We found that both species share a similar genomic composition, with high repeat diversity compared to taxa with small (1C < 10 Gbp) genomes. We also found that, in general, characterized repetitive elements have relatively high heterogeneity scores, indicating ancient diverging evolutionary trajectories. Our results suggest that a whole genome multiplication event, accumulation of repetitive elements, and recent activation of those repeats have all played a role in shaping these genomes. It will be informative to compare these data in the future with data from the giant genome of the angiosperm Paris japonica, to determine if the structures observed here are an emergent property of massive genomic inflation or derived from lineage specific processes.

Hybridisation and chloroplast capture between distinct Themeda triandra lineages in Australia.

Ecotypes are distinct populations within a species that are adapted to specific environmental conditions. Understanding how these ecotypes become established, and how they interact when reunited, is fundamental to elucidating how ecological adaptations are maintained. This study focuses on Themeda triandra, a dominant grassland species across Asia, Africa and Australia. It is the most widespread plant in Australia, where it has distinct ecotypes that are usually restricted to either wetter and cooler coastal regions or the drier and hotter interior. We generate a reference genome for T. triandra and use whole genome sequencing for over 80 Themeda accessions to reconstruct the evolutionary history of T. triandra and related taxa. Organelle phylogenies confirm that Australia was colonized by T. triandra twice, with the division between ecotypes predating their arrival in Australia. The nuclear genome provides evidence of differences in the dominant ploidal level and gene-flow among the ecotypes. In northern Queensland there appears to be a hybrid zone between ecotypes with admixed nuclear genomes and shared chloroplast haplotypes. Conversely, in the cracking claypans of Western Australia, there is cytonuclear discordance with individuals possessing the coastal chloroplast and interior clade nuclear genome. This chloroplast capture is potentially a result of adaptive introgression, with selection detected in the rpoC2 gene which is associated with water use efficiency. The reason that T. triandra is the most widespread plant in Australia appears to be a result of distinct ecotypic genetic variation and genome duplication, with the importance of each depending on the geographic scale considered.

Low dispersal and ploidy differences in a grass maintain photosynthetic diversity despite gene flow and habitat overlap.

Geographical isolation facilitates the emergence of distinct phenotypes within a single species, but reproductive barriers or selection are needed to maintain the polymorphism after secondary contact. Here, we explore the processes that maintain intraspecific variation of C4 photosynthesis, a complex trait that results from the combined action of multiple genes. The grass Alloteropsis semialata includes C4 and non-C4 populations, which have coexisted as a polyploid series for more than 1 million years in the miombo woodlands of Africa. Using population genomics, we show that there is genome-wide divergence for the photosynthetic types, but the current geographical distribution does not reflect a simple habitat displacement scenario as the genetic clusters overlap, being occasionally mixed within a given habitat. Despite evidence of recurrent introgression between non-C4 and C4 groups, in both diploids and polyploids, the distinct genetic lineages retain their identity, potentially because of selection against hybrids. Coupled with strong isolation by distance within each genetic group, this selection created a geographical mosaic of photosynthetic types. Diploid C4 and non-C4 types never grew together, and the C4 type from mixed populations constantly belonged to the hexaploid lineage. By limiting reproductive interactions between photosynthetic types, the ploidy difference probably allows their co-occurrence, reinforcing the functional diversity within this species. Together, these factors enabled the persistence of divergent physiological traits of ecological importance within a single species despite gene flow and habitat overlap.

Biogeography and genome size evolution of the oldest extant vascular plant genus, Equisetum (Equisetaceae).

BACKGROUND AND AIMS: Extant plant groups with a long fossil history are key elements in understanding vascular plant evolution. Horsetails (Equisetum, Equisetaceae) have a nearly continuous fossil record dating back to the Carboniferous, but their phylogenetic and biogeographic patterns are still poorly understood. We use here the most extensive phylogenetic analysis to date as a framework to evaluate their age, biogeography and genome size evolution. METHODS: DNA sequences of four plastid loci were used to estimate divergence times and investigate the biogeographic history of all extant species of Equisetum. Flow cytometry was used to study genome size evolution against the framework of phylogenetic relationships in Equisetum. KEY RESULTS: On a well-supported phylogenetic tree including all extant Equisetum species, a molecular clock calibrated with multiple fossils places the node at which the outgroup and Equisetum diverged at 343 Mya (Early Carboniferous), with the first major split among extant species occurring 170 Mya (Middle Jurassic). These dates are older than those reported in some other recent molecular clock studies but are largely in agreement with a timeline established by fossil appearance in the geological record. Representatives of evergreen subgenus Hippochaete have much larger genome sizes than those of deciduous subgenus Equisetum, despite their shared conserved chromosome number. Subgenus Paramochaete has an intermediate genome size and maintains the same number of chromosomes. CONCLUSIONS: The first divergences among extant members of the genus coincided with the break-up of Pangaea and the resulting more humid, warmer climate. Subsequent tectonic activity most likely involved vicariance events that led to species divergences combined with some more recent, long-distance dispersal events. We hypothesize that differences in genome size between subgenera may be related to the number of sperm flagellae.

Contrasted histories of organelle and nuclear genomes underlying physiological diversification in a grass species.

C4 photosynthesis evolved multiple times independently in angiosperms, but most origins are relatively old so that the early events linked to photosynthetic diversification are blurred. The grass Alloteropsis semialata is an exception, as this species encompasses C4 and non-C4 populations. Using phylogenomics and population genomics, we infer the history of dispersal and secondary gene flow before, during and after photosynthetic divergence in A. semialata. We further analyse the genome composition of individuals with varied ploidy levels to establish the origins of polyploids in this species. Detailed organelle phylogenies indicate limited seed dispersal within the mountainous region of origin and the emergence of a C4 lineage after dispersal to warmer areas of lower elevation. Nuclear genome analyses highlight repeated secondary gene flow. In particular, the nuclear genome associated with the C4 phenotype was swept into a distantly related maternal lineage probably via unidirectional pollen flow. Multiple intraspecific allopolyploidy events mediated additional secondary genetic exchanges between photosynthetic types. Overall, our results show that limited dispersal and isolation allowed lineage divergence, with photosynthetic innovation happening after migration to new environments, and pollen-mediated gene flow led to the rapid spread of the derived C4 physiology away from its region of origin.

Polyploidy in gymnosperms - Insights into the genomic and evolutionary consequences of polyploidy in Ephedra.

While polyploidization is recognized as a major evolutionary driver for ferns and angiosperms, little is known about its impact in gymnosperms, where polyploidy is much less frequent. We explore Ephedra to evaluate (i) the extent of genome size diversity in the genus and the influence polyploidy has had on the evolution of nuclear DNA contents, and (ii) identify where shifts in genome size and polyploidy have occurred both temporally and spatially. A phylogenetic framework of all Ephedra species together with genome sizes and karyotypes for 87% and 67% of them respectively, were used to explore ploidy evolution and its global distribution patterns. Polyploidy was shown to be extremely common, with 41 species (83%) being polyploid (up to 8×) or having polyploid cytotypes - the highest frequency and level reported for any gymnosperm. Genome size was also diverse, with values ranging ~5-fold (8.09-38.34 pg/1C) - the largest range for any gymnosperm family - and increasing in proportion to ploidy level (i.e. no genome downsizing). Our findings provide novel data which support the view that gymnosperms have a more conserved mode of genomic evolution compared with angiosperms.

Genome size variation at constant chromosome number is not correlated with repetitive DNA dynamism in Anacyclus (Asteraceae).

BACKGROUND AND AIMS: Changes in the amount of repetitive DNA (dispersed and tandem repeats) are considered the main contributors to genome size variation across plant species in the absence of polyploidy. However, the study of repeatome dynamism in groups showing contrasting genomic features and complex evolutionary histories is needed to determine whether other processes underlying genome size variation may have been overlooked. The main aim here was to elucidate which mechanism best explains genome size evolution in Anacyclus (Asteraceae). METHODS: Using data from Illumina sequencing, we analysed the repetitive DNA in all species of Anacyclus, a genus with a reticulate evolutionary history, which displays significant genome size and karyotype diversity albeit presenting a stable chromosome number. KEY RESULTS: By reconstructing ancestral genome size values, we inferred independent episodes of genome size expansions and contractions during the evolution of the genus. However, analysis of the repeatome revealed a similar DNA repeat composition across species, both qualitative and quantitative. Using comparative methods to study repeatome dynamics in the genus, we found no evidence for repeat activity causing genome size variation among species. CONCLUSIONS: Our results, combined with previous cytogenetic data, suggest that genome size differences in Anacyclus are probably related to chromosome rearrangements involving losses or gains of chromosome fragments, possibly associated with homoploid hybridization. These could represent balanced rearrangements that do not disrupt gene dosage in merged genomes, for example via chromosome segment exchanges.

Satellite DNA in Paphiopedilum subgenus Parvisepalum as revealed by high-throughput sequencing and fluorescent in situ hybridization.

BACKGROUND: Satellite DNA is a rapidly diverging, largely repetitive DNA component of many eukaryotic genomes. Here we analyse the evolutionary dynamics of a satellite DNA repeat in the genomes of a group of Asian subtropical lady slipper orchids (Paphiopedilum subgenus Parvisepalum and representative species in the other subgenera/sections across the genus). A new satellite repeat in Paphiopedilum subgenus Parvisepalum, SatA, was identified and characterized using the RepeatExplorer pipeline in HiSeq Illumina reads from P. armeniacum (2n = 26). Reconstructed monomers were used to design a satellite-specific fluorescent in situ hybridization (FISH) probe. The data were also analysed within a phylogenetic framework built using the internal transcribed spacer (ITS) sequences of 45S nuclear ribosomal DNA. RESULTS: SatA comprises c. 14.5% of the P. armeniacum genome and is specific to subgenus Parvisepalum. It is composed of four primary monomers that range from 230 to 359 bp and contains multiple inverted repeat regions with hairpin-loop motifs. A new karyotype of P. vietnamense (2n = 28) is presented and shows that the chromosome number in subgenus Parvisepalum is not conserved at 2n = 26, as previously reported. The physical locations of SatA sequences were visualised on the chromosomes of all seven Paphiopedilum species of subgenus Parvisepalum (2n = 26-28), together with the 5S and 45S rDNA loci using FISH. The SatA repeats were predominantly localisedin the centromeric, peri-centromeric and sub-telocentric chromosome regions, but the exact distribution pattern was species-specific. CONCLUSIONS: We conclude that the newly discovered, highly abundant and rapidly evolving satellite sequence SatA is specific to Paphiopedilum subgenus Parvisepalum. SatA and rDNA chromosomal distributions are characteristic of species, and comparisons between species reveal that the distribution patterns generate a strong phylogenetic signal. We also conclude that the ancestral chromosome number of subgenus Parvisepalum and indeed of all Paphiopedilum could be either 2n = 26 or 28, if P. vietnamense is sister to all species in the subgenus as suggested by the ITS data.

Evolutionary diversification of CYC/TB1-like TCP homologs and their recruitment for the control of branching and floral morphology in Papaveraceae (basal eudicots).

Angiosperms possess enormous morphological variation in plant architectures and floral forms. Previous studies in Pentapetalae and monocots have demonstrated the involvement of TCP domain CYCLOIDEA/TEOSINTE BRANCHED1-like (CYC/TB1) genes in the control of floral symmetry and shoot branching. However, how TCP/CYC-like (CYL) genes originated, evolved and functionally diversified remain unclear. We conducted a comparative functional study in Ranunculales, the sister lineage to all other eudicots, between Eschscholzia californica and Cysticapnos vesicaria, two species of Papaveraceae with actinomorphic and zygomorphic flowers, respectively. Phylogenetic analysis indicates that CYL genes in Papaveraceae form two paralogous lineages, PapaCYL1 and PapaCYL2. Papaveraceae CYL genes show highly diversified expression patterns as well as functions. Enhanced branching by silencing of EscaCYL1 suggests that the role of CYC/TB1-like genes in branching control is conserved in Papaveraceae. In contrast to the arrest of stamen development in Pentapetalae, PapaCYL genes promote stamen initiation and growth. In addition, we demonstrate that CyveCYLs are involved in perianth development, specifying sepal and petal identity in Cysticapnos by regulating the B-class floral organ identity genes. Our data also suggest the involvement of CyveCYL genes in the regulation of flower symmetry in Cysticapnos. Our work provides evidence of the importance of TCP/CYC-like genes in the promotion of morphological diversity across angiosperms.

Genome evolution of ferns: evidence for relative stasis of genome size across the fern phylogeny.

The genome evolution of ferns has been considered to be relatively static compared with angiosperms. In this study, we analyse genome size data and chromosome numbers in a phylogenetic framework to explore three hypotheses: the correlation of genome size and chromosome number, the origin of modern ferns from ancestors with high chromosome numbers, and the occurrence of several whole-genome duplications during the evolution of ferns. To achieve this, we generated new genome size data, increasing the percentage of fern species with genome sizes estimated to 2.8% of extant diversity, and ensuring a comprehensive phylogenetic coverage including at least three species from each fern order. Genome size was correlated with chromosome number across all ferns despite some substantial variation in both traits. We observed a trend towards conservation of the amount of DNA per chromosome, although Osmundaceae and Psilotaceae have substantially larger chromosomes. Reconstruction of the ancestral genome traits suggested that the earliest ferns were already characterized by possessing high chromosome numbers and that the earliest divergences in ferns were correlated with substantial karyological changes. Evidence for repeated whole-genome duplications was found across the phylogeny. Fern genomes tend to evolve slowly, albeit genome rearrangements occur in some clades.

Genome biogeography reveals the intraspecific spread of adaptive mutations for a complex trait.

Physiological novelties are often studied at macro-evolutionary scales such that their micro-evolutionary origins remain poorly understood. Here, we test the hypothesis that key components of a complex trait can evolve in isolation and later be combined by gene flow. We use C4 photosynthesis as a study system, a derived physiology that increases plant productivity in warm, dry conditions. The grass Alloteropsis semialata includes C4 and non-C4 genotypes, with some populations using laterally acquired C4 -adaptive loci, providing an outstanding system to track the spread of novel adaptive mutations. Using genome data from C4 and non-C4 A. semialata individuals spanning the species' range, we infer and date past migrations of different parts of the genome. Our results show that photosynthetic types initially diverged in isolated populations, where key C4 components were acquired. However, rare but recurrent subsequent gene flow allowed the spread of adaptive loci across genetic pools. Indeed, laterally acquired genes for key C4 functions were rapidly passed between populations with otherwise distinct genomic backgrounds. Thus, our intraspecific study of C4 -related genomic variation indicates that components of adaptive traits can evolve separately and later be combined through secondary gene flow, leading to the assembly and optimization of evolutionary innovations.

Recent updates and developments to plant genome size databases.

Two plant genome size databases have been recently updated and/or extended: the Plant DNA C-values database (http://data.kew.org/cvalues), and GSAD, the Genome Size in Asteraceae database (http://www.asteraceaegenomesize.com). While the first provides information on nuclear DNA contents across land plants and some algal groups, the second is focused on one of the largest and most economically important angiosperm families, Asteraceae. Genome size data have numerous applications: they can be used in comparative studies on genome evolution, or as a tool to appraise the cost of whole-genome sequencing programs. The growing interest in genome size and increasing rate of data accumulation has necessitated the continued update of these databases. Currently, the Plant DNA C-values database (Release 6.0, Dec. 2012) contains data for 8510 species, while GSAD has 1219 species (Release 2.0, June 2013), representing increases of 17 and 51%, respectively, in the number of species with genome size data, compared with previous releases. Here we provide overviews of the most recent releases of each database, and outline new features of GSAD. The latter include (i) a tool to visually compare genome size data between species, (ii) the option to export data and (iii) a webpage containing information about flow cytometry protocols.

Genomic incongruence accompanies the evolution of flower symmetry in Eudicots: a case study in the poppy family (Papaveraceae, Ranunculales).

Reconstructing evolutionary trajectories and transitions that have shaped floral diversity relies heavily on the phylogenetic framework on which traits are modelled. In this study, we focus on the angiosperm order Ranunculales, sister to all other eudicots, to unravel higher-level relationships, especially those tied to evolutionary transitions in flower symmetry within the family Papaveraceae. This family presents an astonishing array of floral diversity, with actinomorphic, disymmetric (two perpendicular symmetry axes), and zygomorphic flowers. We generated nuclear and plastid datasets using the Angiosperms353 universal probe set for target capture sequencing (of 353 single-copy nuclear ortholog genes), together with publicly available transcriptome and plastome data mined from open-access online repositories. We relied on the fossil record of the order Ranunculales to date our phylogenies and to establish a timeline of events. Our phylogenomic workflow shows that nuclear-plastid incongruence accompanies topological uncertainties in Ranunculales. A cocktail of incomplete lineage sorting, post-hybridization introgression, and extinction following rapid speciation most likely explain the observed knots in the topology. These knots coincide with major floral symmetry transitions and thus obscure the order of evolutionary events.

The Plant DNA C-Values Database: A One-Stop Shop for Plant Genome Size Data.

Genome size is a plant character with far-reaching implications, ranging from impacts on the financial and computing feasibility of sequencing and assembling genomes all the way to influencing the very ecology and evolution of species. The increasing recognition of the role of genome size in plant science has led to a rising demand for comprehensive and easily accessible sources of genome size data. The Plant DNA C-values database has established itself as a trusted and widely used central hub for users needing to access available plant genome size data, complemented with related cytogenetic (ploidy level) and karyological (chromosome number) information where available. Since its inception in 2001, the database has undergone six major updates to incorporate newly available genome size information, leading to the most recent release (Release 7.1), which comprises data for 12,273 species across all the major land plant and some algal lineages. Here we describe how to use the database efficiently, making use of its different query and filtering settings.

A Bioinformatic Pipeline to Estimate Ploidy Level from Target Capture Sequence Data Obtained from Herbarium Specimens.

Whole genome duplications (WGD) are frequent in many plant lineages; however, ploidy level variation is unknown in most species. The most widely used methods to estimate ploidy levels in plants are chromosome counts, which require living specimens, and flow cytometry estimates, which necessitate living or relatively recently collected samples. Newly described bioinformatic methods have been developed to estimate ploidy levels using high-throughput sequencing data, and these have been optimized in plants by calculating allelic ratio values from target capture data. This method relies on the maintenance of allelic ratios from the genome to the sequence data. For example, diploid organisms will generate allelic data in a 1:1 proportion, with an increasing number of possible allelic ratio combinations occurring in individuals with higher ploidy levels. In this chapter, we explain step-by-step this bioinformatic approach for the estimation of ploidy level.

Virus-induced gene silencing (VIGS) in Cysticapnos vesicaria, a zygomorphic-flowered Papaveraceae (Ranunculales, basal eudicots).

BACKGROUND AND AIMS: Studies of evolutionary diversification in the basal eudicot family Papaveraceae, such as the transition from actinomorphy to zygomorphy, are hampered by the lack of comparative functional studies. So far, gene silencing methods are only available in the actinomorphic species Eschscholzia californica and Papaver somniferum. This study addresses the amenability of Cysticapnos vesicaria, a derived fumitory with zygomorphic flowers, to virus-induced gene silencing (VIGS), and describes vegetative and reproductive traits in this species. METHODS: VIGS-mediated downregulation of the C. vesicaria PHYTOENE DESATURASE gene (CvPDS) and of the FLORICAULA gene CvFLO was carried out using Agrobacterium tumefaciens transfer of Tobacco rattle virus (TRV)-based vectors. Wild-type and vector-treated plants were characterized using reverse transcription-PCR (RT-PCR), in situ hybridization, and macroscopic and scanning electron microscopic imaging. KEY RESULTS: Cysticapnos vesicaria germinates rapidly, can be grown at high density, has a short life cycle and is self-compatible. Inoculation of C. vesicaria with a CvPDS-VIGS vector resulted in strong photobleaching of green parts and reduction of endogenous CvPDS transcript levels. Gene silencing persisted during inflorescence development until fruit set. Inoculation of plants with CvFLO-VIGS affected floral phyllotaxis, symmetry and floral organ identities. CONCLUSIONS: The high penetrance, severity and stability of pTRV-mediated silencing, including the induction of meristem-related phenotypes, make C. vesicaria a very promising new focus species for evolutionary-developmental (evo-devo) studies in the Papaveraceae. This now enables comparative studies of flower symmetry, inflorescence determinacy and other traits that diversified in the Papaveraceae.

Morphological and Genome-Wide Evidence of Homoploid Hybridisation in Urospermum (Asteraceae).

The genus Urospermum is distributed in the Mediterranean region and Macaronesia, and has been introduced to other extra-Mediterranean regions. Although the two species constituting the genus, U. dalechampii and U. picroides, are frequently found together, hybrids have so far only been reported once, from Morocco. However, we found certain individuals in Catalonia, whose intermediate morphology suggested a potential hybrid origin. In this study, we applied morphological and molecular methods to investigate the origin of those individuals. Intermediate features at phenotype, karyological, cytogenetic, and genomic levels were identified in morphologically intermediate individuals, supporting their homoploid hybrid origin. Chloroplast sequence data suggest that U. dalechampii is the maternal progenitor of the hybrid. Together with the intermediate traits displayed, the lack of fertile seeds suggests that hybrids are probably F1. Future monitoring studies will be, nonetheless, needed to evaluate the extent of hybridisation and its potential impact on the biology of the genus.

The rise of grasslands is linked to atmospheric CO2 decline in the late Palaeogene.

Grasslands are predicted to experience a major biodiversity change by the year 2100. A better understanding of how grasslands have responded to past environmental changes will help predict the outcome of current and future environmental changes. Here, we explore the relationship between past atmospheric CO2 and temperature fluctuations and the shifts in diversification rate of Poaceae (grasses) and Asteraceae (daisies), two exceptionally species-rich grassland families (~11,000 and ~23,000 species, respectively). To this end, we develop a Bayesian approach that simultaneously estimates diversification rates through time from time-calibrated phylogenies and correlations between environmental variables and diversification rates. Additionally, we present a statistical approach that incorporates the information of the distribution of missing species in the phylogeny. We find strong evidence supporting a simultaneous increase in diversification rates for grasses and daisies after the most significant reduction of atmospheric CO2 in the Cenozoic (~34 Mya). The fluctuations of paleo-temperatures, however, appear not to have had a significant relationship with the diversification of these grassland families. Overall, our results shed new light on our understanding of the origin of grasslands in the context of past environmental changes.