Rootstock-Mediated Genetic Variance in Cadmium Uptake by Juvenile Cacao (Theobroma cacao L.) Genotypes, and Its Effect on Growth and Physiology

Grafting typically offers a shortcut to breed tree orchards throughout a multidimensional space of traits. Despite an overwhelming spectrum of rootstock-mediated effects on scion traits observed across several species, the exact nature and mechanisms underlying the rootstock-mediated effects on scion traits in cacao (Theobroma cacao L.) plants often remain overlooked. Therefore, we aimed to explicitly quantify rootstock-mediated genetic contributions in recombinant juvenile cacao plants across target traits, specifically cadmium (Cd) uptake, and its correlation with growth and physiological traits. Content of chloroplast pigments, fluorescence of chlorophyll a, leaf gas exchange, nutrient uptake, and plant biomass were examined across ungrafted saplings and target rootstock × scion combinations in soils with contrasting levels of Cd. This panel considered a total of 320 progenies from open-pollinated half-sib families and reciprocal full-sib progenies (derived from controlled crosses between the reference genotypes IMC67 and PA121). Both family types were used as rootstocks in grafts with two commercial clones (ICS95 and CCN51) commonly grown in Colombia. A pedigree-based best linear unbiased prediction (A-BLUP) mixed model was implemented to quantify rootstock-mediated narrow-sense heritability (h²) for target traits. A Cd effect measured on rootstocks before grafting was observed in plant biomass, nutrient uptake, and content of chloroplast pigments. After grafting, damage to the Photosystem II (PSII) was also evident in some rootstock × scion combinations. Differences in the specific combining ability for Cd uptake were mostly detected in ungrafted rootstocks, or 2 months after grafting with the clonal CCN51 scion. Moderate rootstock effects (h²> 0.1) were detected before grafting for five growth traits, four nutrient uptake properties, and chlorophylls and carotenoids content (h² = 0.19, 95% CI 0.05–0.61, r = 0.7). Such rootstock effects faded (h²< 0.1) when rootstock genotypes were examined in soils without Cd, or 4 months after grafting. These results suggest a pervasive genetic conflict between the rootstock and the scion genotypes, involving the triple rootstock × scion × soil interaction when it refers to Cd and nutrient uptake, early growth, and photosynthetic process in juvenile cacao plants. Overall, deepening on these findings will harness early breeding schemes of cacao rootstock genotypes compatible with commercial clonal scions and adapted to soils enriched with toxic levels of Cd.

Introduction

Grafting is an ancient technique used to propagate plants vegetatively by combining desirable agronomic traits of the rootstock with those of the scion. It has been used in several plant species targeting a wide spectrum of traits, for instance, resistance to pathogens (Boughalleb et al., 2007; Rivard and Louws, 2008; Spanò et al., 2020); tolerance to abiotic stress factors, such as water deficit, (Liu et al., 2016), heavy metals (Rouphael et al., 2008), and salinity (Yanyan et al., 2018; Suárez-Hernández et al., 2019); improved fruit quality (Davis and Perkins-Veazie, 2005); higher yields (Cardinal et al., 2007); and architectural changes in scions (Eltayb et al., 2014). However, breeding rootstocks for tree crops is slower than scion breeding for the same species. This is due to the long generation times and strong testing requirements of rootstocks, which reduce the opportunity for comprehensively testing their performance against multiple scions and environments.

Conventional propagation by rootstocks and grafting has routinely been used to expand cacao (Theobroma cacao L.) cultivation in tropical areas (Ríos, 1957). Seedling rootstocks resistant to Ceratocystis spp. and tolerant to acidic soils are typically obtained by open pollination (OP) of the IMC67, PA121, and PA46 reference genotypes in seedling orchards, and are grafted with susceptible productive clones (Palencia Calderon et al., 2007). Yet, only a few studies have evaluated the effect of the rootstock genotype on key agronomic traits such as yield and disease resistance (Yin, 2004; Ribeiro et al., 2016; Asman et al., 2021). For instance, Yin (2004) detected a significant effect of the rootstock on the vigor of cacao scions, but was unable to capture an influence on yield components such as bean weight and number of beans per pod. In line with these results, Asman et al. (2021) observed little rootstock effects on the scion’s resistance to vascular streak dieback, caused by Ceratobasidium theobromae. On the other hand, Ribeiro et al. (2016) observed a significant influence of the rootstock × scion interaction on the scion’s resistance to witches’s broom disease caused by Moniliophthora perniciosa, which allowed the identification of an elite rootstock genotype with a positive effect on scion disease resistance trait.

Despite the long time during which grafting has been used in the vegetative propagation of T. cacao, the rootstock effects on the expression of key scion’s agronomic traits remains poorly understood, among these, cadmium (Cd) accumulation. Content of Cd in cacao products is one of the most limiting factors for cocoa sale in international markets. Cadmium, a heavy metal that causes health problems, accumulates in the seeds of T. cacao, which are the raw material to produce chocolate (Barraza et al., 2017). Considering that Cd can accumulate in the human body, starting on January 1, 2019 the European Union (EU) began to control the maximum Cd concentrations allowed in chocolate and cocoa products imported to the EU (European Commission, 2014). Other countries are expected to implement similar regulations to Cd concentration in cocoa products (Vanderschueren et al., 2021), which has generated concern in countries where levels of bioavailable Cd for the plant have been detected in soil. Particularly, high Cd concentrations in cacao beans have been reported in plantations from South America (Chavez et al., 2015; Arévalo-Gardini et al., 2017; Vanderschueren et al., 2019), which has been correlated with the naturally high Cd content in the young soils of this region (Argüello et al., 2019). The Cd concentration in cacao bean is related to soil properties such as total soil Cd, pH, percentage of organic carbon, and oxalate extractable manganese (Argüello et al., 2019). Cd adversely affects growth, photosynthetic process, nutrient uptake, content of chloroplast pigments, cell structure, antioxidative metabolism, and gene expression (Dias et al., 2012; Jiao et al., 2012; Borek et al., 2013; Saidi et al., 2013; Castro et al., 2015; Farooq et al., 2016; Pereira de Araújo et al., 2017).

Cadmium enters root cells through ion channels for Zn, Fe, Mn, and Ca (Song et al., 2017). Transporter gene families like zinc-iron permease (ZIP), natural resistance-associated macrophage proteins (NRAMPs), and heavy metal transporting ATPases (HMAs) have been associated with uptake and translocation of Cd in plants (Vanderschueren et al., 2021). According to some studies, TcNRAMP5 may play a role in the regulation of Cd uptake in cacao plants (Ullah et al., 2018; Moore et al., 2020), whereas HMA-family proteins may contribute to Cd sequestration (Moore et al., 2020). However, the information regarding the specific role of transporter genes in cacao is still limited, as well as the regulatory mechanisms of Cd translocation to the shoot. Within the roots, Cd can be transported in plants through a symplasmic pathway. Cd is later loaded from the symplasm into the xylem (Lux et al., 2011). Movement of Cd from the roots to above-ground tissues depends on the plant mechanisms for vacuolar sequestration, xylem loading, and xylem to phloem transfer (Vanderschueren et al., 2021). In soybean and eggplant, the effect of the root system on the accumulation of Cd in the aerial part of the plant has been demonstrated by grafting cultivars on genotypes that differ in their accumulation of Cd (Sugiyama et al., 2007; Arao et al., 2008). The feasibility of mitigating Cd uptake and accumulation through grafting of highly productive scions on top of rootstocks with low Cd uptake capability has been suggested for T. cacao (Lewis et al., 2018; Engbersen et al., 2019). However, as far as we are aware, such effect has not yet been experimentally demonstrated. Further investigation is needed to elucidate the effect of the rootstock in Cd accumulation in cacao. Rootstock genotype may produce unpredicted responses in grafted plants exposed to stressful Cd conditions compared with self-rooted plants (Savvas et al., 2010).

Assessing the effect of rootstock genotypes or specific rootstock × scion combinations on the accumulation of Cd in plant tissues, and the overall plant tolerance to heavy elements, is therefore key for the establishment of new cacao plantations in regions with high levels of available Cd in soil. Therefore, the aim of this study was to quantify the inheritance of the rootstock effects on the scion’s Cd accumulation, as well as associated growth and physiological traits of juvenile cacao plants. Both open-pollinated (OP) and reciprocal full-sib progenies (derived from controlled crosses between the reference genotypes IMC67 and PA121) were used as rootstocks in grafts with two commercial and widely adapted clones (ICS95 and CCN51) in Colombia. Content of chloroplast pigments, fluorescence of chlorophyll a, leaf gas exchange, ion leakage, protein content, nutrient uptake, and plant biomass were also examined across rootstock × scion combinations in soils with contrasting levels of Cd. We hypothesize that due to the complex acquisition and transportation of Cd from the root system throughout the plant, there is a potential to harness rootstock-driven genetic variance for Cd uptake as part of early selection schemes at seed orchards, nurseries, and ungrafted saplings.

Materials and Methods

Location

The experiment was established under greenhouse conditions at the Palmira Research Station of AGROSAVIA, located in Palmira (3°31′12″N, 76°19′50″W), province of Valle del Cauca, Colombia, at an altitude of 1,001 masl. The climatic conditions during the experiment were recorded using a weather station WatchDog 1000 series Micro Station. Average annual temperature was 26°C, and average relative humidity was 65%.

Soil Substrate

Soil collected in a cacao growing region of southwest Colombia (pH 4.5) had a natural total Cd content of 0.43 mg kg^–1. The collected soil was air-dried and passed through a 2-mm mesh sieve. Cd was added to half of the sieved soil using an aqueous solution of Cd(NO₃)₂. Cd-spiked soil was incubated at field capacity for 1 month in a greenhouse. During the incubation time, the soil was constantly mixed and maintained at field capacity. The other half of the sieved soil was kept untouched as control treatment.

Greenhouse substrate was prepared by mixing Cd-spiked soil with rice husk and sand in a ratio of 3:1:1 to reach a final content of 7.49 mg Cd kg^–1, as determined by inductively coupled plasma optical emission spectrometry (Thermo Scientific ICAP 6500). Greenhouse substrate was also prepared using soil without the addition of Cd. According to soil analysis, carried out before the establishment of the experiment, the total nitrogen content in the substrate was 0.4 and 0.5% for soil without addition of Cd and soil enriched with Cd, respectively, which shows that there was no advantage in the concentration of nitrogen in the treatment with addition of Cd.

Plant Material

Plant material was obtained from the Colombian Cacao Germplasm Bank. Full-sib progenies obtained from the crossing between IMC67 × PA121, and its reciprocal cross, were established in Cd-spiked soil, as well as in soil without addition of Cd. Progenies were subsequently evaluated for plant growth, and physiological and nutritional parameters. Meanwhile, a total of two OP half-sib seedling families of IMC67 and PA121 were also considered. These accessions, PA121 and IMC67, corresponded to genotypes recommended as rootstocks in Colombia due to their resistance to pathogens and adaptation to soils with acidic pH (Palencia Calderon et al., 2007). The mucilage of the seed was removed by gently rubbing with sawdust. Seeds were sown in plastic conical tubes containing sand and grown under these conditions for 2 months. After this period, well-developed seedlings were transplanted into black polyethylene bags of 45 cm high and 29 cm wide, containing 10 kg greenhouse substrate with and without the addition of Cd. A commercial fertilizer containing N, P, K, and Mg was applied.

