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

Papa Mamadou Sitor Ndour, Christine Hatté, Wafa Achouak, Thierry Heulin, Laurent Cournac

► To cite this version:

Papa Mamadou Sitor Ndour, Christine Hatté, Wafa Achouak, Thierry Heulin, Laurent Cournac. Rhizodeposition efficiency of pearl millet genotypes assessed on a short growing period by carbon isotopes (δ13C and F14C). Soil, 2022, 8 (1), pp.49 - 57. 10.5194/soil-8-49-2022. hal-03537131

HAL Id: hal-03537131

https://hal.science/hal-03537131

Submitted on 20 Jan 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished 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_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Rhizodeposition efficiency of pearl millet genotypes assessed on a short growing period by carbon isotopes (δ¹³C and F¹⁴C)

Papa Mamadou Sitor Ndour¹,⁵, Christine Hatté²,⁴, Wafa Achouak³, Thierry Heulin³, and Laurent Cournac¹

¹Eco&Sols, Université de Montpellier, IRD, CIRAD, INRAE, Institut Agro, Montpellier, France

²LSCE, CEA, CNRS, UVSQ, Université Paris-Saclay, 91191 Gif-sur-Yvette Cedex, France

³LEMIRE BIAM, Aix Marseille University, CEA, CNRS, FR ECCOREV 3098,
13115 Saint-Paul-Lez-Durance, France

⁴Institute of Physics – CSE, Silesian University of Technology, Gliwice, Poland

⁵PMI Laboratory, Agrobiosciences, Mohammed VI Polytechnic University,
Ben Guerir, Morocco

**Correspondence:** Papa Mamadou Sitor Ndour (sitndour@yahoo.fr)

Received: 19 September 2021 – Discussion started: 28 September 2021

Revised: 10 December 2021 – Accepted: 14 December 2021 – Published: 20 January 2022

**Abstract.** 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 (C₄ plant type: δ¹³C of -12.8‰, F¹⁴C = 1.012) inbred lines showing contrasting rhizosheath sizes in a C₃ soil type (organic matter with δ¹³C of -22.3‰, F¹⁴C = 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 δ¹³C and F¹⁴C 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 δ¹³C 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 CO₂ 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 \text{ 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 ($^{13}$C and $^{14}$C) 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 $C_4$ plant organic matter is enriched in $^{13}$C compared to $C_3$ plant organic matter, He et al. (2019) used carbon isotopic

shift ($\delta^{13}$C) in a wetland with initial $C_3$ soil type colonized by a $C_4$ plant (*Spartina alterniflora*) in order to determine SOC deposition and its carbon sequestration potential. Likewise, analysis of soil $^{14}$C 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 ($^{13}$C and $^{14}$C) 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 (14°29'21.8 N 16°35'36.4 W), in an uncultivated area essentially populated by *Guiera senegalensis* ($C_3$ 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 $\delta^{13}$C 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 δ¹³C 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 δ¹³C 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: δ¹³C = -13.4‰, oxalic acid: δ¹³C = -16.7‰ and GCL with δ¹³C = -26.7‰) were inserted every five samples. Each in-house standard was regularly checked against international standards. The results are reported in the δ notation:

$$ \delta^{13}C = (R_{\text{sample}}/R_{\text{standard}} - 1), $$

where Rsample and Rstandard are the ¹³C/¹²C 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.3 ¹⁴C geochronology

The amount of soil needed to yield 1 mg of carbon was weighed in a tin capsule. Soil carbon was successively converted into CO₂ and reduced in C in the presence of H₂, 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 F¹⁴C 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):

$$ p = \frac{\delta^{13}C_{\text{RAS}} - \delta^{13}C_{\text{control}}}{\delta^{13}C_{\text{mil}} - \delta^{13}C_{\text{control}}}, \quad (1) $$

where δRAS represents the δ¹³C of the root-adhering soil organic carbon, δcontrol represents the δ¹³C of the unplanted soil organic carbon, and δmil represents the δ¹³C 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:

$$ \text{PDCC}_{\text{RAS}} (\text{g C kg}^{-1} \text{ soil}) = p \times \text{SOC} (\text{g C kg}^{-1} \text{ soil}). $$

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:

$$ \text{PDCD}_{\text{RAS}} (\text{g C}) = \text{PDCC}_{\text{RAS}} (\text{g C kg}^{-1} \text{ soil}) \times \text{mass of RAS (kg)}. $$

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

$$ \text{plant-derived C / plant biomass} = \text{PDCD}_{\text{RAS}} / \text{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, δ¹³C, F¹⁴C), 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:

$$ C_{\text{tot}} = C_{\text{old}} + C_{\text{new}} - C_{\text{lost}}, \quad (2) $$

$$ C_{\text{tot}} \cdot \delta^{13}C_{\text{tot}} = C_{\text{old}} \cdot \delta^{13}C_{\text{old}} + C_{\text{new}} \cdot \delta^{13}C_{\text{new}} - C_{\text{lost}} \cdot \delta^{13}C_{\text{lost}}, \quad (3) $$

$$ C_{\text{tot}} \cdot F^{14}C_{\text{tot}} = C_{\text{old}} \cdot F^{14}C_{\text{old}} + C_{\text{new}} \cdot F^{14}C_{\text{new}} - C_{\text{lost}} \cdot F^{14}C_{\text{lost}}, \quad (4) $$

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 $C_{lost} = 0$

H2 Steady state: replacement of old carbon with new carbon $C_{lost} = C_{new}$

H3 Negative net soil C balance: replacement and loss of old carbon $C_{lost} > C_{new}$

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 ($\delta^{13}C_{new}$) was considered equal to millet leaf $\delta^{13}C$ and $\delta^{13}C_{lost}$ was considered equal to $\delta^{13}C_{old}$. Solving the equation provided us with information on newly inputted carbon ($C_{new}$) and released carbon ($C_{lost}$, $F^{14}C_{lost}$). The derived parameters have been gathered in Table 2 in the Results section.

# 3 Results

## 3.1 SOC, $\delta^{13}C$ and $F^{14}C$ 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 $^{13}C$ abundances ($\delta^{13}C$), significant increases were noticed in the root-adhering soil fraction ($\delta^{13}C_{RAS}$) of all four pearl millet lines compared to the unplanted control soil ($p<0.05$), in agreement with the higher $\delta^{13}C$ value of pearl millet leaves and probably also root exudates ($-12.8\%e$) compared to that of soil with $C_3$ plant history ($-22.3\%e$). On the contrary, there is no significant difference in $\delta^{13}C$ between all four bulk soils and the control soil, showing that there is no significant input of fresh $C_4$-derived carbon outside the rhizosheath. For all four pearl millet lines, the $\delta^{13}C$ 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 $\delta^{13}C$ measured in the RAS between the four pearl millet lines (Fig. 1b).

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)

