# Cultivated and wild pearl millet display contrasting patterns of abundance and co-occurrence in their root mycobiome

Marie-Thérèse Mofini, Abdala Diedhiou, Marie Simonin, Donald Tchouomo

Dondjou, Sarah Pignoly, Cheikh Ndiaye, Doohong Min, Yves Vigouroux,

Laurent Laplaze, Aboubacry Kane

► To cite this version:

Marie-Thérèse Mofini, Abdala Diedhiou, Marie Simonin, Donald Tchouomo Dondjou, Sarah Pignoly, et al. Cultivated and wild pearl millet display contrasting patterns of abundance and co-occurrence in their root mycobiome. Scientific Reports, 2022, 12 (1), pp.207. 10.1038/s41598-021-04097-8. hal-03576717

HAL Id: hal-03576717

https://hal.inrae.fr/hal-03576717

Submitted on 12 Jul 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

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

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

OPEN

# Cultivated and wild pearl millet display contrasting patterns of abundance and co-occurrence in their root mycobiome

Marie-Thérèse Mofini<fnref n="1" />,²,³,⁴,⁵, Abdala G. Diedhiou<fnref n="1" />,²,³,⁴✉, Marie Simonin⁶,⁹, Donald Tchouomo Dondjou<fnref n="1" />,²,³,⁴,⁵, Sarah Pignoly²,³,⁷, Cheikh Ndiaye<fnref n="1" />,²,³,⁴, Doohong Min⁸, Yves Vigouroux⁷, Laurent Laplaze²,³,⁴,⁷✉ & Aboubacry Kane<fnref n="1" />,²,³,⁴✉

Fungal communities associated with roots play a key role in nutrient uptake and in mitigating the abiotic and biotic stress of their host. In this study, we characterized the roots mycobiome of wild and cultivated pearl millet [Pennisetum glaucum (L.) R. Br., synonym: Cenchrus americanus (L.) Morrone] in three agro-ecological areas of Senegal following a rainfall gradient. We hypothesized that wild pearl millet could serve as a reservoir of endophytes for cultivated pearl millet. We therefore analyzed the soil factors influencing fungal community structure and whether cultivated and wild millet shared the same fungal communities. The fungal communities associated with pearl millet were significantly structured according to sites and plant type (wild vs cultivated). Besides, soil pH and phosphorus were the main factors influencing the fungal community structure. We observed a higher fungal diversity in cultivated compared to wild pearl millet. Interestingly, we detected higher relative abundance of putative pathotrophs, especially plant pathogen, in cultivated than in wild millet in semi-arid and semi-humid zones, and higher relative abundance of saprotrophs in wild millet in arid and semi-humid zones. A network analysis based on taxa co-occurrence patterns in the core mycobiome revealed that cultivated millet and wild relatives had dissimilar groups of hub taxa. The identification of the core mycobiome and hub taxa of cultivated and wild pearl millet could be an important step in developing microbiome engineering approaches for more sustainable management practices in pearl millet agroecosystems.

Climate change poses new threats to food security as the human population continues to grow<fnref n="1" />, and the price of fertilizers continues to rise. A paradigm shift is needed to address future challenges faced by agriculture. One such potential change would be to harness the potential of plant microbiomes and especially root-associated fungi to improve plant nutrition and health. Root-associated fungi are a diverse group of microorganisms that develop on the surface or in the tissues of the host. Through the interaction with their host, some of these fungi play beneficial roles that contribute to the health of their host. For instance, they can produce bioactive secondary metabolites that can protect their host against pathogen and insect attacks and enhance abiotic stress tolerance²⁻⁴. They also promote the water and minerals absorption by their host⁵. On the other hand, some of these fungi, especially pathogens, can strongly reduce their host productivity under certain conditions⁶,⁷.

<footnote n="1">¹Département de Biologie Végétale, Faculté des Sciences et Techniques, Université Cheikh Anta Diop (UCAD), BP 5005, Dakar Fann, Senegal. ²Laboratoire Mixte International Adaptation des Plantes et Microorganismes associés aux Stress Environnementaux (LAPSE), Centre de recherche de Bel-Air, Dakar, Sénégal. ³Laboratoire Commun de Microbiologie (LCM), Centre de Recherche de Bel-Air, Dakar, Sénégal. ⁴Centre d'Excellence Africain en Agriculture pour la Sécurité Alimentaire et Nutritionnelle (CEA-AGRISAN), UCAD, Dakar, Sénégal. ⁵Centre d'Etude Régional pour l'Amélioration de l'Adaptation à la Sécheresse (CERAAS), Institut Sénégalais de Recherches Agricoles (ISRA), Route de Khombole, Thiès, Sénégal. ⁶IPME, IRD, Cirad, Université de Montpellier, Montpellier, France. ⁷DIADE, Université de Montpellier, IRD, Cirad, 911 Avenue Agropolis, 34394 Montpellier cedex 5, France. ⁸Department of Agronomy, Kansas State University, Manhattan, KS, USA. ⁹Present address: Université d'Angers, Institut Agro, INRAE, IRHS, SFR 4207 QuaSaV, 49000 Angers, France.  email: abdala.diedhiou@ucad.edu.sn; laurent.laplaze@ird.fr; aboubacry.kane@ucad.edu.sn</footnote>
In cereals, root-associated fungi belong to taxonomically diverse groups including mycorrhizal fungi, mainly arbuscular mycorrhizal fungi (AMF), and non-mycorrhizal fungi¹⁻¹¹. The latter group is dominated by endophytic fungi which have different lifestyles (mutualistic, latent pathogen and latent saprophyte) depending on host genotype and physiology (senescence, flowering, fruition, vegetative period and age¹²⁻¹⁴). According to Nguyen et al.¹⁵ these root-associated fungi can be classified into three main groups based on their trophic modes: (1) symbiotroph, receiving nutrients by exchanging resources with host cells; (2) pathotroph, receiving nutrients by harming host cells; and (3) saprotroph, receiving nutrients by breaking down dead host cells. It is important to note, however, that many fungal genera contain more than one trophic strategy¹⁶. These fungal genera have been considered as multi-trophic mode groups¹⁷. Through their different trophic modes, the root-associated fungi can interact directly or indirectly with each other positively or negatively, and thereby influence the fitness of their host plant¹⁸⁻²⁰. On the other hand, it has been established that host genotype and environment factors shape microbiome assembly in natural environments²⁰⁻²².

Pearl millet [Pennisetum glaucum (L.) R. Br., synonym: Cenchrus americanus (L.) Morrone] is a major staple food and source of fodder and fuel in the arid and semi-arid regions of sub-Saharan Africa and India. The vegetative, reproductive and physiological characteristics of pearl millet make it a crop well suited for growth under difficult conditions, including low soil fertility, high pH, low soil moisture, high temperature, high salinity and limited rainfall where other cereals such as maize, rice, sorghum or durum wheat would fail²³. Pearl millet grain content is highly nutritious, with 8–19% protein, low starch levels, high fiber content (1.2 g/100 g²⁴), and higher micronutrient concentrations (iron and zinc) than rice, wheat, maize and sorghum²⁵. It was domesticated in the central Sahel (Mali-Niger) about 4900 years ago as corroborated by archaeological and genomic studies²⁶.

In Senegal, cultivated pearl millet still coexists with wild millet around or even inside some farmers' plots. Gene flows with wild relatives in the western and eastern Sahel are important and have contributed to increase the diversity of cultivated pearl millet in Africa²⁶. However, the comparison of the microbiomes of cultivated and wild plants has not yet been performed. We hypothesized that wild pearl millet could act as a potential microbial reservoir for cultivated plants. In addition, understanding the complex fungal community assemblage of cultivated and wild millet could help to develop more sustainable management practices in pearl millet agroecosystems.

Here, we employed DNA metabarcoding techniques to characterize fungal communities associated with cultivated and wild pearl millet roots in three agro-ecological zones of Senegal following a rainfall gradient. We examined the influence of soil parameters on the structure of fungal communities and identified the pearl millet core mycobiome common to wild and cultivated plants across the three zones. Identification of this core mycobiome might provide information on fungi potentially involved in maintaining community stability and which may play an important role in plant health and productivity. Network analysis was then performed to investigate co-occurrence patterns and identify potential hub taxa in the core mycobiome of cultivated and wild plants. Once again, this approach should provide potential target hub groups for further studies and for engineering approaches.

## Results

**Experimental sites have contrasted soil properties.** In order to analyze the fungal communities associated with millet roots and the factors driving their structure, we first characterized the soil properties of the experimental plots from three agro-ecological zones [Darou-Mousty (arid zone), Dya (semi-arid zone) and Nioro (semi-humid zone)] in Senegal. Statistically significant differences in soil properties [pH, total nitrogen (N), total phosphorus (P), carbon to nitrogen ratio (C:N) and ammonium (NH₄⁺)] were found between the three zones (Supplementary Table S1). When the properties of soils collected under cultivated plants were compared to those of soils collected under wild plants, we observed significant differences for the carbon:nitrogen (C:N) ratio which was higher in soil collected under wild plants in Dya (p = 0.036), and Nioro (p = 0.022), total P which was higher in soil collected under cultivated plants in Nioro (p = 0.002), and pH H₂O which was lower in soil collected under cultivated plants in Nioro (p = 0.005).