The experiment was established in a randomized complete block design with four replicates arranged in a 4 × 2 factorial design and 10 plants per experimental unit. The treatments corresponded to OP families IMC67 and PA121, as well as IMC67 × PA121 and PA121 × IMC67 progenies from controlled crosses established in soil with and without the addition of Cd. Seedlings were grown on the Cd treatments for 5 months and then evaluated for plant growth, physiological, and nutritional parameters at the end of this period.

Five months after Cd treatment, budwoods of the commercial cultivars ICS95 and CCN51 were grafted onto the seedling rootstocks using the top grafting technique as described by Isele et al. (2020). Taking into consideration that each experimental unit had 10 plants, two of each were evaluated without grafting, four of the remaining plants were grafted with ICS95, and the other four with CCN51. Two and 4 months after grafting, rootstock × scion combinations were evaluated for plant growth, physiological, and nutritional parameters.

Chlorophyll Fluorescence

Chlorophyll fluorescence was determined by means of a portable optical pulse fluorometer (Opti-Sciences OS30P+). The measurements were made between 7:00 and 11:00 on young, fully expanded, healthy, and photosynthetically active leaves of two plants per experimental unit, usually the third or fourth leaf from the apex of the plants. The leaf was adapted to darkness, using suitable clips, for a period of 30 min and subsequently illuminated with a saturating actinic light pulse of 3,500 μmol m^–2 s^–1 for 1 s. The following parameters were evaluated: initial fluorescence (F₀), maximum fluorescence (F_m), and the maximum quantum yield of PSII (F_v/F_m), obtaining two data per plant.

Leaf Gas Exchange

Net photosynthetic rate (P_N), stomatal conductance to water vapor (g_s), leaf transpiration (E), internal carbon dioxide concentration (Ci), and instantaneous water use efficiency (WUE_INST) were measured on the same leaf used to determined chlorophyll fluorescence. Gas exchange measurements were performed between 7:00 and 11:00 using an open system portable gas analyzer (ADC model LCpro+), natural light intensity between 250 and 450 μmol m^–2 s^–1 and atmospheric CO₂ of 390 μmol (CO₂) mol^–1. The leaf area sampled was 6.25 cm². Two data points per plant were obtained.

Chloroplast Pigments

Chlorophyll a, chlorophyll b, their ratio, total chlorophyll, and carotenoids were quantified in the same leaves used to measure gas exchange and fluorescence emission. Pigments were extracted from three disk-shaped leaf segments per plant using cold 80% ammoniacal acetone (4°C), following the protocol described by Melgarejo et al. (2010) and calculated according to Lichtenthaler (1987). The area of each disk was 0.78 cm². The absorbance of the extracts was determined with a spectrophotometer (PerkinElmer Lambda 25) using the following wavelengths: 663, 647, and 470 nm to calculate the concentrations of chlorophyll a, chlorophyll b, and total carotenoids, respectively (Lichtenthaler, 1987).

Pigment concentration was calculated from the following equations:

$\mathrm{Ca}=12.25\mathrm{A}{}_{663-2.79}\mathrm{A}{}_{647}$

$\mathrm{Cb}=21.5\mathrm{A}{}_{647-5.1}\mathrm{A}{}_{663}$

$\mathrm{C}(\mathrm{a}+\mathrm{b})=7,15\mathrm{A}{}_{663+18.7}\mathrm{A}{}_{647}$

$\mathrm{Cx}+\mathrm{c}=1000\mathrm{A}{}_{470-1.82}\mathrm{Ca}-85,02\mathrm{Cb}\mathrm{/}198$

where Ca stands for chlorophyll a, Cb stands for chlorophyll b, C (a + b) stands for total chlorophyll, and Cx + c stands for carotenoids.

Electrolyte Leakage

To determine electrolyte leakage (EL), 0.78 cm² leaf disks were taken from the same leaf used to determine the content of chloroplast pigments. Three disks were obtained from each leaf sample, washed with distilled water, and then arranged in 15 mL conical Falcon tubes containing 3 mL of deionized water. The samples were incubated at room temperature (22°C) on a shaker for a period of 5 h and the electrical conductivity was recorded at the end of this period (EC1). The Falcon tubes were then taken to a water bath at 80°C for 10 min and electrical conductivity (EC2) recorded again after cooling the bathing solution to 22°C. The loss of electrolytes (EL) was calculated using the following formula:

$EL=\frac{EC1}{EC2}\times 100$

EC1 corresponds to the electrical conductivity measured after 5 h in deionized water. EC2 corresponds to the electrical conductivity measured after incubation at 80°C.

Plant Growth and Biomass

Once the physiological parameters were evaluated during morning hours, seedling rootstocks, and also grafting combinations, were harvested for shoot and root lengths measures. Plants were then divided into roots, stem, and leaves. Plant tissue was washed with distilled water and gently dried with absorbent towels. Fresh weight of each plant portion was then determined. To wash the excess of Cd adhered to the roots, these were immersed for 3 min in a solution containing 5 mM EDTA and 20 mM TRIS (pH 8.0) with constant agitation followed by three washes with distilled water for 3 min (Nguyen et al., 2016). The reason why roots were washed was to prevent externally adhered Cd from being quantified as if it were inside the root system. Samples were placed in paper bags and oven-dried at 60°C for 72 h, and after this period dry weight was determined.

Leaf Area

Leaf area was measured using an electronic area meter (LICOR-3000) immediately after harvesting leaf tissue.

Protein Content

Approximately 100 mg of leaf tissue were macerated in liquid nitrogen. Phosphate buffer (0.05 M) was then added, and the samples centrifuged for 30 min at 12,000 × g and 4°C. Total soluble protein content was measured by the Coomassie blue method according to Bradford (1976) using BSA as standard in a microplate spectrophotometer (BioTek EPOCH) at 595 nm.

Cadmium and Mineral Nutrients Content

Dried plant material was ground in a mill (Thomas Wiley) and submitted to nitric–perchloric digestion assisted by microwave (Milestone UltraWave). Mineral elements such as Cd, Ca, Mg, Fe, Zn, Cu, and Mn were determined by atomic absorption spectrophotometry, P by colorimetry, and K by flame emission photometry (Isaac and Kerber, 1971). Nitrogen content was determined by the Kjeldahl method (Bradstreet, 1954).

Data Analyses and Rootstock-Mediated Heritability Scores

Trait variation across families was explicitly compared using a mixed linear model in which family identity and Cd treatment were indexed as fixed effect and repetition as random effect. Second, rootstock-mediated additive genetic variance (i.e., narrow sense hereditability) was obtained for each trait using a best linear unbiased prediction (BLUP) mixed model that relied on the pedigree information (i.e., a pedigree relationship matrix for an A-BLUP model) from the rootstocks (half-sib families for OP, and full-sib families for controlled crossing schemes, allowing for a 5% of maternal effects). The corresponding linear mixed model, following Arenas et al. (2021), was fitted as in:

$y=Xb+Za+e$

where y is the phenotypic trait vector, a is a vector of individual random additive genetic effects with a normal distribution so that a ∼ N(0, Aσ_a²) with A being the pedigree relationship matrix among rootstocks and σ_a² the additive genetic variance (Müller et al., 2017), b is a vector of fixed effects (i.e., intercept or general mean, and experimental site effects), e is the vector of residual error effect, and X and Z are the corresponding incidence matrices for fixed effects and additive genetic effects (Chen et al., 2018; Gutierrez et al., 2018), respectively.

To aid interpretability while accounting for clonal scion differences and in situ environmental (soil) variation (of major importance in Cd fixation), genetic parameters for all phenotypic traits were individually estimated for each environment and scion, rather than including the latter as explicit fixed effects within the model. All models implemented the reproducing kernel Hilbert space (RKHS) method throughout the BGLR (Bayesian Generalized Linear Regression) package (Pérez-Rodríguez and de los Campos, 2014) in R v.3.4.4 (R Core Team). RKHS is a semiparametric approach to infer a given function without making a strong a priori assumption about the distribution of effects (Cuevas et al., 2016). RKHS was executed using a Gibbs sampling with 10,000 iterations after discarding the 5,000 initial steps as burn-in. A thinning interval of 10 was implemented for data recording. Convergence of posterior distributions was verified using trace plots, whereas rootstock-mediated specific combining abilities, original trait correlations (Supplementary Figure 1), and A input pedigree relationship matrices (Supplementary Figure 2) were respectively checked using the R (R Core Team) functions lme (from the nlme package, treating families as fixed effects and repetitions as random effects), cor.test, and heatmap.

Narrow sense rootstock-mediated heritability (h²) was then computed using the additive (σ_a²) and residual (σ_e²) variances, following de los Campos et al. (2015), as in:

$h^2={\mathrm{\sigma }}_a^2\mathrm{/}\left({\mathrm{\sigma }}_a^2+{\mathrm{\sigma }}_e^2\right)$

Variance within experimental units was captured by the residual variance term because it did not correspond to an additive variance component. Narrow sense rootstock-mediated heritability (h²) estimates were summarized using the median and the 95% confidence interval from the BGLR’s posterior distribution. Overall models’ fits were examined by computing the prediction ability (r_y) estimated for each trait as the Pearson’s correlation between the observed and the predicted (i.e., breeding value deviation from the overall mean) trait values (Müller et al., 2017; Zhang et al., 2019).

Results

Cadmium Accumulation

Significant differences in the accumulation of Cd in roots before grafting were observed between the seedling rootstocks established in Cd-spiked soil. Accumulation of Cd in the roots of PA121 × IMC67 full-sib progenies was significantly higher (p < 0.05) than in IMC67 OP half-sib and IMC67 × PA121 full-sib families (Figure 1A). On the other hand, no significant differences were observed in the accumulation of Cd in leaves between the family rootstocks subjected to stress by Cd (Figure 1B).

FIGURE 1

Figure 1. Cadmium accumulation across four ungrafted seedling rootstocks grown in soil enriched with Cd (7.49 mg Cd kg^–1 soil) 5 months after treatment. (A) Cadmium accumulation in roots. (B) Cadmium accumulation in leaves. Data are mean values ± SE (n = 8).

Cadmium accumulation in both rootstocks and scions was quantified 2 and 4 months after grafting the ICS95 and CCN51 clonal scions on the seedling rootstock families (Figure 2). At 2 months after grafting, lower accumulation of Cd in leaf tissue was observed when the ICS95 clone was grafted on the IMC67 OP half-sib progenies (Figure 2A). However, no significant differences in Cd accumulation were observed between any of the rootstock families on which ICS95 was grafted (Figure 2C). Four months after grafting, a greater accumulation of Cd in leaf tissue of the ICS95 clonal scion grafted on PA121 × IMC67 full-sib rootstocks was observed (Figure 2A). For this case, lower accumulation of Cd in roots of the IMC67 OP half-sib and IMC67 × PA121 full-sub families was observed when grafted with the ICS95 clonal scion (Figure 2C).

FIGURE 2

Figure 2. Cadmium accumulation in rootstock × scion combinations 2 and 4 months after grafting. Seedling rootstocks families were established in soil enriched with Cd (7.49 mg Cd kg^–1 soil). Panels (A,C) depict cadmium accumulation in grafting combinations having ICS95 as scion. Panels (B,D) show Cd accumulation in grafting combinations having CCN51 as scion. Data are mean values ± SE (n = 8).

