Population Structure in the Model Grass <i>Brachypodium distachyon</i> Is Highly Correlated with Flowering Differences across Broad Geographic Areas

The small, annual grass (L.) Beauv., a close relative of wheat ( L.) and barley ( 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 by the lack of large numbers of well-characterized and sufficiently diverse accessions. Here, we report on genotyping-by-sequencing (GBS) of 84 , seven , and three 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 genome, which appears as a mosaic of -like and -like sequences. Analysis of more than 50,000 SNPs for the 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 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 populations. Our characterization elucidates genes underlying population differences, expands the germplasm resources available for , and illustrates the feasibility and limitations of GWAS in this model grass.