- 1Département de Sciences Biologiques, Institut de Recherche en Biologie Végétale, Université de Montréal, Montréal, QC, Canada
- 2Quebec Research and Development Centre, Agriculture and Agri-Food Canada, Québec, QC, Canada
- 3Lacombe Research and Development Centre, Agriculture and Agri-Food Canada, Lacombe, AB, Canada
- 4African Genome Center, Mohammed VI Polytechnic University (UM6P), Ben Guerir, Morocco
Arbuscular mycorrhizal fungi (AMF) are obligate plant symbionts that improve the nutrition and health of their host. Most, but not all the crops form a symbiosis with AMF. It is the case for canola (Brassica napus), an important crop in the Canadian Prairies that is known to not form this association. From 2008 to 2018, an experiment was replicated at three locations of the Canadian Prairies and it was used to assess the impact of canola on the community of AMF naturally occurring in three cropping systems, canola monoculture, or canola in two different rotation systems (2-years, canola-wheat and 3-years, barley-pea-canola). We sampled canola rhizosphere and bulk soils to: (i) determine diversity and community structure of AMF, we expected that canola will negatively impact AMF communities in function of its frequency in crop rotations and (ii) wanted to assess how these AMF communities interact with other fungi and bacteria. We detected 49 AMF amplicon sequence variants (ASVs) in canola rhizosphere and bulk soils, confirming the persistence of a diversified AMF community in canola-planted soil, even after 10 years of canola monoculture, which was unexpected considering that canola is among non-mycorrhizal plants. Network analysis revealed a broad range of potential interactions between canola-associated AMF and some fungal and bacterial taxa. We report for the first time that two AMF, Funneliformis mosseae and Rhizophagus iranicus, shared their bacterial cohort almost entirely in bulk soil. Our results suggest the existence of non-species-specific AMF-bacteria or AMF-fungi relationships that could benefit AMF in absence of host plants. The persistence of an AMF community in canola rhizosphere and bulk soils brings a new light on AMF ecology and leads to new perspectives for further studies about AMF and soil microbes interactions and AMF subsistence without mycotrophic host plants.
Introduction
Arbuscular mycorrhizal fungi (AMF) are obligate plant symbionts. 400 million years ago, plants were already in association with AMF (Selosse et al., 2015). The arbuscular mycorrhizal symbiosis has evolved with land plants and, as a result, AMF associate with the roots of 90% of the plant species, including most important crops, such as cereals, legumes, and members of Solanaceae (Zhu et al., 2010). The arbuscular mycorrhizal symbiosis brings several benefits for the plant. The AMF facilitate plant nutrition in providing soil nutrients, notably P and N (Smith and Read, 2008) and can mitigate abiotic stresses, such as drought (da Silva Folli-Pereira et al., 2013; Begum et al., 2019), and protect plant roots against soil-borne pathogen attacks (St-Arnaud and Vujanovic, 2007; Lewandowski et al., 2013; Del Fabbro and Prati, 2014). The association of crop plants with AMF often increases crop yield (Hijri, 2016; Séry et al., 2016; Buczkowska and Sałata, 2020).
In bulk and rhizosphere soils, AMF play an important role in nutrient cycling and in structuring the microbial communities (Smith and Smith, 2011; Iffis et al., 2016; Dagher et al., 2020). AMF interact with the other members of the plant and soil microbiome in various ways. With their hyphae, AMF can be viewed as the bacterial highway of the soil: the labile interface between the hyphae and the soil facilitates bacterial migration (Abeysinghe et al., 2020; de Novais et al., 2020). AMF can also interact with phosphorus solubilizing bacteria (Taktek et al., 2017) or can form a tripartite interaction with other fungi to facilitate symbiosis in plant (Liu et al., 2020). But, there is still little knowledge about the interaction of AMF and other microbes in soil and plant rhizosphere.
Certain vascular plants do not host AMF in their roots and it is the case of canola (Brassica napus), an important crop in the Canadian Prairies. Canola produces glucosinolates that transform into isocyanates in soil (Smith et al., 2004; Ma et al., 2015). Among others, glucosinolates and isocyanates are toxic, which may explain the non-host plant status of canola. But, from an analysis of canola-associated soil fungi based on fungal internal transcribed spacer (ITS), Floc’h et al. (2021) reported the presence of AMF from the Glomeraceae family in the rhizosphere and bulk soil samples from fields that were in canola monoculture for 10 years. The AMF depend on a host for their carbon needs (Rich et al., 2017) and the persistence of AMF for such a long period of time in absence of mycotrophic host plants is unexpected. However, axenic AMF growth is known to be stimulated in vitro by coculture with plant growth-promoting bacteria (Abdellatif et al., 2019). Thus, the report of Floc’h et al. (2021) motivated the verification of this hypothesis: We expect that canola negatively impacts the AMF community in function of its frequency in crop rotations. Furthermore, by using network analysis, we will explore the potential interactions that AMF develop with other microbes, such as fungi and bacteria, in the canola rhizosphere and bulk soil.
In this study, we used a long-term assay held by Agriculture and Agri-Food Canada in three locations of the Canadian Prairies to identify the AMF living in canola rhizosphere and bulk soil and the samples we used are those of Floc’h et al. (2021) from which we extracted and amplified the 18S region of AMF ribosomal DNA (rDNA) to only target AMF. We used network analysis to explore the possible interactions of these AMF with the fungi and bacteria also residing in the canola rhizosphere and bulk soil.
Materials and Methods
A subset of plots from a long-term experiment initiated in 2008 and replicated at three locations in the Canadian Prairies was sampled in 2018. The experiment tested the effects of diversification in canola-based crop rotation systems. The cultivar of canola used in this study is the canola glufosinate-resistant L241C. Glufosinate amendment was 900 g ha–1. The three crop diversification treatments used in this study were: (1) monoculture of canola, (2) wheat-canola, and (3) pea-barley-canola (Table 1). These treatments were applied in a randomized complete block design with four blocks. Plot dimension was 3.7 × 15.2 m. All the rotation phases were present each year at each of the three locations where the experiment was replicated, but only the canola phases of the rotation were used in this study.
Table 1. Results of the Kruskal–Wallis tests of the effect of crop rotation diversification on indices of arbuscular mycorrhizal fungi (AMF) alpha diversity in canola rhizosphere and bulk soil (α = 0.05).
The three experiment sites were located in three pedoclimatic zones of the canola-producing regions of Canada. Two sites were in Alberta in the sub-humid brown soil zone: one site in Lacombe (latitude 52.5°N, longitude 113.7°W) and the other site in Lethbridge (latitude 49.7°N, longitude 112.8°W) and a third site was in Swift Current Saskatchewan, in the semi-arid brown soil zone (latitude 50.3°N, longitude 107.7°W). The fourth site (Melfort), used in the study by Floc’h et al. (2021), was discarded. Crops were grown according to best management practices, as described in a study by Harker et al. (2015). The growing season at all the sites was characterized by more frequent rain events in July; just before sampling, Lethbridge was drier (Supplementary Figure 1).
Rhizosphere and bulk soil samples were collected at the mid-bloom stage (50% of flowers opened) of canola development. This occurred in the fourth week of July 2018. Three to four plants within each plot were randomly selected and uprooted with a shovel. The shoots were removed and roots were placed in plastic bags and brought to the laboratory on ice in a cooler. The soil tightly attached to roots was considered as rhizosphere soil. About 5 g of rhizosphere soil per plot was collected by gently brushing the roots. The brushed roots were then gently washed with sterile distilled water. Bulk soil was taken from the top 0–7 cm soil layer by using a 2-cm diameter soil probe, exactly in between two plant rows. The top 0–7 cm of the soil is the “plow layer” in the Canadian Prairie agriculture, in contrast to many soils in humid and subhumid climates, where the plow layer is the top 0–15 or 0–20 cm. The samples were kept at 4°C before being shipped on dry ice to the laboratory in Quebec City, Quebec, where they were preserved at −80°C until DNA extraction.
More details on site description, experimental design, and sampling methods are provided in a study by Floc’h et al. (2020a).
Deoxyribonucleic Acid Extraction and Amplification
As DNA does not persist well in soil (Romanowski et al., 1993; Widmer et al., 1997; Nielsen et al., 2000), using 18S rDNA to target the AMF community in the rhizosphere and bulk soil of canola was the most adequate approach.
The bulk and rhizosphere samples DNA were extracted by using the PowerSoil™ DNA Isolation Kit (Qiagen, Montreal, Quebec, Canada) for the characterization of resident AMF. The manufacturer’s instructions for both kits were followed, except that soil and rhizosphere DNA were eluted in 50 μl. The DNA extraction of each sample was performed in duplicate and the duplicates were pooled. The quantity and quality of the DNA extracts were first verified on 1.5% agarose gel stained with Biotium GelRed® diluted at ratio 1:10,000 (VWR, Montreal, Quebec, Canada), run at 70 V for 60 min, and visualized by using the Gel Documentation System (Bio-Rad Laboratories, Mississauga, Ontario, Canada). The quantity and quality of the DNA extracts were also verified by using the Qubit Fluorometer version 2.0 (Life Technologies, Burlington, Ontario, Canada) and the Qubit double-stranded DNA (dsDNA) HS Assay Kit. DNA extracts were stored at –20°C until use.
A partial sequence of approximately 800 bp of the nuclear 18S small subunit (SSU) ribosomal RNA gene of AMF was first amplified by using the primer pair AML1/AML2 (Lee et al., 2008). The amplification was performed in 20 μl of reaction mixture in triplicate as follows: 1 μl of genomic DNA (gDNA), 200 μM of each deoxynucleoside triphosphate (dNTP), 2 mM of Mg2+, 0.8 μM of each primer, and 2.5 U of Q5 High-Fidelity DNA Polymerase NEBNext® Q5 Hot Start HiFi PCR Master Mix (BioLabs, Whitby, Ontario, Canada). The thermocycling conditions were as follows: initial denaturation at 98°C for 30 s, 20 cycles at 98°C for 10 s, 64°C for 30 s, 65°C for 60 s, and final extension performed at 65°C for 5 min. The DNA was amplified in the Biometra TProfessional Thermocycler (Biometra GmbH, Goettingen, Germany, United Kingdom). The three amplicon replicates were pooled and purified by using the QIAquick PCR Purification Kit (Qiagen, Montreal, Quebec, Canada) and eluted in 50 μl of elution buffer. This step is important to prevent interactions between the remaining primers during nested PCR. PCR products were visualized in a GelRed-stained 1.5% agarose gel.
To comply with the sequencing length capacity of the Illumina MiSeq® Reagent Kit version 3 (2 bp × 300 bp), a new primer pair from Dr. F Stefani laboratory (Ottawa Research and Development Centre, 960 Carling Avenue, Ottawa, ON K1A 06C, Canada) yielding a 490-bp length amplicon (including primers) was used to target the V3-V4 region of the nuclear 18S ribosomal RNA (rRNA) gene: nu-SSU-0450-5′ (5′-CGCAAATTACCCAATCCC-3′) and nu-SSU-0899-3′ (5′-ATAAATCCAAGAATTTCACCTC-3′). Primers were named according to the primer nomenclature system of Gargas and DePriest (1996). The number in the primer name refers to the 5′ end position of the primer on the 18S sequence standard of Saccharomyces cerevisiae (GenBank accession Z75578). Primers were designed based on the guidelines provided by the Integrated DNA Technologies (IDT Incorporation, San Diego, California, United States). Purified PCR products amplified with AML1/AML2 were used as templates for nested PCR. A 1–3 bp “heterogeneity spacer” was introduced between the 3′ end of the adapter and the 5′ end of the primer pair nu-SSU-0450-5′/nu-SSU-0899-3′ to dampen the effect of the low sequence diversity issue of the MiSeq platform. The amplification reaction mixture was the same as for the first PCR, except for the primer concentration which was 0.5 μM. The thermocycling conditions were as for the first PCR except for the number of cycles which was reduced to 15 and the annealing temperature which was 59°C. The nested PCR was performed in triplicate. Products were verified by electrophoresis on a GelRed-stained 1.5% agarose gel and replicates were pooled.
Library preparation followed the protocol described in a study by Stefani et al. (2020). Briefly, the PCR products from the nested PCR were purified by using the Agencourt AMPure® XP Beads (Beckman Coulter Incorporation, Indianapolis, Indiana, United States), normalized to 1–2 ng/μl with the SequalPrep™ Normalization Plate Kit (Thermo Fisher Scientific), and indexed by using the Nextera Index Kit (Illumina, San Diego, California, United States). Indexed amplicons were then purified and normalized. Purified indexed amplicons were quantified by quantitative PCR (qPCR) by using the LightCycler® 480 System (Roche Molecular Systems Incorporation, Branchburg, New Jersey, United States) with the KAPA Library Quantification Kit for Illumina platforms (KAPA Biosystems, Massachusetts, United States) to determine the volume of each sample to makeup a 1-nM amplicon pool for library preparation. Paired-end sequencing (2 × 300 bp) was carried out by using the Illumina MiSeq® Sequencer for 500 cycles at the Centre d’expertise et de service Génome Québec (Montreal, Quebec, Canada).
Amplicon Sequence Variant Determination and Bioinformatic Pipeline
Bioinformatics used quantitative insights inoto microbial eology 2 (QIIME2) version 2021.4 (Bolyen et al., 2019). The bioinformatic pipeline used for the processing of nu-SSU-0450-5′ and nu-SSU-0899-3′ 18S small subunit rRNA gene sequences was DADA2 version 1.18.0 (Callahan et al., 2016). First, we used Cutadapt version 3.4 to remove the primer part of the nu-SSU-0450-5′ and nu-SSU-0899-3′ RNA gene sequences with “minimum-length” at 50 and “p-error-rate” at 0.1, “—p-times” at 2 and “—p-overlap” at 6. Then, we excluded the sequences with less than 193 bp on the forward sequences and 195 bp on the reverse sequences with the command “—p-trunc-len” with “—p-mas-ee” set to 2, as the base quality of the sequences tended to diminish below that threshold in our data. Next, the Amplicon Sequence Variant (ASV) table was calculated and chimeras were eliminated. A total of 1,905 ASVs were identified by using the naïve Bayesian classifier method on the National Center for Biotechnology Information (NCBI)/nt database. Only the 62 ASVs belonging to the Mucoromycota and to the Glomeromycota or to an unknown order were kept into a phylogenetic analysis conducted to identify with certainty the ASVs that belong to the group AMF.
The phylogenetic tree was constructed from the alignment of the 62 retained ASVs, 144 sequences from Krüger et al. (2012), and 8 sequences from NCBI/nt database with multiple alignment using fast fourrier transform (MAFFT) (default settings) with the software UGENE version 39.0 (Okonechnikov et al., 2012). A maximum-likelihood tree was calculated in RAxML version 8.2.10 (Stamatakis, 2014) via cyberinfrastructure for phylogenetic research (CIPRES) (Miller et al., 2011) with bootstrap resampling set to 1,000 and the GTRGAMMA sequence evolutionary model chosen. The ASVs that were not identified as AMF in the tree were removed to produce a tree with 49 AMF ASVs.
The MiSeq sequencing data generated as part of this study are publicly available on Zenodo.1
Data Processing and Statistical Analyses
We were not able to perform PERMANOVA to assess the effect of crop diversification on AMF community structure due to data structure and neither we were able to do ANOVA to test for the effect of crop diversification on AMF alpha diversity: AMF were sometimes not found in samples leading to an empty row in our matrices. We, thus, used discriminant analysis with the software JMP version 16 (Gonzales, 2021) for community structure and replacement from indicator species analysis and the Kruskal–Wallis test for alpha diversity.
As we wanted to know with what microbes AMF could potentially interact with, we used the exact same samples from a study by Floc’h et al. (2021), in which we reused the ITS and 16S ASV tables of this study in our analysis. The protocols of DNA extraction and amplification, sequencing, and bioinformatic processing of ITS and 16S sequences are described in a study by Floc’h et al. (2021).
To assess the interactions between AMF and fungi and between AMF and bacteria, we created a co-occurrence interkingdom network by using the package SpiecEasi version 1.1.0 in R version 4.1.0 (Kurtz et al., 2015). The analysis considered rhizosphere and bulk soil fungal, bacterial, and AMF communities. The input data consisted of the raw abundance matrices of the AMF, ITS, and 16S ASVs. We first filtered the ITS and 16S datasets to remove ASVs with a frequency lower than 20% to avoid rare species. The SpiecEasi run was conducted with the algorithm “mb” with the lambda min ratio set at 10–2 and 50 repetitions. We then imported the networks into Cytoscape version 3.8.2 (Smoot et al., 2011) for plotting and used the “organic” layout to draw the networks. Edges were defined as co-occurrences or mutual exclusion based on the positive or negative values of inverse covariance linking the nodes.
Results
Raw Sequencing Datasets
The 7,513,787 reads obtained from sequencing were inputted in the pipeline yielding 4,253,351 non-chimeric reads with a mean of 50,635 reads per sample. Rarefaction curves reached saturation for all the samples (Supplementary Figure 2). These reads were assigned to 1,205 ASVs, which were classified into 49 AMF ASV after phylogenetic filtering and clustering (Supplementary Figure 3), totalizing 222,628 reads distributed heterogeneously across our 72 samples.
Taxonomic Profiles and Crop Rotation Effect on the Arbuscular Mycorrhizal Fungi Communities of Canola Rhizosphere and Bulk Soils
The taxonomic profiles of the AMF community differed in canola rhizosphere and bulk soils and varied with crop rotations (Figure 1). However, no significant effect of crop rotations on the alpha diversity of AMF in canola rhizosphere or bulk soil could be detected (Table 1).
Figure 1. Taxonomic profile of the canola arbuscular mycorrhizal fungi (AMF) community at genus level in function of the crop rotations and biotope. Three degrees of crop rotation systems were used: canola monoculture, wheat-canola rotation, and pea-barley-canola rotation.
In the rhizosphere of monocultured canola, AMF were only represented by two genera: Claroideoglomus (99.1%) and Archaeospora (0.9%). The 2-crop system (wheat-canola) showed three genera: Claroideoglomus (17.7%), Diversispora (59.7%), and Paraglomus (22.6%). Finally, the 3-crop system (pea, barley, and canola) showed six genera: Claroideoglomus (57.2%), Funneliformis (0.8%), Diversispora (2.2%), Archaeospora (28.8%), Glomus (3.1%), and one unknown AMF genus (7.7%).
In bulk soil, canola monoculture showed again only two genera of Glomeromycota, Funneliformis (65.1%) and Paraglomus (34.9%). The 2-crop system (wheat-canola) showed two genera: Funneliformis (90.39%) and Rhizophagus (9.61%) and the 3-crop system (pea, barley, and canola) showed four genera: Claroideoglomus (19.54%), Funneliformis (11.42%), Diversispora (48.64%), and Paraglomus (20.4%).
Crop rotations influenced canola AMF community structure both in the rhizosphere and bulk soil as shown by canonical analyses (Figure 2 and Supplementary Figure 4). In the rhizosphere, two AMF ASV (ASV40 Claroideoglomus spp. and ASV42 Archaeospora spp.) were associated with canola monoculture, five ASVs (ASV86 and ASV1487 belong to Diversispora, ASV358 and ASV1667 belong to Claroideoglomus, and ASV711 belongs to Paraglomus) were associated with the two-crop rotation system and six ASVs (two Claroideoglomus, two Glomus, one Paraglomus, and an unidentified AMF) were associated with the three-crop rotation system (Supplementary Figure 4).
Figure 2. Canonical analysis of AMF community from the canola rhizosphere. Circles in red, blue, and green represent the different crop rotations, while vectors show different AMF amplicon sequence variant (ASV), which are likely to be preferentially associated with each of the crop rotations. The more separated the circles are, the more the community structure between the crop rotations is different. The more the arrow is directed to the center of a circle, the more the ASV was associated with a certain crop rotation.
In bulk soil under canola monoculture, we found ASV149 (Paraglomus occultum) and ASV18 [Funneliformis mosseae (F. mosseae)], under the two-crop system, we found two Funneliformis: ASV1254 and ASV10, one Rhizophagus ASV1809, and one Glomus ASV1396, and under the three-crop system, we found three different Diversispora: ASV1774, ASV47, and ASV22, two Claroideoglomus: ASV79 and ASV203, and, finally, two Paraglomus: ASV61 and ASV1462.
Network Analysis and the Microbial Cohorts of Canola-Associated Arbuscular Mycorrhizal Fungi
Network analysis revealed that 24 AMF ASVs are putatively interacting in the rhizosphere and 26 AMF ASVs are putatively interacting in the bulk soil (Figures 3–7). Potential interactions between AMF were always co-occurrences, while no mutual exclusions were found.
Figure 3. Fungal cohorts of canola AMF in the rhizosphere. The nodes in blue are the AMF ASV and the nodes in green are the fungal ASV (FASV). Red edges signify mutual exclusion, whereas green edges signify positive co-occurrence of species.
Figure 4. Fungal cohorts of canola AMF in bulk soil. The nodes in blue are the AMF ASV and the nodes in green are the FASV. Red edges signify mutual exclusion, whereas green edges signify positive co-occurrence of species.
Figure 5. Bacterial cohort of canola AMF in the rhizosphere. The nodes in blue are the AMF ASV and the nodes in orange are the bacterial ASV (BASV). Red edges signify mutual exclusion, whereas green edges signify positive co-occurrence of species.
Figure 6. Bacterial cohort of canola AMF in bulk soil. The nodes in blue are the AMF ASV and the nodes in orange are the BASV. Red edges signify mutual exclusion, whereas green edges signify positive co-occurrence of species.
Figure 7. Bacterial cohort shared between AMF ASV58 (Funneliformis mosseae) and AMF ASV59 (Rhizophagus iranicus). The nodes in blue are the AMF ASV and the nodes in orange are the BASV. Red edges signify mutual exclusion, whereas green edges signify positive co-occurrence of species.
Microbial networks were broadly different in terms of connectivity. Fungi-AMF networks showed relatively few potential interactions: 41 nodes and 33 edges in the rhizosphere (Figure 3) and 50 nodes and 48 edges in the bulk soil (Figure 4), whereas potential interactions in the bacteria-AMF network were more abundant with 130 nodes and 142 edges in the rhizosphere (Figure 5) and 175 nodes and 227 edges in the bulk soil (Figure 6). The occurrence of AMF ASVs rarely corresponded with the occurrence of fungal ASV (FASV). In contrast, all the AMF ASV, except two (AMF ASV1722 and AMF ASV1724, both identified as Glomus indicum), have potential relationships with at least one bacterial ASV (BASV), in each of the interkingdom networks (Figures 5, 6).
In the rhizosphere, AMF cohorts rarely shared common genera with the same type of relationship, but in the bulk soil (Supplementary Table 1), the cohorts were dominated by Vicinamibacteraceae, as this clade was present 21 times across all the AMF cohorts (Supplementary Table 2). In the bulk soil, AMF ASVs rarely shared an important fraction of their bacterial or fungal cohort (Supplementary Tables 1, 2), but AMF ASV58 (F. mosseae) and AMF ASV59 [Rhizophagus iranicus (R. iranicus)] almost entirely shared their microbial cohort (Figure 7). All the BASV shared between the two AMF ASVs were shared with the same type of relationship.
Discussion
How an Arbuscular Mycorrhizal Fungi Community Persists After 10 Years of Canola Monoculture?
Arbuscular mycorrhizal fungi depend entirely on a host plant, as they cannot produce palmitic acid by themselves (Trépanier et al., 2005). It is believed that they cannot complete their life cycle, reproduce, and maintain their community in soil without a living host plant. As the experiment was conducted with a canola cultivar possessing the gene of resistance to the herbicide Liberty® and was grown with the Liberty® herbicide technology, weeds could have hardly supported AMF in this study. Weeds were for all-purpose absent, especially in canola plots. Harker et al. (2015) reported weeds 4 years from the onset of this experiment and their weed count, which was done before preseeding burnout, reflects the size of the seed bank at that time. The glufosinate is a powerful non-selective herbicide that killed 90% of the common weeds in North America at a rate of 70 g ha–1 (Hoss et al., 2003). In the present experiment, glufosinate was applied at a rate of 900 g ha–1; thus, after 10 years, an impact of weeds on the AMF community is very unlikely. Canola, as other Brassicaceae plants, is well documented to be non-mycorrhizal plants (Tommerup, 1984; Anthony et al., 2020) and it was assumed that its successive monoculture leads to a poor presence or absence of AMF in its soil. This study clearly demonstrated the persistence of AMF communities in the soil in which canola was cultivated as a sole crop for a long period of 10 years, even if the community is very poor (Supplementary Table 3). More importantly, we retrieved sequences of several AMF taxa from a narrow zone of influence (rhizosphere) of this presumably non-host crop plant. This could be explained by the following hypotheses:
(i) Arbuscular mycorrhizal fungi may colonize epidermis and external cell layers of canola roots without establishment of functional mycorrhizal symbiosis. Such superficial colonization of canola roots could be sufficient for AMF to undergo limited growth and spore production contributing to their maintenance. Mycorrhizal colonization of the superficial layer of the root cortex of Arabidopsis thaliana (A. thaliana), another non-host plant of the same family as canola, was previously reported (Veiga et al., 2013; Cosme et al., 2018). AMF can be detrimental to several non-AMF host plant genera, such as Stelaria and Pinus, by infecting their root interior (Wagg et al., 2011; Veiga et al., 2012); it is possible that AMF could infect canola in a similar way. Fonseca et al. (2001) reported that when Brassica rapa was inoculated with AMF, it reacts and increases the concentration of δ15N and δ13C in the shoot, but the influence of AMF in this case is to be relativized as the well-watered regime of the experiment emphasized the dry weight of the plants. Da̧browska et al. (2012) also reported that canola gene BnMT2 is involved with mycorrhizal symbiosis and changed its expression patterns when canola was inoculated with Glomus spores, leading to longer shoots and lower fresh biomass on 10-week-old canola seedlings that suggested a negative effect of the inoculation of AMF spores on canola. As canola roots are not a beneficial environment for AMF, their presence in the canola rhizosphere and soil could be from another reason.
(ii) Bacteria are closely or loosely interacting with AMF mycelia forming biofilms (Seneviratne et al., 2008; Lecomte et al., 2011; Iffis et al., 2014; Guennoc et al., 2017). These biofilms can help the AMF acquire phosphate through phosphate rock solubilization (Taktek et al., 2017) or allow AMF to resist harsh environmental conditions (Iffis et al., 2016). Root endophyte bacteria could be actively recruited by AMF after their penetration into the root interior and have beneficial effects on the plant host (Ujvári et al., 2021). They can also increase AMF root colonization (Bourles et al., 2020; Sangwan and Prasanna, 2021; Ujvári et al., 2021). According to Scheublin et al. (2010) and Agnolucci et al. (2015), the most common culturable bacteria ever identified in the biofilms forming on AMF spores are Pseudomonas, Streptomyces, Arthrobacter, and Oxalobacteraceae. In a study by Taktek et al. (2015), the most frequent bacteria encountered were Burkholderia; Lecomte et al. (2011) reported the dominance of Variovorax, Bacillus, Kocuria, Microbacterium, and Sphingomonas; Iffis et al. (2016) reported that Sphingomonas, Pseudomonas, Massilia, and Methylobacterium as their most abundant species associated to AMF vesicles in roots. These taxonomic profiles are very different from what we could identify in AMF cohorts by network analysis, as these genera were not detected in this study. With our more inclusive method, we found a high frequency of Vicinamibacterales in the cohorts of AMF living in canola rhizosphere and bulk soil, in the Canadian Prairie. Vicinamibacterales was associated with tolerance to trace metal contamination of soil, in particular Cu, but its ecological role still remains obscure (Chun et al., 2021). The fact that almost all the AMF in canola rhizosphere and bulk soil had a bacterial cohort in the network analysis comfort the possibility that bacteria could bring advantage to AMF in a hostile environment, facilitating interface with non-host plant. Bacteria acting as host for AMF remains another possibility (Hildebrandt et al., 2006; Horii and Ishii, 2014; Abdellatif et al., 2019).
Arbuscular mycorrhizal fungi may get some palmitate from canola through complex interactions with other microbes, such as mycorrhiza helper bacteria or Trichoderma spp. and soil fungi could also be one of the reasons. AMF were able to maintain their presence in canola rhizosphere and bulk soil through the years. However, the number of interactions of AMF with fungi recorded here was much lower than the number of interactions with bacteria (Figures 3–6). It is possible that soil fungus increases yield and overall production of host plants when co-inoculated with AMF and in presence of pathogenic fungi (Martínez-Medina et al., 2009; Nzanza et al., 2012; Commatteo et al., 2019). It can also allow AMF to colonize the rhizosphere of Brassicaceae, such as A. thaliana or canola (B. napus) (Poveda et al., 2019; Jatana et al., 2021). For example, Nectria and Leptosphaeria were found in association with AMF spores (Hijri et al., 2002). However, soil fungi ecology remains largely unknown and it is difficult to assess their real ecological roles. The fungi associated with AMF in network analysis could be AMF-helpers candidates, such as FASV592 (unknown Nectriaceae), found in the microbial cohort of AMF ASV82 (Claroideoglomus spp.) in the canola rhizosphere. The fact that AMF ASV47 (Diversispora spp.) had FASV377 Olpidium brassicae in its cohort is also of interest. O. brassicae is the most common fungus found in canola roots and its influence on canola is not yet fully understood (Lay et al., 2018a,b; Floc’h et al., 2021). The fact that O. brassicae and Diversispora spp. in canola bulk soil could share a relationship of co-occurrence is to be noted and cooperation between AMF and O. brassicae remains a possibility.
The presence of AMF in canola rhizosphere and bulk soil could also be explained by a combination of these scenarios. However, we do not rule out that dispersion of AMF propagules by wind, rain, and animals may contribute to some extent AMF to persist in canola monoculture overtime (Addy et al., 1998; Smith and Read, 2008; Cornejo et al., 2009; Pepe et al., 2018). Further investigations are needed to shed light on the mechanisms by which AMF persist in soils in the absence of non-host plants. Microbial complexity in soil and plant roots is still poorly understood and microbial interactions in plant microbiomes are just being investigated with mathematical tools since two decades (Friedman et al., 2008; Kurtz et al., 2015). It is also possible that the AMF community retrieved from canola soil and rhizosphere was a remnant of the AMF community of the previous crop from the rotations, other studies report the presence of AMF taxa in the canola rhizosphere (Poveda et al., 2019; Floc’h et al., 2021). Since we found AMF in 10-year-old canola monoculture, AMF coexistence with canola is very likely to be reality.
Effects of Rotation Systems on Arbuscular Mycorrhizal Fungi Communities and Interactions With Their Associated Microbes
Crop rotation is often used in agriculture to mitigate the accumulation of pathogens that occurs in monoculture and that is the case for canola (Hummel et al., 2009; Harker et al., 2015). Crop rotation is known to be a tool that affects the subterranean microbiota of plants (Suzuki et al., 2012; Xuan et al., 2012; Souza et al., 2013; Detheridge et al., 2016; Fan et al., 2020). However, crop rotation, depending on their diversification, impacts the community structure of canola fungal microbiota, but not its bacterial microbiota in this experiment (Floc’h et al., 2020a,b). Crop rotation is also known to have an impact on the AMF community in host plants (Higo et al., 2014, 2018; Detheridge et al., 2016; Bakhshandeh et al., 2017; Hontoria et al., 2019). In this study, we were not able to test for the effect of crop rotation on AMF community structure due to data scarcity. However, our canonical analysis show a clear differentiation between the AMF community of monoculture, 2-crop, and 3-crop systems (Figure 2 and Supplementary Figure 4). With this, we can hypothesize that AMF community composition, as the composition of fungal community in canola rhizosphere and bulk soil, is influenced by crop rotation diversification.
Arbuscular Mycorrhizal Fungi ASV58 (F. mosseae) and ASV59 (R. iranicus) Share Their Microbial Cohort in Canola Bulk Soil
Funneliformis mosseae is a widely spread AMF taxon in natural and disturbed environments. It is able to colonize a large range of hosts (Tanwar et al., 2013; Yang et al., 2017; Feitosa de Souza et al., 2018; Jie et al., 2019). It was commonly found in highly polluted environments (Hassan et al., 2011; Li et al., 2020). Rhizophagus is also a genus of AMF known for its adaptability and for being an R-strategist (Campagnac and Khasa, 2013; Sêdzielewska Toro and Brachmann, 2016; Buysens et al., 2017; Zhang et al., 2018; Cruz-Paredes et al., 2020; Lee et al., 2020). Members of the Rhizophagus genus have been commercialized and utilized as biofertilizers for more than three decades (Badri et al., 2016; Hijri, 2016; Basiru et al., 2021). Rhizophagus is also known for its ability to interact with plant growth-promoting bacteria (Battini et al., 2016; Loján et al., 2017). These two AMF ASV seem to share the same ecological niche and interact with the same bacteria. To the best of our knowledge, it is the first time that an overlap of 90% of bacterial cohort between two AMF species is reported. We know very little about the potential ecology and functions of the members shared in F. mosseae (AMF ASV58) and R. iranicus (AMF ASV59) cohorts. Among 11 taxa shared between AMF ASV58 and AMF ASV59, seven taxa were identified at genus level and among those, only four genera were the subject of more than two scientific publications: Sphingomonas, Altererythrobacter, Luteolibacter, and Gemmatimonas (Supplementary Table 2). Sphingomonas (BASV234) is the only one with a relation of mutual exclusion with the two AMF ASVs. Sphingomonas was isolated from spores of Rhizophagus irregularis and showed biofilm-like formation on hyphae (Lecomte et al., 2011). It was also found in leek rhizosphere (Nunes da Rocha et al., 2011) and diverse environments including biological soil crusts and freshwater (Ko et al., 2017; Lee et al., 2017; Zhang et al., 2017). The other three taxa were also found in diverse environments (Jiang et al., 2012; Yuan et al., 2017; Kang et al., 2019; Meng et al., 2019; Dahal et al., 2021) and Gemmatimonas was reported as a potential denitrifying bacterium (Chee-Sanford et al., 2019). Potential interactions between F. mosseae and R. iranicus and the bacterial taxa represented in their cohorts should deserve further research attention.
Conclusion
A community of AMF can be found in the rhizosphere and bulk soils of canola (B. napus), a non-host plant, even after 10 years of canola monoculture. This finding puts a new light on the ecology of AMF in the rhizosphere and bulk soil in absence of a mycotrophic host plant. AMF form a broad range of potential interactions with diverse bacterial and fungal species that may be important for AMF, especially in absence of host plants. How AMF interact with other microbes is unclear, but their associations in network analysis indicate the possibility of a direct or indirect interaction with other fungi and bacteria, which need to be clarified.
Data Availability Statement
The original contributions presented in the study are publicly available. This data can be found here: The MiSeq sequencing data generated as part of this work are publicly available on Zenodo (https://zenodo.org/record/5639078).
Author Contributions
J-BF performed the experiment, analyzed the data, and wrote the manuscript with input from all co-authors. BT conducted field experiments and collected the data. CH, MS-A, and MH supervised the work, designed the experiment, obtained fundings, and revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This study received funding from the following sources: the Natural Sciences and Engineering Research Council of Canada, CRD Fund (Grant No. CRDPJ 500507-16), the Canola Council of Canada, Saskatchewan Pulse Growers, the Agricultural Bioproducts Innovation Program (ABIP), the Prairie Canola Agronomic Research Program (PCARP), the Agriculture and Agri-Food Canada Canola Cluster Initiative, the Alberta Canola Producers Commission, the Saskatchewan Canola Development Commission, the Manitoba Canola Growers Association, and the Western Grains Research Foundation.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We are grateful to the technical staff of the Agriculture and Agrifood Canada (AAFC) research centers in Lacombe, Lethbridge, Melfort, and Swift Current for carrying out the sampling and providing useful advice.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.828145/full#supplementary-material
Footnotes
References
Abdellatif, L., Lokuruge, P., and Hamel, C. (2019). Axenic growth of the arbuscular mycorrhizal fungus Rhizophagus irregularis and growth stimulation by coculture with plant growth-promoting rhizobacteria. Mycorrhiza 29, 591–598. doi: 10.1007/s00572-019-00924-z
Abeysinghe, G., Kuchira, M., Kudo, G., Masuo, S., Ninomiya, A., Takahashi, K., et al. (2020). Fungal mycelia and bacterial thiamine establish a mutualistic growth mechanism. Life Sci. Alliance 3:e202000878. doi: 10.26508/lsa.202000878
Addy, H. D., Boswell, E. P., and Koide, R. T. (1998). Low temperature acclimation and freezing resistance of extraradical VA mycorrhizal hyphae. Mycol. Res. 102, 582–586. doi: 10.1017/S0953756297005376
Agnolucci, M., Battini, F., Cristani, C., and Giovannetti, M. (2015). Diverse bacterial communities are recruited on spores of different arbuscular mycorrhizal fungal isolates. Biol. Fertil. Soils 51, 379–389. doi: 10.1007/s00374-014-0989-5
Anthony, M. A., Celenza, J. L., Armstrong, A., and Frey, S. D. (2020). Indolic glucosinolate pathway provides resistance to mycorrhizal fungal colonization in a non-host Brassicaceae. Ecosphere 11:e03100. doi: 10.1002/ecs2.3100
Badri, A., Stefani, F. O. P., Lachance, G., Roy-Arcand, L., Beaudet, D., Vialle, A., et al. (2016). Molecular diagnostic toolkit for Rhizophagus irregularis isolate DAOM-197198 using quantitative PCR assay targeting the mitochondrial genome. Mycorrhiza 26, 721–733. doi: 10.1007/s00572-016-0708-1
Bakhshandeh, S., Corneo, P. E., Mariotte, P., Kertesz, M. A., and Dijkstra, F. A. (2017). Effect of crop rotation on mycorrhizal colonization and wheat yield under different fertilizer treatments. Agric. Ecosyst. Environ. 247, 130–136. doi: 10.1016/j.agee.2017.06.027
Basiru, S., Mwanza, H. P., and Hijri, M. (2021). Analysis of Arbuscular Mycorrhizal Fungal Inoculant Benchmarks. Microorganisms 9:81. doi: 10.3390/microorganisms9010081
Battini, F., Cristani, C., Giovannetti, M., and Agnolucci, M. (2016). Multifunctionality and diversity of culturable bacterial communities strictly associated with spores of the plant beneficial symbiont Rhizophagus intraradices. Microbiol. Res. 183, 68–79. doi: 10.1016/j.micres.2015.11.012
Begum, N., Qin, C., Ahanger, M. A., Raza, S., Khan, M. I., Ashraf, M., et al. (2019). Role of Arbuscular Mycorrhizal Fungi in Plant Growth Regulation: Implications in Abiotic Stress Tolerance. Front. Plant Sci. 10:1068. doi: 10.3389/fpls.2019.01068
Bolyen, E., Rideout, J. R., Dillon, M. R., Bokulich, N. A., Abnet, C. C., Al-Ghalith, G. A., et al. (2019). Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 37, 852–857. doi: 10.1038/s41587-019-0209-9
Bourles, A., Guentas, L., Charvis, C., Gensous, S., Majorel, C., Crossay, T., et al. (2020). Co-inoculation with a bacterium and arbuscular mycorrhizal fungi improves root colonization, plant mineral nutrition, and plant growth of a Cyperaceae plant in an ultramafic soil. Mycorrhiza 30, 121–131. doi: 10.1007/s00572-019-00929-8
Buczkowska, H. B., and Sałata, A. (2020). Inoculation with arbuscular mycorrhizal fungi (AMF) and plant irrigation with yield – forming factors in organic sweet pepper (Capsicum annuum L.) cultivation. Acta Sci. Pol. Hortorum Cultus 19, 125–138. doi: 10.24326/asphc.2020.6.11
Buysens, C., Alaux, P. L., César, V., Huret, S., Declerck, S., and Cranenbrouck, S. (2017). Tracing native and inoculated Rhizophagus irregularis in three potato cultivars (Charlotte, Nicola and Bintje) grown under field conditions. Appl. Soil Ecol. 115, 1–9. doi: 10.1016/j.apsoil.2017.03.007
Callahan, B. J., McMurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J. A., and Holmes, S. P. (2016). DADA2: High resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583. doi: 10.1038/nmeth.3869
Campagnac, E., and Khasa, D. P. (2013). Relationship between genetic variability in Rhizophagus irregularis and tolerance to saline conditions. Mycorrhiza 24, 121–129. doi: 10.1007/s00572-013-0517-8
Chee-Sanford, J., Tian, D., and Sanford, R. (2019). Consumption of N2O and other N-cycle intermediates by Gemmatimonas aurantiaca strain T-27. Microbiology 165, 1345–1354. doi: 10.1099/mic.0.000847
Chun, S. J., Kim, Y. J., Cui, Y., and Nam, K. H. (2021). Ecological network analysis reveals distinctive microbial modules associated with heavy metal contamination of abandoned mine soils in Korea. Environ. Pollut. 289:117851. doi: 10.1016/j.envpol.2021.117851
Commatteo, J. G., Consolo, V. F., Barbieri, P. A., and Covacevich, F. (2019). Indigenous arbuscular mycorrhiza and Trichoderma from systems with soybean predominance can improve tomato growth. Soil Environ. 38, 151–161. doi: 10.25252/SE/19/91805
Cornejo, P., Rubio, R., and Borie, F. (2009). Mycorrhizal Propagule Persistence in a Succession of Cereals in a Disturbed and Undisturbed Andisol Fertilized with Two Nitrogen Sources. Chil. J. Agric. Res. 69, 426–434. doi: 10.4067/S0718-58392009000300016
Cosme, M., Fernández, I., Van der Heijden, M. G. A., and Pieterse, C. M. J. (2018). Non-Mycorrhizal Plants: The Exceptions that Prove the Rule. Trends Plant Sci. 23, 577–587. doi: 10.1016/j.tplants.2018.04.004
Cruz-Paredes, C., Jakobsen, I., and Nybroe, O. (2020). Different sensitivity of a panel of Rhizophagus isolates to AMF-suppressive soils. Appl. Soil Ecol. 155:103662. doi: 10.1016/j.apsoil.2020.103662
da Silva Folli-Pereira, M., Sant’Anna Meira-Haddad, L., Vilhena da Cruz Houghton, C. M. N. S., and Megumi Kasuya, M. C. (2013). “Plant-Microorganism Interactions: Effects on the Tolerance of Plants to Biotic and Abiotic Stresses,” in Crop Improvement, eds K. R. Hakeem, P. Ahmad, and M. Ozturk (Boston: Springer US), 209–238. doi: 10.1007/978-1-4614-7028-1_6
Da̧browska, G., Hrynkiewicz, K., and Trejgell, A. (2012). Do Arbuscular Mycorrhizal Fungi Affect Metallothionein Mt2 Expression In Brassica Napus L. Roots? Acta Biol. Cracoviensia Ser. Bot. 54, 34–39. doi: 10.2478/v10182-012-0003-1
Dagher, D. J., de la Providencia, I. E., Pitre, F. E., St-Arnaud, M., and Hijri, M. (2020). Arbuscular Mycorrhizal Fungal Assemblages Significantly Shifted upon Bacterial Inoculation in Non-Contaminated and Petroleum-Contaminated Environments. Microorganisms 8:602. doi: 10.3390/microorganisms8040602
Dahal, R. H., Chaudhary, D. K., Kim, D. U., and Kim, J. (2021). Luteolibacter luteus sp. nov., isolated from stream bank soil. Arch. Microbiol. 203, 377–382. doi: 10.1007/s00203-020-02048-x
de Novais, C. B., Sbrana, C., da Conceição Jesus, E., Rouws, L. F. M., Giovannetti, M., Avio, L., et al. (2020). Mycorrhizal networks facilitate the colonization of legume roots by a symbiotic nitrogen-fixing bacterium. Mycorrhiza 30, 389–396. doi: 10.1007/s00572-020-00948-w
Del Fabbro, C., and Prati, D. (2014). Early responses of wild plant seedlings to arbuscular mycorrhizal fungi and pathogens. Basic Appl. Ecol. 15, 534–542. doi: 10.1016/j.baae.2014.08.004
Detheridge, A. P., Brand, G., Fychan, R., Crotty, F. V., Sanderson, R., Griffith, G. W., et al. (2016). The legacy effect of cover crops on soil fungal populations in a cereal rotation. Agric. Ecosyst. Environ. 228, 49–61. doi: 10.1016/j.agee.2016.04.022
Fan, P., Lai, C., Yang, J., Hong, S., Yang, Y., Wang, Q., et al. (2020). Crop rotation suppresses soil-borne Fusarium wilt of banana and alters microbial communities. Arch. Agron. Soil Sci. 0, 1–13. doi: 10.1080/03650340.2020.1839058
Feitosa de Souza, T. A., Rodriguez-Echeverria, S., Freitas, H., Alves de Andrade, L., and Santos, D. (2018). Funneliformis mosseae and Invasion by Exotic Legumes in a Brazilian Tropical Seasonal Dry Forest. Russ. J. Ecol. 49, 500–506. doi: 10.1134/S1067413618060127
Floc’h, J. B., Hamel, C., Harker, K. N., and St-Arnaud, M. (2020a). Fungal Communities of the Canola Rhizosphere: Keystone Species and Substantial Between-Year Variation of the Rhizosphere Microbiome. Microb. Ecol. 80, 762-777. doi: 10.1007/s00248-019-01475-8
Floc’h, J. B., Hamel, C., Lupwayi, N., Harker, K. N., Hijri, M., and St-Arnaud, M. (2020b). Bacterial Communities of the Canola Rhizosphere: Network Analysis Reveals a Core Bacterium Shaping Microbial Interactions. Front. Microbiol. 11:1587. doi: 10.3389/fmicb.2020.01587
Floc’h, J. B., Hamel, C., Laterrière, M., Tidemann, B., St-Arnaud, M., and Hijri, M. (2021). Inter-Kingdom Networks of Canola Microbiome Reveal Bradyrhizobium as Keystone Species and Underline the Importance of Bulk Soil in Microbial Studies to Enhance Canola Production. Microb. Ecol. [Epub online ahead of print]. doi: 10.1007/s00248-021-01905-6
Fonseca, H. M., Berbara, R. L., and Daft, M. J. (2001). Shoot δ15N and δ13C values of non-host Brassica rapa change when exposed to ±Glomus etunicatum inoculum and three levels of phosphorus and nitrogen. Mycorrhiza 11, 151–158. doi: 10.1007/s005720100125
Friedman, J., Hastie, T., and Tibshirani, R. (2008). Sparse inverse covariance estimation with the graphical lasso. Biostatistics 9, 432–441. doi: 10.1093/biostatistics/kxm045
Gargas, A., and DePriest, P. T. (1996). A nomenclature for fungal PCR primers with examples from intron-containing SSU rDNA. Mycologia 88, 745–748. doi: 10.1080/00275514.1996.12026712
Gonzales, J. E. (2021). Structural Equation Modeling with JMP§Pro. Meas. Interdiscip. Res. Perspect. 19, 80–92. doi: 10.1080/15366367.2020.1809231
Guennoc, C. M., Rose, C., Labbé, J., and Deveau, A. (2017). Bacterial biofilm formation on soil fungi: a widespread ability under controls. Microbiology 94:fiy093. doi: 10.1101/130740
Harker, K. N., O’Donovan, J. T., Turkington, T. K., Blackshaw, R. E., Lupwayi, N. Z., Smith, E. G., et al. (2015). Canola rotation frequency impacts canola yield and associated pest species. Can. J. Plant Sci. 95, 9–20. doi: 10.4141/cjps-2014-289
Hassan, S. E. D., Boon, E., St-Arnaud, M., and Hijri, M. (2011). Molecular biodiversity of arbuscular mycorrhizal fungi in trace metal-polluted soils. Mol. Ecol. 20, 3469–3483. doi: 10.1111/j.1365-294X.2011.05142.x
Higo, M., Isobe, K., Drijber, R. A., Kondo, T., Yamaguchi, M., Takeyama, S., et al. (2014). Impact of a 5-year winter cover crop rotational system on the molecular diversity of arbuscular mycorrhizal fungi colonizing roots of subsequent soybean. Biol. Fertil. Soils 50, 913–926. doi: 10.1007/s00374-014-0912-0
Higo, M., Takahashi, Y., Gunji, K., and Isobe, K. (2018). How are arbuscular mycorrhizal associations related to maize growth performance during short-term cover crop rotation? J. Sci. Food Agric. 98, 1388–1396. doi: 10.1002/jsfa.8606
Hijri, M. (2016). Analysis of a large dataset of mycorrhiza inoculation field trials on potato shows highly significant increases in yield. Mycorrhiza 26, 209–214. doi: 10.1007/s00572-015-0661-4
Hijri, M., Redecker, D., Petetot, J. A. M. C., Voigt, K., Wöstemeyer, J., and Sanders, I. R. (2002). Identification and Isolation of Two Ascomycete Fungi from Spores of the Arbuscular Mycorrhizal Fungus Scutellospora castanea. Appl. Environ. Microbiol. 68, 4567–4573. doi: 10.1128/AEM.68.9.4567-4573.2002
Hildebrandt, U., Ouziad, F., Marner, F. J., and Bothe, H. (2006). The bacterium Paenibacillus validus stimulates growth of the arbuscular mycorrhizal fungus Glomus intraradices up to the formation of fertile spores. FEMS Microbiol. Lett. 254, 258–267. doi: 10.1111/j.1574-6968.2005.00027.x
Hontoria, C., García-González, I., Quemada, M., Roldán, A., and Alguacil, M. M. (2019). The cover crop determines the AMF community composition in soil and in roots of maize after a ten-year continuous crop rotation. Sci. Total Environ. 660, 913–922. doi: 10.1016/j.scitotenv.2019.01.095
Horii, S., and Ishii, T. (2014). Effect of Arbuscular Mycorrhizal Fungi and Their Partner Bacteria on the Growth of Sesame Plants and the Concentration of Sesamin in the Seeds. Am. J. Plant Sci. 5, 3066–3072. doi: 10.4236/ajps.2014.520323
Hoss, N. E., Al-Khatib, K., Peterson, D. E., and Loughin, T. M. (2003). Efficacy of glyphosate, glufosinate, and imazethapyr on selected weed species. Weed Sci. 51, 110–117. doi: 10.1614/0043-17452003051
Hummel, J. D., Dosdall, L. M., Clayton, G. W., Turkington, T. K., Lupwayi, N. Z., Harker, K. N., et al. (2009). Canola–Wheat Intercrops for Improved Agronomic Performance and Integrated Pest Management. Agron. J. 101, 1190–1197. doi: 10.2134/agronj2009.0032
Iffis, B., St-Arnaud, M., and Hijri, M. (2014). Bacteria associated with arbuscular mycorrhizal fungi within roots of plants growing in a soil highly contaminated with aliphatic and aromatic petroleum hydrocarbons. FEMS Microbiol. Lett. 358, 44–54. doi: 10.1111/1574-6968.12533
Iffis, B., St-Arnaud, M., and Hijri, M. (2016). Petroleum hydrocarbon contamination, plant identity and arbuscular mycorrhizal fungal (AMF) community determine assemblages of the AMF spore-associated microbes: Assemblages of AMF and their spore-associated microbes. Environ. Microbiol. 18, 2689–2704. doi: 10.1111/1462-2920.13438
Jatana, B. S., Kitchens, C., Ray, C., Gerard, P., and Tharayil, N. (2021). Dual Inoculation With Arbuscular Mycorrhizal Fungi and Phosphorus Solubilizing Fungi Synergistically Enhances the Mobilization and Plant Uptake of Phosphorus From Meat and Bone Meal. Front. Soil Sci. 1:757839. doi: 10.3389/fsoil.2021.757839
Jiang, F., Li, W., Xiao, M., Dai, J., Kan, W., Chen, L., et al. (2012). Luteolibacter luojiensis sp. nov., isolated from Arctic tundra soil, and emended description of the genus Luteolibacter. Int. J. Syst. Evol. Microbiol. 62, 2259–2263. doi: 10.1099/ijs.0.037309-0
Jie, W., Lin, J., Guo, N., Cai, B., Yan, X., Jie, W., et al. (2019). Community composition of rhizosphere fungi as affected by Funneliformis mosseae in soybean continuous cropping soil during seedling period. Chil. J. Agric. Res. 79, 356–365. doi: 10.4067/S0718-58392019000300356
Kang, H., Kim, H., and Joh, K. (2019). Altererythrobacter maritimus sp. nov., isolated from seawater. Int. J. Syst. Evol. Microbiol. 69, 1566–1572. doi: 10.1099/ijsem.0.003336
Ko, Y., Hwang, W. M., Kim, M., Kang, K., and Ahn, T. Y. (2017). Sphingomonas silvisoli sp. nov., isolated from forest soil. Int. J. Syst. Evol. Microbiol. 67, 2704–2710. doi: 10.1099/ijsem.0.002001
Krüger, M., Krüger, C., Walker, C., Stockinger, H., and Schüßler, A. (2012). Phylogenetic reference data for systematics and phylotaxonomy of arbuscular mycorrhizal fungi from phylum to species level. New Phytol. 193, 970–984. doi: 10.1111/j.1469-8137.2011.03962.x
Kurtz, Z. D., Müller, C. L., Miraldi, E. R., Littman, D. R., Blaser, M. J., and Bonneau, R. A. (2015). Sparse and Compositionally Robust Inference of Microbial Ecological Networks. PLoS Comput. Biol. 11:e1004226. doi: 10.1371/journal.pcbi.1004226
Lay, C. Y., Bell, T. H., Hamel, C., Harker, K. N., Mohr, R., Greer, C. W., et al. (2018a). Canola Root–Associated Microbiomes in the Canadian Prairies. Front. Microbiol. 9:1188. doi: 10.3389/fmicb.2018.01188
Lay, C. Y., Hamel, C., and St-Arnaud, M. (2018b). Taxonomy and pathogenicity of Olpidium brassicae and its allied species. Fungal Biol. 122, 837–846. doi: 10.1016/j.funbio.2018.04.012
Lecomte, J., St-Arnaud, M., and Hijri, M. (2011). Isolation and identification of soil bacteria growing at the expense of arbuscular mycorrhizal fungi: Bacteria associated with arbuscular mycorrhizal fungi. FEMS Microbiol. Lett. 317, 43–51. doi: 10.1111/j.1574-6968.2011.02209.x
Lee, J., Lee, S., and Young, J. P. W. (2008). Improved PCR primers for the detection and identification of arbuscular mycorrhizal fungi. FEMS Microbiol. Ecol. 65, 339–349. doi: 10.1111/j.1574-6941.2008.00531.x
Lee, J. H., Kim, D. I., Choe, H. N., Lee, S. D., and Seong, C. N. (2017). Sphingomonas limnosediminicola sp. nov. and Sphingomonas palustris sp. nov., isolated from freshwater environments. Int. J. Syst. Evol. Microbiol. 67, 2834–2841. doi: 10.1099/ijsem.0.002029
Lee, S. J., Kong, M., St-Arnaud, M., and Hijri, M. (2020). Arbuscular Mycorrhizal Fungal Communities of Native Plant Species under High Petroleum Hydrocarbon Contamination Highlights Rhizophagus as a Key Tolerant Genus. Microorganisms 8:872. doi: 10.3390/microorganisms8060872
Lewandowski, T. J., Dunfield, K. E., and Antunes, P. M. (2013). Isolate Identity Determines Plant Tolerance to Pathogen Attack in Assembled Mycorrhizal Communities. PLoS One 8:e61329. doi: 10.1371/journal.pone.0061329
Li, Y., Zeng, J., Wang, S., Lin, Q., Ruan, D., Chi, H., et al. (2020). Effects of cadmium-resistant plant growth-promoting rhizobacteria and Funneliformis mosseae on the cadmium tolerance of tomato (Lycopersicon esculentum L.). Int. J. Phytoremediation 22, 451–458. doi: 10.1080/15226514.2019.1671796
Liu, H., Wu, M., Liu, J., Qu, Y., Gao, Y., and Ren, A. (2020). Tripartite Interactions Between Endophytic Fungi, Arbuscular Mycorrhizal Fungi, and Leymus chinensis. Microb. Ecol. 79, 98–109. doi: 10.1007/s00248-019-01394-8
Loján, P., Demortier, M., Velivelli, S. L. S., Pfeiffer, S., Suárez, J. P., de Vos, P., et al. (2017). Impact of plant growth-promoting rhizobacteria on root colonization potential and life cycle of Rhizophagus irregularis following co-entrapment into alginate beads. J. Appl. Microbiol. 122, 429–440. doi: 10.1111/jam.13355
Ma, Y., Gentry, T., Hu, P., Pierson, E., Gu, M., and Yin, S. (2015). Impact of brassicaceous seed meals on the composition of the soil fungal community and the incidence of Fusarium wilt on chili pepper. Appl. Soil Ecol. 90, 41–48. doi: 10.1016/j.apsoil.2015.01.016
Martínez-Medina, A., Pascual, J. A., Lloret, E., and Roldán, A. (2009). Interactions between arbuscular mycorrhizal fungi and Trichoderma harzianum and their effects on Fusarium wilt in melon plants grown in seedling nurseries. J. Sci. Food Agric. 89, 1843–1850. doi: 10.1002/jsfa.3660
Meng, F. X., Li, G., Fang, C., Wu, Y. H., Cheng, H., Chen, Y., et al. (2019). Altererythrobacter aerophilus sp. nov., isolated from deep-sea water of the north-west Pacific. Int. J. Syst. Evol. Microbiol. 69, 1689–1695. doi: 10.1099/ijsem.0.003377
Miller, M. A., Pfeiffer, W., and Schwartz, T. (2011). The CIPRES Science Gateway: a Community Resource for Phylogenetic Analyses. in Proceedings of the 2011 TeraGrid Conference: Extreme Digital Discovery TG’11. New York: Association for Computing Machinery, 1–8.
Nielsen, K. M., Smalla, K., and van Elsas, J. D. (2000). Natural Transformation of Acinetobactersp. Strain BD413 with Cell Lysates of Acinetobacter sp.,Pseudomonas fluorescens, and Burkholderia cepaciain Soil Microcosms. Appl. Environ. Microbiol. 66, 206–212. doi: 10.1128/AEM.66.1.206-212.2000
Nunes da Rocha, U., van Elsas, J. D., and van Overbeek, L. S. (2011). Verrucomicrobia subdivision 1 strains display a difference in the colonization of the leek (Allium porrum) rhizosphere: Rhizosphere colonization by Verrucomicrobia strains. FEMS Microbiol. Ecol. 78, 297–305. doi: 10.1111/j.1574-6941.2011.01156.x
Nzanza, B., Marais, D., and Soundy, P. (2012). Response of tomato (Solanum lycopersicum L.) to nursery inoculation with Trichoderma harzianum and arbuscular mycorrhizal fungi under field conditions. Acta Agric. Scand. Sect. B Soil Plant Sci. 62, 209–215. doi: 10.1080/09064710.2011.598544
Okonechnikov, K., Golosova, O., and Fursov, M. Ugene team. (2012). Unipro UGENE: a unified bioinformatics toolkit. Bioinformatics 28, 1166–1167. doi: 10.1093/bioinformatics/bts091
Pepe, A., Giovannetti, M., and Sbrana, C. (2018). Lifespan and functionality of mycorrhizal fungal mycelium are uncoupled from host plant lifespan. Sci. Rep. 8:10235. doi: 10.1038/s41598-018-28354-5
Poveda, J., Hermosa, R., Monte, E., and Nicolás, C. (2019). Trichoderma harzianum favours the access of arbuscular mycorrhizal fungi to non-host Brassicaceae roots and increases plant productivity. Sci. Rep. 9:11650. doi: 10.1038/s41598-019-48269-z
Rich, M. K., Nouri, E., Courty, P. E., and Reinhardt, D. (2017). Diet of Arbuscular Mycorrhizal Fungi: Bread and Butter? Trends Plant Sci. 22, 652–660. doi: 10.1016/j.tplants.2017.05.008
Romanowski, G., Lorenz, M. G., and Wackernagel, W. (1993). Use of polymerase chain reaction and electroporation of Escherichia coli to monitor the persistence of extracellular plasmid DNA introduced into natural soils. Appl. Environ. Microbiol. 59, 3438–46. doi: 10.1128/aem.59.10.3438-3446.1993
Sangwan, S., and Prasanna, R. (2021). Mycorrhizae Helper Bacteria: Unlocking Their Potential as Bioenhancers of Plant–Arbuscular Mycorrhizal Fungal Associations. Microb. Ecol. [Epub online ahead of print]. doi: 10.1007/s00248-021-01831-7
Scheublin, T. R., Sanders, I. R., Keel, C., and van der Meer, J. R. (2010). Characterisation of microbial communities colonising the hyphal surfaces of arbuscular mycorrhizal fungi. ISME J. 4, 752–763. doi: 10.1038/ismej.2010.5
Sêdzielewska Toro, K., and Brachmann, A. (2016). The effector candidate repertoire of the arbuscular mycorrhizal fungus Rhizophagus clarus. BMC Genom. 17:101. doi: 10.1186/s12864-016-2422-y
Selosse, M. A., Strullu-Derrien, C., Martin, F. M., Kamoun, S., and Kenrick, P. (2015). Plants, fungi and oomycetes: a 400-million year affair that shapes the biosphere. New Phytol. 206, 501–506. doi: 10.1111/nph.13371
Seneviratne, G., Zavahir, J. S., Bandara, W. M. M. S., and Weerasekara, M. L. M. A. W. (2008). Fungal-bacterial biofilms: their development for novel biotechnological applications. World J. Microbiol. Biotechnol. 24:739. doi: 10.1007/s11274-007-9539-8
Séry, D. J. M., Kouadjo, Z. G. C., Voko, B. R. R., and Zézé, A. (2016). Selecting Native Arbuscular Mycorrhizal Fungi to Promote Cassava Growth and Increase Yield under Field Conditions. Front. Microbiol. 7:2063. doi: 10.3389/fmicb.2016.02063
Smith, B. J., Kirkegaard, J. A., and Howe, G. N. (2004). Impacts of Brassica break-crops on soil biology and yield of following wheat crops. Aust. J. Agric. Res. 55:1. doi: 10.1071/AR03104
Smith, S. E., and Smith, F. A. (2011). Roles of Arbuscular Mycorrhizas in Plant Nutrition and Growth: New Paradigms from Cellular to Ecosystem Scales. Annu. Rev. Plant Biol. 62, 227–250. doi: 10.1146/annurev-arplant-042110-103846
Smoot, M. E., Ono, K., Ruscheinski, J., Wang, P. L., and Ideker, T. (2011). Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics 27, 431–432. doi: 10.1093/bioinformatics/btq675
Souza, R. C., Cantão, M. E., Vasconcelos, A. T. R., Nogueira, M. A., and Hungria, M. (2013). Soil metagenomics reveals differences under conventional and no-tillage with crop rotation or succession. Appl. Soil Ecol. 72, 49–61. doi: 10.1016/j.apsoil.2013.05.021
Stamatakis, A. (2014). RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30, 1312–1313. doi: 10.1093/bioinformatics/btu033
St-Arnaud, M., and Vujanovic, V. (2007). “Effect of the arbuscular mycorrhizal symbiosis on plant diseases and pests,” in Mycorrhizae in Crop Production, (Eds) H. Chantal and P. Christian (Binghamton: Haworth Food & Agricultural Products Press), 67–122.
Stefani, F., Bencherif, K., Sabourin, S., Hadj-Sahraoui, A. L., Banchini, C., Séguin, S., et al. (2020). Taxonomic assignment of arbuscular mycorrhizal fungi in an 18S metagenomic dataset: a case study with saltcedar (Tamarix aphylla). Mycorrhiza 30, 243–255. doi: 10.1007/s00572-020-00946-y
Suzuki, C., Takenaka, M., Oka, N., Nagaoka, K., and Karasawa, T. (2012). A DGGE analysis shows that crop rotation systems influence the bacterial and fungal communities in soils. Soil Sci. Plant Nutr. 58, 288–296. doi: 10.1080/00380768.2012.694119
Taktek, S., St-Arnaud, M., Piché, Y., Fortin, J. A., and Antoun, H. (2017). Igneous phosphate rock solubilization by biofilm-forming mycorrhizobacteria and hyphobacteria associated with Rhizoglomus irregulare DAOM 197198. Mycorrhiza 27, 13–22. doi: 10.1007/s00572-016-0726-z
Taktek, S., Trépanier, M., Servin, P. M., St-Arnaud, M., Piché, Y., Fortin, J. A., et al. (2015). Trapping of phosphate solubilizing bacteria on hyphae of the arbuscular mycorrhizal fungus Rhizophagus irregularis DAOM 197198. Soil Biol. Biochem. 90, 1–9. doi: 10.1016/j.soilbio.2015.07.016
Tanwar, A., Aggarwal, A., Yadav, A., and Parkash, V. (2013). Screening and selection of efficient host and sugarcane bagasse as substrate for mass multiplication of Funneliformis mosseae. Biol. Agric. Hortic. 29, 107–117. doi: 10.1080/01448765.2013.771955
Tommerup, I. C. (1984). Development of Infection by a Vesicular–Arbuscular Mycorrhizal Fungus in Brassica Napus L. and Trifolium Subterraneum L. New Phytol. 98, 487–495. doi: 10.1111/j.1469-8137.1984.tb04142.x
Trépanier, M., Bécard, G., Moutoglis, P., Willemot, C., Gagné, S., Avis, T. J., et al. (2005). Dependence of Arbuscular-Mycorrhizal Fungi on Their Plant Host for Palmitic Acid Synthesis. Appl. Environ. Microbiol. 71, 5341–5347. doi: 10.1128/AEM.71.9.5341-5347.2005
Ujvári, G., Turrini, A., Avio, L., and Agnolucci, M. (2021). Possible role of arbuscular mycorrhizal fungi and associated bacteria in the recruitment of endophytic bacterial communities by plant roots. Mycorrhiza 31, 527–544. doi: 10.1007/s00572-021-01040-7
Veiga, R. S. L., Faccio, A., Genre, A., Pieterse, C. M. J., Bonfante, P., and van der Heijden, M. G. (2013). Arbuscular mycorrhizal fungi reduce growth and infect roots of the non-host plant Arabidopsis thaliana. Plant Cell Environ. 36, 1926–1937. doi: 10.1111/pce.12102
Veiga, R. S. L., Howard, K., and van der Heijden, M. G. A. (2012). No evidence for allelopathic effects of arbuscular mycorrhizal fungi on the non-host plant Stellaria media. Plant Soil 360, 319–331. doi: 10.1007/s11104-012-1256-x
Wagg, C., Antunes, P. M., and Peterson, R. L. (2011). Arbuscular mycorrhizal fungal phylogeny-related interactions with a non-host. Symbiosis 53, 41–46. doi: 10.1007/s13199-011-0107-5
Widmer, F., Seidler, R. J., Donegan, K. K., and Reed, G. L. (1997). Quantification of transgenic plant marker gene persistence in the field. Mol. Ecol. 6, 1–7. doi: 10.1046/j.1365-294X.1997.00145.x
Xuan, D. T., Guong, V. T., Rosling, A., Alström, S., Chai, B., and Högberg, N. (2012). Different crop rotation systems as drivers of change in soil bacterial community structure and yield of rice, Oryza sativa. Biol. Fertil. Soils 48, 217–225. doi: 10.1007/s00374-011-0618-5
Yang, Y., Chen, Y., Li, W., Zhu, C., and Waqas, M. (2017). Inoculation of Funneliformis mosseae to Enhance Desiccation Tolerance of Populus euphratica Seedlings in Hyper-arid Region. Int. J. Agric. Biol. 19, 983–991. doi: 10.17957/IJAB/15.0369
Yuan, C. G., Chen, X., Jiang, Z., Chen, W., Liu, L., Xian, W. D., et al. (2017). Altererythrobacter lauratis sp. nov. and Altererythrobacter palmitatis sp. nov., isolated from a Tibetan hot spring. Antonie Van Leeuwenhoek 110, 1077–1086. doi: 10.1007/s10482-017-0882-y
Zhang, K., Tang, K., Feng, F., Yuan, B., Zhang, X., and Meng, J. (2017). Sphingomonas crusticola sp. nov., isolated from biological soil crusts. Int. J. Syst. Evol. Microbiol. 67, 3033–3037. doi: 10.1099/ijsem.0.002074
Zhang, X., Ren, B., Wu, S., Sun, Y., Chen, B., and Li, R. (2018). Rhizophagus irregularis influences As and P uptake by alfafa and the neighboring non-host pepperweed growing in an As-contaminated soil. J. Environ. Sci. 67, 36–44. doi: 10.1016/j.jes.2017.07.005
Keywords: arbuscular mycorrhizal fungi, rhizosphere, canola, microbiome, crop rotation, amplicon sequencing
Citation: Floc’h J-B, Hamel C, Laterrière M, Tidemann B, St-Arnaud M and Hijri M (2022) Long-Term Persistence of Arbuscular Mycorrhizal Fungi in the Rhizosphere and Bulk Soils of Non-host Brassica napus and Their Networks of Co-occurring Microbes. Front. Plant Sci. 13:828145. doi: 10.3389/fpls.2022.828145
Received: 03 December 2021; Accepted: 25 January 2022;
Published: 25 February 2022.
Edited by:
Raffaella Balestrini, Institute for Sustainable Plant Protection (CNR), ItalyReviewed by:
Alessandra Turrini, University of Pisa, ItalyLuciano Avio, University of Pisa, Italy
Katalin Posta, Szent István University, Hungary
Copyright © 2022 Floc’h, Hamel, Laterrière, Tidemann, St-Arnaud and Hijri. 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.
*Correspondence: Mohamed Hijri, mohamed.hijri@umontreal.ca