Diversity and association of phenotypic and metabolomic traits in the close model grasses Brachypodium distachyon, B. stacei and B. hybridum

Repeated migration, interbreeding and bottlenecking shaped the phylogeography of the selfing grass Brachypodium stacei

Brachypodium stacei is the most ancestral lineage in the genus Brachypodium, a model system for grass functional genomics. B. stacei shows striking and sometimes contradictory biological and evolutionary features, including a high selfing rate yet extensive admixture, an ancient Miocene origin yet with recent evolutionary radiation, and adaptation to different dry climate conditions in its narrow distribution range. Therefore, it constitutes an ideal system to study these life history traits. We studied the phylogeography of 17 native circum-Mediterranean B. stacei populations (39 individuals) using genome-wide RADseq SNP data and complete plastome sequences. Nuclear SNP data revealed the existence of six distinct genetic clusters, low levels of intra-population genetic diversity and high selfing rates, albeit with signatures of admixture. Coalescence-based dating analysis detected a recent split between crown lineages in the Late Quaternary. Plastome sequences showed incongruent evolutionary relationships with those recovered by the nuclear data, suggesting interbreeding and chloroplast capture events between genetically distant populations. Demographic and population dispersal coalescent models identified an ancestral origin of B. stacei in the western-central Mediterranean islands, followed by an early colonization of the Canary Islands and two independent colonization events of the eastern Mediterranean region through long-distance dispersal and bottleneck events as the most likely evolutionary history. Climate niche data identified three arid niches of B. stacei in the southern Mediterranean region. Our findings indicate that the phylogeography of B. stacei populations was shaped by recent radiations, frequent extinctions, long-distance dispersal events, occasional interbreeding, and adaptation to local climates.

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.

Gradual polyploid genome evolution revealed by pan-genomic analysis of Brachypodium hybridum and its diploid progenitors

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.

Multiple founder events explain the genetic diversity and structure of the model allopolyploid grass Brachypodium hybridum in the Iberian Peninsula hotspot

BACKGROUND AND AIMS

It is accepted that contemporary allopolyploid species have originated recurrently, but very few cases have been documented using multiple natural formations of the same species. To extend our knowledge, we have investigated the multiple origins, genetic variation and structure of the allotetraploid grass Brachypodium hybridum with respect to its progenitor diploid species B. distachyon (D genome) and B. stacei (S genome). For this, our primary focus is the Iberian Peninsula, an evolutionary hotspot for the genus Brachypodium.

METHODS

We analysed 342 B. hybridum individuals from 36 populations using ten nuclear SSR loci and two plastid loci. The B. hybridum genetic profiles were compared with those previously reported for B. stacei and B. distachyon. In addition, phylogenetic analysis of the plastid data was performed for a reduced subset of individuals.

KEY RESULTS

The nuclear simple sequence repeat (SSR) genetic analysis detected medium to high genetic diversity, with a strong south-to-north genetic structure cline, and a high selfing rate in B. hybridum. Comparative genetic analysis showed a close relatedness of current B. hybridum D allelic profiles with those of B. distachyon, but a lack of similarity with those of B. stacei, suggesting another B. stacei source for the B. hybridum S alleles. Plastid analysis detected three different bidirectional allopolyploidization events: two involved distinct B. distachyon-like ancestors and one involved a B. stacei-like ancestor. The south-eastern Iberian Peninsula B. hybridum populations were more genetically diverse and could have originated from at least two hybridization events whereas north-eastern/north-western Iberian Peninsula B. hybridum populations were less diverse and may have derived from at least one hybridization event.

CONCLUSIONS

The genetic and evolutionary evidence supports the plausible in situ origin of the south-eastern and northern Iberian Peninsula B. hybridum allopolyploids from their respective local B. distachyon and unknown B. stacei ancestors. The untapped multiple origins and genetic variation detected in these B. hybridum populations opens the way to future evolutionary analysis of allopolyploid formation and genomic dominance and expression in the B. hybridum-B. distachyon-B. stacei grass model complex.

Sustainable bioenergy for climate mitigation: developing drought-tolerant trees and grasses

BACKGROUND AND AIMS

