Changes in the Diversity of Soil Arbuscular Mycorrhizal Fungi after Cultivation for Biofuel Production in a Guantanamo (Cuba) Tropical System

Abstract

The arbuscular mycorrhizal fungi (AMF) are a key, integral component of the stability, sustainability and functioning of ecosystems. In this study, we characterised the AMF biodiversity in a native vegetation soil and in a soil cultivated with Jatropha curcas or Ricinus communis, in a tropical system in Guantanamo (Cuba), in order to verify if a change of land use to biofuel plant production had any effect on the AMF communities. We also asses whether some soil properties related with the soil fertility (total N, Organic C, microbial biomass C, aggregate stability percentage, pH and electrical conductivity) were changed with the cultivation of both crop species. The AM fungal small sub-unit (SSU) rRNA genes were subjected to PCR, cloning, sequencing and phylogenetic analyses. Twenty AM fungal sequence types were identified: 19 belong to the Glomeraceae and one to the Paraglomeraceae. Two AMF sequence types related to cultured AMF species (Glo G3 for Glomus sinuosum and Glo G6 for Glomus intraradices-G. fasciculatum-G. irregulare) did not occur in the soil cultivated with J. curcas and R. communis. The soil properties (total N, Organic C and microbial biomass C) were higher in the soil cultivated with the two plant species. The diversity of the AMF community decreased in the soil of both crops, with respect to the native vegetation soil, and varied significantly depending on the crop species planted. Thus, R. communis soil showed higher AMF diversity than J. curcas soil. In conclusion, R. communis could be more suitable for the long-term conservation and sustainable management of these tropical ecosytems.

About the Authors:

Maria del Mar Alguacil

* E-mail: [email protected]

Affiliation: CSIC-Centro de Edafología y Biología Aplicada del Segura, Department of Soil and Water Conservation, Campus de Espinardo, Murcia, Spain

Emma Torrecillas

Affiliation: CSIC-Centro de Edafología y Biología Aplicada del Segura, Department of Soil and Water Conservation, Campus de Espinardo, Murcia, Spain

Guillermina Hernández

Affiliation: Instituto de Ecología y Sistemática-CITMA, Grupo de Ecología del Suelo, Capdevila Rancho Boyeros, La Habana, Cuba

Antonio Roldán

Affiliation: CSIC-Centro de Edafología y Biología Aplicada del Segura, Department of Soil and Water Conservation, Campus de Espinardo, Murcia, Spain

Introduction

The exhaustion of oil reserves and requirements for the protection of the environment, together with increasing diesel fuel prices, have made it necessary to find new sources of energy. An alternative that is currently booming worldwide is fuel obtained from renewable resources or “biofuel”. In several countries around the world, plant species such as Jatropha curcas and Ricinus communis are being used to produce biodiesel [1], [2], [3], [4]. These species, native to Africa and Central America, respectively, both belong to the Euphorbiaceae family and are widely distributed in wild or semi-cultivated areas in tropical and sub-tropical countries such as Central and South America, Africa, India and South East Asia [5]. They are small trees/shrubs, well-adapted to arid conditions and able to withstand long periods of drought. Both species are multi-purpose, since besides biodiesel production they are also used for traditional medicine, hedging fences and preventing soil erosion.

Productive plant communities require a supportive soil microbial community; this is especially important for sustainable, low-cost biofuel production on marginal lands. Within the soil microbial community, the arbuscular mycorrhizal fungi (AMF) are the most-important symbionts, being a key, integral component of plant communities in both natural and agricultural ecosystems. The AMF form symbiotic relationships with the majority of land plants, including many crops [6]. They play a vital role in plant growth and productivity since they provide nutrients, especially phosphorus, to plants and, in return, receive plant carbon assimilates [6]. Moreover, AMF protect their hosts from pathogens [7], [8] and drought [9]. Also, it has been shown that AMF improve the soil structure [10], [11], thus helping to enhance the sustainability of ecosystems. Moreover, plant diversity and productivity in ecosystems are influenced significantly by the AM fungal diversity in the soil [12]. Therefore, knowledge of the natural diversity of AMF in the root-associated soil of established biofuel-crop plantations is essential for better management, sustainability and productivity of these tropical ecosystems.

Until now very few studies have been carried out with AMF and these two crop species. It has been reported that J. curcas crops are highly colonised by AMF in northern Thailand [13] and that inoculation of J. curcas seedlings with AMF can promote their establishment under NaCl-induced stress [14]. In other study a positive response to mycorrhizal inoculation in the rhizosphere of J. curcas has been reported [15]. However, in spite of their relevance to the sustainability of the ecosystems in a global change basis, the effects of biofuel plant production on the populations and diversity of AMF remain unknown. Therefore, the objective of this study was to assess, through a molecular approach, the changes in the biodiversity of AMF mediated by the cultivation of two biofuel crops, J. curcas and R. communis, in a tropical system.

Results

Soil properties

The pH and electrical conductivity did not differ significantly between the native vegetation soil and soils planted with J. curcas or R. communis (Table 1). The soil total N, total organic C and concentration of microbial biomass C were significantly higher in the soil planted with both crop species with respect to native vegetation soil; however, the aggregate stability was significantly lower.

[Figure omitted. See PDF.]

Table 1. Biological, chemical and physical properties of native vegetation soil and Ricinus communis and Jatropa curcas soil (n = 5).

https://doi.org/10.1371/journal.pone.0034887.t001

The soil used to grow J. curcas showed a significantly-higher concentration of microbial biomass C than R. communis soil, but lower organic C.

The percentage of root colonisation for J. curcas and R. communis were average values 54 and 42%, respectively, without significant differences between the species. The roots present in the native vegetation soil showed less colonisation than those in soil planted with either crop species, but without significant differences (Table 1).

PCR amplification of AMF SSU rRNA gene sequences from soil and diversity analyses