<figure><img src="image_2.png" /><figcaption>**Figure 1.** (a) Soil organic carbon (SOC) content of the RAS (root-adhering soil), bulk soil and unplanted control soil. (b) $\delta^{13}C\%$ in the RAS, bulk soil and unplanted control soil. (c) $F^{14}C$ in the RAS, bulk soil and unplanted control soil. The dotted line indicates the $F^{14}C$ value of pearl millet tissue, i.e. the new carbon $F^{14}C$ ($F^{14}C = 1.012 \pm 0.002$). Different letters indicate significant difference using a Tukey post hoc test ($p<0.05$).</figcaption></figure>
lines (Fig. 2b). The same trend was found for carbon deposited in the root-adhering soil expressed by plant total biomass (Fig. 2c).

The F¹⁴C 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 ¹⁴C dilution away from the rhizosphere. The low-aggregation line (L220) and intermediary-aggregation line (L3) significantly diluted the ¹⁴C abundance in the RAS fraction compared to the bulk soil, in contrast to the two high-aggregation lines (L132 and L253) characterized by F¹⁴C values in RAS non-significantly different from those of the bulk soil (Fig. 1c).

The F¹⁴C values of bulk soil and unplanted soil control were significantly higher (F¹⁴C > 1.04) than those of pearl millet leaves (F¹⁴C = 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 F¹⁴C in root-adhering soil fractions, due to the root exudation of C compounds with the same F¹⁴C 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).

## 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 C<sub>lost</sub> / C<sub>new</sub> 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 ¹⁴C isotope allowed the lost carbon to be specified. The F¹⁴C 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.

## 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 ¹³C abundances, we demonstrated a significant increase in the δ¹³C of the root-adhering soil fraction of the four pearl

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

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

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

