Quantitative trait loci associated with natural diversity in water-use efficiency and response to soil drying in Brachypodium distachyon

Patterns of pan-genome occupancy and gene coexpression under water-deficit in Brachypodium distachyon

Natural populations are characterized by abundant genetic diversity driven by a range of different types of mutation. The tractability of sequencing complete genomes has allowed new insights into the variable composition of genomes, summarized as a species pan-genome. These analyses demonstrate that many genes are absent from the first reference genomes, whose analysis dominated the initial years of the genomic era. Our field now turns towards understanding the functional consequence of these highly variable genomes. Here, we analysed weighted gene coexpression networks from leaf transcriptome data for drought response in the purple false brome Brachypodium distachyon and the differential expression of genes putatively involved in adaptation to this stressor. We specifically asked whether genes with variable "occupancy" in the pan-genome - genes which are either present in all studied genotypes or missing in some genotypes - show different distributions among coexpression modules. Coexpression analysis united genes expressed in drought-stressed plants into nine modules covering 72 hub genes (87 hub isoforms), and genes expressed under controlled water conditions into 13 modules, covering 190 hub genes (251 hub isoforms). We find that low occupancy pan-genes are under-represented among several modules, while other modules are over-enriched for low-occupancy pan-genes. We also provide new insight into the regulation of drought response in B. distachyon, specifically identifying one module with an apparent role in primary metabolism that is strongly responsive to drought. Our work shows the power of integrating pan-genomic analysis with transcriptomic data using factorial experiments to understand the functional genomics of environmental response.

Variation in climatic tolerance, but not stomatal traits, partially explains Pooideae grass species distributions

BACKGROUND AND AIMS

Grasses in subfamily Pooideae live in some of the world's harshest terrestrial environments, from frigid boreal zones to the arid windswept steppe. It is hypothesized that the climate distribution of species within this group is driven by differences in climatic tolerance, and that tolerance can be partially explained by variation in stomatal traits.

METHODS

We determined the aridity index (AI) and minimum temperature of the coldest month (MTCM) for 22 diverse Pooideae accessions and one outgroup, and used comparative methods to assess predicted relationships for climate traits versus fitness traits, stomatal diffusive conductance to water (gw) and speed of stomatal closure following drought and/or cold.

KEY RESULTS

Results demonstrate that AI and MTCM predict variation in survival/regreening following drought/cold, and gw under drought/cold is positively correlated with δ 13C-measured water use efficiency (WUE). However, the relationship between climate traits and fitness under drought/cold was not explained by gw or speed of stomatal closure.

CONCLUSIONS

These findings suggest that Pooideae distributions are at least partly determined by tolerance to aridity and above-freezing cold, but that variation in tolerance is not uniformly explained by variation in stomatal traits.

A genetic link between leaf carbon isotope composition and whole-plant water use efficiency in the C4 grass Setaria

Genetic selection for whole-plant water use efficiency (yield per transpiration; WUE_plant ) in any crop-breeding programme requires high-throughput phenotyping of component traits of WUE_plant such as intrinsic water use efficiency (WUE_i ; CO₂ assimilation rate per stomatal conductance). Measuring WUE_i by gas exchange measurements is laborious and time consuming and may not reflect an integrated WUE_i over the life of the leaf. Alternatively, leaf carbon stable isotope composition (δ¹³ C_leaf ) has been suggested as a potential time-integrated proxy for WUE_i that may provide a tool to screen for WUE_plant . However, a genetic link between δ¹³ C_leaf and WUE_plant in a C₄ species has not been well established. Therefore, to determine if there is a genetic relationship in a C₄ plant between δ¹³ C_leaf and WUE_plant under well watered and water-limited growth conditions, a high-throughput phenotyping facility was used to measure WUE_plant in a recombinant inbred line (RIL) population created between the C₄ grasses Setaria viridis and S. italica. Three quantitative trait loci (QTL) for δ¹³ C_leaf were found and co-localized with transpiration, biomass accumulation, and WUE_plant . Additionally, WUE_plant for each of the δ¹³ C_leaf QTL allele classes was negatively correlated with δ¹³ C_leaf , as would be predicted when WUE_i influences WUE_plant . These results demonstrate that δ¹³ C_leaf is genetically linked to WUE_plant , likely to be through their relationship with WUE_i , and can be used as a high-throughput proxy to screen for WUE_plant in these C₄ species.

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.

Generating Plants with Improved Water Use Efficiency