**Wild and cultivated pearl millet have different taxonomic diversity and composition of root-associated fungal communities.** We used metabarcoding targeting the internal transcribed spacer 1 (ITS1) region to characterize the fungal communities associated with cultivated and wild millet roots. After sequencing, the normalized dataset accounted for 5524 fungal operational taxonomic units (OTUs) for wild and cultivated pearl millet. The site had no effect on fungal species richness, Shannon index and Simpson diversity index (Fig. 1A–C, Supplementary Table S2). However, the plant type (cultivated vs wild) had a significant effect on all the indices studied (species richness, *p* = 0.000; Shannon, *p* = 0.000; Simpson, *p* = 0.009). In addition, the variance explained by the random factor (plot, R²c – R²m) was low compared to those explained by the fixed factors (plant type and site, R²m) for the different α-diversity indices (Supplementary Table S2). Hence, the various α-diversity indices revealed that cultivated pearl millet display a higher fungal diversity than wild pearl millet.

The non-metric multidimensional scaling (NMDS) allowed us to visualize the distribution of our samples in the two-dimensional space (stress value of 0.19). We observed that site and plant type and their interaction were the major factors structuring fungal communities in the data set, while the effects of plant type and site-plant type interaction were minor (Table 2, Supplementary Table S3, Fig. 2). The root fungal communities were distributed along the first two dimensions of the NMDS based on the site where their plant host inhabits, with a separation between the samples from Darou-Mousty and those from Nioro across the axis 1, and the samples from Dya and those from Nioro across the axis 2 (Fig. 2). The analysis indicated that variation in fungal community structure was mainly explained by soil pH and phosphorus (Fig. 2; Supplementary Table S4).

In Darou-Mousty, 19 OTUs were differentially abundant based on plant type, of which 9 OTUs were significantly enriched in wild plants and 10 enriched in cultivated plants (Fig. 3). In Nioro, 12 OTUs were enriched in
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Species richness (A) and Shannon diversity index (B) of the fungal communities associated with cultivated and wild pearl millet in the three studied sites (Darou-Mousty, Dya and Nioro). In the linear mixed effects (LME) model used to test the effect of plant type and site, plot was included as a random factor (see Table 1).</figcaption></figure>
<table><caption>**Table 1.** Results from ANOVA of the linear mixed effects (LME) model testing the effect of plant type, site and their interaction on species richness, Shannon and Simpson diversity indexes of the fungal communities associated with cultivated and wild Pearl millet. R²m (marginal r squared) represents the variance explained by the fixed factors, and R²c (conditional r squared) represents the variance explained by the both fixed and random factors.</caption><thead><tr><th>Factors</th><th>numDF</th><th>denDF</th><th>F</th><th>p</th><th>R<sup>2</sup>m</th><th>R<sup>2</sup>c</th></tr></thead><tbody><tr><td colspan="7"><strong>Richness</strong></td></tr><tr><td>Intercept</td><td>1</td><td>51</td><td>6450.426</td><td>&lt;0.0001</td><td>0.310</td><td>0.310</td></tr><tr><td>Plants</td><td>1</td><td>51</td><td>19.153</td><td>0.000</td><td></td><td></td></tr><tr><td>Sites</td><td>2</td><td>3</td><td>0.936</td><td>0.483</td><td></td><td></td></tr><tr><td>Plants:Sites</td><td>2</td><td>51</td><td>2.736</td><td>0.074</td><td></td><td></td></tr><tr><td colspan="7"><strong>Shannon</strong></td></tr><tr><td>Intercept</td><td>1</td><td>51</td><td>5988.246</td><td>&lt;0.0001</td><td>0.216</td><td>0.258</td></tr><tr><td>Plants</td><td>1</td><td>51</td><td>14.182</td><td>0.000</td><td></td><td></td></tr><tr><td>Sites</td><td>2</td><td>3</td><td>0.755</td><td>0.542</td><td></td><td></td></tr><tr><td>Plants:Sites</td><td>2</td><td>51</td><td>0.321</td><td>0.727</td><td></td><td></td></tr><tr><td colspan="7"><strong>Simpson</strong></td></tr><tr><td>Intercept</td><td>1</td><td>51</td><td>32,490.240</td><td>&lt;0.0001</td><td>0.162</td><td>0.209</td></tr><tr><td>Plants</td><td>1</td><td>51</td><td>7.470</td><td>0.009</td><td></td><td></td></tr><tr><td>Sites</td><td>2</td><td>3</td><td>1.100</td><td>0.439</td><td></td><td></td></tr><tr><td>Plants:Sites</td><td>2</td><td>51</td><td>0.530</td><td>0.591</td><td></td><td></td></tr></tbody></table>

<table><caption>**Table 2.** Summary of non-parametric permutational multivariate analysis of variance (PERMANOVA) based on Bray-Curtis distance to test the effects of site and plant type on the structure of fungal communities associated to pearl millet roots. *df* degrees of freedom, *SS* sum of squares, *MS* mean sum of squares, *F model F* statistics, *R*² partial R-squared, based on 999 permutations.</caption><thead><tr><th>Factors</th><th>df</th><th>SS</th><th>MS</th><th>F.Model</th><th>R<sup>2</sup></th><th>p value</th></tr></thead><tbody><tr><td>Sites</td><td>2</td><td>3.1885</td><td>1.59425</td><td>7.9021</td><td>0.21097</td><td>0.001</td></tr><tr><td>Plants</td><td>1</td><td>0.4361</td><td>0.43606</td><td>2.1614</td><td>0.02885</td><td>0.001</td></tr><tr><td>Sites × plants</td><td>2</td><td>0.5947</td><td>0.29737</td><td>1.4739</td><td>0.03935</td><td>0.005</td></tr><tr><td>Residuals</td><td>54</td><td>10.8945</td><td>0.20175</td><td>0.72083</td><td></td><td></td></tr><tr><td>Total</td><td>59</td><td>15.1138</td><td>1.00000</td><td></td><td></td><td></td></tr></tbody></table>

wild plants and 9 were enriched in cultivated plants whereas only 3 OTUs were significantly enriched in cultivated plants in Dya (Fig. 3). On the other hand, only the genera *Fusarium* and *Chaetomium* had differentially abundant OTUs in more than one site (Darou-Mousty and Nioro). The mean abundances and log twofold change values of differentially abundant OTUs in each site and plant type are given in Supplementary Table S5.

**Wild and cultivated pearl millet show differences in hosted fungal functional groups.** A total of 4048 out of 5524 OTUs (73%) were assigned to functional guilds using the FUNGuild database. Among them, 2990 OTUs (54%) classified as either pathotroph, saprotroph or symbiotroph (Fig. 4A), while 1058 OTUs (19%) belong to fungi with more than one trophic strategy (Fig. 4B). Putative pathotrophs (38% of relative abundance) and saprotrophs (28% of relative abundance) were the dominant groups in terms of relative abundance in cultivated and wild millet taken together. They were followed by saprotroph-symbiotroph (14% of relative abundance) and pathotroph-saprotroph-symbiotroph (11% of relative abundance).

The relative abundance of putative pathotrophs and those of saprotrophs was significantly influenced by plant type and site-plant type interaction (LME, Supplementary Table S6). Putative pathotrophs had a higher relative abundance in cultivated than in wild pearl millet in Dya (*p* = 0.038) and Nioro (*p* < 0.0001) sites. By contrast, saprotrophs had a higher relative abundance in wild than in cultivated millet in Darou-Mousty (*p* = 0.000) and Nioro (*p* = 0.003) sites (Fig. 4A). However, the relative abundance of symbiotrophs was not influenced by plant type, site or their interaction (LME, Supplementary Table S6). Yet, we observed significant site-plant type interaction effects on the relative abundance of pathotroph-saprotroph and those of pathotroph-saprotroph-symbiotroph (Supplementary Table S6, Fig. 4B).

At a finer scale, putative plant pathogens (32% of relative abundance), undefined saprotrophs (17% of relative abundance) and clavicipitaceous endophyte-saprophitic fungi (13% of relative abundance) were the dominant root-associated fungal guilds of cultivated and wild millet taken together. The effects of site and plant type on the 15 most abundant fungal guilds are shown in Supplementary Table S7. The relative abundance of putative plant pathogens was significantly influenced by plant type, site and site-plant type interaction, while those of undefined saprotrophs was only influenced by plant type (Supplementary Table S7). Specifically, undefined
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Non-Metric Multidimensional Scale (NMDS) plots depicting the similarity of fungal communities according to sites and plant types, and the relative importance of soil properties (arrows) explaining the variation in fungal communities according to sites. Only significant factors (*p* < 0.05) correlated with fungal community dissimilarity are presented. Each point represents a single sample.</figcaption></figure>