On the other hand, no significant differences were observed in Cd accumulation of the CCN51 clonal scion when grafted on any of the evaluated rootstocks 2 months after grafting (Figure 2B). However, a greater accumulation of Cd was observed in PA121 × IMC67 full-sib progenies when they were grafted with the CCN51 scion (Figure 2D). Four months after grafting, a greater accumulation of Cd was observed in CCN51 clonal scion grafted on IMC67 × PA121 full-sib rootstocks, and its reciprocal progenies (Figure 2B). For this case, no significant differences in Cd accumulation were observed between rootstock families having CCN51 as scion (Figure 2D).

Cadmium accumulation in the leaf tissue of the CCN51 scion grafted on IMC67 OP half-sib rootstock progenies, 2 and 4 months after grafting, was significantly higher than for the ICS95 clonal scion (Figure 3A). No significant differences in Cd accumulation were observed between ICS95 and CCN51 when grafted onto any of the other rootstock families at 2 or 4 months after grafting (Figure 3A). Significant differences in Cd accumulation in the roots were observed only for the PA121 × IMC67 full-sib rootstock family. A higher accumulation of Cd was observed when this rootstock family was grafted with the CCN51 scion 2 months after grafting. However, 4 months after grafting a higher accumulation of Cd in the roots was observed when using ICS95 as scion (Figure 3B).

FIGURE 3

Figure 3. Cadmium accumulation of ICS95 and CCN51 clonal scions grafted on four seedling rootstock families 2 and 4 months after grafting. Seedling rootstocks (x-axis) were grown in soil enriched with Cd (7.49 mg Cd kg^–1 soil). (A) Cadmium accumulation in leaf tissue of the ICS95 and CCN51 clonal scions. (B) Cadmium accumulation in the roots of the four OP rootstock families when grafted with the ICS95 and CCN51 clones. Blue bars represent grafting combinations with ICS95 as scion 2 months after grafting. Green bars represent grafting combinations with CCN51 as scion 2 months after grafting. Data are mean values ± SE (n = 8). Dashed bars show Cd accumulation 4 months after grafting.

Plant Growth

A significant decrease in plant biomass, shoot length, and leaf area was observed 5 months after establishment of the ungrafted seedlings in Cd-spiked soil (Figures 4A–F). Fresh weight of the stem decreased by 13.7, 37.6, 8.3, and 13.2% for the OP half-sib families IMC67 and PA121, and for the full-sib IMC67 × PA121 and PA121 × IMC67 progenies, respectively, as compared with the control (without Cd addition) (Figure 4A). A greater decrease in stems’ (37.4%) and roots’ (19.5%) dry weight was also observed for PA121 OP half-sib rootstocks subjected to stress by Cd addition in the soil (Figures 4B,D). A similar response was observed for shoot length, for which a decrease of 15.8, 32, 13.6, and 2.8% was observed respectively, for OP half-sib progenies IMC67 and PA121, and full-sib IMC67 × PA121 and PA121 × IMC67 families (Figure 4E).

FIGURE 4

Figure 4. Plant growth of four families of ungrafted seedlings under greenhouse conditions at two different concentrations of Cd in soil, 5 months after cadmium treatment. (A) Stem fresh weight, (B) stem dry weight, (C) leaf fresh weight, (D) root dry weight, (E) shoot length, and (F) leaf area. White bars represent rootstocks grown in soil without addition of Cd (0.43 mg kg^–1). Gray bars represent rootstocks grown in Cd spiked soil (7.49 mg kg^–1). Data are mean values ± SE (n = 8).

Two months after grafting, significant changes in root biomass were observed for grafting combinations having both ICS95 and CCN51 as scion (Figures 5A–D). Combinations of ICS95 as scion and PA121 half-sib progenies as rootstock showed a significant decrease in root biomass. Specifically fresh and dry root weight fell in 25.5 and 24.6%, respectively, compared with the control (Figures 5A,C). Likewise, 2 months after grafting, PA121 half-sib rootstock families showed a significant decrease in root biomass in combination with the CCN51 clonal scion. For this case, fresh and dry root weight respectively, decreased by 35.13 and 36.15% compared with the control (without Cd addition) (Figures 5B,D).

FIGURE 5

Figure 5. Root biomass of rootstock × scion combinations 2 months after grafting seedling rootstocks subjected to two different concentrations of Cd in soil under greenhouse conditions. Panel (A,B) depict fresh root weight of rootstock families having ICS95 and CCN51 as scions. Panel (C,D) show dry root weight of OP and full-sib rootstock progenies having ICS95 and CCN51 as scion. White bars represent rootstocks grown in soil without addition of Cd (0.43 mg kg^–1). Blue bars represent rootstocks grown in Cd spiked soil (7.49 mg kg^–1) grafted with the ICS95 clonal scion. Green bars represent rootstocks grown in Cd spiked soil (7.49 mg kg^–1) grafted with the CCN51 clonal scion. Data are mean values ± SE (n = 8).

Four months after grafting, a significant decrease in root biomass was also observed in all rootstock families grafted with the ICS95 scion (Figures 6A,C). On the other hand, 4 months after grafting, seedling rootstocks grafted with the CCN51 clonal scion showed significant changes not only in root growth parameters but also in the scion growth. A significant decrease in leaf area and in fresh and dry weight of leaf tissue was observed in the CCN51 scion grafted on all the rootstock families, except when it was grafted on PA121 × IMC67 full-sib rootstocks (Figures 6B,D). A decrease in dry root weight of 23.5 and 37.2% was respectively, observed in the OP half-sib rootstock families IMC67 and PA121, both grafted with the CCN51 clonal scion (to be presented in section Rootstock-Mediated Specific Combining Abilities and Heritability Scores in Full-Sib Families).

FIGURE 6

Figure 6. Plant growth of rootstock × scion combinations 4 months after grafting at two different concentrations of Cd in soil under greenhouse conditions. Panels (A,C) show fresh and dry root weight of seedling rootstocks grafted with the ICS95 clonal scion. Panels (B,D) depict fresh and dry leaf weight of the CCN51 clonal scion across rootstock families. White bars represent rootstocks grown in soil without addition of Cd (0.43 mg kg⁻¹). Blue bars represent rootstock families grafted with the ICS95 clonal scion established in Cd spiked soil (7.49 mg kg^–1). Green bars represent rootstock families grafted with the CCN51 clonal scion established in Cd spiked soil (7.49 mg kg^–1). Data are mean values ± SE (n = 8).

Chlorophyll Fluorescence and Leaf Gas Exchange

A decrease in the net photosynthetic rate (P_N) and in instantaneous water use efficiency (WUE_INST) was observed in response to Cd toxicity before grafting (Figure 7). A significant decrease (p < 0.05) in the carbon assimilation rate equivalent to 40.09%, compared with the control, was observed when using PA121 OP half-sib seedling rootstocks (Figure 7A). Water use efficiency respectively, decreased by 30.46 and 40.6% when using IMC67 and PA121 OP half-sib seedling rootstock families, compared with control plants (Figure 7B). There was no effect of Cd on leaf gas exchange variables when using ICS95 and CCN51 as clonal scions, neither at 2 nor at 4 months after grafting.

FIGURE 7

Figure 7. Leaf gas exchange in four seedling rootstock families subjected to two different concentrations of Cd in soil under greenhouse conditions 5 months after cadmium treatment. (A) Photosynthetic rate (P_N), and (B) instantaneous water use efficiency (WUE_INST). White bars represent rootstock families grown in soil without addition of Cd (0.43 mg kg^–1). Gray bars represent rootstock families grown in Cd spiked soil (7.49 mg kg^–1). Data are mean values ± SE (n = 8).

Variables of chlorophyll a fluorescence emission did not show significant differences (p < 0.05) between seedling rootstock families established in soil with and without the addition of Cd before grafting (to be presented in section Rootstock-Mediated Specific Combining Abilities and Heritability Scores in Full-Sib Families). However, 2 months after grafting, a significant decrease in F_m was observed in ICS95 clonal scion grafted on all the rootstock families. A greater decrease of 19.4% in F_m was observed in this clonal scion grafted on IMC67 × PA121 full-sib rootstock progenies (Figure 8A). An effect of Cd on F_m resulted in a change in F_v/F_m ratio scores. A significant decrease in F_v/F_m was observed in the ICS95 clonal scion grafted on PA121 OP half-sib and IMC67 × PA121 full-sib rootstock families 2 months after grafting (Figure 8B). On the other hand, an increase in F_v/F_m ratio was observed in the same clone when grafted on the reciprocal PA121 × IMC67 full-sib progenies 2 and 4 months after grafting (Figures 8B,C), which suggests for a maternal effect. There was no significant difference when the clonal scion CCN51 was grafted with the different seedling rootstock families for the variables of chlorophyll fluorescence at 2 or at 4 months after grafting (to be presented in section Rootstock-Mediated Specific Combining Abilities and Heritability Scores in Full-Sib Families).

FIGURE 8

Figure 8. Chlorophyll fluorescence in the ICS95 clonal scion when grafted with 4 different rootstock families at two different concentrations of Cd in soil under greenhouse conditions. (A) F_m 2 months after grafting. (B) F_v/F_m 2 months after grafting. (C) F_v/F_m 4 months after grafting. White bars represent rootstock families grown in soil without addition of Cd (0.43 mg kg^–1). Blue bars represent rootstock families grown in Cd spiked soil (7.49 mg kg^–1). Dashed bars show Cd accumulation 4 months after grafting. Data are mean values ± SE (n = 8).

Chloroplast Pigments

Significant changes in the content of chloroplast pigments were observed in the ungrafted seedlings 5 months after establishment in Cd-spiked soil (Figure 9). A different effect of Cd was observed depending on the rootstock. A significant increase in chlorophyll a equivalent to 41.6 and 21.8% was respectively, observed in IMC67 OP half-sib and PA121 × IMC67 full-sib progenies, compared with the control (Figure 9A). IMC67 OP half-sib and PA121 × IMC67 full-sib families also showed a significant increase in chlorophyll b content respectively, equivalent to 32 and 29% (Figure 9B). Changes in chlorophyll a and b resulted in changes in total chlorophyll, accordingly. A significant increase of total chlorophyll was observed in IMC67 OP half-sib and PA121 × IMC67 full-sib families (Figure 9C). The carotenoid content also increased by 17% in IMC67 OP half-sib progenies 5 months after Cd addition. The content of carotenoids in the remaining families decreased (Figure 9D).

FIGURE 9

Figure 9. Content of chloroplast pigments of ungrafted seedling families subjected to two different concentrations of Cd in soil under greenhouse conditions 5 months after cadmium addition. (A) Content of chlorophyll a. (B) Content of Chlorophyll b. (C) Content of total Chlorophyll. (D) Content of carotenoids. White bars represent seedling families grown in soil without addition of Cd (0.43 mg kg^–1). Gray bars represent seedling progenies grown in Cd spiked soil (7.49 mg kg^–1). Data are mean values ± SE (n = 8).

Two months after grafting, a significant decrease in all chloroplast pigments was observed across seedling rootstock families grafted with both ICS95 and CCN51 clonal scions (Figure 10).

FIGURE 10

