Global Diversity of the Brachypodium Species Complex as a Resource for Genome-Wide Association Studies Demonstrated for Agronomic Traits in Response to Climate

Coevolution with hosts underpins speciation in brood-parasitic cuckoos

Coevolution between interacting species is thought to increase biodiversity, but evidence linking microevolutionary processes to macroevolutionary patterns is scarce. We leveraged two decades of behavioral research coupled with historical DNA analysis to reveal that coevolution with hosts underpins speciation in brood-parasitic bronze-cuckoos. At a macroevolutionary scale, we show that highly virulent brood-parasitic taxa have higher speciation rates and are more likely to speciate in sympatry than less-virulent and nonparasitic relatives. We reveal the microevolutionary process underlying speciation: Hosts reject cuckoo nestlings, which selects for mimetic cuckoo nestling morphology. Where cuckoos exploit multiple hosts, selection for mimicry drives genetic and phenotypic divergence corresponding to host preference, even in sympatry. Our work elucidates perhaps the most common, but poorly characterized, evolutionary process driving biological diversification.

Polygenic architecture of flowering time and its relationship with local environments in the grass Brachypodium distachyon

Synchronizing the timing of reproduction with the environment is crucial in the wild. Among the multiple mechanisms, annual plants evolved to sense their environment, the requirement of cold-mediated vernalization is a major process that prevents individuals from flowering during winter. In many annual plants including crops, both a long and short vernalization requirement can be observed within species, resulting in so-called early-(spring) and late-(winter) flowering genotypes. Here, using the grass model Brachypodium distachyon, we explored the link between flowering-time-related traits (vernalization requirement and flowering time), environmental variation, and diversity at flowering-time genes by combining measurements under greenhouse and outdoor conditions. These experiments confirmed that B. distachyon natural accessions display large differences regarding vernalization requirements and ultimately flowering time. We underline significant, albeit quantitative effects of current environmental conditions on flowering-time-related traits. While disentangling the confounding effects of population structure on flowering-time-related traits remains challenging, population genomics analyses indicate that well-characterized flowering-time genes may contribute significantly to flowering-time variation and display signs of polygenic selection. Flowering-time genes, however, do not colocalize with genome-wide association peaks obtained with outdoor measurements, suggesting that additional genetic factors contribute to flowering-time variation in the wild. Altogether, our study fosters our understanding of the polygenic architecture of flowering time in a natural grass system and opens new avenues of research to investigate the gene-by-environment interaction at play for this trait.

Natural variation in Brachypodium distachyon responses to combined abiotic stresses

The demand for agricultural production is becoming more challenging as climate change increases global temperature and the frequency of extreme weather events. This study examines the phenotypic variation of 149 accessions of Brachypodium distachyon under drought, heat, and the combination of stresses. Heat alone causes the largest amounts of tissue damage while the combination of stresses causes the largest decrease in biomass compared to other treatments. Notably, Bd21-0, the reference line for B. distachyon, did not have robust growth under stress conditions, especially the heat and combined drought and heat treatments. The climate of origin was significantly associated with B. distachyon responses to the assessed stress conditions. Additionally, a GWAS found loci associated with changes in plant height and the amount of damaged tissue under stress. Some of these SNPs were closely located to genes known to be involved in responses to abiotic stresses and point to potential causative loci in plant stress response. However, SNPs found to be significantly associated with a response to heat or drought individually are not also significantly associated with the combination of stresses. This, with the phenotypic data, suggests that the effects of these abiotic stresses are not simply additive, and the responses to the combined stresses differ from drought and heat alone.

High-throughput, dynamic, multi-dimensional: an expanding repertoire of plant respiration measurements