To improve sustainability of agriculture, high yielding crop varieties with improved water use efficiency (WUE) are needed. Despite the feasibility of assessing WUE using different measurement techniques, breeding for WUE and high yield is a major challenge. Factors influencing the trait under field conditions are complex, including different scenarios of water availability. Plants with C3 photosynthesis are able to moderately increase WUE by restricting transpiration, resulting in higher intrinsic WUE (iWUE) at the leaf level. However, reduced CO2 uptake negatively influences photosynthesis and possibly growth and yield as well. The negative correlation of growth and WUE could be partly disconnected in model plant species with implications for crops. In this paper, we discuss recent insights obtained for Arabidopsis thaliana (L.) and the potential to translate the findings to C3 and C4 crops. Our data on Zea mays (L.) lines subjected to progressive drought show that there is potential for improvements in WUE of the maize line B73 at the whole plant level (WUEplant). However, changes in iWUE of B73 and Arabidopsis reduced the assimilation rate relatively more in maize. The trade-off observed in the C4 crop possibly limits the effectiveness of approaches aimed at improving iWUE but not necessarily efforts to improve WUEplant.

A glycosyl transferase family 43 protein involved in xylan biosynthesis is associated with straw digestibility in Brachypodium distachyon

The recalcitrance of secondary plant cell walls to digestion constrains biomass use for the production of sustainable bioproducts and for animal feed. We screened a population of Brachypodium recombinant inbred lines (RILs) for cell wall digestibility using commercial cellulases and detected a quantitative trait locus (QTL) associated with this trait. Examination of the chromosomal region associated with this QTL revealed a candidate gene that encodes a putative glycosyl transferase family (GT) 43 protein, orthologue of IRX14 in Arabidopsis, and hence predicted to be involved in the biosynthesis of xylan. Arabinoxylans form the major matrix polysaccharides in cell walls of grasses, such as Brachypodium. The parental lines of the RIL population carry alternative nonsynonymous polymorphisms in the BdGT43A gene, which were inherited in the RIL progeny in a manner compatible with a causative role in the variation in straw digestibility. In order to validate the implied role of our candidate gene in affecting straw digestibility, we used RNA interference to lower the expression levels of the BdGT43A gene in Brachypodium. The biomass of the silenced lines showed higher digestibility supporting a causative role of the BdGT43A gene, suggesting that it might form a good target for improving straw digestibility in crops.

Construction of a high-density genetic map for grape using specific length amplified fragment (SLAF) sequencing

Genetic maps are important tools in plant genomics and breeding. We report a large-scale discovery of single nucleotide polymorphisms (SNPs) using the specific length amplified fragment sequencing (SLAF-seq) technique for the construction of high-density genetic maps for two elite wine grape cultivars, ‘Chardonnay’ and ‘Beibinghong’, and their 130 F₁ plants. A total of 372.53 M paired-end reads were obtained after preprocessing. The average sequencing depth was 33.81 for ‘Chardonnay’ (the female parent), 48.20 for ‘Beibinghong’ (the male parent), and 12.66 for the F₁ offspring. We detected 202,349 high-quality SLAFs of which 144,972 were polymorphic; 10,042 SNPs were used to construct a genetic map that spanned 1,969.95 cM, with an average genetic distance of 0.23 cM between adjacent markers. This genetic map contains the largest molecular marker number of the grape maps so far reported. We thus demonstrate that SLAF-seq is a promising strategy for the construction of high-density genetic maps; the map that we report here is a good potential resource for QTL mapping of genes linked to major economic and agronomic traits, map-based cloning, and marker-assisted selection of grape.

Quantitative Trait Loci Associated with Drought Tolerance in Brachypodium distachyon

