Molecular Barcoding Reveals the Genus Streptomyces as Associated Root Endophytes of Apple (Malus domestica) Plants Grown in Soils Affected by Apple Replant Disease

Apple replant disease (ARD) is a complex phenomenon which affects apple tree nurseries and orchards worldwide, causing growth reductions and losses in fruit yield and quality (Manici et al. 2013; Mazzola und Manici 2012; Winkelmann et al. 2019). ARD occurs when apple or a closely related species is repeatedly planted at the same site and is described as a “harmfully disturbed physiological and morphological reaction of apple plants to soils that faced alterations in their (micro-)biome due to previous apple cultures” (Winkelmann et al. 2019). This disease is species specific and can persist for decades (Savory 1966). Because disinfection of the soil leads to better growth, it is generally accepted that biotic factors are the primary cause (Mahnkopp et al. 2018; Mai and Abawi 1981; Yim et al. 2013). Next to fungi belonging to the genera Fusarium, Cylindrocarpon, and Rhizoctonia, a number of other taxa, including oomycetes such as Pythium and Phytophthora spp., nematodes such as Pratylenchus spp., and various bacterial species such as members of the genera Pseudomonas and Bacillus as well as the phylum Actinobacteria, have been reported to contribute to ARD (Čatská et al. 1982; Manici et al. 2017; Mazzola 1998; Otto and Winkler 1993; Tewoldemedhin et al. 2011; Utkhede and Li 1988). However, despite decades of research, the etiology of ARD is still poorly known.

Based on molecular barcoding approaches in the last decade, many studies confirmed not only changes in the abundance of specific pathogens in ARD-affected soils but also significant shifts in the overall structure of the microbiome of the bulk soil and the rhizosphere (Winkelmann et al. 2019). These microbiome shifts also affect major functional properties, including the potential to degrade aromatic compounds and functions for biocontrol (Radl et al. 2019).

However, despite their close interaction with host cells, surprisingly, there is still a lack of knowledge of how microbes colonizing the root interior (root endophytes) are affected by ARD. The interior of roots can be colonized mostly by rhizosphere microbiota (e.g., through cracks formed during lateral root emergence and at root tips) (Bulgarelli et al. 2013; Hardoim et al. 2008). Positive effects of endophytes on plants include direct or indirect provision of nutrients (Gaiero et al. 2013; White et al. 2019); production of plant hormones such as auxin, cytokinins, or gibberellins (Hardoim et al. 2015; Santoyo et al. 2016); increased tolerance against abiotic stress (Hardoim et al. 2015); and biocontrol due to competitive mechanisms or production of antimicrobial substances (Haas and Keel 2003). Yet there are also endophytes known for their negative effects on plant health. Some of these facultative pathogens can shift their lifestyle, depending on several factors such as host and endophyte development stage, plant defense reactions, or environmental conditions (Schulz and Boyle 2005). Rosenblueth and Martínez-Romero (2006) put forward the hypothesis of an equilibrium between endophytes and plants that, under certain conditions, gets unbalanced to the detriment of one of the partners.

Only a few studies have investigated the role of endophytes in ARD, focusing on potential fungal root pathogens. Manici et al. (2013) found the root endophytic Cylindrocarpon-like fungi (Ilyonectria and Thelonectria spp.) and Pythium spp. to be main causal agents of growth reduction in the rootstock M9 growing in ARD-affected soil. Cylindrocarpon spp. were also identified next to Rhizoctonia spp. as a pathogenic root endophyte by Kelderer et al. (2012) in row-planted (ARD-affected) and interrow-planted (control) apple trees. In addition, Fusarium solani and F. oxysporum were most abundant in roots but not considered pathogenic. Popp et al. (2019) isolated several fungal endophytes from ARD-affected apple roots and reinoculated them in a soil-free biotest. Cadophora, Calonectria, Dactylonectria, Ilyonectria, and Leptosphaeria spp. were reported to have negative effects on plant health. In contrast, studies on the effects of ARD on bacterial root endophytes are scarce. Thus far, only a targeted, cultivation-dependent approach has been published, where the focus was on the biocontrol properties of Actinobacteria isolates (mostly belonging to the genus Streptomyces) from the root interior of apple trees. However, no effect was observed when coinoculated to apple seedlings with Pythium irregulare and Cylindrocarpon macrodidymum (Tewoldemedhin et al. 2011).

Therefore, the aim of this study was to investigate the bacterial root endophytic community structure in plants growing in ARD-affected soils compared with ARD-unaffected soils based on greenhouse biotests using a molecular barcoding approach. In these biotests, the ARD-susceptible genotype M26 was planted into ARD-affected soil (untreated or γ-sterilized) and grass control soil (untreated or γ-sterilized). To go beyond a local response pattern, we used soils from three different sites from northern Germany in the frame of this study. For generating more robust data, we performed our study in two subsequent years to exclude specific effects of the used soils based on one particular season. We propose that a possible causal agent of ARD should be present in roots from all ARD-affected soils in all three sites.

MATERIALS AND METHODS

Experimental setting.

Soil for these experiments was sampled from three different sites: Heidgraben (x-coordinate 53.699199; y-coordinate 9.683171; WGS 84, Schleswig-Holstein, northern Germany), Ellerhoop (x-coordinate 53.71435; y-coordinate 9.770143; WGS 84, Schleswig-Holstein, northern Germany), and Ruthe (x-coordinate 52.243668; y-coordinate 9.819700; WGS 84, Lower Saxony, Germany). These sites differed in their climatic conditions and soil properties. The upper soil textures of the three sites were defined (based on World Reference Base for soil resources) as sand (Heidgraben), loamy sand (Ellerhoop), and silt loam (Ruthe) (Mahnkopp et al. 2018). Every site contained two different plot variants: (i) ARD plots, where ARD was successfully induced by repeatedly replanting Bittenfelder apple seedlings since 2009 in a 2-year cycle, and (ii) control plots, which were only covered with grass since then. ARD plots in Ruthe and Ellerhoop were replanted for the last time in spring 2015 and in Heidgraben in spring 2016.

