<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Agron.</journal-id>
<journal-title>Frontiers in Agronomy</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Agron.</abbrev-journal-title>
<issn pub-type="epub">2673-3218</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fagro.2021.754646</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Agronomy</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Interactive and Dynamic Effects of Rootstock and Rhizobiome on Scion Nutrition in Cacao Seedlings</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name><surname>Schmidt</surname> <given-names>Jennifer E.</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x0002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/314527/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>DuVal</surname> <given-names>Ashley</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/461782/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Puig</surname> <given-names>Alina</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/916130/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Tempeleu</surname> <given-names>Alexandra</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1532705/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Crow</surname> <given-names>Taylor</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1517992/overview"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Breeding and Genetics, Cocoa Plant Sciences, Mars Wrigley</institution>, <addr-line>Davis, CA</addr-line>, <country>United States</country></aff>
<aff id="aff2"><sup>2</sup><institution>Subtropical Horticulture Research Station, United States Department of Agriculture (USDA), Agricultural Research Service (ARS)</institution>, <addr-line>Miami, FL</addr-line>, <country>United States</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Mohamed Hijri, Universit&#x000E9; de Montr&#x000E9;al, Canada</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Durgesh Kumar Jaiswal, Banaras Hindu University, India; Dirceu Mattos Jr., Instituto Agron&#x000F4;mico de Campinas (IAC), Brazil</p></fn>
<corresp id="c001">&#x0002A;Correspondence: Jennifer E. Schmidt <email>jennifer.schmidt&#x00040;effem.com</email></corresp>
<fn fn-type="other" id="fn001"><p>This article was submitted to Plant-Soil Interactions, a section of the journal Frontiers in Agronomy</p></fn></author-notes>
<pub-date pub-type="epub">
<day>18</day>
<month>11</month>
<year>2021</year>
</pub-date>
<pub-date pub-type="collection">
<year>2021</year>
</pub-date>
<volume>3</volume>
<elocation-id>754646</elocation-id>
<history>
<date date-type="received">
<day>06</day>
<month>08</month>
<year>2021</year>
</date>
<date date-type="accepted">
<day>18</day>
<month>10</month>
<year>2021</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2021 Schmidt, DuVal, Puig, Tempeleu and Crow.</copyright-statement>
<copyright-year>2021</copyright-year>
<copyright-holder>Schmidt, DuVal, Puig, Tempeleu and Crow</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license>
</permissions>
<abstract><p>Perennial agroecosystems often seek to optimize productivity by breeding nutrient-efficient, disease-resistant rootstocks. In cacao (<italic>Theobroma cacao</italic> L.), however, rootstock selection has traditionally relied on locally available open pollinated populations with limited data on performance. Furthermore, rootstock associations with the rhizobiome, or rhizosphere microbiome, have been neglected. Better understanding of rootstock and scion effects on cacao-specific traits, particularly those involved in root-microbe interactions and nutrient acquisition, could contribute to more efficient rootstock selection and breeding. A rootstock-scion interaction study was conducted using three scion genotypes and eight rootstock populations under greenhouse conditions to better understand the relationships among rootstock and scion identities, soil fertility, and rhizobiome composition and the impacts of these factors on plant uptake of macro- and micronutrients. We show that rootstock genotype has a stronger influence than scion on nutrient uptake, bacterial and fungal diversity, and rhizobiome composition, and that the relative contributions of rootstock and scion genotype to foliar nutrient status are dynamic over time. Correlation analysis and stepwise regression revealed complex relationships of soil physicochemical parameters and the rhizobiome to plant nutrition and emphasized strong impacts of microbial diversity and composition on specific nutrients. Linear discriminant analysis effect size estimation identified rootstock-responsive taxa potentially related to plant nutrition. This study highlights the importance of considering root-associated microbial communities as a factor in cacao rootstock breeding and the need for further investigation into mechanisms underlying nutrient acquisition and microbial interactions in grafted plants.</p></abstract>
<kwd-group>
<kwd>rootstock</kwd>
<kwd><italic>Theobroma cacao</italic></kwd>
<kwd>rhizobiome</kwd>
<kwd>nutrient uptake</kwd>
<kwd>rhizosphere</kwd>
<kwd>microbiome</kwd>
<kwd>grafting</kwd>
</kwd-group>
<counts>
<fig-count count="5"/>
<table-count count="2"/>
<equation-count count="0"/>
<ref-count count="56"/>
<page-count count="13"/>
<word-count count="7602"/>
</counts>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<title>Introduction</title>
<p>Grafting is a common practice in perennial crops to combine desirable combinations of rootstock and scion traits. Although grafting has been called a &#x0201C;surgical alternative to breeding&#x0201D; for its utility in creating desirable combinations of traits (Albacete et al., <xref ref-type="bibr" rid="B3">2015</xref>), it is better understood as a synergistic counterpart to breeding. Physical attachment of rootstocks bred for disease resistance and nutrient use efficiency to scions bred for fruit quality and quantity enables divergent sets of traits to be combined without incurring linkage drag and yield penalties. However, breeding programs have not yet optimized both rootstocks and scions in many grafted perennial crops, including <italic>Theobroma cacao</italic>.</p>
<p>Valued for its pods, which are used in chocolate, cocoa butter, and other economically important commodities, cacao is grown in lowland tropical regions such as West Africa, Latin America, and Southeast Asia. Modern cacao breeding programs have only existed since the 1930s, despite the popularity of this crop, which translates to a limited number of generations given the timeframe required for breeding perennials (DuVal et al., <xref ref-type="bibr" rid="B20">2017</xref>). To date, these programs have focused primarily on generating elite scions, leaving the development of nutrient-efficient, stress-tolerant, disease-resistant rootstocks as a significant research gap. One reason may be that at least two of the three primary breeding objectives (bean quality, yield, and disease resistance) are controlled by scions, and soil-borne pathogens that can be controlled with resistant rootstocks are not as significant as foliar and pod pathogens in cacao (Bailey and Meinhardt, <xref ref-type="bibr" rid="B6">2016</xref>). Another reason for the relative disinterest in rootstocks may be the dominance of seedlings over vegetative propagation methods in West Africa, which produces 70% of the world&#x00027;s cacao (Yin, <xref ref-type="bibr" rid="B55">2004</xref>). Difficulties inherent in phenotyping belowground traits may also contribute. Yet rootstock breeding could play a critical role in the future of cacao. Meeting rising consumer demand without conversion of additional land to cacao production will require increasing resource use efficiency, especially at the root-soil interface.</p>
<p>Studies of nutrient use efficiency in cacao rootstocks are rare, but have shown potential that should be exploited to develop rootstocks that effectively capture soil nutrients, even in depleted soils. A comparison of two cacao rootstocks found differences in nitrogen use efficiency and growth response that were related to root dry weight (Ribeiro et al., <xref ref-type="bibr" rid="B45">2008</xref>). Evidence from other tropical perennials such as mango supports a role for rootstock as well as scion identity in leaf nutrient content (Sarkhosh et al., <xref ref-type="bibr" rid="B48">2021</xref>). This can be traced in some cases to soil resource acquisition. In grapevine, for example, contrasting effects on scion vigor were traced to differences in responsiveness to heterogeneous soil nitrogen and in transcriptomic reprogramming of nitrogen-related genes (Cochetel et al., <xref ref-type="bibr" rid="B15">2017</xref>). Further exploration of homologous root and rhizosphere traits in cacao is warranted.</p>
<p>Root associations with the rhizosphere microbiome, or rhizobiome, may also play an important role in soil resource acquisition and nutrient use efficiency in perennial systems but have not been fully characterized in cacao. Grafting is known to affect rhizobiome composition: Small but significant rootstock effects on rhizosphere bacterial communities have been identified in grafted tomato (Poudel et al., <xref ref-type="bibr" rid="B41">2019</xref>) and grapevine (D&#x00027;Amico et al., <xref ref-type="bibr" rid="B18">2018</xref>; Marasco et al., <xref ref-type="bibr" rid="B33">2018</xref>). These effects may be small [around 3% of variation in tomato rhizosphere bacterial communities, far lower than soil compartment and site effects (Poudel et al., <xref ref-type="bibr" rid="B41">2019</xref>)] but stronger than scion effects (Liu et al., <xref ref-type="bibr" rid="B31">2018</xref>). Furthermore, these root-rhizobiome interactions can affect nutrient uptake (D&#x00027;Amico et al., <xref ref-type="bibr" rid="B18">2018</xref>). Research in cacao has shown that management practices impact fungal communities (Ar&#x000E9;valo-Gardini et al., <xref ref-type="bibr" rid="B4">2020</xref>) and total microbial communities (Buyer et al., <xref ref-type="bibr" rid="B7">2017</xref>), but the effect of rootstock identity on root-soil-microbe interactions affecting plant nutrition has not been examined.</p>
<p>We grafted seedlings of eight rootstock populations with three scion genotypes under greenhouse conditions to test the individual and interactive effects of rootstock and scion identity on cacao seedling nutrition and rhizobiome composition. Furthermore, we investigated the nature of plant-soil-microbe interactions affecting nutrient uptake using a variety of statistical techniques. We hypothesized that (a) rootstock and scion identity would both affect leaf nutrient content; (b) rootstock would have a stronger influence than scion on rhizobiome composition; and (c) soil physicochemical parameters and microbial community composition would interactively affect plant nutrition. Because nutritionally relevant plant-microbe interactions could occur at many scales, we explored potential impacts of individual taxa, the entire rhizobiome, and subsets of bacterial and fungal communities.</p>
</sec>
<sec sec-type="methods" id="s2">
<title>Methods</title>
<sec>
<title>Grafting and Seedling Growth</title>
<p>Eighteen-month-old <italic>T. cacao</italic> seedlings representing eight different rootstock populations (BYNC, EQX3348, IMC14, PA107, SCA6, T294, T384, T484) were cleft grafted with budwood from three scions (Criollo 22, Matina 1-6, Pound 7) (<xref ref-type="table" rid="T1">Table 1</xref>). Seedlings were grafted at this developmental stage because graft take is highest when rootstock diameter is close to the diameter of the scion budwood and the seedling has matured enough to flush. These diverse rootstocks were chosen to represent a range of genetic backgrounds, but detailed information about growth and developmental characteristics was not available due to limited research. Similarly, scions were selected from widely planted elite commercial cultivars representing different cacao genetic groups and origins (<xref ref-type="table" rid="T1">Table 1</xref>) (Motamayor et al., <xref ref-type="bibr" rid="B35">2008</xref>; Cornejo et al., <xref ref-type="bibr" rid="B16">2018</xref>). The decision to use a diverse range of rootstocks and these specific scions was made to reflect commercial production, where a wide variety of rootstocks (often locally adapted seedlings) are commercially planted and grafted with a smaller pool of elite scions, which are chosen primarily for bean quality and yield. In addition, the larger number of rootstocks was chosen because we hypothesized that rootstock effects would be stronger than scion effects.</p>
<table-wrap position="float" id="T1">
<label>Table 1</label>
<caption><p>Cacao germplasm used in the study.</p></caption>
<table frame="hsides" rules="groups">
<thead><tr>
<th valign="top" align="left"><bold>Accession</bold></th>
<th valign="top" align="left"><bold>Graft material</bold></th>
<th valign="top" align="left"><bold>Provenance</bold></th>
<th valign="top" align="left"><bold>Pedigree or structure group</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Criollo 22</td>
<td valign="top" align="left">Scion</td>
<td valign="top" align="left">Costa Rica</td>
<td valign="top" align="left">Criollo</td>
</tr>
<tr>
<td valign="top" align="left">Matina 1-6</td>
<td valign="top" align="left">Scion</td>
<td valign="top" align="left">Costa Rica</td>
<td valign="top" align="left">Amelonado</td>
</tr>
<tr>
<td valign="top" align="left">Pound 7</td>
<td valign="top" align="left">Scion</td>
<td valign="top" align="left">R. Nanay, Loreto, Peru</td>
<td valign="top" align="left">Nanay</td>
</tr>
<tr>
<td valign="top" align="left">BYNC</td>
<td valign="top" align="left">Rootstock</td>
<td valign="top" align="left">Reassigned offtype, SHRS</td>
<td valign="top" align="left">Iquitos admixed</td>
</tr>
<tr>
<td valign="top" align="left">EQX 3348</td>
<td valign="top" align="left">Rootstock</td>
<td valign="top" align="left">R. Coca/R. Napo, Oriente, Ecuador</td>
<td valign="top" align="left">Iquitos</td>
</tr>
<tr>
<td valign="top" align="left">IMC 14</td>
<td valign="top" align="left">Rootstock</td>
<td valign="top" align="left">Marper farm, Trinidad</td>
<td valign="top" align="left">P 18 x P 21, Iquitos</td>
</tr>
<tr>
<td valign="top" align="left">PA 107</td>
<td valign="top" align="left">Rootstock</td>
<td valign="top" align="left">R. Maranon, Parinari, Peru</td>
<td valign="top" align="left">Iquitos/Maranon</td>
</tr>
<tr>
<td valign="top" align="left">SCA 6</td>
<td valign="top" align="left">Rootstock</td>
<td valign="top" align="left">R. Nanay, Iquitos, Peru</td>
<td valign="top" align="left">Contamana</td>
</tr>
<tr>
<td valign="top" align="left">T-294</td>
<td valign="top" align="left">Rootstock</td>
<td valign="top" align="left">Tecas Farm, INIAP, Ecuador</td>
<td valign="top" align="left">CCN-51 x B-60</td>
</tr>
<tr>
<td valign="top" align="left">T-384</td>
<td valign="top" align="left">Rootstock</td>
<td valign="top" align="left">Tecas Farm, INIAP, Ecuador</td>
<td valign="top" align="left">CCN-51 x TAP-3</td>
</tr>
<tr>
<td valign="top" align="left">T-484</td>
<td valign="top" align="left">Rootstock</td>
<td valign="top" align="left">Tecas Farm, INIAP, Ecuador</td>
<td valign="top" align="left">AMAZ-11 x CCN-51</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p><italic>SHRS, Subtropical Horticultural Research Station, U.S. Department of Agriculture; INIAP, Instituto Nacional de Investigaciones Agropecuarias</italic>.</p>
</table-wrap-foot>
</table-wrap>
<p>Grafted seedlings were maintained under controlled greenhouse conditions in the Mars Breeding Greenhouse at the USDA SHRS in Miami, Florida, USA. Seedlings were maintained in 2 gallon polybags in a substrate of 37% pine bark, 25% perlite, 8% sand, and 30% coir with 6 lbs/yard dolomite and 1.66 lbs/yard micro max. This substrate mixture was autoclaved when first mixed to ensure homogeneous initial conditions, then transferred into polybags. Physical, chemical, and microbiological characteristics were not measured in each polybag at this time because the substrate had been thoroughly mixed and was assumed to be the same in all polybags. At 18 months, the plants were transferred to 7 gallon pots. The seedlings received an application of 15 ml of NUTRICOTE 18-6-18 fertilizer (Arysta LifeScience America Inc., New York, NY, USA) every 6 months and applications of 5 ml of K-Mag 0-0-22 (Mosaic Crop Nutrition, Tampa, FL, USA), which contains 10.8% magnesium and 22% sulfur, at planting and transplanting. Greenhouse temperatures were maintained between 68 and 85&#x000B0;F and humidity was between 30 and 95%. Photosynthetically active radiation ranged from 161 to 485 umol/s/m<sup>2</sup> and there was no supplemental lighting.</p>
</sec>
<sec>
<title>Leaf and Soil Nutrient Analysis</title>
<p>At 4 and 16 months after grafting, leaf tissue was sampled for analysis to assess how grafting affected scion nutrition soon after graft take and a year after establishment. Five leaves were sampled from recently hardened flushes of each plant. Leaves were dried, ground, and analyzed for nutrient concentrations by Waypoint Analytical (Richmond, VA, USA).</p>
<p>To allow time for rootstock effects on soil and microbial communities to develop, soil was collected at 20 months after grafting for analysis of macro- and micronutrients as well as rhizosphere microbial communities. Soil for nutrient analysis was sieved to pass a 2 mm mesh. pH was measured in a 1:1 soil:water slurry, soil organic matter (SOM) was measured via loss-on-ignition, and plant-available nutrient concentrations were measured with the Mehlich 3 method (Mehlich, <xref ref-type="bibr" rid="B34">1984</xref>). Cation exchange capacity (CEC) was calculated from measured cation concentrations and soil pH.</p>
</sec>
<sec>
<title>Bioinformatics</title>
<p>Genomic DNA was extracted from 0.25 g of rhizosphere soil using DNeasy PowerSoil Pro kits (QIAGEN, Germantown, MD, USA). Bacterial and fungal rhizosphere communities were assessed by sequencing of the V4 and ITS1 regions, respectively, on an Illumina Novaseq platform. The original 515F/806R primer pair was used to amplify the V4 region (Caporaso et al., <xref ref-type="bibr" rid="B10">2011</xref>) and the ITS1F/ITS2 primer pair was used for the ITS1 region (White et al., <xref ref-type="bibr" rid="B53">1990</xref>; Gardes and Bruns, <xref ref-type="bibr" rid="B22">1993</xref>).</p>
<p>Primers and adapters were removed and low-quality reads (length &#x0003C; 60 bp, Q score &#x0003C; 5, or unidentified bases &#x0003E; 10%) were filtered out using fastp software (Chen et al., <xref ref-type="bibr" rid="B12">2018</xref>). Subsequent analyses were done in R software v4.0.3 (R Core Team, <xref ref-type="bibr" rid="B44">2020</xref>). Reads were filtered using the dada2 package v1.18 (Callahan et al., <xref ref-type="bibr" rid="B8">2016</xref>) with a maximum of two expected errors. Paired-end reads were denoised, merged, and used to generate amplicon sequence variant (ASV) tables separately for 16S and ITS sequences, containing 44,150 unique prokaryotic ASVs and 5,153 fungal ASVs. Chimeras were removed for both 16S and ITS sequences. For 16S sequences, products of non-specific priming that fell outside the anticipated 250-256 bp length for the V4 region, chloroplast sequences, and mitochondrial sequences were additionally removed. Taxonomy was assigned using the SILVA reference database v138.1 for prokaryotes (Quast et al., <xref ref-type="bibr" rid="B43">2013</xref>; Yilmaz et al., <xref ref-type="bibr" rid="B54">2014</xref>) and the UNITE database v8.2 for fungi (Abarenkov et al., <xref ref-type="bibr" rid="B1">2020</xref>), and species-level identification was done only in the case of 100% sequence identity. Chloroplasts and mitochondria were removed and subsequent analyses of 16S sequences were restricted to bacteria only.</p>
</sec>
<sec>
<title>Statistical Analysis</title>
<p>All analyses were done in R software v4.0.3 (R Core Team, <xref ref-type="bibr" rid="B44">2020</xref>). Leaf nutrient data were scaled to have unit variance and principal components analysis (PCA) was used to ordinate data from each sampling date separately [prcomp function of stats package (R Core Team, <xref ref-type="bibr" rid="B44">2020</xref>)].</p>
<p>Microbial diversity and community composition were investigated separately for bacterial and fungal communities using the vegan package v.2.5-7 (Oksanen et al., <xref ref-type="bibr" rid="B38">2020</xref>). Alpha diversity metrics included richness, Shannon index, and Pielou index (taxonomy-based methods) and Faith&#x00027;s PD (a phylogeny-based method). The effects of rootstock, scion, and the rootstock:scion interaction on microbial alpha diversity metrics were tested with ANOVA followed by <italic>post-hoc</italic> Tukey tests. Microbial communities were ordinated with principal coordinates analysis (PCoA) based on Bray-Curtis distance matrices of log-transformed count data. Effects of rootstock, scion, and the rootstock:scion interaction on community composition were tested with permutational multivariate analysis of variance (PERMANOVA) with 1,000 permutations (Adonis function of vegan package; Oksanen et al., <xref ref-type="bibr" rid="B38">2020</xref>).</p>
<p>Because the functional relevance of the microbiome for plant nutrition was of interest, we sought to identify individual taxa or groups of taxa associated with variation in leaf nutrients. Kendall correlations between leaf nutrients and microbial ASVs were calculated to test potential direct effects of individual bacterial and fungal taxa on plant nutrition (cor function of stats package). Full Mantel tests were carried out on Bray-Curtis dissimilarity matrices of leaf nutrients, bacterial communities, and fungal communities to assess direct relationships among these datasets (mantel function of vegan package). Partial Mantel tests were used to test correlations between plant and microbial communities when soil nutrients were held constant. Relationships among soil nutrients, microbial communities, and leaf nutrients were then investigated with stepwise regression. Individual models were established for each leaf nutrient from the 16-month sampling point. Predictor variables included all soil nutrients, all measures of microbial diversity discussed previously, and the first two eigenvectors of bacterial and fungal communities. Stepwise regression was implemented and variable importance was calculated for each predictor variable using the leapSeq method in the train function of the caret package v6.0.86 with 10-fold cross-validation and all possible variables (Kuhn, <xref ref-type="bibr" rid="B29">2020</xref>).</p>
<p>Linear discriminant analysis effect size (LEfSe), an ensemble learning approach that identifies features responsible for differentiation among groups, was used to identify rootstock-responsive bacterial and fungal ASVs with potential impacts on plant nutrition. Rootstock-responsive taxa were defined as those that differ significantly between rootstock genotypes. LEfSE identifies microbiome features that may explain differences among groups, tests their significance with Kruskal-Wallis and Wilcoxon rank-sum tests, then uses redundancy analysis to estimate the effect size associated with each feature (Segata et al., <xref ref-type="bibr" rid="B49">2011</xref>). This method was implemented with the lefse function of the microbiomeMarker package v0.0.1.9000 (Cao, <xref ref-type="bibr" rid="B9">2021</xref>). Functionally relevant taxa were prioritized by analyzing the 50 bacterial and fungal ASVs with the highest loadings on the first two eigenvectors, which were shown to be important for plant nutrition in the stepwise regression analysis.</p>
</sec>
</sec>
<sec sec-type="results" id="s3">
<title>Results</title>
<sec>
<title>Relationships Among Leaf Nutrients</title>
<p>Ranges of leaf nutrient values for grafted seedlings 4 and 16 months after grafting and normal reference ranges are given in <xref ref-type="table" rid="T2">Table 2</xref>. Four months after grafting, PCA revealed that 46% of variation in leaf nutrients was accounted for by the first two eigenvectors (<xref ref-type="fig" rid="F1">Figure 1A</xref>). Ca, Mg, and Zn had the highest loadings on the first eigenvector, while N, Cu, and P had the highest loadings on the second eigenvector. Significant correlations were also observed between Ca &#x02013; Mg (<italic>p</italic> &#x0003C; 0.001), Ca &#x02013; Zn (<italic>p</italic> &#x0003C; 0.05), and P &#x02013; Cu (<italic>p</italic> &#x0003C; 0.05; <xref ref-type="supplementary-material" rid="SM1">Supplementary Figure 1</xref>). At 16 months, variation among samples was driven by K, P, Ca, and B in the first eigenvector and Fe, Cu, and Al in the second eigenvector (<xref ref-type="fig" rid="F1">Figure 1B</xref>). Correlations among nutrients were stronger at this point, with three main nutrient groupings that were positively correlated with one another and negatively correlated with other groupings: P-K, S-Ca-Mg-B-Zn-Mn, Fe-Al (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure 2</xref>). Some clustering by rootstock was apparent at 4 months, and two groups of rootstocks (BYNC, EQX3348, IMC14, PA107 vs. SCA6, T294, T384, T484) were visible at 16 months (<xref ref-type="fig" rid="F1">Figure 1</xref>).</p>
<table-wrap position="float" id="T2">
<label>Table 2</label>
<caption><p>Ranges of leaf nutrient values in grafted T. cacao seedlings.</p></caption>
<table frame="hsides" rules="groups">
<thead><tr>
<th valign="top" align="left"><bold>Nutrient</bold></th>
<th valign="top" align="center"><bold>4 months</bold></th>
<th valign="top" align="center"><bold>16 months</bold></th>
<th valign="top" align="center"><bold>Normal reference range</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Aluminum</td>
<td valign="top" align="center">18&#x02013;59 ppm</td>
<td valign="top" align="center">1&#x02013;80 ppm</td>
<td valign="top" align="center">0&#x02013;100 ppm</td>
</tr>
<tr>
<td valign="top" align="left">Boron</td>
<td valign="top" align="center">36&#x02013;104 ppm</td>
<td valign="top" align="center">19&#x02013;69 ppm</td>
<td valign="top" align="center">25&#x02013;70 ppm</td>
</tr>
<tr>
<td valign="top" align="left">Calcium</td>
<td valign="top" align="center">0.84&#x02013;2.25%</td>
<td valign="top" align="center">0.27&#x02013;2.29%</td>
<td valign="top" align="center">0.3&#x02013;0.59%</td>
</tr>
<tr>
<td valign="top" align="left">Copper</td>
<td valign="top" align="center">7&#x02013;15 ppm</td>
<td valign="top" align="center">6&#x02013;20 ppm</td>
<td valign="top" align="center">8&#x02013;12 ppm</td>
</tr>
<tr>
<td valign="top" align="left">Iron</td>
<td valign="top" align="center">39&#x02013;357 ppm</td>
<td valign="top" align="center">9&#x02013;400 ppm</td>
<td valign="top" align="center">60&#x02013;200 ppm</td>
</tr>
<tr>
<td valign="top" align="left">Magnesium</td>
<td valign="top" align="center">0.59&#x02013;1.48%</td>
<td valign="top" align="center">0.38&#x02013;1.12%</td>
<td valign="top" align="center">0.2&#x02013;0.49%</td>
</tr>
<tr>
<td valign="top" align="left">Manganese</td>
<td valign="top" align="center">150&#x02013;761 ppm</td>
<td valign="top" align="center">58&#x02013;841 ppm</td>
<td valign="top" align="center">50&#x02013;300 ppm</td>
</tr>
<tr>
<td valign="top" align="left">Nitrogen</td>
<td valign="top" align="center">1.14&#x02013;2.22%</td>
<td valign="top" align="center">1.58&#x02013;3.21%</td>
<td valign="top" align="center">2&#x02013;2.49%</td>
</tr>
<tr>
<td valign="top" align="left">Phosphorus</td>
<td valign="top" align="center">0.16&#x02013;0.33%</td>
<td valign="top" align="center">0.11&#x02013;0.35%</td>
<td valign="top" align="center">0.13&#x02013;0.24%</td>
</tr>
<tr>
<td valign="top" align="left">Potassium</td>
<td valign="top" align="center">1.03&#x02013;2.67%</td>
<td valign="top" align="center">1.31&#x02013;3.94%</td>
<td valign="top" align="center">1.3&#x02013;2.19%</td>
</tr>
<tr>
<td valign="top" align="left">Sodium</td>
<td valign="top" align="center">0.02&#x02013;0.04%</td>
<td valign="top" align="center">0.01&#x02013;0.19%</td>
<td valign="top" align="center">0&#x02013;0.19%</td>
</tr>
<tr>
<td valign="top" align="left">Sulfur</td>
<td valign="top" align="center">0.18&#x02013;0.48%</td>
<td valign="top" align="center">0.21&#x02013;0.7%</td>
<td valign="top" align="center">0.16&#x02013;0.26%</td>
</tr>
<tr>
<td valign="top" align="left">Zinc</td>
<td valign="top" align="center">70&#x02013;335 ppm</td>
<td valign="top" align="center">33&#x02013;168 ppm</td>
<td valign="top" align="center">20&#x02013;200 ppm</td>
</tr>
</tbody>
</table>
</table-wrap>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption><p>Principal components analysis (PCA) of cacao leaf nutrients at two sampling dates. Values were scaled to have unit variance prior to analysis. <bold>(A)</bold> At 4 months, 46% of variation in leaf nutrients was accounted for by the first two eigenvectors. Ca, Mg, and Zn had the highest loadings for the first eigenvector (Dim1 above) and N, Cu, and P had the highest loadings for the second eigenvector (Dim2). <bold>(B)</bold> At 16 months, 55% of variation was accounted for by the first two eigenvectors, which were driven by K, P, Ca, and B (Dim1) and Fe, Cu, and Al (Dim2).</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fagro-03-754646-g0001.tif"/>
</fig>
</sec>
<sec>
<title>Rootstock and Scion Effects on Leaf Nutrients</title>
<p>Rootstock, scion, and the rootstock:scion interaction significantly affected leaf nutrients at both sampling dates. At 4 months, the rootstock:scion interaction was significant for B (<italic>p</italic> &#x0003C; 0.05) and rootstock effects were significant for N (<italic>p</italic> &#x0003C; 0.001), P (<italic>p</italic> &#x0003C; 0.01), K (<italic>p</italic> &#x0003C; 0.05), Mg (<italic>p</italic> &#x0003C; 0.05), Ca (<italic>p</italic> &#x0003C; 0.001), Zn (<italic>p</italic> &#x0003C; 0.001), Mn (<italic>p</italic> &#x0003C; 0.01), Fe (<italic>p</italic> &#x0003C; 0.001), Cu (<italic>p</italic> &#x0003C; 0.05), and Al (<italic>p</italic> &#x0003C; 0.001). Only Na responded to scion genotype (<italic>p</italic> &#x0003C; 0.01). At 16 months, fewer nutrients were affected by rootstock or scion identity. B (<italic>p</italic> &#x0003C; 0.05), Fe (<italic>p</italic> &#x0003C; 0.001), and Al (<italic>p</italic> &#x0003C; 0.001) were affected by the rootstock:scion interaction, while only Na (<italic>p</italic> &#x0003C; 0.01) and Cu (<italic>p</italic> &#x0003C; 0.001) differed among rootstocks. No nutrients were affected by scion identity at 16 months.</p>
</sec>
<sec>
<title>Microbial Alpha Diversity</title>
<p>Rootstock identity (but not scion or the rootstock:scion interaction) significantly affected both phylogeny-based and taxonomy-based measures of bacterial alpha diversity (all <italic>p</italic> &#x0003C; 0.001; <xref ref-type="fig" rid="F2">Figure 2A</xref>). A trend was observed for all measures in which bacterial communities associated with SCA6, T294, T384, and T484 were less diverse than bacterial communities associated with BYNC, EQX3348, IMC14, and PA107, although this distinction was significant only for richness.</p>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption><p>Cacao rootstock effects on microbial alpha diversity. <bold>(A)</bold> In bacterial communities, Faith&#x00027;s PD (a diversity metric that accounts for phylogeny), the Pielou and Shannon indices (integrating richness and evenness), and taxonomic richness were significantly affected by rootstock (all <italic>p</italic> &#x0003C; 0.05). Lower diversity was observed in bacterial communities associated with SCA6, T294, T384, and T484. <bold>(B)</bold> In fungal communities, rootstock identity significantly affected all measures of alpha diversity (all <italic>p</italic> &#x0003C; 0.05). Fungi associated with T384 showed the lowest alpha diversity.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fagro-03-754646-g0002.tif"/>
</fig>
<p>Rootstock identity (but not scion or rootstock:scion interaction) similarly affected all measures of fungal alpha diversity (richness <italic>p</italic> &#x0003C; 0.05, Shannon index <italic>p</italic> &#x0003C; 0.01, Pielou index <italic>p</italic> &#x0003C; 0.05, Faith&#x00027;s PD <italic>p</italic> &#x0003C; 0.05; <xref ref-type="fig" rid="F2">Figure 2B</xref>). Unlike for bacteria, however, decreased diversity was observed only in rootstock T384.</p>
</sec>
<sec>
<title>Microbial Community Composition</title>
<p>PCoA showed clustering of bacterial and fungal communities by two groups of rootstocks (<xref ref-type="fig" rid="F3">Figure 3</xref>). Rhizosphere communities associated with BYNC, EQX3348, IMC14, and PA107 formed a distinct cluster separated from communities associated with SCA6, T294, T384, and T48. PERMANOVA revealed significant effects of rootstock (bacterial <italic>R</italic><sup>2</sup> = 0.43, fungal <italic>R</italic><sup>2</sup> = 0.19, both <italic>p</italic> &#x0003C; 0.001) but not scion or the rootstock:scion interaction on bacterial and fungal communities. The rootstock effect was also observed when the ordination was based on weighted UniFrac distance, which accounts for phylogeny, rather than Bray-Curtis distance (both PERMANOVA <italic>p</italic> &#x0003C; 0.001, data not shown).</p>
<fig id="F3" position="float">
<label>Figure 3</label>
<caption><p>Principal coordinates analysis (PCoA) ordination of rhizobiomes. Separation of <bold>(A)</bold> bacterial communities and <bold>(B)</bold> fungal communities by two groups of rootstocks occurred primarily along the first axis, which accounted for 26.9% of variation in bacterial communities and 13.7% of variation in fungal communities. Further clustering by rootstock within each group was evident along the second axis, which accounted for 5% of variation in both bacterial and fungal communities.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fagro-03-754646-g0003.tif"/>
</fig>
<p>High bacterial richness and evenness was observed at the order level, with 303 unique orders represented by at least one ASV (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure 3</xref>). The Rhizobiales were the most abundant order with 2,486 ASVs across all samples, followed by the Burkholderiales (1,124) and Chitinophagales (848). Fungal ASVs were distributed across 74 orders, but samples were dominated by fewer taxonomic groups (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure 4</xref>). Although over half of fungal sequences were not identifiable to the order level, the Agaricales and Hypocreales were the highest-abundance identifiable orders with 505 and 259 ASVs, respectively.</p>
</sec>
<sec>
<title>Plant-Soil-Microbe Interactions Affecting Leaf Nutrient Content</title>
<p>Correlation analysis did not identify individual bacterial or fungal ASVs whose abundance was strongly and significantly correlated with the concentration of any leaf nutrient (Kendall correlation &#x0003E; 0.75, <italic>p</italic> &#x0003C; 0.05). Full and partial Mantel tests based on Bray-Curtis dissimilarity matrices were used to test for significant relationships among leaf nutrients, soil nutrients, and bacterial and fungal communities. Leaf nutrients did not vary significantly with either bacterial or fungal communities (full Mantel test, both <italic>p</italic> &#x0003E; 0.05), even if soil nutrients were held constant (partial Mantel test, both <italic>p</italic> &#x0003E; 0.05). However, bacterial and fungal communities were correlated with one another (full Mantel test, <italic>p</italic> &#x0003C; 0.05).</p>
<p>Stepwise regression was implemented for each leaf nutrient individually to investigate the combinations of soil physicochemical properties and microbial parameters that affect foliar nutrition (<xref ref-type="supplementary-material" rid="SM1">Supplementary Tables 1&#x02013;5</xref>). Analysis of the relative variable importance in each model showed that measures of microbial diversity and community composition contributed substantially to variation in many leaf nutrients (<xref ref-type="fig" rid="F4">Figure 4</xref>). For example, richness and Faith&#x00027;s PD of bacterial communities were the most important predictor variables in the regression models defined for leaf N, Ca, Zn, and Al, and the second eigenvector of fungal communities (PC2F) was the most important predictor variable for leaf Na. Among soil physicochemical properties, only SOM and Zn had the highest variable importance in multiple models (SOM: leaf Mg, Mn; soil Zn: leaf P, K, and Cu).</p>
<fig id="F4" position="float">
<label>Figure 4</label>
<caption><p>Variable importance heatmap and dendrograms for leaf nutrient models. Stepwise regression analysis identified linear models with the combination of soil physicochemical properties and microbial diversity and composition metrics that best explained the concentrations of each leaf nutrient individually. The heatmap portrays the relative importance of each explanatory variable (columns) in each linear model (rows by response variable). Hierarchical clustering dendrograms along the &#x0201C;x&#x0201D; and &#x0201C;y&#x0201D; axes show similarity among variables and models, respectively. PC1B, PC2B: first and second principal components of bacterial community; PC1F, PC2F: first and second principal components of fungal community; CEC: cation exchange capacity; [B]: bacteria; [F]: fungi.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fagro-03-754646-g0004.tif"/>
</fig>
</sec>
<sec>
<title>Rootstock-Responsive Bacteria and Fungi</title>
<p>The Rhizobiales were the most-represented order among rootstock-responsive bacteria potentially related to plant nutrition, and these were enriched in BYNC, EQX3348, and T384 (<xref ref-type="fig" rid="F5">Figure 5</xref>). BYNC, EQX3348, T294, and T384 were all enriched in multiple bacterial orders. In contrast, IMC14 was enriched only in Pseudomonodales, SCA6 in the Bacillales, and PA107 and T484 in groups not identified at the order level.</p>
<fig id="F5" position="float">
<label>Figure 5</label>
<caption><p>Distribution of rootstock-responsive taxa with potential impacts on plant nutrition. The linear discriminant analysis effect size (LEfSe) method identified <bold>(A)</bold> bacteria and <bold>(B)</bold> fungi responsible for differentiation among rootstocks (LDA &#x0003E; 2, <italic>p</italic> &#x0003C; 0.05). To focus the analysis on taxa potentially relevant for plant nutrition, LEfSe was applied to bacterial and fungal amplicon sequence variants (ASVs) with the highest loadings on the first two eigenvectors used in leaf nutrient models (PC1B, PC2B, PC1F, and PC2F in <xref ref-type="fig" rid="F4">Figure 4</xref>). ASVs are identified to order and classified as &#x0201C;Unidentified&#x0201D; if only higher taxonomic levels were available. &#x0201C;Other&#x0201D; is used to represent orders with fewer than three members. Dots represent the number of ASVs from a given order whose relative abundance was enriched in one rootstock.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fagro-03-754646-g0005.tif"/>
</fig>
<p>BYNC and IMC14 were enriched in the Agaricales, the most-represented order among root-responsive fungi (<xref ref-type="fig" rid="F5">Figure 5</xref>). Other enriched orders were the Pezizales (in BYNC), Glomerellales (in IMC14), and Chaetothyriales (in T294). Enriched fungal groups in PA107 and T484 were not identified at the order level. EQX3348, SCA6, and T384 were not enriched in any rootstock-responsive taxa.</p>
<p>Rootstock-responsive bacteria and fungi were similar to the overall rhizobiome composition (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figures 3, 4</xref>). The most abundant bacterial orders were consistent between the two datasets, but the composition of fungal communities was more variable. In addition, a large number of fungal sequences could be identified only at the kingdom, phylum, or class levels, making higher-resolution analysis of community composition difficult.</p>
</sec>
</sec>
<sec sec-type="discussion" id="s4">
<title>Discussion</title>
<p>Our findings of dominant rootstock effects on plant nutrition and microbial interactions in cacao suggest that rootstock breeding could be employed to increase nutrient use efficiency. In addition, this strategy would not compromise cacao bean quality, which is under scion control (Prawoto, <xref ref-type="bibr" rid="B42">1990</xref>). Differences in plant nutrition among rootstocks and scions appear to be important primarily for the first few months after grafting. If differences had been observed among rootstocks in nutrient accumulation at both time points, it would have warranted investigation of the genetic basis of differences in nutrient uptake and partitioning as well as rootstock-specific microbial taxa that may have contributed. As the plant matures, plant physiological processes may demand stoichiometric nutrient ratios regardless of plant genotype, an explanation supported by the stronger correlations among nutrients at sixteen as compared to 4 months after grafting (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figures 1, 2</xref>). Nonetheless, juvenile vigor of rootstocks and young grafts is a relevant trait for field production and a better understanding of the genes responsible for differentiation among rootstocks and scions in this study could help inform cacao breeding programs. A useful next step would be to validate these results by planting the same genotypes in the field under fertile and nutrient-limited conditions.</p>
<p>In accordance with our second hypothesis, rootstock effects outweighed scion effects on root-associated microbes. Bacterial diversity and community composition of both bacteria and fungi followed a similar trend in which rootstocks clustered in two groups: BYNC/EQX3348/IMC14/PA107 and SCA6/T294/T384/T484. Identifying the mechanisms responsible for the higher bacterial diversity in the first group and differences in rhizobiome composition was beyond the scope of this study, but we suggest that genome sequencing and analysis of root exudates from all eight rootstock populations would be useful starting points. The functional implications of reduced bacterial and fungal diversity for nutrient uptake by rootstocks SCA6, T294, T384, and T484 also deserve further investigation. D&#x00027;Amico et al. (<xref ref-type="bibr" rid="B18">2018</xref>) found reduced rhizosphere bacterial diversity in one <italic>Vitis vinifera</italic> rootstock, including lower abundance of potassium-solubilizing taxa, that was hypothesized to be related to frequent potassium deficiencies in that rootstock.</p>
<p>We used multiple statistical methods to test our hypothesis that soil physicochemical parameters and microbial communities would impact plant nutrition, and to understand the nature of these relationships. The results of stepwise regression show that it may be important to develop new approaches to cacao nutrition that factor in soil biology as well as fertility. However, because this was a greenhouse study, limited variation in soil physicochemical parameters may have led to underestimation of the importance of these variables for plant nutrition. Applying the statistical methods used in this study to field data would be a valuable area for future research. Soil pH, for example, strongly affects nutrient availability and microbial communities in tropical agroecosystems (Juo and Franzluebbers, <xref ref-type="bibr" rid="B27">2003</xref>), but did not show up as particularly important in stepwise regression, likely because of low variability in this study (<xref ref-type="fig" rid="F4">Figure 4</xref>). Interestingly, despite limited variation, soil organic matter (SOM, also a key management concern in tropical soils and particularly in intensive cacao production systems) was the most important variable in multiple plant nutrient models (<xref ref-type="fig" rid="F4">Figure 4</xref>). Increasing SOM in soils under cacao should remain a high priority given its apparent role in plant nutrition as well as significance for microbially-mediated soil health and C sequestration.</p>
<p>We interpret the insignificant nutrient-microbe correlations and Mantel tests and the results of stepwise regression as a call for more complex understanding of plant-soil-microbe interactions at intermediate scales, e.g., between binary associations and entire communities. Identification of functional guilds, statistical methods to identify the most pertinent features of high-dimensionality datasets, and synthetic-community-based studies could prove useful here among other methods (recently reviewed by Gupta et al., <xref ref-type="bibr" rid="B25">2021</xref>). Plant growth promotion or enhanced plant nutrient content mediated by these intermediate-scale groups of microorganisms could result from direct mechanisms (e.g., direct nutrient transfer by symbionts) or indirect mechanisms (e.g., increased bioavailability of soil nutrients), but in either case would be an emergent property of the group rather than attributable to a single organism. Still, we gain insight into potentially beneficial members of the rootstock-responsive taxa identified in this study from single-strain or low-complexity inoculation studies common in the literature.</p>
<p>Many of the rootstock-responsive bacterial taxa identified belonged to orders associated with plant growth promotion, although incomplete taxonomic assignments limit our ability to confirm the presence of plant-growth-promoting rhizobacteria (PGPR) here. For example, the Rhizobiales, Frankiales, and Xanthomonodales, identified as enriched in BYNC, EQX3384, T294, and T384 (<xref ref-type="fig" rid="F5">Figure 5</xref>) include free-living and associative nitrogen-fixing bacteria. Unfortunately, these could be identified to the genus level only in some cases (<italic>Pedomicrobium, Devosia, Mycoplana, Hyphomicrobium, Bauldia, Rhodanobacter, Chujaibacter</italic>, and <italic>Jatrophihabitans</italic>), and N fixation is not a conserved trait for all of these genera (Rivas et al., <xref ref-type="bibr" rid="B46">2002</xref>; DeAngelis and Pold, <xref ref-type="bibr" rid="B19">2020</xref>). The Rhizobiales, which include both diazotrophs and pathogens (Carvalho et al., <xref ref-type="bibr" rid="B11">2010</xref>), were the most abundant order across all samples (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure 3</xref>). Future metagenomic studies of the cacao rhizobiome could screen for gene clusters involved in N fixation that might become part of biologically-based nutrient management strategies.</p>
<p>The Pseudomonodales and Bacillales, present at low relative abundance in all rhizobiome samples and identified in rootstock-responsive bacteria (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure 3</xref>, &#x0201C;Other&#x0201D;), contain taxa capable of increasing plant productivity in many plant hosts by both direct and indirect methods, e.g., producing growth-promoting auxins or suppressing pathogens (Kloepper et al., <xref ref-type="bibr" rid="B28">2004</xref>; Lugtenberg and Kamilova, <xref ref-type="bibr" rid="B32">2009</xref>). Furthermore, strains belonging to the Bacillales have shown promise in cacao. <italic>Bacillus subtilis</italic> inoculation of germinating seeds increased leaf area by 14%, shoot height by 8%, and graft success rate by 24% (Falc&#x000E4;o et al., <xref ref-type="bibr" rid="B21">2014</xref>), and repeated inoculations of seeds and seedlings increased shoot height by 22% (Leite et al., <xref ref-type="bibr" rid="B30">2013</xref>). Inoculation of 1-month-old cacao seedlings with other <italic>Bacillus</italic> spp. increased height and biomass by up to 37 and 73%, respectively (Gupta et al., <xref ref-type="bibr" rid="B24">2011</xref>). However, inoculation with <italic>Pseudomonas</italic> sp. did not promote cacao growth in another study (Guillermo Ram&#x000ED;rez et al., <xref ref-type="bibr" rid="B23">2016</xref>).</p>
<p>As with bacteria, incomplete taxonomic assignments limit our ability to clarify the ecological roles of many of the rootstock-responsive fungi. Unfortunately, many of the rootstock-responsive fungi could not be identified by genus or species, the taxonomic levels used for functional guild assignment in fungal trait databases such as Fun<sup>Fun</sup> (Zanne et al., <xref ref-type="bibr" rid="B56">2020</xref>), FUNGuild (Nguyen et al., <xref ref-type="bibr" rid="B37">2016</xref>), and FungalTraits (P&#x000F5;lme et al., <xref ref-type="bibr" rid="B40">2020</xref>). Two of the four fungal orders detected, the Glomerellales and the Chaetothyriales, have been found at low abundance (2.0&#x02013;2.3%) in the cacao phyllosphere, but their roles there were not defined (Wemheuer et al., <xref ref-type="bibr" rid="B52">2020</xref>). The majority of studies investigating fungi as potential plant-growth-promoting agents in cacao have focused on mycorrhizae (Chulan and Ragu, <xref ref-type="bibr" rid="B13">1986</xref>; Cuenca et al., <xref ref-type="bibr" rid="B17">1990</xref>; Chulan and Martin, <xref ref-type="bibr" rid="B14">1992</xref>; Tchameni et al., <xref ref-type="bibr" rid="B51">2012</xref>; Guillermo Ram&#x000ED;rez et al., <xref ref-type="bibr" rid="B23">2016</xref>) or <italic>Trichoderma</italic> sp. (Bae et al., <xref ref-type="bibr" rid="B5">2009</xref>; Tchameni et al., <xref ref-type="bibr" rid="B50">2011</xref>). While 112 ASVs belonging to the genus <italic>Trichoderma</italic> and 8 belonging to the mycorrhizal order Glomerales were found in fungal rhizobiomes in this study, none were among the rootstock-responsive taxa (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure 4</xref>).</p>
<p>The possible presence of pathogenic taxa among rootstock-responsive fungi could have had detrimental impacts on leaf nutrition in diseased plants. The Glomerellales (not to be confused with the Glomerales) contain <italic>Colletotrichum</italic> sp., cacao pathogens responsible for anthracnose (Nascimento et al., <xref ref-type="bibr" rid="B36">2019</xref>). While two of the rootstock-responsive Glomerellales strains were identified as <italic>Musicillium</italic> sp., the other two were not identified to the genus level and the possible presence of <italic>Colletotrichum</italic> sp. could not be ruled out. A member of the Hypocreales enriched in T484 (<xref ref-type="fig" rid="F5">Figure 5</xref>, &#x0201C;Other&#x0201D;) was identified as <italic>Fusarium</italic> sp. This genus is present at low relative abundance in cacao agroforestry systems (e.g., 1.3&#x02013;6.0%, Ar&#x000E9;valo-Gardini et al., <xref ref-type="bibr" rid="B4">2020</xref>). While it contains pathogens that cause cushion gall, dieback, and other diseases in cacao (Ploetz, <xref ref-type="bibr" rid="B39">2006</xref>; Adu-Acheampong et al., <xref ref-type="bibr" rid="B2">2012</xref>; Rosmana et al., <xref ref-type="bibr" rid="B47">2013</xref>), it also contains strains that may provide biocontrol against <italic>Phytophthora palmivora</italic>, the causal agent of black pod rot (Hanada et al., <xref ref-type="bibr" rid="B26">2010</xref>).</p>
<p>Despite limited attention to cacao rootstocks and rhizosphere microbial associations to date, our results highlight the potential of coordinated breeding efforts in these areas, combined with informed soil nutrient management, to enhance plant nutrition in cacao. This will be a challenging endeavor for numerous reasons: time and cost required to generate new perennial germplasm, dispersed smallholder production adding logistical difficulties to participatory cacao research and breeding, and remaining knowledge gaps in the genetic basis of cacao-rhizobiome interactions. Nonetheless, significant interest by consumers and industry in sustainable cacao production will help advance the field and translate research into practice.</p>
</sec>
<sec sec-type="data-availability" id="s5">
<title>Data Availability Statement</title>
<p>The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: NCBI (accession: <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="PRJNA755178">PRJNA755178</ext-link>).</p>
</sec>
<sec id="s6">
<title>Author Contributions</title>
<p>AD, AP, AT, and JS contributed to the study design and data collection. JS, AD, and TC analyzed the data. JS wrote the manuscript with input from all authors. All authors approved the final version for publication.</p>
</sec>
<sec sec-type="COI-statement" id="conf1">
<title>Conflict of Interest</title>
<p>JS, AD, AT, and TC are employed by Mars Wrigley. The remaining author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s7">
<title>Publisher&#x00027;s Note</title>
<p>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.</p>
</sec>
</body>
<back>
<ack><p>The authors would like to acknowledge Alana Firl and David Mackill for bioinformatics and project support.</p>
</ack>
<sec sec-type="supplementary-material" id="s8">
<title>Supplementary Material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fagro.2021.754646/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fagro.2021.754646/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Data_Sheet_1.docx" id="SM1" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="web"><person-group person-group-type="author"><name><surname>Abarenkov</surname> <given-names>K.</given-names></name> <name><surname>Zirk</surname> <given-names>A.</given-names></name> <name><surname>Piirmann</surname> <given-names>T.</given-names></name> <name><surname>P&#x000F6;h&#x000F6;nen</surname> <given-names>R.</given-names></name> <name><surname>Ivanov</surname> <given-names>F.</given-names></name> <name><surname>Nilsson</surname> <given-names>R. H.</given-names></name> <etal/></person-group> (<year>2020</year>). <source>UNITE general FASTA Release for Fungi. Version 04.02.2020. UNITE Community</source>. Available online at: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.15156/BIO/786368">https://doi.org/10.15156/BIO/786368</ext-link> (accessed April 8, 2021).</citation>
</ref>
<ref id="B2">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Adu-Acheampong</surname> <given-names>R.</given-names></name> <name><surname>Archer</surname> <given-names>S.</given-names></name> <name><surname>Leather</surname> <given-names>S.</given-names></name></person-group> (<year>2012</year>). <article-title>Resistance to dieback disease caused by Fusarium and Lasiodiplodia species in cacao (Theobroma cacao L.) genotypes</article-title>. <source>Exp. Agric.</source> <volume>48</volume>, <fpage>85</fpage>&#x02013;<lpage>98</lpage>. <pub-id pub-id-type="doi">10.1017/S0014479711000883</pub-id><pub-id pub-id-type="pmid">30886898</pub-id></citation></ref>
<ref id="B3">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Albacete</surname> <given-names>A.</given-names></name> <name><surname>Mart&#x000ED;nez-And&#x000FA;jar</surname> <given-names>C.</given-names></name> <name><surname>Mart&#x000ED;nez-P&#x000E9;rez</surname> <given-names>A.</given-names></name> <name><surname>Thompson</surname> <given-names>A. J.</given-names></name> <name><surname>Dodd</surname> <given-names>I. C.</given-names></name> <name><surname>P&#x000E9;rez-Alfocea</surname> <given-names>F</given-names></name></person-group>. (<year>2015</year>). <article-title>Unravelling rootstock&#x000D7;scion interactions to improve food security</article-title>. <source>J. Exp. Bot.</source> <volume>66</volume>, <fpage>2211</fpage>&#x02013;<lpage>2226</lpage>. <pub-id pub-id-type="doi">10.1093/jxb/erv027</pub-id><pub-id pub-id-type="pmid">25754404</pub-id></citation></ref>
<ref id="B4">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ar&#x000E9;valo-Gardini</surname> <given-names>E.</given-names></name> <name><surname>Canto</surname> <given-names>M.</given-names></name> <name><surname>Alegre</surname> <given-names>J.</given-names></name> <name><surname>Ar&#x000E9;valo-Hern&#x000E1;ndez</surname> <given-names>C. O.</given-names></name> <name><surname>Loli</surname> <given-names>O.</given-names></name> <name><surname>Julca</surname> <given-names>A.</given-names></name> <etal/></person-group>. (<year>2020</year>). <article-title>Cacao agroforestry management systems effects on soil fungi diversity in the Peruvian Amazon</article-title>. <source>Ecol. Indic.</source> <volume>115</volume>:<fpage>106404</fpage>. <pub-id pub-id-type="doi">10.1016/j.ecolind.2020.106404</pub-id></citation>
</ref>
<ref id="B5">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bae</surname> <given-names>H.</given-names></name> <name><surname>Sicher</surname> <given-names>R. C.</given-names></name> <name><surname>Kim</surname> <given-names>M. S.</given-names></name> <name><surname>Kim</surname> <given-names>S.-H.</given-names></name> <name><surname>Strem</surname> <given-names>M. D.</given-names></name> <name><surname>Melnick</surname> <given-names>R. L.</given-names></name> <etal/></person-group>. (<year>2009</year>). <article-title>The beneficial endophyte Trichoderma hamatum isolate DIS 219b promotes growth and delays the onset of the drought response in Theobroma cacao</article-title>. <source>J. Exp. Bot.</source> <volume>60</volume>, <fpage>3279</fpage>&#x02013;<lpage>3295</lpage>. <pub-id pub-id-type="doi">10.1093/jxb/erp165</pub-id><pub-id pub-id-type="pmid">19564160</pub-id></citation></ref>
<ref id="B6">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Bailey</surname> <given-names>B. A.</given-names></name> <name><surname>Meinhardt</surname> <given-names>L. W.</given-names></name></person-group> (<year>2016</year>). <source>Cacao Diseases: A History of Old Enemies and New Encounters</source>. <publisher-loc>Basel</publisher-loc>: <publisher-name>Springer International Publishing</publisher-name>.</citation>
</ref>
<ref id="B7">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Buyer</surname> <given-names>J. S.</given-names></name> <name><surname>Baligar</surname> <given-names>V. C.</given-names></name> <name><surname>He</surname> <given-names>Z.</given-names></name> <name><surname>Ar&#x000E9;valo-Gardini</surname> <given-names>E.</given-names></name></person-group> (<year>2017</year>). <article-title>Soil microbial communities under cacao agroforestry and cover crop systems in Peru</article-title>. <source>Appl. Soil Ecol.</source> <volume>120</volume>, <fpage>273</fpage>&#x02013;<lpage>280</lpage>. <pub-id pub-id-type="doi">10.1016/j.apsoil.2017.09.009</pub-id></citation>
</ref>
<ref id="B8">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Callahan</surname> <given-names>B. J.</given-names></name> <name><surname>McMurdie</surname> <given-names>P. J.</given-names></name> <name><surname>Rosen</surname> <given-names>M. J.</given-names></name> <name><surname>Han</surname> <given-names>A. W.</given-names></name> <name><surname>Johnson</surname> <given-names>A. J. A.</given-names></name> <name><surname>Holmes</surname> <given-names>S. P.</given-names></name></person-group> (<year>2016</year>). <article-title>DADA2: High-resolution sample inference from Illumina amplicon data</article-title>. <source>Nat. Methods</source> <volume>13</volume>, <fpage>581</fpage>&#x02013;<lpage>583</lpage>. <pub-id pub-id-type="doi">10.1038/nmeth.3869</pub-id><pub-id pub-id-type="pmid">27214047</pub-id></citation></ref>
<ref id="B9">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Cao</surname> <given-names>Y.</given-names></name></person-group> (<year>2021</year>). <source>microbiomeMarker: Microbiome Biomarker Analysis.</source></citation>
</ref>
<ref id="B10">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Caporaso</surname> <given-names>J. G.</given-names></name> <name><surname>Lauber</surname> <given-names>C. L.</given-names></name> <name><surname>Walters</surname> <given-names>W. A.</given-names></name> <name><surname>Berg-Lyons</surname> <given-names>D.</given-names></name> <name><surname>Lozupone</surname> <given-names>C. A.</given-names></name> <name><surname>Turnbaugh</surname> <given-names>P. J.</given-names></name> <etal/></person-group>. (<year>2011</year>). <article-title>Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample</article-title>. <source>Proc. Natl. Acad. Sci. U. S. A.</source> <volume>108</volume>, <fpage>4516</fpage>&#x02013;<lpage>4522</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.1000080107</pub-id><pub-id pub-id-type="pmid">20534432</pub-id></citation></ref>
<ref id="B11">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Carvalho</surname> <given-names>F. M.</given-names></name> <name><surname>Souza</surname> <given-names>R. C.</given-names></name> <name><surname>Barcellos</surname> <given-names>F. G.</given-names></name> <name><surname>Hungria</surname> <given-names>M.</given-names></name> <name><surname>Vasconcelos</surname> <given-names>A. T. R.</given-names></name></person-group> (<year>2010</year>). <article-title>Genomic and evolutionary comparisons of diazotrophic and pathogenic bacteria of the order Rhizobiales</article-title>. <source>BMC Microbiol.</source> <volume>10</volume>, <fpage>1</fpage>&#x02013;<lpage>15</lpage>. <pub-id pub-id-type="doi">10.1186/1471-2180-10-37</pub-id><pub-id pub-id-type="pmid">20144182</pub-id></citation></ref>
<ref id="B12">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname> <given-names>S.</given-names></name> <name><surname>Zhou</surname> <given-names>Y.</given-names></name> <name><surname>Chen</surname> <given-names>Y.</given-names></name> <name><surname>Gu</surname> <given-names>J.</given-names></name></person-group> (<year>2018</year>). <article-title>fastp: an ultra-fast all-in-one FASTQ preprocessor</article-title>. <source>Bioinformatics</source> <volume>34</volume>, <fpage>i884</fpage>&#x02013;<lpage>i890</lpage>. <pub-id pub-id-type="doi">10.1093/bioinformatics/bty560</pub-id><pub-id pub-id-type="pmid">30423086</pub-id></citation></ref>
<ref id="B13">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chulan</surname> <given-names>A.</given-names></name> <name><surname>Ragu</surname> <given-names>P.</given-names></name></person-group> (<year>1986</year>). <article-title>Growth response of <italic>Theobroma cacao</italic> L. seedlings to inoculation with vesicular-arbuscular mycorrhizal fungi</article-title>. <source>Plant Soil</source> <volume>285</volume>, <fpage>279</fpage>&#x02013;<lpage>285</lpage>. <pub-id pub-id-type="doi">10.1007/BF02374771</pub-id></citation>
</ref>
<ref id="B14">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chulan</surname> <given-names>H. A.</given-names></name> <name><surname>Martin</surname> <given-names>K.</given-names></name></person-group> (<year>1992</year>). <article-title>The vesicular-arbuscular (VA) mycorrhiza and its effects on growth of vegetatively propagated <italic>Theobroma cacao</italic> L</article-title>. <source>Plant Soil</source> <volume>144</volume>, <fpage>227</fpage>&#x02013;<lpage>233</lpage>. <pub-id pub-id-type="doi">10.1007/BF00012879</pub-id></citation>
</ref>
<ref id="B15">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cochetel</surname> <given-names>N.</given-names></name> <name><surname>Escudi&#x000E9;</surname> <given-names>F.</given-names></name> <name><surname>Cookson</surname> <given-names>S. J.</given-names></name> <name><surname>Dai</surname> <given-names>Z.</given-names></name> <name><surname>Vivin</surname> <given-names>P.</given-names></name> <name><surname>Bert</surname> <given-names>P. F.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title>Root transcriptomic responses of grafted grapevines to heterogeneous nitrogen availability depend on rootstock genotype</article-title>. <source>J. Exp. Bot.</source> <volume>68</volume>, <fpage>4339</fpage>&#x02013;<lpage>4355</lpage>. <pub-id pub-id-type="doi">10.1093/jxb/erx224</pub-id><pub-id pub-id-type="pmid">28922755</pub-id></citation></ref>
<ref id="B16">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cornejo</surname> <given-names>O. E.</given-names></name> <name><surname>Yee</surname> <given-names>M.-C.</given-names></name> <name><surname>Dominguez</surname> <given-names>V.</given-names></name> <name><surname>Andrews</surname> <given-names>M.</given-names></name> <name><surname>Sockell</surname> <given-names>A.</given-names></name> <name><surname>Strandberg</surname> <given-names>E.</given-names></name> <etal/></person-group>. (<year>2018</year>). <article-title>Population genomic analyses of the chocolate tree, <italic>Theobroma cacao</italic> L., provide insights into its domestication process</article-title>. <source>Commun. Biol.</source> <volume>1</volume>:<fpage>167</fpage>. <pub-id pub-id-type="doi">10.1038/s42003-018-0168-6</pub-id><pub-id pub-id-type="pmid">30345393</pub-id></citation></ref>
<ref id="B17">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cuenca</surname> <given-names>G.</given-names></name> <name><surname>Herrera</surname> <given-names>R.</given-names></name> <name><surname>Meneses</surname> <given-names>E.</given-names></name></person-group> (<year>1990</year>). <article-title>Effects of VA mycorrhiza on the growth of cacao seedlings under nursery conditions in Venezuela</article-title>. <source>Plant Soil</source> <volume>126</volume>, <fpage>71</fpage>&#x02013;<lpage>78</lpage>. <pub-id pub-id-type="doi">10.1007/BF00041370</pub-id></citation>
</ref>
<ref id="B18">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>D&#x00027;Amico</surname> <given-names>F.</given-names></name> <name><surname>Candela</surname> <given-names>M.</given-names></name> <name><surname>Turroni</surname> <given-names>S.</given-names></name> <name><surname>Biagi</surname> <given-names>E.</given-names></name> <name><surname>Brigidi</surname> <given-names>P.</given-names></name> <name><surname>Bega</surname> <given-names>A.</given-names></name> <etal/></person-group>. (<year>2018</year>). <article-title>The rootstock regulates microbiome diversity in root and rhizosphere compartments of vitis vinifera cultivar lambrusco</article-title>. <source>Front. Microbiol.</source> <volume>9</volume>:<fpage>2240</fpage>. <pub-id pub-id-type="doi">10.3389/fmicb.2018.02240</pub-id><pub-id pub-id-type="pmid">30319569</pub-id></citation></ref>
<ref id="B19">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>DeAngelis</surname> <given-names>K. M.</given-names></name> <name><surname>Pold</surname> <given-names>G.</given-names></name></person-group> (<year>2020</year>). <article-title>Genome sequences of Frankineae sp. strain MT45 and Jatrophihabitans sp. strain GAS493, two actinobacteria isolated from forest soil</article-title>. <source>Microbiol. Resour. Announc.</source> <volume>9</volume>:<fpage>e00614</fpage>-<lpage>20</lpage>. <pub-id pub-id-type="doi">10.1128/MRA.00614-20</pub-id><pub-id pub-id-type="pmid">32943557</pub-id></citation></ref>
<ref id="B20">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>DuVal</surname> <given-names>A.</given-names></name> <name><surname>Gezan</surname> <given-names>S. A.</given-names></name> <name><surname>Mustiga</surname> <given-names>G.</given-names></name> <name><surname>Stack</surname> <given-names>C.</given-names></name> <name><surname>Marelli</surname> <given-names>J. P.</given-names></name> <name><surname>Chaparro</surname> <given-names>J.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title>Genetic parameters and the impact of off-types for <italic>Theobroma cacao</italic> L. in a breeding program in Brazil</article-title>. <source>Front. Plant Sci.</source> <volume>8</volume>:<fpage>2059</fpage>. <pub-id pub-id-type="doi">10.3389/fpls.2017.02059</pub-id><pub-id pub-id-type="pmid">29250097</pub-id></citation></ref>
<ref id="B21">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Falc&#x000E4;o</surname> <given-names>L. L.</given-names></name> <name><surname>Silva-Werneck</surname> <given-names>J. O.</given-names></name> <name><surname>Vilarinho</surname> <given-names>B. R.</given-names></name> <name><surname>da Silva</surname> <given-names>J. P.</given-names></name> <name><surname>Pomella</surname> <given-names>A. W. V.</given-names></name> <name><surname>Marcellino</surname> <given-names>L. H.</given-names></name></person-group> (<year>2014</year>). <article-title>Antimicrobial and plant growth-promoting properties of the cacao endophyte <italic>Bacillus subtilis</italic> ALB629</article-title>. <source>J. Appl. Microbiol.</source> <volume>116</volume>, <fpage>1584</fpage>&#x02013;<lpage>1592</lpage>. <pub-id pub-id-type="doi">10.1111/jam.12485</pub-id><pub-id pub-id-type="pmid">24905220</pub-id></citation></ref>
<ref id="B22">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gardes</surname> <given-names>M.</given-names></name> <name><surname>Bruns</surname> <given-names>T. D.</given-names></name></person-group> (<year>1993</year>). <article-title>ITS primers with enhanced specificity for basidiomycetes - application to the identification of mycorrhizae and rusts</article-title>. <source>Mol. Ecol.</source> <volume>2</volume>, <fpage>113</fpage>&#x02013;<lpage>118</lpage>. <pub-id pub-id-type="doi">10.1111/j.1365-294X.1993.tb00005.x</pub-id><pub-id pub-id-type="pmid">8180733</pub-id></citation></ref>
<ref id="B23">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Guillermo Ram&#x000ED;rez</surname> <given-names>J.</given-names></name> <name><surname>Osorno</surname> <given-names>L.</given-names></name> <name><surname>Walter Osorio</surname> <given-names>N.</given-names></name></person-group> (<year>2016</year>). <article-title>Presence of mycorrhizal fungi and a fluorescent Pseudomonas sp. in the rhizosphere of cacao in two agroecosystems and their effects on cacao seedling growth</article-title>. <source>Agron. Colomb.</source> <volume>34</volume>, <fpage>385</fpage>&#x02013;<lpage>392</lpage>. <pub-id pub-id-type="doi">10.15446/agron.colomb.v34n3.57950</pub-id></citation>
</ref>
<ref id="B24">
<citation citation-type="web"><person-group person-group-type="author"><name><surname>Gupta</surname> <given-names>A.</given-names></name> <name><surname>Gopal</surname> <given-names>M.</given-names></name> <name><surname>Thomas</surname> <given-names>G. V.</given-names></name></person-group> (<year>2011</year>). <source>Efficacy of Rhizospheric Bacillus spp. for Growth Promotion in Theobroma cacao L. seedlings</source>. Available online at: <ext-link ext-link-type="uri" xlink:href="https://www.researchgate.net/publication/237862366">https://www.researchgate.net/publication/237862366</ext-link> (accessed July 29, 2020).</citation>
</ref>
<ref id="B25">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gupta</surname> <given-names>G.</given-names></name> <name><surname>Ndiaye</surname> <given-names>A.</given-names></name> <name><surname>Filteau</surname> <given-names>M.</given-names></name></person-group> (<year>2021</year>). <article-title>Leveraging experimental strategies to capture different dimensions of microbial interactions</article-title>. <source>Front. Microbiol.</source> <volume>12</volume>:<fpage>700752</fpage>. <pub-id pub-id-type="doi">10.3389/fmicb.2021.700752</pub-id><pub-id pub-id-type="pmid">34646243</pub-id></citation></ref>
<ref id="B26">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hanada</surname> <given-names>R. E.</given-names></name> <name><surname>Pomella</surname> <given-names>A. W. V.</given-names></name> <name><surname>Costa</surname> <given-names>H. S.</given-names></name> <name><surname>Bezerra</surname> <given-names>J. L.</given-names></name> <name><surname>Loguercio</surname> <given-names>L. L.</given-names></name> <name><surname>Pereira</surname> <given-names>J. O.</given-names></name></person-group> (<year>2010</year>). <article-title>Endophytic fungal diversity in <italic>Theobroma cacao</italic> (cacao) and <italic>T. grandiflorum</italic> (cupua&#x000E7;u) trees and their potential for growth promotion and biocontrol of black-pod disease</article-title>. <source>Fungal Biol.</source> <volume>114</volume>, <fpage>901</fpage>&#x02013;<lpage>910</lpage>. <pub-id pub-id-type="doi">10.1016/j.funbio.2010.08.006</pub-id><pub-id pub-id-type="pmid">21036333</pub-id></citation></ref>
<ref id="B27">
<citation citation-type="web"><person-group person-group-type="author"><name><surname>Juo</surname> <given-names>A. S. R.</given-names></name> <name><surname>Franzluebbers</surname> <given-names>K.</given-names></name></person-group> (<year>2003</year>). <source>Tropical Soils: Properties and Management for Sustainable Agriculture</source>. <publisher-loc>New York, NY</publisher-loc>: <publisher-name>Oxford University Press</publisher-name>. Available online at: <ext-link ext-link-type="uri" xlink:href="https://books.google.com/books?hl=en&#x00026;lr=&#x00026;id=dVGgQcuoibAC&#x00026;oi=fnd&#x00026;pg=PA3&#x00026;dq=franzluebbers&#x0002B;tropical&#x0002B;soil&#x0002B;sustainable&#x00026;ots=I-fie2C3ed&#x00026;sig=C2BfBRHlUTEqHVI1HCUEC1WVHEc&#x00023;v=onepage&#x00026;q=franzluebbers">https://books.google.com/books?hl=en&#x00026;lr=&#x00026;id=dVGgQcuoibAC&#x00026;oi=fnd&#x00026;pg=PA3&#x00026;dq=franzluebbers&#x0002B;tropical&#x0002B;soil&#x0002B;sustainable&#x00026;ots=I-fie2C3ed&#x00026;sig=C2BfBRHlUTEqHVI1HCUEC1WVHEc&#x00023;v=onepage&#x00026;q=franzluebbers</ext-link> tropical soil sustainable&#x00026;f=false (accessed June 25, 2021).<pub-id pub-id-type="pmid">31639144</pub-id></citation></ref>
<ref id="B28">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kloepper</surname> <given-names>J. W.</given-names></name> <name><surname>Ryu</surname> <given-names>C.-M.</given-names></name> <name><surname>Zhang</surname> <given-names>S.</given-names></name></person-group> (<year>2004</year>). <article-title>Induced systemic resistance and promotion of plant growth by Bacillus sp</article-title>. <source>Phytopathology</source> <volume>94</volume>, <fpage>1259</fpage>&#x02013;<lpage>1266</lpage>. <pub-id pub-id-type="doi">10.1094/PHYTO.2004.94.11.1259</pub-id><pub-id pub-id-type="pmid">18944464</pub-id></citation></ref>
<ref id="B29">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Kuhn</surname> <given-names>M.</given-names></name></person-group> (<year>2020</year>). <source>caret: Classification and Regression Training</source>. <pub-id pub-id-type="pmid">28207752</pub-id></citation></ref>
<ref id="B30">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Leite</surname> <given-names>H. A. C.</given-names></name> <name><surname>Silva</surname> <given-names>A. B.</given-names></name> <name><surname>Gomes</surname> <given-names>F. P.</given-names></name> <name><surname>Gramacho</surname> <given-names>K. P.</given-names></name> <name><surname>Faria</surname> <given-names>J. C.</given-names></name> <name><surname>De Souza</surname> <given-names>J. T.</given-names></name> <etal/></person-group>. (<year>2013</year>). <article-title>Bacillus subtilis and Enterobacter cloacae endophytes from healthy <italic>Theobroma cacao</italic> L. trees can systemically colonize seedlings and promote growth</article-title>. <source>Appl. Microbiol. Biotechnol.</source> <volume>97</volume>, <fpage>2639</fpage>&#x02013;<lpage>2651</lpage>. <pub-id pub-id-type="doi">10.1007/s00253-012-4574-2</pub-id><pub-id pub-id-type="pmid">23212670</pub-id></citation></ref>
<ref id="B31">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname> <given-names>J.</given-names></name> <name><surname>Abdelfattah</surname> <given-names>A.</given-names></name> <name><surname>Norelli</surname> <given-names>J.</given-names></name> <name><surname>Burchard</surname> <given-names>E.</given-names></name> <name><surname>Schena</surname> <given-names>L.</given-names></name> <name><surname>Droby</surname> <given-names>S.</given-names></name> <etal/></person-group>. (<year>2018</year>). <article-title>Apple endophytic microbiota of different rootstock/scion combinations suggests a genotype-specific influence</article-title>. <source>Microbiome</source> <volume>6</volume>:<fpage>18</fpage>. <pub-id pub-id-type="doi">10.1186/s40168-018-0403-x</pub-id><pub-id pub-id-type="pmid">29374490</pub-id></citation></ref>
<ref id="B32">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lugtenberg</surname> <given-names>B.</given-names></name> <name><surname>Kamilova</surname> <given-names>F.</given-names></name></person-group> (<year>2009</year>). <article-title>Plant-growth-promoting rhizobacteria</article-title>. <source>Annu. Rev. Microbiol.</source> <volume>63</volume>, <fpage>541</fpage>&#x02013;<lpage>556</lpage>. <pub-id pub-id-type="doi">10.1146/annurev.micro.62.081307.162918</pub-id><pub-id pub-id-type="pmid">19575558</pub-id></citation></ref>
<ref id="B33">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Marasco</surname> <given-names>R.</given-names></name> <name><surname>Rolli</surname> <given-names>E.</given-names></name> <name><surname>Fusi</surname> <given-names>M.</given-names></name> <name><surname>Michoud</surname> <given-names>G.</given-names></name> <name><surname>Daffonchio</surname> <given-names>D.</given-names></name></person-group> (<year>2018</year>). <article-title>Grapevine rootstocks shape underground bacterial microbiome and networking but not potential functionality</article-title>. <source>Microbiome</source> <volume>6</volume>:<fpage>3</fpage>. <pub-id pub-id-type="doi">10.1186/s40168-017-0391-2</pub-id><pub-id pub-id-type="pmid">29298729</pub-id></citation></ref>
<ref id="B34">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mehlich</surname> <given-names>A.</given-names></name></person-group> (<year>1984</year>). <article-title>Mehlich 3 soil test extractant: a modification of mehlich 2 extractant</article-title>. <source>Commun. Soil Sci. Plant Anal.</source> <volume>15</volume>, <fpage>1409</fpage>&#x02013;<lpage>1416</lpage>. <pub-id pub-id-type="doi">10.1080/00103628409367568</pub-id></citation>
</ref>
<ref id="B35">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Motamayor</surname> <given-names>J. C.</given-names></name> <name><surname>Lachenaud</surname> <given-names>P.</given-names></name> <name><surname>da Silva E Mota</surname> <given-names>J. W.</given-names></name> <name><surname>Loor</surname> <given-names>R.</given-names></name> <name><surname>Kuhn</surname> <given-names>D. N.</given-names></name> <name><surname>Brown</surname> <given-names>J. S.</given-names></name> <etal/></person-group>. (<year>2008</year>). <article-title>Geographic and genetic population differentiation of the Amazonian chocolate tree (<italic>Theobroma cacao</italic> L)</article-title>. <source>PLoS ONE</source> <volume>3</volume>:<fpage>e3311</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0003311</pub-id><pub-id pub-id-type="pmid">18827930</pub-id></citation></ref>
<ref id="B36">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nascimento</surname> <given-names>A. D.</given-names></name> <name><surname>Lima</surname> <given-names>M. O.</given-names></name> <name><surname>Feij&#x000F3;</surname> <given-names>F. M.</given-names></name> <name><surname>J&#x000FA;nior</surname> <given-names>J. H.</given-names></name> <name><surname>Sobrinho</surname> <given-names>R. R.</given-names></name> <name><surname>Assun&#x000E7;&#x000E3;o</surname> <given-names>I. P.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>First report of colletotrichum aeschynomenes causing anthracnose in cacao (<italic>Theobroma cacao</italic>) in Brazil</article-title>. <source>Plant Dis.</source> <volume>103</volume>, <fpage>3284</fpage>&#x02013;<lpage>3284</lpage>. <pub-id pub-id-type="doi">10.1094/PDIS-11-18-2047-PDN</pub-id><pub-id pub-id-type="pmid">30812563</pub-id></citation></ref>
<ref id="B37">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nguyen</surname> <given-names>N. H.</given-names></name> <name><surname>Song</surname> <given-names>Z.</given-names></name> <name><surname>Bates</surname> <given-names>S. T.</given-names></name> <name><surname>Branco</surname> <given-names>S.</given-names></name> <name><surname>Tedersoo</surname> <given-names>L.</given-names></name> <name><surname>Menke</surname> <given-names>J.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>FUNGuild: an open annotation tool for parsing fungal community datasets by ecological guild</article-title>. <source>Fungal Ecol.</source> <volume>20</volume>, <fpage>241</fpage>&#x02013;<lpage>248</lpage>. <pub-id pub-id-type="doi">10.1016/j.funeco.2015.06.006</pub-id></citation>
</ref>
<ref id="B38">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Oksanen</surname> <given-names>J.</given-names></name> <name><surname>Blanchet</surname> <given-names>F. G.</given-names></name> <name><surname>Friendly</surname> <given-names>M.</given-names></name> <name><surname>Kindt</surname> <given-names>R.</given-names></name> <name><surname>Legendre</surname> <given-names>P.</given-names></name> <name><surname>Mcglinn</surname> <given-names>D.</given-names></name> <etal/></person-group>. (<year>2020</year>). <source>Package &#x0201C;vegan&#x0201D; Title Community Ecology Package Version 2.5-7</source>.</citation>
</ref>
<ref id="B39">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ploetz</surname> <given-names>R. C.</given-names></name></person-group> (<year>2006</year>). <article-title>Fusarium-induced diseases of tropical perennial crops fusarium-induced diseases of tropical, perennial crops</article-title>. <source>Phytopathology</source> <volume>96</volume>, <fpage>648</fpage>&#x02013;<lpage>652</lpage>. <pub-id pub-id-type="doi">10.1094/PHYTO-96-0648</pub-id><pub-id pub-id-type="pmid">18943183</pub-id></citation></ref>
<ref id="B40">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>P&#x000F5;lme</surname> <given-names>S.</given-names></name> <name><surname>Abarenkov</surname> <given-names>K.</given-names></name> <name><surname>Henrik Nilsson</surname> <given-names>R.</given-names></name> <name><surname>Lindahl</surname> <given-names>B. D.</given-names></name> <name><surname>Clemmensen</surname> <given-names>K. E.</given-names></name> <name><surname>Kauserud</surname> <given-names>H.</given-names></name> <etal/></person-group>. (<year>2020</year>). <article-title>FungalTraits: a user-friendly traits database of fungi and fungus-like stramenopiles</article-title>. <source>Fungal Divers</source> <volume>105</volume>, <fpage>1</fpage>&#x02013;<lpage>16</lpage>. <pub-id pub-id-type="doi">10.1007/s13225-020-00466-2</pub-id></citation>
</ref>
<ref id="B41">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Poudel</surname> <given-names>R.</given-names></name> <name><surname>Jumpponen</surname> <given-names>A.</given-names></name> <name><surname>Kennelly</surname> <given-names>M. M.</given-names></name> <name><surname>Rivard</surname> <given-names>C. L.</given-names></name> <name><surname>Gomez-Montano</surname> <given-names>L.</given-names></name> <name><surname>Garrett</surname> <given-names>K. A.</given-names></name></person-group> (<year>2019</year>). <article-title>Rootstocks shape the rhizobiome: rhizosphere and endosphere bacterial communities in the grafted tomato system</article-title>. <source>Appl. Environ. Microbiol.</source> <volume>85</volume>, <fpage>e01765</fpage>-<lpage>18</lpage>. <pub-id pub-id-type="doi">10.1128/AEM.01765-18</pub-id><pub-id pub-id-type="pmid">30413478</pub-id></citation></ref>
<ref id="B42">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Prawoto</surname> <given-names>A.</given-names></name></person-group> (<year>1990</year>). <article-title>A Study on cocoa (<italic>Theobroma cacao</italic> L.) budding V. The influence of rootstocks on scion-yield quality</article-title>. <source>Pelita Perkeb.</source> <volume>6</volume>, <fpage>38</fpage>&#x02013;<lpage>46</lpage>. <pub-id pub-id-type="doi">10.22302/iccri.jur.pelitaperkebunan.v6i2.366</pub-id></citation>
</ref>
<ref id="B43">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Quast</surname> <given-names>C.</given-names></name> <name><surname>Pruesse</surname> <given-names>E.</given-names></name> <name><surname>Yilmaz</surname> <given-names>P.</given-names></name> <name><surname>Gerken</surname> <given-names>J.</given-names></name> <name><surname>Schweer</surname> <given-names>T.</given-names></name> <name><surname>Yarza</surname> <given-names>P.</given-names></name> <etal/></person-group>. (<year>2013</year>). <article-title>The SILVA ribosomal RNA gene database project: improved data processing and web-based tools</article-title>. <source>Nucleic Acids Res.</source> <volume>41</volume>, <fpage>D590</fpage>&#x02013;<lpage>D596</lpage>. <pub-id pub-id-type="doi">10.1093/nar/gks1219</pub-id><pub-id pub-id-type="pmid">23193283</pub-id></citation></ref>
<ref id="B44">
<citation citation-type="web"><person-group person-group-type="author"><collab>R Core Team</collab></person-group> (<year>2020</year>). <source>R: A Language and Environment for Statistical Computing</source>. Available online at: <ext-link ext-link-type="uri" xlink:href="https://www.r-project.org/">https://www.r-project.org/</ext-link> (accessed March 30, 2021).</citation>
</ref>
<ref id="B45">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ribeiro</surname> <given-names>M. A. Q.</given-names></name> <name><surname>Da Silva</surname> <given-names>J O</given-names></name> <name><surname>Aitken</surname> <given-names>W. M.</given-names></name> <name><surname>Machado</surname> <given-names>R. C. R.</given-names></name> <name><surname>Baligar</surname> <given-names>V. C.</given-names></name></person-group> (<year>2008</year>). <article-title>Nitrogen use efficiency in cacao genotypes</article-title>. <source>J. Plant Nutr.</source> <volume>31</volume>, <fpage>239</fpage>&#x02013;<lpage>249</lpage>. <pub-id pub-id-type="doi">10.1080/01904160701853720</pub-id></citation>
</ref>
<ref id="B46">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rivas</surname> <given-names>R.</given-names></name> <name><surname>Vel&#x000E1;zquez</surname> <given-names>E.</given-names></name> <name><surname>Willems</surname> <given-names>A.</given-names></name> <name><surname>Vizca&#x000ED;no</surname> <given-names>N.</given-names></name> <name><surname>Subba-Rao</surname> <given-names>N. S.</given-names></name> <name><surname>Mateos</surname> <given-names>P. F.</given-names></name> <etal/></person-group>. (<year>2002</year>). <article-title>A new species of Devosia that forms a unique nitrogen-fixing root-nodule symbiosis with the aquatic legume Neptunia natans (L.f.) Druce</article-title>. <source>Appl. Environ. Microbiol.</source> <volume>68</volume>, <fpage>5217</fpage>&#x02013;<lpage>5222</lpage>. <pub-id pub-id-type="doi">10.1128/AEM.68.11.5217-5222.2002</pub-id><pub-id pub-id-type="pmid">12406707</pub-id></citation></ref>
<ref id="B47">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rosmana</surname> <given-names>A.</given-names></name> <name><surname>Hikmawati</surname> <given-names>H.</given-names></name> <name><surname>Asman</surname> <given-names>A.</given-names></name></person-group> (<year>2013</year>). <article-title>Identification of a disease on cocoa caused by Fusarium in Sulawesi</article-title>. <source>Pelita Perkeb.</source> <volume>29</volume>, <fpage>210</fpage>&#x02013;<lpage>219</lpage>. <pub-id pub-id-type="doi">10.22302/iccri.jur.pelitaperkebunan.v29i3.13</pub-id></citation>
</ref>
<ref id="B48">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sarkhosh</surname> <given-names>A.</given-names></name> <name><surname>Shahkoomahally</surname> <given-names>S.</given-names></name> <name><surname>Asis</surname> <given-names>C.</given-names></name> <name><surname>McConchie</surname> <given-names>C.</given-names></name></person-group> (<year>2021</year>). <article-title>Influence of rootstocks on scion leaf mineral content in mango tree (<italic>Mangifera indica</italic> L.)</article-title>. <source>Hortic. Environ. Biotechnol.</source> <volume>62</volume>, <fpage>725</fpage>&#x02013;<lpage>735</lpage>. <pub-id pub-id-type="doi">10.1007/s13580-021-00355-w</pub-id></citation>
</ref>
<ref id="B49">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Segata</surname> <given-names>N.</given-names></name> <name><surname>Izard</surname> <given-names>J.</given-names></name> <name><surname>Waldron</surname> <given-names>L.</given-names></name> <name><surname>Gevers</surname> <given-names>D.</given-names></name> <name><surname>Miropolsky</surname> <given-names>L.</given-names></name> <name><surname>Garrett</surname> <given-names>W. S.</given-names></name> <etal/></person-group>. (<year>2011</year>). <article-title>Metagenomic biomarker discovery and explanation</article-title>. <source>Genome Biol.</source> <volume>12</volume>:<fpage>R60</fpage>. <pub-id pub-id-type="doi">10.1186/gb-2011-12-6-r60</pub-id><pub-id pub-id-type="pmid">21702898</pub-id></citation></ref>
<ref id="B50">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tchameni</surname> <given-names>S. N.</given-names></name> <name><surname>Ngonkeu</surname> <given-names>M. E. L.</given-names></name> <name><surname>Begoude</surname> <given-names>B. A. D.</given-names></name> <name><surname>Wakam Nana</surname> <given-names>L.</given-names></name> <name><surname>Fokom</surname> <given-names>R.</given-names></name> <name><surname>Owona</surname> <given-names>A. D.</given-names></name> <etal/></person-group>. (<year>2011</year>). <article-title>Effect of <italic>Trichoderma asperellum</italic> and arbuscular mycorrhizal fungi on cacao growth and resistance against black pod disease</article-title>. <source>Crop Prot.</source> <volume>30</volume>, <fpage>1321</fpage>&#x02013;<lpage>1327</lpage>. <pub-id pub-id-type="doi">10.1016/j.cropro.2011.05.003</pub-id></citation>
</ref>
<ref id="B51">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tchameni</surname> <given-names>S. N.</given-names></name> <name><surname>Nwaga</surname> <given-names>D.</given-names></name> <name><surname>Wakam</surname> <given-names>L. N.</given-names></name> <name><surname>Mangaptche Ngonkeu</surname> <given-names>E. L.</given-names></name> <name><surname>Fokom</surname> <given-names>R.</given-names></name> <name><surname>Kuat&#x000E9;</surname> <given-names>J.</given-names></name> <etal/></person-group>. (<year>2012</year>). <article-title>Growth enhancement, amino acid synthesis and reduction in susceptibility towards <italic>Phytophthora megakarya</italic> by arbuscular mycorrhizal fungi inoculation in cocoa plants</article-title>. <source>J. Phytopathol.</source> <volume>160</volume>, <fpage>220</fpage>&#x02013;<lpage>228</lpage>. <pub-id pub-id-type="doi">10.1111/j.1439-0434.2012.01888.x</pub-id></citation>
</ref>
<ref id="B52">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wemheuer</surname> <given-names>F.</given-names></name> <name><surname>Berkelmann</surname> <given-names>D.</given-names></name> <name><surname>Wemheuer</surname> <given-names>B.</given-names></name> <name><surname>Daniel</surname> <given-names>R.</given-names></name> <name><surname>Vidal</surname> <given-names>S.</given-names></name> <name><surname>Daghela</surname> <given-names>H. B. B.</given-names></name></person-group> (<year>2020</year>). <article-title>Agroforestry management systems drive the composition, diversity, and function of fungal and bacterial endophyte communities in theobroma cacao leaves</article-title>. <source>Microorganisms</source> <volume>8</volume>:<fpage>405</fpage>. <pub-id pub-id-type="doi">10.3390/microorganisms8030405</pub-id><pub-id pub-id-type="pmid">32183118</pub-id></citation></ref>
<ref id="B53">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>White</surname> <given-names>T. J.</given-names></name> <name><surname>Bruns</surname> <given-names>T.</given-names></name> <name><surname>Lee</surname> <given-names>S.</given-names></name> <name><surname>Taylor</surname> <given-names>J. W.</given-names></name></person-group> (<year>1990</year>). <article-title>Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics</article-title>, in <source>PCR Protocols: A Guide to Methods and Applications</source>, eds <person-group person-group-type="editor"><name><surname>Innis</surname> <given-names>M. A.</given-names></name> <name><surname>Gelfand</surname> <given-names>D. H.</given-names></name> <name><surname>Sninsky</surname> <given-names>J. J.</given-names></name> <name><surname>White</surname> <given-names>T. J.</given-names></name></person-group> (<publisher-loc>New York, NY</publisher-loc>: <publisher-name>Academic Press, Inc.</publisher-name>), <fpage>315</fpage>&#x02013;<lpage>322</lpage>.</citation>
</ref>
<ref id="B54">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yilmaz</surname> <given-names>P.</given-names></name> <name><surname>Parfrey</surname> <given-names>L. W.</given-names></name> <name><surname>Yarza</surname> <given-names>P.</given-names></name> <name><surname>Gerken</surname> <given-names>J.</given-names></name> <name><surname>Pruesse</surname> <given-names>E.</given-names></name> <name><surname>Quast</surname> <given-names>C.</given-names></name> <etal/></person-group>. (<year>2014</year>). <article-title>The SILVA and &#x0201C;all-species Living Tree Project (LTP)&#x0201D; taxonomic frameworks</article-title>. <source>Nucleic Acids Res.</source> <volume>42</volume>, <fpage>D643</fpage>&#x02013;<lpage>D648</lpage>. <pub-id pub-id-type="doi">10.1093/nar/gkt1209</pub-id><pub-id pub-id-type="pmid">24293649</pub-id></citation></ref>
<ref id="B55">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yin</surname> <given-names>J. P. T.</given-names></name></person-group> (<year>2004</year>). <article-title>Rootstock effects on cocoa in Sabah, Malaysia</article-title>. <source>Exp. Agric.</source> <volume>40</volume>, <fpage>445</fpage>&#x02013;<lpage>452</lpage>. <pub-id pub-id-type="doi">10.1017/S0014479704002108</pub-id><pub-id pub-id-type="pmid">30886898</pub-id></citation></ref>
<ref id="B56">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zanne</surname> <given-names>A. E.</given-names></name> <name><surname>Abarenkov</surname> <given-names>K.</given-names></name> <name><surname>Afkhami</surname> <given-names>M. E.</given-names></name> <name><surname>Aguilar-Trigueros</surname> <given-names>C. A.</given-names></name> <name><surname>Bates</surname> <given-names>S.</given-names></name> <name><surname>Bhatnagar</surname> <given-names>J. M.</given-names></name> <etal/></person-group>. (<year>2020</year>). <article-title>Fungal functional ecology: bringing a trait-based approach to plant-associated fungi</article-title>. <source>Biol. Rev.</source> <volume>95</volume>, <fpage>409</fpage>&#x02013;<lpage>433</lpage>. <pub-id pub-id-type="doi">10.1111/brv.12570</pub-id><pub-id pub-id-type="pmid">31763752</pub-id></citation></ref>
</ref-list>
</back>
</article>