The PCR products obtained from the DNA in soil samples of native vegetation soil and soil planted with J. curcas or R. communis were used to construct SSU rDNA libraries. A total of 480 clones from 15 libraries were screened for the presence of the insert (32 clones/sample) and 320 clones resulted positive and were subsequently sequenced. BLAST analyses revealed that 206 sequences had a high similarity (96–100% identity) to AM fungal sequences (Table 2) and 114 to sequences of non-AMF origin, principally ascomycetes and basidiomycetes (data not shown); these were excluded from analyses. No chimeric sequences were detected. After phylogenetic analyses of the 206 sequences, 20 AM fungal sequence types or phylotypes could be distinguished on the basis of bootstrap values of least 85%. Of these, 19 belonged to the Glomeraceae and one to the Paraglomeraceae (Fig. 1). Since identical sequences were detected, the clones producing the same sequence for each host plant species soil sample and native vegetation soil were represented once in the alignment for clarity. A number of AMF sequence types were related to cultured AMF species. These were Para 1 for Paraglomus laccatum-occultum-brasilianum, Glo G3 for Glomus sinuosum, Glo G6 for Glomus intraradices-G. fasciculatum-G. irregulare, Glo G15 for Glomus mosseae and Glo G16 for Glomus indicum. Fifteen other AM fungal sequence types did not cluster with any sequence from AMF culture.

[Figure omitted. See PDF.]

Figure 1. Neighbour-Joining (NJ) phylogenetic tree showing AM fungal sequences isolated from of native vegetation soil and Ricinus communis and Jatropa curcas soil and reference sequences from GeneBank.

All bootstrap values >85% are shown (1000 replicates). Numbers above branches indicate the bootstrap values of the NJ analysis; numbers below branches indicate the bootstrap values of the maximum likelihood analysis. Sequences obtained in the present study are shown in bold type. They are labelled with the database accession number, the soil from the host plant from which they were obtained (C = native vegetation soil, R = R. communis soil and J = J. curcas soil) and the clone identity number. Group identifiers (for example Glo G1) are AM fungal sequences types found in our study. Endogone pisiformis and Mortierella polycephala were used as out-groups.

https://doi.org/10.1371/journal.pone.0034887.g001

[Figure omitted. See PDF.]

Table 2. Number of clones detected for each fungal type in the native vegetation soil and Ricinus communis and Jatropa curcas soil (n = 5).

https://doi.org/10.1371/journal.pone.0034887.t002

The AMF diversity, expressed by the Shannon index, was higher in native vegetation soil (H′ = 2.42) than in J. curcas (H′ = 1.20) or R. communis (H′ = 1.85) soil. The native vegetation soil showed the highest number of AMF sequence types (15), only Glo G7, Glo G8 and Glo G13 being common to the soil of both crop species (Table 2). The R. communis soil had a richer AMF community than the J. curcas soil; of eight AMF sequence types found in the soil planted with R. communis, four (Para 1, Glo G16, Glo G18 and Glo G19) appeared exclusively here and four (Glo G2, Glo G7, Glo G8 and Glo G13) were found also in the native vegetation soil. The AMF sequence type Glo G12 was exclusive to the soil planted with J. curcas.

The sampling effort curves (Fig. 2) showed that, in spite of the high number of non-AMF sequences found in this experiment, the number of samples analysed was sufficient to detect the majority of AMF sequence types present in the soils planted with J. curcas or R. communis and in the native vegetation soil samples, since the curves showed very-stable plateaus.

[Figure omitted. See PDF.]

Figure 2. Sampling effort curves for native vegetation soil, Ricinus communis soil and Jatropa curcas soil.

The sample order was randomized by 100 replications in EstimateS, version 8.0 [56].

https://doi.org/10.1371/journal.pone.0034887.g002

The distribution in the canonical correspondence analysis (CCA) diagram also shows that the AMF communities in the three soils were different: the communities from the soils of both crop species and the native vegetation soil grouped separately (Fig. 3).

[Figure omitted. See PDF.]

Figure 3. Canonical correspondence analysis (CCA) of the different AMF sequence types and soil properties detected in the three soils studied.

The eigenvalues of the first and second axes in the two-dimensional ordination diagrams are as follows: CCA1: 0.59 and CCA2: 0.55. C: native vegetation soil, R: Ricinus communis soil, J: Jatropa curcas soil, MBC:.microbial biomass C, CRL: colonized root length, AS: aggregate stability, OC: organic C, N: total N, EC: electrical conductivity.

https://doi.org/10.1371/journal.pone.0034887.g003

The relationship between the soil parameters measured, the root length colonised and the AMF diversity

The correlation analysis shows that only microbial biomass C was correlated negatively with AMF diversity, measured as the Shannon index (P<0.05) (Table 3). The rest of parameters did not show any significant correlation.

[Figure omitted. See PDF.]

Table 3. Pearson's coefficients of correlation and significance level between the soil parameters measured, the root length colonised and the AMF diversity (n = 5).

https://doi.org/10.1371/journal.pone.0034887.t003

Discussion

This is the first study using molecular approaches of the diversity of the AMF in the root-associated soil after changing the land use to biofuel plant production. Our results clearly show a significant change in the composition of the soil AMF community when the two crop species were established. In the soils where J. curcas and R. communis grew, there were specific AMF communities. Several studies have shown host preferences of AMF in different habitats [16], [17], [18], [19], [20], [21], [22]. The AMF are totally dependent on the photosynthates provided by their host plant and therefore their survival is dependent on the host plant, so the host plant identity could have determined the AMF community composition detected in the soil of each crop species [23].

Although the mechanism involved in the differences in AMF community composition between the native vegetation soil and those of the two crop species is not actually known, it has been documented that various factors could be responsible for this selective interaction.

Several surveys [24], [25] have reported than cultivated soils have a lower AMF richness and diversity than non-cultivated soils (with native vegetation). In accordance with our study, the native vegetation soil had a higher Shannon index and number of AMF sequence types (H′ = 2.42) than J. curcas soil (H′ = 1.20) and R. communis soil (H′ = 1.85).

Other factors which could be influencing the AMF community composition are: competition for plant colonization among different species and/or strains, specificity in the interaction with the host plant and factors related with the soil conditions [26]. It is known that soil-related factors, such as pH, nutrient content, total soil C and N, humidity and temperature, may control and influence the diversity of both the plant and the fungi [27]. Also, Schreiner and Mihara [28] found that the diversity of AMF in grapevine roots was related to soil fertility. In our study, there was a significant improvement in the properties related to soil fertility (total N, aggregate stability and organic and microbial biomass C) when the two plant species were established (Table 1).

