Population structure in the model grass Brachypodium distachyon is highly correlated with flowering differences across broad geographic areas

The small, annual grass Brachypodium distachyon (L.) Beauv., a close relative of wheat (Triticum aestivum L.) and barley (Hordeum vulgare L.), is a powerful model system for cereals and bioenergy grasses. Genome-wide association studies (GWAS) of natural variation can elucidate the genetic basis of complex traits but have been so far limited in B. distachyon by the lack of large numbers of well-characterized and sufficiently diverse accessions. Here, we report on genotyping-by-sequencing (GBS) of 84 B. distachyon, seven B. hybridum, and three B. stacei accessions with diverse geographic origins including Albania, Armenia, Georgia, Italy, Spain, and Turkey. Over 90,000 high-quality single-nucleotide polymorphisms (SNPs) distributed across the Bd21 reference genome were identified. Our results confirm the hybrid nature of the B. hybridum genome, which appears as a mosaic of B. distachyon-like and B. stacei-like sequences. Analysis of more than 50,000 SNPs for the B. distachyon accessions revealed three distinct, genetically defined populations. Surprisingly, these genomic profiles are associated with differences in flowering time rather than with broad geographic origin. High levels of differentiation in loci associated with floral development support the differences in flowering phenology between B. distachyon populations. Genome-wide association studies combining genotypic and phenotypic data also suggest the presence of one or more photoperiodism, circadian clock, and vernalization genes in loci associated with flowering time variation within B. distachyon populations. Our characterization elucidates genes underlying population differences, expands the germplasm resources available for Brachypodium, and illustrates the feasibility and limitations of GWAS in this model grass.