Identification of barley genetic regions influencing plant–microbe interactions and carbon cycling in soil

Rhizodeposition shapes soil microbial communities that perform important processes such as soil C mineralization, but we have limited understanding of the plant genetic regions influencing soil microbes. Here, barley chromosome regions affecting soil microbial biomass-C (MBC), dissolved organic-C (DOC) and root biomass were characterised. A quantitative trait loci analysis approach was applied to identify barley chromosome regions affecting soil MBC, soil DOC and root biomass. This was done using barley Recombinant Chromosome Substitution Lines (RCSLs) developed with a wild accession (Caesarea 26-24) as a donor parent and an elite cultivar (Harrington) as recipient parent. Significant differences in root-derived MBC and DOC and root biomass among these RCSLs were observed. Analysis of variance using single nucleotide polymorphisms genotype classes revealed 16 chromosome regions influencing root-derived MBC and DOC. Of these chromosome regions, five on chromosomes 2H, 3H and 7H were highly significant and two on chromosome 3H influenced both root-derived MBC and DOC. Potential candidate genes influencing root-derived MBC and DOC concentrations in soil were identified. The present findings provide new insights into the barley genetic influence on soil microbial communities. Further work to verify these barley chromosome regions and candidate genes could promote marker assisted selection and breeding of barley varieties that are able to more effectively shape soil microbes and soil processes via rhizodeposition, supporting sustainable crop production systems.