Figure 10. Content of chloroplast pigments in ICS95 and CCN51 clonal scions 2 months after grafting on four different seedling rootstock families subjected to two different concentrations of Cd in soil under greenhouse conditions. Panels (A,B) show content of Chlorophyll a of ICS95 and CCN51 clonal scions. Panels (C,D) depict content of Chlorophyll b of ICS95 and CCN51 clonal scions. Panels (E,F) exhibit content of total Chlorophyll of ICS95 and CCN51. Panels (G,H) show content of carotenoids of ICS95 and CCN51 clonal scions. White bars represent rootstocks × scion combinations grown in soil without addition of cadmium (0.43 mg kg^–1). Blue bars show seedling rootstock families grafted with the ICS95 clonal scion grown in Cd spiked soil (7.49 mg kg^–1). Green bars depict seedling rootstock families grafted with the CCN51 clonal scion grown in Cd spiked soil (7.49 mg kg^–1). Data are mean values ± SE (n = 8).

Electrolyte Leakage

A significant increase of 24.7 and 30.6% in electrolyte leakage was observed in ungrafted IMC67 and PA121 OP half-sib families established in Cd-spiked soil, suggesting cell membrane damage in these rootstocks (Figure 11). No significant difference was observed in electrolyte leakage after grafting with the ICS95 and CCN51 clonal scions in Cd-spiked soil (Supplementary Table 5).

FIGURE 11

Figure 11. Electrolyte leakage in four ungrafted seedling families subjected to two different concentrations of Cd in soil under greenhouse conditions 5 months after Cd addition. White bars represent ungrafted seedling families grown in soil without the addition of Cd (0.43 mg kg^–1). Gray bars represent ungrafted seedling families grown in Cd-spiked soil (7.49 mg kg^–1). Data are mean values ± SE (n = 8).

Protein Content

A decrease of 22.2% in total soluble protein content was observed in PA121 OP half-sib seedlings before grafting (Figure 12A). Two months after grafting, a significant decrease in total soluble protein content of 39.5 and 38.5% was observed respectively in the ICS95 clonal scion grafted on the IMC67 OP half-sib, and IMC67 × PA121 full-sib families, compared to the control plants without Cd addition (Figure 12B). A significant decrease in protein content in leaf tissue of the CCN51 clonal scion grafted on IMC67 OP half-sib progenies, equivalent to 31% compared with the control, was also observed 2 months after grafting (Figure 12C).

FIGURE 12

Figure 12. Protein content in leaf tissue of rootstocks and grafting combinations subjected to two different concentrations of Cd in soil under greenhouse conditions. (A) Protein content in leaf tissue of ungrafted seedling families 5 months after Cd addition. (B) Protein content in leaf tissue of the ICS95 clonal scion 2 months after grafting with four different rootstock families. (C) Protein content in leaf tissue of the CCN51 clonal scion 2 months after grafting with four different rootstock families. White bars represent families grown in soil without addition of Cd (0.43 mg kg^–1). Gray bars represent ungrafted seedling families grown in Cd-spiked soil (7.49 mg kg^–1). Blue bars represent seedling rootstock families grafted with the ICS95 clonal scions in Cd-spiked soil (7.49 mg kg^–1). Green bars represent seedling rootstock families grafted with the CCN51 clonal scion in Cd-spiked soil (7.49 mg kg^–1). Data are mean values ± SE (n = 8).

Mineral Nutrients Content in Leaf Tissue

Cadmium absorption resulted in changes in the accumulation of mineral elements in leaf tissue of ungrafted seedling (Supplementary Table 1). Significant increases of 14 and 17.3% in K content were observed in PA121 OP half-sib and PA121 × IMC67 full-sib families compared with the control. N content decreased in all seedling families, while Mg content increased. However, these changes were significant only in IMC67 × PA121 full-sib progenies, compared with the control (without Cd addition). A significant increase in P content in leaf tissue was observed only in IMC67 × PA121 full-sib family established in Cd-spiked soil. In relation to micronutrients, changes in Fe content in leaf tissue were significant only in IMC67 OP half-sib progenies, for which there was a reduction of 16.3% in relation to the control. Zn and Mn content decreased in all seedling families subjected to stress by Cd. A significantly higher decrease in the Zn content of 25.8% was observed in PA121 OP half-sib progenies, compared with the control. Mn content respectively decreased by 10.9, 24.4, 11.5, and 21.5% for PA121 and IMC67 OP half-sib, and PA121 × IMC67 and IMC67 × PA121 full-sib families, in relation to the control.

Mineral elements concentration in leaf tissue of clonal scions across seedling rootstock families is shown in Supplementary Table 2. Two months after grafting, significant changes were observed in the concentrations of macro- and microelements in both CCN51 and ICS95.

Rootstock-Mediated Specific Combining Abilities and Heritability Scores in Full-Sib Families

In ungrafted full-sib progenies obtained from the crossing between IMC67 × PA121 and its reciprocal cross, specific combining abilities were only significantly different for the Cd treatment (Supplementary Table 3). Specifically, they differed significantly in 23 of the 38 examined traits, including Cd content in roots, eight nutrient uptake traits, six biomass traits, three cationic interchange traits, and five pigments’ content traits (Table 1). Meanwhile, only two traits exhibited heritability (h²) scores above 0.1 under the non-Cd-added treatment, as compared with seven from the Cd-added treatment (Supplementary Table 4). Of the latter, all matched with significantly different specific combining abilities (Table 1), which speaks for a predominant role of additive genetic effects for ungrafted progenies.

TABLE 1

Table 1. Specific combining abilities and narrow-sense heritability (h²) in cacao seedlings from controlled crosses before grafting.

Concerning cacao full-sib seedlings from controlled crosses between IMC67 × PA121 (and reciprocal) 2-months after grafting with the ICS95 and CCN51 scions, specific combining abilities were mostly different under the Cd-added treatment (Supplementary Table 5). Specifically, these differences concerned 27 traits from the overall 35 examined traits (Table 2). Mismatches due to the scion origin were overserved in 14 cases (including Cd content in leaves and roots). Similarly, when examining the very same cacao seedlings from these controlled crosses 4 months after grafting, specific combining abilities differed for nine of the 17 studied traits (Supplementary Table 6). Specifically, these differences only overlapped between the ICS95 and CCN51 clonal scions for Cd content in leaves and dry root weight (Table 2).

TABLE 2

Table 2. Specific combining abilities in cacao seedling rootstocks from controlled crosses 2 and 4 months after grafting with the ICS95 and CCN51 scions.

Most of these specific combining abilities from the full-sib rootstock families fell when looking at the narrow-sense rootstock-mediated heritability (h²) scores 2 months after grafting (Supplementary Table 7), and completely vanished after 4 months (Supplementary Table 8). When grafting the clonal scion ICS95, heritability scores above 0.1 were only observed 2 months after grafting for root weight, dry root weight, and protein content, whereas for the CCN51 scion this only applied for chlorophyll b (Table 3). Four-months after grafting, none h² estimate was above 0.1.

TABLE 3

Table 3. Narrow-sense rootstock-mediated heritability (h²) in cacao seedling rootstocks from controlled crosses 2 months after grafting with the ICS95 and CCN51 scions.

Discussion

Cadmium Toxicity Affects a Wide Spectrum of Growth and Physiological Traits Before and After Grafting

In this study, an inhibition in plant growth and biomass was observed as result of Cd toxicity. Cd showed a greater toxic effect on the growth of PA121 OP seedlings than in the other families used as rootstocks. Also, a significant decrease in root biomass of PA121 OP rootstocks was observed both in combinations with clones ICS95 and CCN51 as scions. An effect of toxic levels of Cd on decreasing plant growth and biomass accumulation has been reported for other plant species such as cotton (Farooq et al., 2016), mustard (Ahmad et al., 2011), peanut (Lu et al., 2013), lettuce (Dias et al., 2012), lupin (Zornoza et al., 2002), and Rorippa globose (Wei et al., 2012), but this is the first explicit quantification in cacao scion × rootstocks combinations. Two months after grafting, changes were observed in growth parameters, although these were significant only for root biomass. Interestingly, 4 months after grafting, rootstock families having ICS95 as scion continued exhibiting significant changes only for root biomass, whereas the rootstock families having CCN51 as scion already started presenting significant changes also for growth parameters measured on the aerial part such as leaf area and leaf biomass. Decreased leaf area may be the result of reduction in cells size or alternatively more condense intercellular spaces (Barceló et al., 1988). Such effect of Cd on plant growth may result from alterations in cell division (Liu et al., 2003).

Meanwhile, a reduced Cd translocation from root to shoot could be explained by a differential loading of Cd into the xylem (Arao et al., 2010). Differences in Cd uptake by a symplastic pathway could also be related to differences between genotypes in Cd translocation to the stem of rootstocks and to scion’s tissues (Arao et al., 2010). The latter would explain a lower accumulation of Cd in ICS95, grafted on IMC67, compared with CCN51 grafted on the same rootstock. Expression of specific transporter proteins could also explain differences in Cd uptake and accumulation in different cacao genotypes (Moore et al., 2020). Still, further studies are needed to evaluate Cd translocation in the xylem of rootstocks × scion combinations to increase knowledge of Cd distribution in plant organs.

Cadmium has an effect on different photosynthetic processes such as decrease in chlorophyll content, inhibition of chlorophyll formation, inhibition of Rubisco activity, inhibition of both photosynthesis reaction centers, PSI and PSII, and an increase in lipoxygenase activity (Rai et al., 2016). The inhibition of the growth of plants subjected to Cd stress could further be due to a toxic effect of the heavy metal on physiological processes, such as reduction of the maximum photochemical efficiency of PSII (Ge et al., 2015; Pereira de Araújo et al., 2017), in addition to an imbalance in the uptake of essential mineral elements (Zornoza et al., 2002; Hédiji et al., 2015) and to an inhibition in the production of sugars due to the decrease in the carbon assimilation rate, P_N (Dias et al., 2012). Ahmad et al. (2011) correlated the toxic effect of Cd on mustard plant growth with a Cd-induced decrease in P_N and WUE. Similarly, a greater decreasing effect of Cd concentration in P_N and WUE_INST was observed in our study for PA121 that, as indicated above, was the rootstock that presented a greater effect of Cd toxicity on plant growth. A decrease in WUE suggests alterations in water balance, which could be due to changes in water uptake and transport as a result of the toxic effect of Cd (Singh and Tewari, 2003).

Cadmium binds to several sites in PSII, affecting both the donor and the acceptor side. Cd also inhibits oxygen evolution in a high affinity site by competition with Ca on the donor side (Sigfridsson et al., 2004). Concerning underlying physiological processes, efficiency of PSII, measured as the fluorescence F_v/F_m ratio, did not appear to be affected in Cd-stressed rootstocks. However, significant changes in F_v/F_m ratio of scion ICS95 were observed at 2 and 4 months after grafting on IMC67 × PA121 full-sib family rootstocks established in Cd-spiked soil compared with the control. Interestingly, no significant differences were observed in the maximum efficiency of the reaction centers of the PSII of scion CCN51 grafted on any of the rootstocks subjected to Cd stress in soil. This result suggested that grafting combinations having CCN51 as scion had a more stable photosynthetic performance than those having ICS95 clone as scion. A similar effect of Cd on photosynthetic performance has also been reported in poplar (Ge et al., 2015; Jiao et al., 2015), pea (Sandalio et al., 2001), and lettuce (Dias et al., 2012). One of the factors that could result in an effect on photosynthetic performance is the decrease in chlorophyll content as a result of Cd toxicity (Sandalio et al., 2001). However, for the case of ungrafted rootstocks, an increase in the content of chloroplast pigments at IMC67 OP seedlings is observed and does not result in significant changes in the F_v/F_m ratio. On the other hand, changes in the photosynthetic performance of scion ICS95 2 months after grafting effectively coincided with a decrease in chlorophyll content. In the case of scion CCN51, although a decrease in chlorophyll content is observed under Cd stress, the photosynthetic performance remains unchanged.