Soils were sampled from all three sites at a depth of 0 to 20 cm at the end of 2015 and 2016. After sampling and sieving (8 mm), soils were either γ irradiated (G) at a minimal dose of 10 kGy or left untreated (UT), resulting in four variants per site: ARD untreated (ARD UT), ARD γ-irradiated (ARD G), grass untreated (grass UT), and grass γ-irradiated (grass G).

The ARD-susceptible apple rootstock M26, which was propagated and rooted in vitro (Weiß et al. 2017a), was acclimatized for 4 weeks and, afterwards, one plantlet each was planted in 1-liter pots containing the different soil variants. Soils were supplemented with Osmocote Exact 3-4M (16 + 9 + 12 + 2 MgO; https://icl-sf.com/de-de/products/ornamental_horticulture/8840-osmocote-exact-standard-3-4m/) at 2 g liter⁻¹ to exclude nutrient effects. Shoot lengths were measured weekly. Plants were grown for 8 weeks in the greenhouse at a mean daily temperature of 21°C and a 16-h photoperiod achieved by additional light (SON-T Philips Master Agro 400 W) with a set point of 25 klx. Plant protection measures and irrigation were done according to Yim et al. (2015). At the first sign of insect pests, 0.3% NeemAzal was sprayed. During the night, a sulfur evaporator was used in order to prevent fungal diseases. The greenhouse experiment was conducted twice, in February 2016 and 2017 (Mahnkopp et al. 2018), with nine replicates per variant.

Acclimatized plants (before planting into the soil variants) were treated as described by Mahnkopp et al. (2018) and served as the source for “timepoint zero” (T0) samples in both years.

Sampling.

After 8 weeks of cultivation in the greenhouse, 4 representative plants per variant were taken (48 per year, 96 in total) as biological replicates. Roots were washed carefully to get rid of the adhering soil. Shoot and root fresh mass were determined. For surface sterilization, roots were rinsed for 30 s in EtOH (70%), followed by stirring in 2% NaOCl for 7.5 min and, finally, washing five times in sterile deionized water. The final washing water was plated on 523 medium (Viss et al. 1991) and incubated at room temperature for 1 week. Plating resulted in <10 CFU/plate in all cases. Roots were stored in sterile 2-ml Eppendorf tubes at −80°C until DNA extraction.

DNA extraction and amplicon sequencing.

Surface-sterilized roots (50 to 100 mg/sample) were homogenized under frozen conditions using steel beads (Ø 6 mm) in a mixer mill (MM400; Retsch, Haan, Germany) with a frequency of 23 Hz for 2 min using sterilized devices. DNA was extracted using the Invisorb Spin Plant Mini Kit (Stratec, Berlin, Germany) according to the provided protocol. DNA quality was checked using a spectrophotometer (Nanodrop 2000c; Peqlab, Erlangen, Germany).

The primer combination 335F (CADACTCCTACGGGAGGC)/769R (ATCCTGTTTGMTMCCCVCRC) (Dorn-In et al. 2015), including overhang adapter sequence, was used to amplify the V3-V4 region of the 16S ribosomal RNA (rRNA) gene. PCR assays contained 2× Phusion High-Fidelity Master Mix (1.5 mM MgCl₂, 200 μM each dNTP, and 0.2 U of Phusion DNA Polymerase [Thermo Fisher Scientific, Waltham, MA, U.S.A.]), 10 pmol of each primer, 5 ng of DNA template, and water to a final volume of 10 µl. The PCR cycling conditions consisted of an initial denaturation step of 98°C for 10 s; followed by 30 cycles involving 1 s of denaturation at 98°C, 5 s of annealing at 59°C, and 45 s of extension at 72°C; with a final extension of 1 min at 72°C. Triplicate PCR assays were pooled and purified using Agencourt AMPure XP kit (Beckman Coulter, Brea, CA, U.S.A.). The purified products were quantitated using the Quant-IT PicoGreen dsDNA assay kit (Life Technologies Europe, Gent, Belgium). Sample indexing was carried out with Nextera XT Index Kit v2 Set A and B (Illumina, San Diego, CA, U.S.A.) in reaction mixtures containing 10 ng of purified PCR product, 2× Phusion High-Fidelity Master Mix (1.5 mM MgCl₂, 200 μM each dNTP, and 0.2 U of Phusion DNA Polymerase [Thermo Fisher Scientific]), 10 pmol of each indexing primer, and water to a final volume of 25 µl. The indexing PCR cycling conditions consisted of an initial denaturation step of 98°C for 30 s; followed by 8 cycles involving 10 s of denaturation at 98°C, 30 s of annealing at 55°C, and 30 s of extension at 72°C; with a final extension of 5 min at 72°C.

Indexed samples were purified as described above. Equimolar concentrations of the purified indexed samples were prepared and diluted to a final concentration of 4 nM. The library was sequenced using the Illumina MiSeq platform with the MiSeq Reagent Kit v3 (600 cycle) (Illumina).

Bioinformatic and statistical analysis.

FASTQ files were trimmed with a minimum read length of 50 and a minimum Phred score of 15 using AdapterRemoval (Schubert et al. 2016) without merging forward and reverse reads. Afterward, sequences were analyzed using the QIIME 2 software package release 2017.11 (Caporaso et al. 2010) with default parameters. The QIIME 2 plugin DADA2 (Callahan et al. 2016) was used for quality control with the following parameters: 10 bp were removed n-terminally, and reads were truncated at position 300 (forward) and 260 (reverse) for universal 16S rRNA genes. Expected error was adjusted to 2.

Taxonomic analysis of the resulting unique amplicon sequence variants (ASVs) was performed using primer-specific pretrained Naive Bayes classifiers of the SILVA_132_QIIME release 99% and the q2-feature-classifier plugin, setting the confidence threshold to 0.9. Because the PCR-negative control showed no ASVs, contamination during sample processing could be excluded. For further data analysis, unassigned reads and singletons (in sum, <0.03% of all reads) were excluded.

Raw sequence data were deposited in GenBank (https://www.ncbi.nlm.nih.gov/genbank/) under the accession PRJNA647245.

16S amplicon sequencing of DNA extracted from the surface-sterilized M26 roots resulted after quality control in a total of 4,132,410 reads with a mean of 72,498 reads/sample in the biotest in 2016. After removal of chloroplast, archaea, and eukaryotic ASVs, sequence data were rarefied at a number of 28,817 reads (2016). In total, 5,898 ASVs were detected. Because rarefaction analysis (Supplementary Fig. S1) indicated that saturation was reached already with read numbers approximately 5,000 reads/sample, for 2017, the sequencing effort was reduced, resulting in in total 884,916 reads and a mean of 17,351 reads/sample (rarefied to 4,213 reads/sample after the removal of the chloroplasts). The total number of ASVs was 4,971, which were nearly all covered after rarefying at 4,813 reads/sample (Supplementary Fig. S2). In order to identify identical ASVs between the 2 years, alignments of the sequences of the 2 years on genus level were done using ClustalW Multiple Alignment (Thompson et al. 1994), with number of bootstraps set to 1,000 using BioEdit v7.2.5 (Hall 1999) followed by calculating a sequence differences count matrix. ASVs from 2017 which were 100% identical to ASVs from 2016 were given the corresponding name of 2016 ASVs to improve comparability of figures and tables.

To calculate the relative abundance, the number of reads per ASV in the samples was divided by the sum of total reads per sample and multiplied by 100. The relative abundances of ASVs belonging to the same phylum or genus were combined to calculate the overall relative abundance of the corresponding phylum or genus. Species diversity (Shannon and Simpson) and richness (Chao1) indices were determined using the “Phyloseq” (McMurdie and Holmes 2013) and “Vegan” (Oksanen et al. 2019) packages of R v3.6.1 (R Development Core Team 2019) (http://www.R-project.org) and tested for normal distribution based on the Shapiro-Wilk test (Shapiro and Wilk 1965) and homogeneity of variance based on Levene’s test (Levene 1960) using the program PAST3 v. 3.20 (Hammer et al. 2001). If the null hypotheses of normal distribution and equal variances were rejected, the Tukey test based on Herberich et al. (2010) was used at P < 0.05 to determine significant differences of the raw diversity and richness scores. In order to compare the relative abundance of different genera of the initial microbiome (T0 plants) between the years, all ASVs belonging to the same genus were merged. Because of unequal sample size and unequal variance, Welch’s two-sample t test was used at P < 0.05 to determine significant differences. Nonmetric multidimensional scaling (NMDS) was performed with the program PAST3 v. 3.20 (Hammer et al. 2001) using the Bray-Curtis similarity index and analysis of similarity (ANOSIM) in order to visualize the community composition of the different samples. To indicate the influence of the different genera, vectors were added which show the correlation between the corresponding genus and the NMDS score. Spearman’s correlation was used in order to correlate ASVs to shoot growth and fresh mass using the program PAST3 v. 3.20 (Hammer et al. 2001). Venn diagrams were designed using the Venn diagram tool of Bioinformatics & Evolutionary Genomics (http://bioinformatics.psb.ugent.be/webtools/Venn).

To identify specific Streptomyces ASVs, their nucleotide sequences (414 bp) were blasted against the NCBI database (https://www.ncbi.nlm.nih.gov/) using BLASTn.

RESULTS

Root and shoot biomass.

As expected, in both years, plants grown in untreated ARD-affected soils showed the lowest increase in shoot length and lowest shoot fresh mass in comparison with the other soil variants (Table 1). The γ irradiation of grass and ARD soil led to increased shoot length and shoot fresh mass compared with the respective untreated soils. Overall, plants grown in the biotest in 2017 showed a higher biomass in comparison with plants grown in 2016 but the response pattern to the different soils variants was comparable between both years. Detailed growth data for the full set of nine plants per variant can be found in Mahnkopp et al. (2018) whereas, in Table 1, only the data of the plants selected for the barcoding approach of this study are presented.

TABLE 1 Shoot fresh mass and increase of shoot length of M26 apple plants grown for 8 weeks in the greenhouse biotest in 2016 and 2017^y

	Heidgraben				Ellerhoop				Ruthe
	ARD		Grass		ARD		Grass		ARD		Grass
Yearz	UT	G	UT	G	UT	G	UT	G	UT	G	UT	G
2016
Length (cm)	5.9 ± 2.3 a	18.5 ± 1.5 c	13.5 ± 2 b	27.7 ± 6.3 c	7.5 ± 1 a	12.4 ± 3.8 ab	19.6 ± 1.6 b	25.9 ± 1.7 c	8.5 ± 1.4 a	25.2 ± 0.2 c	16.7 ± 1.9 b	32.0 ± 1.4 d
Mass (g)	2.8 ± 0.6 a	7.2 ± 0.6 b	5.9 ± 0.6 b	10.7 ± 1 c	3.2 ± 0.7 a	5.9 ± 1.3 ab	7.4 ± 1.1 b	14.4 ± 0.8 c	3.7 ± 0.4 a	11.4 ± 1 c	7.1 ± 0.1 b	17.6 ± 2.3 d
2017
Length (cm)	5.7 ± 2.2 a	31.7 ± 1.3 b	29.8 ± 3.4 b	42.3 ± 3.9 c	13.4 ± 3.4 a	24.7 ± 0.8 b	31.6 ± 3.8 bc	43.0 ± 4.9 c	22.9 ± 1.9 a	34.3 ± 1.5 b	36.3 ± 4.9 bc	41.6 ± 1.7 c
Mass (g)	2.8 ± 0.6 a	11.7 ± 0.7 bc	9.4 ± 1.2 b	19.8 ± 3.9 c	4.8 ± 0.9 a	9.4 ± 0.3 b	9.7 ± 1.9 b	18.7 ± 3.6 b	8.3 ± 1.1 a	12.8 ± 1.4 b	11.8 ± 2.1 ab	15.9 ± 1 b

^y Surface-sterilized roots of these plants were used for DNA extraction and amplicon sequencing. ARD = apple replant disease, UT = untreated, and G = γ irradiated. Shown is the mean and the standard deviation of n plants. Different letters indicate significant differences within the sites (Tukey’s test, P ≤ 0.05).

^z Year and parameter: Length = increase in shoot length and Mass = shoot fresh mass.

View as image HTML

Endophytic bacterial community composition and diversity in plant roots grown in different soil variants.

The mean number of observed ASVs per sample was 244 in 2016 (Table 2) and 201 in 2017 (Table 3). The highest number of observed ASVs were found in plants grown in ARD G soil from Ellerhoop (n = 339 ± 93) and the lowest number in the variant Heidgraben ARD G (n = 148 ± 62) (both in 2016). In both years, no significant differences in diversity or richness indices were recorded within or between the sites (Tukey’s test, P < 0.05).

TABLE 2 Richness and diversity of endophytic bacterial communities based on amplicon sequence variants (ASVs) in roots grown for 8 weeks in soils from different sites and treatments (ARD = apple replant disease, UT = untreated, and G = γ irradiated) of the biotest in 2016^z

Site, soil	Treatment	n	Observed ASVs	Chao1	Shannon	Simpson
T0		4	290 ± 70	292 ± 70	5.12 ± 0.24	0.99 ± 0.00
Heidgraben
ARD	UT	4	255 ± 106	257 ± 108	4.37 ± 0.50	0.97 ± 0.02
	G	4	148 ± 62	149 ± 61	3.73 ± 0.88	0.91 ± 0.08
Grass	UT	4	289 ± 119	294 ± 124	4.42 ± 0.50	0.97 ± 0.01
	G	4	329 ± 88	331 ± 89	4.88 ± 0.27	0.98 ± 0.01
Ellerhoop
ARD	UT	4	303 ± 88	306 ± 89	4.86 ± 0.35	0.98 ± 0.00
	G	4	339 ± 93	344 ± 95	4.66 ± 0.52	0.97 ± 0.02
Grass	UT	4	202 ± 38	204 ± 39	4.15 ± 0.65	0.96 ± 0.03
	G	4	225 ± 110	228 ± 111	3.95 ± 0.65	0.94 ± 0.04
Ruthe
ARD	UT	4	252 ± 118	264 ± 135	4.13 ± 0.72	0.94 ± 0.06
	G	4	211 ± 56	212 ± 56	4.33 ± 0.56	0.96 ± 0.02
Grass	UT	3	205 ± 18	205 ± 18	4.39 ± 0.13	0.97 ± 0.01
	G	4	167 ± 16	168 ± 17	4.00 ± 0.45	0.95 ± 0.02

^z Additionally, acclimatized plants at timepoint zero (T0) before transferring into the soil variants are shown. There was no significant difference within and between the sites according to Tukey’s test at P ≤ 0.05. Shown are mean ± standard deviation of n replicates.

View as image HTML

TABLE 3 Richness and diversity of endophytic bacterial communities based on amplicon sequence variants (ASVs) in roots grown for 8 weeks in soils from different sites and treatments (ARD = apple replant disease, UT = untreated, and G = γ irradiated) of the biotest in 2017^z

Site, soil	Treatment	n	Observed ASVs	Chao1	Shannon	Simpson
T0		5	166 ± 42	170 ± 47	4.3 ± 0.30	0.97 ± 0.02
Heidgraben
ARD	UT	4	149 ± 35	152 ± 36	4.07 ± 0.30	0.96 ± 0.01
	G	4	235 ± 84	239 ± 88	4.96 ± 0.28	0.99 ± 0.00
Grass	UT	4	152 ± 21	156 ± 24	3.94 ± 0.24	0.96 ± 0.02
	G	3	154 ± 49	154 ± 50	4.23 ± 0.49	0.97 ± 0.02
Ellerhoop
ARD	UT	4	263 ± 50	275 ± 55	4.73 ± 0.35	0.98 ± 0.01
	G	2	206 ± 37	212 ± 33	4.84 ± 0.10	0.99 ± 0.00
Grass	UT	4	194 ± 60	200 ± 63	4.22 ± 0.42	0.94 ± 0.04
	G	3	165 ± 64	172 ± 64	4.26 ± 0.62	0.97 ± 0.02
Ruthe
ARD	UT	3	193 ± 6	202 ± 10	4.3 ± 0.11	0.97 ± 0.01
	G	4	200 ± 25	201 ± 24	4.59 ± 0.24	0.98 ± 0.01
Grass	UT	4	253 ± 105	263 ± 117	4.59 ± 0.33	0.97 ± 0.01
	G	4	246 ± 148	261 ± 161	4.32 ± 0.74	0.96 ± 0.02

^z Additionally, acclimatized plants at timepoint zero (T0) plants before transferring into the soil variants are shown. There was no significant difference within and between the sites according to Tukey’s test at P ≤ 0.05. Shown are mean ± standard deviation of n replicates.

View as image HTML

Proteobacteria was clearly the dominant phylum in all variants, with relative abundance ranging from 66.9% (Heidgraben ARD G) to 83.7% (Ruthe grass UT) (Supplementary Fig. S3) for the biotest performed in 2016. Phylum Bacteroidetes was of second most abundance, with a mean value of 10.3%, followed by Actinobacteria (5%) and Firmicutes (4.4%). Members of Actinobacteria appeared in higher relative abundance in roots grown in ARD UT compared with the other variants of the respective site. The strongest difference was observed in Heidgraben, where the abundance of ASVs linked to Actinobacteria in ARD UT variants (15.60%) was significantly higher than in ARD G (5.03%), grass UT (2.75%), and grass G variants (3.47%). The second biotest in 2017 showed similar shares for the different phyla (Supplementary Fig. S4). ASVs assigned to Proteobacteria ranging from 66.7 to 89.3% were dominant, followed by Bacteroidetes, with a mean value of 13.6% relative abundance; Actinobacteria, with 8.1%; and Firmicutes, with 2.2%. However, the higher abundance of ASVs linked to phylum Actinobacteria in the untreated ARD variants was not observed in this year.

Dynamics of endophytic bacteria during the biotest.

T0 plants of 2016 had the highest diversity of all treatments over the years (Table 2). In 2017, T0 plants showed a significantly lower diversity compared with 2016 T0 plants (Supplementary Fig. S8). After cultivation for 8 weeks in the different soils, 9 of 12 variants of 2016 still had higher numbers in the observed ASVs compared with 2017. However, these differences were not significant (Supplementary Fig. S8).

On the phylum level, T0 plants grown in 2016 were dominated by Proteobacteria (79%, Supplementary Fig. S3), which did not change after 8 weeks of cultivation in the different soils. In 2017, this value increased to 89% for T0 plants but, here, this high relative abundance was found reduced by approximately 16% after the plants had been grown for 8 weeks in the different soil variants, regardless of the soil (Supplementary Fig. S4).

On higher phylogenetic levels, proteobacterial groups of the genera Shewanella and Halomonas belonged to the top three genera regarding relative abundance in T0 plants in both years (Supplementary Table S2). However, in total, nearly one-third (31.9%) of the genera showed significant differences in abundance between the years. In T0 plants from 2016, for example, Ralstonia was (with 6.6%) the most abundant genus but was not present in 2017 T0 plants. Even after growing for 8 weeks in different soils, this abundance pattern still remained for ASVs linked to Ralstonia. Similar contrasting abundance patterns were observed for genus Pseudomonas when both years were compared. Here, we could link 12.6% of all ASVs from T0 plants to this genus in 2017, which was more than four times higher than in 2016. However, in contrast to ASVs linked to Ralstonia, after 8 weeks of cultivation in the different soils, these initial differences in abundance of Pseudomonas were no longer detectable.

To analyze β-diversity, three-dimensional nonmetric multidimensional scalings were created for the untreated soils. In general, high variability within variants could be observed. In both years, T0 plants significantly separated from the other variants (ANOSIM with P ≤ 0.05) (Fig. 1; Supplementary Fig. S5). Only for the biotest in 2016, significant differences between other variants were observed, especially for the treatments with soil from Ruthe (Fig. 1). ASVs assigned to genera Pseudomonas, Rhizobium, and especially Streptomyces were closely linked with ARD, whereas ASVs related to Rhodanobacter, Dyella, Bradyrhizobium, Sphingomonas, and Rhizomicrobium pointed to the untreated grass variants (Fig. 1). Most responsible for differentiation of T0 were ASVs which were linked to genera Halomonas, Acinetobacter, and Shewanella. In the biotest in 2017, no clear clustering except for T0 was observed (Supplementary Fig. S5).

Download as PowerPoint

Identification of bacterial responders in the different treatments and correlation with plant growth.

To further investigate ASVs responding to the different treatments of each site, Venn diagrams were designed. On the one hand, the number of ASVs shared by all four different variants per site, which we considered as the core microbiome of a given site, was surprisingly small. In 2016, only six ASVs (relative abundance >0.5%) in soil variants from Heidgraben, three in those from Ellerhoop, and four in those from Ruthe were present in all variants (Fig. 2). In 2017, these numbers were reduced to zero (Heidgraben), two (Ellerhoop), and one (Ruthe) (Fig. 3). On the other hand, the number of ASVs which were unique for each variant was high. In 2016, unique ASVs in untreated ARD variants ranged from 15 for Ruthe and 19 for Heidgraben to 28 for Ellerhoop (Fig. 2). In the untreated grass variants, the number of unique ASVs ranged from 23 for Ruthe to 26 for Ellerhoop and 29 for Heidgraben. Although the overall distribution was very similar in both years, for the soil from Ellerhoop, some variations were observed: the number of unique ASVs changed in untreated grass variants from 26 in 2016 to 9 in 2017. Results for soil variants from Heidgraben showed the lowest variation between the years, except for the unique ASVs for the grass variant sterilized by γ-irradiation, where 16 (2016) and 31 (2017) unique ASVs were observed.

Download as PowerPoint

Download as PowerPoint

In order to identify responders toward ARD, the unique ASVs of the untreated ARD variants were correlated with shoot growth and fresh mass of all variants of the three sites. In 2016, most noticeable were ASVs related to the genus Streptomyces (Fig. 2), which closely linked to plants grown in ARD-affected soils, confirming the overall observation that ASVs related to Actinobacteria were positively responding to the ARD-affected soils with increased levels in relative abundance. In Heidgraben ARD UT, 7 of 19 unique ASVs were linked to the genus Streptomyces, followed by Ellerhoop ARD UT (4 of 28) and Ruthe ARD UT (3 of 15). Most of these ASVs were high in relative abundance. Streptomyces ASV66 in Heidgraben and Streptomyces ASV42 in Ruthe showed the highest relative abundance, with 4.48 and 4.10%, respectively. All Streptomyces ASVs were negatively correlated with the increase of shoot length and shoot fresh mass and some of them were even present in at least two sites. Streptomyces ASV21, which was present in all three sites as a unique ASV, showed the second highest negative correlation with both plant growth parameters (−0.54 and −0.58, respectively). This number was only surpassed by Streptomyces ASV70 and Streptomyces ASV76 (both present in Heidgraben and Ellerhoop), with a correlation of −0.59 to the increase of shoot length and −0.65 to shoot fresh mass.

A high number of other genera harboring unique ASVs were also negatively correlated with plant growth parameters. For example, Novosphingobium ASV92 and Neorhizobium ASV47 (highly abundant in Heidgraben and Ellerhoop) negatively correlated with increase of shoot length (−0.59 and −0.53, respectively) and shoot fresh mass (−0.53 and −0.59, respectively) (Fig. 2).

In 2017, the overall number of negatively correlated ASVs was lower (Fig. 3). Nevertheless, several Streptomyces ASVs (also present in at least two sites) again were negatively correlated with plant growth. In Heidgraben and Ellerhoop, the relative abundance of Streptomyces ASV76 and Streptomyces ASV621 showed a correlation of −0.53 and −0.57 to increase of shoot length and −0.60 and 0.61 to shoot fresh mass, respectively.

We were further interested in unique ASVs of the untreated grass variants to identify possible plant-growth-promoting bacteria (PGPB), which could help to counteract ARD. In 2016, several ASVs of different genera were present in more than two sites (Supplementary Fig. S6). These included ASVs related to genera Dyella, Massilia, Rhizobium, Rhodanobacter, and unclassified species of family Moraxellaceae. In the biotest 2017, only two ASVs (assigned to genera Rhizobium and Sphingobium) were present at at least at two sites (Supplementary Fig. S7) However, none of them showed positive correlations with plant growth. The only positively correlated ASV was related to unclassified members of family Rhizobiaceae and found in 2017 with a relative abundance of 0.95% and a correlation of 0.32 and 0.30 to shoot growth and shoot fresh mass, respectively (Supplementary Fig. S7).

DISCUSSION

In the present study, we characterized the bacterial root endophytic community of apple plants grown in replant and nonreplant soil in order to understand the etiology of ARD and develop countermeasures.

Endophytic bacterial communities in apple roots were dominated by Proteobacteria.

Proteobacteria was the dominant phylum in most studies where bulk soil or rhizosphere samples from ARD-affected sites had been analyzed (Franke-Whittle et al. 2015; Peruzzi et al. 2017; Sun et al. 2014; Tilston et al. 2018; Yim et al. 2015), with an average relative abundance of 35% (Nicola et al. 2018). The same was true for the root endophytes analyzed in our experiments in both years (Supplementary Figs. S3 and S4). However, in comparison with Nicola et al. (2018), the relative abundance of Proteobacteria in roots from plants grown in ARD UT was clearly higher (76% in 2016 and 71% in 2017 in an average of all three sites). This enrichment of Proteobacteria in the endosphere could be explained by selective recruitment or colonization or a higher competitiveness inside the plant. Members of this phylum are known for their various secretion systems (Preston et al. 2005), their fast growth, and their high metabolic activity and, therefore, they mostly predominate the endosphere (Lundberg et al. 2012; Reinhold-Hurek et al. 2015).

In many studies where rhizosphere and bulk soil samples of ARD-affected sites have been analyzed, Actinobacteria was among the most abundant phyla, with an average relative abundance of 16% (Nicola et al. 2018), which was slightly higher than in our study when we focused on root endophytes (10 and 11% in 2016 and 2017, respectively). In other studies, Bacteroidetes showed an average relative abundance of 14% (Nicola et al. 2018) and 13% (Tilston et al. 2018). In our study, values ranged between 8% (2016) and 15% (2017). As expected, members of Acidobacteria, which were also highly abundant in the rhizosphere and bulk soil in the abovementioned studies, were low in relative abundance in the root interior, due to the ecophysiological properties of these bacteria, including the use of complex organic compounds and their slow growth.

Differences between the outcome of the biotest of 2016 and 2017.

Significant differences occurred in the results comparing biotests between 2016 and 2017. In 2016, Actinobacteria were significantly higher in relative abundance in root samples from ARD UT compared with the grass or γ-sterilized variants (Supplementary Fig. S3). Surprisingly, this was not the case in 2017 (Supplementary Fig. S4). Furthermore, there was a clustering of ARD variants apart from the grass variants in 2016 (Fig. 1) but not in 2017 (Supplementary Fig. S5). These differences in the outcome of the biotests could be due to various factors (e.g., higher shoot lengths of T0 plants of 2017 or variation in environmental factors). Another reason could be related to the soil. First of all, soil collection might have resulted in samples of different microbial composition due to patchy appearance of ARD in the field (Simon et al. 2020). Furthermore, at our reference sites, replanting takes place every second year, and was carried out in 2015 and 2017 at Ruthe and Ellerhoop and in 2014 and 2016 at Heidgraben. Soil for the first greenhouse experiment in 2016 was sampled at the end of 2015, where plants at Ruthe and Ellerhoop had been replanted for the fourth time in spring 2015, whereas plants in Heidgraben had just been uprooted. For the experiment in 2017, soil was collected at the end of 2016, when plants at Ellerhoop and Ruthe had been uprooted and at Heidgraben had been replanted for the fifth time in spring. It is known that the microbial community composition in the rhizosphere of different apple genotypes varies seasonally and among years (Rumberger et al. 2007). Also, replanting is known to have an influence of the rhizosphere community composition (Sun et al. 2014).

Another reason for these different outcomes of the biotest in the 2 years could be differences in the initial endophytic bacterial community composition of the starting material (T0) plants. In this study, a higher diversity and number of ASVs of the initial microbiome in the plant roots was found in 2016 compared with that in 2017. However, this difference between the 2 years was not significant after 8 weeks of growth in the soil. Therefore, a higher diversity and number of ASVs of the initial microbiome did not lead to higher number of ASVs and diversity in plant roots after 8 weeks of growth in the soil. Because soil is the main reservoir of microorganisms for the plant microbiome (Berg and Smalla 2009; Bonito et al. 2014; Hartman and Tringe 2019; Lareen et al. 2016), it is one of the major factors influencing the number of ASVs and diversity. However, some genera (e.g., Ralstonia) were present in the 2017 T0 plants and were still present in the plant roots after growth for 8 weeks in the soil (Supplementary Table S2). With this in mind, one strategy to help to overcome ARD could be to inoculate apple plants with PGPB before transferring them into the soil. This so-called microbiome engineering of plants was recently reviewed by Orozco-Mosqueda et al. (2018). Johnston-Monje and Raizada (2011) could show that green fluorescent protein-tagged Enterobacter asburiaes could systemically colonize the roots of maize and even the rhizosphere. This means that some genera of the initial endophytic microbiome may possess the ability to colonize not only the plant roots but also the rhizosphere. In order to have plants which are preinoculated with PGPB that can influence not only the endosphere but also the rhizosphere and, therefore, may be able to reduce the ARD effect, inoculation studies are needed.

Next to differences between the years, variations within the soil variants were observed. The observed ASVs and the different diversity indices within the variants showed high standard deviations (Tables 2 and 3). With the reanalysis of several studies of microbiomes of ARD-affected soils, Nicola et al. (2018) determined that the strongest factor for bacterial community variation were environmental variables. In our study, several factors responsible for variations were reduced to a minimum (soils mixed, clonally propagated plants, and the same greenhouse conditions). However, here, we were investigating the root endophytic community. Its selection is strongly controlled by the host plant and dependent on soil and several other factors such as stress and environmental conditions (Afzal et al. 2019). Although all plants and soils within a variant were treated in the same way, individual differences in the soil microbiome and, therefore, differences in root colonization cannot be excluded. To reduce these variations, future biotests should increase the number of analyzed plants.

Are streptomycetes part of the ARD complex?

Our results show that nearly all Streptomyces ASVs were negatively correlated with increase in shoot length and shoot fresh mass (Figs. 2 and 3). Interestingly, the same ASVs were also identified in the roots of the rootstock cultivar Bittenfelder grown in ARD-affected soil in the three reference field sites Heidgraben, Ellerhoop, and Ruthe (results not shown). Therefore, regardless of the year, whether greenhouse biotest or field experiment, the site, or the apple rootstock genotype, Streptomyces ASVs were associated with apple roots grown in ARD soils. This raises the question whether Streptomyces spp. are a causative part of the ARD complex or just opportunistic or secondary colonizers.

Streptomyces is a well-studied genus and most famous for its production of antibiotics, with 80% of today’s antibiotics being derived from Streptomyces spp. (de Lima Procópio et al. 2012). Next to this, traits such as production of antifungal substances and siderophores, solubilization of phosphate, synthesis of plant growth regulators, secretion of volatile compounds, biocontrol (competition for nutrients), and degradation of phytotoxins makes this genus a potent PGPB, intensively reviewed by Olanrewaju and Babalola (2019), Sousa and Olivares (2016), Viaene et al. (2016), and Vurukonda et al. (2018). These reviews also highlight that genus Streptomyces is able to colonize a broad range of plant hosts. It is further believed that these plants can selectively recruit Streptomyces spp. (Viaene et al. 2016). However, the signals which attract them or the way of their entering and colonizing the roots are still unknown (Viaene et al. 2016; Vurukonda et al. 2018).

However, these various plant-growth-promoting effects of Streptomyces were not affirmed by our findings. Roots growing in soil from the grass variants showed better growth than those in ARD soils. Only two Streptomyces ASVs unique for at least two grass variants were found which had no correlation with increase in shoot length and shoot fresh mass (Supplementary Figs. S6 and S7). In contrast, in ARD variants, a clear negative correlation of the relative abundance of Streptomyces spp. and plant growth was shown (Figs. 2 and 3). On the one hand, this could indicate that genus Streptomyces is pathogenic and part of the replant disease. On the other hand, because Streptomyces has a saprophytic lifestyle, it could be an opportunist and degrade dead or damaged root material. Structurally damaged and partially necrotic root systems are typical symptoms for ARD-affected plants (Grunewaldt-Stöcker et al. 2019). Streptomyces spp. are able to break down organic remains of plants using several hydrolytic exoenzymes such as cellulases, lignocellulases, pectinases, xylanases, and cutinases (Chater 2016; Chater et al. 2010). Streptomyces was also shown to appear in higher abundance in the rhizosphere of Arabidopsis thaliana when plant exudated phenolic-related compounds such as salicylic acid were present (Badri et al. 2013; Lebeis et al. 2015) and can even grow on minimal media with only salicylic acid as a carbon source (Lebeis et al. 2015). Due to tyrosinase activity, some isolates were partially protected against plant-produced phenols, leading to increased colonization rates on A. thaliana roots (Chewning et al. 2019). Gene expression studies revealed that genes responsible for the production of phytoalexins (some of which belong to polyphenols) are upregulated in M26 roots growing in ARD-affected soils (Weiß et al. 2017a,b).

Overall, these reasons make it seem likely that genus Streptomyces finds favorable conditions and occurs in higher abundance in ARD-affected roots and, hence, is opportunistic. Yet pathogenicity cannot be excluded.

Of the 672 known Streptomyces spp. (Euzéby 1997; Parte 2018) (number as of 3 December 2020), only 10 have pathogenic features (Viaene et al. 2016). Most known are Streptomyces scabies, S. acidiscabies, and S. turgidiscabies, which cause common scab on roots and tuber crops. These species are able to directly penetrate plant cells and, in addition to necrotic scab lesions, lead to reduced growth, root stunting and browning, and a reduction of the complexity of the root system (Loria et al. 2003, 2006; Seipke et al. 2012) (i.e., symptoms that resemble the phenotype of ARD-affected roots). However, despite the large host range, none of these species were reported to infect woody plants, although the host range likely includes all higher plants, because dicot and monocot seedlings of several plant species have shown symptoms after inoculation with S. scabies (Leiner et al. 1996; Loria et al. 2006). One reason for this large host range is based on the assumption that genus Streptomyces is believed to have originated 400 million years ago, when green plants started to colonize the land (Chater 2016). Another reason for this flexibility is the fact that Streptomyces virulence genes are clustered on a pathogenic island which can be mobilized and, via conjugation, transferred to nonpathogenic relatives, which leads to the emergence of new plant-pathogenic streptomycetes (Lerat et al. 2009).

A closer look at the genus Streptomyces from our greenhouse experiment revealed that the Streptomyces ASVs which occurred in at least two sites (Figs. 2 and 3) shared a high similarity with the pathogen S. turgidiscabies. Blasting the sequences against the NCBI database (https://www.ncbi.nlm.nih.gov/) showed a similarity of 99.76% (ASV21 and ASV70), 99.51% (ASV76), and 99.52% (ASV621) (Supplementary Table S1). All of these ASVs showed a negative correlation with shoot fresh mass of approximately −0.60 or more, whereas ASV121, which showed the lowest negative correlation of −0.34 also shared the lowest identity with S. turgidiscabies (97.32%). However, for further comparisons to Streptomyces spp., the complete 16S rRNA sequence of the apple root endophytes identified in this study is necessary.

Nevertheless, the high similarity to pathogenic Streptomyces spp., the broad host range, and the ability for horizontal gene transfer of virulence genes may be arguments in favor of genus Streptomyces as a possible causative organism of ARD.

Role of genus Streptomyces in ARD.

Several previous studies investigated Streptomyces in relation to ARD. However, they resulted in controversial conclusions. Genus Streptomyces is part of the order of Actinomycetales, members of which were first mentioned as a possible cause of ARD by Otto and Winkler (1977). The authors at that time could only identify the bacteria by their morphology at the level of the phylum, which was called “Actinomycetes” in those days. In their histological analysis, “Actinomycetes” species were found in damaged roots of apple seedlings with a frequency of 47.3% in replant affected soil but not (0.3%) in steamed soil (Otto and Winkler 1977). Also, in plants from our greenhouse experiments, Actinobacteria were histologically observed more frequently in roots in untreated ARD soils than in non-ARD soils (Grunewaldt-Stöcker et al. 2019). The so called “root-pathogenic Actinomycetes” (Otto et al. 1993) were observed in ARD-affected roots of apple seedlings. Thereafter, the degree of infestation increased with increasing shoot growth and decreased with stagnating growth. This led to the assumption that root exudates, which are influenced qualitatively and quantitatively by the growing buds, triggered the germination of persistent spores (Otto et al. 1993).

In contrast to a pathogenic role, Streptomyces spp. were considered to be plant growth promoters in other studies dealing with ARD. 16S rRNA pyrosequencing revealed that the genus Streptomyces was positively (0.64) correlated with shoot growth in plants grown in fumigated ARD soil (Nicola et al. 2017). A function in disease suppression was also associated with Streptomyces (Cohen and Mazzola 2006; Cohen et al. 2005; Mazzola et al. 2007), when the effect of seed meal amendments on the putative ARD-causing pathogens Rhizoctonia solani or Pythium spp. was investigated. Seed meal amendments resulted in increased populations of Streptomyces, which were able to suppress infections by R. solani. Disease suppression was attributed to a transformation of bacterial community structure and the production of nitric oxide (Cohen and Mazzola 2006; Cohen et al. 2005), which plays a role in the induction of plant systemic resistance. Most Streptomyces isolates recovered from the apple rhizosphere were able to produce nitric oxide (Cohen et al. 2005). By adding any of several Streptomyces strains, Cohen and Mazzola (2006) could restore disease suppressiveness in previously pasteurized soil. Next to disease suppression, promotion of root infection by Streptomyces spp. was also observed in apple (Zhao et al. 2009) and Picea abies (Lehr et al. 2007). Root infections were significantly elevated in the presence of Streptomyces spp. This may be a negative side effect, because genus Streptomyces is known to promote mycorrhizal formation by promoting fungal growth and by decreasing the plant defense response (Lehr et al. 2007; Tarkka et al. 2008; Vurukonda et al. 2018). Streptomyces sp. AcH 505 was shown to downregulate the peroxidase activity and pathogenesis-related peroxidase gene (Spi2) expression (Lehr et al. 2007) of the host plant, thus promoting fungal root colonization.

Furthermore, two more traits of genus Streptomyces match the characteristics of ARD. First, like ARD, Streptomyces can persist for a very long time in soil. Due to no or minimal metabolic activity, spores can survive harsh conditions for years (Bobek et al. 2017). Second, Streptomyces is very sensitive to waterlogged conditions. Streptomyces is more abundant in drained soils (sandy loam) than in heavy soils (Gowdar et al. 2018) and, similarly, ARD is usually more severe in light soil compared with heavy soils (Mahnkopp et al. 2018; Winkelmann et al. 2019).

All of these findings indicate that Streptomyces spp. could be responsible for ARD or be part of it. However, to prove this, inoculation experiments are necessary, as was done by Tewoldemedhin et al. (2011). They isolated 92 Streptomyces strains from surface-sterilized roots from six ARD-affected sites in South Africa and inoculated 37 of them to 4-week-old apple seedlings to test pathogenicity. Moreover, 11 were coinoculated with the pathogens Pythium irregulare and C. macrodidymum. All tested streptomycetes had no effect on plant growth. At first, this seems to be a clear sign that these Streptomyces isolates were not pathogenic (directly or indirectly). However, these isolates had low identity (less than 98%) to known Streptomyces spp. and none showed close similarity to S. turgidiscabies, which had a high identity to our ASVs, with a negative correlation with plant growth. Also, in our experiments, not all Streptomyces spp. were negatively correlated with plant growth. In 2017, of 61 detected ASVs, only 6 showed negative correlations with shoot fresh mass (15 of 32 in 2016). Furthermore, inoculation trials were done in artificial soil (bark medium and sand [2:1]) (Tewoldemedhin et al. 2011), which means that potential “copathogens” were not present, unlike in ARD soil.

Streptomyces not present in T0 plants.

Only very few of the detected ASVs assigned to genus Streptomyces negatively linked to plant growth were present in T0 plants (data not shown). In 2017, only Streptomyces ASV611 was present in one of four replicates with a relative abundance of 0.38%. In 2016, ASV121 was detected in three of four replicates, one with a relative abundance of 1.39%. All other Streptomyces ASVs were not present in T0 plants. Because genus Streptomyces is widely distributed in soils (Ferrer et al. 2018; Olanrewaju and Babalola 2019; Seipke et al. 2012), plants in our experiments were most likely colonized after planting in the different soil variants. Based on molecular fingerprints of rhizosphere and bulk soil, Lucas et al. (2018) confirmed that Streptomyces is more abundant in ARD compared with grass control soil.

PGPB in plants grown in non-ARD soils.

In order to find possible PGPB that may be used to overcome ARD, we also looked at the unique ASVs in the grass variants to find ASVs positively correlated with plant growth. However, in 2016, no ASV showed any significant positive correlation (Supplementary Fig. S6). In 2017, only one ASV (NA_ASV4691 [Rhizobiaceae]) showed a positive correlation with an increase in shoot length and shoot fresh mass, with 0.32 and 0.30, respectively.

Conclusion.

Here, we showed, for the first time, the apple root endophytic community composition in plants grown on three replant-affected soils in comparison with non-replant-affected soils based on next-generation sequencing in 2 years. Although no PGPB to counteract ARD could be found, several ASVs with negative correlations with plant growth were associated with ARD. With Streptomyces spp. showing strong negative correlations and being present in all soils over the years, a potential key player for the cause of ARD may have been found. However, it remains to be clarified in future studies whether genus Streptomyces as root endophyte in ARD situations acts opportunistically or is pathogenic. Streptomyces spp. can grow saprophytically and just degrade plant material and metabolize plant exudates; however, they can also play an essential role in the ARD complex by suppressing plant defense responses and, thereby, promote infection of fungal pathogens. Further inoculation studies with Streptomyces isolates in combination with fungal pathogens as coinoculants will help to answer this question.

ACKNOWLEDGMENTS

We thank several colleagues and students at Leibniz Universität Hannover, especially E. Schneider and C. Popp, for their help; and B. Fösel from Helmholtz Center Munich for providing microbiome-associated next-generation sequencing service and proofreading the manuscript.