Articles | Volume 8, issue 1
SOIL, 8, 49–57, 2022
SOIL, 8, 49–57, 2022

Short communication 20 Jan 2022

Short communication | 20 Jan 2022

Rhizodeposition efficiency of pearl millet genotypes assessed on a short growing period by carbon isotopes (δ13C and F14C)

Rhizodeposition efficiency of pearl millet genotypes assessed on a short growing period by carbon isotopes (δ13C and F14C)
Papa Mamadou Sitor Ndour1,5, Christine Hatté2,4, Wafa Achouak3, Thierry Heulin3, and Laurent Cournac1 Papa Mamadou Sitor Ndour et al.
  • 1Eco&Sols, Université de Montpellier, IRD, CIRAD, INRAE, Institut Agro, Montpellier, France
  • 2LSCE, CEA, CNRS, UVSQ, Université Paris-Saclay, 91191 Gif-sur-Yvette Cedex, France
  • 3LEMIRE BIAM, Aix Marseille University, CEA, CNRS, FR ECCOREV 3098, 13115 Saint-Paul-Lez-Durance, France
  • 4Institute of Physics – CSE, Silesian University of Technology, Gliwice, Poland
  • 5PMI Laboratory, Agrobiosciences, Mohammed VI Polytechnic University, Ben Guerir, Morocco