Introduction
Soil microbes mediate carbon (C) and organic matter cycling in soil, contributing a vital role for the regulation of CO 2 emissions from soil (Prentice et al. 2001;Li et al. 2013) and nutrient release from soil organic matter (SOM) (Fontaine et al. 2011;Dijkstra et al. 2013;Alegria Terrazas et al. 2016). The soil microbial communities and their interactions with plants are impacted by the release of a range of compounds from living roots through root exudation, sloughed cells, mucilage and so on, collectively defined as rhizodeposition (Jones et al. 2004). These rhizodeposit compounds, in particular root exudates, are utilized by microbes as C sources to derive energy for their activity (Paterson 2003;Cheng and Kuzyakov 2005), with this resulting in the decrease or increase of SOM decomposition (Jenkinson et al. 1985;Kuzyakov et al. 2000;Yin et al. 2019). Indeed, it is known that the growth of plants can alter (via rhizodeposition) this microbially mediated SOM decomposition to varying extents (e.g. Cheng et al. 2003;Mwafulirwa et al. 2016Mwafulirwa et al. , 2021, with increases of up to 380% relative to unplanted soil reported by Cheng et al. (2003).
There is growing evidence that plant influences on soil microbial communities, and the functions they undertake, vary not only between plant species but also between individual genotypes within a single plant species. For example, studies by Aira et al. (2010), Bouffaud et al. (2012), Peiffer et al. (2013) and Walters et al. (2018) suggest that rhizosphere microbial community composition under maize is related to plant genotype. In barley, our previous findings (Mwafulirwa et al. 2016(Mwafulirwa et al. , 2017 and those of Pausch et al. (2016) are indicative that soil microbial activity and, in turn, the decomposition of SOM are also impacted by plant genotype. However, there is a lack of knowledge about the plant genetic regions and plant genes influencing these plant-microbe interactions. Identifying the plant genes influencing these interactions could, in particular, underpin the breeding of crop varieties to control microbially mediated soil processes (such as C mineralization). Hence, better understanding of the plant genes influencing (via rhizodeposition) soil microbes and SOM dynamics could help inform crop breeding to support sustainable agricultural production.
A major limitation for crop breeding to control soil processes, and in turn agricultural and/or environmental sustainability, is the current loss of beneficial plant traits associated with soil microbial interactions in the elite gene pool. This is because the development of modern crop cultivars through selection for yield and other beneficial crop plant traits (such as tolerance to biotic and abiotic stresses), usually under intensive chemical fertilizer applications to soil, have resulted in the loss of the plant genetic variation influencing plant-soil interactions (Tanksley and McCouch 1997;Wissuwa et al. 2009). To overcome this problem, others have proposed the use of wild relatives of crop species as donors of exotic germplasm to improve elite varieties. For example, Matus et al. (2003) developed a population of Recombinant Chromosome Substitution Lines (RCSLs) using wild barley Hordeum vulgare subsp. spontaneum (Caesarea 26-24) as a donor and Hordeum vulgare subsp. vulgare (Harrington, a North American malting cultivar) as the recurrent parent. Importantly, Caesarea 26-24 is adapted to specific soil conditions, in the view that it was collected in a dry and saline environment in Israel (Matus et al. 2003). This suggests that this accession could provide vital genes for regulating plant-soil interactions, especially for stress adaptation/tolerance promoting sustainable production. Our previous work investigating the impacts of plant intraspecific variation on SOM decomposition using a small number of these RCSLs showed that plant genotype influenced microbial activity and soil functioning (Mwafulirwa et al. 2016). The genetic composition of these RCSLs (Close et al. 2009;Comadran et al. 2012;de la Fuente Cantó et al. 2018) and a minimum set representing the entire genome of the wild donor parent (de la Fuente Cantó et al. 2018) have been previously determined using mapped single nucleotide polymorphisms (SNPs). The usefulness of minimum sets of introgression lines representing the entire genome of the donor parent for the detection of quantitative trait loci (QTL) effects and identification of novel exotic alleles is demonstrated by several QTL studies in crop plants, facilitating rapid screening of genetic variation for traits requiring detailed or complex phenotypic evaluations (Schmalenbach et al. 2009;de la Fuente Cantó et al. 2018). For example, Prudent et al. (2009) andTripodi et al. (2020)  Here we used these 28 barley RCSLs, that together represent the whole genome of the wild donor Caesarea 26-24 in the genetic background of the elite variety Harrington, to assess the variation in selected root and soil microbe related traits, i.e. root biomass, dissolved organic-C (DOC) and microbial biomass-C (MBC). Previous studies (Blagodatskaya et al. 2009;Tian et al. 2012;Mwafulirwa et al. 2016) showed that soil DOC and MBC are strongly related to soil respiration rates, and thus to soil functioning. The specific objectives of the present study were to (i) determine the range of variation of the impacts of the barley RCSLs on DOC in soil solution and soil MBC using a 13 C-CO 2 isotopic labelling approach, and (ii) apply a QTL analysis approach with these RCSLs to identify the barley chromosome regions and potential candidate genes influencing soil DOC and MBC.

Soil type
The soil was sampled from a conventionally managed field at Balruddery farm (56.4837° N, 3.1314° W) near Dundee, Scotland, from a depth of 0-10 cm and was sieved to < 6 mm onsite before storing at 4 °C for one week. The soil was a sandy loam of Balrownie Series, Balrownie Association, as identified by Bell et al. (2014, unpublished), and had an organic matter content of 5.8% (muffle furnace, 450 °C, 24 h), pH of 6.0 (H 2 O) and water content (w/w) of 16.9%. The field was planted with barley that was at vegetative (i.e. stem elongation) stage during soil sampling.

Plant materials and genotyping
Twenty-eight barley RCSLs developed using a wild donor (Caesarea 26-24, from a dry and saline region in Israel) and an elite cultivar (Harrington, a North American malting variety) as recipient parent via an advanced backcrossing strategy were used (Matus et al. 2003). These lines represent the minimum number covering the entire wild donor genome, and each line carries a small introgression of the wild barley genome in the predominantly elite background (Generation Challenge Program, unpublished). These lines were chosen because (i) they represent a unique source of genetic diversity to study plant-soil interactions (since they were derived using a wild barley accession adapted to unique soil and environmental conditions), and (ii) selected lines from this population showed differences in rhizodeposition-derived C and the respective impacts on MBC, DOC and SOM mineralization, as observed in our earlier studies (Mwafulirwa et al. 2016(Mwafulirwa et al. , 2017. The genetic architecture of each line was determined from earlier work (de la Fuente Cantó et al. 2018), using the barley iSelect SNP chip (Comadran et al. 2012). Both parent genotypes (Caesarea 26-24 and Harrington) were used for phenotype evaluation, while only the elite parent Harrington was used as a control for genetic analysis.
Experimental setup and 13 C labelling Soil was packed in 93 pots (22.5 cm × 5.5 cm) to a bulk density of 1 g cm −3 and adjusted to 60% water holding capacity (the soil packing volume in each pot was 20.0 cm × 5.5 cm). After one week of soil stabilization in pots, soil solution samplers (Rhizon™ SMS, Rhizosphere Research Products, Wageningen, Netherlands) were inserted to 10 cm depth for DOC measurements. The system was left to stabilize to conditions used in the experiment for a further week before planting.
The barley plants were grown over 39 days without fertilizer addition to soil. Each pot was planted with one of the 30 genotypes (28 RCSLs plus parental genotypes), with one plant per pot, and fallow pots with soil only were included as a control treatment. These were replicated (n = 3) in a randomized complete block design under controlled environment conditions within a plant growth chamber (Conviron CG90; Winnipeg, Canada) at 22 °C and 70% relative humidity. Soil water content was maintained by adding deionized water on a mass basis twice a week. A 12 h daily photoperiod was set with 512 µmol m −2 s −1 PAR within the chamber. Watering was done during the dark period to avoid disruption to the labelling atmosphere.
Labelling plants with 13 C-CO 2 started at the seedling growth stage, one week after sowing seeds. This was achieved by passing a continuous flow of 13 C-enriched CO 2 (20 atom% 13 C) through the plant growth chamber continuously over the experiment period, which was derived by blending CO 2 -free air routed via pressure swing adsorption CO 2 scrubber unit (Parker Balston, Haverhill, USA) with 99 atom% excess 13 C-CO 2 (Sercon Ltd., Cheshire, UK) and CO 2 from a standard CO 2 cylinder (BOC, Worsley, UK) via Brooks thermal mass flow controllers (Flotech Solutions Ltd., Stockport, UK).

Phenotyping
Measurement of DOC was done at 30d and 39d (harvest point), MBC was measured at 39d, and plant biomass (separated into root and shoot biomass) was quantified at 39d. For DOC sampling, vacuum sealed 10 mL bottles were connected to Rhizon soil solution samplers via needle ends and left overnight. The vacuum allowed the soil solution to be drawn into the bottle. Bottles were detached and the collected solutions were kept frozen until they could be analysed for total organic C (TOC) concentration (Shimadzu TOC analyser, Japan) and 13 C-enrichment of this TOC. Soil solution TOC was assumed to represent total DOC. The 13 C-enrichment of DOC was determined using a method described by Garcia-Pausas and Paterson (2011). In brief, the frozen soil solutions were defrosted and sufficient solutions to contain a minimum of 20 µg C and a maximum of 60 µg C per sample, as established from previous TOC analysis of the solutions, were dispensed into mufflefurnanced 12 mL Exetainer vials (Labco Ltd., High Wycombe, UK). Maximum volume in the Exetainers was 4 mL, and deionized water was added to reach the 4 mL volume if a lesser volume was taken. An aliquot of 100 µL of 1.3 M phosphoric acid was added to each sample to remove dissolved inorganic C from the solution. The soil solutions were left with caps off for one hour, following which the vials were capped and flushed with CO 2 -free air on the gas bench. To evolve organic C as CO 2 , an aliquot of 100 µL 1.05 M sodium persulphate was injected into each solution sample through the rubber septum of the vial-cap. The samples were then heated on a dry block at 90 °C for 30 min to release CO 2 . The CO 2 released from the soil solution was then transferred to pre-evacuated N 2 flush-filled Exetainer vials using a syringe with a flow control valve. The 13 C-enrichment of the CO 2 sample was determined on a Delta PLUS Advantage isotope ratio mass spectrometer via an interfaced Gasbench II unit (both Thermo Finnigan, Bremen, Germany). The measured 13 C-enrichment (atom% 13 C) of soil solutions were used to separate root-derived DOC (DOC plant ) and SOM-derived DOC (DOC soil ) proportions of total DOC (DOC total ) following Eqs. 1 and 2 (Garcia-Pausas and Paterson 2011).
where atom% 13 C control is the mean atom% 13 C of DOC measured in the unlabelled fallow control treatments and atom% 13 C total is the measured atom% 13 C value of sample total DOC. Atom% 13 C plant is the atom% 13 C value of the plant tissue (described below).
At harvest, soils were re-wetted to initial moisture level by adding deionized water. Plant shoots were harvested by cutting at the soil surface. Roots were carefully removed from soil (shaking off most of the adhering soil) and put in separate containers then washed with deionized water. The harvested plant shoot and root fractions were taken for freeze-drying. The soil was harvested as one fraction, considering that the soil was densely colonised by roots. The fresh harvested soil was thoroughly mixed by hand and immediately stored at 4 °C for subsequent analyses of soil MBC. The dry weights of root and shoot fractions were used to quantify root and shoot biomass, and total plant biomass was calculated as root (1) DOC plant = DOC total atom% 13 C control −atom% 13 C total ∕ atom% 13 C control −atom% 13 C plant (2) DOC soil = DOC total − DOC plant plus shoot biomass. Dried root and shoot samples were ball-milled (Retsch Ball Mill, model MM2000) and analysed for 13 C-enrichment on a Flash EA 1112 Series Elemental Analyser connected via a Conflo III to a Delta Plus XP isotope ratio mass spectrometer (all Thermo Finnigan, Bremen, Germany).
The fresh soil samples were used to determine the soil MBC concentration, within 24 h following harvest. Two sub samples of the fresh harvested soil were used for MBC analysis. MBC was determined by the chloroform fumigation-extraction method (Vance et al. 1987), where fresh fumigated and nonfumigated soil samples (equivalent 12.5 g dry soil) were extracted with 50 mL of 0.5 M K 2 SO 4 solution. Organic C of the extracts was analysed on a TOC Analyser 700 (Corporation College Station, TX). MBC was calculated as the difference between organic C in the paired fumigated and non-fumigated extracts using a conversion factor k EC of 0.45 (Eq. 3) (Joergensen 1996).
where TOC Fumigated is TOC of fumigated soil sample and TOC Non-fumigated is TOC of the paired nonfumigated soil sample. We further determined the 13 C-enrichment of MBC using a method outlined above for determining the 13 C-enrichment of DOC (Garcia-Pausas and Paterson 2011) and calculated the fractions of MBC derived from plant and SOM using Eqs. 1 and 2.

Statistical and genetic analyses
The software package GenStat (Eighteenth Edition, VSN International Ltd) was used for all statistical analyses. Repeated-measures analysis of variance (ANOVA) was used to assess the effects of barley genotype and sampling date on soil solution DOC (i.e. total DOC, root-derived DOC and SOM-derived DOC), with barley genotype as the fixed factor and sampling date as the repeated factor. In addition, oneway ANOVA was used to test for differences in soil MBC and plant biomass (including root and shoot biomass) among genotypes at harvest point. Where statistically significant (p < 0.05) genotype effects were found, the least significant difference (LSD) was used to assess differences between individual means.
(3) MBC = TOC Fumigated − TOC Non−fumigated ∕k EC For identification of barley chromosome regions influencing the observed phenotypes (i.e. root-derived MBC, root-derived DOC and root biomass, all at 39d), one-way ANOVA was used to analyse the genotype means for each trait to test whether these were related to the genotype for each SNP on the map, with 235 SNPs used in the analysis. Each of these markers represents a block of contiguous SNPs that were found polymorphic for the same RCSLs. We removed redundant markers from the initial group of 1848 SNPs used to characterise the RCSLs. Nevertheless, the mapping information for the group of SNPs within each block was considered to define the QTL regions as explained in de la Fuente Cantó et al. (2018). We present all marker effects that are significant with p < 0.05, but we focus on the most significant ones with p < 0.01 to take into account that multiple markers have been tested. The relative performance of individual RCSLs (RP RCSL ), in comparison to the lines with the genotype of the elite parent Harrington, was calculated following Eq. 4. Caesarea 26-24 (donor parent) was not used in the genetic analysis.
where M RCSL is the trait mean of the RCSL genotypes and M Harrington is the trait mean of the lines with the genotype of the elite parent.
Where two phenotypes or traits were found to be influenced by the same chromosome region (i.e. multiple marker effects), Pearson correlation was used to determine the relationship between those two phenotypes. Pearson correlation was also applied to determine the relationship between root biomass and root-derived MBC.

Identification of candidate genes
Genes which are located within the introgressed chromosome regions that span the observed significant markers and their physical positions were determined from the gene database Barlex (https:// apex. ipk-gater sleben. de/ apex/f? p= 284: 10) and a map-based barley genome assembly with high-confidence genes (Mascher et al. 2017). Based on gene functional annotation, in relation to root-derived MBC and DOC subsets of potential candidates were identified. (4)

Phenotypes
The barley plants were at vegetative (i.e. tillering) stage at the harvest date, 39d after planting, and showed no signs of stress (i.e. pest or pathogen infestation, nor nutrient or water deficiency). Plant tissue 13 C-enrichment ranged from 2.06 to 2.17 atom% 13 C, and did not significantly differ among genotypes. Statistically significant (p < 0.05) differences in root biomass (Fig. 1) and shoot biomass (Fig. S1) among the genotypes were observed.
For soil related characteristics, ANOVA (Tables 1 and S1) showed significant (p < 0.05)  variation among genotypes for root-derived MBC (Fig. 2a) and root-derived DOC (Fig. 2b). The largest root-derived MBC, measured at 39d, was determined as 3.1 mg C kg −1 soil, while the smallest root-derived MBC was estimated as 1.2 mg C kg −1 soil. The estimated root-derived DOC among genotypes ranged from 0.01 to 0.15 mg C L −1 soil solution (0.05 mg C L −1 soil solution average) at 30d and 0.05-0.12 mg C L −1 soil solution (0.08 mg C L −1 soil solution average) at 39d, showing an increasing trend over time. There were no significant differences (Table 1) in total or SOM-derived MBC and total or SOM-derived DOC among the genotypes. Average SOM-derived MBC and average total MBC were 68.20 and 70.30 mg C kg −1 soil, respectively. SOM-derived DOC decreased from an average of 19.97 mg C L −1 soil solution at 30d to 17.11 mg C L −1 soil solution at 39d, while total DOC decreased from an average of 20.12 mg C L −1 soil solution at 30d to 17.20 mg C L −1 soil solution at 39d. Pearson correlation showed a moderate, positive relationship between root-derived DOC and rootderived MBC (r = 0.52; p < 0.001) (Fig. 3a) and a low, positive correlation between root biomass and root-derived MBC (r = 0.44; p < 0.001) (Fig. 3b).

Genetic marker effects and annotated genes
In total, 16 statistically significant (p < 0.05) marker effects (including two multiple marker effects) for root-derived MBC and root-derived DOC, both at . Bars show ± one standard error of the mean. Significant differences (p < 0.05) among barley genotypes and between sampling times were estimated using ANOVA 39d, were detected (Table 2). Root-derived MBC had the highest number of significant marker effects. These marker effects were found on chromosomes 2H, 3H, 4H and 7H. Five of the markers (12_31293, SCRI_RS_129857 and SCRI_RS_171032 on 2H, SCRI_RS_183659 on 3H and 11_21437 on 7H, representing block regions containing 39, 64, 33, 248 and 118 annotated and unknown genes, respectively) were highly significant (p ≤ 0.01) (   were also detected for root-derived MBC (Tables 2  and S2). For these two multiple marker effects, exotic introgressions increased root-derived DOC and rootderived MBC by 38 and 44%, respectively. Seven statistically significant (p < 0.05) markers were identified for root biomass (Table 4). These marker effects were found on chromosomes 2H, 3H, 5H and 7H, with the exotic introgressions increasing root biomass by 29-36%.
A barley genetic map showing statistically significant markers identified for root-derived MBC, rootderived DOC and root biomass is presented in Fig. 4.

Phenotypic evaluation
Analysis of soil microbial biomass showed that concentration of root-derived MBC varied among genotypes. Total MBC and its proportion derived from SOM did not significantly differ among the genotypes. The lack of significant differences in total MBC among genotypes (notwithstanding that rootderived MBC varied among the genotypes) may suggest that barley rhizodeposition impacted microbial activity, including microbial use of root-derived C, rather than microbial growth. This assumption agrees with Shahzad et al. (2015) who showed that root exudation by grassland species could stimulate microbial activity without increasing microbial biomass. It also concurs with Paterson et al. (2008) who suggested that smaller C inputs to soil may not cause significant Markers marked with a dagger ( ‡) were significant for both traits. Those marked with a star (*) showed marker-trait associations with p   changes in the microbial population, and would be processed by those organisms already present in soil. Moreover, this result is consistent with our previous study (Mwafulirwa et al. 2016) where a small number of genotypes from the same barley population did not vary in total MBC and SOM-derived MBC but varied in root-derived MBC, root-derived CO 2 -C, SOMderived CO 2 -C and total CO 2 -C surface soil fluxes. The present results also did not show genotype effects on SOM-derived DOC in soil solutions sampled at 30d and 39d. However, the fraction of total DOC derived from roots varied among the genotypes at both sampling dates. Likewise, these results corroborate our earlier work (Mwafulirwa et al. 2016) where root-derived DOC varied among genotypes (that were selected from the barley population used in the present study) while SOMderived DOC also did not vary among the genotypes. Root-derived DOC increased over time, consistent with plant growth increasing root inputs to soil. In contrast, SOM-derived DOC and total DOC decreased over time, in line with depletion of the available SOM stock. In the present study, it was also noted that concentrations of root-derived DOC among genotypes at 30d did not correlate with concentrations at 39d. This could be explained by restricted root growth because of the use of small pots, in which root growth rates may change with changes in soil nutrients (Marschner 1995) or when all the readily available soil has been explored (Garnett et al. 2009). As such, it is likely that the root-to-soil ratio for individual genotypes varied between the two sampling dates (i.e. 30d and 39d). This potentially affected exudation amounts and thus root-derived DOC concentrations in the soil solutions between those time points. Nevertheless, the pattern of root-derived DOC among genotypes at 30d was consistent with that observed in our previous work using a small number of RCSLs (Mwafulirwa et al. 2016).
Wild barley Caesarea 26-24 had a larger root biomass relative to modern barley Harrington. Indeed, detection of marker effects in the RCSLs (discussed below) showed that the wild barley genome increased root biomass (Table 4). These findings are in agreement with the general understanding that wild barley accessions have an inherent ability to develop vigorous or extensive rooting systems (White et al. 2009;Naz et al. 2014     that lead to greater contact between roots and soil, which in turn enhance water uptake, and thus tolerance to drought (Gahoonia and Nielsen 2004). Naz et al. (2014) (Table 3) were associated with smaller root-derived MBC, relative to modern barley (Harrington). Therefore, this result suggests that Caesarea 26-24 harbours alleles on chromosome 4H that influence microbial use of root-derived C and its flow through the microbial biomass. All other detected marker effects for root-derived MBC showed that the wild barley genome introgressions increased root-derived MBC, but this was in contrast to the observed phenotype (discussed above). This, however, may be indicative of interactive allele or loci effects (Li et al. 2010), or that those effects on MBC were conferred by the modern barley (Harrington) genome. For DOC, the detected marker effects showed that the wild barley genome introgressions increased root-derived DOC. Here, all the detected marker effects (on positions 95.9-100.9 cM on chromosome 1H and 58.6-61.8 cM and 62.5-62.7 cM on chromosome 3H) were in agreement with the phenotype data. This is because the wild genotype Caesarea 26-24 was associated with larger root-derived DOC concentration (at 39d) in comparison to the modern genotype Harrington. As such, these results suggest that the wild barley genotype Caesarea 26-24 carries alleles that influence root-derived DOC on those three chromosome regions. This assumption is in line with the gene encoding T6P found on chromosome 3H, which could be a candidate gene influencing root C deposition and, in turn, root-derived DOC and MBC concentrations in soil. The T6P signalling system is a major regulator of resource allocation (e.g. C allocation and utilization in plants) and has been implicated in several processes in crop plants including assimilate partitioning and source-sink relationships (Paul et al. 2018(Paul et al. , 2020. To our knowledge, this study is the first to demonstrate the localization of the barley genetic influence on soil MBC and DOC, and thus on soil C cycling. The two genome regions on chromosome 3H (discussed above, although this could be considered a single QTL) affected both root-derived DOC and root-derived MBC. This is consistent with the relationship between these two traits, as determined by Pearson correlation analysis, being significant. Pearson correlation also showed a significant but low, positive relationship between root biomass and rootderived MBC. It is known that soil microbes acquire C (including root-derived C) in the form of DOC in the soil solution (Smolander and Kitunen 2002;Montaño et al. 2007), and that quantity of root exudation in barley plants is closely correlated with root biomass production (Darwent et al. 2003). Thus, our results support the likely links between root-derived MBC and root biomass or root-derived DOC.
However, moderate and low correlations between root-derived DOC and MBC and root biomass and root-derived MBC, respectively, suggest that rootderived MBC size was not mainly affected by root biomass size or DOC concentration. In previous work (Mwafulirwa et al. 2016), we showed that the activity (i.e. functional diversity) of soil microbes was affected by barley variety, likely due to differences in rhizodeposit quality (i.e. chemical composition) among the varieties or genotypes, which included genotypes used in the present study. This may lead us to assume that the variation in root-derived MBC observed here among the barley genotypes was also likely due to differences in rhizodeposit chemical composition among the genotypes, and that the marker effects for root-derived MBC or DOC revealed in this work were also associated with rhizodeposit composition. There are reports for QTLs and candidate genes influencing rhizodeposit chemical composition in other crop plants (Yan et al. 2004;Hongni et al. 2011;Qiu et al. 2014;Ramongolalaina et al. 2018). For instance, Hongni et al. (2011) mapped QTLs for three root exudates related to phosphorus efficiency in maize, while Yan et al. (2004) identified QTLs influencing root exudation of organic acids in common bean.

Chromosome regions affecting root biomass
Seven significant marker effects, on chromosomes 2H, 3H, 5H and 7H, were identified for root biomass. Others have also reported the localization of QTLs influencing root biomass in barley on these chromosomes (reviewed by Sallam et al. 2019). For example, Naz et al. (2014) used the 1536-SNP barley BOPA1 set and detected a total of 13 QTLs for root dry weight and root volume using wild barley introgression lines (developed by crossing a German spring cultivar Scarlett and a wild accession ISR42-8), of which seven QTLs were located on chromosomes 2H, 3H, 5H and 7H. Another study using barley lines genotyped using the Illumina 1536-SNP array ) also identified seven QTLs for root dry weight, with four of the QTLs localized on chromosomes 2H, 3H, 5H and 7H. In particular, Arifuzzaman et al. (2014) detected a QTL on 2H spanning from 38.9 to 66.0 cM, corroborating our present study showing significant marker effects on the same chromosome region (38.1-53.8 cM in the present study). This correspondence with previous studies adds additional confidence in our identification of barley chromosome regions influencing MBC and DOC pools in soil.

Conclusions
Our study is the first to identify barley chromosome regions and potential candidate genes influencing MBC and DOC concentrations in soil, and thus microbial populations and C cycling in soil. In total, 16 chromosome regions were identified for rootderived MBC and DOC, mainly on chromosomes 2H, 3H and 7H. Two locations on 3H influenced both traits (i.e. root-derived MBC and DOC). In addition, one chromosome region on 4H with the wild barley genome decreased root-derived MBC, and three on chromosomes 1H and 3H increased root-derived DOC, corroborating the observed phenotypes. This work also supports the previously detected QTLs for root biomass in barley, and therefore validating the chromosome regions influencing root-derived MBC and DOC. Thus, these findings are an important step towards better understanding of the plant genetic influences on soil microbial communities, which could be vital to control C and N cycling in soil through plant breeding approaches. Further work, for example using more replications, is needed to verify the barley chromosome regions or marker effects for MBC and DOC and potential candidate genes revealed in this study, and to sequence the soil microbes to examine the interaction between plant genotype and the diversity or structure of rhizosphere microbiomes. otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/. phosphorus uptake in common bean. Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.