A recent burst of technological innovation and adaptation has greatly improved our ability to capture respiration rate data from plant sources. At the tissue level, several independent respiration measurement options are now available, each with distinct advantages and suitability, including high-throughput sampling capacity. These advancements facilitate the inclusion of respiration rate data into large-scale biological studies such as genetic screens, ecological surveys, crop breeding trials, and multi-omics molecular studies. As a result, our understanding of the correlations of respiration with other biological and biochemical measurements is rapidly increasing. Difficult questions persist concerning the interpretation and utilization of respiration data; concepts such as allocation of respiration to growth versus maintenance, the unnecessary or inefficient use of carbon and energy by respiration, and predictions of future respiration rates in response to environmental change are all insufficiently grounded in empirical data. However, we emphasize that new experimental designs involving novel combinations of respiration rate data with other measurements will flesh-out our current theories of respiration. Furthermore, dynamic recordings of respiration rate, which have long been used at the scale of mitochondria, are increasingly being used at larger scales of size and time to reflect processes of cellular signal transduction and physiological response to the environment. We also highlight how respiratory methods are being better adapted to different plant tissues including roots and seeds, which have been somewhat neglected historically.

Brachypodium: 20 years as a grass biology model system; the way forward?

It has been 20 years since Brachypodium distachyon was suggested as a model grass species, but ongoing research now encompasses the entire genus. Extensive Brachypodium genome sequencing programmes have provided resources to explore the determinants and drivers of population diversity. This has been accompanied by cytomolecular studies to make Brachypodium a platform to investigate speciation, polyploidisation, perenniality, and various aspects of chromosome and interphase nucleus organisation. The value of Brachypodium as a functional genomic platform has been underscored by the identification of key genes for development, biotic and abiotic stress, and cell wall structure and function. While Brachypodium is relevant to the biofuel industry, its impact goes far beyond that as an intriguing model to study climate change and combinatorial stress.

Enriched root bacterial microbiome in invaded vs native ranges of the model grass allotetraploid Brachypodium hybridum

Invasive species can shift the composition of key soil microbial groups, thus creating novel soil microbial communities. To better understand the biological drivers of invasion, we studied plant-microbial interactions in species of the Brachypodium distachyon complex, a model system for functional genomic studies of temperate grasses and bioenergy crops. While Brachypodium hybridum invasion in California is in an incipient stage, threatening natural and agricultural systems, its diploid progenitor species B. distachyon is not invasive in California. We investigated the root, soil, and rhizosphere bacterial composition of Brachypodium hybridum in both its native and invaded range, and of B. distachyon in the native range. We used high-throughput, amplicon sequencing to evaluate if the bacteria associated with these plants differ, and whether biotic controls may be driving B. hybridum invasion. Bacterial community composition of B. hybridum differed based on provenance (native or invaded range) for root, rhizosphere, and bulk soils, as did the abundance of dominant bacterial taxa. Bacteroidetes, Cyanobacteria and Bacillus spp. (species) were significantly more abundant in B. hybridum roots from the invaded range, whereas Proteobacteria, Firmicutes, Erwinia and Pseudomonas were more abundant in the native range roots. Brachypodium hybridum forms novel biotic interactions with a diverse suite of rhizosphere microbes from the invaded range, which may not exert a similar influence within its native range, ostensibly contributing to B. hybridum’s invasiveness. These associated plant microbiomes could inform future management approaches for B. hybridum in its invaded range and could be key to understanding, predicting, and preventing future plant invasions.

Pan-Genome miRNomics in Brachypodium

Pan-genomes are efficient tools for the identification of conserved and varying genomic sequences within lineages of a species. Investigating genetic variations might lead to the discovery of genes present in a subset of lineages, which might contribute into beneficial agronomic traits such as stress resistance or yield. The content of varying genomic regions in the pan-genome could include protein-coding genes as well as microRNA(miRNAs), small non-coding RNAs playing key roles in the regulation of gene expression. In this study, we performed in silico miRNA identification from the genomic sequences of 54 lineages of Brachypodium distachyon, aiming to explore varying miRNA contents and their functional interactions. A total of 115 miRNA families were identified in 54 lineages, 56 of which were found to be present in all lineages. The miRNA families were classified based on their conservation among lineages and potential mRNA targets were identified. Obtaining information about regulatory mechanisms stemming from these miRNAs offers strong potential to provide a better insight into the complex traits that were potentially present in some lineages. Future work could lead us to introduce these traits to different lineages or other economically important plant species in order to promote their survival in different environmental conditions.