Similarly, a significant decrease in the concentration of chloroplast pigments has been reported for other species such as cotton (Farooq et al., 2016), beans (Saidi et al., 2013), Phragmites australis (Pietrini et al., 2003), poplar (Jiao et al., 2015), maize (Ekmekçi et al., 2008), pea (Sandalio et al., 2001), and cucumber plants exposed to Cd (Zhang et al., 2002). In this regard, carotenoids play a role as antioxidants, and increases in carotenoid content in some species have been explained as an attempt to protect chlorophyll from the photooxidative damage caused by Cd stress (Mishra et al., 2006). Therefore, the decrease in carotenoid content observed in grafting combinations having ICS95 as scion may have jeopardized the detoxification of radicals formed in response to Cd stress (Ekmekçi et al., 2008).

Cadmium absorption by plants can also result in changes in the accumulation of essential mineral elements in plant tissues (Castro et al., 2015; Pereira de Araújo et al., 2017). Cd may affect the transport of mineral elements by disturbing the radial movement of transporters in the root, loading into the xylem vessels or into the leaves, thereby promoting morphological changes of the xylem tissue, changes in H⁺-ATPase activity, and alterations in IRT1 transporter selectivity (Sandalio et al., 2001). Changes in nutritional elements may indicate alterations in ionic homeostasis (Saidi et al., 2013). In our work, changes in mineral elements were observed in rootstock families before grafting, and in both ICS95 and CCN51 leaf tissue 2 months after grafting. Cd translocation in the plant is carried out using the same transporters than some nutritional elements such as Ca, Mn, and Zn (Rahman et al., 2016). Competition for these same transporters may explain the decrease in the content of Zn and Mn in ungrafted rootstocks under Cd stress.

On the other hand, according to Küpper et al. (1998), Cd could replace Mg in chlorophylls, which may decrease chlorophyll content in plants under Cd stress. A direct proportional relationship between toxicity of heavy metals and Mg substitution in chlorophylls has been reported (Küpper et al., 1996). In our study, a significant increase in Mg concentration in leaf tissue was observed 2 months after grafting in Cd-spiked soil. These results are in agreement with the observation made by Ciećko et al. (2005), in which the content of Mg in oat increased when it was grown in Cd-contaminated soil. An effect of Cd on electron transport on the reducing side of photosystem I was observed on isolated chloroplasts of maize plants grown in nutrient solution containing the heavy metal. A reduction in electron transport was associated with a decrease in ferredoxin content, which was then correlated with a low Fe concentration, suggesting that Cd induced Fe deficiency (Siedlecka and Baszyński, 1993). In our work, significant changes in Fe content in leaf tissue of ungrafted rootstocks were observed only for IMC67 seedlings. Accordingly, a decrease in F_v/F_m ratio scores was observed in IMC67 seedlings before grafting. A significant decrease in Fe content has also been observed in bean plants (Saidi et al., 2013). On the other hand, a significant increase in Fe content for cacao in grafting combinations having full-sib IMC67 × PA121 progenies as rootstocks was observed too. Despite this, a significant decrease in F_v/F_m was observed in the ICS95 clone when grafted on this rootstock family, whereas no significant changes were observed in this parameter when using CCN51 as scion.

A decrease in the content of total soluble protein was observed in PA121 seedlings before grafting. Significant decreases were also observed in both CCN51 and ICS95 scions 2 months after grafting on the IMC67 OP rootstock family. An effect in total soluble protein content could be explained by Cd-induced protein degradation and an increase in proteolytic activity (Mishra et al., 2006). In some cases, a slight boost in total soluble protein was observed, but it was not significant compared with the control. Slight increases in protein content may be explained by the induction of stress proteins as part of the plant defense system to Cd toxicity (Di Toppi and Gabbrielli, 1999).

Finally, a significant increase in electrolytes leakage was observed in IMC67 and PA121 established in Cd-spiked soil, suggesting a Cd effect on cell membrane integrity. An increase in electrolyte loss has been observed in cucumber (Gonçalves et al., 2007), maize (Ekmekçi et al., 2008), and Bacopa monnieri (Mishra et al., 2006). Heavy metals induce an alteration in the lipid composition of thylakoid membranes. Reactive oxygen species (ROS) induced by Cd results in lipid peroxidation, which implies the degradation of polyunsaturated fatty acids of membrane lipids. The later causes distortion of the lipid bilayer and alters membrane ion channels, resulting in leakage of ions (Mishra et al., 2006). The apparent lack of cell membrane damages in the CCN51 and ICS95 scions, 2 months after grafting, may be due to an insufficient exposure time to Cd as to cause a perceptible effect.

Rootstock Versus Scion Genetic Conflict

This work has enlightened some major trends regarding the complexity at the rootstock–scion interface. First, phenotypic differences due to rootstock effects are more notorious early after grafting (i.e., 2-months), in concert with expectations observed at the ungrafted seedlings. Such differences tend to vanish at older grafted seedlings (i.e., 4 months after grafting). The tendencies observed when computing specific combining abilities across rootstock families are aligned with this conclusion, in the sense that significantly different specific combining abilities among rootstock families were more common at 2 months after grafting than at 4 months. Second, heritability estimates also speak for the dilution of the rootstock effects through genotypes and time. Several of the heritabilities that were calculated in ungrafted families were not significant in the grafted portion of the experiment, which may speak for a genuine physiological communication gap at the graft interface. Alternatively, it could reflect an intrinsic limitation in degrees of freedom due to the complexity of the rootstock × scion factorial design. Still, rootstock-mediated heritability scores above 0.1 were more commonly observed at 2 months after grafting than at 4 months, regardless of the specific rootstock × scion combination. Such instability ultimately suggests an underlying and unavoidable conflict between the two genomes that shape the chimeric grafted organism (Warschefsky et al., 2016; Gautier et al., 2019).

Genetic conflict is known to be pervasive at multiple nested evolutionary scales, for example, among genes, among chromosomes, between chromosomes and cytoplasmic organelles, and between sexes in dioicous species (Ågren and Clark, 2018). These multiple scenarios are evidenced by the recurrent segregation distortion due to gene drives, transposons, and unconventional sex determination systems (Renner and Müller, 2021). However, our study is the first in suggesting an analogous mechanism at the rootstock × scion interface, typically regarded as leading to emergent heterotic properties (Reyes-Herrera et al., 2020). This conflict may be due to underlying additive and combined physiological drivers (Loupit and Cookson, 2020), such as water and nutrients uptake and transport, hormone production and transport, and large-scale movement of molecules during grafting and through time (Rasool et al., 2020).

Foreseeing the Complexity of the Rootstock × Scion Interaction

The novelty of this work lays on the explicit comparison of rootstock families in terms of Cd uptake. However, after quantifying an additive component such as the inheritance of rootstock effects (i.e., rootstock-mediated genetic variance) across recombinant cacao saplings, a next step is to consider more thoroughly the complexity of the rootstock × scion interaction. After all, rootstock metabolites transcend the root system and could reach the grafted scion, i.e., rootstock’s additive contribution (Loupit and Cookson, 2020; Rasool et al., 2020), which in turn may have contrasting consequences on rootstock traits, i.e., scion’s additive contribution (Shu et al., 2017). These concurrent effects would ultimately feedback an emergent rootstock × scion interaction, a statistical interaction in the strict sense.

To be able to accurately estimate the rootstock × scion component, it will be necessary to further validate the rootstock-mediated pedigree-estimated heritability scores via controlled experiments across an expanded panel of clonal scion genotypes grafted on rootstocks with a more continuous gradient of pairwise relatedness values (Reyes-Herrera et al., 2020). Specifically, future experimental assessments might rely on factorial designs of diverse clonal scions grafted on clonal rootstocks. This would allow to reduce recombination uncertainty and optimize statistical power to estimate the interaction term.

From an analytical point of view, the A-BLUP model implemented here would be capable to condition pedigree-based rootstock-mediated heritability scores as a function of the scion’s pedigree. Alternatively, the rootstock × scion interaction may also be quantified via indirect genetic effect (IGE) models (Bijma, 2010, 2013; Fisher and Mcadam, 2019) at relatively low phenotyping costs. It would also be desirable to extend IGE monitoring through time as a way to validate whether some of the significant rootstock × scion effects may persist in adult trees, even several years after grafting.

Perspectives

A recurrent caveat of pedigree-based heritability estimates concerns the potential fortuitous unbalanced between each phenotypic vector and the A pedigree-based relatedness matrix within the “genetic prediction” model. However, we did not detect this trend for any of the significantly rootstock-inherited traits. On the contrary, after examining for 2 years, as part of a parallel experiment, eight traits in full-sib families obtained by cross pollination of clonal accessions from the Cacao Germplasm at greenhouse conditions and two water regimes, we have been able to identify three candidate families for further testing. We recommend using this promissory dataset as reference (i.e., training) population to calibrate explicit eco-physiological mechanistic models (Loìpez-Hernández and Cortés, 2019) and last-generation machine learning algorithms (Cortés and López-Hernández, 2021; Montesinos-Loìpez et al., 2021) as innovative alternatives beyond A-BLUP models (Guevara-Escudero et al., 2021). As part of this task, we envision the following pipeline: (1) developing explicit eco-physiological indices for cacao targeting neo-tropical localities, (2) calibrating last-generation predictive breeding models aiming to forecast such indices based on extensive genealogical information, (3) extending the previous models to account for the complexity of the rootstock-scion interaction (i.e., G × S × E term, as expanded in the previous section), and (4) validating the corresponding predictions across seed orchards and cacao saplings at local nurseries to leverage natural variation for early selection (i.e., before grafting) of low Cd uptake. Ultimately, this combined strategy promises speeding up breeding of polygenic trait variation in a perennial tree crop, while accounting for the interaction of multiple genotypes at the rootstock–scion interface. In parallel, high throughput genotyping of the rootstock families will enable a more accurate description of the underlying genetic architecture.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Author Contributions

CR-M conceived the original sampling and experiment and compiled the datasets. JF-P, CH-V, and CR-M collected phenotypic data at the greenhouses. JF-P and CH-V carried out lab work procedures. CH-V, JF-P, and AC prepared input datasets for statistical softwares and carried out data analyses. CR-M and AC drafted a first version of this manuscript and edited by the other co-authors. CH-V, JF-P, AC, MM-D-T, CR-M, and VB interpreted results, contributed to the manuscript, and approved the submitted version.

Funding

This research was funded by USAID/USDA Cacao for Peace initiative (Collaboration agreement No. 1939, project 1001136). JF-P received financial support to carry out her master’s studies from AGROSAVIA, which is financed with public resources from the Ministry of Agriculture and Rural Development from Colombia (MADR). AGROSAVIA also financed this publication.

Author Disclaimer

The findings and conclusion in this publication are those of the authors and should not be construed to represent any official USDA or United States Government determination or policy.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We would like to thank the National Germplasm Banks (SBGNAA) for supplying plant material of the cacao genotypes included in this study. We also would like to acknowledge technical support from Diego Fernando Morales Villegas AGROSAVIA, Palmira. We extend our thank to Marshall Elson for commenting an early version of this manuscript, as well as to the reviewers and the editor for their insightful recommendations.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.777842/full#supplementary-material