**Figure 2.** (a) Plant-derived carbon content (PDCC; g C kg⁻¹ 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).
<table><caption>**Table 1.** Data 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 ($C_{new}$), the amount of carbon lost ($C_{lost}$) and the radiocarbon fraction of this lost carbon ($F^{14}C_{lost}$) are then evaluated. Carbon content is expressed in % wt, $\delta^{13}C$ as ‰, and $F^{14}C$ is unitless. The values are given with a digit number in accordance with the number of significant figures.</caption><thead><tr><th></th><th></th><th>L220</th><th>L3</th><th>L132</th><th>L253</th></tr></thead><tbody><tr><td colspan="6">Hypothesis 1: C<sub>lost</sub> = 0</td></tr><tr><td>RAS</td><td>C<sub>lost</sub></td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>Control</td><td>C<sub>old</sub></td><td>0.43</td><td>0.43</td><td>0.43</td><td>0.43</td></tr><tr><td>C<sub>tot</sub>–C<sub>old</sub></td><td>C<sub>new</sub></td><td>No solution</td><td>No solution</td><td>No solution</td><td>No solution</td></tr><tr><td colspan="6">Hypothesis 2: C<sub>lost</sub> = C<sub>new</sub></td></tr><tr><td>Control</td><td>C<sub>old</sub></td><td>0.43</td><td>0.43</td><td>0.43</td><td>0.43</td></tr><tr><td>RAS</td><td>C<sub>tot</sub></td><td>0.39</td><td>0.38</td><td>0.40</td><td>0.40</td></tr><tr><td>C<sub>tot</sub> = C<sub>old</sub></td><td></td><td>Not fulfilled</td><td>Not fulfilled</td><td>Not fulfilled</td><td>Not fulfilled</td></tr><tr><td colspan="6">Hypothesis 3: C<sub>lost</sub>–C<sub>new</sub> > 0</td></tr><tr><td>Control</td><td>C<sub>old</sub></td><td>0.43</td><td>0.43</td><td>0.43</td><td>0.43</td></tr><tr><td>RAS</td><td>C<sub>tot</sub></td><td>0.39</td><td>0.38</td><td>0.40</td><td>0.40</td></tr><tr><td></td><td>δ<sup>13</sup>C<sub>old</sub></td><td>-22.3</td><td>-22.3</td><td>-22.3</td><td>-22.3</td></tr><tr><td></td><td>δ<sup>13</sup>C<sub>tot</sub></td><td>-22.1</td><td>-21.8</td><td>-22.0</td><td>-22.0</td></tr><tr><td></td><td>δ<sup>13</sup>C<sub>new</sub></td><td>-12.8</td><td>-12.8</td><td>-12.8</td><td>-12.8</td></tr><tr><td></td><td>δ<sup>13</sup>C<sub>lost</sub></td><td>-22.3</td><td>-22.3</td><td>-22.3</td><td>-22.3</td></tr><tr><td>(2) + (3) →</td><td>C<sub>lost</sub></td><td>0.040</td><td>0.066</td><td>0.038</td><td>0.039</td></tr><tr><td>(2) →</td><td>C<sub>new</sub></td><td>0.009</td><td>0.019</td><td>0.013</td><td>0.012</td></tr><tr><td></td><td>C<sub>lost</sub> / C<sub>new</sub></td><td>4.4</td><td>3.5</td><td>2.9</td><td>3.3</td></tr><tr><td></td><td>F<sup>14</sup>C<sub>old</sub></td><td>1.045</td><td>1.045</td><td>1.045</td><td>1.045</td></tr><tr><td></td><td>F<sup>14</sup>C<sub>tot</sub></td><td>1.017</td><td>1.014</td><td>1.024</td><td>1.028</td></tr><tr><td></td><td>F<sup>14</sup>C<sub>new</sub></td><td>1.012</td><td>1.012</td><td>1.012</td><td>1.012</td></tr><tr><td>(4) →</td><td>F<sup>14</sup>C<sub>lost</sub></td><td>1.311</td><td>1.212</td><td>1.257</td><td>1.213</td></tr></tbody></table>

millet lines compared to their respective bulk soil fractions and to the unplanted control soil ($p<0.05$). This shift in $\delta^{13}C$ is in agreement with the higher $\delta^{13}C$ value of pearl millet leaves ($-12.8\‰$) and probably also of root exudates compared to that of soil with a $C_3$ 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 $\delta^{13}C$-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 mil-

let 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 $F^{14}C$ 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 $F^{14}C$ 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 $F^{14}C$ 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 $^{14}$C 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 trac-
ing 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 no-
ticed 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 ex-
opolysaccharides 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, phe-
nolics) and high molecular weight molecules (proteins and
mucilage) (Bais et al., 2006), which are probably subject to
differential degradation capacity by soil microbiota. There-
fore, 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 re-
ported in different plant species (Liu et al., 2019; Semchenko
et al., 2021).

# 5 Conclusions

