While prokaryotic pan-genomes have been shown to contain many more genes than any individual organism, the prevalence and functional significance of differentially present genes in eukaryotes remains poorly understood. Whole-genome de novo assembly and annotation of 54 lines of the grass Brachypodium distachyon yield a pan-genome containing nearly twice the number of genes found in any individual genome. Genes present in all lines are enriched for essential biological functions, while genes present in only some lines are enriched for conditionally beneficial functions (e.g., defense and development), display faster evolutionary rates, lie closer to transposable elements and are less likely to be syntenic with orthologous genes in other grasses. Our data suggest that differentially present genes contribute substantially to phenotypic variation within a eukaryote species, these genes have a major influence in population genetics, and transposable elements play a key role in pan-genome evolution.

Brachypodium distachyon is being widely investigated across the world as a model plant for temperate cereals. This annual plant has three cytotypes (2n = 10, 20, 30) that are still regarded as part of a single species. Here, a multidisciplinary study has been conducted on a representative sampling of the three cytotypes to investigate their evolutionary relationships and origins, and to elucidate if they represent separate species. Statistical analyses of 15 selected phenotypic traits were conducted in individuals from 36 lines or populations. Cytogenetic analyses were performed through flow cytometry, fluorescence in situ hybridization (FISH) with genomic (GISH) and multiple DNA sequences as probes, and comparative chromosome painting (CCP). Phylogenetic analyses were based on two plastid (ndhF, trnLF) and five nuclear (ITS, ETS, CAL, DGAT, GI) genes from different Brachypodium lineages, whose divergence times and evolutionary rates were estimated. The phenotypic analyses detected significant differences between the three cytotypes and demonstrated stability of characters in natural populations. Genome size estimations, GISH, FISH and CCP confirmed that the 2n = 10 and 2n = 20 cytotypes represent two different diploid taxa, whereas the 2n = 30 cytotype represents the allotetraploid derived from them. Phylogenetic analysis demonstrated that the 2n = 20 and 2n = 10 cytotypes emerged from two independent lineages that were, respectively, the maternal and paternal genome donors of the 2n = 30 cytotype. The 2n = 20 lineage was older and mutated significantly faster than the 2n = 10 lineage and all the core perennial Brachypodium species. The substantial phenotypic, cytogenetic and molecular differences detected among the three B. distachyon sensu lato cytotypes are indicative of major speciation processes within this complex that allow their taxonomic separation into three distinct species. We have kept the name B. distachyon for the 2n = 10 cytotype and have described two novel species as B. stacei and B. hybridum for, respectively, the 2n = 20 and 2n = 30 cytotypes.

The genus Brachypodium represents a model system that is advancing our knowledge of the biology of grasses, including small grains, in the postgenomics era. The most widely used species, Brachypodium distachyon, is a C₃ plant that is distributed worldwide. B. distachyon has a small genome, short life cycle, and small stature and is amenable to genetic transformation. Due to the intensive and thoughtful development of this grass as a model organism, it is well-suited for laboratory and field experimentation. The intent of this review is to introduce this model system genus and describe some key outcomes of nearly a decade of research since the first draft genome sequence of the flagship species, B. distachyon, was completed. We discuss characteristics and features of B. distachyon and its congeners that make the genus a valuable model system for studies in ecology, evolution, genetics, and genomics in the grasses, review current hot topics in Brachypodium research, and highlight the potential for future analysis using this system in the coming years.

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.