An increase in the soil fertility could have been due to the higher input of root exudates as a consequence of the establishment of crop plants. It has been documented that the soil organic C released by roots promotes a dense microbial community in the immediate environment of the root that can stimulate microbial activity and, in turn, reactivate microbial populations [29]. In fact, microbial biomass C has been used frequently as an indicator of soil microbial activity [30], [31] and is related closely to soil fertility [32]. Microbial biomass C was increased significantly in our study when the two crops were planted, the soil planted with J. curcas showing the highest value (Table 1). It is noteworthy that there was a negative correlation between microbial biomass C and AMF diversity measured as the Shannon index (P<0.05). This might be due to nutritional competition between non-AM fungi and AM fungi when the microbial populations are increased in their rhizospheres, since the DNA extracted from soil planted with the crops generated clones which provided non-AMF sequences, mainly from ascomycetes and basidiomycetes. It has been reported that these fungal groups colonise as many plant species and are as abundant as the AMF [33]. In this regard, the J. curcas soil had a higher number of non-AMF sequences and lower AMF diversity than the R. communis soil (Table 2).

With regard to AMF species distribution, several studies showed that G. intraradices is considered a cosmopolitan fungus that dominates the AMF community in different environments [18], [19], [24], [25], [26], [27], [34], [35], [36], [37]. In our study, surprisingly, Glo G6, that showed very-high similarity (99% identity) to G. intraradices, was not found in the soils planted with either crop species studied. It was found only in the native vegetation soil samples. We do not know of any study where the AMF diversity in J. curcas or R. communis has been investigated. In a study carried out with spores isolated from a soil in Thailand, whose aim was to check the mycorrhizal colonization of indigenous AMF of J. curcas, Charoenpakdee et al. [13] also did not find the presence of G. intraradices. We have no previous information about the AMF which colonise these crop species; G. intraradices may even be incompatible with them. Also, it might be, together with Glo G3 (98% identity with Glomus sinuosum), ruderal species [27].

It is interesting to note that only clones detected in R. communis soil showed high homology with Glomus indicum, a new arbuscular mycorrhizal fungus that was described recently by Blaszkowski et al. [38] in the rhizosphere of Euphorbia heterophylla growing naturally in coastal sands of Alappuzha (India) and in Lactuca sativa cultivated in Asmara (Africa). Another cultured AMF species found exclusively in R. communis soil was Para 1, that showed high similarity (97%) with Paraglomus laccatum-occultum brasilianum group. Glo G15, related to the Glomus mosseae group, was found in both native vegetation soil and soil planted with J. curcas, showing its greatest abundance in the latter case.

The rest of the AMF sequence types found in this study are related to database sequences from uncultured AMF detected in different host plant roots and ecosystems or are sequences which have been found for the first time in our survey.

In our study, the soil cultivated with R. communis had higher AMF diversity than the J. curcas soil. Since the diversity of AMF in the soil can affect the diversity and productivity of plants and, therefore, the stability, sustainability and functioning of the ecosystem [12], a major concern in biofuel plant production should be to ascertain whether the crop species affect the natural AMF diversity in the soil and thus the long-term sustainability of these agroecosystems.

In conclusion, our results show that some soil properties related to soil fertility showed higher values in J. curcas and R. communis soils, cultivated as biofuel crops in a tropical system. The diversity of the AMF community in the soils of both crops decreased with respect to the native vegetation soil and varied significantly depending on the crop species planted. Thus, R. communis soil showed higher AMF diversity than J. curcas soil, suggesting that this biofuel crop could be more suitable in long-term conservation and sustainable management of these tropical ecosytems; although more studies are needed in order to confirm this hypothesis.

Materials and Methods

Ethics Statement

No specific permits were required for the described field studies since these locations/activities belong to the “Instituto de Ecología y Sistemática”. The location is not privately-owned or protected in any way. Field studies did not involve endangered or protected species.

Study Site

This research was conducted in the experimental area called “Paraguay” located in the province of Guantánamo (eastern Cuba) (coordinates 20°03′N, and 75°08′W). The climate is tropical, with an average annual rainfall of 800 mm and a mean annual temperature of 25°C. The soil type in the experimental area is an Inceptisol [39], and the natural vegetation in the zone was grassland dominated by Jatropha gossypifolia L., Dichantium caricosum L. A. Camus, Cynodon nlenfuensis Vanderyst and Cynodon plectostachyus (K. Schum) Pilg.

Experimental design

Experimental plantations were conducted using a randomized factorial design with five replication blocks of 25×25 m. Three levels were established: two crop types: Jatropha curcas L. and Ricinus communis L., and non-planted areas (conserving the natural vegetation mentioned above) used as control. Initial seedlings spacing within plantations was 1 mand at least 40 seedlings per plant type and replication block were established. The soil was not tilled (seedlings were planted in individual holes). The native plants were removed when J. curcas and R. communis were established. The plants were manually pulled up trying not to leave roots in the soil. There was not management in the soil (fertilization, herbicide, ploughing), only weeding. The experiment started in November 2007. The weeding was performed whenever seedlings germinated appeared on the surface.

Sampling

Two years and a half after plantation, in June 2010, soil samples for each treatment and replication block were taken (15 samples in total), each sample was a composite of 4 subsamples randomly taken at a depth of 10–30 cm and with a diameter of ca. 2.5 cm. The soil samples were taken from rhizospheres of biofuel plants, exactly below the plants (inside plants shade area). The soil samples were sieved through 2-mm pores to eliminate large particles and separate roots and then divided into two subsamples: one soil subsample was stored in plastic bags at −20°C for molecular analysis and another soil subsample was stored at 2°C to carry out the rest of analysis. At the same time, five root samples (aprox. 0.5 g fresh weight) per treatment and replication block (15 in total) were also harvested. In the case of native vegetation soil, the sampling was done in the grassland rhizosphere at random. These samples contained rhizospheric soil and roots.

Mycorrhizal determinations