Supplementary Figure 1 | Pearson correlations coefficients among examined traits from half-sib and full-sib seedling rootstocks (A) before grafting, (B) two, and (C) four months after grafting. Estimates were gathered in two OP half-sib families (IMC67 and PA121), and two full-sib rootstock progenies obtained from the crossing between IMC67 × PA12 and its reciprocal cross, grafted with the ICS95 and CCN51 clonal scions under two Cd treatments. Correlation estimates are presented above the diagonal and below diagonal circles are sized and coloured accordingly.

Supplementary Figure 2 | Heat plots of A pedigree relatedness matrices from two OP half-sib families (IMC67 and PA121), and two full-sib rootstock families derived from controlled crosses between IMC67 × PA121 and its reciprocal cross, inputed in a “genetic prediction” additive mixed linear model, according to de los Campos et al. (2009). Heat plots are depicted for verification purposes (A) before grafting, and (B) two and (C) four months after grafting with the ICS95 and CCN51 clonal scions under two Cd treatments.

Supplementary Table 1 | Concentration of mineral nutrients in leaf tissue of ungrafted rootstocks five months after cadmium treatment. Different lowercase letters indicate significant differences in mineral elements concentration between ungrafted seedling families at the same Cd concentration in soil. Different capital letters indicate significant differences in mineral elements concentration of the same ungrafted seedling at different Cd concentrations in soil. Results are according to Duncan test (p < 0.05). Data are mean values ± SE (n = 8).

Supplementary Table 2 | Concentration of mineral nutrients in leaf tissue of scions ICS95 and CCN51 two months after grafting onto rootstocks. Different lowercase letters indicate significant differences in mineral elements concentration in leaf tissue of the ICS95 or CCN5 clonal scions grafted on four different seedling rootstock families at the same Cd concentration in soil. Different capital letters indicate significant differences in mineral elements concentration in leaf tissue of ICS95 or CCN51 clonal scions grafted on the same seedling rootstock family at different Cd concentrations in soil. Results are according to Duncan test (p < 0.05). Data are mean values ± SE (n = 8).

Supplementary Table 3 | Specific combining abilities in cacao seedlings from controlled crosses before grafting. Estimates were gathered in two OP half-sib families (IMC67 and PA121), and two full-sib progenies obtained from the crossing between IMC67 × PA121 and its reciprocal cross, under two Cd treatments. Specific combining abilities are indexed to a general mean (intercept under the first column after the p-value), and the corresponding positive or negative contributions (three last columns). Significantly different specific combining abilities are highlighted.

Supplementary Table 4 | Narrow-sense heritability (h2) in cacao seedlings from controlled crosses before grafting. Estimates were gathered in two OP half-sib families (IMC67 and PA121), and two full-sib progenies obtained from the crossing between IMC67 × PA121 and its reciprocal cross, grafted with the ICS95 and CCN51 clonal scions under two Cd treatments. Heritability (h²) and model fits (r) estimates were computed using an A pedigree relatedness matrix (Supplementary Figure 2A) inputted in a “genetic prediction” additive mixed linear model, according to de los Campos et al. (2009). Heritability (h²) scores above 0.1 are highlighted.

Supplementary Table 5 | Specific combining abilities in cacao seedling rootstocks from controlled crosses two-months after grafting. Estimates were gathered in full-sib rootstock progenies obtained from the crossing between IMC67 × PA121, and its reciprocal cross, grafted with the ICS95 and CCN51 clonal scions under two Cd treatments. Specific combining abilities are indexed to a general mean (intercept under the first column after the p-value), and the corresponding positive or negative contributions (three last columns). Significantly different specific combining abilities are highlighted.

Supplementary Table 6 | Specific combining abilities in cacao seedling rootstocks from controlled crosses four-months after grafting. Estimates were gathered in two OP half-sib families (IMC67 and PA121), and two full-sib progenies obtained from the crossing between IMC67 × PA121 and its reciprocal cross, grafted with the ICS95 and CCN51 clonal scions under two Cd treatments. Specific combining abilities are indexed to a general mean (intercept under the first column after the p-value), and the corresponding positive or negative contributions (three last columns). Significantly different specific combining abilities are highlighted.

Supplementary Table 7 | Narrow-sense rootstock-mediated heritability (h2) in cacao seedling rootstocks from controlled-crosses two-months after grafting. Estimates were gathered in two OP half-sib families (IMC67 and PA121), and two full-sib progenies obtained from the crossing between IMC67 × PA121 and its reciprocal cross, grafted with the ICS95 and CCN51 clonal scions under two Cd treatments. Heritability (h²) and model fits (r) estimates were computed using an A pedigree relatedness matrix (Supplementary Figure 2B) inputted in a “genetic prediction” additive mixed linear model, according to de los Campos et al. (2009). Heritability (h²) scores above 0.1 are highlighted.

Supplementary Table 8 | Narrow-sense rootstock-mediated heritability (h2) in cacao seedling rootstocks from controlled-crosses four-months after grafting. Estimates were gathered in two OP half-sib families (IMC67 and PA121), and two full-sib progenies obtained from the crossing between IMC67 × PA121 and its reciprocal cross, grafted with the ICS95 and CCN51 clonal scions under two Cd treatments. Heritability (h²) and model fits (r) estimates were computed using an A pedigree relatedness matrix (Supplementary Figure 2C) inputted in a “genetic prediction” additive mixed linear model, according to de los Campos et al. (2009). None h² estimate was above 0.1.

References

Ågren, J. A., and Clark, A. G. (2018). Selfish genetic elements. PLoS Genet. 14:e1007700. doi: 10.1371/journal.pgen.1007700

PubMed Abstract | CrossRef Full Text | Google Scholar

Ahmad, P., Nabi, G., and Ashraf, M. (2011). Cadmium-induced oxidative damage in mustard [Brassica juncea (L.) Czern. & Coss.] plants can be alleviated by salicylic acid. South Afr. J. Bot. 77, 36–44. doi: 10.1016/j.sajb.2010.05.003

CrossRef Full Text | Google Scholar

Arao, T., Ishikawa, S., Murakami, M., Abe, K., Maejima, Y., and Makino, T. (2010). Heavy metal contamination of agricultural soil and counter measures in Japan. Paddy and water Environ. 8, 247–257. doi: 10.1007/s10333-010-0205-7

CrossRef Full Text | Google Scholar

Arao, T., Takeda, H., and Nishihara, E. (2008). Reduction of cadmium translocation from roots to shoots in eggplant (Solanum melongena) by grafting onto Solanum torvum rootstock. Soil Sci. Plant Nutr. 54, 555–559. doi: 10.1111/j.1747-0765.2008.00269.x

CrossRef Full Text | Google Scholar

Arenas, S., Cortés, A. J., Mastretta, A., and Jaramillo, J. P. (2021). Evaluating the accuracy of genomic prediction for the management and conservation of relictual natural tree populations. Tree Genet. Genomes 17, 1–19. doi: 10.1007/s11295-020-01489-1

CrossRef Full Text | Google Scholar

Arévalo-Gardini, E., Arévalo-Hernández, C. O., Baligar, V. C., and He, Z. L. (2017). Heavy metal accumulation in leaves and beans of cacao (Theobroma cacao L.) in major cacao growing regions in Peru. Sci. Total Environ. 605, 792–800. doi: 10.1016/j.scitotenv.2017.06.122

PubMed Abstract | CrossRef Full Text | Google Scholar

Argüello, D., Chavez, E., Lauryssen, F., Vanderschueren, R., Smolders, E., and Montalvo, D. (2019). Soil properties and agronomic factors affecting cadmium concentrations in cacao beans: a nationwide survey in Ecuador. Sci. Total Environ. 649, 120–127. doi: 10.1016/j.scitotenv.2018.08.292

PubMed Abstract | CrossRef Full Text | Google Scholar

Asman, A., bin Purung, M. H., Lambert, S., Amiruddin, A., and Rosmana, A. (2021). Effect of rootstock and scion on resistance of cocoa clones to vascular streak dieback caused by Ceratobasidium theobromae. Ann. Agricult. Sci. 66, 25–30. doi: 10.1016/j.aoas.2021.02.005

CrossRef Full Text | Google Scholar

Barceló, J., Vázquez, M. D., and Poschenrieder, C. H. (1988). Structural and ultrastructural disorders in cadmium-treated bush bean plants (Phaseolus vulgaris L.). New Phytol. 108, 37–49. doi: 10.1111/j.1469-81371988.tb00202.x

CrossRef Full Text | Google Scholar

Barraza, F., Schreck, E., Lévêque, T., Uzu, G., López, F., Ruales, J., et al. (2017). Cadmium bioaccumulation and gastric bioaccessibility in cacao: a field study in areas impacted by oil activities in Ecuador. Environ. Pollut. 229, 950–963. doi: 10.1016/j.envpol.2017.07.080

PubMed Abstract | CrossRef Full Text | Google Scholar

Bijma, P. (2010). Estimating indirect genetic effects: precision of estimates and optimum designs. Genetics 186, 1013–1028. doi: 10.1534/genetics.110.120493

PubMed Abstract | CrossRef Full Text | Google Scholar

Bijma, P. (2013). The quantitative genetics of indirect genetic effects: a selective review of modelling issues. Heredity 112, 61–69. doi: 10.1038/hdy.2013.15

PubMed Abstract | CrossRef Full Text | Google Scholar

Borek, M., Ba̧czek-Kwinta, R., and Rapacz, M. (2013). “Chlorophyll fluorescence imaging of cadmium-treated white cabbage plants,” in Proceedings of the 16th International Conference on Heavy Metals in the Environment E3S Web of Conferences, Vol. 1, (Les Ulis: EDP Sciences), 39004.

Google Scholar

Boughalleb, N., Tarchoun, N., El-Mbarki, A., and El-Mahjoub, M. (2007). Resistance evaluation of nine cucurbit rootstocks and grafted watermelon (Citrullus lanatus L.) varieties against Fusarium wilt and Fusarium crown and root rot. J. Plant Sci. 2, 102–107. doi: 10.3923/jps.2007.102.107

CrossRef Full Text | Google Scholar

Bradford, M. M. (1976). A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal. Biochem. 72, 248–254. doi: 10.1006/abio.1976.9999

PubMed Abstract | CrossRef Full Text | Google Scholar

Bradstreet, R. B. (1954). Kjeldahl method for organic nitrogen. Anal. Chem. 26, 185–187.

Google Scholar

Cardinal, A. B. B., Gonçalves, P. D. S., and Martins, A. L. M. (2007). Stock-scion interactions on growth and rubber yield of Hevea brasiliensis. Sci. Agric. 64, 235–240. doi: 10.1590/S0103-90162007000300004

CrossRef Full Text | Google Scholar

Castro, A. V., de Almeida, A. A. F., Pirovani, C. P., Reis, G. S., Almeida, N. M., and Mangabeira, P. A. (2015). Morphological, biochemical, molecular and ultrastructural changes induced by Cd toxicity in seedlings of Theobroma cacao L. Ecotoxicol. Environ. Safety 115, 174–186. doi: 10.1016/j.ecoenv.2015.02.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Chavez, E., He, Z. L., Stoffella, P. J., Mylavarapu, R. S., Li, Y. C., Moyano, B., et al. (2015). Concentration of cadmium in cacao beans and its relationship with soil cadmium in southern Ecuador. Sci. Total Environ. 533, 205–214. doi: 10.1016/j.scitotenv.2015.06.106

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Z. Q., Baison, J., Pan, J., Karlsson, B., Andersson, B., Westin, J., et al. (2018). Accuracy of genomic selection for growth and wood quality traits in two control-pollinated progeny trials using exome capture as the genotyping platform in Norway spruce. BMC Genomics 19:946.

