# Domestication-driven changes in plant traits associated with changes in the assembly of the rhizosphere microbiota in tetraploid wheat

Aymé Spor, Agathe Roucou, Arnaud Mounier, David Bru, Marie-Christine Breuil, Florian Fort, Denis Vile, Pierre Roumet, Laurent Philippot, Cyrille Violle

<figure><img src="image_1.png" /></figure>

## To cite this version:

Aymé Spor, Agathe Roucou, Arnaud Mounier, David Bru, Marie-Christine Breuil, et al. Domestication-driven changes in plant traits associated with changes in the assembly of the rhizosphere microbiota in tetraploid wheat. Scientific Reports, 2020, 10 (1), 10.1038/s41598-020-69175-9. hal-02918224

HAL Id: hal-02918224

https://hal.inrae.fr/hal-02918224

Submitted on 20 Aug 2020

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

OPEN

# Domestication-driven changes in plant traits associated with changes in the assembly of the rhizosphere microbiota in tetraploid wheat

Aymé Spor<fnref n="1" />,⁷✉, Agathe Roucou²,³,⁷, Arnaud Mounier⁴, David Bru<fnref n="1" />, Marie-Christine Breuil<fnref n="1" />, Florian Fort⁵, Denis Vile³, Pierre Roumet⁶, Laurent Philippot<fnref n="1" />✉, & Cyrille Violle²

Despite the large morphological and physiological changes that plants have undergone through domestication, little is known about their impact on their microbiome. Here we characterized rhizospheric bacterial and fungal communities as well as the abundance of N-cycling microbial guilds across thirty-nine accessions of tetraploid wheat, *Triticum turgidum*, from four domestication groups ranging from the wild subspecies to the semi dwarf elite cultivars. We identified several microbial phylotypes displaying significant variation in their relative abundance depending on the wheat domestication group with a stronger impact of domestication on fungi. The relative abundance of potential fungal plant pathogens belonging to the Sordariomycetes class decreased in domesticated compared to wild emmer while the opposite was found for members of the Glomeromycetes, which are obligate plant symbionts. The depletion of nitrifiers and of arbuscular mycorrhizal fungi in elite wheat cultivars compared to primitive domesticated forms suggests that the Green Revolution has decreased the coupling between plant and rhizosphere microbes that are potentially important for plant nutrient availability. Both plant diameter and fine root percentage exhibited the highest number of associations with microbial taxa, highlighting their putative role in shaping the rhizosphere microbiota during domestication. Aside from domestication, significant variation of bacterial and fungal community composition was found among accessions within each domestication group. In particular, the relative abundances of Ophiostomataceae and of Rhizobiales were strongly dependent on the host accession, with heritability estimates of ~ 27% and ~ 25%, indicating that there might be room for genetic improvement via introgression of ancestral plant rhizosphere-beneficial microbe associations.

Crop domestication events led to spectacular modifications of plant traits increasing their suitability to human requirements. Cereals and especially wheat have been domesticated in the Middle East<fnref n="1" />. Domestication process led to dramatic phenotypic changes in cultivated species in relation to cultivation conditions and human needs. The Domestication syndrome² refers to the whole set of phenotypic changes occurred during this process, including: loss of dormancy, increasing seed size, modifying seed dispersal mode and apical dominance as well photoperiod sensitivity.

<footnote n="1">¹Université Bourgogne Franche-Comté, INRAE, AgroSup Dijon, Agroécologie, 21000 Dijon, France. ²CEFE, Univ Montpellier, CNRS, EPHE, IRD, Université Paul Valéry, Montpellier, France. ³LEPSE, INRAE, Institut Agro, Univ Montpellier, Montpellier, France. ⁴Agroécologie, AgroSup Dijon, CNRS, INRAE, Univ. Bourgogne, Univ. Bourgogne Franche-Comté, 21000 Dijon, France. ⁵CEFE, Univ Montpellier, CNRS, EPHE, Institut Agro, IRD, Université Paul Valéry, Montpellier, France. ⁶AGAP, Univ Montpellier, CIRAD, INRAE, Institut Agro, Montpellier, France. ⁷These authors contributed equally: Aymé Spor and Agathe Roucou. ✉email: ayme.spor@inrae.fr; laurent.philippot@inrae.fr</footnote>
Domestication of tetraploid wheat (wild emmer, *Triticum turgidum* ssp. *Dicoccoides* (Körn.), *Thell.*) firstly led to a non-brittle hulled subspecies (cultivated emmer, *T. t.* ssp. *Dicoccum Schrank*) which was one of the first cereals domesticated¹. Then, human selection of mutations at multiple loci, notably at *Tg*—tenacious glume—and *Q* loci³⁻⁵, led to the evolution of traits such as soft glumes, non-hulled grains improving threshing efficiency⁶ and facilitating its cultivation in various environments. This phenotypic evolution together with hybridization between different forms⁷ led to free-threshing subspecies (*T. t.* ssp. *polonicum*, *T. t.* ssp. *carthicum*, *T. t.* ssp. *turgidum*, *T. t.* ssp. *Durum* (*Desf.*) *MacKey*) that were widely spread out with the early agricultural movements. Nowadays, elite durum varieties and landraces from *T. t. durum* grown in different environments are of major importance for grain production in the Mediterranean basin⁸. Although domestication is a continuous process including modern revolution and modern plant breeding, Gioia et al.⁹ suggested to distinguish a primary domestication period (from wild emmer to cultivated emmer) from a secondary one's (from cultivated emmer to the free-threshing subspecies including durum wheat). Molecular markers-based studies showed that *durum* wheat history was associated with a significant decrease of the level of genetic diversity, from wild ancestors to the most recent modern varieties¹⁰,¹¹. It has also been speculated that domestication-driven phenotypic changes have profoundly altered above and belowground biotic interactions, notably because they induce changes in the associated microbiota habitat¹². This is also strengthened by the fact that domesticated crops during the Green Revolution were selected for agroecosystems characterized by an extensive use of pesticides and fertilizers, and are therefore less prone to rely on their associated microbiota for their growth.