The percentage of mycorrhizal root colonization was estimated by visual observation of fungal colonization after clearing washed roots in 10% KOH and staining with 0.05% trypan blue in lactic acid (v/v), according to Phillips and Hayman [40]. The extent of mycorrhizal colonization was calculated according to the gridline intersect method [41]. At least 100 fields per sample were checked at 80× magnification.

Soil biological, chemical and physical analyses.

Soil microbial biomass C was determined using a fumigation-extraction method [42]. Ten grams of soil at 60% of its field capacity are fumigated in a 125-ml Erlenmeyer flask with purified CHCl₃ (Panreac) for 24 h placed in a glass desiccator. After removal of residual CHCl₃, 40 ml of 0.5 M K₂SO₄ (Panreac) solution is added and the sample is shaken for 1 h before filtration of the mixture. The K₂SO₄-extracted C was determined with an automatic carbon analyser for liquid samples (Shimadzu TOC) and microbial biomass C is calculated as the difference between fumigated and non-fumigated samples.

Soil pH and electrical conductivity were measured in a 1∶5 (w/v) aqueous solution. Total N was determined with the Kjeldahl method, consisting of titration after distillation and sample digestion and the total organic C according to Yeomans and Bremner [43].

The percentage of stable aggregates was determined by the method described by Lax et al. [44]. A 4 g aliquot of sieved (0.2–4 mm) soil was placed on a small 0.250 mm sieve and wetted by spray. After 15 min the soil was subjected to an artificial rainfall of 150 ml with energy of 270 Jm⁻². The remaining soil on the sieve was placed in a previously weighed capsule (T), dried at 105°C and weighed (P1). Then, the soil was soaked in distilled water and, after 2 h, passed through the same 0.250 mm sieve with the assistance of a small stick to break the remaining aggregates. The residue remaining on the sieve, which was made up of plant debris and sand particles, was dried at 105°C and weighed (P2). The percentage of stable aggregates with regard to the total aggregates was calculated by (P1-P2)×100/(4-P2+T).

Soil DNA extraction and PCR

All PCR experiments were run using DNA preparations consisting of pooled soil extracts. DNA extractions from 15 soil samples were carried out.

For each of the 15 soil samples, genomic DNA was extracted from 0.5 g of soil fresh weight using a FastDNA™ Spin kit for soil according to the recommendations of the manufacturer (Q-BIOgene, Heidelberg, Germany). DNA extracts were stored at −20°C.

Several dilutions of extracted DNA (1/10, 1/50, 1/100) were prepared and 2 µl were used as template. Partial small subunit (SSU) ribosomal RNA gene fragments were amplified using nested PCR with the universal eukaryotic primers NS1 and NS4 [45]. PCR was carried out in a final volume of 25 µl using the “ready to go” PCR beads (Amersham Pharmacia Biotech), 0.2 µM dNTPs and 0.5 µM of each primer (PCR conditions: 94°C for 3 min, then 30 cycles at 94°C for 30 s, 40°C for 1 min, 72°C for 1 min, followed by a final extension period at 72°C for 10 min).

Two µl of several dilutions (1/10, 1/20, 1/50 and 1/100) from the first PCR were used as template DNA in a second PCR reaction performed using the specific primers AML1 and AML2 [46]. PCR reactions were carried out in a final volume of 25 µl using the “ready to go” PCR beads (Amersham Pharmacia Biotech), 0.2 µM dNTPs and 0.5 µM of each primer (PCR conditions: 94°C for 3 min, then 30 cycles of 1 min denaturation at 94°C, 1 min primer annealing at 50°C and 1 min extension at 72°C, followed by a final extension period of 10 min at 72°C). Positive and negative controls using PCR positive samples and sterile water respectively were also included in all amplifications. All the PCR reactions were run on a Perkin Elmer Cetus DNA Thermal Cycler. Reaction yields were estimated by using a 1.2% agarose gel containing ethidium bromide.

Cloning and sequencing

The PCR products were purified using a Gel extraction Kit (Qiagen) cloned into pGEM-T Easy (Promega) and transformed into Escherichia coli (Xl1 blue). Thirty-two putative positive transformants were screened in each resulting SSU rRNA gene library, using 0.7 unit of RedTaq DNA polymerase (Sigma) and the supplied reaction buffer to a final volume of 25 µl and a re-amplification with AML1 and AML2 primers with the same cycling parameters described above. Product quality and size were checked in agarose gels as described above. All clones having inserts of the correct size in each library were sequenced.

Clones were grown in liquid culture and the plasmid extracted using the QIAprep Spin Miniprep Kit (Qiagen). The sequencing was done by Laboratory of Sistemas Genómicos (Valencia, Spain) using the universal primers SP6 and T7. Sequence editing was done using the program Sequencher version 4.1.4 (Gene Codes Corporation).

A search for chimeric sequences was performed using Chimera and Cross-Over Detection and Evaluation (Ccode) [47].

Phylogenetical analysis

Sequence similarities were determined using the Basic Local Alignment Search Tool (BLASTn) sequence similarity search tool [48] provided by GenBank. Phylogenetic analysis was carried out on the sequences obtained in this study and those corresponding to the closest matches from GenBank. Sequences were aligned with other published glomeralean sequences using the program ClustalX [49] and the alignment was adjusted manually in GeneDoc [50]. Neighbour-joining (NJ) [51] and maximum likelihood (ML) phylogenetic analyses were performed with the programs PAUP4.08b [52] and RAxML v.7.0.4 [53], respectively. Distances for the NJ tree were computed using the default parameters. For the ML analysis, a GTR-GAMMA model of evolution was used. A total of 200 independent bootstrap analyses were performed to provide nodal support. The ML bootstrap values were calculated with 1000 replicates using the same substitution model. Endogone pisiformis Link and Mortierella polycephala Coem, were used as the out-groups.

Different AMF sequence types or phylotypes, were defined as groups of closely related sequences, usually with a high level of bootstrap support in the phylogenetic analyses (higher than 85%) and sequence similarity ≥97%. The pairwise analysis within clusters was carried out using MEGA software version 4 [54].

