ORIGINAL RESEARCH article
Sec. Plant Pathogen Interactions
Volume 7 - 2016 | https://doi.org/10.3389/fmicb.2016.00004
Microbial Profiling of a Suppressiveness-Induced Agricultural Soil Amended with Composted Almond Shells
- Departamento de Microbiología, Facultad de Ciencias, Instituto de Hortofruticultura Subtropical y Mediterránea “La Mayora”, Universidad de Málaga, Consejo Superior de Investigaciones Científicas, Málaga, Spain
This study focused on the microbial profile present in an agricultural soil that becomes suppressive after the application of composted almond shells (AS) as organic amendments. For this purpose, we analyzed the functions and composition of the complex communities present in an experimental orchard of 40-year-old avocado trees, many of them historically amended with composted almond shells. The role of microbes in the suppression of Rosellinia necatrix, the causative agent of avocado white root rot, was determined after heat-treatment and complementation experiments with different types of soil. Bacterial and fungal profiles obtained from natural soil samples based on the 16S rRNA gene and ITS sequencing revealed slight differences among the amended (AS) and unamended (CT) soils. When the soil was under the influence of composted almond shells as organic amendments, an increase in Proteobacteria and Ascomycota groups was observed, as well as a reduction in Acidobacteria and Mortierellales. Complementary to these findings, functional analysis by GeoChip 4.6 confirmed these subtle differences, mainly present in the relative abundance of genes involved in the carbon cycle. Interestingly, a group of specific probes included in the “soil benefit” category was present only in AS-amended soils, corresponding to specific microorganisms previously described as potential biocontrol agents, such as Pseudomonas spp., Burkholderia spp., or Actinobacteria. Considering the results of both analyses, we determined that AS-amendments to the soil led to an increase in some orders of Gammaproteobacteria, Betaproteobacteria, and Dothideomycetes, as well as a reduction in the abundance of Xylariales fungi (where R. necatrix is allocated). The combination of microbial action and substrate properties of suppressiveness are discussed.
The enhancement of soil suppressiveness using organic amendments has been widely described, especially for soil-borne diseases (Lazarovits et al., 2001; Bailey and Lazarovits, 2003; van Elsas and Postma, 2007; Bonilla et al., 2012b; Pane et al., 2013). However, this effect can be extremely variable depending on the pathosystem and the environmental conditions, and there are even some examples of the amendment application increasing disease incidence (Termorshuizen et al., 2006; Janvier et al., 2007). The soils that become suppressive soils provide an environment in which plant disease development is reduced, even in the presence of a virulent pathogen and a susceptible host (Hadar and Papadopoulou, 2012). This phenomenon could be induced as a direct result of the activity of microorganism consortia that are naturally established on soil after application of the amendment (Weller et al., 2002).
As such, understanding the diversity, composition, structure, function and interactions of microbial communities is crucial to gain insight into the basis for suppressiveness mediated by this organic amendment (Janvier et al., 2007). Approaches for studying microbial communities in the soil are complex. Thus, employing genomic approaches to understand which changes occur in soil could be a good alternative strategy to decipher the profiling of soil microbiota (Garbeva et al., 2004).
The use of genomic techniques rely on PCR amplification of the conserved and variable regions of the microbial genome, commonly 16S ribosomal RNA (rRNA) for bacteria and 18S rRNA or internal transcribed sequences (ITS) for fungi, allowing for direct sequencing of these PCR amplicons using different high-throughput next-generation sequencing methods. Each group of PCR amplicons that shares a similar or identical variable region is considered an “operational taxonomic unit” (OTU) and is assumed to be equivalent to a microbial species or genus. The analysis of OTUs provide information about the phylogenetic diversity of the soil microbial community (van Elsas et al., 2007, 2008; Hirsch et al., 2013; Koyama et al., 2014).
Moreover, complementary techniques have arisen, such as microarrays, which have considerable potential in environmental microbial ecology, providing novel insights into how environmental factors affect microbial communities in various habitats (Hazen et al., 2010; He et al., 2012; Bai et al., 2013; Zhang et al., 2013; Tu et al., 2014). The GeoChip microarray is a comprehensive functional gene array (FGA) targeting hundreds to thousands of different gene families that play important roles in various biogeochemical processes, enabling researchers to comprehensively analyse the functional diversity, composition, and structure of microbial communities in various environments. It is a powerful FGA-based technology that can be used to survey the functional diversity, composition, structure, metabolic potential/activity, and dynamics of microbial communities, and then link them with ecosystem processes and functions (Xie et al., 2011; Xue et al., 2013; Cong et al., 2015).
Our research interest is focused on the avocado (Persea americana Mill.), for which southern Spain is one of the most relevant zones in the Mediterranean area for this crop. In this part of the world, one of the most limiting soilborne diseases affecting avocado trees is white root rot, caused by the fungus R. necatrix Prill. White root rot is considered to be an emergent threat to many woody crop plants worldwide (Pliego et al., 2009, 2012).
The role of soil microorganisms in the plant protection have been broadly reported. Thus, different microbes can contribute to the biocontrol of avocado white root rot using different weaponry such as antagonism (Pseudomonas chlororaphis PCL1606 or Bacillus subtilis PCL1608; Cazorla et al., 2006, 2007), competition for niches and nutrients (Calderón et al., 2014), or induction of systemic resistance or predation (Trichoderma spp.; Ruano-Rosa and López-Herrera, 2009). These microorganisms can act as single or combined with other biocontrol agents against R. necatrix (Ruano-Rosa et al., 2014). Other studies have reported the positive effect of the application of arbuscular mycorrhizal fungi to soil and the biocontrol activity on avocado (Hass and Menge, 1990; González-Cortés et al., 2012).
During the past decades, several approaches have been implemented to achieve an integrated management of R. necatrix, including physical, chemical and biological control approaches (López-Herrera et al., 1998; López-Herrera and Zea-Bonilla, 2007; González-Sánchez et al., 2013). All of these approaches seem to be effective at the experimental level, and some of them have been proven to be effective under certain conditions. However, at the same time, traditional strategies of land management have improved, and some of these strategies could be considered useful approaches to fight against diseases in avocado management, thus increasing the weaponry available against white root rot (Bonilla et al., 2012b).
One of these approaches is the use of organic amendments or mulches, which have produced beneficial effects for plants, including increasing health and yields in avocado crops (Moore-Gordon et al., 1997; Wolstenholme et al., 1997; Hermoso et al., 2011). It has been previously shown that the application of such organic matter to avocado agricultural soil can affect soil physicochemical properties and microbial communities (Bonilla et al., 2012a; López et al., 2014). Additionally, organic amendments could play a critical role in global biochemical cycles (Bonanomi et al., 2014) and could cause different effects, such as the improvement of soil fertility and the enhancement of natural suppressiveness of the soil against several phytopathogens (Cretoiu et al., 2013). Several organic amendments have shown an obvious suppressive effect against another important avocado soil-borne phytopathogen, Phytophthora cinnamomi (Bender et al., 1992; Downer et al., 2001).
In a previous study, it was shown that different organic matter applied as a mulch to the avocado crop exhibited suppressive effects against white root rot (Bonilla et al., 2015). Composted almond shells were one type of organic matter tested. The application of composted almond shells as a mulch led to an enhancement of the bacterial composition and activities of the soil communities in relation to the observed suppressiveness (Bonilla et al., 2015).
The objective of the present study was to gain insight into the microbial profiling present in the amended soils showing suppressive ability against the avocado soil-borne phytopathogen R. necatrix. The use of different microbial approaches should uncover the microbial communities potentially involved in the suppressive phenotype.
Materials and Methods
Field of Study
Soil samples were obtained from an avocado crop field (cv. Hass avocado trees grafted onto cv. Topa-Topa seedling rootstocks) located at the Experimental Station “La Mayora” (IHSM-UMA-CSIC, Málaga, Spain) on the coast of the Malaga Province (SE Spain). This experimental field of 2.5 km2 (36°75′N, 4°04′O) contains 195 40-year-old avocado trees planted at 8 × 8 m. Selected avocado trees were grouped in pairs to facilitate their management. Sixteen pairs of trees were under ecological management (massive application of composted almond shells in 2002, 2007, and 2012), and another 16 pairs of trees were under conventional management (addition of mineral nutrients twice per year, as well as the application of herbicides and pesticides when necessary, López et al., 2014) and without any organic amendment.
Natural field soil samples allocated underneath of avocado trees unamended (CT) or amended with composted almond shells (AS) were taken to perform the different experiments. Soil samples were collected in April 2013, November 2013 and April 2014. Composite soil samples were taken from four different groups of paired trees with (AS) or without (CT) organic amendment and were randomly selected from throughout the avocado orchard. To obtain a composite soil sample, two sampling distal points at 1.5 m around the trunk base for each tree of a pair of trees under the same treatment were selected; the upper layer of compost was carefully removed, and 5–10 kg of soil samples (15 cm depth) were collected per pair of trees and merged. Samples were placed in cold storage and transported to the laboratory. Samples of each type of soil were sieved through a 20 mm mesh and immediately used for physicochemical and suppressiveness experiments. To perform DNA extractions, three soil samples (1 g each) from composite soil samples per each pair of trees were sieved again (2 mm diameter) and processed independently. The remaining unused soil samples were stored at −80°C.
Physicochemical Analysis of Soil Samples
Physicochemical analysis of both types of soil samples were performed at Laboratorio Caisur S.L. (Granada, Spain) using standardized methodologies. Four samples from each composite field soil sample (AS and CT) were analyzed independently.
To test the potential role of soil microorganisms in suppressiveness, we prepared three types of processed soils using different treatments: Field soils (raw soils), heat-treated soils, and complemented soils (Table 1). We applied a moist heat treatment to the field soil samples as previously described (Weller et al., 2002), with slight modifications. Briefly, the heat treatment consisted of heating the soil in high moisture conditions at 100°C for 20 min in an autoclave. The soil was allowed to recover at 4°C overnight. Then, we performed a second treatment step, heating the soil at 100°C for 10 min in high moisture conditions. After allowing it to cool, the soil was ready to be used (Figure 1). Complemented soils were prepared with the purpose of observing the partial recovery of the microbial characteristics of the natural soil (Weller et al., 2002). The complemented soil consisted of heat-treated soil mixed with natural raw field soil in a 9:1 (w/w) ratio (Table 1).
Figure 1. Processing scheme of the soil heat-treatment and complementation used in this study for the agricultural field soil samples. The same procedure was followed for both unamended soil and soil amended with composted almond shells. T0–3 indicates sampling points to perform bacterial and fungal plate counts.
To evaluate changes in the culturable microbiota fraction during different times of the soil sample processing, counts of cultivable colony forming units (CFUs) of bacteria and fungi per gram of soil were performed. For this, 2 g samples of soil obtained at the different key times during the process were suspended in 20 ml of sterile saline solution (0.85% NaCl) with 0.5 g of sterile gravel and mixed at 150 rpm for 30 min on an orbital shaker at room temperature. Ten-fold serial dilutions of the obtained suspensions were plated on Luria Bertani (LB) agar with 100 mg of cycloheximide per liter, to analyse the heterotrophic bacteria group, and on potato dextrose agar (PDA) with 50 mg of chlortetracycline and 1 ml of tergitol NP-10 (Sigma) per liter (Bonilla et al., 2012a).
Suppressiveness assays against white root rot caused by the virulent strain R. necatrix CH53 (López-Herrera and Zea-Bonilla, 2007) were conducted using two different susceptible pathosystems, avocado (Cazorla et al., 2006) and wheat (Triticum aestivum). The R. necatrix inoculum was produced on wheat seeds (Freeman et al., 1986). The seeds were soaked for 12 h in 250-ml Erlenmeyer flasks filled with distilled water. The flasks were autoclaved after excess water had been drained off. After sterilization, fungal disks of a 1-week-old culture of R. necatrix grown on PDA were placed aseptically in each flask. Flasks were incubated at 25°C for 2–3 weeks and were shaken every 2–3 days to avoid clustering of the seeds.
Avocado/R. necatrix Test System
Six-month-old commercial avocado plants were obtained from Brokaw nurseries (Brokaw España, S.L., Vélez-Málaga, Spain). The roots from the avocado plants were disinfected by immersion in 0.1% NaOCl for 20 min and then washed twice (20 min) with sterile distilled water. Then, avocado plants were placed into square plastic pots (10.5 × 10.5 × 10.5 cm) containing 0.64 L of the sieved CT and AS types of soils. Fungal infection with R. necatrix was performed using wheat grains (4 holes of 2 cm depth were made per pot, 3 infected wheat grains were placed per hole) as previously described (Freeman et al., 1986). Non-infected plants were used as controls. Three sets of 15 avocado plants were tested per type of soil. The plants were grown in a chamber at 25°C with 70% relative humidity and 16 h of daylight and were watered twice per week. Aerial symptoms of avocado white root rot were recorded on a scale of 0–3, and a disease index (DI) was calculated after 5 weeks using the previously described formula (Cazorla et al., 2006).
Wheat/R. necatrix Test System
Wheat seeds were disinfected by immersion in 0.05% NaOCl for 10 min, washed and then placed in darkness between pieces of moist filter paper in a growth chamber for 2–3 days at 25°C to induce germination. Then, germinated seedlings were disinfected again by immersion in 0.1% NaOCl for 20 min and washed (20 min) with sterile distilled water. Seedlings were placed into plastic seedling trays (5 cm diameter × 5.5 cm) containing 0.08 L of different types of soils and either infected with R. necatrix using wheat grains (three grains per slot) or not infected to be used as controls. Three sets of 50 wheat seedlings were tested per type of soil. The seedlings were grown in a chamber at 25°C with 70% relative humidity and 16 h of daylight and were watered twice per week. Aerial symptoms were evaluated, and the disease index percentage was calculated as previously described for the avocado/R. necatrix system (Cazorla et al., 2006). Disease index percentage was recorded after evaluation of symptoms, with values ranging between 0 (healthy plant), 1 (yellowing stem base), 2 (drying stem base), and 3 (dead plant). The number of diseased seedlings was determined 7 weeks after beginning the assay, and the disease index was calculated as previously described (Cazorla et al., 2006).
Soil DNA Extraction
Soil DNA extraction was performed using 1.0 g of soil samples and a PowerSoil® DNA Isolation Kit (MOBIO Laboratories, Inc, Carlsbad, CA, USA). DNA was extracted from three independent soil samples per pair of trees for amended and unamended soil (AS and CT) and checked for quality. To test the DNA quality we performed a DNA digestion using the restriction enzyme EcoRI (New England BioLabs®, Inc., Ipswich, MA, UK) and PCR amplification of the variable region of the bacterial 16S rDNA with the universal bacterial primers 341F and 907R as described by Muyzer et al. (2004). Digestion and PCR products were analyzed for size by agarose gel electrophoresis and ethidium bromide staining. Suitable samples were mixed and DNA quantity and quality (A260/A230 > 1.8 and A260/A280 > 1.7) were evaluated using a NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies Inc., Wilmington, DE, USA).
Three independent DNA extractions were performed per each pair of trees, and then merged to create a composite DNA sample. Three of these composite DNA extractions were independently analyzed for each type of field soil (AS and CT). DNA was stored at −20°C for further analyses.
Analysis of 16S rRNA and its Gene Sequence
Two composite DNA samples from each soil type were sent for sequencing by STAB VIDA (NGS Laboratories, Caparica, Portugal) and sent to ChunLab (Seoul, Korea) to obtain the microbial DNA sequences of the 16S rRNA gene and ITS hypervariable regions. Sequences were analyzed using QIIME software (Caporaso et al., 2010) and CLcommunity™ software (ChunLab). Sequences of a length less than 200 nt were excluded from the analysis. The data were filtered for noisy sequences, checked for the presence of chimeras, and binned into OTUs (Peiffer et al., 2013) at the 97% sequence similarity level. A representative sequence of each OTU was taxonomically classified. The relative abundance of microbial clades at different taxonomic levels was calculated as the average value from two independent analyses and was used to perform the comparative distribution analysis.
Three of the composite samples of purified test DNA (800 ng per sample) from the two different types of soils studied (AS and CT) were sent to Glomics Inc (Norman, Oklahoma) for the sequencing analysis (Tu et al., 2014). Briefly, after the hybridization steps, the arrays were washing, dried and then scanned. The images obtained were analyzed by NimbleScan software (Roche NimbleGen Inc., Madison, WI) using the gridding file containing GeoChip 4.6 probes and NimbleGen control probes to determine the intensity of each spot and to identify low quality spots, which were removed prior to statistical analysis (probe spots with coefficient of variance > 0.8 were removed). Extracted data were then loaded into the GeoChip data analysis pipeline at the Institute for Environmental Genomics (Microarray Data Manager, http://ieg.ou.edu/microarray/; Liang et al., 2010; Deng and Zhou, 2013). First, the average signal intensity of the common oligo reference standard (CORS) was calculated for each array, and the maximum average value was applied to normalize the signal intensity of samples in each array. Second, the sum of the signal intensity of the samples was calculated for each array, and the maximum sum value was applied to normalize the signal intensity of all of the spots on an array, which produced a normalized value for each spot in each array. Spots were scored as positive based on a floating signal-to-noise ratio [SNR = (signal mean–background mean)/background standard deviation] so that hyperthermophile control probes accounted for 5% of positive probes. Spots that were detected in less than two samples were also removed. Before statistical analysis, logarithmic transformation was carried out for the remaining spots, and the signals of all spots were transformed into relative abundances (the sum of the number of hybridized probes for each gene category or gene function between the number of total detected probes).
Data processing was used for further analyses. Genes that overlapped between treatments (AS and CT) were calculated by dividing the number of overlapped genes between the treatments by the number of all genes detected in both treatments. Gene function diversity was calculated using the Shannon-Weiner index (H', alpha diversity) and Simpson's index (1/D, beta diversity). We performed a detrended correspondence analysis (DCA) to measure the differences of community functional gene structure between treatments. For comparing the different gene function communities, a hierarchical clustering analysis using Bray-Curtis distances was also performed. To analyse the unique detected probes in the AS samples, we performed a Venn diagram analysis using an on-line tool (http://bioinfogp.cnb.csic.es/tools/venny/). Previously, we prepared two databases by selecting genes (probes) that hybridize exclusively in each type of soil and compared them. This website provided us with a list of 2766 AS unique detected sequences from suppressive soil, which were selected to perform specific comparative analysis.
For suppressive analysis, the data were statistically analyzed using an analysis of variance (Sokal and Rohlf, 1986), followed by Fisher's least significant difference test (P = 0.05) using SPSS 22 software (SPSS Inc., Chicago). For GeoChip 4.6 analysis, significant differences in relative abundances of the microbial gene diversity between different soils were analyzed by an unpaired Student's t-test. A significance level of P < 0.1 was adopted for all comparisons. Based on the standard error, the 95% confident interval for each response variable was obtained and the significant differences between the soils were estimated.
Characteristics of Avocado Field Soils
The soil samples were taken from the same avocado orchard but from trees under different soil management (AS-amended or unamended). Soil characteristics of the experimental avocado field revealed sandy-loam textures for the amended (AS) and unamended (CT) soils. The pH was not substantially different among these samples and ranged from 7.20 to 7.55 (nearly neutral pH). Some macro- and micro-nutrients, such as potassium, iron and manganese, were also increased in the AS-amended soil (data not shown).
White Root Rot Suppressiveness Assay
Suppressiveness assays against white root rot were performed using the avocado/R. necatrix and the wheat/R. necatrix experimental plant test systems. AS-amended and unamended avocado agricultural soils, after different experimental heat treatments an complementantions were used (Figure 1; Table 1).
Bacterial and fungal counts of AS-amended and CT soil were very similar, with values of 6.5 and 6.6 log10 bacterial cfu/g, respectively, and 5.0 and 5.1 log10 fungal cfu/g, respectively. After the heat treatment of the soil, bacterial counts decreased and stabilized, without any further changes after a second heat treatment in any type of soil (Table 2). There were no differences in the results obtained for fungal count (Table 2).
Table 2. Plate counts of total heterotrophic bacteria and fungi during the soil heat-treatment of the unamended and amended with composted almond shells.
For avocado/Rosellinia test system, the disease incidence was evaluated after 5 weeks and at the end of the assay, and the disease index (DI) was calculated (Figure 2A). In these studies, AS field soil samples displayed better suppressive ability than CT field soil samples. Plants growing in the presence of AS-amended soil samples displayed a significantly lower DI than plants cultivated in the presence of CT soil samples at the end of the experiment (Figure 2A). The disease suppressiveness activity was reduced when AS soil samples were heat-treated (ASt) but showed no changes in CTt soil. Moreover, suppressiveness was complemented by soils ASc and CTt+AS, when incorporating AS soil samples. Complemented soil ASt+CT and CTc did not have a disease-suppressive ability, with levels resembling those for the heat-treated unamended soil (Figure 2A).
Figure 2. Suppressiveness assays using the avocado/R. necatrix (A) and wheat/R. necatrix (B) test systems. AS, agricultural field soil amended with composted almond shells; ASt, AS heat-treated soil; ASc, ASt complemented with AS in 9:1 (w/w) ratio; ASt+CT, ASt complemented with CT in 9:1 (w/w) ratio; CT, Agricultural field soil under conventional management; CTt, CT heat-treated soil; CTc, CTt complemented with CT in 9:1 (w/w) ratio; CTt+AS, CTt complemented with AS in 9:1 (w/w) ratio. Data were analyzed for significance after arcsine square root transformation with analysis of variance, followed by Fisher's least significant difference test (P = 0.05). Values of bars with different letters indications denote a statistically significant difference.
For the wheat/R. necatrix plant test system, disease incidence was tested 7 weeks after inoculation when the disease index (DI) was calculated (Figure 2B). Similar to the results shown by the avocado/R. necatrix test system, the AS-amended soil exhibited better suppressive ability than CT soil. The suppressiveness phenotype was significantly lost in heat-treated soils (ASt and CTt) and was partially recovered when we used amended field soil to complement (ASc and CTt+AS). The soils complemented with unamended soil, CTt and ASt+CT, had a disease-suppressive ability similar to that of heat-treated unamended soil (Figure 2B).
Characterization of the Soil Microbial Community Based on 16S rRNA Gene and its Sequencing
DNA profiling approaches and the sequencing of 16S rRNA and the ITS variable regions of extracted and mixed DNA revealed the relative abundances of microbial clades at different taxonomic levels. However, only the most abundant OTUs were quantified with a level of precision sufficient to perform the comparative distribution analysis due to the high level of OTU richness.
In both samples, Archaea were found in a very low relative abundance (< 0.1%). Thus, the bacterial 16S rRNA gene sequences allowed us to identify 33 different representative phyla in AS soil samples and 26 phyla in CT soil samples, from which 5 and 7 phyla comprise more than 1% of the community in AS and CT, respectively (Figure 3).
Figure 3. Analysis of microbial communities present in field soil samples unamended (CT) and amended with composted almond shells (AS). Relative abundance (percentage) of different prokaryotic groups detected by 16S rRNA gene sequencing analysis of soil DNA. Analysis of microbial groups are marked at the class level (thick bars) and the phylum level (boxed thin bars). < 1%, sum of all detected groups with a relative abundance less than 1%. *Taxonomic characteristics of these groups are uncertain.
In AS soil samples, the 5 most abundant phyla (above 89% of relative abundance) were Proteobacteria (50.08%), Acidobacteria (22.64%), Bacteroidetes (8.05%), Planctomycetes (4.27%), and Actinobacteria (4.09%). In contrast, the analysis of CT soil samples revealed that the most abundant (representing above 95%) phyla were Proteobacteria (45.48%), Acidobacteria (27.39%), Bacteroidetes (8.79%), Planctomycetes (60.99%), Actinobacteria (3.19%), Nitrospirae (1.70%), and Gemmatimonadetes (1.63%).
At the class level, the AS soils presented a high abundance of uncultured bacteria from the groups of Acidobacteria (EU686603, 18.44%), Gammaproteobacteria (17.85%), Alphaproteobacteria (15.28%), and Betaproteobacteria (11.4%) (Figure 3). In CT soil samples, the class analysis resulted in a similar representation of class abundance, including uncultured bacteria EU686603 (22.99%), Alphaproteobacteria (17.7%), and Gammaproteobacteria (10.7%).
In both soil samples, the phylum Proteobacteria is the most abundant (50.08 and 45.48%). Differences in this group have been shown between the two soil samples. In general, diversity is higher in AS soil samples that exhibit a predominance of the classes Gammaproteobacteria (36%) and Alphaproteobacteria (30%) and a low percentage of Deltaproteobacteria. In CT soil samples, a clear predominance of Alphaproteobacteria can be observed (39%). Remarkably, we observed an increase in AS soil samples (almost 2x) of the orders Steroidobacter (28%) and Burkholderiales (13%) and the decrease of Rhodospirales (from 18% in CT to 8% in AS) (Figure S1A).
We observed 76 different classes in AS soil samples and 65 classes in CT soil samples. We detected 24 and 13 specific bacterial classes in AS and CT, respectively, and a slightly higher richness in AS samples (Figure S2A).
The analysis of ITS sequences to reveal the abundance of eukaryotic microbes allowed us to identify a high abundance of fungal microbes. Eukaryotic microbes different from fungi ranged from 7.97% (AS) to 9.52 (CT). Among the fungi detected, the unclassified fungi comprises 8.04% (AS) and 4.28% (CT), and those below 1% represent 2.9% in CT soil samples and 3.4% in AS soil samples.
The most abundant fungal groups (approximately 70%) that are in both soil samples are of the phyla Ascomycota and Basidiomycota and of the group Mortierellales. In AS soil samples, an increase in the relative abundance of Ascomycota can be observed (Figure 4), (35.37% in CT and 45.79% in AS), as well as a reduction in the group of Mortierellales (18.37 in CT and 9.92% in AS).
Figure 4. Analysis of microbial communities present in field soil samples unamended (CT) and amended with composted almond shells (AS). Relative abundance (percentage) of different eukaryotic groups detected by ITS region sequence analysis of soil DNA. Analysis of microbial groups are marked at the class level (thick bars) and at phylum level (boxed thin bars). < 1%, sum of all detected groups with a relative abundance less than 1%. *Taxonomic characteristics of these groups are uncertain.
The analysis of the most abundant group of microorganisms (Ascomycota) revealed that in AS soil samples an increase of the class of Dothideomycetes (from 40% in CT to 54% in AS) was observed. Additionally, a reduction of the class of Sordariomycetes (from 38% in CT to 29% in AS) was observed. Also of note in reference to fungal order in AS soil samples, a huge increase of Pleosporales (from 16% in CT to 48% in AS) was observed. Remarkably, one of the fungal order that decreased in AS soil samples was the order Xylariales (from 8% in CT to 3% in AS), where the pathogen R. necatrix is allocated (Figure S1B).
We observed 39 different classes in AS soil and 50 classes in CT soil. We detected 7 and 18 specific bacterial classes in AS and CT soil, respectively, and observed a slightly higher richness in CT samples (Figure S2B).
GeoChip Analysis in Soil Samples
The number of total genes detected by GeoChip analysis and overlapping genes between treatments were measured to understand the functional diversity and structure of the microbial communities. The number of total genes detected ranged from 27348 to 28491 and from 29311 to 33526 in AS and CT samples, respectively. An unpaired Student's t-test showed that these values were significantly different. The percentage of overlapping genes between samples ranged from 77.18% for AS (77.41, 75.25, and 78.88%) to 73.16% for CT (76.25, 65.70, and 77.52%) (Figure S3). This value fell to 65.43% when we compared overlapping genes between treatments (AS1–3 and CT1–3). DCA (detrended correspondence analysis) and hierarchical clustering (with Bray-Curtis distance) were performed (Figure S3) using all of the detected genes, showing that functional structure of the microbial community was similar in the replicates but different among the soils (AS and CT).
To understand the effects of composted almond shells on the microbial communities and the acquired suppressive capacity, microbial functional genes categorized as participating in biogeochemical cycles and other important soil processes were examined (Figure 5). Gene functions related to the carbon cycle were the gene category most represented in all samples. C cycling probes were significantly more abundant than other categories in AS samples (36.65% in AS and 34.54% in CT), whereas genes related to organic contaminant degradation (12.42% in AS and 12.81% in CT), metal resistance (14.58% in AS and 16.32 in CT) and virulence (1.59% in AS and 1.61% in CT) were significantly more abundant in CT samples. There were no significant differences in N, P, and S cycle genes and other gene categories such as stress, fungi functions, soil benefit and soil borne pathogens (Figure 5).
Figure 5. GeoChip analysis of functional gene categories. Relative abundance of all detected genes from different gene categories analyzed in this study. *Indicates significant statistical differences (p < 0.1) between the two types of soil samples, amended soil (AS) and conventionally managed soil (CT). Standard error bars are shown.
Key genes for acetogenesis, C degradation, C fixation, methane metabolism, and other genes related to the C cycle were detected in the two types of soils (Figure S4A). The relative abundance of genes related to the C degradation category were the highest and exhibited significant differences between the AS samples and the CT samples. In this category, we found the presence of degradative genes of the most abundant C sources derived from plant and animal sources that could be present in soil ecosystems, such as starch, hemicellulose, cellulose, chitin, and lignin. There were few significant differences between samples in these categories of detected genes (Figure S4A).
Of the nitrogen cycle category, only the ammonification subcategory had a higher significant difference for amended soil (Figure S4B). In this subcategory, there are genes that function in the decomposition of organic matter and cycling of accumulated N.
Related to the sulfur cycle, the analyses performed exhibited a higher significant difference (P < 0.1) in only the sulphite reductase genes of AS samples compared to CT samples. These genes encode enzymes that catalyze the reduction of sulphite to sulfide, using iron as cofactor, and provide a source of S to microbiota. The CT samples exhibited a higher significant difference in sulfate reductase, a protein involved in sulfur reduction by anaerobic respiration (Figure S4C).
Statistical analyses showed no significant differences in the relative abundance of genes involved in the phosphorous cycle for these samples.
The analysis of genes in the category of environmental adaptability showed significant differences (P < 0.1) in the subcategories, as shown in Figures S4D–F. Genes involved in the organic degradation of aromatics, such as chlorinated and pesticide-related compounds, had a higher significant relative abundance for amended soil than conventional managed soil. Similar results were obtained for genes related to osmotic and oxygen stress, from the stress category, and metal resistance to cobalt and lead, which had slightly higher significant relative abundance for AS samples than CT samples. On the other hand, unamended soils exhibited significantly higher values of relative abundance for genes related with stress induced by glucose limitation and metal resistance to cadmium and other metals.
The category of plant interaction covers a wide range of different functional genes involved in microbial interactions with plants, including genes related to fungal function, soil benefit, soil borne pathogens, and virulence. The analyses performed showed significant differences (P < 0.1) in some subcategories, as shown in Figures S4G–J. There were not any significant differences in the genes in the categories of soil benefit or fungi function. Nevertheless, CT samples exhibited a higher significant relative abundance of detected genes from the oomycetes subcategory (soil borne pathogen), which included different genes from this pathogenic group. Genes related to virulence processes such as iron oxidation or secretion had a higher significant relative abundance for amended soils; whereas unamended soils exhibited significantly higher values for genes involved in virulence actions such as iron uptake (aerobactin genes) and pilin formation.
Unique DNA Probes Detected in as Suppressive Soil Samples
Results of the GeoChip analysis and the Venn diagram representation allowed us to determine microbial specific gene functions detected exclusively in each treatment and the number of commonly detected probes (27364) (Figure 6A). We found 6674 unique detected probes in CT samples and 2766 unique detected probes in AS samples (approximately 10% of the total AS detected genes) from the gene categories analyzed. Approximately 34.49% of the unique hybridizations were related to the Carbon cycle category (Figure 6B), mainly to starch and chitin degradation (Table S1). The Organic remediation gene category exhibited 14.53% unique hybridizations of genes related to the degradation of aromatic compounds. The Stress category had 13.38% unique hybridized probes and the Metal resistance category had 11.86% unique hybridized probes. The Nitrogen cycle category exhibited 8.57% unique hybridized probes, mostly in genes related with denitrification. The remaining gene categories had lower percentages: Sulfur cycle 5.60%, Fungi function 3.69%, Soil benefit 2.64% [approximately 44% of unique detected probes in this category correspond mainly with antimicrobial genes such as cat (catalase), phzF (phenazine), or pcbC (isopenicillin)], Phosphorus cycle 2.28%, Virulence 1.88%, and Soil borne pathogen 1.08% (Figure 6B). This analysis allowed us to relate different gene functions implicated in the metabolism of different soil compounds with bacterial or fungal classes present in the AS soil (Table S1).
Figure 6. GeoChip analysis of unique detected genes. (A) Number of core and unique detected genes (different gene ID) of amended soil (AS) and conventionally managed soil (CT). (B) Assigned functions of the unique genes detected in the AS sample.
The application of organic amendments to agricultural soils is a longstanding practice, and examples of organic-amendment-mediated suppression of soilborne diseases were reported as early as the late nineteenth century (Stone et al., 2004). Growers have observed that different types of organic materials suppress root rot for varying lengths of time. At present, nursery and greenhouse growers successfully use compost-amended potting mixes to suppress soilborne diseases, such as Pythium and Phytophthora root rots, in container systems (Hoitink et al., 1991). However, limited field studies have been conducted to determine the impact of soil amendments on microbial communities in actual organic and conventional production systems (Drinkwater et al., 1995; Gunapala and Scrow, 1998; Bulluck and Ristaino, 2001). In the case of avocado orchards, organic matter-mediated disease suppression against Phytophthora cinnamomi has been observed in avocado agricultural fields organically managed in Australia. Organic amendments (barley straw, sorghum residues, and native grass) were added to the soil under the trees as a mulch layer resulting in the suppression of Phytophthora root rot of avocado (Malajczuk, 1979, 1983). Additionally, our previous studies also demonstrated that different organic amendments can influence the composition and diversity of soil bacterial communities in avocado plants growing in microcosms after DGGE analysis, showing enhancement of specific populations such as Burkholderia and Frateuria (Bonilla et al., 2012a, 2015). Among different organic matter tested on avocado crops, composted almond shells (AS; commercial almond shells derived from the almond industry were piled and traditionally composted) exhibited enhancement of soil suppressiveness against R. necatrix (Bonilla et al., 2012a, 2015), the causal agent of avocado white root rot (Pliego et al., 2012). Even when soil suppressiveness against R. necatrix is improved after the addition of AS, only subtle changes in the bacterial community and composition and specific enzymatic activities have been reported using DGGE analysis (Bonilla et al., 2015). It must be considered that a wide range of factors can affect soil microorganism communities (van Veen et al., 1997). The soil samples used in our study came from the same orchard (same type of soil, environmental conditions, plant age, and cultivar, etc.), but were under different management, and this was assumed to be the only difference between the samples. The soil influenced by the amendment of AS showed some characteristics that differed from the conventional unamended soil. The almond shells are a high dry matter-containing substrate, composed of approximately 95% organic matter, with poor values of glucose, fructose, or sucrose. The characteristics and composition of AS makes this substrate an acceptable growing media for soilless culture (Valverde et al., 2013). Moreover, it must be taken into account that the avocado is a shallow rooted tree, with most of the feeder roots allocated in the top 15 cm, which needs good aeration. Roots are helped by the presence of a rich surface of organic mulch, as shown by the tendency of healthy feeder roots to grow into any decomposing litter layer (Chanderbali et al., 2013).
In this work, a metagenomic approach to the community composition of amended and unamended avocado soils have been performed for the first time. The use of metabarcoding and GeoChip techniques allowed a better knowledge on the community composition and their potential activities. In first place, an attempt to identify key factors involved in this enhanced suppressivity after the addition of organic amendments revealed the crucial role of the microbiota present in the organic amended soil. The microbiota evolved in the composted almond shells and was crucial for suppressiveness because the reduction of the bacterial population after a heat treatment in the organic amendment resulted in a more conducive phenotype (heat-treated soil samples harbor 105 cfu/g, most likely composed mainly by sporulated bacterial and fungal microorganisms). Moreover, total or partial suppressiveness was recovered when these heat-treated soil samples were complemented with a portion of soil influenced by AS, but it remained conducive when complemented with a portion of conventional soil (CT). This effect has been previously described for different suppressive soils, where sterilization by autoclaving, steam pasteurization, and irradiation rendered soils conducive to the pathogen studied (Malajczuk, 1983; Weller et al., 2002; Mendes et al., 2011). Suppressiveness experiments performed do not excluded the possibility that the disinfected avocado root used could harbor endophytic microorganisms, but our results significantly pointed out the role of the composted almond shells in the plant protection against R. necatrix. Thus, our results support the crucial role of microbes present in AS for turning the conducive CT soil into a more suppressive soil against R. necatrix.
To gain insights into the microbial diversity present in the soil samples, we used several different approaches. Phylogenetic marker analysis based on the sequencing of 16S rDNA and ITSs revealed a relatively similar array of prokaryotic and eukaryotic populations present in the AS and CT soil samples; however, a different response has been described in the literature for other types of organic matter from different sources, such as composted municipal waste (Zaccardelli et al., 2013). It is remarkable that in our model system, the group of fast-growing, easily cultivable Proteobacteria is the dominant group of prokaryotes in both soil samples. These data are similar to those previously observed for other soil and rhizosphere samples with a high presence of organic matter (Lynch and Whipps, 1990; Paul and Clark, 1996; Hawkes et al., 2007; Mendes et al., 2011). Moreover, the representation of the other phyla different than Proteobacteria were quite similar among AS-amended and unamended soils, thus contradicting the idea that a highly specific community is stimulated by the addition of AS. Diversity analysis confirmed the previously obtained results (Bonilla et al., 2015), highlighting the enhancement of specific microbial populations in AS-amended samples, such as Betaproteobacteria (Burkholderiales) and the class of Gammaproteobacteria, which have been reported to protect plants from fungal infections in other suppressive soils (Mendes et al., 2011). It is important to note the clear enhancement in AS-amended soil of the order Steroidobacter, previously reported to play an essential role in the positive interactions with plants; for example, controlling seed germination, stem, and root elongation or stress protection in plants (Zarraonaindia et al., 2015).
In contrast, analysis of eukaryotic ITS revealed a different abundance distribution of microbes among the two types of soil samples. Fungal clones were the most common and dominant microbial eukaryotes in the soil. AS-amended soil samples had an increased relative abundance of Ascomycota. This fact is not surprising considering that Ascomycetes are the largest group on true fungi (Larena et al., 1999). Moreover, the dominance of Ascomycota has been observed during different composting processes (De Gannes et al., 2013; Neher et al., 2013), where most of them are saprophytic and live on dead organic material that they help decompose (Agrios, 1997; Viebahn et al., 2005). This behavior easily explains their higher abundance when composted almond shells are added to the soil as mulch. Within Ascomycota, the group that exhibited the most apparent and highest increase of abundance in AS-amended soil samples was the fungal class of Dothideomycetes. A high abundance of Dothideomycetes in soils with at high hydrocarbon concentrations has been previously reported (Ferrari et al., 2011), suggesting its preference for those habitats with a high concentration of organic matter where it participates in biomass conversion (Shrestha et al., 2011). Moreover, the large increase of the phylum Pleosporales (Dothideomycetes) is also not surprising because this group is very well-known to contain species that chlorinate lignin as a first step of biomass conversion during plant litter degradation (Ortíz-Bermúdez et al., 2007). Interestingly, it has been shown that several genera of Dothideomycetes exhibit an increased presence in suppressive soils because they harbor endohyphal bacteria from groups that are capable of hydrocarbon biodegradation, such as the Xanthomonadales, Pseudomonadales, Burkholderiales, and Sphingomonadales (Hoffman and Arnold, 2010). Dothideomycetes have also been shown to increase slightly in AS-amended soils. However, the group that shows an apparent decrease in AS-amended soils is Mortierellales. This group has a complex phylogeny (Wagner et al., 2013) and is considered to be ubiquitous in the bulk and rhizospheric soil, implying that it could play a role in maintenance of the micro-ecological balance (Miao et al., in press). Interestingly, the group of Glomeromycota, which contains different groups of symbiotic fungi previously detected in avocado (Hass and Menge, 1990; González-Cortés et al., 2012), it is clearly detected in unamended soils, but decreased in the amended ones (below 1%). A possible explanation could be that in the AS amended soils, take place a strong competition with other decomposing fungi, such as the Dothideomycetes, more adapted to an environment with high amount of decomposing organic matter. Finally, it should be noted that members of Xylariaceae, to which R. necatrix belongs (Pliego et al., 2012), are less abundant in AS-amended soils, thus revealing a negative effect on this fungal group. These results indicate that the soil fungal community was affected by the soil amendment with AS.
Phylogenetic markers such as the prokaryotic 16S and eukaryotic ITS region do not carry explicit functional information. For this, the use of GeoChip-based analysis allowed for the analysis of microbial functional genes encoding key enzymes involved in major biogeochemical processes that facilitate linking microbial community structure to potential ecological functions (Torsvik and Ovreas, 2002). Using this technique, we screened potential functional gene diversity among unamended and AS-amended soil samples.
Probe signals and DCA analysis indicated that the microbial community functional structures differed between CT and AS soil samples. The sample sites are very close together, so the differences observed in the microbial communities are thought to be the result of amendment with organic matter.
Generally, similar abundance patterns of functional genes involved in nutrient cycling processes such a nitrogen, phosphorous or sulfur cycling, were found in both types of samples. However, AS-amended samples had higher signal intensities for C degradation (carbon cycle) genes than CT, with some differences being statistically significant. Substrates for this group of genes ranged from labile C to more recalcitrant C (e.g., starch, hemicelluloses, cellulose, chitin, and lignin). These results suggest that AS-amended microbial have a greater capacity for C degradation than CT communities. This suggests, as expected, an important role of carbon cycling in response to the addition of organic matter to the soil. However, no differences in gene abundance for N, P, or S cycling was observed. This can be explained because almond shells are a lignin-rich waste resulting from the almond industry, mostly composed of approximately 27% lignin and 73% holocellulose (Caballero et al., 1996), and those cycles were not compromised. However, statistical differences in the abundance of genes related to organic remediation and metal resistance were observed in AS-amended soil displaying lower levels than CT. This observation may be due to a decrease in the available compounds due to the high sorption ability of the composted almond shells and derivate compounds from its degradation, which have been previously reported to be able to remove such substances from the soil (Pehlivan et al., 2009).
Interestingly, both soil samples shared a core of probes corresponding to approximately 90% of the assayed sequences (27364 probes). However, approximately 10% of the total probes analyzed were unique for AS-amended samples (2766 probes). When the sequence of these probes were analyzed, they resulted in a very similar distribution to that previously shown for the whole GeoChip analysis, with above 34.5% corresponding to C cycling, followed by probes related to organic remediation (14.5%), stress (13.4%), metal resistance (11.9%), or the N cycle (8.6%). These results support the following previously described results: systems associated with organic matter-mediated general suppression; suppression typically occurs as a result of the activation of the indigenous microbial community (Lockwood, 1990); and suppressive activities can be generated by one to few populations of organisms (Gerlagh, 1968; Cook and Baker, 1983; Hoitink and Boehm, 1999; Weller et al., 2002). Postma et al. (2000) found that qualitative rather than quantitative shifts in the bacterial community correlate with disease suppressiveness, and several studies indicated that mechanisms within the microbial activity of the soil are responsible for the suppression of pathogens (Rovira and Wildermuth, 1981; Nitta, 1991; Workneh and van Bruggen, 1994; van Os and van Ginkel, 2001).
Among the specific taxa stimulated, Pseudomonadaceae, Burkholderiaceae, Xanthomonadales, and Actinobacteria, harbor genera and species with activity against plant pathogenic fungi (Postma et al., 2010). Additionally, it is important to note that Pseudomonas, Rhizobium, Bacillus, Variovorax, Phyllobacterium, and Azospirillum, are considered the most efficient plant growth-promoting bacteria (Bertrand et al., 2001).
Sequencing of specific probes present in AS-amended soils revealed the presence in such soil samples of genes for bacterial and fungal catalases, phenazine biosynthetic genes (from Proteobacteria) or the presence of potential antibiotics produced by Actinobacteria (data not shown). Nearly all of these probes corresponded to the GeoChip category “soil benefit,” where the antimicrobials from different groups were analyzed. To the best of our knowledge, no probes from Bacilli were used, so the role of antimicrobials such as iturin or fengicins, produced by Bacillus spp., cannot be discussed based on our results.
It is important to note that the genus Pseudomonas (class Gammaproteobacteria) and Bacillus (class Bacilli) are two of the most prominent bacteria that can be isolated from avocado soil and rhizosphere displaying antifungal activity and plant protection against soil-borne pathogens (Cazorla et al., 2006, 2007; González-Sánchez et al., 2010). Our results reinforce the importance of such microorganisms in the soil and root ecology of the avocado crop. These groups of microorganisms can produce metabolites, such as siderophores and antibiotics, with specific suppressive activity against soilborne pathogens. Antagonistic pseudomonads, including Pseudomonas chlororaphis, play a role in white root rot suppressiveness (Cazorla et al., 2006; Calderón et al., 2014). However, other types of rhizobacterial taxa may differ in prevalence between suppressive and conducive soils, suggesting that the microbial basis of white root rot could be far more complex than solely a Pseudomonas property; it has also been observed for other pathosystems such as Thielaviospsis basicola-mediated black root rot of tobacco (Almario et al., 2014).
In conclusion, and taking together the results obtained in this work and in previous works related, a theoretical model about the role of the microorganisms in enhancing suppressiveness after amendment with composted almond shells can be proposed (Figure S5). Soil amendments with composted almond shells resulted in an extra input of organic matter rich in lignin that could be initially degraded by fungal members of the community (such as Dothideomycetes) and Actinobacterias. Lignin degradation from composting almond shells would produce a progressive release to the soil of more simple compounds. Those compounds, together with others also present in the almond shells, could lead to an increase in carbon sources available, such as cellulose, hemicellulose, and aromatic compounds. At this point, some Proteobacteria already present in the soil (such as Gammaproteobacteria and Betaproteobacteria) could take advantage metabolizing that available organic matter, thus slightly enhancing their population. These groups of microorganisms could harbor, among other, genes involved in antifungal enzymatic activities and production of antimicrobial compounds that could have an effect on the interaction with other microbes. The resulting modified microbiota after addition of composted almond shells could be more active against some groups of phytopathogenic fungi (as Xilariales, where R. necatrix is included) finally showing a phenotype of induced suppressiveness effect.
FC and AD designed and corrected the manuscript. NB contributed to data management. CV performed date management and wrote the manuscript.
Conflict of Interest Statement
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.
We want to thank José M. Hermoso, Jorge González, and Emilio Guirado for assistance in suppressive assay design and management. We would like to thank Irene Linares and Claudia E. Calderón for their assistance during various parts of the projects. We are especially grateful to José M. Farré and Juan A. Torés for their helpful ideas and discussion, which were essential for this project. This work was supported by Plan Nacional I+D+I from Ministerio de Economia (MINECO) (AGL11-30354-C02-01 and AGL14-52518-C2-1-R), co-financed by FEDER funds (EU). CV was supported by a PhD fellowship from the FPI program of MINECO, and NB was supported by a PhD fellowship from the FPU program of MICINN.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2016.00004
Agrios, G. N. (ed.). (1997). Plant Pathology. San Diego, CA: Academic Press.
Almario, J., Muller, D., Défago, G., and Moënne-Loccoz, Y. (2014). Rhizosphere ecology and phytoprotection in soils naturally suppressive to Thielaviopsis black root rot of tobacco. Environ. Microbiol. 16, 1949–1960. doi: 10.1111/1462-2920.12459
Bai, S., Li, J., He, Z., van Nostrand, J. D., Tian, Y., Lin, G., et al. (2013). GeoChip-based analysis of the functional gene diversity and metabolic potential of soil microbial communities of mangroves. Appl. Microbiol. Biotechnol. 97, 7035–7048. doi: 10.1007/s00253-012-4496-z
Bailey, K. L., and Lazarovits, G. (2003). Suppressing soil-borne diseases with residue management and organic amendments. Soil Tillage Res. 72, 169–180. doi: 10.1016/S0167-1987(03)00086-2
Bender, G. S., Casale, W. L., and Rahimian, M. (1992). “Use of worm-composted sludge as a soil amendment for avocados in Phytophthora-infested soil,” in Proceeding of Second World Avocado Congress (Orange, CA), 143.
Bertrand, H., Nalin, R., Balli, E., and Cleyet-Maret, J. C. (2001). Isolation and indentification of the most efficient plant growth-promoting bacteria associated to canole (Brassica napus). Biol. Fertil. Soils 33, 152–156. doi: 10.1007/s003740000305
Bonanomi, G., D'Ascoli, R., Scotti, R., Gaglione, S. A., Caceres, M. G., Sultana, S., et al. (2014). Soil quality recovery and crop yield enhancement by combined application of compost and wood to vegetables grown under plastic tunnels. Agric. Ecosyst. Environ. 192, 1–7. doi: 10.1016/j.agee.2014.03.029
Bonilla, N., Cazorla, F. M., Martínez-Alonso, M., Hermoso, J. M., González-Fernández, J. J., Gaju, N., et al. (2012a). Organic amendments and land management affect bacterial community composition, diversity and biomass in avocado crop soils. Plant Soil 357, 215–226. doi: 10.1007/s11104-012-1155-1
Bonilla, N., Gutiérrez-Barranquero, J., Vicente, A., and Cazorla, F. (2012b). Enhancing soil quality and plant health through suppressive organic amendments. Diversity 4, 475–491. doi: 10.3390/d4040475
Bonilla, N., Vida, C., Martínez-Alonso, M., Landa, B. B., Gaju, N., Cazorla, F. M., et al. (2015). Organic amendments to avocado crops induce suppressiveness and influence the composition and activity of soil microbial communities. Appl. Environ. Microbiol. 81, 3405–3418. doi: 10.1128/AEM.03787-14
Bulluck, L. R. III, and Ristaino, J. B. (2001). Effect of synthetic and organic soil fertility amendments on southern blight, soil microbial communities, and yield of processing tomatoes. Phytopathology 92, 181–189. doi: 10.1094/PHYTO.2002.92.2.181
Caballero, J. A., Font, R., and Marcilla, A. (1996). Comparative study of the pyrolysis of almond shells and their fractions, holocellulose and lignin. Product yields and kinetics. Thermochim. Acta 276, 57–77. doi: 10.1016/0040-6031(95)02794-7
Calderón, C. E., de Vicente, A., and Cazorla, F. M. (2014). Role of 2-hexyl, 5-propyl resorcinol production by Pseudomonas chlororaphis PCL1606 in the multitrophic interactions in the avocado rhizosphere during the biocontrol process. FEMS Microbiol. Ecol. 89, 20–31. doi: 10.1111/1574-6941.12319
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Cazorla, F. M., Duckett, S. B., Bergström, E. T., Noreen, S., Odijk, R., Lugtenberg, B. J. J., et al. (2006). Biocontrol of avocado dematophora root rot by antagonistic Pseudomonas fluorescens PCL1606 correlates with the production of 2-hexyl 5-propyl resorcinol. Mol. Plant Microbe Interact. 19, 418–428. doi: 10.1094/MPMI-19-0418
Cazorla, F. M., Romero, D., Pérez-García, A., Lugtenberg, B. J. J., de Vicente, A., and Bloemberg (2007). Isolation and characterization of antagonistic Bacillus subtilis strains from the avocado rhizoplane displaying biocontrol activity. J. Appl. Microbiol. 103, 1950–1959. doi: 10.1111/j.1365-2672.2007.03433.x
Chanderbali, A. S., Soltis, D. E., Soltis, P. S., and Wolstenholme, B. N. (2013). “Taxonomy and Botany,” in The Avocado: Botany, Production and Uses, 2nd Edn, eds B. A. Schaffer, B. N. Wolstenholme, and A. W. Whiley (Oxfordshire: CABI), 31–50. doi: 10.1079/9781845937010.0031
Cong, J., Yang, Y., Liu, X., Lu, H., Liu, X., Zhou, J., et al. (2015). Analyses of soil microbial community compositions and functional genes reveal potential consequences of natural forest succession. Sci. Rep. 5:10007. doi: 10.1038/srep10007
Cook, R. J., and Baker, F. (1983). The Nature and Practice of Biological Control of Plant Pathogens. Saint Paul, MN: The American Phytopathological Society.
Cretoiu, M. S., Korthals, G. W., Visser, J. H. M., and van Elsas, J. D. (2013). Chitin amendment increases soil suppressiveness toward plant pathogens and modulates the actinobacterial and oxalobacteraceal communities in an experimental agricultural field. Appl. Environ. Microbiol. 79, 5291–5301. doi: 10.1128/AEM.01361-13
De Gannes, V. D., Endoxic, G., and Hickey, W. J. (2013). Insight into fungal communities in compost revealed by 454-pyrosequencing: implications for human health and safety. Front. Microbiol. 4:164. doi: 10.3389/fmicb.2013.00164
Deng, Y., and Zhou, J. (2013). “Molecular ecological network of microbial communities,” in Encyclopedia of Metagenomics, ed K. Nelson (Berlin; Heidelberg: Springer-Verlag), 1–8. doi: 10.1007/978-1-4614-6418-1_758-2
Downer, A. J., Menge, J. A., and Pond, E. (2001). Association of cellulolytic enzyme activities in Eucalyptus mulches with biological control of Phytophthora cinnammomi. Phytopathology 91, 847–855. doi: 10.1094/PHYTO.2001.91.9.847
Drinkwater, L. E., Letourneau, D. K., Workneh, F., van Bruggen, A. H. C., and Shennan, C. (1995). Fundamental differences between conventional and organic tomato agroecosystems in California. Ecol. Appl. 5, 1098–1112. doi: 10.2307/2269357
Ferrari, B. C., Zhang, C., and van Dorst, J. (2011). Recovering greater fungal diversity from pristine and diesel fuel contaminated sub-Antartic soil through cultivation using both a high and a low nutrient media approach. Front. Microbiol. 2:217. doi: 10.3389/fmicb.2011.00217
Freeman, S., Sztejnberg, A., and Chet, I. (1986). Evaluation of Trichoderma as a biocontrol agent for Rosellinia necatrix. Plant Soil 94, 163–170. doi: 10.1007/BF02374340
Garbeva, P., Veen, J. A., and van Elsas, J. D. (2004). Microbial diversity in soil: selection of microbial populations by plant and soil type and implications for disease suppressiveness. Ann. Rev. Phytopathol. 42, 243–270. doi: 10.1146/annurev.phyto.42.012604.135455
Gerlagh, M. (1968). Introduction of Ophiobolus graminis into new polders and its decline. Neth. J. Plant Pathol. 74, 1–97. doi: 10.1007/BF02019999
González-Cortés, J. C., Vega-Fraga, M., Varela-Fregoso, L., Martínez-Trujillo, M., Carreón-Abud, Y., and Gavito, M. E. (2012). Arbuscular mycorrhizal fungal (AMF) communities and land use change: the conversion of temperate forest to avocado plantations and maize fields in central Mexico. Fungal Ecol. 5, 16–23. doi: 10.1016/j.funeco.2011.09.002
González-Sánchez, M. A., de Vicente, A., Pérez-García, A., Pérez-Jiménez, R., Romero, D., and Cazorla, F. M. (2013). Evalation of the effectiveness of biocontrol bacteria against avocado white root rot occurring under commercial greenhouse plant production conditions. Biol. Cont. 67, 94–100. doi: 10.1016/j.biocontrol.2013.08.009
González-Sánchez, M. A., Pérez-Jiménez, R. M., Pliego, C., Ramos, C., de Vicente, A., and Cazorla, F. M. (2010). Biocontrol bacteria selected by a direct plant protection strategy against avocado white root rot show antagonism as a prevalent trait. J. Appl. Microbiol. 109, 65–78. doi: 10.1111/j.1365-2672.2009.04628.x
Gunapala, N., and Scrow, K. M. (1998). Dynamics of soil microbial biomass and activity in conventional and organic farming systems. Soil Biol. Biochem. 30, 805–816. doi: 10.1016/S0038-0717(97)00162-4
Hadar, Y., and Papadopoulou, K. K. (2012). Suppressive composts: microbial ecology links between abiotic environments and healthy plants. Annu. Rev. Phytopathol. 50, 133–153. doi: 10.1146/annurev-phyto-081211-172914
Hass, J. H., and Menge, J. A. (1990). VA-mycorrhizal fungi and soil characteristics in avocado (Persea Americana Mill.) orchard soils. Plant Soil 127, 207–212. doi: 10.1007/BF00014427
Hawkes, C. V., DeAngelis, K. M., and Firestone, M. K. (2007). “Root interactions with soil microbial communities and processes,” in The Rhizosphere, Vol. 1, eds Z. Cardon, and J. Whitbeck (New York, NY: Elsevier), 1–29.
Hazen, T. C., Dubinsky, E. A., DeSantis, T. Z., Andersen, G. L., Piceno, Y. M., Singh, N., et al. (2010). Deep-sea oil plume enriches indigenous oil-degrading bacteria. Science 330, 204–208. doi: 10.1126/science.1195979
He, Z., Deng, Y., and Zhou, J. (2012). Development of functional gene microarrays for microbial community analysis. Curr. Opin. Biotechnol. 23, 49–55. doi: 10.1016/j.copbio.2011.11.001
Hermoso, J. M., Torres, M. D., and Farré, J. M. (2011). “Effects of organic materials on adult Hass trees,” in Proceedings of the Seventh World Avocado Congress (Cairns, QLD).
Hirsch, P. R., Mauchline, T. H., and Clark, I. M. (2013). “Culture-independent molecular approaches to microbial ecology in soil and the rhizosphere,” in Molecular Microbial Ecology of the Rhizosphere, Vol. 5, ed F. J. de Bruinjn (Hoboken, NJ: John Wiley & Sons, Inc). 45–55. doi: 10.1002/9781118297674.ch5
Hoffman, M. T., and Arnold, A. E. (2010). Diverse bacteria inhabit living hyphae of phylogenetically diverse fungal endophytes. Appl. Environ. Microbiol. 76, 4063–4075. doi: 10.1128/AEM.02928-09
Hoitink, H., and Boehm, M. (1999). Biocontrol within the context of soil microbial communities: a substrate-dependent phenomenon. Annu. Rev. Phytopathol. 37, 427–446. doi: 10.1146/annurev.phyto.37.1.427
Hoitink, H., Stone, A. G., and Han, D. Y. (1991). Suppression of plant diseases by composts. HortScience 32, 184–187.
Janvier, C., Villeneuve, F., Alabouvette, C., Edel-Hermann, V., Mateille, T., and Steinberg, C. (2007). Soil health through soil disease suppression: which strategy from descriptors to indicators? Soil Biol. Biochem. 39, 1–23. doi: 10.1016/j.soilbio.2006.07.001
Koyama, A., Wallenstein, M., Simpson, R. T., and Moore, J. C. (2014). Soil bacterial community composition altered by increased nutrient availability in Arctic tundra soils. Front. Microbiol. 5:516. doi: 10.3389/fmicb.2014.00516
Larena, I., Salazar, O., González, V., Julián, M. C., and Rubio, V. (1999). Design of a primer for ribosomal DNA internal transcriber spacer with enhanced specificity for Ascomycetes. J. Biotech. 75, 187–194. doi: 10.1016/S0168-1656(99)00154-6
Lazarovits, G., Tenuta, M., and Conn, K. (2001). Organic amendments as a disease control strategy for soilborne diseases of high-value agricultural crops. Australas. Plant. Pathol. 30, 111–117. doi: 10.1071/AP01009
Liang, Y., He, Z., Wu, L., Deng, Y., Li, G., and Zhou, J. (2010). Development of a common oligonucleotide reference standard for microarray data normalization and comparison across different microbial communities. Appl. Environ. Microbiol. 76, 1088–1094. doi: 10.1128/AEM.02749-09
Lockwood, J. L. (1990). “Relation of energy stress to behavior of soil-borne plant pathogens and to disease development,” in Biological Control of Soil-borne Plant Pathogens, eds D. Hornby and R. J. Cook (Wallinford, CT: CAB International), 197–214.
López, R., Burgos, P., Hermoso, J. M., Hormaza, J. I., and González-Fernández, J. J. (2014). Long term changes in soil properties and enzyme activities after almond shell mulching in avocado organic production. Soil Tillage Res. 143, 155–163. doi: 10.1016/j.still.2014.06.004
López-Herrera, C. J., Pérez-Jiménez, R. M., Zea-Bonilla, T., Basallote-Ureba, M. J., and Melero-Vara, J. M. (1998). Soil solarization in established avocado trees for control of Dematophora necatrix. Plant Dis. 82, 1088–1092. doi: 10.1094/PDIS.1918.104.22.1688
López-Herrera, C. J., and Zea-Bonilla, T. (2007). Effects of benomyl, carbendazim, fluazinam and thiophanate methyl on white root rot of avocado. Crop Prot. 26, 1186–1192. doi: 10.1016/j.cropro.2006.10.015
Lynch, J. M., and Whipps, J. M. (1990). Substrate flow in the rhizosphere. Plant Soil 129, 1–10. doi: 10.1007/BF00011685
Malajczuk, N. (1979). “Biological suppression of Phytophthora cinnamomi in eucalypts and avocados in Australia,” in Soil-borne Plant Pathogens, eds B. Schippers and W. Gams (London: Academic Press), 635–652.
Malajczuk, N. (1983). “Microbial antagonism to Phytophthora,” in Phytophthora, Its Biology, Taxonomy, Ecology and Pathology, eds D. C. Erwin, S. Bartnicki-Garcia, and P. H. Tsao (Saint Paul, MN: American Phytopathological Society), 197–218.
Mendes, R., Kruijt, M., de Bruijn, I., Dekkers, E., van der Voort, M., Shneider, J. H. M., et al. (2011). Deciphering the rhizosphere microbiome for disease-suppressive bacteria. Science 332, 1097–1100. doi: 10.1126/science.1203980
Miao, C. P., Mi, Q. L., Qiao, X. G., Zheng, Y. K., Chen, Y. W., and Xu, L. (in press). Rhizospheric fungi of Panax notoginseng: diversity antagonism to host phytopathogens. J. Ginseng Res. doi: 10.1016/j.jgr.2015.06.004
Moore-Gordon, C., Cowan, A. K., and Wolstenholme, B. N. (1997). “Mulching of avocado orchards to increase Hass yield and fruit size and boost financial rewards a three season summary of research findings,” in South African Avocado Growers' Association Yearbook, Vol. 20 (Tzanee: SAAGA) 46–49.
Muyzer, G., Brinkhoff, T., Nübel, U., Santegoeds, C., Schäfer, H., and Wawer, C. (2004). “Denaturing gradient gel electrophoresis (DGGE) in microbial ecology,” in Molecular Microbial Ecology, Vol. 1, eds G. A. Kowalchuk, F. J. Bruinjn, I. M. Heas, A. D. L. Akkermans, and J. D. van Elsas (Dordrecht: Kluwer Academic Publications), 743–770.
Neher, D. A., Weicht, T. R., Bates, S. T., Left, J. W., and Fierer, N. (2013). Changes in bacterial and fungal communities across compost recipes, preparation method, and composting times. PLoS ONE 8:e79512. doi: 10.1371/journal.pone.0079512
Nitta, T. (1991). Diversity of root fungal floras: its implications for soil-borne diseases and crop growth. Jpn. Agric. Res. 25, 6–11.
Ortíz-Bermúdez, P., Hirth, K. C., Srebotnik, E., and Hammel, K. E. (2007). Chlorination of lignin by ubiquitous fungi has a likely role in global organichloride production. Proc. Natl. Acad. Sci. U.S.A. 104, 3895–3900. doi: 10.1073/pnas.0610074104
Pane, C., Piccolo, A., Spaccini, R., Celano, G., Villecco, D., and Zaccardelli, M. (2013). Agricultural waste-based composts exhibiting suppressivity to diseases caused by the phytopathogenic soil-borne fungi Rhizoctonia solani and Sclerotinia minor. Appl. Soil Ecol. 65, 43–51. doi: 10.1016/j.apsoil.2013.01.002
Paul, E. A., and Clark, F. E. (1996). Soil Microbiology and Biochemistry, 2nd Edn. San Diego, CA: Academic Press, Inc.
Pehlivan, E., Altun, T., Cetin, S., and Bhanger, M. I. (2009). Lead sorption by waste biomass of hazelnut and almond shell. J. Hazard. Mater. 167, 1203–2128. doi: 10.1016/j.jhazmat.2009.01.126
Peiffer, J. A., Spor, A., Koren, O., Jin, Z., Green, S., and Dangl, J. L. (2013). Diversity and heritability of the maize rhizosphere microbiome under field conditions. Proc. Natl. Acad. Sci. U.S.A. 110, 6548–6553. doi: 10.1073/pnas.1302837110
Pliego, C., Kanematsu, S., Ruano-Rosa, D., de Vicente, A., López-Herrera, C., Cazorla, F. M., et al. (2009). GFP sheds light on the infection process of avocado roots by Rosellinia necatrix. Fungal Genet. Biol. 46, 137–145. doi: 10.1016/j.fgb.2008.11.009
Pliego, C., López-herrera, C., Ramos, C., and Cazorla, F. M. (2012). Developing tools to unravel the biological secrets of Rosellinia necatrix, an emergent threat to woody crops. Mol. Plant Pathol. 13, 226–239. doi: 10.1111/j.1364-3703.2011.00753.x
Postma, J., Sheper, R. W. A., and Schilder, M. T. (2010). Effect of successive cauliflower plantings and Rhizoctonia solani AG 2-1 inoculations on disease suppressiveness of a suppressive and conducive soil. Soil Biol. Biochem. 42, 804–812. doi: 10.1016/j.soilbio.2010.01.017
Postma, J., Willemsen-de Klein, M. J., and van Elsas, J. D. (2000). Effect of the indigenous microflora on the development of root and crown rot caused by Pythium aphanidermatum in cucumber grown on rockwool. Phytopathology 90, 125–133. doi: 10.1094/PHYTO.2000.90.2.125
Rovira, A. D., and Wildermuth, G. B. (1981). “The nature and mechanisms of suppression,” in Biology and Control of Take All, eds M. J. C. Asher and P. Shipton (London: Academic), 385–415.
Ruano-Rosa, D., and López-Herrera, C. J. (2009). Evaluation of Trichoderma spp. As biocontrol agents against avocado white root rot. Biol. Cont. 51, 66–71. doi: 10.1016/j.biocontrol.2009.05.005
Ruano-Rosa, D., Martín-Pérez, R., Bonilla, N., Cazorla, F. M., de Vicente, A., and López-Herrera, C. J. (2014). Biological control of avocado white root rot with combined applications of Trichoderma spp. and rhizobacteria. Europ. Plant Pathol. 138, 751–762. doi: 10.1007/s10658-013-0347-8
Shrestha, K., Shrestha, P., Walsh, K. B., Harrower, K. M., and Midmore, D. J. (2011). Microbial enhancement of compost extracts based on cattle rumen content. compost-characterisation of a system. Bioresour. Technol. 102, 8027–8034. doi: 10.1016/j.biortech.2011.06.076
Stone, A. G., Scheuerell, S. J., and Darby, H. M. (2004). “Suppression of soilborne diseases in field agricultural systems: organic matter management, cover cropping, and other cultural practices,” in Soil Organic Matter in Sustainable Agriculture, eds F. Magdoff and R. R. Well (Boca Ratón, FL: CRC Press), 131–177. doi: 10.1201/9780203496374.ch5
Termorshuizen, A. J., van Rijn, E., van der Gaag, D. J., Alabouvette, C., Chen, Y., Lagerföf, J., et al. (2006). Suppressiveness of 18 composts against 7 pathosystems: variability in pathogen response. Soil Biol. Biochem. 38, 2461–2477. doi: 10.1016/j.soilbio.2006.03.002
Torsvik, V., and Ovreas, L. (2002). Microbial diversity and function in soil: from genes to ecosystems. Curr. Opin. Microbiol. 5, 240–224. doi: 10.1016/S1369-5274(02)00324-7
Tu, Q., Yu, H., He, Z., Deng, Y., Wu, L., van Nostrand, J. D., et al. (2014). GeoChip 4: a functional gene-array-based high-throughput environmental technology for microbial community analysis. Mol. Ecol. Resour. 14, 914–928. doi: 10.1111/1755-0998.12239
Valverde, M., Madrid, R., García, A. L., del Amor, F. M., and Rincón, L. (2013). Use of almond shell and almond hull as substrate for sweet pepper cultivation. Effect on fruit yield and mineral content. Span. J. Agric. Res. 11, 164–172. doi: 10.5424/sjar/2013111-3566
van Elsas, J. D., Costa, R., Jansson, J., Sjöling, S., Bailey, M., Nalin, R., et al. (2008). The metagenome of disease suppressive soils—experiences from the METACONTROL project. Trends Biotecnol. 26, 591–601. doi: 10.1016/j.tibtech.2008.07.004
van Elsas, J. D., Jansson, J. K., and Trevors, T. (eds.). (2007). Modern Soil Microbiology II. Boca Raton, FL: CRC Press, Taylor & Francis Group, LLC.
van Elsas, J. D., and Postma, J. (2007). “Suppression of soil-borne phytopathogens by compost,” in Compost Science and Technology, Vol. 8, eds L. F. Diaz, M. de Bertoldi, W. Bidlingmaier, and E. Stentiford (Amsterdam: Elsevier), 201–214.
van Os, G. J., and van Ginkel, J. H. (2001). Suppression of Pythium root rot in bulbous iris in relation to biomass and activity of the soil microflora. Soil Biol. Biochem. 32, 1447–1454. doi: 10.1016/S0038-0717(01)00053-0
van Veen, J. A., van Oberbeek, L. S., and van Elsas, J. D. (1997). Fate and activity of microorganisms introduced into soil. Microbiol. Mol. Biol. Rev. 61, 121–135.
Viebahn, M. C., Veeman, C., Wernars, K., van Loon, L. C., Smith, E., and Bakker, P. A. H. M. (2005). Assessment of differences in ascomycete communities in the rhizosphere of field-grown wheat and potato. FEMS Microbiol. Ecol. 53, 245–253. doi: 10.1016/j.femsec.2004.12.014
Wagner, L., Stielow, B., Hoffmann, K., Petkovits, T., Papp, T., Vágvölgyi, C., et al. (2013). A comprehensive molecular phylogeny of the Mortierellales (Mortierellomycotina) based on nuclear ribosomal DNA. Persoonia 30, 77–93. doi: 10.3767/003158513X666268
Weller, D. M., Raaijmakers, J. M., Gardener, B. B. M., and Thomashow, L. S. (2002). Microbial populations responsible for specific soil suppressiveness to plant pathogens. Annu. Rev. Phytopathol. 40, 309–348. doi: 10.1146/annurev.phyto.40.030402.110010
Wolstenholme, B. N., Moore-Gordon, C. S., and Cowan, A. K. (1997). “Orchard mulching effects on avocado fruiting,” in Proceedings from the Australian Avocaco Growers' Association Conference'97: Searching for Quality, ed I. G. Cutting (Rotorna), 119–130.
Workneh, F., and van Bruggen, A. H. C. (1994). Microbial densitiy, composition, and diversity in organically and conventionally managed rhizosphere soil in relation to suppression of corky root of tomatoes. Appl. Soil Ecol. 1, 219–230. doi: 10.1016/0929-1393(94)90013-2
Xie, J., He, Z., Liu, X., Liu, X., Van Nostrand, J. D., Deng, Y., et al. (2011). GeoChip-based analysis of the functional gene diversity and metabolic potential of microbial communities in acid mine drainage. Appl. Environ. Microbiol. 77, 991–999. doi: 10.1128/AEM.01798-10
Xue, K., Wu, L., Deng, Y., He, Z., van Nostrand, J., Robertson, P. G., et al. (2013). Functional gene differences in soil microbial communities from conventional, low-input, and organic farmlands. Appl. Environ. Microbiol. 79, 1284–1292. doi: 10.1128/AEM.03393-12
Zaccardelli, M., de Nicola, F., Villecco, D., and Scotti, R. (2013). The development and suppressive activity of soil microbial communities under compost amendment. J. Soil Sci. Plant Nutr. 13, 730–742. doi: 10.4067/s0718-95162013005000058
Zarraonaindia, I., Owens, S. M., Weisenhorn, P., West, K., Hampton-Marcell, J., Lax, S., et al. (2015). The soil microbiome influences grapevine-associated microbiota. MBio 6:e02527–14. doi: 10.1128/mBio.02527-14
Zhang, Y., Lu, Z., Liu, S., Yang, Y., He, Z., Ren, Z., et al. (2013). GeoChip-based analysis of microbial communities in alpine meadow soils in the Qinghai-Tibetan plateau. BMC Microbiol. 13:72. doi: 10.1186/1471-2180-13-72
Keywords: soil, amendment, almond shells, microbial profiling, suppressiveness
Citation: Vida C, Bonilla N, de Vicente A and Cazorla FM (2016) Microbial Profiling of a Suppressiveness-Induced Agricultural Soil Amended with Composted Almond Shells. Front. Microbiol. 7:4. doi: 10.3389/fmicb.2016.00004
Received: 22 October 2015; Accepted: 05 January 2016;
Published: 22 January 2016.
Edited by:Aurelio Ciancio, Istituto per la Protezione Sostenibile delle Piante del CNR, Italy
Reviewed by:Raffaella Balestrini, Consiglio Nazionale delle Ricerche, Italy
Giancarlo Roccuzzo, Consiglio per la Ricerca in Agricoltura e l'Analisi dell'Economia Agraria, Italy
Copyright © 2016 Vida, Bonilla, de Vicente and Cazorla. 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) or licensor 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: Francisco M. Cazorla, firstname.lastname@example.org