Rhizosphere microbes are important for plants because they provide benefits to their host. Both bacteria and fungi, in particular symbiotic rhizobia and arbuscular mycorrhizal fungi (AMF), are known to mediate nutrient acquisition such as nitrogen, iron or phosphorus (reviewed in Ref. ¹³), enhance plant tolerance to environmental stresses¹⁴ and protect against pathogens by acting as the first line of defence against soil-borne pathogens¹⁵. In turn, roots secrete large amounts of rhizodeposits that can be utilized as carbon sources by microbes for growth and development, but also secondary metabolites that can be important for establishing symbiosis with microbes or acting against plant pathogens invasion. Recent literature describing the effect of plant genotypes on their associated rhizosphere microbiota revealed in some species that an heritable intraspecific variation driven by the host's genetics exists¹⁶⁻²⁰. However, little is known about the extent of the impact of plant domestication on the rhizosphere microbiota. When analysing a handful of genotypes of wild and domesticated barley or *Phaseolus vulgaris*, Bulgarelli et al.²¹ and Perez-Jaramillo et al.²² reported small but significant differences in rhizosphere bacterial composition associated with differences in host genotypes. On the opposite, in a more comprehensive study based on 33 sunflower genotypes that varied in their extent of domestication, Leff et al. demonstrated that domestication affected only rhizosphere and seed associated fungal communities and neither root nor rhizosphere bacterial communities²³.

Here, we took advantage of a unique collection of genotypes belonging to three main sub-species representative of tetraploid wheat domestication. These genotypes have been characterized for a variety of above and belowground morphological traits, and have been shown to display contrasted ecological strategies depending on domestication²⁴. We characterized the rhizosphere microbiota of 39 selected accessions, homogeneously divided in four representative accession subsets corresponding to four main steps of the evolutionary history of *Triticum turgidum* i.e. wild, first domesticated, free threshing landraces and post Green Revolution elite lines respectively. We hypothesized that wheat domestication, by modifying key crop traits, might have weakened the association between wheat and its rhizosphere microbiota, and therefore led to significant changes in rhizosphere bacterial and fungal communities.

## Materials and methods

**Plant material.** Forty accessions of tetraploid wheat were selected from the four domestication series groups represented by the three subspecies *T. turgidum* ssp. *dicoccoides* (Körn) *Thell*, wild emmer), *T. turgidum* ssp. *dicoccum*, (cultivated emmer), and *T. turgidum* ssp. *Durum* (*Desf*) *MacKey*. This last sub-species was split in two subgroups to distinguish pre *versus* post Green Revolution durum wheat accessions; the first subgroup is issued from landraces and the second one related to semi dwarf elite accessions registered in Europe after the Green Revolution (from 1970s to 1990s) (Supplementary Table S1).

For each group, ten accessions have been sampled to maximize the genetic diversity of each group based on 21 non-linked microsatellite markers mapped on 14 *durum* chromosomes²⁵. The accessions came from seed collections maintained in different stock centers (see Supplementary Table S1). For each accession, seeds came from successive self-fertilizations realized in common gardens to limit residual heterozygosity and to ensure that the material is genetically fixed. Seeds were provided by INRAE UMR AGAP in Montpellier.