Metabolomic Variation Aligns with Two Geographically Distinct Subpopulations of Brachypodium Distachyon before and after Drought Stress

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.

Genetic and Methylome Variation in Turkish Brachypodium Distachyon Accessions Differentiate Two Geographically Distinct Subpopulations

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.

Extending the Genotype in Brachypodium by Including DNA Methylation Reveals a Joint Contribution with Genetics on Adaptive Traits

Epigenomic changes have been considered a potential missing link underlying phenotypic variation in quantitative traits but is potentially confounded with the underlying DNA sequence variation. Although the concept of epigenetic inheritance has been discussed in depth, there have been few studies attempting to directly dissect the amount of epigenomic variation within inbred natural populations while also accounting for genetic diversity. By using known genetic relationships between Brachypodium lines, multiple sets of nearly identical accession families were selected for phenotypic studies and DNA methylome profiling to investigate the dual role of (epi)genetics under simulated natural seasonal climate conditions. Despite reduced genetic diversity, appreciable phenotypic variation was still observable in the measured traits (height, leaf width and length, tiller count, flowering time, ear count) between as well as within the inbred accessions. However, with reduced genetic diversity there was diminished variation in DNA methylation within families. Mixed-effects linear modeling revealed large genetic differences between families and a minor contribution of DNA methylation variation on phenotypic variation in select traits. Taken together, this analysis suggests a limited but significant contribution of DNA methylation toward heritable phenotypic variation relative to genetic differences.

Mechanical stimulation in Brachypodium distachyon: Implications for fitness, productivity, and cell wall properties

Mechanical stimulation, including exposure to wind, is a common environmental variable for plants. However, knowledge about the morphogenetic response of the grasses (Poaceae) to mechanical stimulation and impact on relevant agronomic traits is very limited. Two natural accessions of Brachypodium distachyon were exposed to wind and mechanical treatments. We surveyed a wide range of stem-related traits to determine the effect of the two treatments on plant growth, development, and stem biomass properties. Both treatments induced significant quantitative changes across multiple scales, from the whole plant down to cellular level. The two treatments resulted in shorter stems, reduced biomass, increased tissue rigidity, delayed flowering, and reduced seed yield in both accessions. Among changes in cell wall-related features, a substantial increase in lignin content and pectin methylesterase activity was most notable. Mechanical stimulation also reduced the enzymatic sugar release from the cell wall, thus increasing biomass recalcitrance. Notably, treatments had a distinct and opposite effect on vascular bundle area in the two accessions, suggesting genetic variation in modulating these responses to mechanical stimulation. Our findings highlight that exposure of grasses to mechanical stimulation is a relevant environmental factor affecting multiple traits important for their utilization in food, feed, and bioenergy applications.

Genetic Variation for Seed Metabolite Levels in Brachypodium distachyon

Metabolite composition and concentrations in seed grains are important traits of cereals. To identify the variation in the seed metabolotypes of a model grass, namely Brachypodium distachyon, we applied a widely targeted metabolome analysis to forty inbred lines of B. distachyon and examined the accumulation patterns of 183 compounds in the seeds. By comparing the metabolotypes with the population structure of these lines, we found signature metabolites that represent different accumulation patterns for each of the three B. distachyon subpopulations. Moreover, we found that thirty-seven metabolites exhibited significant differences in their accumulation between the lines Bd21 and Bd3-1. Using a recombinant inbred line (RIL) population from a cross between Bd3-1 and Bd21, we identified the quantitative trait loci (QTLs) linked with this variation in the accumulation of thirteen metabolites. Our metabolite QTL analysis illustrated that different genetic factors may presumably regulate the accumulation of 4-pyridoxate and pyridoxamine in vitamin B6 metabolism. Moreover, we found two QTLs on chromosomes 1 and 4 that affect the accumulation of an anthocyanin, chrysanthemin. These QTLs genetically interacted to regulate the accumulation of this compound. This study demonstrates the potential for metabolite QTL mapping in B. distachyon and provides new insights into the genetic dissection of metabolomic traits in temperate grasses.