Brachypodium distachyon (Brachypodium) is rapidly emerging as a powerful model system to facilitate research aimed at improving grass crops for grain, forage and energy production. To characterize the natural diversity of Brachypodium and provide a valuable new tool to the growing list of resources available to Brachypodium researchers, we created and characterized a large, diverse collection of inbred lines. We developed 84 inbred lines from eight locations in Turkey. To enable genotypic characterization of this collection, we created 398 SSR markers from BAC end and EST sequences. An analysis of 187 diploid lines from 56 locations with 43 SSR markers showed considerable genotypic diversity. There was some correlation between SSR genotypes and broad geographic regions, but there was also a high level of genotypic diversity at individual locations. Phenotypic analysis of this new germplasm resource revealed considerable variation in flowering time, seed size, and plant architecture. The inbreeding nature of Brachypodium was confirmed by an extremely high level of homozygosity in wild plants and a lack of cross-pollination under laboratory conditions. Taken together, the inbreeding nature and genotypic diversity observed at individual locations suggest a significant amount of long-distance seed dispersal. The resources developed in this study are freely available to the research community and will facilitate experimental applications based on natural diversity.

The local environment plays a major role in the spatial distribution of plant populations. Natural plant populations have an extremely poor displacing capacity, so their continued survival in a given environment depends on how well they adapt to local pedoclimatic conditions. Genomic tools can be used to identify adaptive traits at a DNA level and to further our understanding of evolutionary processes. Here we report the use of genotyping-by-sequencing on local groups of the sequenced monocot model species Brachypodium distachyon. Exploiting population genetics, landscape genomics and genome wide association studies, we evaluate B. distachyon role as a natural probe for identifying genomic loci involved in environmental adaptation. Brachypodium distachyon individuals were sampled in nine locations with different ecologies and characterized with 16,697 SNPs. Variations in sequencing depth showed consistent patterns at 8,072 genomic bins, which were significantly enriched in transposable elements. We investigated the structuration and diversity of this collection, and exploited climatic data to identify loci with adaptive significance through i) two different approaches for genome wide association analyses considering climatic variation, ii) an outlier loci approach, and iii) a canonical correlation analysis on differentially sequenced bins. A linkage disequilibrium-corrected Bonferroni method was applied to filter associations. The two association methods jointly identified a set of 15 genes significantly related to environmental adaptation. The outlier loci approach revealed that 5.7% of the loci analysed were under selection. The canonical correlation analysis showed that the distribution of some differentially sequenced regions was associated to environmental variation. We show that the multi-faceted approach used here targeted different components of B. distachyon adaptive variation, and may lead to the discovery of genes related to environmental adaptation in natural populations. Its application to a model species with a fully sequenced genome is a modular strategy that enables the stratification of biological material and thus improves our knowledge of the functional loci determining adaptation in near-crop species. When coupled with population genetics and measures of genomic structuration, methods coming from genome wide association studies may lead to the exploitation of model species as natural probes to identify loci related to environmental adaptation.

Brachypodium distachyon (Brachypodium) is a non-domesticated model grass species that can be used to test if variation in genetic sequence or methylation are linked to environmental differences. To assess this, we collected seeds from 12 sites within five climatically distinct regions of Turkey. Seeds from each region were grown under standardized growth conditions in the UK to preserve methylated sequence variation. At six weeks following germination, leaves were sampled and assessed for genomic and DNA methylation variation. In a follow-up experiment, phenomic approaches were used to describe plant growth and drought responses. Genome sequencing and population structure analysis suggested three ancestral clusters across the Mediterranean, two of which were geographically separated in Turkey into coastal and central subpopulations. Phenotypic analyses showed that the coastal subpopulation tended to exhibit relatively delayed flowering and the central, increased drought tolerance as indicated by reduced yellowing. Genome-wide methylation analyses in GpC, CHG and CHH contexts also showed variation which aligned with the separation into coastal and central subpopulations. The climate niche modelling of both subpopulations showed a significant influence from the “Precipitation in the Driest Quarter” on the central subpopulation and “Temperature of the Coldest Month” on the coastal subpopulation. Our work demonstrates genetic diversity and variation in DNA methylation in Turkish accessions of Brachypodium that may be associated with climate variables and the molecular basis of which will feature in ongoing analyses.