The temperate wild grass Brachypodium distachyon (Brachypodium) serves as model system for studying turf and forage grasses. Brachypodium collections show diverse responses to drought stress, but little is known about the genetic mechanisms of drought tolerance of this species. The objective of this study was to identify quantitative trait loci (QTLs) associated with drought tolerance traits in Brachypodium. We assessed leaf fresh weight (LFW), leaf dry weight (LDW), leaf water content (LWC), leaf wilting (WT), and chlorophyll fluorescence (Fv/Fm) under well-watered and drought conditions on a recombinant inbred line (RIL) population from two parents (Bd3-1 and Bd1-1) known to differ in their drought adaptation. A linkage map of the RIL population was constructed using 467 single nucleotide polymorphism (SNP) markers obtained from genotyping-by-sequencing. The Bd3-1/Bd1-1 map spanned 1,618 cM and had an average distance of 3.5 cM between adjacent single nucleotide polymorphisms (SNPs). Twenty-six QTLs were identified in chromosome 1, 2, and 3 in two experiments, with 14 of the QTLs under well-watered conditions and 12 QTLs under drought stress. In Experiment 1, a QTL located on chromosome 2 with a peak at 182 cM appeared to simultaneously control WT, LWC, and Fv/Fm under drought stress, accounting for 11-18.7% of the phenotypic variation. Allelic diversity of candidate genes DREB2B, MYB, and SPK, which reside in one multi-QTL region, may play a role in the natural variation in whole plant drought tolerance in Brachypodium. Co-localization of QTLs for multiple drought-related traits suggest that the gene(s) involved are important regulators of drought tolerance in Brachypodium.

Interactive effects of water limitation and elevated temperature on the physiology, development and fitness of diverse accessions of Brachypodium distachyon

An enduring question in plant physiology and evolution is how single genotypes of plants optimize performance in diverse, often highly variable, environments. We grew 35 natural accessions of the grass Brachypodium distachyon in four environments in the glasshouse, contrasting soil water deficit, elevated temperature and their interaction. We modeled treatment, genotype and interactive effects on leaf-level and whole-plant traits, including fecundity. We also assessed the relationship between glasshouse-measured traits and parameters related to climate at the place of origin. We found abundant genetic variation in both constitutive and induced traits related to plant-water relations. Most traits showed strong interaction between temperature and water availability, and we observed genotype-by-environment interaction for several traits. Notably, leaf free proline abundance showed a strong effect of genotype × temperature × water. We found strong associations between phenology, biomass and water use efficiency (WUE) with parameters describing climate of origin. Plants respond to multiple stressors in ways not directly predictable from single stressors, underscoring the complex and trait-specific mechanisms of environmental response. Climate-trait correlations support a role for WUE and phenology in local adaptation to climate in B. distachyon.

Leaf Segmentation and Tracking in Arabidopsis thaliana Combined to an Organ-Scale Plant Model for Genotypic Differentiation

A promising method for characterizing the phenotype of a plant as an interaction between its genotype and its environment is to use refined organ-scale plant growth models that use the observation of architectural traits, such as leaf area, containing a lot of information on the whole history of the functioning of the plant. The Phenoscope, a high-throughput automated platform, allowed the acquisition of zenithal images of Arabidopsis thaliana over twenty one days for 4 different genotypes. A novel image processing algorithm involving both segmentation and tracking of the plant leaves allows to extract areas of the latter. First, all the images in the series are segmented independently using a watershed-based approach. A second step based on ellipsoid-shaped leaves is then applied on the segments found to refine the segmentation. Taking into account all the segments at every time, the whole history of each leaf is reconstructed by choosing recursively through time the most probable segment achieving the best score, computed using some characteristics of the segment such as its orientation, its distance to the plant mass center and its area. These results are compared to manually extracted segments, showing a very good accordance in leaf rank and that they therefore provide low-biased data in large quantity for leaf areas. Such data can therefore be exploited to design an organ-scale plant model adapted from the existing GreenLab model for A. thaliana and subsequently parameterize it. This calibration of the model parameters should pave the way for differentiation between the Arabidopsis genotypes.

Serial sequencing of isolength RAD tags for cost-efficient genome-wide profiling of genetic and epigenetic variations

Isolength restriction site-associated DNA (isoRAD) sequencing is a very simple but powerful approach that was originally developed for genome-wide genotyping at minimal labor and cost, and it has recently extended its applicability to allow quantification of DNA methylation levels. The isoRAD method is distinct from other genotyping-by-sequencing (GBS) methods because of its use of special restriction enzymes to produce isolength tags (32-36 bp), and sequencing of these uniform tags can bring many benefits. However, the relatively short tags produced by the original protocol are mostly suited to single-end (SE) sequencing (36-50 bp), and therefore they cannot efficiently match the gradually increased sequencing capacity of next-generation sequencing (NGS) platforms. To address this issue, we describe an advanced protocol that allows the preparation of five concatenated isoRAD tags for Illumina paired-end (PE) sequencing (100-150 bp). The configuration of the five concatenated tags is highly flexible, and can be defined by users to work with a desired combination of samples and/or restriction enzymes to suit specific research purposes. In comparison with the original protocol, the advanced protocol has an additional digestion and ligation step, and library preparation can be completed in ∼8 h.