Zelalem M. Taye1*
Bobbi L. Helgason2
Jennifer K. Bell2
Charlotte E. Norris2Sally Vail3
Stephen J. Robinson3
Isobel A. P. Parkin3
Melissa Arcand2
Steven Mamet2
Matthew G. Links4,5Tanner Dowhy4
Steven Siciliano3
Eric G. Lamb1- 1Department of Plant Sciences, College of Agriculture and Bioresources, University of Saskatchewan, Saskatoon, SK, Canada
- 2Department of Soil Science, College of Agriculture and Bioresources, University of Saskatchewan, Saskatoon, SK, Canada
- 3Saskatoon Research and Development Centre, Agriculture and Agri-Food Canada, Saskatoon, SK, Canada
- 4Department of Computer Science, College of Arts and Science, University of Saskatchewan, Saskatoon, SK, Canada
- 5Department of Animal and Poultry Science, College of Agriculture and Bioresources, University of Saskatchewan, Saskatoon, SK, Canada
Modifying the rhizosphere microbiome through targeted plant breeding is key to harnessing positive plant–microbial interrelationships in cropping agroecosystems. Here, we examine the composition of rhizosphere bacterial communities of diverse Brassica napus genotypes to identify: (1) taxa that preferentially associate with genotypes, (2) core bacterial microbiota associated with B. napus, (3) heritable alpha diversity measures at flowering and whole growing season, and (4) correlation between microbial and plant genetic distance among canola genotypes at different growth stages. Our aim is to identify and describe signature microbiota with potential positive benefits that could be integrated in B. napus breeding and management strategies. Rhizosphere soils of 16 diverse genotypes sampled weekly over a 10-week period at single location as well as at three time points at two additional locations were analyzed using 16S rRNA gene amplicon sequencing. The B. napus rhizosphere microbiome was characterized by diverse bacterial communities with 32 named bacterial phyla. The most abundant phyla were Proteobacteria, Actinobacteria, and Acidobacteria. Overall microbial and plant genetic distances were highly correlated (R = 0.65). Alpha diversity heritability estimates were between 0.16 and 0.41 when evaluated across growth stage and between 0.24 and 0.59 at flowering. Compared with a reference B. napus genotype, a total of 81 genera were significantly more abundant and 71 were significantly less abundant in at least one B. napus genotype out of the total 558 bacterial genera. Most differentially abundant genera were Proteobacteria and Actinobacteria followed by Bacteroidetes and Firmicutes. Here, we also show that B. napus genotypes select an overall core bacterial microbiome with growth-stage-related patterns as to how taxa joined the core membership. In addition, we report that sets of B. napus core taxa were consistent across our three sites and 2 years. Both differential abundance and core analysis implicate numerous bacteria that have been reported to have beneficial effects on plant growth including disease suppression, antifungal properties, and plant growth promotion. Using a multi-site year, temporally intensive field sampling approach, we showed that small plant genetic differences cause predictable changes in canola microbiome and are potential target for direct and indirect selection within breeding programs.
Introduction
The plant microbiome is a promising avenue of exploration to enhance crop productivity and management. Recent studies are revealing that plant breeding can shape the composition of root-associated bacterial communities including enhancing antagonistic potential toward pathogens (Peiffer and Ley, 2013; Bouffaud et al., 2014; Cardinale et al., 2015). Promising results that indicate both microbiome heritability and associations with yield have been reported. For example, a study on 27 maize inbred lines indicated the presence of a small but significant proportion of heritable variation in total bacterial diversity across field environments and substantially more heritable variation between replicates of lines within each field (Peiffer et al., 2013). A study of Brassica napus root-associated microbiomes in the Canadian Prairies identified bacterial taxa which were positively correlated with canola yield (Lay et al., 2018). Identifying genetically controlled positive plant–microbial interactions by comparing lines within breeding programs and across diversity panels is the first step in determining if plant breeders could develop varieties by selecting for genetic factors controlling beneficial plant–microbial interrelationships.
The rhizosphere microbiome has commonly been targeted to identify positive plant–microbial relationships. Past microbiome studies using culture dependent and independent approaches have shed light on the structure and composition of the plant microbiota in many crops (Garland, 1996; Germida et al., 1998; Bulgarelli et al., 2012; Lundberg et al., 2012; Peiffer et al., 2013; Schlaeppi et al., 2013; Edwards et al., 2015; Lay et al., 2018). Continuing advances in biotechnology and bioinformatics are enabling researchers to evaluate the microbiome to a greater depth, incorporating more replications and to account for variables such as genotype, time, and space.
Our goal was to characterize the core rhizosphere microbiome and identify contrasting components of the bacterial microbiota between plant genotypes of B. napus (canola) grown under field conditions. Our specific objectives were to: (1) identify bacterial taxa differentially abundant between multiple canola genotypes, (2) screen taxa that were differentially abundant between genotypes to identify potential beneficial plant–microbial interactions, (3) characterize the core canola rhizosphere bacterial microbiota across the full growing season at a single site, (4) examine the consistency of the core microbiota across years and location, (5) estimate the heritability of alpha diversity measures at flowering and for the whole growing season, and (6) determine correlation between the microbial and plant genetic distances among canola genotypes at different growth stages.
Materials and Methods
Site Description
Experiments were conducted in 2016 at a single location and at three locations in 2017. The 2016 experimental site was near Saskatoon, SK, Canada (latitude 52.181366, longitude -106.502941). The soil is a haplic Kastanozem (Iuss Working Group WRB, 2015), with a clay loam texture. The field was managed with an oilseed, wheat, barley, fallow rotation with fallow in the growing season prior to this experiment. Pre-trial field nutrient assessments guided pre-trial fertilizer applications of 78.5 kg ha–1 NH3 in the fall of 2015 with N21.5, P35 S25 being applied prior to planting in the spring of 2016. EdgeTM (Ethalfluralin – Group 3 – dinitroanaline) (Gowan Company, Yuma, AZ, United States) a pre-emergence herbicide, was also applied to the field at a rate of 19.1 kg ha–1 for weed control. The 2017 locations were at Saskatoon (latitude 52.183149, longitude -106.514904), Melfort (latitude 52.819333, longitude −104.596348), and Scott (latitude 52.365370, longitude −108.875710), SK.
Canola Genotypes
Sixteen diverse B. napus genotypes varying in seed color, glucosinolate, erucic acid, and fiber contents (Table 1) were used in this study to represent genetically diverse germplasm. Five of the B. napus lines were breeding lines from the Agriculture and Agri-Food Canada (AAFC) canola breeding program, whereas the other 11 lines were selected from a larger diversity collection comprised of accessions representing the variation found across spring B. napus housed within various germplasm collections. A genetic similarity matrix for the 16 B. napus genotypes was generated based on single nucleotide polymorphisms (SNPs) across the genome determined using the Brassica 60K Illumina Infinium SNP array (Clarke et al., 2016). Two of the genotypes (Table 1) used in this study have alternate names, NAM-48 or DH27298 and NAM-94 or YN04-C1213sp013.
Table 1. Sixteen diverse Brassica napus genotypes selected for the rhizosphere microbial analysis and their seed quality traits.
Experimental Design
The field design in 2016 was a randomized complete block design with three replicates. The three blocks each were comprised of two strips of eight 2 m by 6 m plots. The 16 B. napus genotypes were randomly assigned to plots and seeded on May 27, 2016. B. napus seeds were pre-treated with HELIX XTra® (active ingredients: Thiamethoxam, difenoconazole, metalaxyl-M, fludioxonil, and sedaxane) (Syngenta Canada Inc., Guelph, ON, United States) and were seeded at a rate of 100 seeds m–2 with HELIX-treated corn grits. Roots with attached rhizosphere and bulk soil were collected weekly for 10 weeks starting 18 days after seeding (DAS) and ending 81 DAS. The corresponding phenotypic stages for the 2016 sampling weeks were 2–3 leaf stage (Week 1, 18 DAS), 4–5 leaf stage (Week 2, 25 DAS), 6–9 leaf stage (Week 3, 32 DAS), flowering (Weeks 4–7, 39–70 DAS), maturity (Weeks 8 and 9, 67 and 74 DAS), and harvest (Week 10, 81 DAS). Plants were harvested at maturation to determine seed yield (g) for the entire plot. Each plot-level sample was comprised of three composited subsamples. A subsample was an individual canola plant taken at random from within a plot. For each subsample, bulk and rhizosphere soils were separated, with rhizosphere soil defined as soil that remained attached to roots after shaking to dislodge the loosely attached (bulk) soil (Supplementary Figure S1). Roots of each sub-sample with the adhered rhizosphere soil were handled aseptically, cut and placed together in an Erlenmeyer flask with 100 mL of buffer (0.05 M NaCl), and shaken (Innova 2100 Platform Shaker at 180 r/min) for 15 min. Following shaking, roots were removed from the flask, and the buffer and rhizosphere soil were transferred into two 50 mL centrifuge tubes. Tubes containing the sample mixtures were then centrifuged at 5000 r/min for 15 min at room temperature. The supernatant was decanted and pellets of the rhizosphere soil were subsampled into two 1.5 mL Eppendorf tubes and stored frozen at −80°C.
In 2017, each of the three sites were established following a similar experimental design and sampling protocol to 2016. Seeding of the 16 diverse B. napus genotypes was done on May 29 at Saskatoon and May 19 at Melfort. At Scott, delayed seeding date of June 20 was due to re-seeding after hail damage. Rhizosphere soil samples were collected three times in 2017 at the 6–9 leaf stage (Week 3), mid-flowering (Week 6), and maturity (Week 9). Plants were harvested at maturation to determine seed yield (g) for the entire plot.
DNA Extraction, Amplification, and Sequencing
Rhizosphere soil pellets were removed from the −80°C freezer and 0.1 g was transferred into 96-well DNA extraction plates. DNA was extracted using the DNeasy PowerSoil Kit (Qiagen, catalogue number 12955-4) following the recommended standard procedure. The extracted DNA quantity was determined using a standard Qubit protocol (Thermo Fisher Scientific, Waltham, MA, United States). 16S rRNA genes were amplified using the primer set 342F-806R with Illumina sequencing adapters (342F: 5′-ACA CTG ACG ACA TGG TTC TAC ACT ACG GGG GGC AGC AG-3′ and 806R: 5′-TAC GGT AGC AGA GAC TTG GTC TGG ACT ACC GGG GTA TCT-3′) (Mori et al., 2014). Amplified PCR products were then sent to the Innovation Centre at Genome Quebec for amplicon barcoding, normalization, library QC, and subsequently sequenced on an Illumina MiSeq. Sequencing of the 2017 samples was done at the University of Saskatchewan sequencing facility following the same procedure using an Illumina MiSeq platform. A total of 887 samples were analyzed. Replicates for DNA extraction were individually sequenced. To compare absolute abundance of bacterial 16S rRNA genes across samples, a known amount of a bacterial species not expected to be present in soil (0.3 ng μL–1 of Aliivibrio fisheri) was spiked into each sample as an internal standard (Smets et al., 2016).
Sequence Data Processing and Statistical Analysis
Adaptors and primers were trimmed from demultiplexed paired reads using Cutadapt version 2.1 (Martin, 2011). The Cutadapt processed reads were further processed using qiime2 version 2019.1 (Bolyen et al., 2018) DADA2 plugin (Callahan et al., 2016) to filter low-quality and chimera errors and generate a final amplicon sequence variants table and corresponding taxonomic table. DADA2 pipeline is designed to resolve exact biological sequences from Illumina sequence data and does not involve sequence clustering. Following Knight et al. (2018), exact sequence variant approach was used. Oligotyping (Eren et al., 2013) improves upon traditional operational taxonomic unit (OTU) picking, by including position-specific information from 16S rRNA sequencing to identify subtle nucleotide variation and by discriminating between closely related but distinct taxa. The sequences were mapped at 99% sequence identity to an optimized version of the Greengenes reference database (version 13.8) containing the target V3–V5 16S region to determine taxonomies at seven different levels.
The amplicon sequence variant table was further processed to: (1) remove 16S rRNA gene sequences identified as chloroplast, mitochondria, or archaea, (2) adjust sequence counts to account for differences in sequencing depth based on the internal standard (A. fisheri counts), and (3) remove A. fisheri sequences. We rounded A. fisheri adjusted values to the nearest integer prior to analysis with negative binomial generalized linear models (GLMs) which require count data. All other summary statistics and analyses were completed using the unrounded values unless stated. Statistical analyses were done in R (version 3.6.0, R Core Team, 2019) using Phyloseq (version 1.22.3) (McMurdie and Holmes, 2013), microbiome (version 1.5.28) (Leo and Sudarshan, 2017), and their associated dependencies.
The B. napus core microbiota, or the set of amplicon sequence variants detected in 50–100% of the samples with a relative abundance threshold value above 0.01%, was identified using the core function in microbiome R package version 1.5.28 (Leo and Sudarshan, 2017) using the whole dataset from Saskatoon in 2016. The dataset was then grouped into three broad phenological stages: vegetative (weeks 1–3), flowering (weeks 4–7), and maturity (weeks 8–10). Core bacterial taxa with at least 75% prevalence were determined within each phenological category and were compared with the core taxa identified for the whole dataset to determine if there were patterns in how the core taxa were recruited. Core bacterial taxa with at least 75% prevalence were determined for each of the three sampling sites in 2017, and then compared between sites and to the core bacterial taxa from the single site 2016 experiment.
Prior to α-diversity calculations and permutational multivariate analysis of variance (PERMANOVA) test (Anderson, 2001), the amplicon sequence variant table was normalized using the edgeR “edgernorm” method in R package microbiomeSeq (Ssekagiri et al., 2017). The α-diversity indices (Richness, Pielou’s evenness) calculations, pairwise ANOVA of the diversity measures between canola genotypes as well as between sampling weeks were done using microbiomeSeq and microbiome (Leo and Sudarshan, 2017) packages. Bray–Curtis distance for each pair of samples was calculated and PEMANOVA and homogeneity of dispersion test were run using the “adonis” and “betadisper” functions in R package Vegan (Oksanen et al., 2019). The 2016 10-week dataset and its subsets vegetative, flowering, and maturity stages were used for α-diversity and vegetative, flowering, vegetative and flowering combined, and flowering and maturity stages combined subsets for PERMANOVA analysis.
We calculated the mean Bray–Curtis distance between canola genotypes by averaging the paired sample values of each group. Person’s correlation coefficient between plant genetic distance and the mean microbial Bray–Curtis distance among the genotypes were then calculated. The datasets that showed difference in microbial community structure among canola genotypes (PERMANOVA test) were used for this analysis.
We considered the diversity indexes calculated as phenotypic records and hence estimated their broad-sense heritability (h2) separately for flowering and whole 2016 datasets. To estimate h2, first variance components of each variable (richness, evenness, and other calculated indices) were generated by fitting a linear mixed effects model using restricted maximum likelihood (REML) with the lmer function of the lme4 R package version 1.1.21 (Bates et al., 2019). The two models where: response (whole data) = canola genotype + Rep (Week) + canola X week + residual, and response (flowering) = canola + canola X week + residual. Due to model over parameterization and singular fit “Rep (week)” (replication nested in sampling week) from whole data model and “canola X week” (canola sampling week interaction) from the flowering stage model were dropped in the final model. The h2 decreased by a maximum of only 0.01 in the simpler models. The proportion of variance explained by canola genotype (h2) was calculated using the formula h2 = σ G2/(σG2 + (σGx W 2)/w + σe2/rw) for the whole 2016 data, or h2 = σG2/(σG2 + σe2/rw) for flowering stage data, where σG2 is genetic variance, σGx W2 is the variance of genotype by sampling week interactions, σe2 is residual variance, r is the number of replications in each sampling week, and w is the number of sampling weeks.
Bacterial genera that were differentially abundant between each of the 15 B. napus genotypes and a reference genotype (NAM-0) were identified using negative binomial models fitted in GLM framework using the edgeR package version 3.8 (Robinson et al., 2010). The 2016 Saskatoon dataset was first agglomerated to the lowest identified taxonomic rank (genus), resulting in a total of 558 genera. First, Upperquartile (Bullard et al., 2010) normalization was used to align the upper quartiles of the count per million within the libraries. The goal was to identify bacterial taxa that were differentially abundant between the reference genotype and each of the 15 remaining B. napus genotypes. Hence, the design matrix contained B. napus genotype with the intercept set as NAM-0, and a significant coefficient for a given genotype identified significant abundance differences from NAM-0. Dispersion estimates were subsequently calculated using estimateDisp for the given design matrix. A quasi-negative-binomial-GLM was fit with the glmQLFit function, and significance tested with the glmQLFtest function. A 1% false discovery rate (FDR; Benjamin and Hochberg adjusted p-value) was used. R scripts for the core bacteria, heritability and differential abundance analysis are provided (Supplementary Data Sheets S1–S3).
Results
Rhizosphere Taxonomic Characterization
A total of 1,571,433,759 reads (ranging from 5925 to 354,234,808 per sample) across 477 samples (Supplementary Table S1) were present in the dataset. We detected 32 named bacterial phyla (Supplementary Table S2 and Supplementary Figure S2), approximately a third of the currently named phyla (Hug et al., 2016). Proteobacteria, Acidobacteria, Actinobacteria, and Chloroflexi were the most abundant phyla, present in 69.4, 13.8, 12.2, and 4.6% of the samples, respectively (Supplementary Tables S3, S4). Four named classes were identified from phylum Proteobacteria which included the common Gammaproteobacteria, Alphaproteobacteria, Betaproteobacteria, and Deltaproteobacteria, 15 classes from Acidobacteria, and 10 from Actinobacteria.
Single Site 2-Year B. napus Rhizosphere Core Microbiome
Six bacterial genera were identified as core at a 75% prevalence threshold; at 50, 60, 65, and 70% prevalence thresholds the number of core genera were 32, 15, 11, and 8, respectively. The B. napus core microbiome primarily included Proteobacteria and Actinobacteria. Arthrobacter (Actinobacteria) was found in 96% of the samples, Bradyrhizobium (Proteobacteria) in 95%, Stenotrophomonas (Proteobacteria), Skermanella (Proteobacteria), and one unclassified Acidobacteria and one unclassified Actinobacteria were present in 75% of the samples (Table 2 and Supplementary Table S5). Three of the four named core taxa in 2016 were consistently observed in the 2017 experiment. A total of 11 core bacterial taxa were identified at 75% in the 2017 experiment (Table 3 and Supplementary Table S5).
Table 2. Overall and growth stage-related B. napus core bacterial taxa identified at Saskatoon 2016 experiment.
Table 3. Cross-year B. napus core bacterial taxa identified in at least 75% of the samples at Saskatoon site.
Core genera associated with each growth stage were generally part of the overall core taxa at the 75% prevalence threshold. Arthrobacter, Bradyrhizobium, and an unclassified Acidobacteria in the class Ellin6075 were present in all growth stages, while two variants of Arthrobacter, Stenotrophomonas, and Skermanella joined the core taxa at flowering. One unclassified Gaiellaceae (Actinobacteria) is observed only at the vegetative and flowering stages. Acinetobacter a Gammaproteobacteria in family Moraxellaceae was only present during flowering while Agromyces an Actinobacteria in family Microbacteriaceae was only present in the maturity stage (Table 2).
2017 Three-Site B. napus Rhizosphere Core Microbiome
The 2017 cross-site core microbiome closely resembled the 2016 core microbiome with 11, 21, and 7 core bacterial genera at 75% prevalence at Saskatoon, Melfort, and Scott sites, respectively. The cross-site B. napus core microbiome primarily included Proteobacteria and Actinobacteria. Arthrobacter, Erwinia spp. (two variants), Kaistobacter, Pedobacter, and Stenotrophomonas (Stenotrophomonas retroflexus) were cross-site conserved. Bradyrhizobium and Microlunatus were conserved at two sites (Saskatoon and Melfort) (Table 4 and Supplementary Table S5). Arthrobacter was the most prevalent core taxa at Saskatoon (both years) and Melfort and Stenotrophomonas was most prevalent at Scott.
Table 4. Brassica napus core bacterial taxa identified in at least 75% of the samples at Saskatoon, Melfort, and Scot sites.
Genotype Variations in Alpha Diversity Measures
Alpha diversity measures (richness and Pielou’s evenness) varied among the canola genotypes. More variability is observed in Pielou’s evenness. Flowering stage followed by whole data analysis resulted in the highest variability among genotypes (Figure 1). No significant difference among genotypes was observed at vegetative stage and only two pairs of genotypes (NAM-48 – NAM-76 and NAM-14 – NAM-43) showed significant difference in richness and Pielou’s evenness at maturity.
Figure 1. Variability in alpha diversity measures (richness and evenness) among canola genotypes. Bars connect significantly different canola genotype pairs and significance level is indicated with an asterisk (∗0.05, ∗∗0.01). Figures present alpha diversity comparisons based on (A) Whole 2016 dataset, (B) vegetative, (C) flowering, and (D) maturity stages.
Microbial Patterns Among B. napus Genotypes
The PERMANOVA and betadisper analyses were used to test whether the rhizosphere bacterial community structure differed among B. napus genotypes. The PERMANOVA partitions the variability between factors based on the Bray–Curtis dissimilarity matrix. To check the contribution of growth stage at which samples were collected, we used the complete 2016 dataset and its subsets. Statistically significant differences were observed in the rhizosphere microbiome among B. napus genotypes when using the complete (F = 1.23, R2 = 0.04, p < 0.001), flowering (F = 1.12, R2 = 0.09, p < 0.05), vegetative and flowering combined (F = 1.14, R2 = 0.05, p < 0.05), and flowering and maturity combined (F = 1.17, R2 = 0.05, p < 0.01) datasets. The test for homogeneity of multivariate dispersions was not significant indicating that canola genotypes have similar dispersion and we can trust the significant variation observed. We did not observe significant variation in the rhizosphere microbiome among B. napus genotypes at the vegetative and maturity stages.
Relationship Between B. napus Genotypes and Rhizosphere Microbial Composition
To check whether the composition of the rhizosphere microbiota could be correlated to the genetic relatedness among the B. napus genotypes, we computed mean microbial Bray–Curtis distances and correlated those with plant genetic distance among the genotypes. In addition, to infer if the correlation is affected by the growth stage of B. napus, we did correlations using datasets from flowering, vegetative and flowering combined, and flowering and maturity stages combined. The vegetative stage was omitted since we did not find significant variation during the PERMANOVA analysis. The significant (p < 0.001) and highest (R = 0.65) correlation was observed when using the dataset across the whole growing season followed by in the vegetative and flowering stage combined (p < 0.001, R = 0.55), and flowering and maturity stage combined datasets (p < 0.001, R = 0.44), respectively (Figure 2). At flowering the least but significant correlation was observed (p < 0.05, R = 0.21).
Figure 2. Correlation between mean microbial Bray–Curtis distance and plant genetic distance among canola genotypes. Change in correlation when considering the (A) whole 2016 dataset, (B) flowering stage, (C) vegetative and flowering stages combined, and (D) flowering and maturity stages combined.
Heritability of Microbial Attributes
The heritability estimates for the alpha diversity measures were higher at flowering than for the whole growing season (Tables 5, 6). At flowering, diversity inverse Simpson and evenness Simpson were with the highest heritability (h2 = 0.59) followed by diversity Shannon and evenness Pielou (h2 = 0.51). Diversity gini Simpson (h2 = 0.37) followed by diversity inverse Simpson and evenness Simpson (h2 = 0.25) were with the highest heritability estimates for the whole growing season.
Table 5. Broad-sense heritability of alpha diversity traits (diversity_inverse_simpson, diversity_gini_simpson, diversity_shannon, evenness_pielou, evenness_simpson) using the full 10 weeks dataset of 2016.
Table 6. Broad-sense heritability of alpha diversity traits (diversity_inverse_simpson, diversity_gini_simpson, diversity_shannon, evenness_pielou, evenness_simpson) using the flowering stage dataset of 2016 experiment.
Differentially Abundant Bacterial Genera
Of the total 558 genera, 81 bacterial genera were significantly (p < 0.01) more and 71 genera were similarly less abundant in at least one B. napus genotype relative to the reference genotype (NAM-0). With the exception of genotypes NAM-14, NAM-30, NAM-32, NAM-43, and NAM-46, there was an observed general trend toward differentially less abundant taxa per genotype (Table 7 and Figure 3). The most genetically distinct genotype (NAM-46) had the greatest number of differentially highly abundant taxa; however, there was not a clear pattern linking plant genetic distance across the diversity panel with number of differentially abundant bacterial taxa (Table 7 and Supplementary Figure S3).
Table 7. Number of differentially less and more abundant bacterial genera at 1% false discovery rate in 15 B. napus lines compared with the reference NAM-0.
Figure 3. Differentially abundant bacterial genera compared with the reference genotype (NAM-0). Change in absolute abundance in B. napus genotypes (log fold change) is shown against average abundance in count per million (CPM). Red dot indicates significantly differentially more and blue less abundant taxa (FDR < 0.01). The non-significant genera are indicated in gray. Names of B. napus genotypes are indicated at the top left corner of individual plots. Plots are arranged left to right from genotype with the highest to the lowest number of differentially abundant genera compared with NAM-0.
A total of 29 bacterial phyla (17–23 phyla per genotype) were represented among the differentially abundant bacterial genera. The majority were Proteobacteria and Actinobacteria, together accounting for 52–61% of differentially abundant genera in each B. napus genotype. Phylum FBP, Tenericutes, and BRC1were unique phyla observed only in NAM-13, NAM-14, and NAM-48, respectively.
Across the 15 genotypes, there were 67 genotype-specific (Table 8) differentially abundant bacterial genera of which 56 (84%) were differentially more abundant and 11 (16%) less abundant compared with the reference NAM-0. Per B. napus genotype, the number of unique genera ranged from one [NAM-23, ∼2% (per total number of genotype-specific genera)] to 11 (NAM-46, 16%). In NAM-76, all unique genera were differentially less abundant whereas in NAM-5, NAM-13, NAM-14, NAM-23, NAM-32, NAM-46, and NAM-94, the unique genera were differentially more abundant.
Table 8. Differentially abundant named bacterial genera unique to each B. napus genotypes in 2016 and their traits or plant beneficial roles.
We ranked differentially abundant bacteria based on abundance changes relative to the reference line [log fold change (logFC)] to identify bacteria with the strongest population responses to each B. napus genotype. Responses in both positive (more abundant) and negative (less abundant) directions were observed. For example, in the genetically distinct genotype NAM-46 there were more genera with strong positive responses, this also seemed the case for most lines (Supplementary Figure S4). The full list of the significantly differentially abundant bacterial genera within each genotype together with their corresponding taxonomic classification and values for logFC, log count per million (logCPM), p-value, and FDR are provided in Supplementary Table S6.
Discussion
B. napus Genotype Regulated Changes in Rhizosphere Bacteria
Here we show fine-scale regulation of the rhizosphere bacterial community by B. napus host genotypes as a total of 67 differentially abundant genera were genotype-specific. This suggests an extensive and selective control by B. napus genotype on associated rhizosphere bacterial genera. Given these controls are genetically based, they may represent potential breeding targets if the associated bacterial shown to be positively associated with yield or positive traits in subsequent work. The most genetically distinct B. napus genotype (NAM-46) contributed a greater proportion (∼16%) of the genotype specific differentially abundant taxa, and almost all of its top 20 differentially abundant genera based on logFC were differentially more abundant (Supplementary Figure S4). This suggests that genetic variability among genotypes might be directed toward investing in recruiting taxa rather than efforts toward exclusion (reduction in abundance). The phyla that were mainly responsive to B. napus line were Proteobacteria and Actinobacteria; similar responsiveness of these phyla to dynamic changes in bacterial community has been reported (Yang et al., 2017).
Several potentially important beneficial B. napus genotype–bacterial interrelationships were evident among the differentially abundant taxa. Differentially abundant bacteria with potential benefits include Cellulosimicrobium, Enterobacter, Klebsiella, Modestobacter, Neisseria, Peredibacter, Pseudomonas, Serratia, and Streptosporangium. These bacteria have been reported to have roles in the biological control of diseases, including production of antibiotics, association with disease suppressive soils, anti-fungal properties, and biocontrol mechanisms. Cellulosimicrobium spp., for example, inhibit growth of the pathogenic fungi Botryis cinerea, Fusarium oxysporum, and Verticillium dahlia in barley (Nabti et al., 2014). Enterobacter spp. produce complex compounds that are antagonistic toward many fungal pathogens and also have antibiotic properties (Chernin et al., 1995, 1996). Neisseria is associated with disease suppressive soils (Almario et al., 2013; Adam et al., 2014), Peredibacter suppresses soil-borne disease (Jaiswal et al., 2017), and Pseudomonas induces systemic resistance against pathogens such as the soil-borne Rhizoctonia solani (Nandakumar et al., 2001). Streptosporangium is an Actinomycete known to produce antimicrobial compounds (Hardoim et al., 2015). Finally, Serratia promotes resistance to fungal pathogens in rapeseed (Alström, 2001; Neupane et al., 2013). It can also parasitize pathogenic Fusarium and reduces the production of areal hyphae and microconidia (Minerdi et al., 2008).
In addition to disease protection, some of the differentially abundant bacteria have beneficial effects including growth promotion. Cellulosimicrobium spp. can stimulate growth of barley seedlings (Nabti et al., 2014), and Enterobacter spp. are important plant growth-promoting bacteria (Chernin et al., 1996; Nie et al., 2002). Strains of Klebsiella isolated from wheat rhizosphere produce indole acetic acid (IAA) with demonstrated promotion of root growth (Sachdev et al., 2009). Serratia is a plant growth-promoting bacteria and is positively correlated with canola yield (Lay et al., 2018) and promotes growth in rapeseed (Alström, 2001; Neupane et al., 2013). Finally, Pseudomonas, a well-studied plant growth-promoting bacterium enhances root elongation of canola and is characterized by its ability to produce phytohormones, IAA, and cytokinin (Pallai et al., 2012).
The potential roles of the differentially abundant bacterial taxa discussed above are limited to genera with functions previously reported in the literature; examination of the uncharacterized and unclassified bacteria associated with high-performing B. napus genotypes will likely reveal additional positive plant–microbial interrelationships. Evaluating genotypes with contrasting differential abundance of bacteria that influence traits of interest such as disease resistance and yield would be the next step toward practical applications of B. napus microbiome manipulation. The further dissection of the highlighted B. napus–bacterial interrelationships is crucial in conceptualizing their use as potential targets for direct or indirect selection within breeding programs.
The observed high correlation between the mean microbial Bray–Curtis distance and plant genetic distance (R = 0.65) indicates that overall microbial variation is a good predictor of plant genetic distance. Such observations under field conditions encourage further dissection of plant genetic control on the Canola microbiome. The high variability in species evenness compared with species richness also indicates as there could be direct processes including changes in rhizosphere soil chemistry or other indirect effects by host genotype that alters proportional diversity through changes in evenness without changes in species composition (Wilsey and Potvin, 2000). These potential host-driven changes showed maximum effect during flowering stage with close to no effects during vegetative and maturity stages to grant significant variation in evenness among the genotypes. This suggests that plant genetics-driven changes start to be apparent around flowering stage.
This study is also among the first to describe alpha diversity measures as phenotypic traits in canola and to estimate heritability both considering the whole growing season and at flowering stage. The heritability range for different traits (0.37–0.59) at flowering stages indicates that 37–59% of the variation in alpha diversity measures in our canola genotypes at flowering was due to genetics. Across growth stages 16–37% of the variation in our canola genotypes was either directly or indirectly due to genetics. Diversity is generally associated with resilience to perturbations and disease suppression in rhizosphere, thus diversity could be a valuable trait to consider for sustainable and improved crop production. Maximum heritability values were observed at flowering, again suggesting maximum genetic effects at this growth stage. Therefore, if the alpha diversity measures are to be used as indicator traits, we suggest selections to be made at the flowering stage.
B. napus Rhizosphere Core Microbiome
Microbes that are consistently present across all cultivars and sites are likely to provide critical ecological functions (Shade and Handelsman, 2012). Our cross-site core taxa analysis revealed taxa that were conserved across sites (Table 4). This suggests that B. napus genotypes could select relatively conserved core bacterial taxa across a variety of soil types and environments. In contrast to most core microbiome studies, which use samples collected at a single time point, our sampling approach (10 weeks, at least 28 samples per canola genotype) provides insights into how the core members associate with the plant through growth and development stages. Separating the core taxa for each of the three growth stages yielded three key insights into the ecology of the core microbiome. First, there were core members that remained present throughout the development of the crop; second, there were taxa that joined the core membership starting at flowering and continued through to maturity and third, there were taxa that joined the core membership only at a particular growth stage but were not retained in other stages. These findings suggest that efforts to determine the core microbiota of crop species that are limited in the volume and/or timing of sampling should focus on the periods of flowering and/or maturity to get the best overall representation. Additional similar study in a different crop is needed to see if the same phenomenon holds. Cross-site and cross-year stability of core members such as Arthrobacter, Bradyrhizobium, and S. retroflexus indicate their close association with canola.
Core microbiota members with potential beneficial effects include Arthrobacter, Stenotrophomonas, and Bradyrhizobium. Arthrobacter can increase canola yield when applied as a bacterial suspension to seeds (Kloepper et al., 1988). Stenotrophomonas spp. including Stenotrophomonas rhizophila have antagonistic activity against fungal pathogens including V. dahlia in oil seed rape (Alström, 2001; Wolf et al., 2002). Arthrobacter and Stenotrophomonas are likely broadly present in the core microbiome of canola, having been identified in other surveys of the core rhizosphere microbiome of field grown canola (Lay et al., 2018) and in related Brassica rapa (Germida et al., 1998; Siciliano and Germida, 1999; Larcher et al., 2008; Croes et al., 2013; Tkacz et al., 2015). Positive correlations between Arthrobacter and Stenotrophomonas with canola yield (Lay et al., 2018) suggest an important beneficial influence of these bacteria on canola. We observed two variants of Arthrobacter joining the core microbiome during flowering, suggesting an important role in the reproductive phase of canola development. Finally, Bradyrhizobium shares characteristics with plant growth-promoting rhizobacteria, producing compounds including phytohormones, siderophores, and hydrogen cyanide, and exhibiting antagonistic effects toward many plant pathogenic fungi (Antoun et al., 1998). Core bacteria with such antagonistic effects toward pathogens can be targets to enhance suppressive conditions in the rhizosphere. For instance, a suppressive environment in the rhizosphere can contribute to the failure of invasion and disease development by Plasmodiophora brassicae, the causative agent of clubroot disease (Mattey and Dixon, 2015; Dixon, 2017). Focusing on potentially co-evolved core bacterial taxa with both antagonistic potential and a competitive advantage in the B. napus rhizosphere could be a strategy to boost the suppressive potential of the rhizosphere to provide a control approach for clubroot or other diseases.
Harnessing Plant–Microbial Interactions for Canola Breeding
Harnessing the potentially beneficial plant–microbe interactions identified here will require research in a number of directions. Our approach using core-microbiota and differential abundance analysis, correlations between microbial Bray–Curtis distance and plant genetic distance and heritability of alpha diversity measures shows that genotype of B. napus may drive very specific associations. Although potentially positive associations are cited throughout literature, testing canola genotype-specific relationships is key to establishing if selecting for microbiome as a trait is feasible. Once specific positive interactions are confirmed in canola, there are several approaches that could be taken to enhance the benefits. For example, one could manipulate the microbiome by breeding plants that have a competitive advantage for particular microbes (Bakker et al., 2012). Our differential abundance analysis identified bacterial taxa that appear selectively favored in each B. napus genotypes; closely examining the genetic differences between those genotypes could shed light on the genetic basis for specific plant–microbe relationships. Once a genetic basis for a plant–microbe relationship is identified, breeders can work to move that trait into new lines. Developing microbial consortia that are evolutionarily adapted to the host plant of interest is another approach to harness beneficial microbial functions in agricultural systems. This approach is a direct microbiome manipulation where inoculated bacterial consortia may serve to reduce the time required for the rhizosphere microbiome to achieve niche saturation and competitive exclusion of pathogens (Bakker et al., 2012). Further developing insight on our growth stage-related core bacteria and their potential benefits, microbial ecologists could design rationally optimized microbial consortia potentially with a high degree of effectiveness and persistence in the environment. Even if such inoculant design is based on the assessment of the close association with the host plant potentially in its target production area or region, unforeseen risks associated with bringing in inoculants on resident biodiversity and ecosystem functions (Hart et al., 2018) should be considered. While it is promising to see relatively high heritability estimate of broad microbial traits, the differential abundance and core bacterial analysis should expand toward identifying co-occurrence patterns, microbial hubs, evaluating the heritability of these taxa and linking them with plant genotype through genome-wide association studies. To conclude, this study has shown that in a realistic field setting plant genetics influences the canola microbiome, opening the promise of revolutionary ways to enhance canola productivity and sustainability.
Data Availability Statement
The datasets generated for this study can be accessed from the NCBI Sequence Read Archive (https://www.ncbi.nlm.nih.gov/Traces/study/?acc=PRJNA575004) under accession number PRJNA575004.
Author Contributions
SS and BH conceived and designed the experiments. SS, BH, ZT, EL, CN, MA, SV, IP, and JB performed the experiments. JB contributed in sequence library preparation and approved the manuscript. ZT, EL, TD, SM, SS, and ML designed bioinformatics pipelines for sequence data processing. ZT completed the final bioinformatics analysis. IP carried out the SNP analyses. SR created the genetic similarity matrix. ZT and EL set the objectives, analyzed the data, and wrote the manuscript. SS, BH, SV, IP, MA, CN, SM, ML, SR, and TD reviewed, commented, and approved the manuscript.
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.
Funding
This study was funded by the Canada First Research Excellence Fund under the Plant Phenotyping and Imaging Research Centre (P2IRC) project 1.3, Plant Pedological Phenotype: Phenotyping the Plant Microbiome at the University of Saskatchewan, and a Saskatchewan Wheat Development Commission Scholarship to ZT.
Acknowledgments
We thank research technician Alix Schebel, graduate students Shanay Williams-Johnson and Zayda Morales, and summer students Cordell Vandenderen, Yolanda Iannucci, Sarah Valentine, Triann Littlepine, Jennifer Dechka, and Kira Blomquist who were part of the lab and field crew of project 1.3. We thank Agriculture and Agri-Food Canada (AAFC) for assistance with DNA sequencing. We also thank Steve Shirtliffe and project 1.2 field crews.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.03007/full#supplementary-material
References
Abiala, M. A., Odebode, A. C., Hsu, S. F., and Blackwood, C. B. (2015). Phytobeneficial properties of bacteria isolated from the rhizosphere of maize in southwestern Nigerian soils. Appl. Environ. Microbiol. 81, 4736–4743. doi: 10.1128/AEM.00570-15
Adam, M., Westphal, A., Hallmann, J., and Heuer, H. (2014). Specific microbial attachment to root knot nematodes in suppressive soil. Appl. Environ. Microbiol. 80, 2679–2686. doi: 10.1128/AEM.03905-13
Almario, J., Kyselková, M., Kopecký, J., Ságová-Marečková, M., Muller, D., Grundmann, G. L., et al. (2013). Assessment of the relationship between geologic origin of soil, rhizobacterial community composition and soil receptivity to tobacco black root rot in Savoie region (France). Plant Soil 371, 397–408. doi: 10.1007/s11104-013-1677-71
Alström, S. (2001). Characteristics of bacteria from oilseed rape in relation to their biocontrol activity against Verticillium dahliae. J. Phytopathol. 149, 57–64. doi: 10.1046/j.1439-0434.2001.00585.x
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
Antoun, H., Beauchamp, C. J., Goussard, N., Chabot, R., and Lalande, R. (1998). Potential of Rhizobium and Bradyrhizobium species as plant growth promoting rhizobacteria on non-legumes: effect on radishes (Raphanus sativus L.). Plant Soil 204, 57–67. doi: 10.1023/A:1004326910584
Bakker, M. G., Manter, D. K., Sheflin, A. M., Weir, T. L., and Vivanco, J. M. (2012). Harnessing the rhizosphere microbiome through plant breeding and agricultural management. Plant Soil 360, 1–13. doi: 10.1007/s11104-012-1361-x
Bates, D., Maechler, M., Bloker, B., Walker, S., Christensen, R., Dai, B., et al. (2019). Package lme4 : Linear Mixed-Effects Models using “Eigen” and S4 Version 1.1-21. Available at: https://github.com/lme4/lme4/ (accessed August 29, 2019).
Bertrand, H., Nalin, R., Bally, R., and Cleyet-Marel, J. C. (2001). Isolation and identification of the most efficient plant growth-promoting bacteria associated with canola (Brassica napus). Biol. Fertil. Soils 33, 152–156. doi: 10.1007/s003740000305
Bolyen, E., Rideout, J. R., Dillon, M. R., Bokulich, N. A., Abnet, C., Al-Ghalith, G. A., et al. (2018). QIIME 2: reproducible, interactive, scalable, and extensible microbiome data science. PeerJ 6:e27295v2. doi: 10.7287/peerj.preprints.27295v2
Bouffaud, M. L., Poirier, M. A., Muller, D., and Moënne-Loccoz, Y. (2014). Root microbiome relates to plant host evolution in maize and other Poaceae. Environ. Microbiol. 16, 2804–2814. doi: 10.1111/1462-2920.12442
Boukaya, N., Goudjal, Y., Zamoum, M., Chaabane Chaouch, F., Sabaou, N., Mathieu, F., et al. (2018). Biocontrol and plant-growth-promoting capacities of actinobacterial strains from the Algerian Sahara and characterisation of Streptosporangium becharense SG1 as a promising biocontrol agent. Biocontrol Sci. Technol. 28, 858–873. doi: 10.1080/09583157.2018.1501466
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
Bullard, J. H., Purdom, E., Hansen, K. D., and Dudoit, S. (2010). Evaluation of statistical methods for normalization and differential expression in mRNA-Seq experiments. BMC Bioinformatics 11:94. doi: 10.1186/1471-2105-11-94
Callahan, B. J., Sankaran, K., Fukuyama, J. A., McMurdie, P. J., and Holmes, S. P. (2016). Bioconductor workflow for microbiome data analysis: from raw reads to community analyses. F1000Research 5:1492. doi: 10.12688/f1000research.8986.2
Cardinale, M., Grube, M., Erlacher, A., Quehenberger, J., and Berg, G. (2015). Bacterial networks and co-occurrence relationships in the lettuce root microbiota. Environ. Microbiol. 17, 239–252. doi: 10.1111/1462-2920.12686
Chen, Y. P., Rekha, P. D., Arun, A. B., Shen, F. T., Lai, W. A., and Young, C. C. (2006). Phosphate solubilizing bacteria from subtropical soil and their tricalcium phosphate solubilizing abilities. Appl. Soil Ecol. 34, 33–41. doi: 10.1016/j.apsoil.2005.12.002
Chernin, L., Brandis, A., Ismailov, Z., and Chet, I. (1996). Pyrrolnitrin production by an Enterobacter agglomerans strain with a broad spectrum of antagonistic activity towards fungal and bacterial phytopathogens. Curr. Microbiol. 32, 208–212. doi: 10.1007/s002849900037
Chernin, L., Ismailov, Z., Haran, S., and Chet, I. (1995). Chitinolytic Enterobacter agglomerans antagonistic to fungal plant pathogens. Appl. Environ. Microbiol. 61, 1720–1726.
Ciancio, A., Pieterse, C. M. J., and Mercado-Blanco, J. (2016). Editorial: harnessing useful rhizosphere microorganisms for pathogen and pest biocontrol. Front. Microbiol. 7:1620. doi: 10.3389/fmicb.2016.01620
Clarke, W. E., Higgins, E. E., Plieske, J., Wieseke, R., Sidebottom, C., Khedikar, Y., et al. (2016). A high-density SNP genotyping array for Brassica napus and its ancestral diploid species based on optimised selection of single-locus markers in the allotetraploid genome. Theor. Appl. Genet. 129, 1887–1899. doi: 10.1007/s00122-016-2746-2747
Compant, S., and Mathieu, F. (2013). “Use of DOPE-FISH tool to better visualize colonization of plants by beneficial bacteria? An example with Saccharothrix algeriensis NRRL b-24137 colonizing grapevine plants,” in Molecular Microbial Ecology of the Rhizosphere, ed. F. J. de Bruijn, (Hoboken, NJ: John Wiley & Sons, Inc), 929–931. doi: 10.1002/9781118297674.ch87
Costerousse, B., Schönholzer-Mauclaire, L., Frossard, E., and Thonar, C. (2018). Identification of heterotrophic zinc mobilization processes among bacterial strains isolated from wheat rhizosphere (Triticum aestivum L.). Appl. Environ. Microbiol. 84, e1715–e1717. doi: 10.1128/AEM.017150-1717
Croes, S., Weyens, N., Janssen, J., Vercampt, H., Colpaert, J. V., Carleer, R., et al. (2013). Bacterial communities associated with Brassica napus L. grown on trace element-contaminated and non-contaminated fields: a genotypic and phenotypic comparison. Microb. Biotechnol. 6, 371–384. doi: 10.1111/1751-7915.12057
Dai, Y., Yan, Z., Jia, L., Zhang, S., Gao, L., Wei, X., et al. (2016). The composition, localization and function of low-temperature-adapted microbial communities involved in methanogenic degradations of cellulose and chitin from Qinghai–Tibetan Plateau wetland soils. J. Appl. Microbiol. 121, 163–176. doi: 10.1111/jam.13164
Dixon, G. R. (2017). Managing clubroot disease (caused by Plasmodiophora brassicae Wor.) by exploiting the interactions between calcium cyanamide fertilizer and soil microorganisms. J. Agric. Sci. 155, 527–543. doi: 10.1017/S0021859616000800
Edwards, J., Johnson, C., Santos-Medellín, C., Lurie, E., Podishetty, N. K., Bhatnagar, S., et al. (2015). Structure, variation, and assembly of the root-associated microbiomes of rice. Proc. Natl. Acad. Sci. U.S.A. 112, E911–E920. doi: 10.1073/pnas.1414592112
Egamberdiyeva, D., and Höflich, G. (2003). Influence of growth-promoting bacteria on the growth of wheat in different soils and temperatures. Soil Biol. Biochem. 35, 973–978. doi: 10.1016/S0038-0717(03)00158-155
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
Eren, A. M., Maignien, L., Sul, W. J., Murphy, L. G., Grim, S. L., Morrison, H. G., et al. (2013). Oligotyping: differentiating between closely related microbial taxa using 16S rRNA gene data. Methods Ecol. Evol. 4, 1111–1119. doi: 10.1111/2041-210X.12114
Fan, Z. Y., Miao, C. P., Qiao, X. G., Zheng, Y. K., Chen, H. H., Chen, Y. W., et al. (2016). Diversity, distribution, and antagonistic activities of rhizobacteria of Panax notoginseng. J. Ginseng Res. 40, 97–104. doi: 10.1016/j.jgr.2015.05.003
Fhoula, I., Najjari, A., Turki, Y., Jaballah, S., Boudabous, A., and Ouzari, H. (2013). Diversity and antimicrobial properties of lactic acid bacteria isolated from rhizosphere of olive trees and desert truffles of tunisia. Biomed Res. Int. 2013:405708. doi: 10.1155/2013/405708
Garland, J. L. (1996). Patterns of potential C source utilization by rhizosphere communities. Soil Biol. Biochem. 28, 223–230. doi: 10.1016/0038-0717(95)001130-111
Germida, J. J., Siciliano, S. D., De Freitas, J. R., and Seib, A. M. (1998). Diversity of root-associated bacteria associated with field-grown canola (Brassica napus L.) and wheat (Triticum aestivum L.). FEMS Microbiol. Ecol. 26, 43–50. doi: 10.1016/S0168-6496(98)00020-28
Grady, E. N., MacDonald, J., Liu, L., Richman, A., and Yuan, Z. C. (2016). Current knowledge and perspectives of Paenibacillus: a review. Microb. Cell Fact. 15:203. doi: 10.1186/s12934-016-0603-607
Hacène, H., Daoudi-Hamdad, F., Bhatnagar, T., Baratti, J. C., and Lefebvre, G. (2000). H107, a new aminoglycoside anti-Pseudomonas antibiotic produced by a new strain of Spirillospora. Microbios 102, 69–77.
Hardoim, P. R., van Overbeek, L. S., Berg, G., Pirttilä, A. M., Compant, S., Campisano, A., et al. (2015). The hidden world within plants: ecological and evolutionary considerations for defining functioning of microbial endophytes. Microbiol. Mol. Biol. Rev. 79, 293–320. doi: 10.1128/mmbr.00050-14
Hart, M. M., Antunes, P. M., Chaudhary, V. B., and Abbott, L. K. (2018). Fungal inoculants in the field: is the reward greater than the risk? Funct. Ecol. 32, 126–135. doi: 10.1111/1365-2435.12976
Hou, J., Liu, W., Wang, B., Wang, Q., Luo, Y., and Franks, A. E. (2015). PGPR enhanced phytoremediation of petroleum contaminated soil and rhizosphere microbial community response. Chemosphere 138, 592–598. doi: 10.1016/j.chemosphere.2015.07.025
Hsieh, T. F., Huang, H. C., Mündel, H. H., Conner, R. L., Erickson, R. S., and Balasubramanian, P. M. (2005). Resistance of common bean (Phaseolus vulgaris) to bacterial wilt caused by Curtobacterium flaccumfaciens pv. flaccumfaciens. J. Phytopathol. 153, 245–249. doi: 10.1111/j.1439-0434.2005.00963.x
Huang, A., Teplitski, M., Rathinasabapathi, B., and Ma, L. (2010). Characterization of arsenic-resistant bacteria from the rhizosphere of arsenic hyperaccumulator Pteris vittata. Can. J. Microbiol. 56, 236–246. doi: 10.1139/W10-005
Huber, B., Herzog, B., Drewes, J. E., Koch, K., and Müller, E. (2016). Characterization of sulfur oxidizing bacteria related to biogenic sulfuric acid corrosion in sludge digesters. BMC Microbiol. 16:153. doi: 10.1186/s12866-016-0767-767
Hug, L. A., Baker, B. J., Anantharaman, K., Brown, C. T., Probst, A. J., Castelle, C. J., et al. (2016). A new view of the tree of life. Nat. Microbiol. 1:16048. doi: 10.1038/nmicrobiol.2016.48
IUSS Working Group WRB (2015). World Reference Base for Soil Resources 2014, Update 2015 International Soil Classification System for Naming Soils and Creating Legends for Soil Maps. World Soil Resources Reports No. 106. Rome: FAO.
Jaiswal, A. K., Elad, Y., Paudel, I., Graber, E. R., Cytryn, E., and Frenkel, O. (2017). Linking the belowground microbial composition, diversity and activity to soilborne disease suppression and growth promotion of tomato amended with biochar. Sci. Rep. 7:44382. doi: 10.1038/srep44382
Janarthine, S. R. S., and Eganathan, P. (2012). Plant growth promoting of endophytic Sporosarcina aquimarina SjAM16103 isolated from the pneumatophores of Avicennia marina L. Int. J. Microbiol. 2012, 1–10. doi: 10.1155/2012/532060
Janssen, J., Weyens, N., Croes, S., Beckers, B., Meiresonne, L., Van Peteghem, P., et al. (2015). Phytoremediation of metal contaminated soil using willow: exploiting plant-associated bacteria to improve biomass production and metal uptake. Int. J. Phytoremediation 17, 1123–1136. doi: 10.1080/15226514.2015.1045129
Jha, B., Gontia, I., and Hartmann, A. (2012). The roots of the halophyte Salicornia brachiata are a source of new halotolerant diazotrophic bacteria with plant growth-promoting potential. Plant Soil 356, 265–277. doi: 10.1007/s11104-011-0877-879
Kloepper, J. W., Hume, D. J., Scher, F. M., Singleton, C., Tipping, B., Leliberté, M., et al. (1988). Plant growth-promoting rhizobacteria on canola (Rapeseed). Plant Dis. 72, 42. doi: 10.1094/PD-72-42
Knight, R., Vrbanac, A., Taylor, B. C., Aksenov, A., Callewaert, C., Debelius, J., et al. (2018). Best practices for analysing microbiomes. Nat. Rev. Microbiol. 16, 410–422. doi: 10.1038/s41579-018-0029-29
Larcher, M., Rapior, S., and Cleyet-Marel, J. C. (2008). Bacteria from the rhizosphere and roots of Brassica napus influence its root growth promotion by Phyllobacterium brassicacearum. Acta Bot. Gall. 155, 355–366. doi: 10.1080/12538078.2008.10516116
Lay, C., Bell, T. H., Hamel, C., Harker, K. N., Mohr, R., Greer, C. W., et al. (2018). Canola root – associated microbiomes in the Canadian Prairies. Front. Microbiol. 9:1188. doi: 10.3389/fmicb.2018.01188
Lee, M., Woo, S.-G., and Ten, L. N. (2011). Shinella daejeonensis sp. nov., a nitrate-reducing bacterium isolated from sludge of a leachate treatment plant. Int. J. Syst. Evol. Microbiol. 61, 2123–2128. doi: 10.1099/ijs.0.026435-0
Leo, L., and Sudarshan, S. (2017). Tools for Microbiome Analysis in R, Version 1.5.28. Available at: http://microbiome.github.com/microbiome
Lu, T., Ke, M., Lavoie, M., Jin, Y., Fan, X., Zhang, Z., et al. (2018). Rhizosphere microorganisms can influence the timing of plant flowering. Microbiome 6:231. doi: 10.1186/s40168-018-0615-10
Lundberg, D. S., Lebeis, S. L., Paredes, S. H., Yourstone, S., Gehring, J., Malfatti, S., et al. (2012). Defining the core Arabidopsis thaliana root microbiome. Nature 488, 86–90. doi: 10.1038/nature11237
Martin, M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 17:10. doi: 10.14806/ej.17.1.200
Mattey, M., and Dixon, G. R. (2015). Premature germination of resting spores as a means of protecting brassica crops from Plasmodiphora brassicae Wor., (clubroot). Crop Prot. 77, 27–30. doi: 10.1016/j.cropro.2015.07.013
McMurdie, P. J., and Holmes, S. (2013). Phyloseq: an R Package for reproducible interactive analysis and graphics of microbiome census data. PLoS One 8:e61217. doi: 10.1371/journal.pone.0061217
Minerdi, D., Moretti, M., Gilardi, G., Barberio, C., Gullino, M. L., and Garibaldi, A. (2008). Bacterial ectosymbionts and virulence silencing in a Fusarium oxysporum strain. Environ. Microbiol. 10, 1725–1741. doi: 10.1111/j.1462-2920.2008.01594.x
Mobasseri, M., Hutchinson, M. C., Afshar, F. J., and Pedram, M. (2019). New evidence of nematode-endosymbiont bacteria coevolution based on one new and one known dagger nematode species of Xiphinema americanum-group (Nematoda, Longidoridae). PLoS One 14:e0217506. doi: 10.1371/journal.pone.0217506
Mori, H., Maruyama, F., Kato, H., Toyoda, A., Dozono, A., Ohtsubo, Y., et al. (2014). Design and experimental application of a novel non-degenerate universal primer set that amplifies prokaryotic 16S rRNA genes with a low possibility to amplify eukaryotic rRNA genes. DNA Res. 21, 217–227. doi: 10.1093/dnares/dst052
Muthukumarasamy, R., Gopalakrishnan, R., Vadivelu, M., and Arun, K. (2017). Isolation of bacterial strains possessing nitrogen-fixation, phosphate and potassium-solubilization and their inoculation effects on sugarcane. Indian J. Exp. Biol. 55, 161–170.
Nabti, E., Bensidhoum, L., Tabli, N., Dahel, D., Weiss, A., Rothballer, M., et al. (2014). Growth stimulation of barley and biocontrol effect on plant pathogenic fungi by a Cellulosimicrobium sp. strain isolated from salt-affected rhizosphere soil in northwestern Algeria. Eur. J. Soil Biol. 61, 20–26. doi: 10.1016/j.ejsobi.2013.12.008
Nandakumar, R., Babu, S., Viswanathan, R., Raguchander, T., and Samiyappan, R. (2001). Induction of systemic resistance in rice against sheath blight disease by Pseudomonas fluorescens. Soil Biol. Biochem. 33, 603–612. doi: 10.1016/j.micres.2016.08.007
Neupane, S., Goodwin, L. A., Högberg, N., Kyrpides, N. C., Alström, S., Bruce, D., et al. (2013). Non-contiguous finished genome sequence of plant-growth promoting Serratia proteamaculan S4. Stand. Genomic Sci. 8, 441–449. doi: 10.4056/sigs.4027757
Nie, L., Shah, S., Rashid, A., Burd, G. I., Dixon, D. G., and Glick, B. R. (2002). Phytoremediation of arsenate contaminated soil by transgenic canola and the plant growth-promoting bacterium Enterobacter cloacae CAL2. Plant Physiol. Biochem. 40, 355–361. doi: 10.1016/S0981-9428(02)01375-X
Oksanen, J., Blanchet, F. G., Friendly, M., Kindt, R., Legendre, P., Mcglinn, D., et al. (2019). Package “vegan” Title Community Ecology Package. Available at: https://cran.r-project.org/web/packages/vegan/ (accessed August 29, 2019).
Pallai, R., Hynes, R. K., Verma, B., and Nelson, L. M. (2012). Phytohormone production and colonization of canola (Brassica napus L.) roots by Pseudomonas fluorescens 6-8 under gnotobiotic conditions. Can. J. Microbiol. 58, 170–178. doi: 10.1139/w11-120
Peiffer, J. A., and Ley, R. E. (2013). Exploring the maize rhizosphere microbiome in the field: a glimpse into a highly complex system. Commun. Integr. Biol. 6:e25177. doi: 10.4161/cib.25177
Peiffer, J. A., Spor, A., Koren, O., Jin, Z., Tringe, S. G., Dangl, J. L., et al. (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
Purkayastha, G. D., Mangar, P., Saha, A., and Saha, D. (2018). Evaluation of the biocontrol efficacy of a Serratia marcescens strain indigenous to tea rhizosphere for the management of root rot disease in tea. PLoS One 13:e0191761. doi: 10.1371/journal.pone.0191761
R Core Team (2019). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.
Rath, K. (2017). Soil Salinity as a Driver of Microbial Community Structure and Functioning. Available at: http://lup.lub.lu.se/search/ws/files/35616924/thesis_Dec12.pdf (accessed August 23, 2019).
Rilling, J. I., Acuña, J. J., Sadowsky, M. J., and Jorquera, M. A. (2018). Putative nitrogen-fixing bacteria associated with the rhizosphere and root endosphere of wheat plants grown in an andisol from southern Chile. Front. Microbiol. 9:2710. doi: 10.3389/fmicb.2018.02710
Robinson, M. D., Mccarthy, D. J., and Smyth, G. K. (2010). edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinforma. Appl. 26, 139–140. doi: 10.1093/bioinformatics/btp616
Sachdev, D. P., Chaudhari, H. G., Kasture, V. M., Dhavale, D. D., and Chopade, B. A. (2009). Isolation and characterization of indole acetic acid (IAA) producing Klebsiella pneumoniae strains from rhizosphere of wheat (Triticum aestivum) and their effect on plant growth. Indian J. Exp. Biol. 47, 993–1000.
Sanyal, S. K., Mou, T. J., Chakrabarty, R. P., Hoque, S., Hossain, M. A., and Sultana, M. (2016). Diversity of arsenite oxidase gene and arsenotrophic bacteria in arsenic affected Bangladesh soils. AMB Express 6:21. doi: 10.1186/s13568-016-0193-190
Schlaeppi, K., Dombrowski, N., Oter, R. G., Ver Loren van Themaat, E., and Schulze-Lefert, P. (2013). Quantitative divergence of the bacterial root microbiota in Arabidopsis thaliana relatives. Proc. Natl. Acad. Sci. U.S.A. 111, 585–592. doi: 10.1073/pnas.1321597111
Schmalenberger, A., Hodge, S., Bryant, A., Hawkesford, M. J., Singh, B. K., and Kertesz, M. A. (2008). The role of Variovorax and other Comamonadaceae in sulfur transformations by microbial wheat rhizosphere communities exposed to different sulfur fertilization regimes. Environ. Microbiol. 10, 1486–1500. doi: 10.1111/j.1462-2920.2007.01564.x
Schmitz-Esser, S., Toenshoff, E. R., Haider, S., Heinz, E., Hoenninger, V. M., Wagner, M., et al. (2008). Diversity of bacterial endosymbionts of environmental Acanthamoeba isolates. Appl. Environ. Microbiol. 74, 5822–5831. doi: 10.1128/AEM.01093-1098
Shade, A., and Handelsman, J. (2012). Beyond the venn diagram: the hunt for a core microbiome. Environ. Microbiol. 14, 4–12. doi: 10.1111/j.1462-2920.2011.02585.x
Shahid, M., Akram, M. S., Khan, M. A., Zubair, M., Shah, S. M., Ismail, M., et al. (2018). A phytobeneficial strain Planomicrobium sp. MSSA-10 triggered oxidative stress responsive mechanisms and regulated the growth of pea plants under induced saline environment. J. Appl. Microbiol. 124, 1566–1579. doi: 10.1111/jam.13732
Siciliano, S. D., and Germida, J. J. (1999). Taxonomic diversity of bacteria associated with the roots of field-grown transgenic Brassica napus cv. Quest, compared to the non-transgenic B. napus cv. Excel and B. rapa cv. Parkland. FEMS Microbiol. Ecol. 29, 263–272.
Smets, W., Leff, J. W., Bradford, M. A., Mcculley, R. L., Lebeer, S., and Fierer, N. (2016). A method for simultaneous measurement of soil bacterial abundances and community composition via 16S rRNA gene sequencing. Soil Biol. Biochem. 96, 145–151. doi: 10.1016/j.soilbio.2016.02.003
Ssekagiri, A. T., Sloan, W., and Zeeshan Ijaz, U. (2017). microbiomeSeq: An R Package for Analysis of Microbial Communities in an Environmental Context. Entebbe: ISCB Africa ASBCB conference.
Tkacz, A., Cheema, J., Chandra, G., Grant, A., and Poole, P. S. (2015). Stability and succession of the rhizosphere microbiota depends upon plant type and soil composition. ISME J. 9, 2349–2359. doi: 10.1038/ismej.2015.41
Vazquez, P., Holguin, G., Puente, M. E., Lopez-Cortes, A., and Bashan, Y. (2000). Phosphate-solubilizing microorganisms associated with the rhizosphere of mangroves in a semiarid coastal lagoon. Biol. Fertil. Soils 30, 460–468. doi: 10.1007/s003740050024
Visioli, G., Sanangelantoni, A. M., Vamerali, T., Dal Cortivo, C., and Blandino, M. (2018). 16S rDNA profiling to reveal the influence of seed-applied biostimulants on the rhizosphere of young maize plants. Molecules 23:E1461. doi: 10.3390/molecules23061461
Wilsey, B. J., and Potvin, C. (2000). Biodiversity and ecosystem functioning: importance of species evenness in an old field. Ecology 81, 887–892. doi: 10.1890/0012-9658(2000)081%5B0887:baefio%5D2.0.co;2
Wolf, A., Fritze, A., Hagemann, M., and Berg, G. (2002). Stenotrophomonas rhizophila sp. nov., a novel plant-associated bacterium with antifungal properties. Int. J. Syst. Evol. Microbiol. 52, 1937–1944. doi: 10.1099/ijs.0.02135-2130
Xiao, X., Fan, M., Wang, E., Chen, W., and Wei, G. (2017). Interactions of plant growth-promoting rhizobacteria and soil factors in two leguminous plants. Appl. Microbiol. Biotechnol. 101, 8485–8497. doi: 10.1007/s00253-017-8550-8558
Yadav, A. N., Sachan, S. G., Verma, P., and Saxena, A. K. (2015). Prospecting cold deserts of north western Himalayas for microbial diversity and plant growth promoting attributes. J. Biosci. Bioeng. 119, 683–693. doi: 10.1016/j.jbiosc.2014.11.006
Yang, E., Sun, L., Ding, X., Sun, D., Liu, J., and Wang, W. (2019). Complete genome sequence of Caulobacter flavus RHGG3 T, a type species of the genus Caulobacter with plant growth-promoting traits and heavy metal resistance. 3 Biotech 9:42. doi: 10.1007/s13205-019-1569-z
Yang, Y., Wang, N., Guo, X., Zhang, Y., and Ye, B. (2017). Comparative analysis of bacterial community structure in the rhizosphere of maize by highthroughput pyrosequencing. PLoS One 12:e0178425. doi: 10.1371/journal.pone.0178425
Ye, J. (2016). Bacterial Diversity and Functions in Soils Amended with Biochar and Organic Fertilizers. Available at: http://unsworks.unsw.edu.au/fapi/datastream/unsworks:42321/SOURCE02?view=true (accessed August 23, 2019).
Yoneyama, T., Terakado-Tonooka, J., and Minamisawa, K. (2017). Exploration of bacterial N2-fixation systems in association with soil-grown sugarcane, sweet potato, and paddy rice: a review and synthesis. Soil Sci. Plant Nutr. 63, 578–590. doi: 10.1080/00380768.2017.1407625
Zahir, I., Houari, A., Iraqui, M., and Ibnsouda, S. (2011). Aerococcus sp. with an antimycobacterial effect. Afr. J. Biotechnol. 10, 19473–19480. doi: 10.5897/AJB11.2473
Keywords: Brassica napus, breeding, canola, core microbiome, differential abundance, microbiome, plant–microbial interactions, rhizosphere
Citation: Taye ZM, Helgason BL, Bell JK, Norris CE, Vail S, Robinson SJ, Parkin IAP, Arcand M, Mamet S, Links MG, Dowhy T, Siciliano S and Lamb EG (2020) Core and Differentially Abundant Bacterial Taxa in the Rhizosphere of Field Grown Brassica napus Genotypes: Implications for Canola Breeding. Front. Microbiol. 10:3007. doi: 10.3389/fmicb.2019.03007
Received: 15 June 2019; Accepted: 13 December 2019;
Published: 15 January 2020.
Edited by:
Aymé Spor, INRA, UMR 1347 Agroécologie, FranceReviewed by:
Kim Yrjälä, Zhejiang Agriculture and Forestry University, ChinaClaudia Bartoli, Institut National de la Recherche Agronomique (INRA), France
Copyright © 2020 Taye, Helgason, Bell, Norris, Vail, Robinson, Parkin, Arcand, Mamet, Links, Dowhy, Siciliano and Lamb. 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: Zelalem M. Taye, zmt059@mail.usask.ca