Seventy-nine representative sequences of the clones generated in this study have been deposited at the National Centre for Biotechnology Information (NCBI) GenBank (http://www.ncbi.nlm.nih.gov) under the accession numbers FR821524–FR821602.

Diversity of AM fungal community

The Shannon (H′) index was calculated as an additional measure of diversity, as it combines two components of diversity, i.e., species richness and evenness. It is calculated from the equation H′ = −∑p_i(ln p_i), where pi is the proportion of individuals found in the _ith species (in a sample, the true value of p_i is unknown but is estimated as n_i/N, [here and throughout, n_i is the number of individuals in the _ith species]).

Statistical analysis

Aggregate stability and percentage colonization were arcsine-transformed, and the other parameters were log-transformed to compensate for variance heterogeneity before analysis of variance. The data (colonization and soil properties) were subjected to analysis of variance, and comparisons among means were made using a Least Significant Difference (LSD) test calculated at p<0.05. Correlation analysis between all the soil parameters measured, colonized root length and the AMF diversity was carried out using Pearson's rank correlation coefficients. Statistical procedures were carried out with the software package IBM SPSS Statistic 19.0 for Windows.

In order to investigate the influence of change of land use on the AM fungal community composition and soil properties, ordination analyses were conducted in CANOCO for Windows v. 4.5 [55] using the relative abundance data for each AMF sequence type. Initial detrended correspondence analysis suggested a unimodal character of the data response to the sample origin (the lengths of gradients were >4); therefore, the canonical-correspondence analysis (CCA) was used. Monte Carlo permutation were conducted using 499 random permutations.

The presence or absence of AMF phylotypes in each soil sample was used to construct the sampling effort curves (with 95% confidence intervals) using the software EstimateS 8.00 [56]. The sample order was randomized by 100 replications.

Acknowledgments

MM Alguacil was supported by the Ramon and Cajal programme (Ministerio de Educación y Ciencia, Spain).

Author Contributions

Conceived and designed the experiments: MMA ET GH AR. Performed the experiments: MMA ET GH AR. Analyzed the data: MMA ET GH AR. Contributed reagents/materials/analysis tools: MMA ET GH AR. Wrote the paper: MMA AR.

Citation: Alguacil MdM, Torrecillas E, Hernández G, Roldán A (2012) Changes in the Diversity of Soil Arbuscular Mycorrhizal Fungi after Cultivation for Biofuel Production in a Guantanamo (Cuba) Tropical System. PLoS ONE 7(4): e34887. https://doi.org/10.1371/journal.pone.0034887

References

1. Flórez JP, Nieto NA, Londoño NY, Barajas CL (2004) Proceso Experimental de Producción y Operación del Biocombustible obtenido partir del Aceite de Ricino. JP FlórezNA NietoNY LondoñoCL Barajas2004Proceso Experimental de Producción y Operación del Biocombustible obtenido partir del Aceite de Ricino.Tesis doctoral. Universidad Francisco de Paula (Santander). Tesis doctoral. Universidad Francisco de Paula (Santander).

2. Prueksakorn K, Shabbir HG, Pomthong M, Sébastien B (2006) K. PrueksakornHG ShabbirM. PomthongB. Sébastien2006Energy analysis of Jatropha plantation systems for biodiesel production in Thailand, Energy for Sustainable Development, The 2nd Joint International Conference on Sustainable Energy and Environment, Bangkok, Thailand. Energy analysis of Jatropha plantation systems for biodiesel production in Thailand, Energy for Sustainable Development, The 2nd Joint International Conference on Sustainable Energy and Environment, Bangkok, Thailand.

3. Recalde ER, Durán JM (2009) Cultivos energéticos alternativos. ER RecaldeJM Durán2009Cultivos energéticos alternativos.CIITTOL. Universidad Católica del Ecuador. CIITTOL. Universidad Católica del Ecuador.

4. David ML, Joerg AP, Alberte B (2009) Modeling the land requirements and potential productivity of sugarcane and Jatropha in Brazil and India using the LPJmL dynamic global vegetation model. Biomass Bioenerg 33: 1087–1095.ML DavidAP JoergB. Alberte2009Modeling the land requirements and potential productivity of sugarcane and Jatropha in Brazil and India using the LPJmL dynamic global vegetation model.Biomass Bioenerg3310871095

5. Cano-Asseleih LM, Plumbly RA, Hylands PJ (1989) Purification and partial characterization of the hemagglutination from seeds of Jatropha curcas. J Food Biochem 13: 1–20.LM Cano-AsseleihRA PlumblyPJ Hylands1989Purification and partial characterization of the hemagglutination from seeds of Jatropha curcas.J Food Biochem13120

6. Smith SE, Read DJ (2008) Mycorrhizal symbiosis. Cambridge, UK: Academia Press. SE SmithDJ Read2008Mycorrhizal symbiosisCambridge, UKAcademia Press

7. Pozo MJ, Azcón-Aguilar C (2007) Unravelling mycorrhiza induced resistance. Curr Opin Plant Biol 10: 393–398.MJ PozoC. Azcón-Aguilar2007Unravelling mycorrhiza induced resistance.Curr Opin Plant Biol10393398

8. Sikes BA, Cottenie C, Klironomos JN (2009) Plant and fungal identity determines pathogen protection of plant roots by arbuscular mycorrhizas. J Ecol 97: 1274–1280.BA SikesC. CottenieJN Klironomos2009Plant and fungal identity determines pathogen protection of plant roots by arbuscular mycorrhizas.J Ecol9712741280

9. Augé RM, Stodola AJW, Tims JE, Saxton AM (2001) Moisture retention properties of a mycorrhizal soil. Plant Soil 230: 87–97.RM AugéAJW StodolaJE TimsAM Saxton2001Moisture retention properties of a mycorrhizal soil.Plant Soil2308797

10. Caravaca F, Barea JM, Figueroa D, Roldán A (2002) Assessing the effectiveness of mycorrhizal inoculation and soil compost addition for reafforestation with Olea europaea subsp. sylvestris through changes in soil biological and physical parameters. Appl Soil Ecol 20: 107–118.F. CaravacaJM BareaD. FigueroaA. Roldán2002Assessing the effectiveness of mycorrhizal inoculation and soil compost addition for reafforestation with Olea europaea subsp. sylvestris through changes in soil biological and physical parameters.Appl Soil Ecol20107118

11. Wilson GWT, Rice CW, Rillig MC, Springer A, Hartnett DC (2009) Soil aggregation and carbon sequestration are tightly correlated with the abundance of arbuscular mycorrhizal fungi: results from long-term field experiments. Ecol Lett 12: 452–461.GWT WilsonCW RiceMC RilligA. SpringerDC Hartnett2009Soil aggregation and carbon sequestration are tightly correlated with the abundance of arbuscular mycorrhizal fungi: results from long-term field experiments.Ecol Lett12452461

12. Van der Heijden MGA, Klironomos JN, Ursic M, Moutoglis P, Streitwolf-Engel R, et al. (1998) Mycorrhizal fungal diversity determines plant biodiversity, ecosystem variability and productivity. Nature 396: 69–72.MGA Van der HeijdenJN KlironomosM. UrsicP. MoutoglisR. Streitwolf-Engel1998Mycorrhizal fungal diversity determines plant biodiversity, ecosystem variability and productivity.Nature3966972

13. Charoenpakdee S, Phosri C, Dell B, Lumyong S (2010) The mycorrhizal status of indigenous arbuscular mycorrhizal fungi of physic nut (Jatropha curcas) in Thailand. Mycosphere 1: 167–181.S. CharoenpakdeeC. PhosriB. DellS. Lumyong2010The mycorrhizal status of indigenous arbuscular mycorrhizal fungi of physic nut (Jatropha curcas) in Thailand.Mycosphere1167181

14. Kumar A, Sharma S, Mishra S (2010) Influence of Arbuscular Mycorrhizal (AM) Fungi and Salinity on Seedling Growth, Solute Accumulation, and Mycorrhizal Dependency of Jatropha curcas L. J Plant Growth Regul 29: 297–306.A. KumarS. SharmaS. Mishra2010Influence of Arbuscular Mycorrhizal (AM) Fungi and Salinity on Seedling Growth, Solute Accumulation, and Mycorrhizal Dependency of Jatropha curcas L.J Plant Growth Regul29297306

15. Singh AK, Jamaluddin (2011) Effect of arbuscular mycorrhizal inoculation on phosphatase activity in the rhizosphere of Jatropha curcas in calcareous mined spoil. National Academy Science Letters-India 34: 27–29.AK SinghJamaluddin2011Effect of arbuscular mycorrhizal inoculation on phosphatase activity in the rhizosphere of Jatropha curcas in calcareous mined spoil.National Academy Science Letters-India342729

16. Vandenkoornhuyse P, Husband R, Daniell TJ, Watson IJ, Duck JM, et al. (2002) Arbuscular mycorrhizal community composition associated with two plant species in a grassland ecosystem. Mol Ecol 11: 1555–1564.P. VandenkoornhuyseR. HusbandTJ DaniellIJ WatsonJM Duck2002Arbuscular mycorrhizal community composition associated with two plant species in a grassland ecosystem.Mol Ecol1115551564

17. Vandenkoornhuyse P, Ridgway K, Watson IJ, Fitter AH, Young JPW (2003) Co-existing grass species have distinctive arbuscular mycorrhizal communities. Mol Ecol 12: 3085–3095.P. VandenkoornhuyseK. RidgwayIJ WatsonAH FitterJPW Young2003Co-existing grass species have distinctive arbuscular mycorrhizal communities.Mol Ecol1230853095

18. Öpik M, Moora M, Liira J, Kõljalg U, Zobel M, Sen R (2003) Divergent arbuscular mycorrhizal communities colonize roots of Pulsatilla spp. in boreal Scots pine forest and grassland soils. New Phytol 160: 581–593.M. ÖpikM. MooraJ. LiiraU. KõljalgM. ZobelR. Sen2003Divergent arbuscular mycorrhizal communities colonize roots of Pulsatilla spp. in boreal Scots pine forest and grassland soils.New Phytol160581593

19. Sýkorovà Z, Wiemken A, Redecker D (2007) Cooccurring Gentiana verna and Gentiana acaulis and their neighboring plants in two Swiss upper montane meadows harbour distinct arbuscular mycorrhizal fungal communities. Appl Environ Microbiol 73: 5426–5434.Z. SýkorovàA. WiemkenD. Redecker2007Cooccurring Gentiana verna and Gentiana acaulis and their neighboring plants in two Swiss upper montane meadows harbour distinct arbuscular mycorrhizal fungal communities.Appl Environ Microbiol7354265434

20. Alguacil MM, Roldán A, Torres MP (2009) Assessing the diversity of AM fungi in arid gypsophilous plant communities. Environ Microbiol 11: 2649–2659.MM AlguacilA. RoldánMP Torres2009Assessing the diversity of AM fungi in arid gypsophilous plant communities.Environ Microbiol1126492659

21. Alguacil MM, Torres MP, Torrecillas E, Díaz G, Roldán A (2011) Plant type differently promotes the arbuscular mycorrhizal fungi biodiversity in their rhizospheres after revegetation of a degraded, semiarid land. Soil Biol Biochem 43: 167–173.MM AlguacilMP TorresE. TorrecillasG. DíazA. Roldán2011Plant type differently promotes the arbuscular mycorrhizal fungi biodiversity in their rhizospheres after revegetation of a degraded, semiarid land.Soil Biol Biochem43167173

22. Li L, Li T, Zhang Y, Zhao Z (2010) Molecular diversity of arbuscular mycorrhizal fungi and their distribution patterns related to host-plants and habitats in a hot and arid ecosystem, southwest China. Fems Microbiol Ecol 71: 418–427.L. LiT. LiY. ZhangZ. Zhao2010Molecular diversity of arbuscular mycorrhizal fungi and their distribution patterns related to host-plants and habitats in a hot and arid ecosystem, southwest China.Fems Microbiol Ecol71418427

23. Hempel S, Renker C, Buscot F (2007) Differences in the species composition of arbuscular mycorrhizal fungi in spore, root and soil communities in a grassland ecosystem. Environ Microbiol 9: 1930–1938.S. HempelC. RenkerF. Buscot2007Differences in the species composition of arbuscular mycorrhizal fungi in spore, root and soil communities in a grassland ecosystem.Environ Microbiol919301938

24. Öpik M, Moora M, Liira J, Zobel M (2006) Composition of root-colonizing arbuscular mycorrhizal fungal communities in different ecosystems around the globe. J Ecol 94: 778–790.M. ÖpikM. MooraJ. LiiraM. Zobel2006Composition of root-colonizing arbuscular mycorrhizal fungal communities in different ecosystems around the globe.J Ecol94778790

25. Öpik M, Moora M, Zobel M, Saks U, Wheatley R, et al. (2008) High diversity of arbuscular mycorrhizal fungi in a boreal herb-rich coniferous forest. New Phytol 179: 867–876.M. ÖpikM. MooraM. ZobelU. SaksR. Wheatley2008High diversity of arbuscular mycorrhizal fungi in a boreal herb-rich coniferous forest.New Phytol179867876

26. Cesaro P, van Tuinen D, Copetta A, Chatagnier O, Berta G, et al. (2008) Preferential colonization of Solanum tuberosum L. roots by the fungus Glomus intraradices in arable soil of a potato farming area. Appl Environ Microbiol 74: 5776–5783.P. CesaroD. van TuinenA. CopettaO. ChatagnierG. Berta2008Preferential colonization of Solanum tuberosum L. roots by the fungus Glomus intraradices in arable soil of a potato farming area.Appl Environ Microbiol7457765783

27. Husband R, Herre EA, Turner SL, Gallery R, Young JPW (2002) Molecular diversity of arbuscular mycorrhizal fungi and patterns of host association over time and space in a tropical forest. Mol Ecol 11: 2669–2678.R. HusbandEA HerreSL TurnerR. GalleryJPW Young2002Molecular diversity of arbuscular mycorrhizal fungi and patterns of host association over time and space in a tropical forest.Mol Ecol1126692678

28. Schreiner RP, Mihara KL (2009) The diversity of arbuscular mycorrhizal fungi amplified from grapevine roots (Vitis vinifera L.) in Oregon vineyards is seasonally stable and influenced by soil and vine age. Mycologia 101: 599–611.RP SchreinerKL Mihara2009The diversity of arbuscular mycorrhizal fungi amplified from grapevine roots (Vitis vinifera L.) in Oregon vineyards is seasonally stable and influenced by soil and vine age.Mycologia101599611

29. Picone C (2002) Managing mycorrhizae for sustainable agriculture in the tropics. In: Vandermeer JH, editor. Tropical agroecosystems. pp. 95–132.C. Picone2002Managing mycorrhizae for sustainable agriculture in the tropics.JH VandermeerTropical agroecosystems95132CRC Press, Boca Raton, FL. CRC Press, Boca Raton, FL.

30. De Luca TH, Keeney DR (1993) Soluble anthrone-reactive carbon in soils: effect of carbon and nitrogen amendments. Soil Sci Soc Am J 57: 1296–1300.TH De LucaDR Keeney1993Soluble anthrone-reactive carbon in soils: effect of carbon and nitrogen amendments.Soil Sci Soc Am J5712961300

31. García C, Hernández MT, Costa F (1997) Potential use of dehydrogenase activity as an index of microbial activity in degraded soils. Commun Soil Sci Plan 28: 123–134.C. GarcíaMT HernándezF. Costa1997Potential use of dehydrogenase activity as an index of microbial activity in degraded soils.Commun Soil Sci Plan28123134

32. Wardle DA, Yeates GW, Nicholson KS, Bonner KI, Watson RN (1999) Response of soil microbial biomass dynamics, activity and plant litter decomposition to agricultural intensification over a seven-year period. Soil Biol Biochem 31: 1707–1720.DA WardleGW YeatesKS NicholsonKI BonnerRN Watson1999Response of soil microbial biomass dynamics, activity and plant litter decomposition to agricultural intensification over a seven-year period.Soil Biol Biochem3117071720

33. Kovács GM, Szigetvári C (2002) Mycorrhizae and other root-associated fungal structures of the plants of sandy grassland on the great Hungarian Plain. Phyton-Ann Rei Bot A 42: 211–223.GM KovácsC. Szigetvári2002Mycorrhizae and other root-associated fungal structures of the plants of sandy grassland on the great Hungarian Plain.Phyton-Ann Rei Bot A42211223

34. Wirsel SGR (2004) Homogenous stands of a wetland grass harbour diverse consortia of arbuscular mycorrhizal fungi. FEMS Microbiol Ecol 48: 129–138.SGR Wirsel2004Homogenous stands of a wetland grass harbour diverse consortia of arbuscular mycorrhizal fungi.FEMS Microbiol Ecol48129138

35. Jumpponen A, Trowbridge J, Mandyam K, Johnson L (2005) Nitrogen enrichment causes minimal changes in arbuscular mycorrhizal colonization but shifts community composition: evidence from rDNA data. Biol Fert Soils 41: 217–224.A. JumpponenJ. TrowbridgeK. MandyamL. Johnson2005Nitrogen enrichment causes minimal changes in arbuscular mycorrhizal colonization but shifts community composition: evidence from rDNA data.Biol Fert Soils41217224

36. Hijri I, Sýkorová Z, Oehl F, Ineichen K, Mäder P, et al. (2006) Communities of arbuscular mycorrhizal fungi in arable soils are not necessarily low in diversity. Mol Ecol 15: 2277–2289.I. HijriZ. SýkorováF. OehlK. IneichenP. Mäder2006Communities of arbuscular mycorrhizal fungi in arable soils are not necessarily low in diversity.Mol Ecol1522772289

37. Lekberg Y, Koide RT, Rohr JR, Aldrich-Wolfe L, Morton JB (2007) Role of niche restrictions and dispersal in the composition of arbuscular mycorrhizal fungal communities. J Ecol 95: 95–105.Y. LekbergRT KoideJR RohrL. Aldrich-WolfeJB Morton2007Role of niche restrictions and dispersal in the composition of arbuscular mycorrhizal fungal communities.J Ecol9595105

38. Błaszkowski J, Wubet T, Harikumar VS, Ryszka P, Buscot F (2010) Glomus indicum, a new arbuscular mycorrhizal fungus. Botany 88: 132–143.J. BłaszkowskiT. WubetVS HarikumarP. RyszkaF. Buscot2010Glomus indicum, a new arbuscular mycorrhizal fungus.Botany88132143

39. FAO (1998) World Reference Base for Soil Resources. FAO1998World Reference Base for Soil Resources.FAO – ISRIC – ISSS, Rome. FAO – ISRIC – ISSS, Rome.

40. Phillips JM, Hayman DS (1970) Improved procedures for clearing roots and staining parasitic and vesicular- arbuscular mycorrhizal fungi for rapid assessment of infection. Trans Br Mycol Soc 55: 158–161.JM PhillipsDS Hayman1970Improved procedures for clearing roots and staining parasitic and vesicular- arbuscular mycorrhizal fungi for rapid assessment of infection.Trans Br Mycol Soc55158161

41. Giovannetti M, Mosse B (1980) An evaluation of techniques for measuring vesicular-arbuscular mycorrhizal infection in roots. New Phytol 84: 489–499.M. GiovannettiB. Mosse1980An evaluation of techniques for measuring vesicular-arbuscular mycorrhizal infection in roots.New Phytol84489499

42. Vance ED, Brookes PC, Jenkinson D (1987) An extraction method for measuring microbial biomass carbon. Soil Biol Biochem 19: 703–707.ED VancePC BrookesD. Jenkinson1987An extraction method for measuring microbial biomass carbon.Soil Biol Biochem19703707

43. Yeomans JC, Bremner JM (1988) A rapid and precise method for routine determination of organic carbon in soil. Commun Soil Sci Plan 19: 1467–1476.JC YeomansJM Bremner1988A rapid and precise method for routine determination of organic carbon in soil.Commun Soil Sci Plan1914671476

44. Lax A, Díaz E, Castillo V, Albaladejo J (1994) Reclamation of physical and chemical properties of a salinized soil by organic amendment. Arid Soil Res Rehabil 8: 9–17.A. LaxE. DíazV. CastilloJ. Albaladejo1994Reclamation of physical and chemical properties of a salinized soil by organic amendment.Arid Soil Res Rehabil8917

45. White TJ, Bruns T, Lee S, Taylor J (1990) Amplification and direct sequencing of fungal ribosomal RNA genes for phylogen-etics. In: Innis MA, Gelfand DH, Sninsky JJ, White TJ, editors. PCR protocols. A guide to methods and applications. San Diego: Academic. pp. 315–322.TJ WhiteT. BrunsS. LeeJ. Taylor1990Amplification and direct sequencing of fungal ribosomal RNA genes for phylogen-etics.MA InnisDH GelfandJJ SninskyTJ WhitePCR protocols. A guide to methods and applicationsSan DiegoAcademic315322

46. Lee J, Lee S, Young JPW (2008) Improved PCR primers for the detection and identification of arbuscular mycorrhizal fungi. Fems Microbiol Ecol 65: 339–349.J. LeeS. LeeJPW Young2008Improved PCR primers for the detection and identification of arbuscular mycorrhizal fungi.Fems Microbiol Ecol65339349

47. Gonzalez JM, Zimmermann J, Saiz-Jimenez C (2005) Evaluating putative chimeric sequences from PCR amplified products and other cross-over events. Bioinformatics 21: 333–337.JM GonzalezJ. ZimmermannC. Saiz-Jimenez2005Evaluating putative chimeric sequences from PCR amplified products and other cross-over events.Bioinformatics21333337

48. Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, et al. (1997) Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25: 3389–3402.SF AltschulTL MaddenAA SchäfferJ. ZhangZ. Zhang1997Gapped BLAST and PSI-BLAST: a new generation of protein database search programs.Nucleic Acids Res2533893402

49. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG (1997) The ClustalX windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 24: 4876–4882.JD ThompsonTJ GibsonF. PlewniakF. JeanmouginDG Higgins1997The ClustalX windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools.Nucleic Acids Res2448764882

50. Nicholas KB, Nicholas HB Jr (1997) GeneDoc: a tool for editing and annotating multiple sequence alignments. Distributed by the authors. KB NicholasHB Nicholas Jr1997GeneDoc: a tool for editing and annotating multiple sequence alignments.Distributed by the authors

51. Saitou N, Nei M (1987) The Neighbor-Joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol 4: 406–425.N. SaitouM. Nei1987The Neighbor-Joining method: a new method for reconstructing phylogenetic trees.Mol Biol Evol4406425

52. Swofford DL (2002) PAUP: Phylogenetic Analysis Using Parsimony (and other methods), Version 4.08b10 for Machintosh. DL Swofford2002PAUP: Phylogenetic Analysis Using Parsimony (and other methods), Version 4.08b10 for Machintosh.Sinauer Associates Inc., Sunderland, MA. Sinauer Associates Inc., Sunderland, MA.

53. Stamatakis A (2006) RaxML-VI-HPC: maximum likelihoodbased phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 22: 2688–2690.A. Stamatakis2006RaxML-VI-HPC: maximum likelihoodbased phylogenetic analyses with thousands of taxa and mixed models.Bioinformatics2226882690

54. Tamura K, Dudley J, Nei M, Kumar S (2007) MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol 24: 1596–1599.K. TamuraJ. DudleyM. NeiS. Kumar2007MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0.Mol Biol Evol2415961599

55. Ter Braak CFJ, Smilauer P (2004) CANOCO reference manual and CanoDraw for Windows user's guide: software for canonical community ordination (version 4.5). CFJ Ter BraakP. Smilauer2004CANOCO reference manual and CanoDraw for Windows user's guide: software for canonical community ordination (version 4.5).Biometris, Wageningen, The Netherlands. Biometris, Wageningen, The Netherlands.

56. Colwell R (2005) EstimateS: statistical estimation of species richness and shared species from samples, version 8.0. R. Colwell2005EstimateS: statistical estimation of species richness and shared species from samples, version 8.0.User's guide. University of Connecticut, Storrs. http://viceroy.eeb.uconn.edu/estimates. User's guide. University of Connecticut, Storrs. http://viceroy.eeb.uconn.edu/estimates.