Google Scholar

Ciećko, Z., Kalembasa, S., Wyszkowski, M., and Rolka, E. (2005). The magnesium content in plants in soil contaminated with cadmium. Polish J. Environ. Stud. 14, 365–370.

Google Scholar

Cortés, A. J., and López-Hernández, F. (2021). Harnessing crop wild diversity for climate adaptation. Genes 12:783.

Google Scholar

Cuevas, J., Crossa, J., Soberanis, V., Peìrez-Elizalde, S., Peìrez-Rodriìguez, P., de los Campos, G., et al. (2016). Genomic prediction with genotype x environment interaction kernel regres- sion models. Plant Genome 9, 1–20.

Google Scholar

Davis, A. R., and Perkins-Veazie, P. (2005). Rootstock effects on plant vigour and watermelon fruit quality. Cucurbit Genet. Coop. Rep. 28, 39–42.

Google Scholar

de los Campos, G., Naya, H., Gianola, D., Crossa, J., Legarra, A., Manfredi, E., et al. (2009). Predicting quantitative traits with regression models for dense molecular markers and pedigree. Genetics 182, 375–385. doi: 10.1534/genetics.109.101501

PubMed Abstract | CrossRef Full Text | Google Scholar

de los Campos, G., Sorensen, D., and Gianola, D. (2015). Genomic heritability: what is it? PLoS Genet. 11:e1005048. doi: 10.1371/journal.pgen.1005048

PubMed Abstract | CrossRef Full Text | Google Scholar

Di Toppi, L. S., and Gabbrielli, R. (1999). Response to cadmium in higher plants. Environ. Exp. Bot. 41, 105–130. doi: 10.1016/S0098-8472(98)00058-6

CrossRef Full Text | Google Scholar

Dias, M. C., Monteiro, C., Moutinho-Pereira, J., Correia, C., Gonçalves, B., and Santos, C. (2012). Cadmium toxicity affects photosynthesis and plant growth at different levels. Acta Physiol. Plant. 35, 1281–1289. doi: 10.1007/s11738-012-1167-8

CrossRef Full Text | Google Scholar

Ekmekçi, Y., Tanyolac, D., and Ayhan, B. (2008). Effects of cadmium on antioxidant enzyme and photosynthetic activities in leaves of two maize cultivars. J. Plant Physiol. 165, 600–611. doi: 10.1016/j.jplph.2007.01.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Eltayb, M. T. A., Magid, T. D. A., Ibrahim, A. A., and Dirar, A. M. A. (2014). Effect of grafting (rootstock) on morphological changes of scions in some Acacia species. J. Forest Prod. Ind. Corpus. 3, 27–36.

Google Scholar

Engbersen, N., Gramlich, A., Lopez, M., Schwarz, G., Hattendorf, B., Gutierrez, O., et al. (2019). Cadmium accumulation and allocation in different cacao cultivars. Sci. Total Environ. 678, 660–670. doi: 10.1016/j.scitotenv.2019.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

European Commission (2014). Commission Regulation (EC) No 488/2014 of 12 May 2014 amending Regulation (EC) No 1881/2006 as regards maximum levels of cadmium in foodstuffs. Official J. Eur. Union 138:75.

Google Scholar

Farooq, M. A., Ali, S., Hameed, A., Bharwana, S. A., Rizwan, M., Ishaque, W., et al. (2016). Cadmium stress in cotton seedlings: physiological, photosynthesis and oxidative damages alleviated by glycinebetaine. South Afr. J. Bot. 104, 61–68. doi: 10.1016/j.sajb.2015.11.006

CrossRef Full Text | Google Scholar

Fisher, D. N., and Mcadam, A. G. (2019). Indirect genetic e ects clarify how traits can evolve even when fitness does not. Evol. Lett. 3, 4–14. doi: 10.1002/evl3.98

CrossRef Full Text | Google Scholar

Gautier, A. T., Chambaud, C., Brocard, L., Ollat, N., Gambetta, G. A., Delrot, S., et al. (2019). Merging genotypes: graft union formation and scion-rootstock interactions. J. Exp. Bot. 70, 747–755. doi: 10.1093/jxb/ery422

PubMed Abstract | CrossRef Full Text | Google Scholar

Ge, W., Jiao, Y., Zou, J., Jiang, W., and Liu, D. (2015). Ultrastructural and photosynthetic response of populus 107 leaves to cadmium stress. Polish J. Environ. Stud. 24, 519–527. doi: 10.15244/pjoes/27814

CrossRef Full Text | Google Scholar

Guevara-Escudero, M., Osorio, A. N., and Cortés, A. J. (2021). Integrative breeding for biotic resistance in forest trees. Plants 10:2022.

Google Scholar

Gonçalves, J. F., Becker, A. G., Cargnelutti, D., Tabaldi, L. A., Pereira, L. B., Battisti, V., et al. (2007). Cadmium toxicity causes oxidative stress and induces response of the antioxidant system in cucumber seedlings. Braz. J. Plant Physiol. 19, 223–232. doi: 10.1590/S1677-04202007000300006

CrossRef Full Text | Google Scholar

Gutierrez, A. P., Matika, O., Bean, T. P., and Houston, R. D. (2018). Genomic selec- tion for growth traits in Pacific oyster (Crassostrea gigas): potential of low-density marker panels for breeding value prediction. Front. Genet. 9:391. doi: 10.3389/fgene.2018.00391

PubMed Abstract | CrossRef Full Text | Google Scholar

Hédiji, H., Djebali, W., Belkadhi, A., Cabasson, C., Moing, A., Rolin, D., et al. (2015). Impact of long-term cadmium exposure on mineral content of Solanum lycopersicum plants: consequences on fruit production. South Afr. J. Bot. 97, 176–181. doi: 10.1016/j.sajb.2015.01.010

CrossRef Full Text | Google Scholar

Isaac, R. A., and Kerber, J. D. (1971). “Atomic absorption and flame photometry: techniques and uses in soil, plant, and water analysis,” in Instrumental Methods for Analysis of Soils and Plant Tissue, ed L. M. Walsh. doi: 10.2136/1971.instrumentalmethods.c2

CrossRef Full Text | Google Scholar

Isele, E., Breen, M., and Galanti, R. (2020). Grafting Cacao. College of Tropical Agriculture and Human Resources (CTAHR).

Google Scholar

Jiao, Y., Ge, W., Qin, R., Sun, B., Jiang, W., and Liu, D. (2012). Influence of cadmium stress on growth, ultrastructure and antioxidative enzymes in Populus 2025. Fresenius Environ. Bull. 21, 1375–1384.

Google Scholar

Jiao, Y., Zou, J., Ge, W., Jiang, W., and Liu, D. (2015). Physiological and ultrastructural effects of cadmium on Poplar (Populus x euramericana) leaves. Baltic Forestry 21, 106–113.

Google Scholar

Küpper, H., Kupper, F., and Spiller, M. (1996). Environmental relevance of heavy metal-substituted chlorophylls using the example of water plants. J. Exp. Bot. 47, 259–266. doi: 10.1093/jxb/47.2.259

CrossRef Full Text | Google Scholar

Küpper, H., Küpper, F., and Spiller, M. (1998). In situ detection of heavy metal substituted chlorophylls in water plants. Photosynth. Res. 58, 123–133.

Google Scholar

Lewis, C., Lennon, A. M., Eudoxie, G., and Umaharan, P. (2018). Genetic variation in bioaccumulation and partitioning of cadmium in Theobroma cacao L. Sci. Total Environ. 640, 696–703. doi: 10.1016/j.scitotenv.2018.05.365

PubMed Abstract | CrossRef Full Text | Google Scholar

Lichtenthaler, H. K. (1987). Chlorophylls and carotenoids: pigments of photosynthetic biomembranes. Methods Enzymol. 148, 350–382.

Google Scholar

Liu, D., Jiang, W., and Gao, X. (2003). Effects of cadmium on root growth, cell division and nucleoli in root tip cells of garlic. Biol. Plant. 47, 79–83. doi: 10.1023/A:1027384932338

CrossRef Full Text | Google Scholar

Liu, S., Li, H., Lv, X., Ahammed, G. J., Xia, X., Zhou, J., et al. (2016). Grafting cucumber onto luffa improves drought tolerance by increasing ABA biosynthesis and sensitivity. Sci. Rep. 6, 1–14. doi: 10.1038/srep20212

PubMed Abstract | CrossRef Full Text | Google Scholar

Loìpez-Hernández, F., and Cortés, A. J. (2019). Last-generation genome–environment associations reveal the genetic basis of heat tolerance in common bean (Phaseolus vulgaris L.). Front. Genet. 10:954. doi: 10.3389/fgene.2019.00954

PubMed Abstract | CrossRef Full Text | Google Scholar

Loupit, G., and Cookson, S. J. (2020). Identifying molecular markers of successful graft union formation and compatibility. Front. Plant Sci. 11:352. doi: 10.3389/fpls.2020.610352

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, Z., Zhang, Z., Su, Y., Liu, C., and Shi, G. (2013). Cultivar variation in morphological response of peanut roots to cadmium stress and its relation to cadmium accumulation. Ecotoxicol. Environ. Safety 91, 147–155. doi: 10.1016/j.ecoenv.2013.01.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Lux, A., Martinka, M., Vaculík, M., and White, P. J. (2011). Root responses to cadmium in the rhizosphere: a review. J. Exp. Bot. 62, 21–37. doi: 10.1093/jxb/erq281

PubMed Abstract | CrossRef Full Text | Google Scholar

Melgarejo, L. M., Romero, M., Hernández, S., Jaime, S. M. E., Suárez, D., and Pérez, W. (2010). Experimentos en Fisiología Vegetal Lab Físiol Bioquímica Veg. Bogotá: Universidad Nacional de Colombia. Available online at: https://www.uv.mx/personal/tcarmona/files/2019/02/Melgarejo-2010.pdf

Google Scholar

Mishra, S., Srivastava, S., Tripathi, R. D., Govindarajan, R., Kuriakose, S. V., and Prasad, M. N. V. (2006). Phytochelatin synthesis and response of antioxidants during cadmium stress in Bacopa monnieri L? Plant Physiol. Biochem. 44, 25–37. doi: 10.1016/j.plaphy.2006.01.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Montesinos-Loìpez, O. A., Montesinos-Loìpez, A., Peìrez-Rodriìguez, P., Barroìn-Loìpez, J. A., Martini, J. W. R., and Fajardo-Flores, S. B. (2021). A review of deep learning applications for genomic selection. BMC Genomics 22:19.

Google Scholar

Moore, R. E., Ullah, I., de Oliveira, V. H., Hammond, S. J., Strekopytov, S., Tibbett, M., et al. (2020). Cadmium isotope fractionation reveals genetic variation in Cd uptake and translocation by Theobroma cacao and role of natural resistance-associated macrophage protein 5 and heavy metal ATPase-family transporters. Horticult. Res. 7, 1–11. doi: 10.1038/s41438-020-0292-6

CrossRef Full Text | Google Scholar