Bioenergy crops are central to climate mitigation strategies that utilize biogenic carbon, such as BECCS (bioenergy with carbon capture and storage), alongside the use of biomass for heat, power, liquid fuels and, in the future, biorefining to chemicals. Several promising lignocellulosic crops are emerging that have no food role - fast-growing trees and grasses - but are well suited as bioenergy feedstocks, including Populus, Salix, Arundo, Miscanthus, Panicum and Sorghum.

SCOPE

These promising crops remain largely undomesticated and, until recently, have had limited germplasm resources. In order to avoid competition with food crops for land and nature conservation, it is likely that future bioenergy crops will be grown on marginal land that is not needed for food production and is of poor quality and subject to drought stress. Thus, here we define an ideotype for drought tolerance that will enable biomass production to be maintained in the face of moderate drought stress. This includes traits that can readily be measured in wide populations of several hundred unique genotypes for genome-wide association studies, alongside traits that are informative but can only easily be assessed in limited numbers or training populations that may be more suitable for genomic selection. Phenotyping, not genotyping, is now the major bottleneck for progress, since in all lignocellulosic crops studied extensive use has been made of next-generation sequencing such that several thousand markers are now available and populations are emerging that will enable rapid progress for drought-tolerance breeding. The emergence of novel technologies for targeted genotyping by sequencing are particularly welcome. Genome editing has already been demonstrated for Populus and offers significant potential for rapid deployment of drought-tolerant crops through manipulation of ABA receptors, as demonstrated in Arabidopsis, with other gene targets yet to be tested.

CONCLUSIONS

Bioenergy is predicted to be the fastest-developing renewable energy over the coming decade and significant investment over the past decade has been made in developing genomic resources and in collecting wild germplasm from within the natural ranges of several tree and grass crops. Harnessing these resources for climate-resilient crops for the future remains a challenge but one that is likely to be successful.

Untargeted metabolic profiling reveals geography as the strongest predictor of metabolic phenotypes of a cosmopolitan weed

Plants produce a multitude of metabolites that contribute to their fitness and survival and play a role in local adaptation to environmental conditions. The effects of environmental variation are particularly well studied within the genus Plantago; however, previous studies have largely focused on targeting specific metabolites. Studies exploring metabolome-wide changes are lacking, and the effects of natural environmental variation and herbivory on the metabolomes of plants growing in situ remain unknown. An untargeted metabolomic approach using ultra-high-performance liquid chromatography-mass spectrometry, coupled with variation partitioning, general linear mixed modeling, and network analysis was used to detect differences in metabolic phenotypes of Plantago major in fifteen natural populations across Denmark. Geographic region, distance, habitat type, phenological stage, soil parameters, light levels, and leaf area were investigated for their relative contributions to explaining differences in foliar metabolomes. Herbivory effects were further investigated by comparing metabolomes from damaged and undamaged leaves from each plant. Geographic region explained the greatest number of significant metabolic differences. Soil pH had the second largest effect, followed by habitat and leaf area, while phenological stage had no effect. No evidence of the induction of metabolic features was found between leaves damaged by herbivores compared to undamaged leaves on the same plant. Differences in metabolic phenotypes explained by geographic factors are attributed to genotypic variation and/or unmeasured environmental factors that differ at the regional level in Denmark. A small number of specialized features in the metabolome may be involved in facilitating the success of a widespread species such as Plantago major into such wide range of environmental conditions, although overall resilience in the metabolome was found in response to environmental parameters tested. Untargeted metabolomic approaches have great potential to improve our understanding of how specialized plant metabolites respond to environmental change and assist in adaptation to local conditions.

Variation in functional responses to water stress and differentiation between natural allopolyploid populations in the Brachypodium distachyon species complex

Background and Aims

Some polyploid species show enhanced physiological tolerance to drought compared with their progenitors. However, very few studies have examined the consistency of physiological drought response between genetically differentiated natural polyploid populations, which is key to evaluation of the importance of adaptive evolution after polyploidization in those systems where drought exerts a selective pressure.

Methods

