Microbial diversity composition of apple tree roots and resistance of apple Valsa canker with different grafting rootstock types
BMC Microbiology volume 22, Article number: 148 (2022)
The composition and diversity of root microbial community are affected by plant genotypes and soil environment, which in turn affect plant growth and development. Grafting rootstock types of the apple tree can affect phenotypes in cultivation practice, but it is not clear whether grafting rootstock types can affect the composition and diversity of root microbial community and the resistance of apple tree to apple Valsa canker.
To explore root microbial differences and the correlation, 16S rRNA and ITS genes were sequenced using Novaseq technology.
The results showed that the influence of grafting rootstock types on the composition of the root fungal community was greater than that of bacteria. And the bacterial community richness was higher in the healthy (OTUs: 1693) and dwarfing rootstock (OTUs: 1526) than in the disease (OTUs: 1181) and standard rootstock (OTUs: 1412), while the fungal community richness was the opposite. Moreover, the bacterial abundance of root zone, rhizosphere, and root endophytic microorganisms with the same grafting rootstock type exhibited a decreasing trend. Results of Nested PCR assay on soil and root tissue of Valsa mali showed that the content of V. mali in dwarfing rootstocks are lower than standard rootstocks. These results suggest that apple trees grafting with dwarfing rootstocks are more resistant to V. mali than standard rootstocks.
Under different grafting types, the effect on the composition of fungal community in apple tree root was greater than that of bacteria. The bacterial community in dwarfing rootstocks is more abundant and diverse, including more beneficial microorganisms. Therefore, dwarfing rootstock is more conducive to the resistance to apple Valsa canker from biological control.
Apple, one of the main economic crops in China, is planted by grafting the scions onto different rootstocks. The type of rootstock has a great influence on the quality of apple trees . Using MAS marker technology, it is possible to select rootstocks genotypes such as disease resistance, dwarfing, and early maturity . Studies have shown that commercial apple rootstocks have different levels of sensitivity to Apple Replant Disease (ARD), and apple rootstock G.935 might have resistance to Pythium ultimum [3, 4]. Studies on the microorganisms in the root zone of apple trees indicated that the diversity of the microbial community was affected by the rootstock and irrelevant with the scion . Root exudates might affect the diversity of soil microorganisms . Phenolics and rhizodeposits secreted by apple rootstocks can affect the composition of the microbial community in the root zone of apple trees .
Standard and dwarfing rootstocks are distinct patterns of grafting. The apple trees in standard rootstock orchards are tall, with strong branches and high drought resistance . In dwarfing rootstock orchards, flowering and fruiting periods are relatively short, with higher planting density, higher yield, better quality fruits, and so on [9, 10]. Since the 1960s, the cultivation of dwarfing rootstocks emerged in China . However, there are few studies on the utilization and cultivation techniques of dwarfing rootstock, it has not become the mainstream of apple tree grafting rootstock .
Apple tree also faces various abiotic and biotic stresses, including drought, low temperature, pests, microbial diseases, and so on [13,14,15]. In China, apple orchards are severely affected by the apple Valsa canker caused by Valsa mali [16, 17]. In some apple orchards, the disease rate is up to 30 percent, and even cause a large number of trees died [18, 19]. For the prevention and control of apple Valsa canker, chemical pesticide spraying and fruit tree peeling are generally used [17, 20]. Nevertheless, chemical control pollutes the environment, and fruit tree peeling affects the growth and development of the plant. Therefore, green biological control has gradually become a new strategy to control apple Valsa canker.
Root microorganisms form a complex network and their interactions are largely determined the beneficial traits of plants . Dominant species with specific traits can perform specific functions within the microbial community . For example, Burkholderia cepacia has a robust biological control function when acting against fungal diseases . Plant growth promoting rhizobacteria (PGPR) is selected from the soil that promotes plant growth and shows a wide range of plant diseases resistance. When plants are invaded by pathogens, PGPR will produce some antibiotics, lyases, volatiles and siderophers in time to inhibit their growth, so as to reduce the damage of pathogens to plants [24,25,26]. Meanwhile, PGPR also help plants to tolerance abiotic stresses like salt, drought, nutrient excess or deficiency . Many studies have shown that PGPR can improve the microbial community structure of plant rhizosphere soil. Wang et al. reported that microbial co-inoculants 1 (Ensifer sp. NYM3, Acinetobacter sp. P16 and Flavobacterium sp. KYM3) and microbial co-inoculants 2 significantly affected the indigenous soil bacterial community; notably, Gammaproteobacteria, Acidobacteria, Nitrospirae, and Armatimonadetes were significantly increased, while Actinobacteria and Firmicutes were significantly decreased by microbial co-inoculations . Based on quantitative PCR and DNA sequencing network analysis, Wu et al. found that Bacillus amyloliquefaciens partially inhibited the nitrification process by significantly reducing ammonium-oxidizing bacteria in soil . In terms of broad-spectrum biological control activity, compared with a single PGPR strain, the combination usage of PGPR strains has a better control effect on plant diseases [30, 31]. Soil microorganisms can help plants cope with the invasion of potato common scab, which is caused by Streptomyces spp. [32,33,34,35]. Biocontrol microorganisms can provide frontline defense against pathogen invasion . Studies have found that the infection of pathogenic bacteria usually causes changes in the soil microbial community, such as that with Oxalobacteriaceae, Burkholderiaceae, and Sphingosine bacteriaceae in the rhizosphere, which indicates that the invading pathogens directly or indirectly influence the root bacteria [37, 38]. As we all know, diseases of woody plants are caused by many factors. And the microbial communities associated with plants are complex and dynamic, which beneficial species coexist with pathogenic species . Understanding these factors and their effect on plant root microbiome will provide effective support for improving crop yield and preventing disease in the future . Do root microorganisms have a certain inhibitory effect on apple Valsa canker?
Based on the above, we hypothesized that different grafting rootstocks could affect the structure of apple root microbiome, which in turn showed different resistance to apple V. mali. Therefore, next-generation sequencing (NGS) technologies (NovaSeq 6000) were used to study the community composition of the root microorganisms of apple trees under different grafting rootstocks and disease conditions. NovaSeq 6000 relies on Illumina’s SBS chemistry and two-color reversible terminator-based method. Combined with patterned flow cell technology , in excess of 3000 Gb of data can be sequenced on an S4 flow cell. The results of this study will help to further clarify the microbial community composition and diversity of different rootstocks and provide a new perspective on the control of V. mali.
Sequencing data summary
The raw tags obtained by Illumina NovaSeq sequencing were spliced and underwent quality control practices to yield clean tags, and then chimera filtering was conducted to obtain effective tags for subsequent analysis. Operational taxonomic units (OTUs) were clustered using a 97% similarity cutoff with UPARSE (version 7.0.1090). A total of 2,131,276 bacterial sequences, 9,305 Operational Taxonomic Units (OTUs), and 2,728,919 fungal sequences, 3,971 OTUs were obtained in 36 DNA samples. Among bacterial, the number of OTUs that could be annotated to the database Silva132 was 8,674 (93.22%). It can be seen from the rarefaction curve of bacteria and fungi that the curves gradually became flat, indicating that the amount of sequencing data was reasonable (Fig. S1).
Analysis of high abundance species of the apple tree root system
The root bacterial communities of apple trees in standard rootstock and dwarfing rootstock orchards were mainly Proteobacteria, Bacteroidetes, and Actinobacteria (Fig. 1a), and the top 10 phyla accounted for 94.73% ~ 98.68% (Table S1). Meanwhile, the phyla of fungi were relatively widespread, mainly composed of Ascomycota, Mortierellomycota and Basidiomycota (Fig. 1b). And especially in root endophytic fungi and its top 10 only accounted for 16.59% ~ 36.71% (Table S1).
Due to the particularity of roots, we further analyzed niches sample communities. Through the analysis of the genus-level heat map, it can be seen that the bacterial flora of the RZ (root zone soil), RS (rhizosphere soil) and R (root endophytes) were quite different (Fig. 2a). There were certain common genera between the RZ and RS, such as, Bacillus, Sphingomonas, etc. Meanwhile, RS and R also contained some common genera, such as Devosia, Novosphingobium, Pseudoxanthomonas, etc. (Fig. 2a). However, the RZ and RS bacteria shared very few genera. In RZ, Metarhizium, Conocybe, and Microthecium were clustered in standard rootstock samples. Neonectria, Pseudogymnoascus, and Gymnoascus were clustered in dwarfing rootstock. In RS, standard rootstock samples mainly clustered Cladosporium, Minimedusa and Aureobasidium, dwarfing rootstock samples mainly clustered Vishniacozyma, Thelebolus, and Acaulium. For the root endophytic fungi, Ceratobasidium was the most clustered fungus in standard rootstock samples. Moreover, there were more clustered fungi in dwarfing rootstocks, including Alternaria, Plectosphaerella, Dactylonectria, and Paraphoma (Fig. 2b). It could be seen that the grafting rootstock types and pathogenesis had a certain influence on the composition of the fungus.
The α diversity analysis
The Chao, ACE, Shannon, and Simpson indices calculated from the fungal OTUs of all the samples indicated that α diversity of standard rootstocks was higher than that of dwarfing rootstocks (Fig. 3). Particularly, ACE and Shannon indices showed that the relative abundance and diversity of fungi in standard rootstocks were significantly higher than those in dwarfing rootstocks (Fig. 3b, 3c). Similarly, the same phenomenon was observed in bacteria, but the difference was not significant (Fig. S2). There were 2496 fungal OTUs in standard rootstock and dwarfing rootstock, and the specific OTU of standard rootstock (865) was higher than that of dwarfing rootstock (573) (Fig. 4a). However, the specific OTUs of standard rootstock (1412) was less than dwarfing rootstock (1526) in bacteria (Fig. S3a).
Based on the α diversity index, diseased and healthy had no significant effect on the richness and diversity of flora (Fig. 3, Fig. S2). The bacteria-specific OTUs of diseased samples were lower than that of healthy samples, while the OTUs of endemic fungi were more than that of the healthy group (Fig. 4b, Fig. S3b). Healthy apple trees may be enriched with certain biocontrol bacteria that increase their resistance to apple Valsa canker. Those results suggested that apple grafting rootstock types had a great impact on community diversity, followed by the whether it was infected by apple Valsa canker.
Remarkably, the richness and diversity of R were significantly lower than those in RZ and RS, as well, the richness of flora in RZ, RS, and R decreased successively (Fig. 3, Fig. S2). And there were 1720 fungal OTUs and 3992 bacterial OTUs in the three ecological niches, the specific OTUs of RZ (fungi: 470, bacteria: 1276) and RS (fungi: 413, bacteria: 858) than that in R (fungi: 197, bacteria: 538) (Fig. 4c, Fig. S3c). In addition, according to the OTUs statistics of root and soil samples, 2135 fungal OTUs and 4882 bacterial OTUs of root were found to be from soil samples, accounting for 91.55% and 90.07% of the total OTUs, respectively (Fig. 4d, Fig. S3d). This indicated that the composition of soil microorganisms has a great influence on the composition of plant root endophytic flora.
The β diversity analysis
Principal co-ordinate analysis (PCoA) based on the Bray-Curtis distance revealed that bacteria mainly accumulated according to different niches (PC1: 39.16%, PC2: 13.32%) (Fig. 5a), and the grafting rootstock types had little effect on the differences between bacterial groups (PC1:39.14%, PC2: 13.33%) (Fig. 5c). However, the fungal samples were clustered in grafting rootstock types to a certain extent (PC1:30.52%, PC2: 17.85%) (Fig. 5d), root and two soil samples were divided into two clusters (PC1: 30.52%, PC2: 17.84%) (Fig. 5b). These results showed that the influence of grafting rootstock types on the clustering of fungi was greater than that of bacteria, while the niches had a greater effect on bacterial.
Effects of planting patterns on fungal flora composition
Since the grafting rootstock types had a significant effect on the fungal community composition, t-test was used to investigate the difference in fungal composition between the two grafting rootstock types. Hypocreales, Dothideales, Cantharellales, Eurotiales, and Auriculariales at the order lever were significantly more abundant in the standard rootstock than in the dwarfing rootstock orchards (p < 0.05) (Fig. 6a). LEfSe analysis based on LDA was performed to further determine species with significant differences in the two grafting rootstock types. In the dwarfing rootstock orchards, the relative abundance of Ascomycota was much higher than that of standard rootstock orchards (Fig. 6b). On the contrary, Mucoromycota, Ascomycota, and Basidiomycota were more prevalent in the standard rootstock orchards (Fig. 6b). Fungi with high abundance in dwarfing rootstock were closer in genetic evolution, while those with high abundance in standard rootstock were farther apart in evolution (Fig. 6c).
Differences in composition of bacteria and fungi of soil and root endophytes.
The t-test of root endophytes and soil bacterial flora showed that Pseudomonas frederiksbergensis, Lysobacter sp., Acidobacteria sp., and some uncultivated bacterium in RZ and RS were significantly more abundant than those in R (p < 0.05) (Fig. 7a). Simultaneously, Bradyrhizobium elkanii, Variovorax paradoxus and Acidobacteria sp. were more likely to be present in roots as endophytes (p < 0.05) (Fig. 7a). Especially, Bradyrhizobium elkanii had been reported as PGPR, which could promote the growth of rice .
Obviously, fungi were more abundant and diverse in soil than roots (Fig. 3). For the root endophytes and soil microorganism, t-test found that Metarhizium robertsii, Aureobasidium leucospermi, Acaulium caviariforme and so on mainly existed in soil samples, and these were rarely present in R (Fig. 7b). Notably, the abundance of Metarhizium Robertsii mainly caused the differences of soil and root microbiome (Fig. S5b). Ramanpreet et.al found that Metarhizium Robertsii was not only rhizosphere competent but also could be associated with beneficial endophytic bacteria in roots to promote plant growth .
Analysis on the difference of bacterial and fungal flora composition between disease and health
LDA effect size (LEfSe) analysis was performed to further determine species with significant differences in disease and health. For bacteria, Acidovorax and Streptococcus were abundant in the health, while Rubrivivax and Luteolibacter were present in the disease (Fig. 8a). Pararhizobium giardini, Streptomyces scabrisporus, and Vibrio ponticus were more prevalent in the health, contributing 5.48%, 1.76% and 9.08% of the differential species, respectively (Fig. S4f). These may be part of the reason why it was not infected.
In the analysis of fungal composition difference, it was found that the fungi mainly existed in the disease (Fig. 8b). In the disease, Vishniacozyma, Mortierellales, Minimedusa, Cantharellales and Sordariales were more abundant, while Microascales were more likely to be present in the health (Fig. 8b). Among them, Metarhizium robertsii belonging to Ascomycota also played a major role in the differential flora (41.47%) (Fig. S4c).
Function prediction of root fungi
Using FunGuild analysis, based on the OTUs of fungi, the corresponding ecological functions of the fungi could be obtained. From the perspective of the three niches, unassigned species in the R (82.21%) more than the RS (56.56%) and RZ (49.00%), while undefined-Saprotroph (R: 7.06%) was found in RS (21.08%) and RZ (22.69%) accounted for a lower proportion (Fig. 9a). This meant that there was little knowledge surrounding the function of endophytic fungi in apple trees, and the functions of many endophytic fungi were still unknown. There were many Plant_Pathogen in the soil (RS: 5.57%; RZ: 8.08%), but the proportion of endophytic fungi in the R was small (2.67%) (Fig. 9a), so it had to do with the apple tree's own immune response, and "reject" some pathogens entering the plant from the soil.
However, owing to the strong invasiveness of fungi, it could be seen that there were still some pathogenic microorganisms in endophytes. Interestingly, in apple orchards with disease, Endophyte-Plant_Pathogen (0.94%) and Plant_Pathogen (1.86%) were lower than those without disease (Endophyte-Plant_Pathogen: 3.58%; Plant_Pathogen: 3.48%) (Fig. 9b).
It could be clearly seen from the PCA that for the soil samples, similar functional fungi were based on two types of grafting rootstocks (standard and dwarfing rootstocks) for aggregation (Fig. 9c). However, in the root endophytic fungi, the functional fungi clustered together in the standard rootstock orchard, but not in the dwarfing rootstock (Fig. S5a). Endophytic fungi could be divided into two clusters according to whether the plant is healthy (Fig. S5b). These showed that the grafting rootstock types had a larger impact on the fungal function of the RS and RZ, but had a little impact on the endophytic fungus function.
Analysis of V. mali in the root system
Apple Valsa canker has always been regarded as a branch disease, and there is no research on its control from the perspective of soil and root endophytes. Previous studies isolated 281 strains of fungi from Xinjiang wood, and found that they were all asexual, including V. mali (asexual form is Cytospora mali) and Valsa sordida (asexual form is Cytospora chrysosperma) . The Illumina NovaSeq sequencing data from this duration showed that Valsa mali (Valsa, Cytospora) was contained in the RZ (OTUs: 27) and RS (OTUs: 20), and the content of R (OTUs: 6) was low (Table 1). In addition, the number of pathogenic fungi OTUs in the diseased (33) was indeed more than that in the health (20) (Table 1). However, its existence did not determine whether the apple Valsa canker occurs. This showed that V. mali exists in the soil and was "forced" to accumulate in the roots, and it would eventually become the key factor of the diseased branches of the apple tree. In this way, the combination of branch control and soil control were used to prevention the apple Valsa canker.
In China, standard rootstocks and dwarfing rootstocks are main grafting rootstock types. In our research, we found that there were great differences in fungi OTUs under the two grafting rootstock types. The OTUs of fungi in the standard rootstock (865) was more than the dwarfing rootstock orchards (573) (Fig. 4a). In the case of Valsa canker, the OTUs of bacteria in the health was more than the diseased, including some bacterial genera with potential biocontrol effects, for example, Pararhizobium (Fig. S3b, Fig. S4f). However, the fungi in the diseased orchards were more than the healthy orchards, containing genera of fungi that had the potential to cause disease, such as Fusarium (Fig. 4b, Fig. 8b). It could be inferred from the above phenomenon that the relationship between the incidence of apple Valsa canker and root microorganisms can be defined according to whether the beneficial microbial species are enriched. Empirical evidence and theoretical predictions suggest that species-rich communities are more resistant to pathogen invasions [45, 46]. These results also confirmed our initial hypothesis.
Differences in the composition of microbes may be related to the advantages of grafting rootstock types. Some endemic microbial species were found in dwarfing rootstock. For example, Rhizobia had a good inhibition on many soilborne plant pathogenic fungi belonging to different genera like Fusarium, Rhizoctonia, Sclerotium and Macrophomina . Lysobacter can inhibit all kinds of pathogenic bacteria and fungi [48, 49]. Variovorax paradoxus and Streptomyces scabrisporus were also effective biocontrol bacteria, which could restrain the invasion of some pathogenic bacteria and have a positive effect on plant growth . Metarhizium robertsii was not only rhizosphere competent but also displayed a beneficial endophytic association with plant roots . These different microorganisms are probably a connection with the many growth advantages of dwarfing rootstock cultivation. For example, dwarfing rootstock orchards was better than standard rootstock orchards for substance absorption and transmission [52,53,54]. The dwarfed rootstock of the M9 variety has better specificity for the transmission of ABA, and also has a higher nitrogen absorption efficiency than other farming methods . The calcium absorption rate of dwarfing rootstocks was higher than that of standard rootstocks [55, 56]. However, the phosphorus absorption efficiency of different planting methods was controversial [57, 58]. Differences in microbial community composition indicated that the dwarfing rootstock might be more resistant to V. mali than the standard rootstock.
The analysis in this study found that the microflora of the three niches of the root zone, rhizosphere and root of apple trees are significant diffidence, and the abundance and richness are decrease successively. Chen et al. also showed that in mulberry trees, the richness and diversity of microbes showed a decreasing trend in root circumference, rhizosphere and roots [59, 60]. Proteobacteria, Actinobacteria, and Acidobacteria were the main bacterial groups in the root zone soil, while Proteobacteria and Actinobacteria in the rhizosphere soil, Proteobacteria and Bacteroidetes in the root. Meanwhile, Ascomycota, Mortierellomycota and Basidiomycota were the main fungal groups in the root zone, rhizosphere and root. Previous studies had shown that these microbial communities were also the dominant communities in soil and plants [61, 62]. The effect of grafting rootstock types on fungi is greater than on bacteria, while the influence of niche on bacteria and fungi is opposite. Hewavitharana et al. had shown that the type of rootstock in the greenhouse (G9.35, G.41, M.9) had a significant impact on the rhizosphere fungal community composition of apple seedlings, but had no remarkable effect on bacteria . It further clarified the understanding of the relationship between grafting rootstock types and rhizosphere microorganisms. We also found V. mali in the soil and root endophytes, which was confirmed by sequencing data and Nested PCR. V. mali in the soil may from the remnants of diseased branches, or may be a pathogenic fungus existing in the soil, thereby causing apple tree disease under appropriate conditions.
Apple Valsa canker, caused by the fungus V. mali, is one of the most important diseases in apples . Chemical pesticide spraying and fruit tree peeling were widely used to prevent and treat apple Valsa canker in China, which are environmentally unfriendly and labor-intensive [13, 14]. Through our findings, we can combine the method of biocontrol of V. mali with traditional control, improving the control effect of V. mali. Further studies will be conducted on the species with antibacterial ability in the root system of apple trees, especially the dominant species, in order to screen potential strains to control apple Valsa canker.
This study mainly analyzed the relationship between the diversity of apple tree root microbial community composition and the resistance of apple Valsa canker with different grafting rootstock types. Making full use of this diversity to select superior rootstocks and effective antagonistic microorganisms are an effective method to improve agricultural value. Theoretically, it provides ideas for studying the occurrence of plant diseases, and also provides a basis for biological and ecological control of apple Valsa canker.
Materials and methods
Sample collection and processing
All of the samples were collected at Mizhi (110°08′75.5″ E; 37°78′42.4″ N), Shaanxi, China in November 2019. Mizhi County features a mid-temperate semi-arid climate zone. The annual average temperature is 8.1℃, the annual average rainfall is 414 mm .
A total of 36 samples (3 niches × 2 grafting rootstock types × 2 healthy conditions × 3 repeats, Table 2) were from "Fuji" apple trees planted for 6 years. Three niches: soil of the root zone (RZ), rhizosphere (RS), and endophyte of the root (R); two grafting rootstock types: standard rootstocks (H) and dwarfing rootstocks (L); two healthy conditions: the trees with (Vm) and without (nVm) apple Valsa canker.
The planting distance of the standard rootstock orchards is 4 m × 5 m, and the dwarfing rootstock orchards were 4 m × 1.5 m. Five points around the diseased (Vm) and healthy (nVm) apple trees were selected as sampling points. First, remove the impurities and topsoil, and then collected the soil at a depth of 5–10 cm from the ground surface in a ziplock bag as the RZ. The soil and lateral roots 10–20 cm away from the ground surface were collected as the samples of RS and R, respectively. Shook the RS off the lateral roots of the apple tree, and then sieved the RS and RZ with a sterile 2 mm screen. Lateral roots were disinfected and used as the samples of apple tree endophyte, as previously described . The differences caused by non-experimental variables such as edge effect, branch damage and other insect pests were excluded.
DNA extraction and high-throughput sequencing
The genomic DNA of the samples were extracted using the TIANamp Soil DNA Kit and cetyltrimethylammonium bromide (CTAB) method, and PCR amplification was performed with 16S V3 region primers (341F: CCTAYGGGRBGCASCAG; 806R: GGACTACNNGGGTATCTAAT) and ITS region primers (ITS5-1737F: GGAAGTAAAAGTCGTAACAAGG; ITS2-2043R: GCTGCGTTCTTCATCGATGC). All PCR reactions were carried out with 15 µL of Phusion® High-Fidelity PCR Master Mix (New England Biolabs, USA), 6 µM of forward and reverse primers, and about 10 ng template DNA, the total volume is 30 µL. Thermal cycling consisted of initial denaturation at 98℃ for 1 min, followed by 30 cycles of denaturation at 98℃ for 10 s, annealing at 50℃ for 30 s, and elongation at 72℃ for 30 s. Finally, 72℃ for 5 min. Equal concentration mixing was performed according to PCR product concentration and 2% agarose gel electrophoresis was applied to purify the PCR products after full mixing. Screened the bands with sequence size of 400–450 bp and reclaimed the PCR products with the GeneJET gel recovery kit (Thermo Scientific, USA). Next, we shipped the recovered DNA to Novogene (Beijing, China) and used the TruSeq® DNA PCR-Free Sample Preparation Kit (Illumina, USA) for library construction. The library quality was preliminary determined by Qubit@ 2.0 Fluorometer (Thermo Scientific, USA), and Q-PCR (real-time PCR by using primers F: 5’- AATGATACGGCGACCACCGA-3’; R: 5’-CAAGCAGAAGACGGCATACGA-3’) was used for accurate and quantitative library detection. After the library was qualified, sequenced with NovaSeq6000. The sequence files of bacteria and fungi have been uploaded to NCBI, and the accession numbers are PRJNA675028 and PRJNA675150 respectively.
Data processing and analysis
To facilitate analyse the differences of microbial community composition and the resistance of apple Valsa canker with different grafting rootstocks, 12 groups were established: 1) H.RZ.Vm; 2) H.RS.Vm; 3) H.R.Vm; 4) H.RZ.nVm; 5) H.RS.nVm; 6) H.R.nVm; 7) L.RZ.Vm; 8) L.RS.Vm; 9) L.R.Vm; 10) L.RZ.nVm; 11) L.RS.nVm; 12) L.R.nVm.
The Illumina Novaseq6000 is a paired-end sequencing instrument, and the reads are paired-end. We separated each sample data from the offline data according to the barcode sequences and PCR amplification primer sequences. After trimming the barcode and primer sequences, FLASH (V1.2.7, http://ccb.jhu.edu/software/FLASH/) was applied to splice, and QIIME (V1.9.1, http://qiime.org/scripts/split libraries fastq.html) was perform to filter, and FastQC (https://www.bioinformatics.babraham.ac.uk) was utilized quality control the reads of each sample. We then performed chimera filtering to obtain effective tags that could be used for subsequent analysis [66,67,68,69]. Next, Uparse (Uparse v7.0.1001, http://drive5.com/uparse) was utilized to cluster the effective tags of all samples, and the sequences were clustered into Operational Taxonomic Units (OTUs) with 97% identity by default. The most frequently occurring sequence was screened as the representative sequence of OTUs, the OTUs sequences were annotated, and the SSUrRNA database of SILVA132 (http://www.arb-silva.de) were used to annotate the bacteria species (set threshold is 0.8 ~ 1); blasted in QIIME and Unit database (v7.2, https://unite.ut.ee) were employed to perform species annotation analysis of fungi. Using MUSCLE (Version 3.8.31, http://www.drive5.com/muscle) to perform rapidly multiple sequence alignment, and obtain the phylogenetic relationship of all OTUs representative sequences. Finally, the data of each sample was normalized, and the samples with the least amount of data was used as the standard. The subsequent α diversity analysis and β diversity analysis were based on the data after the normalization. R (Version 2.15.3) was used to draw the rarefaction curve, visualize the data and draw the PCoA, analyse the species with significant differences among groups to perform the t-test among groups and plot [70, 71]. According to the results of OTUs obtained by clustering, we analyzed the common and unique OTUs among different samples (groups), and generated a Venn diagram. Then, QIIME was applied to calculate Chao, ACE, Shannon, Simpson, and Unifrac distance matrix. The LDA score of the LEfSe analysis was set to 2, 3 or 4. LEfSe is a complex algorithm, non-parametric factorial Kruskal-Wallis (KW) sum-ranktest was used to detect species with significant differences in abundance between different groups, and then Wilcoxon rank sum test was used to judge the differences between groups. Finally, linear discriminant analysis (LDA) was used to reduce and evaluate the impact of species with significant differences (LDA Score). Finally, the FUNGuild was used to predict the functions of ITS.
Availability of data and materials
Sequence data generated and analyzed during the current study are available in the NCBI SRA, BioProject ID: PRJNA675028 and PRJNA675150, https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA675028 and https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA675150.
Soil of the root zone
Endosphere of the root
Diseased apple trees
Healthy apple trees
Root rhizosphere soil and root zone soil
Principal Component Analysis
Internal Transcribed Spacer
Molecular Marker-Assisted Selection
Apple Replant Disease
Plant growth promoting rhizobacteria
Operational Taxonomic Unit
Principal co-ordinate analysis
Linear discriminant analysis
Linear discriminant analysis Effect Size
Xu H, Iwashiro R, Li T, Harada T. Long-distance transport of Gibberellic Acid Insensitive mRNA in Nicotiana benthamiana. Bmc Plant Biology. 2013;13:165. https://doi.org/10.1186/1471-2229-13-165.
Fazio G, Mazzola M. Target traits for the development of marker assisted selection of apple rootstocks-prospects and benefits. Acta Horticulturae. 2004;663:823–8. https://doi.org/10.17660/ActaHortic.2004.663.149.
Reim S, Siewert C, Winkelmann T, Woehner T, Hanke M.V, Flachowsky H. Evaluation of Malus genetic resources for tolerance to apple replant disease (ARD). Scientia Horticulturae. 2019, 256. https://doi.org/10.1016/j.scienta.2019.05.044.
Zhu Y, Shin S, Mazzola M. Genotype responses of two apple rootstocks to infection by Pythium ultimum causing apple replant disease. Can J Plant Pathol. 2016;38:483–91. https://doi.org/10.1080/07060661.2016.1260640.
Bintarti AF, Wilson JK, Quintanilla-Tornel MA, Shade A. Biogeography and diversity of multi-trophic root zone microbiomes in michigan apple orchards: analysis of rootstock, scion, and local growing region. Phytobiomes J. 2020;4:122–32. https://doi.org/10.1094/pbiomes-01-20-0007-r.
Zhalnina K, Louie KB, Hao Z, Mansoori N, Da Rocha UN, Shi S, Cho H, Karaoz U, Loque D, Bowen BP, Firestone MK, Northen TR, Brodie EL. Dynamic root exudate chemistry and microbial substrate preferences drive patterns in rhizosphere microbial community assembly. Nat Microbiol. 2018;3:470–80. https://doi.org/10.1038/s41564-018-0129-3.
Leisso R, Rudell D, Mazzola M. Metabolic composition of apple rootstock rhizodeposits differs in a genotype-specific manner and affects growth of subsequent plantings. Soil Biol Biochem. 2017;113:201–14. https://doi.org/10.1016/j.soilbio.2017.06.011.
Ferree ME, Barden JA. Influence of strains and rootstocks on photosynthesis, respiration, and morphology of delicious apple trees. United States: N. p., 1971.
Costes E, García-Villanueva E. Clarifying the effects of dwarfing rootstock on vegetative and reproductive growth during tree development: a study on apple trees. Ann Bot. 2007;100(2):347–57. https://doi.org/10.1093/aob/mcm114.
Gao Y, Liu F, Wang K, Gong X, Liu LJ, Christopher M, et al. Genetic diversity of Malus cultivars and wild relatives in the Chinese National Repository of Apple Germplasm Resources. Tree Genet Genomes. 2015;11:106. https://doi.org/10.1007/s11295-015-0913-7.
Wang TL. The difference between apple tree stock cultivation and short stock cultivation. Hebei Fruit Tree. 2015;001:9–10. (in Chinese).
Li CY, Yang YZ, Gao JD, Wang Q, Cao HC, Du XM, et al. Research progress on dwarfing mechanism of Apple dwarfing rootstock. Chin Agric Sci Bull. 2017;33(28):86–92. (in Chinese).
Wang P, Sun X, Jia X, Ma F. Apple autophagy-related protein MdATG3s afford tolerance to multiple abiotic stresses. Plant Sci. 2017;256:53–64. https://doi.org/10.1016/j.plantsci.2016.12.003.
Wang N, Guo T, Sun X, Jia X, Wang P, Shao Y, et al. Functions of two Malus hupehensis (Pamp.) Rehd. YTPs (MhYTP1 and MhYTP2) in biotic- and abiotic-stress responses. Plant Sci. 2017;261:18-27. https://doi.org/10.1016/j.plantsci.2017.05.002.
Meng XL, Qi XH, Han ZY, Guo YB, Wang YN, Hu TL, et al. Latent infection of Valsa mali in the seeds, seedlings and twigs of crabapple and apple trees is a potential inoculum source of Valsa canker. Sci Rep. 2019;9(1):7738. https://doi.org/10.1038/s41598-019-44228-w.
Wang H, Fewer DP, Holm L, Rouhiainen L, Sivonen K. Atlas of nonribosomal peptide and polyketide biosynthetic pathways reveals common occurrence of nonmodular enzymes. P Natl Acad Sci USA. 2014;111:9259–64. https://doi.org/10.1073/pnas.1401734111.
Liu MY, Song LQ, Zhao LL. Occurrence and control of apple Valsa canker. Yantai Fruit Trees. 2020;01:43–4. (in Chinese).
Gao KX, Liu XG. Research situation of apple tree canker. J Hebei Forestry Univ. 1995;10(1):87–91. (in Chinese).
Wang JY. Prevention and control of apple canker in northern Shaanxi. Northwest Horticulture. 2006;4:19–20. (in Chinese).
Zhao ZG. Research on the occurrence regularity and prevention methods of apple Valsa canker. S China Agric. 2020. (in Chinese).
Dubuis C, Keel C, Haas D. Dialogues of root-colonizing biocontrol pseudomonads. Eur J Plant Pathol. 2007;119:311–28. https://doi.org/10.1007/s10658-007-9157-1.
Yang T, Wei Z, Friman VP, Xu Y, Shen Q, Kowalchuk GA, et al. Resource availability modulates biodiversity-invasion relationships by altering competitive interactions. Environ Microbiol. 2017;19:2984–91. https://doi.org/10.1111/1462-2920.13708.
Bach E, dos Santos Seger GD, Fernandes GdC, Lisboa BB, Pereira Passaglia LM. Evaluation of biological control and rhizosphere competence of plant growth promoting bacteria. Appl Soil Ecol. 2016;99:141–9. https://doi.org/10.1016/j.apsoil.2015.11.002.
Ali S, Hameed S, Imran A, Iqbal M, Lazarovits G. Genetic, physiological and biochemical characterization of Bacillus sp. strain RMB7 exhibiting plant growth promoting and broad spectrum antifungal activities. Microbial Cell Fact. 2014;13:144. https://doi.org/10.1186/s12934-014-0144-x.
Tjamos SE, Flemetakis E, Paplomatas EJ, Katinakis P. Induction of resistance to Verticillium dahlia in Arabidopsis thaliana by the biocontrol agent K-165 and pathogenesis-related proteins gene expression. Mol Plant Microbe Interact. 2005;18:555–61. https://doi.org/10.1094/MPMI-18-0555.
Yu XM, Ai CX, Xin L, Zhou G. The siderophore-producing bacterium, Bacillus subtilis CAS15, has a biocontrol effect on Fusarium wilt and promotes the growth of pepper. Eur J Soil Biol. 2011;47:138–45.
Cohen AC, Bottini R, Pontin M, Berli FJ, Moreno D, Boccanlandro H, et al. Azospirillum brasilense ameliorates the response of Arabidopsis thaliana to drought mainly via enhancement of ABA levels. Physiol Plant. 2015;153:79–90. https://doi.org/10.1111/ppl.12221.
Wang J, Li Q, Xu S, Zhao W, Lei Y, Song C, et al. Traits-based integration of multi-species inoculants facilitates shifts of indigenous soil bacterial community. Front Microbiol. 2018;9:1692. https://doi.org/10.3389/fmicb.2018.01692.
Wu S, Zhuang G, Bai Z, Cen Y, Xu S, Sun H, et al. Mitigation of nitrous oxide emissions from acidic soils by Bacillus amyloliquefaciens, a plant growth-promoting bacterium. Glob Chang Biol. 2018;24:2352–65. https://doi.org/10.1111/gcb.14025.
Liu K, McInroy JA, Hu CH, Kloepper JW. Mixtures of plant growth-promoting rhizobacteria enhance biological control of multiple plant diseases and plant growth promotion in the presence of pathogens. Plant Dis. 2018;102:67–72. https://doi.org/10.1094/PDIS-04-17-0478-RE.
Fierer N. Embracing the unknown: disentangling the complexities of the soil microbiome. Nat Rev Microbiol. 2017;15:579–90. https://doi.org/10.1038/nrmicro.2017.87.
Schlatter D, Kinkel L, Thomashow L, Weller D, Paulitz T. Disease suppressive soils: new insights from the soil microbiome. Phytopathology. 2017;107:1284–97. https://doi.org/10.1094/phyto-03-17-0111-rvw.
Shi W, Li M, Wei G, Tian R, Li C, Wang B. The occurrence of potato common scab correlates with the community composition and function of the geocaulosphere soil microbiome. Microbiome. 2019;7:1. https://doi.org/10.1186/s40168-019-0629-2.
Macnish GC, Neate SM. Rhizoctonia bare patch of cereals - An Australian perspective. Plant Dis. 1996;80:965–71. https://doi.org/10.1094/pd-80-0965.
Berendsen RL, Pieterse CMJ, Bakker PAHM. The rhizosphere microbiome and plant health. Trends Plant Sci. 2012;17:478–86. https://doi.org/10.1016/j.tplants.2012.04.001.
Chapelle E, Mendes R, Bakker PAHM, Raaijmakers JM. Fungal invasion of the rhizosphere microbiome. Isme J. 2016;10:265–8. https://doi.org/10.1038/ismej.2015.82.
Zhang Y, Xu J, Riera N, Jin T, Li J, Wang N. Huanglongbing impairs the rhizosphere-to-rhizoplane enrichment process of the citrus root-associated microbiome. Microbiome. 2017;5:97. https://doi.org/10.1186/s40168-017-0304-4.
Berendsen RL, Vismans G, Yu K, Song Y, Jonge R, Burgman WP. Disease-induced assemblage of a plant-beneficial bacterial consortium. Isme J. 2018;12:1496–507. https://doi.org/10.1038/s41396-018-0093-1.
Poudel R, Jumpponen A, Kennelly MM, Rivard CL, Gomez-Montano L, Garrett KA. Rootstocks shape the rhizobiome: rhizosphere and endosphere bacterial communities in the grafted tomato system. Appl Environ Microbiol. 2019;85:9. https://doi.org/10.1128/aem.01765-18.
Hao XZ, Gao ZL, Shi CC, Zhang JP, Cao QX. Analysis on dynamic changes of climatic factors in Yulin. J Inn Mongolia Forestry Sci Technol. 2021;47(03):51–4. (in Chinese).
Kumar KR, Cowley MJ, Davis RL. Next-generation sequencing and emerging technologies. Thieme Med Publishers. 2019;45:661–73. https://doi.org/10.1055/s-0039-1688446.
Ahmad I, Jiménez-Gasco MM, S. Luthe DE, Barbercheck M. Systemic colonization by Metarhizium robertsii enhances cover crop growth. J Fungi. 2020;6(2). https://doi.org/10.3390/jof6020064.
Sasan KR, Bidochka JM. The insect-pathogenic fungus Metarhizium robertsii (Clavicipitaceae) is also an endophyte that stimulates plant root development. Am J Bot. 2012;99(1). https://doi.org/10.3732/ajb.1100136.
Guo KF. Identification of forest rot pathogens in Xinjiang and study on genetic differentiation and rapid detection of primary infection sources of Valsa mali. Shihezi University; 2016. (in Chinese).
Van Elsas JD, Chiurazzi M, Mallon CA, Elhottova D, Kristufek V, Salles JF. Microbial diversity determines the invasion of soil by a bacterial pathogen. PNAS. 2012;109:1159–64. https://doi.org/10.1073/pnas.1109326109.
Wei Z, Yang T, Friman VP, Xu Y, Shen Q, Jousset A. Trophic network architecture of root-associated bacterial communities determines pathogen invasion and plant health. Nat Commun. 2015;6:8413. https://doi.org/10.1038/ncomms9413.
Bonaterra A, Mari M, Casalini L, Montesinos E. Biological control of Monilinia laxa and Rhizopus stolonifer in postharvest of stone fruit by Pantoea agglomerans EPS125 and putative mechanisms of antagonism. Int J Food Microbiol. 2003;84:93–104. https://doi.org/10.1016/s0168-1605(02)00403-8.
Ji GH, Wei LF, He YQ, Wu YP, Bai XH. Biological control of rice bacterial blight by Lysobacter antibioticus strain 13–1. Biol Control. 2008;45:288–96. https://doi.org/10.1016/j.biocontrol.2008.01.004.
Choi H, Kim HJ, Lee JH, Kim JS, Park SK, Kim IS, et al. Insight into genes involved in the production of extracellular chitinase in a biocontrol bacterium Lysobacter enzymogenes C-3. Plant Pathol J. 2012;28:439–45. https://doi.org/10.5423/ppj.Nt.07.2012.0115.
Nguyen Thanh H, Tran Quang M, Pham Xuan H, Nguyen Thi Thanh T, Furuya N, Hoang HL. Biological control of potato tuber soft rot using N-acyl-L-homoserine lactone-degrading endophytic bacteria. Curr Sci. 2018;115:1921-1927. https://doi.org/10.18520/cs/v115/i10/1921-1927.
Fravel D, Olivain C, Alabouvette C. Fusarium oxysporum and its biocontrol. New Phytologist. 2003;157:493–502. https://doi.org/10.1046/j.1469-8137.2003.00700x.
Fu ZQ, Dong X. Systemic acquired resistance: turning local infection into global defense. Annu Rev Plant Biol. 2013;64:839–63. https://doi.org/10.1146/annurev-arplant-042811-105606.
Kamboj JS, Browning G, Blake PS, Quinlan JD, Baker DA. GC-MS-SIM analysis of abscisic acid and indole-3-acetic acid in shoot bark of apple rootstocks. Plant Growth Regul. 1999;28:21–7. https://doi.org/10.1023/a:1006299414481.
Else MA, Taylor JM, Young S, Atkinson CJ. The effect of the graft union on hormonal and ionic signalling between rootstocks and scions of grafted apple (Malus pumila L. Mill.). Environ Exp Bot. 2018;156:325–36. https://doi.org/10.1016/j.envexpbot.2018.07.013.
Polat D, Yildirim F, Yildirim AN. Identification of minerals in leaves of 14 apple rootstocks in stool bed condition. Erwerbs-Obstbau. 2020;62:77–81. https://doi.org/10.1007/s10341-019-00461-1.
Amiri ME, Fallahi E, Safi-Songhorabad M. Influence of rootstock on mineral uptake and scion growth of “golden delicious” and “royal gala” apples. J Plant Nutr. 2014;37:16–29. https://doi.org/10.1080/01904167.2013.792838.
Ershadi A, Talaie A. The effect of clonal rootstocks on leaf mineral composition of several apple cultivars. Neutron News. 2001;12:317–20. https://doi.org/10.17660/ActaHortic.2001.564.37.
Kucukyumuk Z, Erdal I. Rootstock and cultivar effect on mineral nutrition, seasonal nutrient variation and correlations among leaf, flower and fruit nutrient concentrations in apple trees. Bulg J Agric Sci. 2011;17:633–41.
Chen LJ, Wu XQ, Xu Y, Wang BL, Liu S, Niu JF, et al. Microbial diversity and community structure changes in the rhizosphere soils of Atractylodes lancea from different planting years. Plant Signal Behav. 2021;16(2):1854507. https://doi.org/10.1080/15592324.2020.1854507.
Zhao Y, Fu W, Hu C, Chen G, Xiao Z, Chen Y, et al. Variation of rhizosphere microbial community in continuous mono-maize seed production. Sci Rep. 2021;11(1):1544. https://doi.org/10.1038/s41598-021-81228-1.
Yu W, Jun T, Xiao L, Wu RY, Mo DXC. Variation in fungal community structures in rhizosphere soil of Coptis chinensis with cropping mode under natural forest and artificial shed. Chin J Chin Mater Med. 2020;45(21):5160–8. (in Chinese).
Lin M, Xiong H, Xiang X, Zhou Z, Liang L, Mei Z. The effect of plant geographical location and developmental stage on root-associated microbiomes of Gymnadenia conopsea. Front Microbiol. 2020;11:1257. https://doi.org/10.3389/fmicb.2020.01257.
Hewavitharana SS, Mazzola M. Influence of rootstock genotype on efficacy of anaerobic soil disinfestation for control of apple nursery replant disease. Eur J Plant Pathol. 2020;157:39–57. https://doi.org/10.1007/s10658-020-01977-z.
Wang S, Hu T, Wang Y, Luo Y, Michailides TJ, Cao K. New understanding on infection processes of Valsa canker of apple in China. Eur J Plant Pathol. 2016;146:531–40. https://doi.org/10.1007/s10658-016-0937-3.
McPherson MR, Wang P, Marsh EL, Mitchell RB, Schachtman DP. Isolation and analysis of microbial communities in soil, rhizosphere, and roots in perennial grass experiments. Jove-J Vis Exp. 2018;137. https://doi.org/10.3791/57932.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6. https://doi.org/10.1038/nmeth.f.303.
Magoc T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27:2957–63. https://doi.org/10.1093/bioinformatics/btr507.
Bokulich NA, Subramanian S, Faith JJ, Gevers D, Gordon JI, Knight R. Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nat Methods. 2013;10:57-U11. https://doi.org/10.1038/nmeth.2276.
Rognes T, Flouri T, Nichols B, Quince C, Mahe F. VSEARCH: a versatile open source tool for metagenomics. Peer J. 2016;4:2584. https://doi.org/10.7717/peerj.2584.
Xu WF, Wang F, Wang RL, Sui Y, Zhou ZY, Xie J, et al. Seasonal characterization of the endophytic fungal microbiome of mulberry (Morus spp.) cultivars resistant and susceptible to sclerotiniosis. Microorganisms. 2021;9(10):2052. https://doi.org/10.3390/microorganisms9102052.
Fang H, Wang H, Cai L, Yu Y. Prevalence of antibiotic resistance genes and bacterial pathogens in long-term manured greenhouse soils as revealed by metagenomic survey. Environ Sci Technol. 2015;49:1095–104. https://doi.org/10.1021/es504157v.
We are grateful to the fund projects: Natural Science Basic Research Program of Shaanxi (2020zdzx03-03-01); National Natural Science Foundation of China (U19032061007919); National Key R&D Program (2017YFD0200602-2); National Natural Science Foundation of China (32072477); Shaanxi Provincial Natural Science Basic Research Program Key Project (2017JZ006).
Ethics approval and consent to participate
Our sampling has been approved by Plant Protection and Inspection Station, and we comply with the policy statement of the World Conservation Union on endangered species research and the Convention on trade in endangered species of Wild Fauna and flora.
Consent for publication
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. All authors have read and consent to the submission and publication of this manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Wang, J., Wang, R., Kang, F. et al. Microbial diversity composition of apple tree roots and resistance of apple Valsa canker with different grafting rootstock types. BMC Microbiol 22, 148 (2022). https://doi.org/10.1186/s12866-022-02517-x