Successful germination and seedling development are crucial steps in the growth of a new plant. In this study, we investigated the course of the cell cycle during germination in relation to grain hydration in the model grass Brachypodium distachyon (Brachypodium) for the first time. Flow cytometry was performed to monitor the cell cycle progression during germination and to estimate DNA content in embryo tissues. The analyses of whole zygotic embryos revealed that the relative DNA content was 2C, 4C, 8C, and 16C. Endoreplicated nuclei were detected in the scutellum and coleorhiza cells, whereas the rest of the embryo tissues only had nuclei with a 2C and 4C DNA content. This study was accompanied by a spatiotemporal profile analysis of the DNA synthetic activity in the organs of Brachypodium embryos during germination using EdU labelling. Upon imbibition, nuclear DNA replication was initiated in the radicle within 11 h and subsequently spread towards the plumule. The first EdU-labelled prophases were observed after 14 h of imbibition. Analysis of selected genes that are involved in the regulation of the cell cycle, such as those encoding cyclin-dependent kinases and cyclins, demonstrated an increase in their expression profiles.

Our understanding of polyploid genome evolution is constrained because we cannot know the exact founders of a particular polyploid. To differentiate between founder effects and post polyploidization evolution, we use a pan-genomic approach to study the allotetraploid Brachypodium hybridum and its diploid progenitors. Comparative analysis suggests that most B. hybridum whole gene presence/absence variation is part of the standing variation in its diploid progenitors. Analysis of nuclear single nucleotide variants, plastomes and k-mers associated with retrotransposons reveals two independent origins for B. hybridum, ~1.4 and ~0.14 million years ago. Examination of gene expression in the younger B. hybridum lineage reveals no bias in overall subgenome expression. Our results are consistent with a gradual accumulation of genomic changes after polyploidization and a lack of subgenome expression dominance. Significantly, if we did not use a pan-genomic approach, we would grossly overestimate the number of genomic changes attributable to post polyploidization evolution.

SUMMARY: Unraveling the evolution of plant polyploids is a challenge when their diploid progenitor species are extinct or unknown or when genome sequences of known progenitors are unavailable. Existing subgenome identification methods cannot adequately infer the homeologous genomes that are present in the allopolyploids if they do not take into account the potential existence of unknown progenitors. We addressed this challenge in the widely distributed dysploid grass genus Brachypodium, which is a model genus for temperate cereals and biofuel grasses. We used a transcriptome‐based phylogeny and newly designed subgenome detection algorithms coupled with a comparative chromosome barcoding analysis. Our phylogenomic subgenome detection pipeline was validated in Triticum allopolyploids, which have known progenitor genomes, and then used to infer the identities of three subgenomes derived from extant diploid species and four subgenomes derived from unknown diploid progenitors (ghost subgenomes) in six Brachypodium polyploids (B. mexicanum, B. boissieri, B. retusum, B. phoenicoides, B. rupestre and B. hybridum), of which five contain undescribed homeologous subgenomes. The existence of the seven Brachypodium progenitor genomes in the polyploids was confirmed by their karyotypic barcode profiles. Comparative phylogenomics of nuclear versus plastid trees allowed us to formulate hypothetical homoploid hybridizations and allo‐ and autopolyploidization scenarios that could have generated the six Brachypodium polyploids.

The “genomic shock” hypothesis posits that unusual challenges to genome integrity such as whole genome duplication may induce chaotic genome restructuring. Decades of research on polyploid genomes have revealed that this is often, but not always the case. While some polyploids show major chromosomal rearrangements and derepression of transposable elements in the immediate aftermath of whole genome duplication, others do not. Nonetheless, all polyploids show gradual diploidization over evolutionary time. To evaluate these hypotheses, we produced a chromosome-scale reference genome for the natural allotetraploid grass Brachypodium hybridum, accession “Bhyb26.” We compared 2 independently derived accessions of B. hybridum and their deeply diverged diploid progenitor species Brachypodium stacei and Brachypodium distachyon. The 2 B. hybridum lineages provide a natural timecourse in genome evolution because one formed 1.4 million years ago, and the other formed 140 thousand years ago. The genome of the older lineage reveals signs of gradual post-whole genome duplication genome evolution including minor gene loss and genome rearrangement that are missing from the younger lineage. In neither B. hybridum lineage do we find signs of homeologous recombination or pronounced transposable element activation, though we find evidence supporting steady post-whole genome duplication transposable element activity in the older lineage. Gene loss in the older lineage was slightly biased toward 1 subgenome, but genome dominance was not observed at the transcriptomic level. We propose that relaxed selection, rather than an abrupt genomic shock, drives evolutionary novelty in B. hybridum, and that the progenitor species’ similarity in transposable element load may account for the subtlety of the observed genome dominance.