These results show that natural carbon isotope tracing al-
lowed us to detect pearl millet root exudation into soil at an
early stage of plant growth (δ¹³C 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
(F¹⁴C). 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 con-
tribute 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’ bal-
ance 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. Neverthe-
less, 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 pro-
vide interesting tools to breeders for the selection of efficient
genotypes for a sustainable production.
Appendix A

<table><caption>**Table A1.** RAS / 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.</caption><thead><tr><th>Lines</th><th>RAS / RT</th><th>RT (g)</th><th>Shoot (g)</th><th>Plant biomass (g)</th><th>RAS (g)</th></tr></thead><tbody><tr><td>L220</td><td>8.4</td><td>0.116 (± 0.07)</td><td>0.286 (± 0.08)</td><td>0.402 (± 0.15)</td><td>0.972 (± 0.6)</td></tr><tr><td>L3</td><td>14.9</td><td>0.096 (± 0.06)</td><td>0.194 (± 0.07)</td><td>0.290 (± 0.13)</td><td>1.433 (± 0.9)</td></tr><tr><td>L132</td><td>23.3</td><td>0.082 (± 0.04)</td><td>0.231 (± 0.03)</td><td>0.313 (± 0.06)</td><td>1.921 (± 0.8)</td></tr><tr><td>L253</td><td>20.8</td><td>0.128 (± 0.06)</td><td>0.244 (± 0.07)</td><td>0.372 (± 0.13)</td><td>2.659 (± 1.2)</td></tr></tbody></table>

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

**Supplement.** The supplement related to this article is available online at: https://doi.org/10.5194/soil-8-49-2022-supplement.

**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.

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

**Acknowledgements.** 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.

References

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, https://doi.org/10.1128/AEM.66.8.3393-3398.2000, 2000.

Badri, D. V. and Vivanco, J. M.: Regulation and function of root exudates, Plant Cell Environ., 32, 666–681, https://doi.org/10.1111/j.1365-3040.2009.01926.x, 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, https://doi.org/10.1146/annurev.arplant.57.032905.105159, 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, https://doi.org/10.1016/j.agee.2013.08.020, 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, https://doi.org/10.1046/j.1462-2920.2000.00114.x, 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, https://doi.org/10.1038/s41559-018-0643-y, 2018.

Debieu, M., Kanfany, G., and Laplaze, L.: Pearl millet genome: lessons from a tough crop, Trends Plant Sci., 22, 911–913, https://doi.org/10.1016/j.tplants.2017.09.006, 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, https://doi.org/10.1111/j.1469-8137.2012.04183.x, 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, https://doi.org/10.1093/jxb/erv223, 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, https://doi.org/10.1111/j.1461-0248.2005.00813.x, 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, https://doi.org/10.1111/nph.12786, 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, https://doi.org/10.1016/0016-7061(93)90128-8, 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 C₄ Plant *Spartina alterniflora*, Ecosystems, 22, 859–872, https://doi.org/10.1007/s10021-018-0308-3, 2019.

Hua, Q., Barbetti, M., and Rakowski, A. Z.: Atmospheric Radiocarbon for the Period 1950–2010, Radiocarbon, 55, 2059–2072, https://doi.org/10.2458/azu_js_rc.v55i2.16177, 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, https://doi.org/10.1016/j.cosust.2015.09.002, 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, https://doi.org/10.1007/s00425-019-03228-w, 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, https://doi.org/10.1111/gcb.13012, 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, https://doi.org/10.1016/j.soilbio.2015.12.011, 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, https://doi.org/10.1007/s11104-020-04700-3, 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, https://doi.org/10.1007/s11104-021-04917-w, 2021.

Nguyen, C.: Rhizodeposition of organic C by plants: mechanisms and controls, Agron. Sustain. Dev., 23, 375–396, https://doi.org/10.1051/agro:2003011, 2003.

Reimer, P. J., Brown, T. A., and Reimer, R. W.: Discussion: Reporting and Calibration of Post-Bomb ¹⁴C Data, Radiocarbon, 46, 1299–1304, https://doi.org/10.1017/S0033822200033154, 2004.

Sasse, J., Martinoia, E., and Northen, T.: Feed Your Friends: Do Plant Exudates Shape the Root Microbiome?, Trends Plant Sci., 23, 25–41, https://doi.org/10.1016/j.tplants.2017.09.003, 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, https://doi.org/10.1111/nph.17604, 2021.

Synal, H.-A., Stocker, M., and Suter, M.: MICADAS: A new compact radiocarbon AMS system, Nucl. Instrum. Meth. B, 259, 7–13, https://doi.org/10.1016/j.nimb.2007.01.138, 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 ¹⁴C 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 ¹³CO₂-labeling study, Biogeosciences, 17, 2971–2986, https://doi.org/10.5194/bg-17-2971-2020, 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 agro-nomic traits in arid environments, Nat. Biotechnol., 35, 969–976, https://doi.org/10.1038/nbt.3943, 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, https://doi.org/10.1371/journal.pone.0019563, 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, https://doi.org/10.1016/j.nimb.2009.10.067, 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, https://doi.org/10.1111/nph.13966, 2016.