Müller, B. S. F., Neves, L. G., de Almeida Filho, J. E., Resende, M. F. R., Munϸoz, P. R., et al. (2017). Genomic prediction in contrast to a genome-wide association study in explaining heritable variation of complex growth traits in breeding populations of Eucalyptus. BMC Genomics 18:1–17. doi: 10.1186/s12864-017-3920-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Nguyen, N. T., McInturf, S. A., and Mendoza-Cózatl, D. G. (2016). Hydroponics: a versatile system to study nutrient allocation and plant responses to nutrient availability and exposure to toxic elements. J. Visual. Exp. JoVE 113:54317. doi: 10.3791/54317

PubMed Abstract | CrossRef Full Text | Google Scholar

Palencia Calderon, G. E., Gómez Santos, R., and Mejia, L. A. (2007). Patrones Para Cacao: Cartilla (No. Doc. 22050) CO-BAC, Bogotá). Bucaramanga: Editorial AGROSAVIA.

Google Scholar

Pereira de Araújo, R., Furtado, de Almeida, A. A., Pereira, L. S., Mangabeira, P. A., Souza, J. O., et al. (2017). Photosynthetic, antioxidative, molecular and ultrastructural responses of young cacao plants to Cd toxicity in the soil. Ecotoxicol. Environ. Safety 144, 148–157. doi: 10.1016/j.ecoenv.2017.06.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Pérez-Rodríguez, P., and de los Campos, G. (2014). Genome-wide regression and prediction with the Bglr statistical package. Genetics 198, 483–495. doi: 10.1534/genetics

CrossRef Full Text | Google Scholar

Pietrini, F., Iannelli, M. A., Pasqualini, S., and Massacci, A. (2003). Interaction of cadmium with glutathione and photosynthesis in developing leaves and chloroplasts of Phragmites australis (Cav.) Trin. ex Steudel. Plant Physiol. 133, 829–837. doi: 10.1104/pp.103.026518

PubMed Abstract | CrossRef Full Text | Google Scholar

Rahman, A., Mostofa, M. G., Nahar, K., Hasanuzzaman, M., and Fujita, M. (2016). Exogenous calcium alleviates cadmium-induced oxidative stress in rice (Oryza sativa L.) seedlings by regulating the antioxidant defense and glyoxalase systems. Braz. J. Bot. 39, 393–407. doi: 10.1016/j.crvi.2016.08.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Rai, R., Agrawal, M., and Agrawal, S. B. (2016). “Impact of heavy metals on physiological processes of plants: with special reference to photosynthetic system,” in Plant Responses to Xenobiotics, (Singapore: Springer), 127–140. doi: 10.1007/978-981-10-2860-1_6

CrossRef Full Text | Google Scholar

Rasool, A., Mansoor, S., Bhat, K. M., Hassan, G. I., Baba, T. R., Alyemeni, M. N., et al. (2020). Mechanisms underlying graft union formation and rootstock scion interaction in horticultural plants. Front. Plant Sci. 11:590847. doi: 10.3389/fpls.2020.590847

PubMed Abstract | CrossRef Full Text | Google Scholar

Renner, S. S., and Müller, N. A. (2021). Plant sex chromosomes defy evolutionary models of expanding. Nat. Plants 7, 392–402. doi: 10.1038/s41477-021-00884-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Reyes-Herrera, P. H., Muñoz-Baena, L., Velásquez-Zapata, V., Patiño, L., Delgado-Paz, O. A., Díaz-Diez, C. A., et al. (2020). Inheritance of rootstock effects in avocado (Persea americana Mill.) cv. Hass. Front. Plant Sci. 11:555071. doi: 10.3389/fpls.2020.555071

PubMed Abstract | CrossRef Full Text | Google Scholar

Ribeiro, M. A., de Almeida, A. A. F., Alves, T. F., Gramacho, K. P., Pirovani, C. P., and Valle, R. R. (2016). Rootstock× scion interactions on Theobroma cacao resistance to witches’ broom: photosynthetic, nutritional and antioxidant metabolism responses. Acta Physiol. Plant. 38:73. doi: 10.1007/s11738-016-2095-9

CrossRef Full Text | Google Scholar

Ríos, A. V. (1957). La relación patrón-injerto en la propagación vegetativa del cacao (Theobroma cacao L.). Acta Agron. 7, 223–232.

Google Scholar

Rivard, C. L., and Louws, F. J. (2008). Grafting to manage soilborne diseases in heirloom tomato production. HortScience 43, 2104–2111. doi: 10.21273/HORTSCI.43.7.2104

CrossRef Full Text | Google Scholar

Rouphael, Y., Cardarelli, M., Rea, E., and Colla, G. (2008). Grafting of cucumber as a means to minimize copper toxicity. Environ. Exp. Bot. 63, 49–58. doi: 10.1016/j.envexpbot.2007.10.015

CrossRef Full Text | Google Scholar

Saidi, I., Ayouni, M., Dhieb, A., Chtourou, Y., Chaïbi, W., and Djebali, W. (2013). Oxidative damages induced by short-term exposure to cadmium in bean plants: protective role of salicylic acid. South Afr. J. Bot. 85, 32–38. doi: 10.1016/j.sajb.2012.12.002

CrossRef Full Text | Google Scholar

Sandalio, L. M., Dalurzo, H. C., Gomez, M., Romero-Puertas, M. C., and Del Rio, L. A. (2001). Cadmium-induced changes in the growth and oxidative metabolism of pea plants. J. Exp. Bot. 52, 2115–2126. doi: 10.1093/jexbot/52.364.2115

PubMed Abstract | CrossRef Full Text | Google Scholar

Savvas, D., Colla, G., Rouphael, Y., and Schwarz, D. (2010). Amelioration of heavy metal and nutrient stress in fruit vegetables by grafting. Sci. Horticult. 127, 156–161. doi: 10.1016/j.scienta.2010.09.011

CrossRef Full Text | Google Scholar

Shu, B., Liu, L., Jue, D., Wang, Y., Wei, Y., and Shi, S. (2017). Effects of avocado (Persea americana Mill.) scion on arbuscular mycorrhizal and root hair development in rootstock. Arch. Agron. Soil Sci. 63, 1951–1962.

Google Scholar

Siedlecka, A., and Baszyński, T. (1993). Inhibition of electron flow around photosystem I in chloroplasts of Cd-treated maize plants is due to Cd-induced iron deficiency. Physiol. Plant. 87, 199–202.

Google Scholar

Sigfridsson, K. G., Bernát, G., Mamedov, F., and Styring, S. (2004). Molecular interference of Cd2+ with Photosystem II. Biochim. Biophys. Acta (BBA)-Bioenergetics 1659, 19–31. doi: 10.1016/j.bbabio.2004.07.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Singh, P. K., and Tewari, R. K. (2003). Cadmium toxicity induced changes in plant water relations and oxidative metabolism of Brassica juncea L. plants. J. Environ. Biol. 1, 107–112.

Google Scholar

Song, Y., Jin, L., and Wang, X. (2017). Cadmium absorption and transportation pathways in plants. Int. J. Phytoremed. 19, 133–141. doi: 10.1080/15226514.2016.1207598

PubMed Abstract | CrossRef Full Text | Google Scholar

Spanò, R., Ferrara, M., Gallitelli, D., and Mascia, T. (2020). The role of grafting in the resistance of tomato to viruses. Plants 9:1042. doi: 10.3390/plants9081042

PubMed Abstract | CrossRef Full Text | Google Scholar

Suárez-Hernández, ÁM., Vázquez-Angulo, J. C., Grimaldo-Juárez, O., Duran, C. C., González-Mendoza, D., Bazante-González, I., et al. (2019). Production and quality of grafted watermelon in saline soil. Horticult. Brasil. 37, 215–220. doi: 10.1590/S0102-053620190213

CrossRef Full Text | Google Scholar

Sugiyama, M., Ae, N., and Arao, T. (2007). Role of roots in differences in seed cadmium concentration among soybean cultivars—proof by grafting experiment. Plant Soil 295, 1–11. doi: 10.1007/s11104-007-9252-2

CrossRef Full Text | Google Scholar

Ullah, I., Wang, Y., Eide, D. J., and Dunwell, J. M. (2018). Evolution, and functional analysis of natural resistance-associated macrophage proteins (NRAMPs) from Theobroma cacao and their role in cadmium accumulation. Sci. Rep. 8, 1–15. doi: 10.1038/s41598-018-32819-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Vanderschueren, R., Montalvo, D., De Ketelaere, B., Delcour, J. A., and Smolders, E. (2019). The elemental composition of chocolates is related to cacao content and origin: a multi-element fingerprinting analysis of single origin chocolates. J. Food Composition Anal. 83:103277. doi: 10.1016/j.jfca.2019.103277

CrossRef Full Text | Google Scholar

Vanderschueren, R., Argüello, D., Blommaert, H., Montalvo, D., Barraza, F., Maurice, L., et al. (2021). Mitigating the level of cadmium in cacao products: reviewing the transfer of cadmium from soil to chocolate bar. Sci. Total Environ. 146779. doi: 10.1016/j.scitotenv.2021.146779

CrossRef Full Text | Google Scholar

Warschefsky, E. J., Klein, L. L., Frank, M. H., Chitwood, D. H., Londo, J. P., Von Wettberg, E. J. B., et al. (2016). Rootstocks: diversity, domestication, and impacts on shoot phenotypes. Trends Plant Sci. 21, 418–437. doi: 10.1016/j.tplants.2015.11.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, S., Li, Y., Zhan, J., Wang, S., and Zhu, J. (2012). Tolerant mechanisms of Rorippa globosa (Turcz.) Thell. hyperaccumulating Cd explored from root morphology. Bioresource Technol. 118, 455–459. doi: 10.1016/j.biortech.2012.05.049

PubMed Abstract | CrossRef Full Text | Google Scholar

Yanyan, Y. A. N., Shuoshuo, W. A. N. G., Min, W. E. I., Biao, G. O. N. G., and Qinghua, S. H. I. (2018). Effect of different rootstocks on the salt stress tolerance in watermelon seedlings. Horticult. Plant J. 4, 239–249. doi: 10.1016/j.hpj.2018.08.003

CrossRef Full Text | Google Scholar

Yin, J. P. T. (2004). Rootstock effects on cocoa in Sabah, Malaysia. Exp. Agricult. 40, 445–452. doi: 10.1017/S0014479704002108

CrossRef Full Text | Google Scholar

Zhang, F., Shi, W., Jin, Z., and Shen, Z. (2002). Response of antioxidative enzymes in cucumber chloroplasts to cadmium toxicity. J. Plant Nutr. 26, 1779–1788. doi: 10.1081/PLN-120023282

CrossRef Full Text | Google Scholar

Zhang, H., Yin, L., Wang, M., Yuan, X., and Liu, X. (2019). Factors affecting the accuracy of genomic selection for agricultural economic traits in maize, cattle, and pig populations. Front. Genet. 10:189. doi: 10.3389/fgene.2019.00189

PubMed Abstract | CrossRef Full Text | Google Scholar

Zornoza, P., Vázquez, S., Esteban, E., Fernández-Pascual, M., and Carpena, R. (2002). Cadmium-stress in nodulated white lupin: strategies to avoid toxicity. Plant Physiol. Biochem. 40, 1003–1009. doi: 10.1016/S0981-9428(02)01464-X

CrossRef Full Text | Google Scholar