The Brachypodium genus represents a useful model system to study grass genome organization. Palaeogenomic analyses (e.g. Murat F, Armero A, Pont C, Klopp C, Salse J. 2017. Reconstructing the genome of the most recent common ancestor of flowering plants. Nature Genetics49: 490–496) have identified polyploidization and dysploidy as the prime mechanisms driving the diversity of plant karyotypes and nested chromosome fusions (NCFs) crucial for shaping grass chromosomes. This study compares the karyotype structure and evolution in B. distachyon (genome Bd), B. stacei (genome Bs) and in their putative allotetraploid B. hybridum (genomes BdBs). Brachypodium chromosomes were measured and identified using multicolour fluorescence in situ hybridization (mcFISH). For higher resolution, comparative chromosome barcoding was developed using sets of low-repeat, physically mapped B. distachyon-derived bacterial artificial chromosome (BAC) clones. All species had rather small chromosomes, and essentially all in the Bs genome were morphometrically indistinguishable. Seven BACs combined with two rDNA-based probes provided unambiguous and reproducible chromosome discrimination. Comparative chromosome barcoding revealed NCFs that contributed to the reduction in the x = 12 chromosome number that has been suggested for the intermediate ancestral grass karyotype. Chromosome Bd3 derives from two NCFs of three ancestral chromosomes (Os2, Os8, Os10). Chromosome Bs6 shows an ancient Os8/Os10 NCF, whilst Bs4 represents Os2 only. Chromosome Bd4 originated from a descending dysploidy that involves two NCFs of Os12, Os9 and Os11. The specific distribution of BACs along Bs9 and Bs5, in both B. stacei and B. hybridum, suggests a Bs genome-specific Robertsonian rearrangement. mcFISH-based karyotyping identifies all chromosomes in Brachypodium annuals. Comparative chromosome barcoding reveals rearrangements responsible for the diverse organization of Bd and Bs genomes and provides new data regarding karyotype evolution since the split of the two diploids. The fact that no chromosome rearrangements were observed in B. hybridum compared with the karyotypes of its phylogenetic ancestors suggests prolonged genome stasis after the formation of the allotetraploid.

The Brachypodium genus is an informative model system for studying grass karyotype organization. Previous studies of a limited number of species and reference chromosomes have not provided a comprehensive picture of the enigmatic phylogenetic relationships in the genus. Comparative chromosome barcoding, which enables the reconstruction of the evolutionary history of individual chromosomes and their segments, allowed us to infer the relationships between putative ancestral karyotypes of extinct species and extant karyotypes of current species. We used over 80 chromosome-specific BAC (bacterial artificial chromosome) clones derived from five reference chromosomes of B. distachyon as probes against the karyotypes of twelve accessions representing five diploid and polyploid Brachypodium perennials. The results showed that descending dysploidy is common in Brachypodium and occurs primarily via nested chromosome fusions. Brachypodium distachyon was rejected as a putative ancestor for allotetraploid perennials and B. stacei for B. mexicanum. We propose two alternative models of perennial polyploid evolution involving either the incorporation of a putative x = 5 ancestral karyotype with different descending dysploidy patterns compared to B. distachyon chromosomes or hybridization of two x = 9 ancestors followed by genome doubling and descending dysploidy. Details of the karyotype structure and evolution in several Brachypodium perennials are revealed for the first time.