saprotroph had a higher relative abundance (*p* = 0.005) in wild than in cultivated millet. Meanwhile, putative plant pathogens had a higher relative abundance in cultivated millet than in their relative wild only in Nioro (*p* = 0.004). In addition, the relative abundance of clavicipitaceous endophyte-saprophitic fungi was not influenced by plant type, site or their interaction (Supplementary Table S7). Altogether, our data indicate that cultivated pearl millet hosts a higher abundance of putative fungal pathogens in semi-arid and semi-humid zones, while wild pearl millet is associated with a higher abundance of fungal saprotrophs in arid and semi-humid zones.

Pearson's correlations between soil properties and the relative abundance of trophic mode groups revealed contrasting patterns within and across sites (Fig. 5). In Darou-Mousty and Dya for instance, the relative abundance of putative pathotrophs was negatively correlated with total P, C/N, total N and total C, while the relative abundance of saprothrophs was positively correlated with these soil properties. In Nioro, the relative abundance of saprotrophs was positively correlated with pH and C/N, and negatively correlated with total P, assimilable P, NO₃⁻ and NH₄⁺, while the relative abundance of pathotrophs followed the opposite pattern for these soil properties (Fig. 5).

**Pearl millet core mycobiome.** In order to determine the core mycobiome of pearl millet roots, we used a 75% prevalence threshold to identify the taxa shared between the majority of samples that were collected at sites that were hundreds of kilometers apart and had different soil characteristics. The core mycobiome of millet roots contained 260 OTUs that represented only 4.7% of the observed OTU richness of the entire fungal community. Interestingly, these 260 OTUs accounted for 28.8% of the total sequences of the dataset (Fig. 6A). Of these 260 core OTUs, 91 taxa were found in all of the samples. The 40 most abundant core OTUs (Supplementary Table S9) had relative abundances ranging from 0.17 to 2.04% in the entire dataset. These core taxa belonged to four fungal phyla: Ascomycota (217 OTUs), Basidiomycota (33 OTUs), Chytridiomycota (2 OTUs) and Glomeromycota (8 OTUs). Ascomycota, the most dominant, represented from 63.8 to 97% of the core OTUs with 22 families (Fig. 6C) in pearl millet roots, both in cultivated and wild plants. The most abundant families were *Pleosporaceae* (31 OTUs) and *Nectriaceae* (31 OTUs) with relative abundance of 4.9% and 4.46%, respectively. The second most dominant phylum was Basidiomycota representing from 1.9 to 35.9% of the core OTUs with seven families. In contrast, the phyla Chytridiomycota and Glomeromycota were less represented in the core mycobiome (0.01–1.06% and 0.09–1.28% respectively; Fig. 6B). Some OTUs were classified as OTUs with a very high level of affiliation (Sordariales, Ascomycota, Basidiomycota, Sebacinales, Dothideomycetes, Pleosporales, Xylariales and Hypocreales). All the seven trophic mode groups were represented in the core mycobiome (Supplementary Table S9).

**Co-occurrence patterns in the core mycobiome of cultivated and wild pearl millet.** We further analyzed the differences in the core mycobiome of cultivated and wild pearl millet by investigating their taxa
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Differentially abundant OTUs detected in cultivated and wild Pearl millet by pairwise comparison (DeSeq2 analysis of 5524 OTUs, with *p* value adjusted to the 1% threshold). Each point represents an individual OTU, which was assigned to genus (y-axis, 63 OTUs) and phylum level (colours). In the x-axis, positive values of log2 fold change indicate higher abundance of OTUs in wild plants and negative values indicate higher abundance of OTUs in cultivated plants.</figcaption></figure>