Correspondence: Papa Mamadou Sitor Ndour (


Rhizosheath size varies significantly with crop genotype, and root exudation is one among its driving factors. Unravelling the relationships between rhizosheath formation, root exudation and soil carbon dynamics may bring interesting perspectives in terms of crop breeding towards sustainable agriculture. Here we grew four pearl millet (C4 plant type: δ13C of 12.8 ‰, F14C = 1.012) inbred lines showing contrasting rhizosheath sizes in a C3 soil type (organic matter with δ13C of 22.3 ‰, F14C = 1.045). We sampled the root-adhering soil (RAS) and bulk soil after 28 d of growth under a semi-controlled condition. The soil organic carbon (SOC) content and δ13C and F14C of soil samples were measured and the plant-derived C amount and Clost/ Cnew ratio in the RAS were calculated. The results showed a significant increase in δ13C in the RAS of the four pearl millet lines compared to the control soil, suggesting that this approach was able to detect plant C input into the soil at an early stage of pearl millet growth. The concentration of plant-derived C in the RAS did not vary significantly between pearl millet lines, but the absolute amount of plant-derived C varied significantly when we considered the RAS mass of these different lines. Using a conceptual model and data from the two carbon isotopes' measurements, we evidenced a priming effect for all pearl millet lines. Importantly, the priming effect amplitude (Clost/ Cnew ratio) was higher for the small rhizosheath (low-aggregation) line than for the large rhizosheath (high-aggregation) ones, indicating a better C sequestration potential of the latter.

1 Introduction

In the context of climate change, a new challenge for agriculture is to sequester more carbon in the soil to mitigate CO2 increase in the atmosphere (Lal et al., 2015). This will be particularly important as increasing soil organic carbon (SOC) content would enhance soil fertility through improving physical and biological properties of the soil and then would have beneficial effects on agriculture and improve food security (Lal et al., 2015). Moreover, this strategy could be particularly relevant in the Sahel region of Africa, where very little above-ground cover remains after harvest, leading to soil carbon depletion (Baudron et al., 2014).

Pearl millet (Pennisetum glaucum (L.) R. Br.) is the main cereal grown in the Sahel and the arid region of India (Debieu et al., 2017). It is naturally adapted to drought conditions of these semiarid and arid regions and is a staple food for nearly 100 million people around the world (Burgarella et al., 2018; Varshney et al., 2017) and 50 million in western Africa. Following rainfall decline in recent decades, farmers have privileged short-cycle varieties that give optimal yield to cope with the short rainy season (Vigouroux et al., 2011). Despite this adaptation, pearl millet yields remain very low (less than 1 t ha−1) compared to other cereals. The main factors explaining this low productivity are (i) low soil fertility associated with low SOC content due in part to the export of crop residue and low carbon fertilization and (ii) climatic uncertainties (mainly rainfall) which could increase with climate change.

To simultaneously improve soil carbon sequestration through crop cultivation and ensure food security in Africa, a combination of plant breeding and the development of improved agricultural practices could be considered. Then, plant breeding for root traits (architecture and root exudation) could be an original strategy to improve water and mineral nutrition of crops and to enhance SOC sequestration. Root exudation is reported to account for 17 % to 40 % of plant photo-assimilates (Nguyen, 2003; Badri and Vivanco, 2009; Sasse et al., 2018). However, the question of whether cultivar variability can be used to significantly increase plant carbon deposition and its sequestration in soil is not fully documented. We have shown in a previous study that the mass of root-adhering soil (as an estimator of the rhizosheath) varied significantly with pearl millet genotype among a panel of inbred lines (Ndour et al., 2021). This emerging functional trait for plant phenotyping is linked to root exudation and other functional and morphological root traits such as root architecture (Delhaize et al., 2012, 2015; George et al., 2014) and to root-associated microbiota (Ndour et al., 2020). In order to evaluate how it may contribute to the carbon sequestration objective, the carbon input of pearl millet into the soil should be directly determined as well. However, due to the heterogeneous nature of the soil and particularly in-field conditions, conventional carbon measurement methods could not answer this question in short-term experiments such as those used for screening plant genotypes. The change in carbon concentration would remain below the variability around the mean value. Furthermore, it is now well known that the input of energetic molecules (such as root exudates) induces an increase in the activity of microorganisms and can thus contribute to extra mineralization of molecules derived from soil old carbon, that is, the so-called priming effect (Fontaine and Barot, 2005). Assessing the fate of the new carbon, i.e., assessing a balance between carbon gain and loss, cannot be achieved by looking only at carbon concentration on a point-by-point basis. The persistence of the carbon injected into the soil must be monitored over several years. As an alternative method, measuring carbon deposition in the rhizosphere using carbon isotopes (13C and 14C) that are much less sensitive to soil heterogeneity could be very interesting to test and indicative of the age of the primed carbon. For instance, based on the fact that C4 plant organic matter is enriched in 13C compared to C3 plant organic matter, He et al. (2019) used carbon isotopic shift (δ13C) in a wetland with initial C3 soil type colonized by a C4 plant (Spartina alterniflora) in order to determine SOC deposition and its carbon sequestration potential. Likewise, analysis of soil 14C profiles has revealed the ability of soil to store carbon and the impact of certain practices on both carbon stock and its persistence (Mathieu et al., 2015; Jreich et al., 2018).

The objective of this work was to apply these alternative methods to measure carbon deposition in the rhizosphere of pearl millet using carbon isotopes (13C and 14C) and to compare the ability of different pearl millet inbred lines to durably deposit carbon in the soil. We used four pearl millet lines with contrasting phenotypes for root-adhering soil aggregation (i.e. rhizosheath size) estimated by RAS / RT ratio (root-adhering soil mass / root tissue mass).

2 Material and methods

2.1 Pearl millet cultivation

Four pearl millet lines with contrasting RAS / RT ratios (i.e. rhizosheath size or rhizosphere aggregation) were selected from a previous study (Ndour et al., 2021): one low-aggregation line: L220, one intermediary-aggregation line: L3, and two high-aggregation lines: L253 and L132. The soil was sampled at Fissel Mbedap, Thiès, Senegal (142921.8 N 163536.4 W), in an uncultivated area essentially populated by Guiera senegalensis (C3 plant species). It is classified as Arenosol according to the FAO/WRB (IUSS Working Group WRB, 2015), and the texture characteristics were 93 % sand, 5.2 % silt and 2.6 % clay. The soil was sampled at the surface horizon (0–20 cm), sieved at 2 mm and homogenized. The δ13C of the soil was 22.3 ‰. The experiment was performed in greenhouse under natural light conditions in August 2018. The soil was filled in WM pots installed in plastic crates (Ndour et al., 2021). Each pot was filled with 1.5 kg of the soil moistened to its water-holding capacity (8 %), which was determined gravimetrically from the graph of water loss kinetics for 18 h. This watering level was adopted to have a good germination rate as the soil was sandy with a low water retention capacity, and we have used bottomless pots which allow water infiltration as in Ndour et al. (2017). The four pearl millet lines were sown in a randomized bloc design with seven replicates. Seven pots without plants were added to serve as controls. After 1 week, we thinned to one plant per pot. The seedlings were watered daily with 20 mL throughout the experiment until 1 d before the soil sampling.

2.2 Soil sampling

Harvesting was performed 28 d after sowing by detaching the two pieces of the pots and shaking gently to keep only the soil adhering closely to the roots (RAS). Roots and RAS were separated by washing the roots in 50 mL Falcon tubes containing 40 mL of distilled water. We also sampled bulk soil for each pot (away from the root system) as well as control soil from unplanted pots. Soil samples were dried by incubation at 65 C for 4 d. All soil samples were stored in haemolysis tubes until laboratory analyses. Shoot and root biomass were weighted, and shoot biomass was also sampled for δ13C measurement of pearl millet tissue.

2.3 Laboratory analyses

2.3.1 Organic carbon content

Approximately 15 to 20 mg of soil samples were weighed in tin cups for measurement (with a precision of 1 to 2 µg). The sample was combusted in a FlashEA1112 and the carbon content determined using the Eager software. One standard was inserted every 10 samples. For the range of carbon content of the soils analysed here, the relative measurement error is of the order of 1.5 %–2 %.

2.3.2δ13C analysis

Isotopic analysis was performed using continuous-flow EA-IRMS coupling, i.e. a FlashEA1112 Elemental Analyzer coupled to a Thermo Finnigan Delta + XP Isotope-Ratio Mass Spectrometer. Three in-house standards (Hobo5 sediment: δ13C =13.4 ‰, oxalic acid: δ13C =16.7 ‰ and GCL with δ13C =26.7 ‰) were inserted every five samples. Each in-house standard was regularly checked against international standards. The results are reported in the δ notation:


where Rsample and Rstandard are the 13C /12C ratios of the sample and the international standard, Vienna Pee Dee Belemnite (VPDB), respectively. The measurements were triplicated for representativeness. The external reproducibility of the analysis was better than 0.1 ‰, typically 0.06 ‰.

2.3.314C geochronology

The amount of soil needed to yield 1 mg of carbon was weighed in a tin capsule. Soil carbon was successively converted into CO2 and reduced in C in the presence of H2, using automated graphitization equipment, AGE3 (Wacker et al., 2010). The pure graphite was then pressed in the presence of ultrapure iron into a target to be introduced in the solid source of ECHoMICADAS (Synal et al., 2007; Tisnérat-Laborde et al., 2015), a compact radiocarbon system able to run very small samples. Results were reported as F14C as recommended by Reimer et al. (2004).

2.4 Calculations and statistical analyses

The proportion of plant-derived carbon (p) in the root-adhering soil was calculated according to the following Eq. (1):

(1) p = δ 13 C RAS - δ 13 C control δ 13 C mil - δ 13 C control ,

where δRAS represents the δ13C of the root-adhering soil organic carbon, δcontrol represents the δ13C of the unplanted soil organic carbon, and δmil represents the δ13C of pearl millet tissues we measured (12.8 ‰).

The plant-derived carbon content (PDCCRAS) of the root-adhering soil was calculated by multiplying the proportion of plant-derived carbon (p) by the SOC content:


The plant-derived carbon deposited in the root-adhering soil volume (PDCDRAS) was calculated by multiplying the PDCC by the average RAS of the different pearl millet lines in Table A1. This is expressed on a per-plant basis:

PDCDRAS(gC)=PDCCRAS(gCkg-1soil)×mass of RAS (kg).

The plant-derived C deposited / plant biomass was calculated as follows:

plant-derivedC/plant biomass=PDCDRAS/plant biomass.

Statistical analysis was performed using the R statistical environment (version 4.0.3). The normality of the data was tested using the Shapiro test (p<0.05). To test the effects of the pearl millet line on the different parameters (SOC, δ13C, F14C), general linear models (GLMs) were constructed using the “quick linear regression” (glm) function in R. Each model was fitted by considering the distribution mode and using the corresponding link function, i.e. Gaussian (link = “identity”) for normally distributed data and Poisson (link = “log”) for not normally distributed variables. Thereafter, analyses of variance (ANOVAs) were fitted to these models using the chi-squared (Chisq) test, and Tukey honest significance differences (Tukey HSDs) post hoc tests were performed (p<0.05) to compare the mean of the different parameters for the four pearl millet lines using the multicomp library available in R.

2.5 Processed data: integrated model

Three equations that derive from our analyses help in conceptualizing soil carbon dynamics:


with subscript “tot” for total organic carbon after the experiment, “old” for organic carbon before the experiment, “new” for the new carbon input during the experiment and “lost” for organic carbon that was present before the experiment and was lost during the experiment. “C” is for the carbon concentration.

Three hypotheses were investigated.

  • H1

    Input of new carbon without change in original carbon cortege Clost=0

  • H2

    Steady state: replacement of old carbon with new carbon Clost= Cnew

  • H3

    Negative net soil C balance: replacement and loss of old carbon Clost> Cnew

The isotopic information for “old” carbon was retrieved from the control soil. The “new” carbon input was derived from millet rhizodeposition. As a first approximation, the new C in the RAS fractions (δ13Cnew) was considered equal to millet leaf δ13C and δ13Clost was considered equal to δ13Cold. Solving the equation provided us with information on newly inputted carbon (Cnew) and released carbon (Clost, F14Clost). The derived parameters have been gathered in Table 2 in the Results section.

3 Results

3.1 SOC, δ13C and F14C of soil samples

SOC content was not significantly different between RAS, bulk soil and unplanted control fractions in any of the pearl millet lines, which was expected given the method of measuring C (combustion) and the short duration of the experiment (28 d) (Fig. 1a). Using natural 13C abundances (δ13C), significant increases were noticed in the root-adhering soil fraction (δ13CRAS) of all four pearl millet lines compared to the unplanted control soil (p<0.05), in agreement with the higher δ13C value of pearl millet leaves and probably also root exudates (12.8 ‰) compared to that of soil with C3 plant history (22.3 ‰). On the contrary, there is no significant difference in δ13C between all four bulk soils and the control soil, showing that there is no significant input of fresh C4-derived carbon outside the rhizosheath. For all four pearl millet lines, the δ13C of the RAS fractions was always higher than that of the bulk soil, with the difference being significant for only two lines (L3 and L253). No significant difference was found in the δ13C measured in the RAS between the four pearl millet lines (Fig. 1b).

Figure 1(a) Soil organic carbon (SOC) content of the RAS (root-adhering soil), bulk soil and unplanted control soil. (b) δ13C ‰ in the RAS, bulk soil and unplanted control soil. (c) F14C in the RAS, bulk soil and unplanted control soil. The dotted line indicates the F14C value of pearl millet tissue, i.e. the new carbon F14C (F14C =1.012±0.002). Different letters indicate significant difference using a Tukey post hoc test (p<0.05).


The plant-derived carbon content of the RAS varied from 0.124 g C kg−1 soil (L220) to 0.192 g C kg−1 soil (L3), but no significant difference was found between the four pearl millet lines (Fig. 2a). Considering the mass of root-adhering soil of the four pearl millet lines (Table S1), the amount of plant-derived carbon deposited in the RAS varied significantly according to pearl millet line and was significantly higher for the high-aggregation line (L253) compared to low-aggregation (L220) and intermediary-aggregation (L3) lines (Fig. 2b). The same trend was found for carbon deposited in the root-adhering soil expressed by plant total biomass (Fig. 2c).

The F14C values of the bulk soil of the four pearl millet lines were not significantly different from those of the unplanted control soil (Fig. 1c), showing the absence of 14C dilution away from the rhizosphere. The low-aggregation line (L220) and intermediary-aggregation line (L3) significantly diluted the 14C abundance in the RAS fraction compared to the bulk soil, in contrast to the two high-aggregation lines (L132 and L253) characterized by F14C values in RAS non-significantly different from those of the bulk soil (Fig. 1c).

The F14C values of bulk soil and unplanted soil control were significantly higher (F14C > 1.04) than those of pearl millet leaves (F14C =1.012±0.002, dotted line in Fig. 1c), illustrating a significant incorporation of carbon contemporaneous with the bomb peak into the “old” soil carbon pool. The dilution of the F14C in root-adhering soil fractions, due to the root exudation of C compounds with the same F14C value as leaves, was significantly different for the two high-aggregation lines (L253 and L132) compared to the two other lines (L220 and L3) (Fig. 1c).

Figure 2(a) Plant-derived carbon content (PDCC; g C kg−1 soil) in the root-adhering soil (RAS) of the four pearl millet lines. (b) Plant-derived carbon mass deposited (PDCD; g C) in the RAS of the four pearl millet lines. (c) Ratio between the plant-derived carbon mass and the total plant biomass produced by the four pearl millet lines. Black points within the boxes represent the means. Different letters indicate significant differences using an ANOVA and the Tukey post hoc test (p<0.05).


3.2 Conceptual modelling

Solving Eqs. (2), (3), and (4) did not provide a solution under either hypothesis 1 (input of new carbon without change in original carbon content) or hypothesis 2 (steady state). Only hypothesis 3 (priming) made it possible to solve the equations and provided consistent parameters (Table 1). The validation of this hypothesis showed that the experiment (4 weeks of pearl millet growth) is in a priming effect phase. For all pearl millet lines, new carbon has been deposited and some of the “old” carbon released. At this stage of plant growth, the trends showed a higher carbon loss vs. gain (high Clost/ Cnew ratio) for the low-aggregation and intermediary-aggregation lines L220 and L3 (4.4 and 3.5) compared to the two high-aggregation lines L132 and L253 (2.9 and 3.3) (Table 1). Investigation with the 14C isotope allowed the lost carbon to be specified. The F14C value of lost carbon was between 1.2 and 1.3, i.e. carbon which is consistent with the 1980s (Hua et al., 2013). This means that the injection of new carbon also implies loss of some decades-old carbon.

Table 1Data processing for the root-adhering soil (RAS) of pearl millet lines. Hypotheses 1 and 2 are not solvable with the measured parameters. Hypothesis 3 is the only one that has a solution. The amount of carbon added (Cnew), the amount of carbon lost (Clost) and the radiocarbon fraction of this lost carbon (F14Clost) are then evaluated. Carbon content is expressed in % wt, δ13C as ‰, and F14C is unitless. The values are given with a digit number in accordance with the number of significant figures.

Download Print Version | Download XLSX

4 Discussion

Our data showed no statistically significant changes in soil carbon content between the RAS of the different pearl millet lines, neither with the bulk soil nor with the control soil, after 28 d of growth. Similar results were obtained in other studies in which no significant differences were noticed in the rhizosphere SOC content compared to the control soil (Wang et al., 2016; Van de Broek et al., 2020). In contrast, using natural 13C abundances, we demonstrated a significant increase in the δ13C of the root-adhering soil fraction of the four pearl millet lines compared to their respective bulk soil fractions and to the unplanted control soil (p<0.05). This shift in δ13C is in agreement with the higher δ13C value of pearl millet leaves (12.8 ‰) and probably also of root exudates compared to that of soil with a C3 plant history (22.3 ‰), evidencing the carbon deposition into the soil through root exudation. Nevertheless, plant-derived C content (g C kg−1 soil) in the rhizosphere (RAS) did not vary significantly between pearl millet lines. The trend for increasing aggregation efficiency between pearl millet lines was not found in the δ13C-derived evaluation of input carbon content. This result is consistent with those obtained by Van de Broek et al. (2020), who found no significant differences in net carbon rhizodeposition between four wheat cultivars with contrasting root biomass (including two old and two recent varieties). However, the absolute amount of plant-derived C in RAS (obtained by multiplying plant-derived C content by the mean RAS mass of the different pearl millet lines) was significantly higher for the high-aggregation pearl millet line (L253) compared to the low- and intermediary-aggregation lines (L220 and L3, respectively). The ratio between this plant-derived C deposited in the RAS and total plant biomass (expressed in percent) was also higher for the high-aggregation pearl millet line (L253) compared to the L220 and L3 lines, allowing for the detection of about 0.1 % of plant-derived C per unit of plant biomass. This indicates that the significant difference in root exudation among these pearl millet lines would be related to a difference in carbon allocation rather than a difference in biomass quantity production.

The distinct evolution between the F14C contents of the RAS of the different pearl millet lines and the bulk soil clearly showed a different distribution of soil carbon age between the two fractions. The decrease in F14C in the RAS fraction of the L220 and L3 lines may result from an increase in the proportion of “new” carbon from root exudation with F14C equal to that of the millet leaf (see Fig. 12 in Jreich et al., 2018) compared to that of the bulk soil (equivalent to the control soil). This would be a dilution effect by incorporation of “new” carbon. Here we are probably faced with contrasting situations: a loss of “old” carbon due to the priming effect of root exudates with a partial replacement of the lost “old” carbon by “new carbon” from root exudates and probably their associated microbial metabolites, this replacement tending to be of lower amplitude in low-aggregation lines (case of L220 and L3) than in high-aggregation lines (case of L253 and L132). The 14C isotope analyses allow us to estimate both the ratio of carbon input to output and the mean age of the lost carbon, which is a few decades (1980s).

Mwafulirwa et al. (2016) showed a genotypic influence of barley on the stabilization of rhizodeposited C and the SOM mineralization. Using a combination of natural isotope tracing and conceptual modelling, we confirmed these findings in pearl millet, as our results show that lines L253 and L132 better preserve the soil-existing C stock (“old” C) after plant deposition of “new” C than lines L220 and L3. In addition, better C preservation (lower priming) was importantly noticed with the high-aggregation genotypes (L132 and L253). This suggests that the better soil C preservation in the root vicinity of these pearl millet lines could be related to its transformation into microbial macromolecules such as exopolysaccharides which make C more stable and contribute to increasing rhizosphere soil aggregation, as reported in previous studies (Gouzou et al., 1993; Alami et al., 2000; Bezzate et al., 2000). Root exudates include low molecular weight molecules (sugars, amino acids, organic acids, phenolics) and high molecular weight molecules (proteins and mucilage) (Bais et al., 2006), which are probably subject to differential degradation capacity by soil microbiota. Therefore, the variation of soil C stabilization in these four pearl millet lines could also be related to a genotypic variation in the root exudate quality (biochemical composition), as reported in different plant species (Liu et al., 2019; Semchenko et al., 2021).

5 Conclusions

These results show that natural carbon isotope tracing allowed us to detect pearl millet root exudation into soil at an early stage of plant growth (δ13C at 28 d) and to study the dynamics of “old” and “new” C, demonstrating the impact of the rhizosphere priming effect at this stage of plant growth (F14C). With the use of conceptual modelling applied to data from these two isotopes, we pinpoint a way of assessing how the root C metabolism of some plant lines/varieties can contribute effectively to the storage of carbon in the soils with an enrichment in – or a better preservation of – SOM and therefore contribute to a reduction of greenhouse gases' balance in agriculture. Further studies should focus on studying the root exudate profile among these pearl millet genotypes to see how it might help explain the observed variability in the dynamic of plant-derived C in the rhizosphere. Nevertheless, whatever the relative contribution of the quantity and the quality of root exudation in shaping the rhizosheath size, the latter is now recognized as a relevant root trait, and genetic studies are being carried out to detect QTLs and to identify their controlling genes in pearl millet, and this would provide interesting tools to breeders for the selection of efficient genotypes for a sustainable production.

Appendix A

Table A1RAS / RT means (rhizosheath size) of the four pearl millet lines measured in Ndour et al. (2021), root (RT) mass, shoot mass, total plant biomass measured for the four pearl millet lines, and the mean values of the root-adhering soil (RAS) mass calculated for the four pearl millet lines. All these parameters are measured/calculated at 28 d of plant growth.

Download Print Version | Download XLSX

Data availability

The data generated in this study are available from the corresponding authors upon reasonable request.


The supplement related to this article is available online at:

Author contributions

PMSN, CH, WA, TH and LC conceptualized the study. PMSN performed the experiment and the soil sampling. CH performed laboratory analyses and the conceptual modelling. PMSN wrote the initial draft. PMSN, CH, WA, TH and LC contributed to generating and reviewing the subsequent versions of the manuscript.

Competing interests

The contact author has declared that neither they nor their co-authors have any competing interests.


Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


Papa Mamadou Sitor Ndour would like to thank the French National Research Institute for Sustainable Development (IRD) and the Make Our Planet Great Again (MOPGA) initiative of the French government for co-funding his postdoctoral fellowship. We thank Oscar Pascal Malou for communicating data on soil texture.

Financial support

This research has been supported by the Agence Nationale de la Recherche (grant no. ANR-17-CE20-0022) and the Commissariat à l'Énergie Atomique et aux Énergies Alternatives (CarbonSeq, DRF impulsion).

Review statement

This paper was edited by Paul Hallett and reviewed by Paul Hallett and one anonymous referee.


Alami, Y., Achouak, W., Marol, C., and Heulin, T.: Rhizosphere soil aggregation and plant growth promotion of sunflowers by an exopolysaccharide-producing Rhizobium sp. strain isolated from sunflower roots, Appl. Environ. Microbiol., 66, 3393–3398,, 2000. 

Badri, D. V. and Vivanco, J. M.: Regulation and function of root exudates, Plant Cell Environ., 32, 666–681,, 2009. 

Bais, H. P., Weir, T. L., Perry, L. G., Gilroy, S., and Vivanco, J. M.: The role of root exudates in rhizosphere interactions with plants and other organisms, Ann. Rev. Plant Biol., 57, 233–266,, 2006. 

Baudron, F., Jaleta, M., Okitoi, O., and Tegegn, A.: Conservation agriculture in African mixed crop-livestock systems: Expanding the niche, Agr. Ecosyst. Environ., 187, 171–182,, 2014. 

Bezzate, S., Aymerich, S., Chambert, R., Czarnes, S., Berge, O., and Heulin, T.: Disruption of the Paenibacillus polymyxa levansucrase gene impairs its ability to aggregate soil in the wheat rhizosphere, Environ. Microbiol., 2, 333–342,, 2000. 

Burgarella, C., Cubry, P., Kane, N. A., Varshney, R. K., Mariac, C., Liu, X., Shi, C., Thudi, M., Couderc, M., Xu, X., Chitikineni, A., Scarcelli, N., Barnaud, A., Rhoné, B., Dupuy, C., François, O., Berthouly-Salazar, C., and Vigouroux, Y.: A western Sahara centre of domestication inferred from pearl millet genomes, Nat. Ecol. Evol., 2, 1377–1380,, 2018. 

Debieu, M., Kanfany, G., and Laplaze, L.: Pearl millet genome: lessons from a tough crop, Trends Plant Sci., 22, 911–913,, 2017. 

Delhaize, E., James, R. A., and Ryan, P. R.: Aluminium tolerance of root hairs underlies genotypic differences in rhizosheath size of wheat (Triticum aestivum) grown on acid soil, New Phytol., 195, 609–619,, 2012. 

Delhaize, E., Rathjen, T. M., and Cavanagh, C. R.: The genetics of rhizosheath size in a multiparent mapping population of wheat, J. Exp. Bot., 66, 4527–4536,, 2015. 

Fontaine, S. and Barot, S.: Size and functional diversity of microbe populations control plant persistence and long-term soil carbon accumulation, Ecol. Lett., 8, 1075–1087,, 2005. 

George, T. S., Brown, L. K., Ramsay, L., White, P. J., Newton, A. C., Bengough, A. G., Russell, J., and Thomas, W. T. B.: Understanding the genetic control and physiological traits associated with rhizosheath production by barley (Hordeum vulgare), New Phytol., 203, 195–205,, 2014. 

Gouzou, L., Burtin, G., Philippy, R., Bartoli, F., and Heulin, T.: Effect of inoculation with Bacillus polymyxa on soil aggregation in the wheat rhizosphere: preliminary examination, Geoderma, 56, 479–491,, 1993. 

He, Y., Zhou, X., Cheng, W., Zhou, L., Zhang, G., Zhou, G., Liu, R., Shao, J., Zhu, K., and Cheng, W.: Linking Improvement of Soil Structure to Soil Carbon Storage Following Invasion by a C4 Plant Spartina alterniflora, Ecosystems, 22, 859–872,, 2019. 

Hua, Q., Barbetti, M., and Rakowski, A. Z.: Atmospheric Radiocarbon for the Period 1950–2010, Radiocarbon, 55, 2059–2072,, 2013. 

IUSS Working Group WRB: World Reference Base for Soil Resources 2014, Update 2015 International Soil Classification System for Naming Soils and Creating Legends for Soil Maps, FAO, Rome, ISBN 978-92-5-108369-7, 2015. 

Jreich, R., Hatte, C., Balesdent, J., and Parent, É.: Bayesian selection of mixed covariates from a latent layer: application to hierarchical modeling of soil carbon dynamics, J. Soc. Fr. Stat., 159, 128–155, 2018. 

Lal, R., Negassa, W., and Lorenz, K.: Carbon sequestration in soil, Curr. Opin. Environ. Sustain., 15, 79–86,, 2015. 

Liu, T.-Y., Chen, M.-X., Zhang, Y., Zhu, F.-Y., Liu, Y.-G., Tian, Y., Fernie, A. R., Ye, N., and Zhang, J.: Comparative metabolite profiling of two switchgrass ecotypes reveals differences in drought stress responses and rhizosheath weight, Planta, 250, 1355–1369,, 2019. 

Mathieu, J. A., Hatté, C., Balesdent, J., and Parent, É.: Deep soil carbon dynamics are driven more by soil type than by climate: a worldwide meta-analysis of radiocarbon profiles, Glob. Change Biol., 21, 4278–4292,, 2015. 

Mwafulirwa, L., Baggs, E. M., Russell, J., George, T., Morley, N., Sim, A., de la Fuente Cantó, C., and Paterson, E.: Barley genotype influences stabilization of rhizodeposition-derived C and soil organic matter mineralization, Soil Biol. Biochem., 95, 60–69,, 2016. 

Ndour, P. M. S., Heulin, T., Achouak, W., Laplaze, L., and Cournac, L.: The rhizosheath: from desert plants adaptation to crop breeding, Plant Soil, 456, 1–13,, 2020. 

Ndour, P. M. S., Barry, C. M., Tine, D., De la Fuente Cantó, C., Gueye, M., Barakat, M., Ortet, P., Achouak, W., Ndoye, I., Sine, B., Laplaze, L., Heulin, T., and Cournac, L.: Pearl millet genotype impacts microbial diversity and enzymatic activities in relation to root-adhering soil aggregation, Plant Soil, 464, 109–129,, 2021. 

Nguyen, C.: Rhizodeposition of organic C by plants: mechanisms and controls, Agron. Sustain. Dev., 23, 375–396,, 2003. 

Reimer, P. J., Brown, T. A., and Reimer, R. W.: Discussion: Reporting and Calibration of Post-Bomb 14C Data, Radiocarbon, 46, 1299–1304,, 2004. 

Sasse, J., Martinoia, E., and Northen, T.: Feed Your Friends: Do Plant Exudates Shape the Root Microbiome?, Trends Plant Sci., 23, 25–41,, 2018. 

Semchenko, M., Xue, P., and Leigh, T.: Functional diversity and identity of plant genotypes regulate rhizodeposition and soil microbial activity, New Phytol., 232, 776–787,, 2021. 

Synal, H.-A., Stocker, M., and Suter, M.: MICADAS: A new compact radiocarbon AMS system, Nucl. Instrum. Meth. B, 259, 7–13,, 2007. 

Tisnérat-Laborde, N., Thil, F., Synal, H. A., Cersoy, S., Hatté, C., Gauthier, C., Massault, M., Michelot, J. L., Noret, A., Siani, G., Tombret, O., Vigne, J. D., and Zazzo, D.: ECHoMICADAS: A new compact AMS system to measuring 14C for Environment, 22nd International Radiocarbon Conference, Dakar, 16–20, 2015. 

Van de Broek, M., Ghiasi, S., Decock, C., Hund, A., Abiven, S., Friedli, C., Werner, R. A., and Six, J.: The soil organic carbon stabilization potential of old and new wheat cultivars: a 13CO2-labeling study, Biogeosciences, 17, 2971–2986,, 2020. 

Varshney, R. K., Shi, C., Thudi, M., Mariac, C., Wallace, J., Qi, P., Zhang, H., Zhao, Y., Wang, X., Rathore, A., Srivastava, R. K., Chitikineni, A., Fan, G., Bajaj, P., Punnuri, S., Gupta, S. K., Wang, H., Jiang, Y., Couderc, M., Katta, M. A. V. S. K., Paudel, D. R., Mungra, K. D., Chen, W., Harris-Shultz, K. R., Garg, V., Desai, N., Doddamani, D., Kane, N. A., Conner, J. A., Ghatak, A., Chaturvedi, P., Subramaniam, S., Yadav, O. P., Berthouly-Salazar, C., Hamidou, F., Wang, J., Liang, X., Clotault, J., Upadhyaya, H. D., Cubry, P., Rhoné, B., Gueye, M. C., Sunkar, R., Dupuy, C., Sparvoli, F., Cheng, S., Mahala, R. S., Singh, B., Yadav, R. S., Lyons, E., Datta, S. K., Hash, C. T., Devos, K. M., Buckler, E., Bennetzen, J. L., Paterson, A. H., Ozias-Akins, P., Grando, S., Wang, J., Mohapatra, T., Weckwerth, W., Reif, J. C., Liu, X., Vigouroux, Y., and Xu, X.: Pearl millet genome sequence provides a resource to improve agronomic traits in arid environments, Nat. Biotechnol., 35, 969–976,, 2017. 

Vigouroux, Y., Mariac, C., De Mita, S., Pham, J.-L., Gérard, B., Kapran, I., Sagnard, F., Deu, M., Chantereau, J., Ali, A., Ndjeunga, J., Luong, V., Thuillet, A.-C., Saïdou, A.-A., and Bezançon, G.: Selection for earlier flowering crop associated with climatic variations in the Sahel, PLoS One, 6, e19563,, 2011. 

Wacker, L., Němec, M., and Bourquin, J.: A revolutionary graphitisation system: Fully automated, compact and simple, Nucl. Instrum. Meth. B, 268, 931–934,, 2010. 

Wang, X., Tang, C., Severi, J., Butterly, C. R., and Baldock, J. A.: Rhizosphere priming effect on soil organic carbon decomposition under plant species differing in soil acidification and root exudation, New Phytol., 211, 864–873,, 2016. 

Short summary
Unravelling relationships between plant rhizosheath, root exudation and soil C dynamic may bring interesting perspectives in breeding for sustainable agriculture. Using four pearl millet lines with contrasting rhizosheaths, we found that δ13C and F14C of root-adhering soil differed from those of bulk and control soil, indicating C exudation in the rhizosphere. This C exudation varied according to the genotype, and conceptual modelling performed with data showed a genotypic effect on the RPE.