Araceli Lamelas1†
Damaris Desgarennes1†
Daniel López-Lima1
Luc Villain2
Alexandro Alonso-Sánchez1
Alejandro Artacho3
Amparo Latorre3,4,5
Andrés Moya3,4,5
Gloria Carrión1*- 1Red de Estudios Moleculares Avanzados and Red de Biodiversidad y Sistemática, Instituto de Ecología A. C., Xalapa, Mexico
- 2CIRAD, UMR IPME, Montpellier, France
- 3Joint Unit of Research in Genomics and Health, Foundation for the Promotion of Health and Biomedical Research in the Valencian Community (FISABIO) and Cavanilles Institute of Biodiversity and Evolutionary Biology, Universitat de València, Valencia, Spain
- 4CIBER en Epidemiología y Salud Pública, Madrid, Spain
- 5Institute for Integrative Systems Biology (I2SysBio), University of Valencia, Spanish National Research Council (CSIC-UVEG), Valencia, Spain
The Meloidogyne-based disease complexes (MDCs) are caused by the interaction of different root-knot nematode species and phytopathogenic fungi. These complexes are devastating several important crops worldwide including tomato and coffee. Despite their relevance, little is known about the role of the bacterial communities in the MDCs. In this study 16s rDNA gene sequencing was used to analyze the bacterial microbiome associated with healthy and infested roots, as well with females and eggs of Meloidogyne enterolobii and M. paranaensis, the causal agents of MDC in tomato and coffee, respectively. Each MDC pathosystems displayed a specific taxonomic diversity and relative abundances constituting a very complex system. The main bacterial drivers of the MDC infection process were identified for both crops at order level. While corky-root coffee samples presented an enrichment of Bacillales and Burkholderiales, the corcky-root tomato samples presented an enrichment on Saprospirales, Chthoniobacterales, Alteromonadales, and Xanthomonadales. At genus level, Nocardia was common to both systems, and it could be related to the development of tumor symptoms by altering both nematode and plant systems. Furthermore, we predicted the healthy metabolic profile of the roots microbiome and a shift that may result in an increment of activity of central metabolism and the presence of pathogenic genes in both crops.
Introduction
Plants harbor a myriad of microorganisms out and inside their tissues, constituting specific microbiomes. Plant-associated microorganisms, mainly endophytes, play an important role in plant health and productivity, since they provide their hosts with essential functions (Kiers et al., 2003: Kiers and van der Heijden, 2006; Bulgarelli et al., 2012; Vandenkoornhuyse et al., 2015). Pathogens attack, results in the disturbance of the plant microbiome balance, possibly promoting diseases through cooperation or coinfection of different pathogens (Tollenaere et al., 2016; Vannier et al., 2019). The idea that disease could be caused by the combined action of more than one pathogen set the bases for the pathobiome concept, which describes the pathogen as a component of a complex microbial community (Ryan, 2013; Vayssier-Taussat et al., 2014; Lamichhane and Venturi, 2015; Jakuschkin et al., 2016; Brader et al., 2017).
The pathobiome theory considers the interaction between well adapted pathogens, known as keystone pathogens that modulate plant defenses and strongly affect plant microbiome composition; accessory pathogens, which find their niche in the pathobiome by providing nutritional or colonization support to keystone pathogens; and pathobionts, whose commensal or mutualistic relationship switches to pathogenic, in response to changes in the plant microbiome (Brader et al., 2017).
Considering this definition, the pathobiome concept contradicts the Koch and Hill's fundamental postulate “one microbe-one disease” (Ryan, 2013; Vayssier-Taussat et al., 2014), stating that disease is driven by a highly complex network in which microorganisms interact with each other influenced by diverse external factors (Lamichhane and Venturi, 2015; Brader et al., 2017).
Sedentary endoparasitic nematodes such as Meloidogyne spp. are very harmful plant parasites due to their high pathogenicity and their ability to overcome resistance in many hosts (Jones et al., 2013; Elling, 2013). Interactions of root-knot nematodes (RKN) with other phytopathogenic microorganisms can cause syndromes known as Meloidogyne-based disease complex (MDC), severely affecting world agricultural production (Wolfgang et al., 2019). MDC symptoms are characterized by the initial infection of RKN, causing the typical root swellings and knots. As the infection progresses, the concomitant presence of fungi and bacteria is observed along with nematodes in root tissues, which present more severe symptoms such as hyperplasia with deep and cracked cortical root tissues, having a cork-like appearance, leading to necrosis and atrophy of the root system (Koike et al., 2006; Taher et al., 2017; López-Lima et al., in press; Wolfgang et al., 2019). These symptoms, in the advanced stages of infection, are also commonly known as corky-root disease in coffee (Bertrand et al., 2000; Lopez-Lima et al., 2015; Villain et al., 2018). Consequently, water and nutrient flows are altered, causing significant damage to the plant including chlorosis, defoliation, necrosis of tip branches, and decreased production (Hajji-Hedfi et al., 2018).
The bacterial community is an important component of MDC, perhaps favoring disease development. Some groups of bacteria like Herbaspirillum, Rickettsia, Chitinophaga, and Pedobacter have been reported as symbionts of nematodes (Tian et al., 2011; Baquiran et al., 2013; Cheng et al., 2013). In contrast, other bacterial species such as Pseudomonas and Commamonas have been reported reducing nematode populations by inhibiting hatching or by nematicidal action (Siddiqui et al., 2006; Wolfgang et al., 2019). Root bacterial microbiome structure has been studied in some relevant crops affected by RKN and other pathogens with the aim of understanding their function in the agroecosystems and developing better management strategies (Busby et al., 2017; Finkel et al., 2017; Toju et al., 2018; Toju and Tanaka, 2019).
Particularly, severe effects of MDC have been reported in two crops of world importance: tomato (Sasanelli et al., 2008; Ganaie and Khan, 2011; Rizvi et al., 2015; Mwangi et al., 2019) and coffee (Alvarado, 1935; Carneiro et al., 1996; Carneiro et al., 2004; Villain et al., 2013; Lopez-Lima et al., 2015; Monteiro et al., 2016; Villain et al., 2018). In tomato MDC is caused by M. incognita and/or M. enterolobii interacting with F. oxysporum, F. solani, Pyrenochaeta lycopersici, and Rhizoctonia solani (Sasanelli et al., 2008; Ganaie and Khan, 2011; Rizvi et al., 2015; van Bruggen et al., 2016; Mwangi et al., 2019). In Mexico, M. enterolobii causes significant losses in due to its high virulence, which leads to plant death before completing the crop cycle. This is the case of cultivars carrying the resistant gene Mi-1, which confers resistance to M. javanica, M. arenaria and M. incognita (Kiewnick et al., 2008; Rosa et al., 2014; Martínez-Gallardo et al., 2015; Martínez-Gallardo et al., 2019).
In coffee, MDC is caused by the RKN Meloidogyne paranaensis or M. arabicida interacting with phytopathogenic fungi such as Fusarium oxysporum and F. solani (Bertrand et al., 2000; López-Lima et al., in press). In Mexico, Brazil and Guatemala, M. paranaensis is the RKN responsible for MDC, affecting both Arabica and Robusta crops (Carneiro et al., 1996; Carneiro et al., 2004; Villain et al., 2013; Lopez-Lima et al., 2015; Villain et al., 2018).
Notwithstanding, little is known about the bacteria associated with these MDCs and their role in the development of the disease in tomato and coffee crops. M. enterolobii and M. paranaesis are emerging species causing severe damage in tomato and coffee crops, respectively (Elling, 2013). These two pathosystems differ greatly, since tomato is a short-cycle crop mainly grown under greenhouse conditions and subjected to intensive management, while coffee is a woody perennial crop mainly grown within agroforestry systems.
In order to understand the pathogenicity process of a plant disease, the pathobiome components must be characterized. We are particularly interested in investigating the bacterial microbiome associated with the diseases caused by a root-knot nematode-fungi complex: Meloidogyne-based disease complex in tomato and coffee. The study aims were to: i) investigate the bacterial community of tomato and coffee, infested by MDC, and ii) determine whether bacterial community composition and function differ between healthy and infected tissues, as well as in different life stages involved in RKN and rhizospheric soil. This work will allow us to establish a fundamental understanding of the dynamics and ecological role of microbes in the pathobiome of Meloidogyne-based disease complex.
Materials and Methods
Study System and Sample Collection
Tomato (Solanum lycopersicum) samples were collected from a five months-old crop, cultivar Red Grape, located in Culiacan, Sinaloa, Mexico (24°55'59” N, 107°26'34” W at 67m.a.s.l). Coffee (Coffea arabica) samples were collected from a four years-old plantation, cultivar Costa Rica 95 in Jilotepec, Veracruz, Mexico (19°35'42” N, 93°53'01” W at 996 m.a.s.l). Plants were classified according to the progress of MDC symptoms and the root galling index (RGI, Coyne and Ross, 2014) from 1: no visible galling damage to 5: severe galling damage. 1) No symptoms: roots with no galling damage (RGI level 1, healthy root tissues), 2) Initial symptoms: roots with slight galling damage (RGI level 2), separated knots but without rot or necrosis (root nodule tissue), and 3) advanced symptoms: roots with heavy galling damage (RGI level 4–5), hyperplasia, continuous nodules, with deep and cracked cortical tissues and presence of rot and necrosis (corky-root tissues). Five plants were collected per plant species, from them healthy root tissue, root nodule tissue, corky-root tissue (in this study, corky-root refers to severely damaged root tissues in the last stage of the disease), female nematodes, egg masses, and rhizospheric soil were obtained (Figure S1, Table S1). Female nematodes and egg masses were obtained by dissecting corky roots. Nematodes were identified by species-specific SCAR markers (Randig et al., 2002; Tigano et al., 2010). Meloidogyne enterolobii and M. paranaensis were identified from tomato and coffee samples, respectively.
Females and eggs were washed with PBS (137 mM NaCl, 2.7 mM KCl, 10 mM Na2HPO4, pH 7.4) by slow pipetting. PBS females and eggs washing solution, containing external bacteria, was recovered (Figure S1, Table S1). All samples were stored at −20°C until processing.
DNA Extraction, 16S rRNA Gene Amplification, and Sequencing
Female nematodes or eggs were ground using a sterile pestle in 1.5 mL tube after adding 1 mL of cold Extraction Buffer I (10 mM Tris, 60 mM NaCl, 5% Sucrose, 10 mM EDTA, pH = 7.8, sterilized and stored at 4°C). Then, 200 μL of Extraction Buffer II (300 mM Tris, 1.25% SDS, 5% Sucrose, 10 mM EDTA pH = 8.0) was added and mixed. Samples were then incubated for 30 min at 65°C. Afterward, 60 μL 3 M Potassium Acetate (pH = 4.8 to 5.2) were added, and the samples were incubated at −20 °C for 20 min, thawed at room temperature for 5 min and centrifuged at maximum speed. The recovered pellet was washed with 70% ethanol. Samples were air dried and resuspended in 35 μL of sterile molecular grade water. Root and soil DNA extraction was conducted using the Quick-DNA ™Fecal/Soil Microbe Miniprep kit (Zymo Research, Cat. No. D6010) following the manufacturer protocol.
To determine bacterial diversity, 16S ribosomal primers were designed from variable regions three and four and used for sequencing in an Illumina platform (Illumina, 2013) in the FISABIO foundation (Valencia, Spain). PCR reactions were conducted using total genomic DNA and Kappa polymerase (Kappa HiFi Hotstart Ready Mix, cat. No. KK2602, Kappa Biosystems). PCR products were cleaned up using magnetic beads (BeckmanCoulter Agencourt Ampure XP, cat. No. A63881) and resuspended in 10 mM Tris pH 8.5 or Resuspension buffer (RSB, Illumina Inc.).
Libraries were constructed using Nextera XT adapters (Illumina Inc.) using the same polymerase and cleaned again using magnetic beads.
Libraries were quantified by fluorescence using a Qubit DNA HS kit, library quality was analyzed using a Tapestation 2200 or Bioanalyzer 2100 and pooled to an equimolar ratio to be sequenced in a 2X300 bp MiSeq Run. The raw data was deposited in NCBI´s sra archive under bioproject accession number PRJNA591631.
Bioinformatic and Statistical Analysis
All the reads with a quality below 20 were filtrated by Prinseq-lite program (Schmieder and Edwards, 2011). R1 and R2 fastq files were joined by fastq-join from the ea-tools suit. Sequences were then checked for chimeras with Mothur (Schloss et al., 2009) using Greengenes 13-8 database (http://greengenes.lbl.gov/). The reads were clustered at 97% in OTUs (Operational Taxonomic Units) using uclust algorithm as implemented in the script pick_open_reference_otus.py from QIIME 1.9.1 (Caporaso et al., 2010; Lan et al., 2012). The taxonomy assignment was done by RDP method (Lan et al., 2012) against Greengenes 13-8 database. Each sample was further manually filtered to remove OTUs missing in at least four out of the five samples in the same treatment, with a frequency lower than 0.01% and with homology to the chloroplast, archaea, and mitochondrial sequences (Supporting Information Table S2). Alpha diversity estimations were calculated using the next scripts: multiple_rarefactions.py, alpha_diversity.py, collate_alpha.py and make_rarefaction_plots.py from QIIME 1.9.1. Rarefaction curves were generated by subsampling the OTU table with step increments of 10 and the sampling depth of 1,500 reads. The alpha diversity indexes Chao1, Shannon, and Simpson and the number of OTUs by Treatment, were plotted by “vegan” and “stats” for R v.3.3.1 (Oksanen et al., 2016; R Core Team, 2016). We applied the Pairwise Wilcox test tested to determine significant differences between the diversity indexes of the treatments.
The relative taxonomic abundance of the Samples and the Treatment was represented by collapsed histogram plotted by “RColorBrewer” and “ggplot2” libraries for R v.3.3.1 (Wickham, 2009; Neuwirth, 2014). The data table contained the relative frequency of the most frequent OTUs and the collapse frequency of the others by its closest common taxa. To visualize the relative abundance of the most frequent OTUs a heat map per treatment was constructed by “ape”, “gplots”, “vegan” and “RColorBrewer” libraries for R v.3.3.1 (Paradis et al., 2004; Warnes et al., 2016). The OTUs with a frequencies higher than 0.1% were represented and the soil samples were not considered.
To represent dissimilarities across the samples, we performed a nonmetric multidimensional scaling (NMDS) using the Bray-Curtis distances matrix. In the NMDS plot, the distance between objects represents the relative dissimilarity between them. The measure of the goodness of fit of the final NMDS plot was represented by the stress value, which indicates the match between inter-object distance and dissimilarity. The closer the stress value is to zero, the better the ordination (Manly, 1986; Quinn and Keough, 2002). Finally, the effect of crop and sample type was evaluated with a permutational multivariate ANOVA (PERMANOVA) using the Bray-Curtis dissimilarities (Anderson, 2001; Anderson and Walsh, 2013), the significance threshold for PERMANOVA was set at p < 0.001. All the analyses were performed using “vegan” and “mass” for R v.3.3.1 (Oksanen et al., 2016; Venables and Ripley, 2002).
We predict the metagenome functional profiling of the different samples base on the Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway by PICRUSt. To assess the accuracy of the PICRUSt predictions, the Nearest Sequenced Taxon Index was calculated for each sample and treatment (Langille et al., 2013; Tables S3 and S4). The effect of crop and sample type on KEGG orthology (KOs) was evaluated with a PERMANOVA using a Bray-Curtis dissimilarities matrix, which was previously calculated considering the relative abundance of functional categories in all samples. The significance threshold for PERMANOVA was set at p ≤ 0.001. When PERMANOVA was significant, differences between samples were determined with multiple pairwise comparisons using a Wilcoxon test with FDR correction set at p ≤ 0.05. All the analyses were performed using “vegan” and “mass” for R v.3.3.1.
We generated two cluster dendrograms based on the Bray-Curtis distance matrix estimating of the OTUs taxonomic abundance filtered at 0.01%, the KOs abundance predicted by PICRUSt by “vegan” for R v.3.3.1.
We made a graphical representation of the main drivers of the healthy and Corky/knot conditions using relative abundance data at order level.
Based on the functional dendrogram results we defined three main clusters, and we built Venn diagrams of the specific or shared KOs between corky clusters and the healthy cluster, and among coffee and tomato corky clusters and the soil cluster.
Results
Diversity and Species Richness of Bacterial Community
A total of 2,131,578 reads were sequenced by Illumina PE platform. The quality check removed 20% of the reads obtaining an average of 53,846 quality read filters per sample (max = 2,002,726, min = 1,426) (Table S1). The sequencing depth quality was confirmed by estimating the Chao, Shannon, and Simpson alpha diversity indexes and rarefaction curves of the number of observed OTUs, where all the curves reflected a saturated sampling (Figure S2).
The OTUs richness and diversity were calculated using Simpson and Shannon indexes. In general, in both crops the number of observed OTUs differed significantly among sample types (ANOVA, Tomato: F(7, 40) = 8.72, p = 5.77e-06; Coffee: F(7, 40) = 25.97, p = 1.73e-11), with higher values in soil samples than in root samples, and higher values in root samples than in nematode samples (Supporting Information Figure S3). Moreover, the coffee samples exhibited significantly higher values than tomato samples, with a maximum value of 310.45 ± 19.45 (observed OTUs) for coffee soil samples and a minimum value of 81.8 ± 51.9 (observed OTUs) for female wash solution from tomato samples (F(7, 40) = 4.43, p = 0.0394). Regarding the diversity indexes, significant differences were only detected across coffee samples (Shannon: F(7, 40) = 12.58, p = 1.25e-07, Simpson: F(7, 40) = 2.83, p = 0.206), registering the maximum diversity values in soil samples (7.147 ± 0.277) and the minimum in female wash solution samples (4.252 ± 0.609) (Supporting Information Figure S3).
The high-quality reads were clustered over 881 OTUs, most of them belonging to the phylum Proteobacteria (49%), Bacteroidetes (14%), and Actinobacteria (10%). The relative abundance of major genera and families by treatment is shown in Figure 1 (Supporting Information Figure S4). The genera Pseudomonas and Flavobacterium had the highest relative abundance with a total of 8.07% and 5.67%, respectively. The genera shared for all samples were Pseudomonas and Bacillus, while at OTUs level 560886 and 961783 (Pseudomonas) were the only ones common to all samples. The core microbiome of the tomato samples consisted of the genera Cellvibrio and Rhodoplanes, while Microbacterium, Achromobacter, Pedobacter, and Flavobacterium were the genera present in core microbiome of coffee. The OTUs 558264 (Achromobacter) and NCR.OTU928383 (Pedobacter) were shared by all the coffee samples (Figure 2). No common core microbiome was found among nematode samples. Root samples however, shared the genera Novosphingobium, Sphingobacterium, and Janthinobacterium. In the case of Janthinobacterium, we found 10 OTUs with a relative frequency lower than 0.1%. Specifically, the corky root had the genera Nocardia in common, although in low frequencies, 0.001% in corky root tomato and 0.009% in corky root coffee. Soil and Corky tomato samples shared 22 OTUs belonged to the genera Pseudomonas, Sphingobacterium, Flavobacterium, Rhodoplanes, Cellvibrio, Shinella, Luteolibacter, Brevundimonas, Janthinobacterium, Chitinophaga, Bacillus, and Agrobacterium, while Soil and Corky coffee samples shared 10 OTUs belonged to the genera Pedobacter, Janthinobacterium, and Candidatus Solibacter.
Figure 1 Relative taxonomic abundance by treatment represented by the most frequent OTUs, the least frequent taxa are collapsed in others OTUs.
Figure 2 Relative abundance heat map of the OTUs with a frequency higher than 0.1% by treatment. The soil samples are not taken into account. The dendrogram on the left shows the phylogenetic relationship between the OTUs, the dendrogram on the top illustrate the relationship between the samples based on Ray distance. Yellow, low relative abundance; orange, high relative abundance.
Permutational ANOVA (PERMANOVA) and NMDS analysis (NMDS) of the bacterial community revealed that both factors, crop species and sample type, clearly determine the bacterial community assemblage (Figure 3, Table 1). The microbiome structure of all coffee samples differed significantly from the all the tomato samples (PERMANOVA, F(1,80) = 213.987, R2 = 0.27180, p = 0.001; Figure 3). Moreover, the soil samples associated to each crop, differed from the other samples (Tomato: F(1,40) = 27.584, R2 = 0.85783, p = 0.001; Coffee: F(1,40) = 25.753, R2 = 0.84925, p = 0.001; Figure 3). A split analysis of all tomato or coffee samples showed significant differences between root and nematode samples (Figure S5) and among the different disease developmental stages (knots vs. corky tissues) and the nematode stages (eggs vs. females) (Figure S5).
Figure 3 Non-metric multidimensional scaling (NMDS) analysis of all the samples. NMDS plot for Bray–Curtis dissimilarities of the bacterial communities associated with tomato and coffee.
Table 1 Permutational multivariate ANOVA of the bacterial communities associated with different sample types of coffee and tomato.
The microbiome structure was modified by nematode infection, leading to the depletion of Pseudomonadales (Gammaproteobacteria) and Sphingobacteriales (Sphingobacteria) and the enrichment of Actinomycetales (Actinobacteria), in both crops. The microbiome structure variation of both crops differed in the taxonomic orders. The corky-root coffee samples exhibited an enrichment of Bacillales (Bacilli) and Burkholderiales (Betaproteobacteria) in comparison with healthy samples, while in corky-root tomato samples Saprospirales (Saprospirae), Chthoniobacterales (Spartobacteria), Alteromonadales (Gammaproteobacteria), and Xanthomonadales (Gammaproteobacteria) were enriched in comparison with healthy samples and Rhodobacterales were slightly depleted in corky-root samples. The pairwise comparisons of relative abundance between Healthy, Knot, and Corky coffee samples within the order Rhodobacterales were significant with a p-value of 0.03 (Figure S6).
Predictive Functional Analysis of the Bacterial Microbiome
The association of the microbial community with a metabolic profile was assessed by PICRUST analysis. We predicted the gene content of the bacterial microbiome and construct a dendrogram based on taxonomic and metabolic dissimilarities (Figure 4, Tables 1 and 2). The results showed that, while at a taxonomic level the samples grouped primarily based on the crop species (tomato vs. coffee plant) and within each of both groups, based on the origin of the sample (nematode stage vs. root source); at a functional level (KO), in both cases (tomato and coffee plants), the PBS wash solution eggs and the corky root samples grouped together, with the particularity that within the coffee cluster, knot samples were also present in this cluster and that in the tomato samples, egg samples belonged to this cluster. The healthy root samples of both crops clustered together along with knot tomato samples. Regarding coffee and tomato soil samples, at a taxonomic level they both clearly separated and clustered with their relative root and nematode crop samples, whereas at a functional level (KO), they clustered together and with the coffee infested roots and nematode stages samples.
Figure 4 Cluster dendrogram of dissimilarity matrix of (A) OTUs relative abundance and (B) KO's abundance predicted by PICRUSt.
Table 2 Permutational multivariate ANOVA of the metabolic functionality Kyoto Encyclopaedia of Genes and Genomes categories of the bacterial communities associated with different sample types of coffee and tomato crops.
The Healthy Root Cluster (healthy tomato roots, healthy coffee roots, and tomato root-knot) had 6,071 KOs, the Coffee Cork Cluster (corky roots, root-knots, wash solution egg) had 6,211 KOs and the Tomato Cork Cluster (corky roots, eggs, wash solution egg) had 5,264 KOs shared 3,147 KOs and had 135, 220, and 9 specific KOs, respectively. The metabolic biomarkers that defined the healthy root stage were defined by the presence of 135 KOs, which were classified principally in the following functional categories: Membrane transport (42 KOs), Carbohydrate metabolism (13 KOs), Cellular processes (eight KOs), Glycan biosynthesis and metabolism (eight KOs), Metabolism of cofactors and vitamins (seven KOs), and Poorly characterized functions (26 KOs) (Figure 5).
Figure 5 Venn diagram of the KO numbers specific to ot shared by the coffee and tomato corky root clusters and the healthy root cluster.
Moreover, the samples of the Healthy Cluster had 285 KOs with higher relative frequency and 135 KOs with a lower significant relative abundance compared to the Corky Coffee Cluster, while the Healthy Cluster samples had 113 KOs with significantly higher relative abundance and 50 KOs with a significantly lower relative abundance compared to the samples of the Corky Tomato Cluster (Table S5, Figure 6). The 285 mentioned KOs with higher relative frequency in the healthy cluster vs. in the corky coffee cluster classified principally in the following functional categories: Membrane transport (94 KOs), Amino acid metabolism (36 KOs), Transcription (36 KOs), Glycan biosynthesis and metabolism (26 KOs), Carbohydrate metabolism (25 KOs), Cellular processes and signaling (11 KOs), and Unclassified (61 KOs). The other 135 with lower relative frequency classified in almost the same functional categories: Membrane transport (105 KOs), Carbohydrate metabolism (61 KOs), Amino acid metabolism (50 KOs), Transcription (45 KOs), Cellular processes and signaling (13 KOs), and Unclassified (83 KOs). Only the metabolic functions Replication and repair and Glycan biosynthesis and metabolism showed higher numbers in the Healthy Cluster samples (15 and 26 KOs, respectively) than in the Corky Coffee Cluster samples (9 and 14 KOs, respectively). Regarding the 113 KOs with higher frequency in Healthy Cluster compared to Corky Tomato Cluster, they classified mainly in the functional categories: Membrane transport (20 KOs), Carbohydrate metabolism (20 KOs), Amino acid metabolism (19 KOs), Xenobiotics biodegradation and metabolism (19 KOs), Cellular processes and signaling (four KOs), and Unclassified (19 KOs), while the other 50 KOs with lower frequency classified on Metabolism of cofactors and vitamins (six KOs), Amino acid metabolism (five KOs) Cellular processes and signaling (six KOs), and Unclassified (six KOs). While the metabolic profile of the Corky Coffee Cluster samples was defined by the increase in the KOs number of the different metabolic categories compared to the Healthy Cluster samples, the Corky Tomato Cluster samples were defined by a decrease of the same functional categories.
Figure 6 Venn diagram of the KO numbers (A) with higher significant relative frequency, and (B) with lower significant relative frequency, on the healthy root cluster relative to coffee corky root and tomato corky root clusters.
Corky stage has been characterized by the presence of 10 specific KOs, which were classified principally in the following functional categories: Cellular processes and signaling (KOs 4), Membrane transport (3), Unclassified (3). Compared to Corky Coffee stages, only 220 specific KOs were found and classified in: Membrane transport (34 KOs), Carbohydrate metabolism (23 KOs), Amino acid metabolism (15 KOs), Cellular processes and signaling (14 KOs), Enzyme families (13 KOs), Xenobiotics biodegradation and metabolism (seven KOs), and Unclassified (57 KOs). In Corky Tomato stages, only nine specific KOs were found, involved principally in Membrane transport (three KOs), Cellular processes (signaling, 2) and Unclassified (four KOs).
A specific analysis of the metabolic profile shared by the Soil Clusters and both Cork Cluster, showed that Cork coffee Cluster and Soil Cluster shared 290 KOs, with 14.14% (67 KOs) involved in Environmental Information Processing and 60.34% (286 KOs) involved in Metabolism. The Cluster Cork Tomato and Soil Cluster shared 128KO's which were involved in Environmental Information Processing 35,96% (64 KOs) and 35,96% (64 KOs) involved on Metabolism (Figure S7).
Discussion
In this study, we used 16S rDNA gene Illumina sequencing to analyze the microbial communities associated with M. enterolobii and M. paranaensis and different disease stages of the MDC in tomato and coffee field plants. The microbiome isolated from soil, healthy roots, infested roots, and different nematode developmental stages of each crop showed a very complex and distinct structure.
The two crop systems (tomato vs. coffee) differed in their bacterial community, alpha diversity and OTUs richness, all these variables were higher in coffee than in tomato samples. Similarly, the crop type clearly affected the composition of MDC bacterial microbiome, as revealed by NMDS and PERMANOVA. Coffee samples came from perennial plantations, grown under agroforestry systems, holding several tree shade species and with minor use of fertilizers and pesticides (Gordon et al., 2007; González-Zamora et al., 2016). Moreover, tillage is not employed, leading to the establishment of many soil microbiota species (Legrand et al., 2018; Dube et al., 2019). By contrast, tomato, an annual crop with short cycle, is known by its intensive cropping systems, with high demand of agrochemicals (both fertilizers and pesticides), applied before and during the growing cycle (Pérez-Olvera et al., 2011; Pratt and Ortega, 2019). Previous studies have shown that plant management or status has a clear effect on the associated microbiome (Coleman-Derr et al., 2016; Hartman et al., 2018; Compant et al., 2019) and that intensively grown crops have a lower diversity of associated plant and soil microbiota than agroecosystems subjected to less aggressive agricultural practices (Rathore et al., 2019; Saleem et al., 2019).
We consider that the significant differences of microbiota diversity found between the two studied crop systems could be related with the management practices, leading to the establishment of a more diverse bacterial community in the coffee agroecosystem.
Moreover, we found significantly higher alpha diversity and OTUs richness in root samples than in nematode samples. However, for both crops, there was no significant differences between disease stages or nematode life cycle stages. These results differ with those observed by Wolfgang et al. (2019), who describe decreasing alpha diversity during the root disease process. Conversely, we observed a non-significant increasing diversity in disease stages in both crops.
We identified Pseudomonadales, Sphingobacteriales, and Actinomycetales as the main common drivers of the MDC implicated in the infection process for both crops. Additionally, we observed crop-specific drivers of the MDC infection process: Bacillales and Burkholderiales in coffee samples and Saprospirales, Chthoniobacterales, Alteromonadales, and Xanthomonadales in tomato samples. Few similarities with the results published by Tian et al. (2015) and Wolfgang et al. (2019) were observed. Our results and the two cited studies only agree in the richness reduction of Pseudomonadales in the corky-root samples in comparison with the healthy ones. The order Pseudomonadales contains species listed as opportunistic parasites (Castillo et al., 2017); however, other species in the same order have shown nematicide effects by the production of compounds such as cyanide (Siddiqui et al., 2006). The discrepancies regarding the main drivers observed in these studies may be due to cultivar differences (Rybakova et al., 2017). These differences could be explained by the diverse sets of available microorganisms with varied abilities to attach to the nematode cuticle. We also consider that the rhizobiome composition could be affected by the physicochemical conditions of the soil (Bonito et al., 2019; Saad et al., 2019). Wolfgang et al. (2019) proposed that changes in physiological conditions of plant cells is the cause by a microbial structure shift, and this is influenced by the ability of some bacteria to adhere to the nematode cuticle (Elhady et al., 2017). In addition, it has been shown that the rhizosphere microbiome may influence the severity of RKN infection by promoting pathogens in infested soils and suppressing pathogens in uninfected soils (Zhou et al., 2019).
At the genera level, only Nocardia was shared by the corky-root samples. Species of this genus are found worldwide in soils rich in organic matter, some of which are plant pathogenic (Xing et al., 2011; Bai et al., 2016; Conville et al., 2017). Nocardia sp. is the causal agent of the disease known as tobacco false broomrape, which causes cracked tumors in the main and secondary roots, resembling the hyperplasia produced in the advanced stages of the MDC in tomato and coffee (Morán et al., 2013; Morán-Gómez et al., 2018). Moreover, Agrobacterium was externally recovered from M. enterolobii female and tomato corky-roots. Species of this genus have virulence genes that affect the response of the plant cells to growth regulators (auxins and cytokinins) and induce uncontrolled cell division, which results in tissue proliferation and gall formation (Bourras et al., 2015; Lacroix and Citovsky, 2019). Thus, the genera Nocardia and Agrobacterium, associated with RKN microbiome reported in this and other studies (Cao et al., 2015), may contribute to the formation of corky root symptoms in the MDC.
Additionally, M. enterolobii produces the MeTCTP (translationally controlled tumor protein) effector that promotes parasitism by suppressing programmed cell death in host plants (Zhuo et al., 2017). The combined effect of the tumor caused by the nematode and the bacteria could explain the remarkably wide range of hosts and severity of symptoms observed in plants infested with this nematode (Elling, 2013; Jones et al., 2013); however, this requires further investigation.
The grouping pattern based on the microbiome metabolic prediction differs from the grouping pattern based on the taxonomic information. At taxonomic level we observed two big clusters based on the crop of origin (Tomato vs. Coffee) and the type of sample (Roots vs. Nematodes), leaving the soil samples alone as an external group. By contrast, we observed three metabolic clusters: the Healthy Cluster (Healthy coffee and tomato roots plus Tomato Root-Knot samples); the Coffee Corky cluster (Coffee Corky-Roots and Root-Knots plus Coffee PBS washout nematode Eggs samples); and the Tomato Corky cluster (Tomato Corky-Roots plus Tomato nematode Eggs and PBS washout Eggs samples).
Metabolic functions of the microbiota associated with the soil samples of both crops were comparable, as well as the metabolic profile of healthy roots of both crops. The metabolic functional categories enriched in Corky samples versus soil samples, considering the KOs analysis, were associated to pathogenicity, and include cell signaling and regulation and membrane transport, and other categories like carbohydrate metabolism, amino acid metabolism, cofactor and vitamins metabolism, family enzymes, and transcription factors related with the increase in the metabolic activity in the nematode feeding site. Our metagenomics analyses conducted in the two phathosystems (Tomato- M. enterolobii, and Coffee- M. paranaensis) suggest that root infection by the nematode is, on the contrary, not linked to a specific microbiome community, but is associated with a metabolic profile shift characterized by an increase in pathogenicity and central metabolism functional categories, which is consistent with Tian et al. (2015) observations.
Additionally, one of the defense mechanisms of the coffee plant toward the attack of Meloidogyne spp. is the accumulation of phenolic compounds at the feeding site to limit their reproduction (Albuquerque et al., 2010; Silva et al., 2010). Our results, indicate that corky-roots of coffee plants had more metabolic profiles that encode enzymes involved in the pathway of benzoate degradation (aminobenzoate and benzoate) than healthy coffee roots. Benzoate is an intermediate in the pathway of anaerobic metabolism of aromatic compounds, i.e. phenolic compounds (Carmona et al., 2009). It is possible that some of the bacteria associated with coffee corky-root prevent the accumulation of phenolic compounds at the feeding site of M. paranaensis and contribute to the infection process, as occurs in other pathosystems (Cheng et al., 2013).
Our results indicate that, during nematode pathogenesis and nematode feeding site induction in plants, the MDC infection process could be mediated by particular metabolic microbial categories at a functional level, even though the two MDCs studied here involved different RKN species and unspecific associated bacterial taxa. Nonetheless, further studies are required to gain a deeper understanding of the bacterial functional roles in the MDC infection process.
Data Availability Statement
The datasets generated for this study can be found in NCBI using accession number PRJNA591631.
Author Contributions
GC, DL-L, ArL designed the experiments. DL-L and LV collected and prepared the samples. DL-L, ArL extracted the DNA. AM and AmL performed the sequencing. ArL and AA performed the bioinformatic analysis. DD performed the statistical analyses. GC, ArL, DD, DL-L, LV, and AA-S wrote the manuscript. GC, ArL, DD, DL-L, LV, AA-S, AM, AmL read, revised, and approved the manuscript.
Funding
The research leading to these results has received funding from the Research Projects of High Strategic Value for Society (Project 2003530920) founded by Institute of Ecology A. C (INECOL).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank to Martín Aluja Schuneman-Hofer for his support in the management of the resources that financed this project. We would to thanks to Luisa Cruz for the English writing revision.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2020.00136/full#supplementary-material
References
Albuquerque, E. V. S., Carneiro, R. M. D. G., Costa, P. M., Gomes, A. C. M. M., Santos, M., Pereira, A. A., et al. (2010). Resistance to Meloidogyne incognita expresses a hypersensitive-like response in Coffea arabica. Eur. J. Plant Pathol. 127, 365–373. doi: 10.1007/s10658-010-9603-3
Anderson, M. J., Walsh, D. C. I. (2013). PERMANOVA, ANOSIM, and the Mantel test in the face of heterogeneous dispersions: What null hypothesis are you testing? Ecol. Monogr. 83, 557–574. doi: 10.1890/12-2010.1
Anderson, M. J. (2001). A new method for non-parametric multivariate analysis of variance. Austral Ecol. 26, 32–46. doi: 10.1111/j.1442-9993.2001.01070.pp.x
Bai, J. L., Wang, Y., Qin, S., Ding, P., Xing, K., Yuan, ,. B., et al. (2016). ). Nocardia jiangsuensis sp. nov., an actinomycete isolated from coastal soil. Int. J. Syst. Evol. Micr. 66, 4633– 4638. doi: 10.1099/ijsem.0.001402
Baquiran, J.-P., Thater, B., Sedky, S., De Ley, P., Crowley, D., Orwin, P. M. (2013). Culture-independent investigation of the microbiome associated with the nematode Acrobeloides maximus. PloS One 8, e67425. doi: 10.1371/journal.pone.0067425
Bertrand, B., Nuñez, C., Sarah, J. L. (2000). Disease complex in coffee involving Meloidogyne arabicida and Fusarium oxysporum. Plant Pathol. 49, 383–388. doi: 10.1046/j.1365-3059.2000.00456.x
Bonito, G., Benucci, G. M. N., Hameed, K., Weighill, D., Jones, P., Chen, K.-H., et al. (2019). Fungal-bacterial networks in the Populus rhizobiome are impacted by soil properties and host genotype. Front. Microbiol. 10, 481. doi: 10.3389/fmicb.2019.00481
Bourras, S., Rouxel, T., Meyer, M. (2015). Agrobacterium tumefaciens gene transfer: how a plant pathogen hacks the nuclei of plant and nonplant organisms. Phytopathology 105, 1288–1301. doi: 10.1094/PHYTO-12-14-0380-RVW
Brader, G., Compant, S., Vescio, K., Mitter, B., Trognitz, F., Ma, L.-J., et al. (2017). Ecology and genomic insights into plant-pathogenic and plant-nonpathogenic endophytes. Annu. Rev. Phytopathol. 55, 61–83. doi: 10.1146/annurev-phyto-080516-035641
Bulgarelli, D., Rott, M., Schlaeppi, K., Ver Loren van Themaat, E., Ahmadinejad, N., Assenza, F., et al. (2012). Revealing structure and assembly cues for Arabidopsis root-inhabiting bacterial microbiota. Nature 488, 91–95. doi: 10.1038/nature11336
Busby, P. E., Soman, C., Wagner, M. R., Friesen, M. L., Kremer, J., Bennett, A., et al. (2017). Research priorities for harnessing plant microbiomes in sustainable agriculture. PloS Biol. 15, e2001793. doi: 10.1371/journal.pbio.2001793
Cao, Y., Tian, B., Ji, X., Shang, S., Lu, C., Zhang, K. (2015). Associated bacteria of different life stages of Meloidogyne incognita using pyrosequencing-based analysis. J. Basic Microb. 55, 950–960. doi: 10.1002/jobm.201400816
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
Carmona, M., Zamarro, M. T., Blazquez, ,. B., Durante-Rodríguez, G., Juarez, J. F., Valderrama, J. A., et al. (2009). Anaerobic catabolism of aromatic compounds: a genetic and genomic view. Microbiol. Mol. Biol. R. 74, 71–133. doi: 10.1128/MMBR.00021-08
Carneiro, R. M. D. G., Carneiro, R. G., Abrantes, I. M. O., Santos, M. S. N. A., Almeida, M. R. A. (1996). Meloidogyne paranaensis n.sp. (Nemata: Meloidogynidae), a root knot nematode parasitizing coffee in Brazil. J. Nematol. 28, 177–189.
Carneiro, R., Tigano, M., Randig, O., Almeida, M. R., Sarah, J. L. (2004). Identification and genetic diversity of Meloidogyne spp. (Tylenchida: Meloidogynidae) on coffee from Brazil, Central America and Hawaii. Nematology 6, 287–298. doi: 10.1163/1568541041217942
Castillo, J. D., Vivanco, J. M., Manter, D. K. (2017). Bacterial Microbiome and Nematode Occurrence in Different Potato Agricultural Soils. Microb. Ecol. 74, 888–900. doi: 10.1007/s00248-017-0990-2
Cheng, X. Y., Tian, X. L., Wang, Y. S., Lin, R. M., Mao, Z. C., Chen, N., et al. (2013). Metagenomic analysis of the pinewood nematode microbiome reveals a symbiotic relationship critical for xenobiotics degradation. Sci. Rep. 3, 1869. doi: 10.1038/srep01869
Coleman-Derr, D., Desgarennes, D., Fonseca-Garcia, C., Gross, S., Clingenpeel, S., Woyke, T., et al. (2016). Biogeography and cultivation affect microbiome composition in drought-adapted plants. New Phythol. 209, 798–811. doi: 10.1111/nph.13697
Compant, S., Samad, A., Faist, H., Sessitsch, A. (2019). A review on the plant microbiome: Ecology, functions, and emerging trends in microbial application. J. Adv. Res. 19, 29–37. doi: 10.1016/j.jare.2019.03.004
Conville, P. S., Brown-Elliott, B. A., Smith, T., Zelazny, A. M. (2017). The complexities of Nocardia taxonomy and identification. J. Clin. Microbiol. 56, e01419–e01417. doi: 10.1128/JCM.01419-17
Coyne, D. L., Ross, J. L. (2014). Protocol for Nematode Resistance Screening: Root Knot Nematodes, Meloidogyne spp (Ibadan, Nigeria: International Institute of Tropical Agriculture (IITA)), 27 pp.
Dube, J. P., Valverde, A., Steyn, J. M., Cowan, D. A., van der Waals, J. (2019). Differences in bacterial diversity, composition and function due to long-term agriculture in soils in the eastern free State of South Africa. Diversity 11, 61. doi: 10.3390/d11040061
Elhady, A., Giné, A., Topalovic, O., Jacquiod, S., Sørensen, S. J., Sorribas, F. J., et al. (2017). Microbiomes associated with infective stages of root-knot and lesion nematodes in soil. PloS One 12, e0177145. doi: 10.1371/journal.pone.0177145
Elling, A. A. (2013). Major emerging problems with minor Meloidogyne species. Phytopathology 103, 1092–1102. doi: 10.1094/PHYTO-01-13-0019-RVW
Finkel, O. M., Castrillo, G., Herrera Paredes, S., Salas González, I., Dangl, J. L. (2017). Understanding and exploiting plant beneficial microbes. Curr. Opin. Plant Biol. 38, 155–163. doi: 10.1016/j.pbi.2017.04.018
Ganaie, M. A., Khan, T. A. (2011). Studies on the interactive effect of Meloidogyne incognita and Fusarium solani on Lycopersicon esculentum, Mill. Int. J. Bot. 7, 205–208. doi: 10.3923/ijb.2011.205.208
González-Zamora, A., Esperón-Rodríguez, M., Barradas, V. (2016). Mountain cloud forest and grown-shade coffee plantations: A comparison of tree biodiversity in central Veracruz, Mexico. For. Syst. 25 (1), e055. doi: 10.5424/fs/2016251-07538
Gordon, C., Manson, R., Sundberg, J., Cruz-Angón, A. (2007). Biodiversity, profitability, and vegetation structure in a Mexican coffee agroecosystem. Agric. Ecosyst. Environ. 118, 256–266. doi: 10.1016/j.agee.2006.05.023
Hajji-Hedfi, L., Regaieg, H., Larayedh, A., Chihani, N., Horrigue- Raouani, N. (2018). Biological control of wilt disease complex on tomato crop caused by Meloidogyne javanica and Fusarium oxysporum f.sp. lycopersici by Verticillium leptobactrum. Environ. Sci. Pollut. R. 25, 18297–18302. doi: 10.1007/s11356-017-0233-6
Hartman, K., van der Heijden, M. G. A., Wittwer, R. A., Banerjee, S., Walser, J. C., Schlaeppi, K. (2018). Cropping practices manipulate abundance patterns of root and soil microbiome members paving the way to smart farming. Microbiome 6, 14. doi: 10.1186/s40168-017-0389-9
Illumina (2013). 16S Metagenomic Sequencing Library Preparation. Illumina.com, (B), pp. 1–28. Available at: http://support.illumina.com/content/dam/illumina-support/documents/documentation/chemistry_documentation/16s/16s-metagenomic-library-prep-guide-15044223-b.pdf.
Jakuschkin, B., Fievet, V., Schwaller, L., Fort, T., Robin, C., Vacher, C. (2016). Deciphering the pathobiome: Intra- and interkingdom interactions involving the pathogen Erysiphe alphitoides. Microb. Ecol. 72, 870–880. doi: 10.1007/s00248-016-0777-x
Jones, J. T., Haegeman, A., Danchin, E. G., Gaur, H. S., Helder, J., Jones, M. G., et al. (2013). Top 10 plant-parasitic nematodes in molecular plant pathology. Mol. Plant Pathol. 14, 946– 961. doi: 10.1111/mpp.12057
Kiers, E. T., van der Heijden, M. G. A. (2006). Mutualistic stability in the arbuscular mycorrhizal symbiosis: exploring hypotheses of evolutionary cooperation. Ecology 87, 1627–1636. doi: 10.1890/0012-9658(2006)87[1627:msitam]2.0.co;2
Kiers, E. T., Rousseau, R. A., West, S. A., Denison, R. F. (2003). Host sanctions and the legume–rhizobium mutualism. Nature 425, 78–81. doi: 10.1038/nature01931
Kiewnick, S., Karssen, G., Brito, J. A., Oggenfuss, M., Frey, J. E. (2008). First report of root-knot nematode Meloidogyne enterolobii on tomato and cucumber. Plant Dis. 92, 1370. doi: 10.1094/PDIS-92-9-1370A
Koike, S. T., Gladders, P., Paulus, A. O. (2006). Vegetable diseases: a color handbook (San Diego: Academic Press), 448p.
Lacroix, B., Citovsky, V. (2019). Pathways of DNA transfer to plants from Agrobacterium tumefaciens and related bacterial species. Annu. Rev. Phytopathol. 57, 231–251. doi: 10.1146/annurev-phyto-082718-100101
Lamichhane, J. R., Venturi, V. (2015). Synergisms between microbial pathogens in plant disease complexes: a growing trend. Front. Plant Sci. 6, 385. doi: 10.3389/fpls.2015.00385
Lan, Y., Wang, Q., Cole, J. R., Rosen, G. L. (2012). Using the RDP classifier to predict taxonomic novelty and reduce the search space for finding novel organisms. PloS One 7, e32491. doi: 10.1371/journal.pone.0032491
Langille, M. G. I., Zaneveld, J., Caporaso, J. G., McDonald, D., Knights, D., Reyes, J., et al. (2013). Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat. Biotechnol. 31, 814–823. doi: 10.1038/nbt.2676
Legrand, F., Picot, A., Cobo-Díaz, J. F., Carof, M., Chen, W., Le Floch, G. (2018). Effect of tillage and static abiotic soil properties on microbial diversity. Appl. Soil Ecol. 132, 135–145. doi: 10.1016/j.apsoil.2018.08.016
López-Lima, D., Carrión, G., Sánchez-Nava, P., Desgarennes, D., Villain, L. (in press). Fungal diversity and Fusarium oxysporum pathogenicity associated with coffee corky-root disease in Mexico. Rev. FCA UNCUYO. Available at: http://revista.fca.uncu.edu.ar/index.php?option=com_content&view=article&id=614:fungal-diversity-and-fusarium-oxysporum-pathogenicity-associated-with-coffee-corky-root-disease-in-mexico&catid=20:en-prensa&Itemid=23%20.
Lopez-Lima, D., Sánchez-Nava, P., Carrión, G., Espinosa De Los Monteros, A., Villain, L. (2015). Corky-root symptoms for coffee in central Veracruz are linked to the root-knot nematode Meloidogyne paranaensis, a new report for Mexico. Eur. J. Plant Pathol. 141, 623–629. doi: 10.1007/s10658-014-0564-9
Manly, B. F. (1986). Multivariate Statistical Methods: A Primer (London: Chapman and Hall, Ltd), 159.
Martínez-Gallardo, J. A., Díaz-Valdés, T., Allende-Molar, R., García-Estrada, R. S., Carrillo-Fasio, J. A. (2015). Primer reporte de Meloidogyne enterolobii parasitando tomate en Culiacán, Sinaloa, México. REMEXCA 11, 2165–2168. doi: 10.29312/remexca.v0i11.786
Martínez-Gallardo, J. A., Díaz-Valdes, T., Allende-Molar, R., Retes-Manjarrez, J. E., Carrillo-Fasio, J. A. (2019). Identification and distribution of Meloidogyne spp. In tomato in Sinaloa México. REMEXCA 10, 453–459. doi: 10.29312/remexca.v10i2.392
Monteiro, J. M. S., Cares, J. E., Gomes, A. C. M. M., Correa, V. R., Mattos, V. S., Santos, M. F. A., et al. (2016). First report of, and additional information on, Meloidogyne konaensis (Nematoda: Meloidogyninae) Parasitizing various crops in Brazil. Nematology 18, 831–844. doi: 10.1163/15685411-00002997
Morán, Y. M., Chacón, O., Córdoba-Selles, M. C., Domínguez, R., Herrera Isla, L., Borrás-Hidalgo, O. (2013). Identification and molecular characterization of Nocardia sp. as a new causal agent of tobacco false broomrape. J. Phytopathol. 161, 86–91. doi: 10.1111/jph.12029
Morán-Gómez, Y. M., Trémols-González, A. J., Domínguez-Larrinaga, R., Carrillo-Benites, ,. M. G., Cabrera-Alfonso, J. R. (2018). Assessment of antimicrobial potential of iron on phytopathogenic isolates of Nocardia sp. and their effect on tobacco false broomrape symptom expression. Trop. Plant Pathol. 43, 333–340. doi: 10.1007/s40858-018-0213-5
Mwangi, M. W., Muiru, W. M., Narla, R. D., Kimenju, J. W., Kariuki, G. M. (2019). Management of Fusarium oxysporum f. sp. Lycopersici and root-knot nematode disease complex in tomato by use of antagonistic fungi, plant resistance and neem. Biocontrol Sci. Techn. 29, 229–238. doi: 10.1080/09583157.2018.1545219
Neuwirth, E. (2014). RcolorBrewer: ColorBrewer Palettes. R package version 1.1-2. Available at : https://CRAN.R-project.org/package=RColorBrewer.
Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., Minchin, P. R., O'Hara, R. B., et al. (2016). Vegan: Community Ecology Package: Ordination, Diversity and Dissimilarities. Package. R package.
Pérez-Olvera, M. A., Navarro-Garza, H., Miranda-Cruz, E. (2011). “Use of pesticides for vegetable crops in Mexico,” in Pesticides in the modern world – Pesticides use and management. Ed. Stoytcheva, M. (Croacia: InTech Rijeka), 97–118. doi: 10.5772/18510
Paradis, E., Claude, J., Strimmer, K. (2004). APE: analyses of phylogenetics and evolution in R language. Bioinformatics 20, 289–290. doi: 10.1093/bioinformatics/btg412
Pratt, L., Ortega, J. M. (2019). Protected agriculture in Mexico: building the methodology for the first certified agricultural green bond. Inter-Ameriacn Developmental Bank Technical Note No. 1668. doi: 10.18235/0001705
Quinn, G. P., Keough, M. J. (2002). Experimental Design and Data Analysis for Biologists (New York, NY: Cambridge University Press), 537.
R Core Team (2016). R: A language and environment for statistical computing (Vienna, Austriag/: R Foundation for Statistical Computing). URL https://www.R-project.org.
Randig, O., Bongiovanni, M., Carneiro, R., Castagnone-Sereno, P. (2002). Genetic diversity of root-knot nematodes from Brazil and development of SCAR markers specific for the coffee-damaging species. Genome 45, 862–870. doi: 10.1139/g02-054
Rathore, R., Germaine, K. J., Forristal, P. D., Spink, J., Dowling, D. N. (2019). “Meta-omics approach to unravel the endophytic bacterial communities of Brassica napus and other,” in Endophytes for a growing world. Eds. Hodkinson, T. R., Doohan, F. M., Saunders, M. J., Murphy, B. R. (New York, Ny: Cambridge University Press). doi: 10.1017/9781108607667
Rizvi, R., Ansari, R. A., Iqbal, A., Ansari, S., Sumbul, A., Mahmood, I., et al. (2015). Dynamic role of organic matter and bioagent for the management of Meloidogyne incognita–Rhizoctonia solani disease complex on tomato in relation to some growth attributes. Cogent Food Agric. 1, 1. doi: 10.1080/23311932.2015.1068523
Rosa, J. M. O., Westerich, J. N., Wilcken, S. R. S. (2014). Reaction of tomato genotypes and hybrids to Meloidogyne enterolobii. Cienc. Rural 44, 1166–1171. doi: 10.1590/0103-8478cr20130041
Ryan, E. T. (2013). The intestinal pathobiome: its reality and consequences among infants and young children in resource-limited settings. J. Infect. Dis. 208, 1732–1733. doi: 10.1093/infdis/jit509
Rybakova, D., Mancinelli, R., Wikström, M., Birch-Jensen, A. S., Postma, J., Ehlers, R. U., et al. (2017). The structure of the Brassica napus seed microbiome is cultivar-dependent and affects the interactions of symbionts and pathogens. Microbiome 5, 104. doi: 10.1186/s40168-017-0310-316
Saad, R. F., Kobaissi, A., Echevarria, G., Kidd, P., Calusinska, M., Goux, X., et al. (2019). Influence of new agromining cropping systems on soil bacterial diversity and the physico-chemical characteristics of an ultramafic soil. Sci. Total Environ. 645, 380–392. doi: 10.1016/j.scitotenv.2018.07.106
Saleem, M., Hu, J., Jousset, A. (2019). More than the sum of its parts: microbiome biodiversity as a driver of plant growth and soil health. Annu. Rev. Ecol. Evol. Sys. 50, 6.1–6.24. doi: 10.1146/annurev-ecolsys-110617-062606
Sasanelli, N., Ciccarese, F., Papajova, I. (2008). Aphanocladium album by via sub-irrigation in the control of Pyrenochaeta lycopersici and Meloidogyne incognita on tomato in a plastic-house. Helminthologia 45, 137–142. doi: 10.2478/s11687-008-0027-y
Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/AEM.01541-09
Schmieder, R., Edwards, R. (2011). Quality control and preprocessing of metagenomic datasets. Bioinformatics 27, 863–864. doi: 10.1093/bioinformatics/btr026
Siddiqui, I. A., Shaukat, S. S., Sheikh, I. H., Khan, A. (2006). Role of cyanide production by Pseudomonas fluorescens CHA0 in the suppression of root-knot nematode, Meloidogyne javanica in tomato. World J. Microbiol. Biotechnol. 22, 641–650. doi: 10.1007/s11274-005-9084-2
Silva, R. V., Oliveira, R. D. L., Nascimento, K. J. T., Rodrogues, F. A. (2010). Biochemical responses of coffee resistance against Meloidogyne exigua mediated by silicon. Plant Pathol. 59 (3), 586–593. doi: 10.1111/j.1365-3059.2009.02228.x
Taher, I. Y., Kassab, A. S., Mahgoob, A. E., Taha, E. H., Banora, M. (2017). Disease complex in tomato cv. 035 involving Meloidogyne incognita and Fusarium oxysporum f. sp. lycopersici. Arab. Univ. J. Agric. Sci. 25, 349–352. doi: 10.21608/ajs.2017.13617
Tian, X., Cheng, X., Mao, Z., Chen, G., Yang, J., Xie, B. (2011). Composition of bacterial communities associated with a plant parasitic nematode Bursaphelenchus mucronatus. Curr. Biol. 62, 117–125. doi: 10.1007/s00284-010-9681-7
Tian, B. Y., Cao, Y., Zhang, K. Q. (2015). Metagenomic insights into communities, functions of endophytes, and their associates with infection by root-knot nematode, Meloidogyne incognita, in tomato roots. Nature 5, 17087. doi: 10.1038/srep17087
Tigano, M., de Sequeira, K., Castagnone-Sereno, P., Mulet, K., Queiroz, P., dos Santos, M., et al. (2010). Genetic diversity of the root-knot nematode Meloidogyne enterolobii and development of a SCAR marker for this guava-damaging species. Plant Pathol. 59, 1054–1061. doi: 10.1111/j.1365-3059.2010.02350.x
Toju, H., Tanaka, Y. (2019). Consortia of anti-nematode fungi and bacteria in the rhizosphere of soybean plants attacked by root-knot nematodes. R. Soc Open Sci. 6, 181693. doi: 10.1098/rsos.181693
Toju, H., Peay, K. G., Yamamichi, M., Narisawa, K., Hiruma, K., Naito, K., et al. (2018). Core microbiomes for sustainable agroecosystems. Nat. Plants 4, 247–257. doi: 10.1038/s41477-018-0139-4
Tollenaere, C., Susi, H., Laine, A. L. (2016). Evolutionary and epidemiological implications of multiple infection in plants. Trends Plant Sci. 21, 80–90. doi: 10.1016/j.tplants.2015.10.014
van Bruggen, A. H. C., Narouei-Khandan, H. A., Gravel, V., Blok, J. W. (2016). Corky root severity, root knot nematode galling and microbial communities in soil, rhizosphere and rhizoplane in organic and conventional greenhouse compartments. Appl. Soil Ecol. 100, 112–123. doi: 10.1016/j.apsoil.2015.11.015
Vandenkoornhuyse, P., Quaiser, A., Duhamel, M., Le Van, A., Dufresne, A. (2015). The importance of the microbiome of the plant holobiont. New Phythol. 206, 1196–1206. doi: 10.1111/nph.13312
Vannier, N., Agler, M., Hacquard, S. (2019). Microbiota-mediated disease resistance in plants. PloS Pathog. 15, e1007740. doi: 10.1371/journal.ppat.1007740
Vayssier-Taussat, M., Albina, E., Citti, C., Cosson, J. F., Jacques, M.-A., Lebrun, M.-H., et al. (2014). Shifting the paradigm from pathogens to pathobiome: new concepts in the light of meta-omics. Front. Cell. Infect. Mi. 4, 29. doi: 10.3389/fcimb.2014.00029
Venables, W. N., Ripley, B. D. (2002). Modern Applied Statistics with S (New York, NY: Springer), 497.
Villain, L., Sarah, J. L., Hernández, A., Bertrand, B., Anthony, F., Lashermes, P., et al. (2013). Diversity of root-knot nematodes parasitizing coffee in Central America. Nematropica 43, 194–206.
Villain, L., Lima-Salgada, S. M., Trinh-Phap, Q. (2018). “Plant parasitic nematodes in subtropical and tropical agriculture,” in Nematode parasites of coffee and cocoa. Eds. Sikora, R. A., Coyne, D., Hallmann, J., Timper, P. (Wallingford: CABI), 536–583.
Warnes, G. R., Bolker, B., Bonebakker, L., Gentleman, R., Liaw, W. H. A., Lumley, T., et al. (2016). gplots: Various R Programming Tools for Plotting Data. R package version 3.0.1. Available at: https://CRAN.R-project.org/package=gplots.
Wickham, H. (2009). ggplot2: Elegant Graphics for Data Analysis, (New York, Ny: Springer Verlag), 213.
Wolfgang, A., Taffner, J., Guimarães, R. A., Coyne, D., Berg, G. (2019). Novel strategies for soil-borne diseases: exploiting the microbiome and volatile-based mechanisms toward controlling Meloidogyne-based disease complexes. Front. Microbiol. 10, 1296. doi: 10.3389/fmicb.2019.01296
Xing, K., Qin, S., Fei, S. M., Lin, Q., Bian, G. K., Miao, Q., et al. (2011). Nocardia endophytica sp. nov., an endophytic actinomycete isolated from the oil-seed plant Jatropha curcas L. Int. J. Syst. Evol. Microbiol. 61, 1854–1858. doi: 10.1099/ijs.0.027391-0
Zhou, D., Feng, H., Schuelke, T., De Santiago, A., Zhang, Q., Zhang, J., et al. (2019). Rhizosphere microbiomes from root knot nematode non-infested plants suppress nematode infection. Microbial Ecol. 78, 470–481. doi: 10.1007/s00248-019-01319-5
Keywords: corky root, pathobiome, Meloidogyne enterolobii, Meloidogyne paranaensis, functional profile
Citation: Lamelas A, Desgarennes D, López-Lima D, Villain L, Alonso-Sánchez A, Artacho A, Latorre A, Moya A and Carrión G (2020) The Bacterial Microbiome of Meloidogyne-Based Disease Complex in Coffee and Tomato. Front. Plant Sci. 11:136. doi: 10.3389/fpls.2020.00136
Received: 07 November 2019; Accepted: 29 January 2020;
Published: 27 February 2020.
Edited by:
Amir Szitenberg, Dead Sea and Arava Science Center, IsraelReviewed by:
Tariq Mukhtar, Pir Mehr Ali Shah Arid Agriculture University, PakistanXingzhong Liu, Chinese Academy of Sciences, China
Copyright © 2020 Lamelas, Desgarennes, López-Lima, Villain, Alonso-Sánchez, Artacho, Latorre, Moya and Carrión. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Gloria Carrión, Z2xvcmlhLmNhcnJpb25AaW5lY29sLm14
†These authors have contributed equally to this work