A comparative functional approach was used to investigate differentiation of drought-tolerance-related traits in the Brachypodium species complex, a model system for grass polyploid adaptive speciation and functional genomics that comprises three closely related annual species: the two diploid parents, B. distachyon and B. stacei, and the allotetraploid derived from them, B. hybridum. Differentiation of drought-tolerance-related traits between ten genetically distinct B. hybridum populations and its ecological correlates was further analysed.

Key Results

The functional drought response is overall well differentiated between Brachypodium species. Brachypodium hybridum allotetraploids showed a transgressive expression pattern in leaf phytohormone content in response to drought. In contrast, other B. hybridum physiological traits correlated to B. stacei ones. Particularly, proline and water content were the traits that best discriminated these species from B. distachyon under drought.

Conclusions

After polyploid formation and/or colonization, B. hybridum populations have adaptively diverged physiologically and genetically in response to variations in aridity.

Neutral molecular markers support common origin of aluminium tolerance in three congeneric grass species growing in acidic soils

Aluminium (Al) toxicity is the main abiotic stress limiting plant productivity in acidic soils that are widely distributed among arable lands. Plant species differ in the level of Al resistance showing intraspecific and interspecific variation in many crop species. However, the origin of Al-tolerance is not well known. Three annual species, difficult to distinguish phenotypically and that were until recently misinterpreted as a single complex species under Brachypodium distachyon, have been recently separated into three distinct species: the diploids B. distachyon (2n = 10) and B. stacei (2n = 20), and B. hybridum (2n = 30), the allotetraploid derived from the two diploid species. The aims of this work were to know the origin of Al-tolerance in acidic soil conditions within these three Brachypodium species and to develop new DNA markers for species discrimination. Two multiplex SSR-PCRs allowed to genotype a group of 94 accessions for 17 pentanucleotide microsatellite (SSRs) loci. The variability for 139 inter-microsatellite (ISSRs) markers was also examined. The genetic relationships obtained using those neutral molecular markers (SSRs and ISSRs) support that all Al-tolerant allotetraploid accessions of B. hybridum have a common origin that is related with both geographic location and acidic soils. The possibility that the adaptation to acidic soils caused the isolation of the tolerant B. hybridum populations from the others is discussed. We finally describe a new, easy, DNA barcoding method based in the upstream-intron 1 region of the ALMT1 gene, a tool that is 100 % effective to distinguish among these three Brachypodium species.

Environmental isolation explains Iberian genetic diversity in the highly homozygous model grass Brachypodium distachyon

BACKGROUND

Brachypodium distachyon (Poaceae), an annual Mediterranean Aluminum (Al)-sensitive grass, is currently being used as a model species to provide new information on cereals and biofuel crops. The plant has a short life cycle and one of the smallest genomes in the grasses being well suited to experimental manipulation. Its genome has been fully sequenced and several genomic resources are being developed to elucidate key traits and gene functions. A reliable germplasm collection that reflects the natural diversity of this species is therefore needed for all these genomic resources. However, despite being a model plant, we still know very little about its genetic diversity. As a first step to overcome this gap, we used nuclear Simple Sequence Repeats (nSSR) to study the patterns of genetic diversity and population structure of B. distachyon in 14 populations sampled across the Iberian Peninsula (Spain), one of its best known areas.

RESULTS

We found very low levels of genetic diversity, allelic number and heterozygosity in B. distachyon, congruent with a highly selfing system. Our results indicate the existence of at least three genetic clusters providing additional evidence for the existence of a significant genetic structure in the Iberian Peninsula and supporting this geographical area as an important genetic reservoir. Several hotspots of genetic diversity were detected and populations growing on basic soils were significantly more diverse than those growing in acidic soils. A partial Mantel test confirmed a statistically significant Isolation-By-Distance (IBD) among all studied populations, as well as a statistically significant Isolation-By-Environment (IBE) revealing the presence of environmental-driven isolation as one explanation for the genetic patterns found in the Iberian Peninsula.

CONCLUSIONS

The finding of higher genetic diversity in eastern Iberian populations occurring in basic soils suggests that these populations can be better adapted than those occurring in western areas of the Iberian Peninsula where the soils are more acidic and accumulate toxic Al ions. This suggests that the western Iberian acidic soils might prevent the establishment of Al-sensitive B. distachyon populations, potentially causing the existence of more genetically depauperated individuals.