co-occurrence patterns using ecological network analysis. The correlation-based network of cultivated plants contained 214 nodes and 807 edges, while that of the wild plants had 195 nodes and 680 edges (Fig. 7A,B). The networks of cultivated and wild plants had 84.51% and 84.12% of strong positive correlations, respectively. Except for average degree, all estimated network-level topological features (average path length, network diameter, graph density, modularity and clustering coefficient) were higher in wild than in cultivated pearl millet (Supplementary Table S10). The betweenness centrality and eigenvector centrality were higher in the network of cultivated plants compared to that of wild plants, while node degree and closeness centrality did not differ significantly between the two networks (Supplementary Fig. S1). In addition, we identified four potential hubs belonging to three genera which include plant pathogen (Bipolaris, Cochliobolus and Curvularia) and one dung saprotroph-undefined saprotroph-wood saprotroph (Penicillium #1) in the ecological network of cultivated pearl millet. By contrast, four potential hubs belonging to two dung saprotroph-undefined saprotroph-wood saprotroph (Penicillium #1 and Penicillium #2), one clavic pitaceous endophyte-saprophitic fungi (Paecilomyces) and one fungal parasite-plant pathogen (Helminthosporium) were identified in the ecological network of wild pearl millet (Fig. 7C). The two ecological networks shared only one potential hub (Penicillium #1). Furthermore, these seven potential hub taxa belonged to three trophic mode groups including saprotroph, pathotroph and saprotroph-symbiotroph (Fig. 7C).

Discussion

Here, we analyzed the root mycobiome of cultivated and wild pearl millet and the environmental factors struc-
turing its diversity across three agro-ecological zones in Senegal. Our results show that root fungal communities
were mainly structured by sites and plant type. Soil pH and phosphorus were the main factors explaining the
site effect on pearl millet fungal communities as reported in previous studies²⁷⁻³¹. We observed that, in each
site, the relationship between the relative abundance of putative pathotrophs and soil properties (especially total
P, C/N, total N and total C) and the relationship between the relative abundance of saprotrophs and those soil
properties followed opposite trends. On the other hand, contrasting patterns were observed across sites. This
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Relative OTU abundance of fungal functional groups associated with cultivated and wild millet across the three sites (Darou-Mousty, Dya and Nioro). The abundances of fungi with one trophic strategy are given in (A), and those of fungi with more than one trophic strategy (saprotroph-symbiotroph, pathotroph-symbiotroph, pathotroph-saprotroph, pathotroph-saprotroph-symbiotroph) are in (B). In the linear mixed effects (LME) model used to test the effect of plant type and site, plot was included as a random factor (see Table 3).</figcaption></figure>
<table><caption>**Table 3.** Results from ANOVA of the linear mixed effects (LME) model testing the effect of plant type, site and their interaction on the relative abundance of fungal functional groups associated with cultivated and wild millet across the three sites (Darou-Mousty, Dya and Nioro). R<sup>2</sup>m (marginal r squared) represents the variance explained by the fixed factors, and R<sup>2</sup>c (conditional r squared) represents the variance explained by the both fixed and random factors.</caption><thead><tr><th>Factors</th><th>numDF</th><th>denDF</th><th>F</th><th>p</th><th>R<sup>2</sup>m</th><th>R<sup>2</sup>c</th></tr></thead><tbody><tr><td colspan="7"><strong>Pathotrophs</strong></td></tr><tr><td>Intercept</td><td>1</td><td>51</td><td>193.340</td><td>&lt;0.0001</td><td>0.653</td><td>0.811</td></tr><tr><td>Plants</td><td>1</td><td>51</td><td>40.770</td><td>&lt;0.0001</td><td></td><td></td></tr><tr><td>Sites</td><td>2</td><td>3</td><td>7.331</td><td>0.07</td><td></td><td></td></tr><tr><td>Plants:sites</td><td>2</td><td>51</td><td>12.901</td><td>&lt;0.0001</td><td></td><td></td></tr><tr><td colspan="7"><strong>Saprotrophs</strong></td></tr><tr><td>Intercept</td><td>1</td><td>51</td><td>114.884</td><td>&lt;0.0001</td><td>0.393</td><td>0.560</td></tr><tr><td>Plants</td><td>1</td><td>51</td><td>31.554</td><td>&lt;0.0001</td><td></td><td></td></tr><tr><td>Sites</td><td>2</td><td>3</td><td>0.808</td><td>0.524</td><td></td><td></td></tr><tr><td>Plants:sites</td><td>2</td><td>51</td><td>6.710</td><td>0.003</td><td></td><td></td></tr><tr><td colspan="7"><strong>Symbiotrophs</strong></td></tr><tr><td>Intercept</td><td>1</td><td>51</td><td>7.335</td><td>0.009</td><td>0.184</td><td>0.740</td></tr><tr><td>Plants</td><td>1</td><td>51</td><td>2.350</td><td>0.131</td><td></td><td></td></tr><tr><td>Sites</td><td>2</td><td>3</td><td>0.821</td><td>0.520</td><td></td><td></td></tr><tr><td>Plants:sites</td><td>2</td><td>51</td><td>1.301</td><td>0.281</td><td></td><td></td></tr><tr><td colspan="7"><strong>Pathotrophs_Saprotrophs</strong></td></tr><tr><td>Intercept</td><td>1</td><td>51</td><td>11.818</td><td>0.001</td><td>0.224</td><td>0.465</td></tr><tr><td>Plants</td><td>1</td><td>51</td><td>0.012</td><td>0.913</td><td></td><td></td></tr><tr><td>Sites</td><td>2</td><td>3</td><td>1.631</td><td>0.332</td><td></td><td></td></tr><tr><td>Plants:sites</td><td>2</td><td>51</td><td>3.379</td><td>0.042</td><td></td><td></td></tr><tr><td colspan="7"><strong>Pathotrophs_Symbiotrophs</strong></td></tr><tr><td>Intercept</td><td>1</td><td>51</td><td>17.482</td><td>0.0001</td><td>0.053</td><td>0.053</td></tr><tr><td>Plants</td><td>1</td><td>51</td><td>1.610</td><td>0.210</td><td></td><td></td></tr><tr><td>Sites</td><td>2</td><td>3</td><td>0.771</td><td>0.537</td><td></td><td></td></tr><tr><td>Plants:sites</td><td>2</td><td>51</td><td>0.079</td><td>0.924</td><td></td><td></td></tr><tr><td colspan="7"><strong>Saprotrophs_Symbiotrophs</strong></td></tr><tr><td>Intercept</td><td>1</td><td>51</td><td>46.181</td><td>&lt;0.0001</td><td>0.413</td><td>0.557</td></tr><tr><td>Plants</td><td>1</td><td>51</td><td>0.023</td><td>0.879</td><td></td><td></td></tr><tr><td>Sites</td><td>2</td><td>3</td><td>5.822</td><td>0.093</td><td></td><td></td></tr><tr><td>Plants:sites</td><td>2</td><td>51</td><td>2.724</td><td>0.075</td><td></td><td></td></tr><tr><td colspan="7"><strong>Pat_Sap_Sym</strong></td></tr><tr><td>Intercept</td><td>1</td><td>51</td><td>40.341</td><td>&lt;0.0001</td><td>0.189</td><td>0.287</td></tr><tr><td>Plants</td><td>1</td><td>51</td><td>0.771</td><td>0.384</td><td></td><td></td></tr><tr><td>Sites</td><td>2</td><td>3</td><td>1.744</td><td>0.314</td><td></td><td></td></tr><tr><td>Plants:sites</td><td>2</td><td>51</td><td>3.292</td><td>0.045</td><td></td><td></td></tr></tbody></table>

finding contrasts with others that showed that abundance of pathogenic fungi followed the same trend as that of saprotrophic fungi<sup>32,33</sup>. Although we relate the site effect observed on the relative abundance of most guilds and trophic mode groups to differences in soil properties, we cannot rule out the contribution of other factors such as host plant density, temperature, rainfall and moisture<sup>34–37</sup>. We acknowledge that the fungal trophic modes were identified by FUNGuild which is based on fungal taxonomy, and therefore studies about the fungal lifestyles and trophic modes are needed to better understand the factors influencing the contrasting patterns we observed. In addition, the differences in soil properties observed between and within sites could be related to different factors including the quantity and type of inputs used, cropping system, land use history and plant functional trait effects<sup>38–41</sup>.

We observed a higher fungal diversity in cultivated compared to wild pearl millet. We speculate that during pearl millet domestication, genotypic and phenotypic changes resulting from selection and demographic bottlenecks could have affected root traits related to nutrient foraging and plant–microbe interaction in response to resource availability and heterogeneity<sup>42–44</sup>. In addition, we suggest that additional adaptive changes in the root morphological traits and quantitative and/or qualitative root exudate composition of domesticated pearl millet due to constant clogging effects over time may have facilitated colonization by a broader range of fungi<sup>45–48</sup>. Interestingly, we detected higher relative abundance of putative pathotrophs, especially potential plant pathogens, in cultivated pearl millet than in wild pearl millet in semi-arid and semi-humid zones, while having higher relative
<figure><img src="image_7.png" /></figure>

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

<figure><img src="image_9.png" /><figcaption>**Figure 5.** Pearson's correlation between soil properties and the relative abundance of the 7 trophic mode groups of cultivated and wild millet across the three sites. Stars indicate significant correlation: \*p<0.05, \*\*p<0.01, \*\*\*p<0.001. The color key indicates the Pearson correlation coefficient values. For the trophic mode groups, Pat_Sym, Pat_Sap_Sym, Pat_Sap and Sap_Sym refer to pathotroph-symbiotroph, pathotroph-saprotroph-symbiotroph, pathotroph-saprotroph and saprotroph-symbiotroph, respectively. For the soil properties, Ass P and tot P refer to assimilable and total phosphorus respectively. The correlation coefficient r and p-values used to produce heatmap figures are given in Supplementary Table S8.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>B</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>C</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 6.** Fungal community composition of the core mycobiome. (A) core mycobiome diversity vs all taxa diversity (left) and core microbiome relative abundance vs all taxa (right), (B) relative abundance of core mycobiome phyla, (C) 30 most abundant family of the core mycobiome.</figcaption></figure>

abundance of saprotrophs in wild pearl millet in arid and semi-humid zones. This suggests that cultivated pearl millet plants were more colonized by pathotrophs than wild relatives, while the latter provided a more appropriate ecological niche for saprotrophs across the three studied agro-ecological zones. Besides genotypic and phenotypic differences between cultivated millet and wild relatives, resource demand, exudate and litter quality, interaction with other microbes, agriculture practices, and moisture may account for this observation⁹,⁴²,⁴⁹. Although further studies are needed to better understand factors behind the selection of fungal communities by wild and cultivated pearl millet, the fact that the relative abundance of putative pathotrophs was lower in wild millet suggests that it could be a potential resource of genetic resistance to root pathogens⁵⁰,⁵¹. On the other hand, although we believe that our results reflected intrinsic differences between cultivated and wild pearl millet in terms of diversity of their root-associated fungi, they should be interpreted cautiously since our sequencing depth was insufficient to cover the entire fungal diversity in both cultivated and wild pearl millet (Supplementary Fig. S2).

We hypothesized that wild pearl millet could serve as a reservoir of endophytes for cultivated pearl millet. We therefore characterized a pearl millet core mycobiome across three contrasted regions. Interestingly, it was composed of diverse groups including 260 taxa that accounted for almost 28.8% of total abundance of fungal community in our dataset indicating that these taxa are both prevalent and abundant on millet roots. In our study, Ascomycota were dominant in the core mycobiome. Basidiomycota was the second most dominant phylum, and Glomeromycota was the least represented in the core mycobiome. The lower proportion of Glomeromycota could be related to the low levels of pearl millet root colonization by AMF (root length colonization ≤ 5%⁵²). The ecological network analysis we performed revealed a predominance of strong positive correlations between
<figure><img src="image_13.png" /><figcaption>**Figure 7.** Co-occurrence-based network of the core mycobiome of cultivated (A) and wild (B) pearl millet. Each node corresponds to an OTU, and edges correspond to either positive (yellow) or negative (red) correlations inferred from OTU abundances. Node size reflects their proportional abundance (Pro_abundance) and their color reflects their trophic mode. Nodes with fungal taxon names represent the potential hub OTUs of cultivated and wild plants inferred from node degree and betweenness centrality values (C). Symbols after taxon names indicate their guilds: Ø = plant pathogens, # = dung saprotroph-undefined saprotroph-wood saprotrophs, § = clavicipitaceous endophyte-saprophitic fungi, and ¶ = fungal parasite-plant pathogens. Dashed lines indicate the threshold estimated of the top 2% of degree and betweenness centrality values.</figcaption></figure>
fungal OTUs in both cultivated and wild millet (84.51% and 84.12%, respectively), suggesting a potential for extensive cooperative interactions or sharing similar ecological niches between most fungal taxa in pearl millet roots<sup>53,54</sup>. Our results also showed that wild plants had lower betweenness centrality, but higher network-level topological features (e.g. network diameter, graph density, modularity and clustering coefficient) compared to cultivated plants, indicating a more complex and connected network<sup>55,56</sup>. This higher network complexity in wild pearl millet may result in a more stable mycobiome that could contribute to higher plant resilience to environmental perturbations<sup>57,58</sup>. On the other hand, plant genotype and/or agriculture practices may account for the lower network complexity in cultivated pearl millet<sup>59,60</sup>.

We identified four different potential hub taxa in the core mycobiome of cultivated (*Bipolaris*, *Cochliobolus*, *Curvularia* and *Penicillium* #1) and wild (*Helminthosporium*, *Paecilomyces*, *Penicillium* #1 and *Penicillium* #2) plants. These hub taxa may act as keystone taxa<sup>56,61</sup> in the mycobiome with a strong structuring or recruiting role in the community. Surprisingly, the cultivated millet and wild relatives had dissimilar groups of potential hub taxa (only one was shared), suggesting that they may develop dissimilar interactions with their mycobiome for biotic and abiotic stress tolerance<sup>48,62</sup>. Similarly, analyses of seed microbiome of domesticated and wild rice showed that they had different hub OTUs<sup>61</sup>. In our study, two potential hub taxa including the shared one, belonged to *Penicillium*, that was also identified as hub taxon in the seed microbiome of wild rice<sup>61</sup>. Besides its role in nutrient cycling through the decomposition of organic matter, *Penicillium* is known to solubilize phosphate<sup>63-65</sup> and to excrete antimicrobial substances<sup>66,67</sup>. In this respect, Murali and Amruthesh<sup>68</sup> have shown that *Penicillium oxalicum*, significantly reduced the downy mildew disease and enhanced plant growth in pearl millet. We also identified *Paecilomyces*, a clavicipitaceous endophyte-saprophytic fungus, as a potential hub taxon in the core mycobiome of wild pearl millet. *Paecilomyces* fungi are known to produce bioactive metabolites and therefore have been used as plant growth-promoting fungi for horticultural crops<sup>69</sup>. Moreover, it has been reported that extracts of some *Paecilomyces* strains combined with urea or phosphate increased the yield of rice and maize<sup>70,71</sup>. Most of clavicipitaceous fungi are members of dark septate endophyte (DSE) which can be in some cases more frequent than mycorrhizal fungi<sup>5,72,73</sup>. They play a significant ecological and physiological role in different ecosystems by impacting plant growth and nutrition. In the pearl millet core mycobiome, we identified at least six families (*Cladosporiaceae*, *Herpotrichiellaceae*, *Didymellaceae*, *Didymosphaeriaceae*, *Leptosphaeriaceae*, *Pleosporaceae* and *Trichosphaeriaceae*) and two unassigned taxa at family level (Sordariales and Xylariales), containing DSE fungi. The high relative abundance and ubiquity of members of DSE suggest a mutualistic lifestyle that might compensate for low levels of AMF colonization in pearl millet<sup>74</sup>. Interestingly, the other potential hub taxa identified in this study and assigned as plant pathogens (*Bipolaris*, *Cochliobolus* and *Curvularia*) and fungal parasite-plant pathogen (*Helminthosporium*), have also been reported to be phosphate-solubilizing fungi<sup>75-79</sup>. Their impact on the fitness of cultivated and wild pearl millet remains to be elucidated.

In conclusion, this study allowed us to assess the factors influencing fungal communities associated with wild and cultivated pearl millet roots. We found that soil pH and available phosphorus were the main factors explaining the effect of sites on fungal communities. On the other hand, we observed a higher fungal diversity in cultivated compared to wild pearl millet. Interestingly, our data indicate that cultivated pearl millet hosts a higher relative abundance of potential fungal pathogens while wild pearl millet is associated with a higher relative abundance of fungal saprotrophs. In addition, the cultivated millet and wild relatives had dissimilar groups of hub taxa (only one out of eight, was shared), suggesting that they may develop dissimilar sophisticated dialogues with their mycobiome for biotic and abiotic stress tolerance. Future research should focus on these core and hub taxa, which are likely to play an important role in pearl millet fitness and will allow the development of microbiome engineering approaches for agriculture.

## Methods

**Site description.** The study was conducted in Senegal, West Africa, in the North–South region of the groundnut basin (13° 45' N, 15° 47' W) at 18 m above sea level, specifically in the district of Darou-Mousty (15° 02' 31" N, 16° 02' 53" W) with two plots, Dya (14° 23' 60" N, 16° 10' W, 11 m altitude) with three plots and Nioro (13° 45' N, 15° 48' W, 19 m altitude) with one plot (Supplementary Fig. S3). The two main crops grown in these areas are millet (*P. glaucum*) and groundnuts (*Arachis hypogaea*). The mean annual precipitation is 750 mm and mainly comes between July and September, and the annual average temperature is 30 °C<sup>80</sup>. The soil is a Deck-Dior<sup>81</sup> silty-sand fine sandy, mixed haplic ferric lixisol<sup>82</sup>, a leached ferruginous tropical soil. Top soil (0–10 cm) has a sand content > 90%, organic matter and total nitrogen content of 0.52 and 0.03% respectively, total P content of 70 mg kg<sup>-1</sup> and an average pH (water) of 6.2<sup>83</sup>. The sites were chosen on the basis of soil homogeneity (Dior type) following a gradient of rainfall from North to South.

**Experimental design and soil/root sampling.** The experiment was conducted during the 2016 rainy season, with a mean precipitation of 476, 518 and 911 mm in Darou-Mousty (arid zone), Dya (semi-arid zone) and Nioro (semi-humid zone), respectively. Field trials were conducted in collaboration with local farmers' associations and sampling was performed with the help of farmers. Permission to collect plants and soils was obtained from farmers before settling the trials. Trials did not involve endangered or protected species and complied with relevant national and international guidelines and legislation. Seeds from a single pearl millet variety (Souna 3) were provided to all farmers participating to the trials on all sites in order to limit plant genotype effect. The farmers followed their traditional agricultural practices (fertilization or not, previous crops, etc.) that were recorded as well as the cropping history and the vegetation around (Supplementary Table S11).

At the end of the grain-filling stage, soil and root samples were taken. From each cultivated plot, a set of five replicates was collected. Each biological replicate for sequencing analysis came from the roots of five plants harvested from each cultivated plot. Sampled plants were located 10 m from each other to cover the whole plot.
Soils and roots of wild relatives were collected within or around the cultivated plots. For each individual plant, approximately 500 g of soil influenced by the roots was collected at a depth of 0–20 cm. A total of 60 root samples and 60 soil samples were collected: two plant types (cultivated and wild), six plots and five replicates per plot. Soils and roots were placed in plastic bags in ice and transported to the laboratory where they were stored at 4 °C for 24 h before processing. Total DNA extraction was only done from the roots while the soils samples were used for physicochemical analyses. Wild pearl millet is morphologically dissimilar to cultivated pearl millet. Supplementary Figure S4 shows the differences between cultivated and wild pearl millet.

**Soil chemical analysis.** Soil pH was determined with a soil-to-water ratio of 1: 2.5. Total carbon (C) and total nitrogen (N) contents were quantified using Elemental Analyzer (Flash EA 1112 series, ThermoFinnigan, France). Soil nitrate ($NO_3^-$) and ammonium ($NH_4^+$) were extracted with 2 M KCl and were quantified by Bran + Luebbe GmbH AutoAnalyzer 3. Soil available phosphorus (AP) was extracted using sodium bicarbonate and then measured by the molybdenum-blue method⁸⁴. The P concentration was determined after dry mineralization by inductively coupled plasma atomic emission spectrometry (ICP-AES)⁸⁵.

**DNA-metabarcoding.** Root samples were cleaned of soil particles by manual agitation and about 100 mg of fine roots from each replicate were collected and ground with 270 mg of Fontainebleau sand (minimum quantity allowing the crushing of 100 mg of roots) with a mortar and pestle. 185 mg of the obtained grind (corresponding to approximately 50 mg of root material) was used for DNA extraction. No specific treatment was applied to the roots, so the mycobiome considered here included both fungi inhabiting the root surface and root endosphere of pearl millet. Total DNA was extracted using the DNeasy Plant Mini Kit (Qiagen') according to the manufacturer's instructions. A total of 60 DNA samples (30 cultivated and 30 wild) were obtained. The concentration and purity of isolated DNA were determined using a Nanodrop ND-2000 UV-VIS spectrophotometer (NanoDrop Technologies, Wilmington) and DNA samples were stored at -20 °C before sequencing.

DNA amplification and sequencing of fungal rDNA were performed at MR DNA (www.mrdnalab.com, Shallowater, TX, USA). Briefly, the ITS1 region of fungal rDNA was amplified by PCR using 5 µL of Q5 Reaction Buffer (5x), 5 µL of Q5 GC high Enhancer (5x), 1 µL (10 µM) of ITS1F primer (5'-CTTGGTCATTTAGAGGAA GTAA-3';⁸⁶), 1 µL (10 µM) of ITS2 primer (5'-GCTGCGTTCTTCATCGATGC-3';⁸⁷), 2 µL of dNTP (2.5 µM), 1 µL of DNA template (20 ng µL⁻¹), 0.25 µL of Q5 Polymerase (5U µL⁻¹) and 9.75 µL of ddH₂O to a final volume of 25 µL. PCR was performed in triplicate using the following conditions: 5 min at 98 °C for initial denaturing, followed by 27 cycles of 98 °C for 30 s, primer annealing at 56 °C for 30 s and extension at 72 °C for 30 s, followed by a final extension for 5 min at 72 °C.

After amplification, the quality and relative concentration of the amplicons were checked by migration on 2% agarose gel. Multiple replicates were pooled together in equal proportions based on their molecular weight and DNA concentrations. Pooled DNA samples were purified using calibrated Ampure XP beads. Then the pooled and purified amplicons were used to prepare DNA libraries following Illumina Truseq DNA library preparation protocol. Sequencing was performed on a MiSeq Illumina platform (2x300) following the manufacturer's guidelines.

**Sequence analysis.** First, raw Illumina MiSeq paired-end reads were assembled using MR DNA pipeline for ITS1 region fungal libraries. Subsequently, sequences were demultiplexed and formatted for processing using a Python script (http://drive5.com/usearch/manual/uparse_pipeline.html). Next, fungal sequences were separately quality-filtered and clustered into operational taxonomic units (OTUs) using UPARSE pipeline and UPARSE algorithm⁸⁸. Briefly, sequences were quality-filtered allowing a maximum e value of 0.5. Subsequently, reads were trimmed to 240-bp (base pairs) length as well as dereplicated and sorted by abundance, removing singletons (sequences which appeared once) prior OTU determination at 97% sequence similarity threshold. Then, chimeric sequences were detected and removed using UCHIME⁸⁹ and Gold database as reference. Finally, reads from the entire dataset were mapped back to the representative fungal databases to generate one OTU table. The taxonomic affiliation of each OTU was obtained using the UNITE database (version 7.2, https://unite.ut.ee,⁹⁰). OTUs were then assigned to functional groups using the FUNGuild database (https://github.com/UMNFuN/FUNGuild¹⁵). We only accepted the guild assignment that confidence rankings were “highly probable” or “probable”. For taxa with more than one trophic strategy, we subsequently used manual evaluation as recommended by Nilsson et al.¹⁶. The raw sequence data has been deposited in figshare (https://figshare.com/articles/Pearl_Millet_Fungus/11277950).

**Data analyses.** The number of sequences per sample varied from 25,556 to 154,517. Because the library sizes were very uneven, the data were normalized to the same number of counts per sample (25,556 sequences) using the “rarefy” function in the “Vegan” package⁹¹ in R (v 4.1.0; https://cran.r-project.org) as recommended by Weiss et al.⁹². Species richness, Shannon and Simpson diversity indices were used to assess α-diversity of fungal communities associated with cultivated and wild pearl millet. Before testing the differences between the means, we first visualized the distribution of the data using a histogram to determine the distribution of errors. Then, the Shapiro test was used to check normality followed by the Bartlett test to check for homoscedasticity of variances. As normality was verified, the data were subjected to a linear mixed-effects (LME) model fit by restricted maximum likelihood (REML) with plant type, site, and their interaction as fixed factors and plot as random factor using the “lme” function in the “nlme” package⁹³. The significance of fixed effects was assessed using the “anova.lme” function. If there was a significant effect of the site or the site-plant type interaction, pairwise comparisons were conducted using the “emmeans” package with Tukey’s adjusted p values⁹⁴. The variance explained by the fixed factors (marginal R²) and those explained by both the fixed and random factors (conditional R²)
were calculated with the “r.squaredGLMM” function in the “MuMln” package⁹⁵. Differences in the community structure between samples (β-diversity) were visualized using a non-metric multidimensional scaling (NMDS) based on Bray–Curtis distances calculated from the "metaMDS" function implemented in the "Vegan" package⁹². The “betadisper” function in “Vegan” was used to compare group dispersions. The "envfit" function was used to determine the relationships between community structure and soil properties. The PERMANOVA test (non-parametric Permutational Multivariate Analysis of Variance) was used to test significant differences in the structure of fungal communities by the "adonis" function⁹⁶.

Differentially abundant fungal OTUs between plant types (cultivated vs wild pearl millet) were identified in each site by the “DESeq2” package⁹⁷ using the default values of the "Test (Wald)" and "FitType (parametric)" options of the function. The differential abundance calculation is based on a modeling of the OTU distribution by a negative binomial law. The adjusted *p* value cutoff was set to alpha = 0.01, and the list of OTUs declared differentially abundant was extracted. Considering only the assigned OTUs using FUNGuild, we calculated the relative abundance of each guild and trophic mode group, and a LME model fit by REML was used to assess the effects of plant type, site, and their interaction. The model included plant type, site, and their interaction as fixed factors and plot as random factor. The significance of fixed effects was assessed using the “anova.lme” function. If there was a significant effect of the site or the site-plant type interaction, pairwise comparisons were conducted as described above. The variance explained by the fixed factors and those explained by both fixed and random factors were calculated with the “r.squaredGLMM” function. For each site, we used Pearson’s correlation coefficients to determine the relationship between soil properties and the relative abundance of each trophic mode group using the “rcorr” function in the "Hmisc" package⁹⁸, and a heatmap was then built using the "heatmap.2" function in the "gplots" package⁹⁹.

We characterized the core fungal microbiome (mycobiome) of pearl millet roots on a 75% prevalence threshold using the "microbiome" package¹⁰⁰ to identify highly prevalent taxa on pearl millet roots that are present in the majority of samples (across all sites and plant types). We then constructed co-occurrence networks to infer interactions among OTUs in the core mycobiome of cultivated plants and wild plants, based on Spearman's correlation inferred from OTU abundances. We considered only positive correlations (with r > 0.6) and negative correlations (with r < -0.6) associated with FDR-adjusted *p* values < 0.01⁵⁵,¹⁰¹. Gephi software (v0.9.2; https://gephi.org) was used to visualize the ecological networks and estimate node-level topological features (degree, betweenness centrality, closeness centrality and eigenvector centrality) and network-level topological features (average degree, average path length, network diameter, graph density, modularity and clustering coefficient) for the cultivated and wild plants⁵⁵. In each network, nodes correspond to OTUs, and edges correspond to strong correlations inferred from OTU abundances. For each ecological network, the OTUs belonging to the top 2% of degree and betweenness centrality were identified as potential hub OTUs⁶¹. We then used the Wilcoxon test to compare the estimated node-level topological features between cultivated and wild plants, and graphics were made using the 'ggplot' package¹⁰².

Received: 10 September 2021; Accepted: 15 December 2021

Published online: 07 January 2022

## References

1. FAO. Food and Agriculture Organization. The State of Food and Agriculture. Livestock in the balance (FAO, 2009).

2. Zhang, H. W., Song, Y. C. & Tan, R. X. Biology and chemistry of endophytes. *Nat. Prod. Rep.* **23**, 753–771 (2006).

3. Rodriguez, R. J. et al. Fungal endophytes: Diversity and functional roles. *New Phytol.* **182**, 314–330 (2009).

4. Khan, A. R. et al. Plant growth-promoting potential of endophytic fungi isolated from Solanum nigrum leaves. *World J. Microbiol. Biotechnol.* **31**(9), 1461–1466 (2015).

5. Jumpponen, A. Dark septate endophytes-are they mycorrhizal? *Mycorrhiza* **11**(4), 207–211 (2001).

6. Maron, J. L., Marler, M., Klironomos, J. N. & Cleveland, C. C. Soil fungal pathogens and the relationship between plant diversity and productivity. *Ecol. Lett.* **14**, 36–41 (2011).

7. van Ruijven, J., Ampt, E., Francioli, D. & Mommer, L. Do soil-borne fungal pathogens mediate plant diversity-productivity relationships? Evidence and future opportunities. *J Ecol.* **108**, 1810–1821 (2020).

8. Casini, G. et al. Endophytic fungal communities of ancient wheat varieties. *Phytopathol. Mediterr.* **58**(1), 151–162 (2019).

9. Xu, S. Q., Tian, L., Chang, C. L., Li, X. J. & Tian, C. J. Cultivated rice rhizomicrobiome is more sensitive to environmental shifts than that of wild rice in natural environments. *Appl. Soil Ecol.* **140**, 68–77 (2019).

10. Brisson, V. L., Schmidt, J. E., Northen, T. R., Vogel, J. P. & Gaudin, A. C. M. Impacts of maize domestication and breeding on rhizosphere microbial community recruitment from a nutrient depleted agricultural soil. *Sci. Rep.* **9**, 15611 (2019).

11. Gao, C. et al. Fungal community assembly in drought-stressed sorghum shows stochasticity, selection, and universal ecological dynamics. *Nat. Commun.* **11**, 34 (2020).

12. Promputtha, I. et al. A phylogenetic evaluation of whether endophytes become saprotrophs at host senescence. *Microb. Ecol.* **53**, 579–590 (2007).

13. Hyde, K. D. & Soytong, K. The fungal endophyte dilemma. *Fungal Diversity* **33**, 163–173 (2008).

14. Yuan, Z. L. et al. Identity, diversity, and molecular phylogeny of the endophytic mycobiota in the roots of rare wild rice (*Oryza granulate*) from a nature reserve in Yunnan, China. *Appl. Environ. Microbiol.* **76**, 1642–1652 (2010).

15. Nguyen, N. H. et al. FUNGuild: An open annotation tool for parsing fungal community datasets by ecological guild. *Fungal Ecol.* **20**, 241–248 (2016).

16. Nilsson, R. H. et al. Mycobiome diversity: High-throughput sequencing and identification of fungi. *Nat. Rev. Microbiol.* **17**, 95–109 (2019).

17. George, P. B. L. et al. Primer and database choice affect fungal functional but not biological diversity findings in a national soil survey. *Front. Environ. Sci.* **7**, 173 (2019).

18. Barea, J., Azcón, R. & Azcón-Aguilar, C. Mycorrhizosphere interactions to improve plant fitness and soil quality. *Antonie Van Leeuwenhoek* **81**, 343–351 (2002).

19. Hassani, M. A., Durán, P. & Hacquard, S. Microbial interactions within the plant holobiont. *Microbiome* **6**(1), 58 (2018).
20. de la Fuente Cantó, C. et al. An extended root phenotype: The rhizosphere, its formation and impacts on plant fitness. *Plant J.* **103**, 951–964 (2020).

21. Wagner, M. R. et al. Host genotype and age shape the leaf and root microbiomes of a wild perennial plant. *Nat. Commun.* **7**, 1–15 (2016).

22. Freedman, Z. B. et al. Environment-host-microbial interactions shape the *Sarracenia purpurea* microbiome at the continental scale. *Ecology* **102**(5), e03308 (2021).

23. Vadez, V. et al. Phenotyping pearl millet for adaptation to drought. *Front Physiol.* **3**, 386 (2012).

24. Nambiar, V. S. et al. Potential functional implications of pearl millet (*Pennisetum glaucum*) in health and disease. *J. Appl. Pharm. Sci.* **1**, 62–67 (2011).

25. Tako, E. et al. Higher iron pearl millet (*Pennisetum glaucum* L.) provides more absorbable iron that is limited by increased polyphenolic content. *Nutr. J.* **14**, 11 (2015).

26. Burgarella, C. et al. A Western Sahara Centre of domestication inferred from pearl millet genomes. *Nat. Ecol. Evol.* **2**, 1377–1380 (2018).

27. Tedersoo, L. et al. Global diversity and geography of soil fungi. *Science* **346**, 4168–4183 (2014).

28. Siles, J. A. & Margesin, R. Abundance and diversity of bacterial, archaeal, and fungal communities along an altitudinal gradient in alpine forest soils: What are the driving factors?. *Microb. Ecol.* **72**, 207–220 (2016).

29. Ding, J. et al. Influence of inorganic fertilizer and organic manure application on fungal communities in a long-term field experiment of Chinese Mollisols. *Appl. Soil Ecol.* **111**, 114–122 (2017).

30. Yu, P. et al. Root type and soil phosphate determine the taxonomic landscape of colonizing fungi and the transcriptome of field-grown maize roots. *New Phytol.* **217**, 1240–1253 (2018).

31. Li, Y. et al. Studies on fungal communities and functional guilds shift in tea continuous cropping soils by high-throughput sequencing. *Ann. Microbiol.* **70**, 7 (2020).

32. Veach, A. M., Stokes, C. E., Knoepp, J., Jumpponen, A. & Baird, R. Fungal communities and functional guilds shift along an elevational gradient in the southern Appalachian Mountains. *Microb. Ecol.* **76**, 156–168 (2018).

33. Zhao, A. et al. Soil fungal community is more sensitive to nitrogen deposition than increased rainfall in a mixed deciduous forest of China. *Soil Ecol. Lett.* **2**, 20–32 (2020).

34. Hiiesalu, I., Bahram, M. & Tedersoo, L. Plant species richness and productivity determine the diversity of soil fungal guilds in temperate coniferous forest and bog habitats. *Mol. Ecol.* **26**, 4846–4858 (2017).

35. Anthony, M. A., Frey, S. D. & Stinson, K. A. Fungal community homogenization, shift in dominant trophic guild, and appearance of novel taxa with biotic invasion. *Ecosphere* **8**(9), e01951. https://doi.org/10.1002/eecs2.1951 (2017).

36. Makiola, A. et al. Land use is a determinant of plant pathogen alpha—but not beta-diversity. *Mol. Ecol.* **28**, 3786–3798 (2019).

37. Cheng, Y. T., Zhang, L. & He, S. Y. Plant–microbe interactions facing environmental challenge. *Cell Host Microbe* **26**, 183–192 (2019).

38. McClintock, N. C. & Diop, A. M. Soil fertility management and compost use in Senegal's peanut basin. *Int. J. Agric. Sustain.* **3**, 79–91 (2005).

39. Hobbie, S. E. Plant species effects on nutrient cycling: Revisiting litter feedbacks. *Trends Ecol. Evol.* **30**, 357–363 (2015).

40. Faucon, M. P., Houben, D. & Lambers, H. Plant functional traits: Soil and ecosystem services. *Trends Plant Sci.* **22**, 385–394 (2017).

41. Laekemariam, F. & Kibret, K. Explaining soil fertility heterogeneity in smallholder farms of southern Ethiopia. *Appl. Environ. Soil Sci.* **20**, 6161059 (2020).

42. Perez-Jaramillo, J. E., Mendes, R. & Raaijmakers, J. M. Impact of plant domestication on rhizosphere microbiome assembly and functions. *Plant Mol. Biol.* **90**, 635–644 (2015).

43. Leff, J. W., Lynch, R. C., Kane, N. C. & Fierer, N. Plant domestication and the assembly of bacterial and fungal communities associated with strains of the common sunflower, *Helianthus annuus*. *New Phytol.* **214**, 412–423 (2016).

44. Schmidt, J. E., Bowles, T. M. & Gaudin, A. C. M. Using ancient traits to convert soil health into crop yield: Impact of selection on maize root and rhizosphere function. *Front. Plant Sci.* **7**, 1–11 (2016).

45. Berg, G. & Smalla, K. Plant species and soil type cooperatively shape the structure and function of microbial communities in the rhizosphere. *FEMS Microbiol. Ecol.* **68**(1), 1–13 (2009).

46. Badri, D. V., Chaparro, J. M., Zhang, R., Shen, Q. & Vivanco, J. M. Application of natural blends of phytochemicals derived from the root exudates of Arabidopsis to the soil reveal that phenolic related compounds predominantly modulate the soil microbiome. *J. Biol. Chem.* **288**, 4502–4512 (2013).

47. Lebeis, S. L. et al. Salicylic acid modulates colonization of the root microbiome by specific bacterial taxa. *Science* **349**, 860–864 (2015).

48. Lareen, A., Burton, F. & Schäfer, P. Plant root-microbe communication in shaping root microbiomes. *Plant Mol. Biol.* **90**, 575–587 (2016).

49. Nguyen, D. Q. et al. Soil and root nutrient chemistry structure root-associated fungal assemblages in temperate forests. *Environ. Microbiol.* **00**(00), 00–00 (2020).

50. Dempwolf, H. et al. Past and future use of wild relatives in crop breeding. *Crop Sci.* **57**(3), 1070–1082 (2017).

51. Radić, T. et al. Roots-associated community composition and co-occurrence patterns of fungi in wild grapevine. *Fugal Ecol.* **50**, 101034 (2021).

52. Mofini, M. T. et al. Mycorrhizal status of cultivated and wild pearl millet (*Pennisetum glaucum*) in three agro-ecological zones of Senegal. *IJDR* **10**(12), 42550–42556 (2020).

53. Zhang, B., Zhang, J., Liu, Y., Shi, P. & Wei, G. Co-occurrence patterns of soybean rhizosphere microbiome at a continental scale. *Soil Biol. Biochem.* **118**, 178–186 (2018).

54. Qian, X. et al. Leaf and root endospheres harbor lower fungal diversity and less complex fungal co-occurrence patterns than rhizosphere. *Front. Microbiol.* **10**, 1015 (2019).

55. Jiao, S., Chen, W. M. & Wei, G. H. Biogeography and ecological diversity patterns of rare and abundant bacteria in oil-contaminated soils. *Mol. Ecol.* **26**, 5305–5317 (2017).

56. Tipton, L. et al. Fungi stabilize connectivity in the lung and skin microbial ecosystems. *Microbiome* **6**, 12 (2018).

57. de Vries, F. T. et al. Soil bacterial networks are less stable under drought than fungal networks. *Nat. Commun.* **9**, 3033 (2018).

58. Wagg, C., Dudenhöffer, J. H., Widmer, F. & van der Heijden, M. G. Linking diversity, synchrony and stability in soil microbial communities. *Funct. Ecol.* **32**, 1280–1292 (2018).

59. Banerjee, S. et al. Agricultural intensification reduces microbial network complexity and the abundance of keystone taxa in roots. *ISME J.* **13**, 1722–1736 (2019).

60. Tian, L. et al. Comparison between wild and related cultivated rice species reveals strong impacts of crop domestication on methane metabolism of the rhizomicrobiome. *Res. Square* https://doi.org/10.21203/rs.3.rs-28865/v1 (2020).

61. Kim, H., Lee, K. K., Jeon, J., Harrie, W. A. J. & Ly, Y. H. Domestication of Oryza species eco-evolutionarily shapes bacterial and fungal communities in rice seed. *Microbiome* **8**, 20 (2020).

62. Badri, D. V., Weir, T. L., van der Lelie, D. & Vivanco, J. M. Rhizosphere chemical dialogues: Plant–microbe interactions. *Curr. Opin. Biotechnol.* **20**, 642–650 (2009).
63. Wakelin, S. A., Warren, R. A., Harvey, P. R. & Ryder, M. H. Phosphate solubilization by *Penicillium* spp. closely associated with wheat roots. *Biol. Fertil. Soils* **40**, 36–43 (2004).

64. Pandey, A. et al. Phosphate solubilization by *Penicillium* spp. isolated from soil samples of Indian Himalayan region. *World J. Microbiol Biotechnol.* **24**, 97–102 (2008).

65. Qiao, H. et al. The phosphate-solubilizing ability of *Penicillium guanacastense* and its effects on the growth of *Pinus massoniana* in phosphate-limiting conditions. *Biology Open* **8**, e046797 (2019).

66. Yang, L. et al. Antifungal substances produced by *Penicillium oxalicum* strain PY-1-potential antibiotics against plant pathogenic fungi. *World J. Microbiol. Biotechnol.* **24**, 909–915 (2008).

67. Sindhu, G. M., Murali, M., Thriveni, M. C., Anupama, N. & Amruthesh, K. N. Growth promotion and disease resistance in muskmelon induced by crude proteins of *Penicillium verruculosum* against gummy stem blight disease. *Asian J. Crop Sci.* **10**, 160–167 (2018).

68. Murali, M. & Amruthesh, K. N. Plant growth-promoting fungus *Penicillium oxalicum* enhances plant growth and induces resistance in pearl millet against downy mildew disease. *J. Phytopathol.* **163**(9), 743–754 (2015).

69. Moreno-Gavira, A., Dianez, F., Sanchez-Montesinos, B. & Santos, M. Paecilomyces variotii as a plant-growth promoter in horticulture. *Agonomy* **10**(4), 597 (2020).

70. Lu, C. et al. Paecilomyces variotii extracts (ZNC) enhance plant immunity and promote plant growth. *Plant Soil* **441**, 383–397 (2019).

71. Wang, X. et al. Paecilomyces variotii extracts and controlled-release urea synergistically increased nitrogen use efficiency and rice yield. *ACS Omega* **5**, 13303–13311 (2020).

72. Wilberforce, E. M. et al. The widespread occurrence of dark septate endophyte fungi in grassland communities. *Abstr. 7th Int. Mycol. Cong. Oslo Norway* **319**, 10 (2002).

73. Yuan, Z. L., Zhang, C. L. & Lin, F. C. Role of diverse non-systemic fungal endophytes in plant performance and response to stress: Progress and approaches. *J. Plant Growth Regul.* **29**, 116–126 (2010).

74. Mandyam, K. & Jumponen, A.., Seeking the elusive function of the root colonising dark septate endophytic fungi. *Stud. Mycol.* **53**, 173–189 (2005).

75. Fatima, F. et al. Antimicrobial and immunomodulatory efficacy of extracellularly synthesized silver and gold nanoparticles by a novel phosphate solubilizing fungus *Bipolaris tetramera*. *BMC Microbiol.* **15**, 52 (2015).

76. Priyadharsini, P. & Muthukumar, T. The root endophytic fungus *Curvularia geniculata* from *Parthenium hysterophorus* roots improves plant growth through phosphate solubilization and phytohormone production. *Fungal Ecol.* **27**, 69–77 (2017).

77. Mehta, P., Sharma, R., Putatunda, C. & Walia, A. Endophytic fungi: Role in phosphate solubilization. In *Advances in Endophytic Fungal Research* (ed. Singh, B. P.) 183–209 (Springer, 2019). https://doi.org/10.1007/978-3-030-03589-1_9

78. Etesami, H. Enhanced phosphorus fertilizer use efficiency with microorganisms. In *Nutrient Dynamics for Sustainable Crop Production* (ed. Meena, R.) 215–245 (Springer, 2020).

79. Abdalla, M. A. et al. Isolation of endophytic fungi from South African plants, and screening for their antimicrobial and extracel- lular enzymatic activities and presence of type I polyketide synthases. *S. Afr. J. Bot.* **134**, 336–342 (2020).

80. Diedhiou, S. et al. Succession of soil microbial communities during decomposition of native shrub litter of semi-arid senegal. *Pedobiologia* **52**, 273–286 (2009).

81. Badiane, A. N., Khouma, M. & Sene, M. Région de Diourbel: Gestion des sols. *Drylands Res. Work. Pap.* 15, Somerset, UK (2000). www.drylandsresearch.org.uk/pdfs/WP_Badiane_Soils.pdf. Accessed 31 May 2017.

82. FAO. Food and Agriculture Organization. World Reference Base for Soil Resources: A Framework for International Classification, Correlation and Communication World Soil Resources; Reports 103 (FAO, 2006).

83. Dossa, E. L. et al. Crop productivity and nutrient dynamics in a shrub (*Guiera senegalensis*)-based farming system of the Sahel. *Agron. J.* **104**, 1255–1264 (2013).

84. Murphy, J. & Riley, J. P. A modified single solution method for, determination of phosphate in natural waters. *Anal. Chem. Acta* **27**, 31–36 (1962).

85. Pansu, M. & Gautheyrou, J. *Handbook of Soil Analysis. Mineralogical, Organic and Inorganic Methods* 983 (Springer, 2006).

86. Gardes, M. & Bruns, T. D. ITS primers with enhanced specificity for basidiomycetes—application to identification of mycorrhizae and rusts. *Mol. Ecol.* **2**, 113–118 (1993).

87. Baldwin, M. W. Relational schemas and the processing of social information. *Psychol. Bull.* **112**, 461–484 (1992).

88. Edgar, R. C. UPARSE: Highly accurate OTU sequences from microbial amplicon reads. *Nat. Methods* **10**, 996–998 (2013).

89. Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C. & Knight, R. UCHIME improves sensitivity and speed of chimera detection. *Bioinformatics* **27**, 2194–2200 (2011).

90. Köljalg, U. et al. Towards a unified paradigm for sequence-based identification of fungi. *Mol. Ecol.* **22**, 5271–5277 (2013).

91. Oksanen, J. et al. Vegan: Community ecology package. R package version 2.4-4. https://cran.r-project.org/package=vegan (2017).

92. Weiss, S. et al. Normailzation and microbial differential abundance strategies depend upon data characteristics. *Microbiome* **5**, 1–18. https://doi.org/10.1186/s40168-017-0237-y (2017).

93. Pinheiro, J., Bates, D., DebRoy, S. & Sarkar, D. R Core Team. nlme: Linear and nonlinear mixed effects models. R package version 3.1-143. https://CRAN.R-project.org/package=nlme (2019).

94. Lenth R. Emmeans: Estimated Marginal Means, aka Least-Squares Means. https://cran.r-project.org/package=emmeans (2020).

95. Barton, K. Multi-Model Inference. Version 1.43.17 MuMln.pdf (r-project.org) (2020).

96. Anderson, M. J. A new method for non-parametric multivariate analysis of variance. *Austral. Ecol.* **26**, 32–46 (2001).

97. Anders, S. & Wolfgang, H. Differential expression analysis for sequence count data. *Genome Biol.* **11**, 106 (2010).

98. Harrell, F. E. et al. Hmisc: Harrell Miscellaneous. R package version 3.15-0 (2015). http://CRAN.R-project.org/package=Hmisc.

99. Warnes, G. R. et al. Gplots: Various R Programming Tools for Plotting Data. R Package Version 301 2. (2016).

100. Lahti, L. & Shetty, S. Microbiome analytics. Version 1.7.21. (2019).

101. Barberán, A., Bates, S. T., Casamayor, E. O. & Fierer, N. Using network analysis to explore co-occurrence patterns in soil microbial communities. *ISME J.* **6**, 343–351 (2012).

102. Wickham, H. ggplot2: Elegant Graphics for Data Analysis. Springer. http://ggplot2.tidyverse.org/. (2016).

## Acknowledgements

The authors thank Dr Lionel Moulin (IRD, Montpellier) for critical reading of the manuscript. This work was supported by the Deutscher Akademischer Austauschdienst (DAAD, PhD Grant to MTM and DT through the CERAAS), the IRD (postdoctoral Grant to MS), and the Sustainable Intensification Innovation Lab project (SIIL / Feed the Future) through the United States Agency for International Development (USAID, Grant 929773554 to DM), and the CGIAR Research Program on Grain Legumes and Dryland Cereals (support to YV and LL), and the Senegalese Ministry of Higher Education, Research and Innovation through a Grant from the FIRST program (support to AGD and AK).
**Author contributions**

A.G.D., L.L. and A.K. designed the study. M.T.M., S.P. and C.N. performed the experiments. M.T.M., A.G.D., M.S., D.T.D., D.M., Y.V., L.L. and A.K. analyzed the data. M.T.M., A.G.D., M.S., L.L. and A.K. wrote the article. All authors revised and edited the paper.

**Competing interests**

The authors declare no competing interests.

**Additional information**

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-021-04097-8.

**Correspondence** and requests for materials should be addressed to A.G.D., L.L. or A.K.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022, corrected publication 2022