The evolution of chromosome number and ribosomal DNA (rDNA) loci number and localisation were studied in Onobrychis Mill. Diploid and tetraploid species, as well as two basic chromosome numbers, x = 7 and x = 8, were observed among analysed taxa. The chromosomal distribution of rDNA loci was presented here for the first time using fluorescence in situ hybridisation (FISH) with 5S and 35S rDNA probes. Onobrychis species showed a high polymorphism in the number and localisation of rDNA loci among diploids, whereas the rDNA loci pattern was very similar in polyploids. Phylogenetic relationships among the species, inferred from nrITS sequences, were used as a framework to reconstruct the patterns of basic chromosome number and rDNA loci evolution. Analysis of the evolution of the basic chromosome numbers allowed the inference of x = 8 as the ancestral number and the descending dysploidy and polyploidisation as the major mechanisms of the chromosome number evolution. Analyses of chromosomal patterns of rRNA gene loci in a phylogenetic context resulted in the reconstruction of one locus of 5S rDNA and one locus of 35S rDNA in the interstitial chromosomal position as the ancestral state in this genus.

Brachypodium distachyon (Brachypodium) is a non-domesticated model grass that has been used to assess population level genomic variation. We have previously established a collection of 55 Brachypodium accessions that were sampled to reflect five different climatic regions of Turkey; designated 1a, 1c, 2, 3 and 4. Genomic and methylomic variation differentiated the collection into two subpopulations designated as coastal and central (respectively from regions 1a, 1c and the other from 2, 3 and 4) which were linked to environmental variables such as relative precipitation. Here, we assessed how far genomic variation would be reflected in the metabolomes and if this could be linked to an adaptive trait. Metabolites were extracted from eight-week-old seedlings from each accession and assessed using flow infusion high-resolution mass spectrometry (FIE-HRMS). Principal Component Analysis (PCA) of the derived metabolomes differentiated between samples from coastal and central subpopulations. The major sources of variation between seedling from the coastal and central subpopulations were identified. The central subpopulation was typified by significant increases in alanine, aspartate and glutamate metabolism and the tricarboxylic acid (TCA) cycle. Coastal subpopulation exhibited elevated levels of the auxin, indolacetic acid and rhamnose. The metabolomes of the seedling were also determined following the imposition of drought stress for seven days. The central subpopulation exhibited a metabolomic shift in response to drought, but no significant changes were seen in the coastal one. The drought responses in the central subpopulation were typified by changes in amino acids, increasing the glutamine that could be functioning as a stress signal. There were also changes in sugars that were likely to be an osmotic counter to drought, and changes in bioenergetic metabolism. These data indicate that genomic variation in our Turkish Brachypodium collection is largely reflected as distinctive metabolomes (“metabolotypes”) through which drought tolerance might be mediated.

Summary: Nucleolar dominance (ND) consists of the reversible silencing of 35S/45S rDNA loci inherited from one of the ancestors of an allopolyploid. The molecular mechanisms by which one ancestral rDNA set is selected for silencing remain unclear. We applied a combination of molecular (Southern blot hybridization and reverse‐transcription cleaved amplified polymorphic sequence analysis), genomic (analysis of variants) and cytogenetic (fluorescence in situ hybridization) approaches to study the structure, expression and epigenetic landscape of 35S rDNA in an allotetraploid grass that exhibits ND, Brachypodium hybridum (genome composition DDSS), and its putative progenitors, Brachypodium distachyon (DD) and Brachypodium stacei (SS). In progenitor genomes, B. stacei showed a higher intragenomic heterogeneity of rDNA compared with B. distachyon. In all studied accessions of B. hybridum, there was a reduction in the copy number of S homoeologues, which was accompanied by their inactive transcriptional status. The involvement of DNA methylation in CG and CHG contexts in the silencing of the S‐genome rDNA loci was revealed. In the B. hybridum allotetraploid, ND is stabilized towards the D‐genome units, irrespective of the polyphyletic origin of the species, and does not seem to be influenced by homoeologous 35S rDNA ratios and developmental stage.