**Experimental design.** We set up the experiment at the INRAE-SupAgro campus (France, 43° 37' 02" N, 3° 51' 18" E) for a duration of 90 days in outdoor conditions (minimal temperature = 2 °C and maximal temperature = 19 °C). As described in Ref.²⁴, plants were grown in 9 L pots (17.5 cm diameter and 38.5 cm depth) filled with 2 mm sieved topsoil taken from a natural grassland (France, 44° 26' 25" N, 3° 55' 58" E; 0.994 g N kg⁻¹ soil, 4.17 g K kg⁻¹ soil and 0.166 g P kg⁻¹ soil). To synchronize germination, we put seeds into Petri dishes for 3 days at 4 °C, and 4 days at ambient temperature. After germination, two seedlings were transplanted (at 2 cm depth) in each pot and thinned to one plant per pot at the first leaf stage. We performed a randomized complete block design using five blocks, where each accession was represented by one plant in each block. To avoid water stress, plants were watered once or three times per week with 150 mL of tap water, according to climatic conditions. One accession from *T. turgidum* ssp. *dicoccoides* did not grow up during the experiment. We therefore removed it from the following analyses (Supplementary Table S1).
At vegetative stage, we determined traits on each plant. Aboveground traits were measured at the leaf level: leaf nitrogen content (LNC, %) and leaf longevity (LL, °Cd), and at the plant level: plant diameter (PD, cm) and shoot dry mass (SDM, g). At harvest, belowground parts were extracted from soil; architectural trait, crow root angle (RA, °), and morphological traits, percent of fine roots (%), mean root diameter (RD, mm) and root dry mass (RDM, g) were determined after rhizosphere soil samples collection (see below) according to the detailed acquisition protocol reported in Ref.²⁴.

**Rhizosphere soil samples collection.** Plants were manually taken off the pots and shaken vigorously to remove loose soil. Roots were then cut and washed twice with 100 mL KCL 1 N in 500 mL sterile centrifuge pots shaken at 350 rpm. The resulting slurry was then centrifuged at 4,000g. The supernatant was removed and used for quantification of nitrogen pools. The same procedure was used for bulk soil samples. Bulk soil samples are samples coming from pots without any plants (five replicates) and from planted pots from eight genotypes (two per subgroup) replicated five times. We therefore have 195 rhizosphere samples and 45 bulk soil samples. The soil pellet was used for subsequent DNA extractions.

**Amplicon generation and MiSeq sequencing.** Amplicons were generated in two steps according to Ref.²⁶. In the first step, the bacterial 16S rRNA gene V3–V4 hypervariable region was amplified by polymerase chain reaction (PCR) using the fusion primers U341F (5'-CCTACGGGRSGCAGCAG-3') and 805R (5'-GAC TACCAGGGTATCTAAT-3')²⁷, with overhang adapters (forward: TCGTCGGCAGCGTCAGATGTGTATAAG AGACAG, adapter: GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG) to allow the subsequent addition of multiplexing index-sequences. PCR was carried out in duplicate 15 µL reactions containing 7.5 µL Phusion High-Fidelity PCR Master Mix (THERMO FISHER SCIENTIFIC), 0.25 µM of each primer, 250 ng T4 gp32 (MPBIO) and 1 ng template DNA. Thermal cycling conditions were 98 °C for 3 min followed by 25 cycles of 98 °C for 30 s, 55 °C for 30 s and 72 °C for 30 s, with a final extension at 72 °C for 10 min. Duplicate first step PCR products were pooled then used as template for the second step PCR. In the second step, PCR amplification added multiplexing index-sequences to the overhang adapters using a unique multiplex primer pair combination for each sample. The reaction was carried out in duplicate 30 µL volumes containing 15 µL Phusion High-Fidelity PCR Master Mix (THERMO FISHER SCIENTIFIC), 1 µM of one forward and one reverse multiplex primer and 6 µL of first step PCR product. Thermal cycling conditions were 98 °C for 3 min followed by 8 cycles of 98 °C for 30 s, 55 °C for 30 s and 72 °C for 30 s, with a final extension at 72 °C for 10 min. Duplicate second step PCR products were pooled and then visualized in 2% agarose gel to verify amplification and their sizes. The amplicons were cleaned-up and pooled using SequalPrep Normalization plate kit 96-well (INVITROGEN). Sequencing was performed on MiSeq (ILLUMINA, 2×250 bp) using the MiSeq reagent kit v2 (500cycles). Demultiplexing and trimming of ILLUMINA adapters and barcodes was done with ILLUMINA MiSeq Reporter software (version 2.5.1.3). Fungal ITS rRNA region was amplified similarly, using the primers ITS3F (5'-GCA TCGATGAAGAACGCAGC-3') and ITS4R (5'-TCCTCSSCTTATTGATATGC-3'), modified from Ref.²⁹, with 30 cycles for the first step PCR and 10 cycles for the second step PCR.

**Quantification of total and N-cycle related microbial communities.** DNA was extracted from 250 mg dry-weight soil for all 240 samples using the DNeasy PowerSoil-htp 96 well DNA isolation kit (QIAGEN, France). Total bacterial communities were quantified using 16S rDNA and ITS primers-based qPCR assays as described in Refs.²⁸,²⁹. Quantification of the bacterial and archaeal ammonia-oxidizers (AOB and AOA respectively) was performed according to Refs.³⁰,³¹ whereas quantification of denitrifiers was performed according to Refs.³²,³³,³⁴. For this purpose, the genes encoding catalytic enzymes involved in nitrogen fixation (*nifH*), ammonia-oxidation (bacterial and archaeal *amoA*), nitrite reduction (*nirK* and *nirS*) and nitrous oxide reduction (*nosZI* and *nosZII*) were used as molecular markers. Reactions were carried out in a ViiA7 (LIFE TECHNOLOGIES, USA). Quantification was based on the increasing fluorescence intensity of the SYBR Green dye during amplification. The real-time PCR assays were carried out in a 15 µL reaction volume containing SYBR green PCR Master Mix (Absolute Blue QPCR SYBR Green Low Rox Mix, THERMO FISHER SCIENTIFIC, France), 1 µM of each primer, 250 ng of T4 gene 32 (QBIOGENE, France) and 0.5 ng of DNA as previously described³⁵. Three independent replicates were considered for real time PCR assays. Standard curves were done using serial dilutions of linearized plasmids containing appropriate cloned targeted genes from bacterial strains or environmental clones. PCR efficiency for the different assays ranged from 70 to 99%. No template controls gave null values. The presence of PCR inhibitors in DNA extracted from soil and rhizosphere samples was estimated by mixing a known amount of standard DNA with soil DNA extract before qPCR. We did not detect inhibition in any case.

**Bioinformatic analysis of the 16S rRNA and ITS amplicons.** Sequence data were analysed using an in-house developed Python notebook piping together different bioinformatics tools (available upon request). Briefly, for both 16S and ITS data, sequences were assembled using PEAR³⁶ with default settings. Further quality checks were conducted using the QIIME pipeline³⁷ and short sequences were removed (<335 bp for 16S and <310 bp for ITS). Reference based and de novo chimera detection, as well as clustering in OTUs were performed using VSEARCH³⁸ and the adequate reference databases (Greengenes' representative set of sequences for 16S and UNITE's ITS2 reference dynamic dataset for ITS). The identity thresholds were set at 97% for 16S and ITS data. For 16S rRNA data, representative sequences for each OTU were aligned using PYNAST³⁹ and a 16S rRNA phylogenetic tree was constructed using FASTTREE⁴⁰ Taxonomy was assigned using UCLUST⁴¹ and the latest released Greengenes database (v.05/2013⁴²) for 16S rRNA. For ITS, the taxonomy assignment was performed using BLAST⁴³ and the UNITE reference database (v.7-08/2016⁴⁴). Raw sequences were deposited at the
NCBI under the accession number SUB6424930. The process of raw sequence submission was performed using the make.sra command of MOTHUR software.

Bacterial and fungal α-diversity metrics describing the structure of microbial communities included indices pertaining to richness (observed species, chao1), evenness (Simpson reciprocal), both (Shannon), and relatedness (Faith's Phylogenetic Diversity PD; (Faith, 1992) for bacterial community data only) and were calculated based on rarefied OTU tables (12,000 sequences and 20,000 sequences for 16S rRNA and ITS respectively). Bray-Curtis dissimilarity matrices were also computed to detect global variations in the composition of microbial communities for 16S rRNA and ITS. Bacterial co-occurrence networks were constructed for bulk soil, and the four subspecies separately in order to detect associations between bacteria that might have been lost or might have emerged during the course of domestication. Co-occurrence or exclusion between two OTUs was considered as relevant when Spearman correlation coefficient was significant and > 0.6 or < -0.6. Because Spearman correlations are known to be sensitive to zeros, only the most abundant OTUs were kept for this analysis. Also, several other metrics could have been used to detect correlations between OTUs, however, the aim of the present analysis being to compare networks across the four domestication groups, we believe that the choice of the metrics is less important than if it was to compare it to existing knowledge.

**Statistical analyses.** The analyses of the different response variables (diversity indices, QPCRs abundance data, OTUs relative abundance data) was done using general mixed models. Response variables were log10-transformed prior to analysis. To detect differences in response variables related to the environmental compartment (rhizosphere versus bulk soil), the following model was used:

$$Y_{ijk} = \mu + Block_{i} + compartment_{j} + \varepsilon_{ijk}$$

where $Y_{ijk}$ is the response variable, Block is the random Block effect (i=1, ..., 5), compartment is a fixed compartment effect (j=1, 2) and $\varepsilon_{ijk}$ is the residuals.

To detect differences in response variables related to domestication, or to the host genotypic variation, the following model was used on rhizosphere data:

$$Y_{ijkl} = \mu + Block_{i} + subspecies_{j} + Genotype(subspecies)_{jk} + \varepsilon_{ijkl}$$

where $Y_{ijkl}$ is the response variable, Block is the random Block effect (i=1, ..., 5), subspecies is a fixed subspecies effect (j=1, ..., 4) depicting the domestication effect, Genotype(subspecies) is the random genotypic effect nested within each subspecies (k=1, ..., 10) and $\varepsilon_{ijkl}$ is the residuals. For all studied response variables, normality and homoscedasticity of the residuals distribution was visually inspected. Note that we opted for a custom mixed model analysis to take into account the relative complexity of our experimental design that cannot be handled by classical analysis packages such as edgeR or DESeq2. A False Discovery Rate procedure was used to decrease substantially the probability to detect false positive. Broad sense heritability estimates were calculated as the part of variance explained by the Genotype (subspecies) effect divided by the total variance.

Constrained Principal Coordinate Analyses were performed to retrieve the amount of variation of total bacterial and fungal communities explained by the different explaining variables (compartment, subspecies, Genotype) when conditioning for Block using the vegan *capscale* function. Finally, Mantel tests were performed to detect correlations between microbial communities data (Bray-Curtis distance matrices) and plant phenotypic data (Euclidean distance matrix).

Integration and combined visualization of microbial and plant phenotypic data sets were realized using the R package *mixOmics*^45 using DIABLO (Data Integration Analysis for Biomarker discovery using a Latent component method for Omics studies) in order to identify correlated key variables in both datasets^46 explaining differences between subspecies.

## Results and discussion

**Wheat rhizosphere shapes the soil microbiota.** Like many other plant species, wheat affected the microbiome assembly in the rhizosphere. Quantification of the abundance of the bacterial communities showed significant differences between the wheat rhizosphere and the soil compartments (p < 0.05). We also observed a significant effect of the compartment (rhizosphere soil vs bulk soil) on the number of bacterial species (observed OTUs, p < 0.01) as well as on the Simpson reciprocal index for both bacterial and fungal communities (p < 0.01). In addition to the changes in α-diversity,

clear differences between

compartments in microbial community composition were also found (p < 0.05, Supplementary Fig S1a,b). At the class level, several bacterial taxa exhibited contrasted relative abundances between the soil and the rhizosphere (Supplementary Fig S1c,d). Rhizosphere samples were enriched in members of the α- and γ-proteobacteria, flavobacteria, cytophagia and actinobacteria classes and depleted in members of the gemmatimonadetes, acidobacteria phylum and thermoleophilia class compared to bulk soil. Among fungi, members of the Euriomycetes class were more abundant in the rhizosphere while members of the Leotiomycetes, Pezizomycetes, Tremellomycetes and Mortierellomycotina were over-represented in the soil. These results confirm previous findings that the rhizosphere is a selective environment which favours growth of potentially adapted or fast growing microorganisms^17,47.

**Domestication affected specific bacterial and fungal taxa in the rhizosphere.** We detected no changes in bacterial and fungal α-diversity indices related to domestication contrary to what is classically observed in the literature, probably related to larger number of accessions within each domestication group. A
small but significant proportion of the total variance in the structure of the bacterial and fungal communities was attributable to domestication (4.1%, p = 0.001 and 3.4%, p = 0.002 respectively for bacteria and fungi; Fig. 1a,b). The sizes of the total bacterial and fungal communities were equivalent across the four different wheat groups. AOB was the only N-cycle guild whose abundance was significantly impacted by domestication with a decrease along the domestication gradient (wild emmer > cultivated emmer > durum wheat elite genotypes, Fig. 1c), with the durum landraces genotypes being in-between wild emmer and elite genotypes. When analysing the relative abundance of members of the fungal community, we found two clear contrasted patterns across the domestication series. The first pattern, consisting in a strong decrease in relative abundance between wild emmer and the three domestication groups, was detected for four members of the Sordariomycetes class: the *Togniniella*, the *Chaetosphaeria*, the *Ophiostoma* and the *Raffaela* genera, as well as for one member of the Chytridiomycetes class: the *Clydaea* genus (Fig. 1d, Supplementary Fig. S2a). Interestingly, several members of these genera have been found to be pathogens associated to plant diseases (dutch elm disease for *Ophiostoma*<sup>48</sup>, oak and laurel wilt disease for *Raffaela*<sup>49,50</sup>, or grapevine trunk disease for *Togniniella*<sup>51</sup>). On the opposite, most of the members of the Glomeromycetes class showed the second pattern with increased abundance in cultivated emmer compared to wild emmer. This was the cases of the *Clareidoglomus*, the *Archaeospora* and an unclassified genera belonging to the Glomerales order (Fig. 1f, Supplementary Fig. S2b). All members of this order are obligate plant symbionts, known as arbuscular mycorrhizal fungi (AMF), which are colonizing host plant root for reciprocal nutrient exchange<sup>52</sup>. Our findings suggest that the association between AMF and cultivated emmer increased during the first domestication events, at the expense of potentially pathogenic fungi. Accordingly, previous studies suggested that AMF can provide systemic protection against a wide range of pathogens in several plants through "mycorrhiza-induced resistance"<sup>14,53</sup>. Interestingly, another member of the Glomeromycetes class, *Paraglomus*, displayed a strong abundance decrease between wild emmer compared to the three domestication groups (Supplementary Fig. S2a), which supports previous findings suggesting that *Paraglomus* spp. have distinct ecological characteristics compared to members of the glomerales order<sup>54</sup>.

About 20 bacterial phylotypes displayed significant variation in their relative abundance across the four tetraploid wheat domestication groups. Among those, two out of the three most abundant ones were members of the Flavobacteriaceae and Comamonadaceae families whose abundance were higher in wild emmer, cultivated emmer and elite lines than in landraces accessions. The third phylotype, the most abundant one, was classified as similar to the Gemm 1 clone (Fig. 1e) of the Gemmatimonadetes phylum. The relative abundance of this phylotype displayed a significantly increased abundance in modern durum cultivars compared to durum landraces. In contrast to Gemmatimonadetes, Comamonadaceae and Flavobacteria are often described as bacterial phylotypes exhibiting a significant rhizosphere effect<sup>55-57</sup>.

**Assembly co-occurrence networks across wheat domestication groups of the rhizosphere bacterial community.** While no significant associations between fungal OTUs have been observed, a comparison of the co-occurrence networks inferred for bacterial community from the four domestication groups and bulk soil revealed that a large part of those network members are shared between the four groups (Fig. 2 and Supplementary Fig. S3). Specifically, 24 bacterial nodes were present across all four subspecies networks amongst which 8 were unique to the rhizosphere (i.e. absent in the bulk soil network, Supplementary Fig. S3). Interestingly, most of the bacterial phylotypes whose abundances have been shown to be dependent on the host genotype are present in those co-occurrence networks. While the majority of bacterial nodes composing the bulk soil network are bulk soil-specific, very few OTUs are domestication groups-network specific (none for *T. turgidum* spp. *dicoccoides*, 7 out of 41 for *T. turgidum* spp. *dicoccum*, 7 out of 51 for durum landraces and 1 out of 64 for durum elite accessions) which tends to indicate that domestication did not drastically change bacterial association patterns in the rhizosphere (Supplementary Fig. S3b).

**Aside from domestication, genotypic variability is a strong driver of rhizosphere bacterial and fungal community composition.** A significant part of the observed variations in bacterial and fungal α-diversity indices was attributable to genotype variation within each subspecies. The broad sense heritability of the bacterial observed species and Simpson reciprocal indices were respectively 16.8% and 17.6% (p = 0.016 and p = 0.011). In contrast, estimates were not significant for the fungal observed species and Simpson reciprocal index. However, the composition of both bacterial and fungal communities in the rhizosphere was significantly under the influence of the host genotype. About 30% of the total variation in rhizosphere microbial communities is attributable to genotype variation (CAP analyses of Bray-Curtis distance matrices, 28.5%, p = 0.001 and 30.3%, p = 0.001 respectively for bacteria and fungi; Fig. 3a,b). At the global scale, within-domestication group variation in rhizosphere microbial communities is therefore stronger than variation observed between domestication groups. The values observed in this study were higher than those previously reported<sup>16,17</sup>, probably because of the larger genotypic variance existing in our panel, as well as the homogeneous conditions over the experiments. Seven bacterial phylotypes displayed strong and significant broad sense heritability after correcting for multiple testing (FDR procedure, q<0.05). Those phylotypes belong to the Actinobacteria (two of them), the Gemmatimonadetes¹, the Nitrospirae¹, the Proteobacteria² and the Verrucomicrobia¹ phyla. Among those, members of the Rhizobiales order exhibited the highest heritability estimates with 28.2% of their relative abundance explained by differences in host genotype. Interestingly, Rhizobiales are well-known beneficial partners in plant-microbe interactions<sup>58</sup>. Members of the Nitrospiraceae family also showed high heritability estimates (24.8%). Of note, there are several studies describing the role of Nitrospiraceae in nitrogen cycling but,
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Impact of domestication on bacterial and fungal communities composition. Constrained analyses of principal coordinates of bacterial (a) and fungal (b) communities. Black, red, green and blue dots correspond to *T. turgidum* spp. *dicoccoides*, *T. turgidum* spp. *dicoccum*, *T. turgidum* spp. *durum* landraces and *T. turgidum* spp. *durum* elite cultivars, respectively. (c) Mean gene copy numbers (and standard deviations) of ammonia-oxidizing bacteria across the four domestication groups expressed in number of copies per g of DNA. Same colors as above. (d–f) Distribution of counts assigned to Ascomycota:Ophiostoma, Gemmatimonadetes:Gemm1 and Glomeromycota:Claroideoglomus, respectively, across the four domestication groups. Mean and standard errors are represented. Same colors as above.</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>d</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 2.** Bacterial cooccurrence networks across the four domescation groups. (a–d) Correspond to *T. turgidum* spp. *dicoccoides*, *T. turgidum* spp. *dicoccum*, *T. turgidum* spp. *durum* landraces and *T. turgidum* spp. *durum* elite cultivars rhizosphere bacterial cooccurrence networks, respectively. Edges indicate positive (green) or negative correlations (red) as defined by Spearmann's correlation coefficients *ρ* > 0.6 or *ρ* < −0.6, respectively. Node colors differ according to OTUs phylum assignment.</figcaption></figure>

to our knowledge, none has described this taxa as significantly impacted by plant roots. Finally, the Gemm 1 clone within the Gemmatimonadetes phylum is of particular interest because it showed both variations between domescation groups but also significant genotypic variation within each group with a heritability estimate of 23.2%. Unfortunately, despite being found in a variety of soils, little is known about its ecology.

Several fungal phylotypes also displayed significant variation related to the identity of the host genotype. Probably the most interesting one is the Ophiostomataceae phylotype containing the *Ophiostoma* and *Raffaela* genera, which were also impacted by domestication. The relative abundance of Ophiostomataceae family members was strongly dependent on the host genotype (26.6%), especially within the *dicoccoides* subspecies. In this subspecies, the counts variation attributed to this phylotype ranged between ~50 for genotype 46309 to ~1,200 for genotype 46499 while it showed a much smaller range of variation in the three other domestication groups.
<figure><img src="image_9.png" /><figcaption>**Figure 3.** Impact of within-domestication group genotypic variation on bacterial and fungal communities composition. Constrained analyses of principal coordinates of bacterial (a) and fungal (b) communities. Black, red, green and blue dots correspond to *T. turgidum* spp. *dicoccoides*, *T. turgidum* spp. *dicoccum*, *T. turgidum* spp. *durum* landraces and *T. turgidum* spp. *durum* elite cultivars, respectively. Genotypes are encoded with different symbols within each domestication group.</figcaption></figure>

The abundance of the total rhizosphere associated microbial community, as well as those of the N-cycle microbial guilds were not significantly influenced by the host genotypic variation.

**Relationships between plant phenotypes and microbial communities.** We sought to determine whether the variations in the rhizosphere microbiome were related to changes in wheat phenotypic traits during domestication. The phenotypic distance between plants was positively correlated to differences in composition of rhizosphere bacterial communities (Mantel R test, R = 0.25, p = 0.02), but not to the overall fungal community composition (Mantel R test, R = 0.12, p = 0.14). We then combined the microbial community composition data with the plant phenotypic trait dataset in order to detect covariation in those datasets that might be related to domestication using a DIABLO approach. When looking at both microbial communities, a correlation between the most discriminant bacterial and fungal OTUs, and plant phenotypic traits was observed on the 1st component. This 1st component was mainly discriminating *T. t. spp. dicoccoides* accessions from the others while the 2nd one discriminated most of the *T. t. spp. dicoccum* accessions from the others (Fig. 4). We identified nine bacterial and six fungal OTUs that correlated with seven plant phenotypic traits (Table 1). Figure 4 highlights links between those OTUs and plant phenotypic traits indicating negative and positive correlations (|r|>0.4, assessed via DIABLO similarity scores that are analogous to Pearson's correlation coefficients). Plant diameter was the trait exhibiting the strongest associations with microbial taxa, being both positively correlated with the relative abundance of unknown members the Basidiomycota phylum (r = 0.61), and negatively correlated with the relative abundance of members C114 of the Gemmatimonadetes phylum (r = -0.57). Previous study suggested that differences in root architecture can also determine how effective Glomeraceae are at reducing fungal pathogen abundance⁵⁹. This is consistent with our findings showing that potential plant pathogen fungi were both negatively correlated to the percentage of fine roots, namely here Togniniella and Ophiostoma, and to the proportion of beneficial Glomeraceae, namely here Clareidoglomus, Glomus and an unidentified member of the Glomeromycota (Fig. 4). Taken together, our results showed that the abundances of several fungal taxa impacted by domestication were also negatively related to the percentage of fine roots, which indicate that changes in plant morphology, especially in the root architecture during domestication, might have driven changes in the composition of the rhizosphere microbiota.

**Conclusion**

Altogether, our findings indicate that domestication has shaped the rhizosphere associated microbiota in wheat. The impact of domestication was stronger on fungi with a large decrease of potential pathogens associated with an increase of obligate plant symbionts belonging to the glomeromycota phylum, that were related to a change in the morphology of the root system between wild and cultivated emmer. The depletion of microbes potentially important for plant nutrient availability, e.g. nitrifying bacteria as well as arbuscular mycorrhizal fungi, between primitive domesticated forms and modern accessions suggests that the Green revolution, which relied on chemical fertilizers, might have caused a decrease in the association between tetraploid wheat and its rhizosphere microbiota. Nevertheless, the strong genotypic variation that we observed within each domestication group, indicates that there might be room for genetic improvement via introgression of ancestral plant rhizosphere-beneficial microbe associations.
<figure><img src="image_10.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_11.png" /></figure>

**Figure 4.** Network visualization of associations between bacterial and fungal OTUs abundances, and ten above and belowground plant traits across the four domestication groups. Green and blue nodes correspond to plant phenotypic traits and bacterial (a) or fungal (b) OTUs, respectively. Edges indicate positive (red) or negative correlations (blue) as defined by Pearson's correlation coefficients $r > 0.4$ or $r < -0.4$, respectively. Different modules are represented by different colors (green, red, and blue). Left panels are showing discrimination of the samples based only on the sub-selected OTUs and plant traits that correlate.

Received: 25 February 2020; Accepted: 8 July 2020

Published online: 22 July 2020

References

1. Nesbitt, M. & Samuel, D. Wheat domestication: Archaeobotanical evidence. *Science* **279**, 1431–1431 (1998).

2. Hammer, K. Das Domestikationssyndrom. *Kult*. **32**, 11–34 (1984).

3. Peleg, Z., Fahima, T., Korol, A. B., Abbo, S. & Saranga, Y. Genetic analysis of wheat domestication and evolution under domestication. *J. Exp. Bot*. **62**, 5051–5061 (2011).

4. Simons, K. J. et al. Molecular characterization of the major wheat domestication gene Q. *Genetics* **172**, 547 (2006).

5. Simonetti, M. C. et al. Quantitative trait loci influencing free-threshing habit in tetraploid wheats. *Genet. Resour. Crop Evol*. **46**, 267–271 (1999).
<table><caption>**Table 1.** Correlations between plant traits and the relative abundances of microbes related to domestication.</caption><thead><tr><th>Similarity score<sup>a</sup></th><th>Plant diameter</th><th>Diameter</th><th>Leaf longevity</th><th>Root angle</th><th>N mass</th><th>Percent fine roots</th><th>Dry shoot</th></tr></thead><tbody><tr><td>Actinobacteria:Actinomy cetales</td><td>-</td><td>-0.54</td><td>-</td><td>-</td><td>-</td><td>0.49</td><td>-</td></tr><tr><td>Actinobacteria:Streptomyc etaceae</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>0.49</td></tr><tr><td>Ascomycota:Ophiostoma</td><td>0.42</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-0.42</td><td>-</td></tr><tr><td>Ascomycota:Tetrapisispora</td><td>-</td><td>0.42</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Ascomycota:Togniniella</td><td>0.53</td><td>0.40</td><td>-</td><td>-</td><td>-</td><td>-0.52</td><td>-</td></tr><tr><td>Basidiomycota:Sebacina</td><td>-</td><td>0.53</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Basidiomycota:unidentified</td><td>0.61</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Gemmatimonadetes:C114</td><td>-0.57</td><td>-</td><td>-0.43</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Glomeromycota:Paraglomus</td><td>0.47</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-0.49</td><td>-</td></tr><tr><td>Proteobacteria:Legionellales</td><td>-</td><td>-</td><td>0.53</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Proteobacteria:Methylop hilales</td><td>0.41</td><td>-</td><td>-</td><td>-</td><td>-0.55</td><td>-0.41</td><td>-</td></tr><tr><td>Proteobacteria:Oxalobact eraceae</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-0.40</td><td>0.44</td></tr><tr><td>Proteobacteria:Phyllobact riacae</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-0.47</td><td>-</td></tr><tr><td>Proteobacteria:Rhodoferax</td><td>-</td><td>-</td><td>-</td><td>0.41</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Proteobacteria:Sinobacte raceae</td><td>-0.52</td><td>-</td><td>-0.47</td><td>-</td><td>-</td><td>-</td><td>-</td></tr></tbody></table>

6. Tzarfati, R. et al. Threshing efficiency as an incentive for rapid domestication of emmer wheat. *Ann. Bot.* **112**, 829–837 (2013).

7. Matsuoka, Y. Evolution of polyploid triticum wheats under cultivation: The role of domestication, natural hybridization and allopolyploid speciation in their diversification. *Plant Cell Physiol.* **52**, 750–764 (2011).

8. Araus, J. L., Bort, J., Steduto, P., Villegas, D. & Royo, C. Breeding cereals for Mediterranean conditions: Ecophysiological clues for biotechnology application. *Ann. Appl. Biol.* **142**, 129–141 (2003).

9. Gioia, T. et al. Impact of domestication on the phenotypic architecture of durum wheat under contrasting nitrogen fertilization. *J. Exp. Bot.* **66**, 5519–5530 (2015).

10. Thuillet, A.-C., Bataillon, T., Poirier, S., Santoni, S. & David, J. L. Estimation of long-term effective population sizes through the history of durum wheat using microsatellite data. *Genetics* **169**, 1589 (2005).

11. Haudry, A. et al. Grinding up wheat: A massive loss of nucleotide diversity since domestication. *Mol. Biol. Evol.* **24**, 1506–1517 (2007).

12. Pérez-Jaramillo, J. E., Mendes, R. & Raaijmakers, J. M. Impact of plant domestication on rhizosphere microbiome assembly and functions. *Plant Mol. Biol.* **90**, 635–644 (2016).

13. Philippot, L., Raaijmakers, J. M., Lemanceau, P. & van der Putten, W. H. Going back to the roots: The microbial ecology of the rhizosphere. *Nat. Rev. Microbiol.* **11**, 789 (2013).

14. Pozo, M. J. & Azcón-Aguilar, C. Unraveling mycorrhiza-induced resistance. *Curr. Opin. Plant Biol.* **10**, 393–398 (2007).

15. Cook, R. J. et al. Molecular mechanisms of defense by rhizobacteria against root disease. *Proc. Natl. Acad. Sci.* **92**, 4197 (1995).

16. Lundberg, D. S. et al. Defining the core *Arabidopsis thaliana* root microbiome. *Nature* **488**, 86 (2012).

17. Peiffer, J. A. et al. Diversity and heritability of the maize rhizosphere microbiome under field conditions. *Proc. Natl. Acad. Sci.* **110**, 6548 (2013).

18. İnceoğlu, Ö, Salles, J. F., van Overbeek, L. & van Elsas, J. D. Effects of plant genotype and growth stage on the betaproteobacterial communities associated with different potato cultivars in two fields. *Appl. Environ. Microbiol.* **76**, 3675 (2010).

19. Bulgarelli, D. et al. Revealing structure and assembly cues for Arabidopsis root-inhabiting bacterial microbiota. *Nature* **488**, 91 (2012).

20. Hardoim, P. R. et al. Rice root-associated bacteria: Insights into community structures across 10 cultivars. *FEMS Microbiol. Ecol.* **77**, 154–164 (2011).

21. Pérez-Jaramillo, J. E. et al. Linking rhizosphere microbiome composition of wild and domesticated Phaseolus vulgaris to genotypic and root phenotypic traits. *ISME J.* **11**, 2244 (2017).

22. Bulgarelli, D. et al. Structure and function of the bacterial root microbiota in wild and domesticated barley. *Cell Host Microbe* **17**, 392–403 (2015).

23. Leff, J. W., Lynch, R. C., Kane, N. C. & Fierer, N. Plant domestication and the assembly of bacterial and fungal communities associated with strains of the common sunflower, *Helianthus annuus*. *New Phytol.* **214**, 412–423 (2017).

24. Roucou, A. et al. Shifts in plant functional strategies over the course of wheat domestication. *J. Appl. Ecol.* **55**, 25–37 (2018).

25. Röder, M. S. et al. A microsatellite map of wheat. *Genetics* **149**, 2007 (1998).

26. Berry, D., Ben Mahfoudh, K., Wagner, M. & Loy, A. Barcoded Primers used in multiplex amplicon pyrosequencing bias amplification. *Appl. Environ. Microbiol.* **77**, 7846 (2011).

27. Takahashi, S., Tomita, J., Nishioka, K., Hisada, T. & Nishijima, M. Development of a prokaryotic universal primer for simultaneous analysis of bacteria and archaea using next-generation sequencing. *PLoS ONE* **9**, e105592 (2014).

28. Muyzer, G., de Waal, E. C. & Uitterlinden, A. G. Profiling of complex microbial populations by denaturing gradient gel electrophoresis analysis of polymerase chain reaction-amplified genes coding for 16S rRNA. *Appl. Environ. Microbiol.* **59**, 695 (1993).

29. White, T. J., Bruns, T., Lee, S. & Taylor, J. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In *PCR Protocols* (eds Innis, M. A. et al.) 315–322 (Academic Press, Cambridge, 1990).

30. Tourna, M., Freitag, T. E., Nicol, G. W. & Prosser, J. I. Growth, activity and temperature responses of ammonia-oxidizing archaea and bacteria in soil microcosms. *Environ. Microbiol.* **10**, 1357–1364 (2008).

31. Leininger, S. et al. Archaea predominate among ammonia-oxidizing prokaryotes in soils. *Nature* **442**, 806–809 (2006).

<footnote><sup>a</sup>Note that similarity scores in DIABLO are analogous to Pearson's correlation coefficients.</footnote>
32. Henry, S. et al. Quantification of denitrifying bacteria in soils by nirK gene targeted real-time PCR. *J. Microbiol. Methods* **59**, 327–335 (2004).

33. Henry, S., Bru, D., Stres, B., Hallet, S. & Philippot, L. Quantitative detection of the *nosZ* Gene, encoding nitrous oxide reductase, and comparison of the abundances of 16S rRNA, *narG*, *nirK*, and *nosZ* genes in soils. *Appl. Environ. Microbiol.* **72**, 5181 (2006).

34. Jones, C. M., Graf, D. R., Bru, D., Philippot, L. & Hallin, S. The unaccounted yet abundant nitrous oxide-reducing microbial community: A potential nitrous oxide sink. *ISME J.* **7**, 417 (2012).

35. Bru, D. et al. Determinants of the distribution of nitrogen-cycling microbial communities at the landscape scale. *ISME J.* **5**, 532 (2010).

36. Zhang, J., Kobert, K., Flouri, T. & Stamatakis, A. PEAR: A fast and accurate Illumina Paired-End reAd mergeR. *Bioinformatics* **30**, 614–620 (2013).

37. Caporaso, J. G. et al. QIIME allows analysis of high-throughput community sequencing data. *Nat. Methods* **7**, 335 (2010).

38. Rognes, T., Flouri, T., Nichols, B., Quince, C. & Mahé, F. VSEARCH: A versatile open source tool for metagenomics. *PeerJ* **4**, e2584 (2016).

39. Caporaso, J. G. et al. PyNAST: A flexible tool for aligning sequences to a template alignment. *Bioinformatics* **26**, 266–267 (2009).

40. Price, M. N., Dehal, P. S. & Arkin, A. P. FastTree 2—Approximately maximum-likelihood trees for large alignments. *PLoS ONE* **5**, e9490 (2010).

41. Edgar, R. C. Search and clustering orders of magnitude faster than BLAST. *Bioinformatics* **26**, 2460–2461 (2010).

42. McDonald, D. et al. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. *ISME J.* **6**, 610 (2011).

43. Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. Basic local alignment search tool. *J. Mol. Biol.* **215**, 403–410 (1990).

44. Nilsson, R. H. et al. The UNITE database for molecular identification of fungi: Handling dark taxa and parallel taxonomic classifications. *Nucleic Acids Res.* **47**, D259–D264 (2019).

45. Rohart, F., Gautier, B., Singh, A. & Cao, K.-A.L. mixOmics: An R package for 'omics feature selection and multiple data integration. *PLoS Comput. Biol.* **13**, e1005752 (2017).

46. A. Singh, et al., DIABLO—An integrative, multi-omics, multivariate method for multi-group classification. *bioRxiv*, 067611 (2016).

47. Fierer, N., Bradford, M. A. & Jackson, R. B. Toward an ecological classification of soil bacteria. *Ecology* **88**, 1354–1364 (2007).

48. Anderson, P. K. et al. Emerging infectious diseases of plants: Pathogen pollution, climate change and agrotechnology drivers. *Trends Ecol. Evol.* **19**, 535–544 (2004).

49. Dreaden, T. J. et al. Development of multilocus PCR assays for *Raffaelea lauricola*, causal agent of laurel wilt disease. *Plant Dis.* **98**, 379–383 (2013).

50. Takahashi, Y., Matsushita, N. & Hogetsu, T. Spatial distribution of *Raffaelea quercivora* in xylem of naturally infested and inoculated oak trees. *Phytopathology* **100**, 747–755 (2010).

51. Gramaje, D., Armengol, J. & Ridgway, H. J. Genetic and virulence diversity, and mating type distribution of *Togninia minima* causing grapevine trunk diseases in Spain. *Eur. J. Plant Pathol.* **135**, 727–743 (2013).

52. Corradi, N. & Bonfante, P. The arbuscular mycorrhizal symbiosis: Origin and evolution of a beneficial plant infection. *PLoS Pathog.* **8**, e1002600 (2012).

53. Maherali, H. & Klironomos, J. N. Influence of phylogeny on fungal community assembly and ecosystem functioning. *Science* **316**, 1746–1748 (2007).

54. Gosling, P., Proctor, M., Jones, J. & Bending, G. D. Distribution and diversity of *Paraglomus* spp. in tilled agricultural soils. *Mycorrhiza* **24**, 1–11 (2014).

55. B. D. Emmett, N. D. Youngblut, D. H. Buckley, L. E. Drinkwater, plant phylogeny and life history shape rhizosphere bacterial microbiome of summer annuals in an agricultural field. *Front. Microbiol.* **8**, 2414 (2017).

56. Hernández, M., Dumont, M. G., Yuan, Q. & Conrad, R. Different bacterial populations associated with the roots and rhizosphere of rice incorporate plant-derived carbon. *Appl. Environ. Microbiol.* **81**, 2244–2253 (2015).

57. DeBruyn, J. M., Nixon, L. T., Fawaz, M. N., Johnson, A. M. & Radosevich, M. Global biogeography and quantitative seasonal dynamics of gemmatimonadetes in soil. *Appl. Environ. Microbiol.* **77**, 6295 (2011).

58. Oldroyd, G. E. D. Speak, friend, and enter: Signalling systems that promote beneficial symbiotic associations in plants. *Nat. Rev. Microbiol.* **11**, 252 (2013).

59. Sikes, B. A. When do arbuscular mycorrhizal fungi protect plant roots from pathogens?. *Plant Signal. Behav.* **5**, 763–765 (2010).

## Acknowledgements

This study was supported by the European Research Council (ERC) Starting Grant Project "Ecophysiological and biophysical constraints on domestication in crop plants" (Grant ERC-StG-2014-639706-CONSTRAINTS) and by an INRAE Pari Scientifique (IDIOME-Blé).

## Author contributions

A.S., D.V., P.R., L.P. and C.V. designed the experiment. A.R., D.B., M.C.B., F.F., D.V., P.R. and C.V. performed the experiment. A.S. and A.M. analysed the data. A.S. and L.P. wrote the manuscript with inputs from A.R., F.F., D.V., P.R. and C.V. All authors reviewed the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

Supplementary information is available for this paper at https://doi.org/10.1038/s41598-020-69175-9.

Correspondence and requests for materials should be addressed to A.S. or L.P.

Reprints and permissions information is available at www.nature.com/reprints.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
<figure><img src="image_12.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or
format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the
Creative Commons license, and indicate if changes were made. The images or other third party material in this
article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the
material. If material is not included in the article's Creative Commons license and your intended use is not
permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from
the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020