Comparison and interpretation of characteristics of Rhizosphere microbiomes of three blueberry varieties

Background Studies on the rhizosphere microbiome of various plants proved that rhizosphere microbiota carries out various vital functions and can regulate the growth and improve the yield of plants. However, the rhizosphere microbiome of commercial blueberry was only reported by a few studies and remains elusive. Comparison and interpretation of the characteristics of the rhizosphere microbiome of blueberry are critical important to maintain its health. Results In this study, a total of 20 rhizosphere soil samples, including 15 rhizosphere soil samples from three different blueberry varieties and five bulk soil samples, were sequenced with a high-throughput sequencing strategy. Based on these sequencing datasets, we profiled the taxonomical, functional, and phenotypic compositions of rhizosphere microbial communities for three different blueberry varieties and compared our results with a previous study focused on the rhizosphere microbiome of blueberry varieties. Our results demonstrated significant differences in alpha diversity and beta diversity of rhizosphere microbial communities of different blueberry varieties and bulk soil. The distribution patterns of taxonomical, functional, and phenotypic compositions of rhizosphere microbiome differ across the blueberry varieties. The rhizosphere microbial communities of three different blueberry varieties could be distinctly separated, and 28 discriminative biomarkers were selected to distinguish these three blueberry varieties. Core rhizosphere microbiota for blueberry was identified, and it contained 201 OTUs, which were mainly affiliated with Proteobacteria, Actinobacteria, and Acidobacteria. Moreover, the interactions between OTUs of blueberry rhizosphere microbial communities were explored by a co-occurrence network of OTUs from an ecological perspective. Conclusions This pilot study explored the characteristics of blueberry’s rhizosphere microbial community, such as the beneficial microorganisms and core microbiome, and provided an integrative perspective on blueberry’s rhizosphere microbiome, which beneficial to blueberry health and production. Supplementary Information The online version contains supplementary material available at 10.1186/s12866-021-02092-7.


Background
The rhizosphere of plants harbors diverse microorganisms in soil, which evolve alongside plants and environments and form an integral part of the life cycle of plants. Rhizosphere microbiota carry out a variety of vital functions and play a critical role in biogeochemical cycles involving soil formation and carbon cycling [1]. For example, many rhizosphere microorganisms provide nutrients to plants from soil [2] and prevent plants from being infested by pathogens [3]. The complex and dynamic interactions between plants and microbiota, especially between microorganisms, are related to the growth of plants [4,5]. Hence, understanding the taxonomical and functional compositions of rhizosphere microbial community is bene cial to the growth and yield of plants. In recent decades, many studies have been conducted to characterize rhizosphere microbiome in speci c crop plant species, including rice [6], soybean [7], corn [8], barley [9], and wheat [10], and vegetable and fruit crops, including sugarcane [11], cucumber [12], grapevine [13] and citrus [14,15]. A majority of these studies were performed through high-throughput sequencing of the microbial 16S rRNA to fully explore and characterize the role of microbiota in the rhizosphere microbial community. Several consistent trends and speci c traits were demonstrated on the basis of many studies on rhizosphere microbiome of plants. For example, the number of bacteria a liated with Alphaproteobacteria in rhizosphere microbial communities of various plants increases [9,16,17]. However, the current studies on the rhizosphere microbiome are primarily on model plants, and relatively few studies related to blueberry have been carried out to explore the taxonomical and functional compositions of the blueberry rhizosphere microbial community [18], especially for the rhizosphere microbiome of different blueberry varieties [19].
Blueberries are perennial owering plants known for their blue or purple berries. In taxonomy, the species of blueberry are classi ed into the Vaccinium genus. The commercial blueberries are all native to North America and different kinds of blueberries were later introduced to Asia and Europe [20]. In recent years, numerous studies have investigated the effects of blueberry for consumer's health based on their composition in avonoids, polyphenols, anthocyanins, pro-anthocyanidins, phenolic acids and stilbenes, and demonstrated that the anti-oxidant and anti-in ammatory activities of blueberry [21,22]. Moreover, previous studies have explored the dynamic changes of human or mice gut microbes with the consumption of blueberry or its extracts [23,24]. Six-week regular consumption of wild blueberry drink can positively modulate the composition of human gut microbiota and increase the content of Bi dobacteria [23], which have been shown to exert positive bene ts to humans health [25]. Additionally, growing evidence suggested that avonoids of blueberry have the potential to restrict the development and severity of certain cancers and vascular diseases [26]. Given these bene ts, more blueberries are needed and consumed. However, the diseases of blueberry, such as stem and leaf diseases, including phomopsis leaf spot and fruit rot and septoria leaf spot, reduce the yield of blueberry [27][28][29]. Previous studies have suggested that the rhizosphere microbiome can in uence plant susceptibility to diseases and tness response to the changes of environmental factors [30,31], and several diseases of plants are related to rhizosphere microbiota in soil and can be controlled by related microbes [32,33]. Therefore, understanding the blueberry rhizosphere microbiome and comparing the differences in rhizosphere microbial communities of different blueberry varieties, including the universal microbiota (shared microbiota) between different kinds of blueberry varieties and speci c microbiota of each blueberry, are favorable to the cultivation and agricultural management of blueberries. However, only a few study has explored and illustrated the composition of rhizosphere microbiome of blueberry to date [34] and the compositions of blueberry rhizosphere microbial community remains elusive.
In this present study, we collected 15 rhizosphere soil samples of three kinds of blueberry varieties and ve adjacent soil samples (bulk soil) from a blueberry plantation in Hefei City, China, on 13 April 2018. To pro le the structure of rhizosphere microbial community of blueberries, we performed 16S rRNA amplicon sequencing for these samples and analyzed the sequencing data. In our work, we focused on the following scienti c questions: (i) How does the microbial diversity differ between rhizosphere microbial communities of different blueberry varieties? (ii) What are differences in taxonomical, functional, and phenotypic compositions between rhizosphere microbial communities of different blueberry varieties?
(iii) What are the core microbiota of rhizosphere microbial communities in blueberry? (iv) How are the cooccurrence relationships between microbiota in different blueberry varieties? Notably, our study aims to compare and interpret the characterization of the blueberry rhizosphere microbial community and explore the patterns of the blueberry rhizosphere microbial community, which could provide an integrative view on the blueberry rhizosphere microbiome and provide insights on keeping blueberry health to improve the production of blueberry.

Methods
Collection of rhizosphere soil samples Three blueberry varieties, namely Rabbiteye Blueberry, Northern Highbush Blueberry and Southern Highbush Blueberry, were selected from a blueberry plantation in Hefei City, Anhui province, China, to investigate the structure of blueberry rhizosphere microbial community and explore the differences among three different blueberry varieties. The selected plants of three blueberry varieties have been planted for 6-7 years. The rhizosphere soil samples of these three blueberry varieties were collected according to the sampling procedure [6,14]  assessed on agarose gels, respectively. The V3-V4 hypervariable regions of the 16S rRNA gene of microbes for each rhizosphere soil sample were ampli ed and sequenced to pro le the structure of the blueberry rhizosphere microbial community. Speci cally, approximately 50 ng DNA was used as PCR template, and the forward primer 347F 5¢-CCTACGGRRBGCASCAGKVRVGAAT-3¢ and reverse primer 802R 5¢-GGACTACNVGGGTWTCTAATCC-3¢ were used to amplify the V3-V4 amplicons [35]. Indexed adapters were added to the ends of 16S rDNA amplicons and the sequencing library was constructed. The sequencing library was veri ed, quanti ed, and sequenced on an Illumina MiSeq platform (San Diego, CA, USA) using the paired-end sequencing strategy.
16S rRNA amplicon data processing and taxonomical pro les The paired-end reads of 16S rDNA amplicons of each sample were spliced by using the Fast Length Adjustment of Short reads (FLASH, v1.2.11) software [36] with default settings. The spliced reads containing ambiguous base calls (N) were removed, and the lengths of spliced reads ranging from 220 bp to 550 bp were chosen by using "trim.seqs" command in the mothur platform [37] (version 1.25.0). The putative chimeras were identi ed against the SILVA database [38] (release 123) and removed in the mothur platform. The high-quality sequences were used for taxonomical analysis against the Greengenes database [39] (version 13_5) in QIIME (Quantitative Insights Into Microbial Ecology, Boulder, CO, USA, v1.9.1) [40]. The operational taxonomic units (OTUs) were clustered at the 97% nucleotide identity threshold by using the "pick_closed_reference_otus.py" script, and the singletons of OTUs were removed. The nal OTU table was rare ed to 18,652 reads per sample prior to downstream analysis to eliminate the effect of sequencing depth.

Functional and phenotypic compositions of blueberry rhizosphere microbiome
To compare the differences in functional and phenotypic compositions of rhizosphere microbial communities of different blueberry varieties, two popular tools in current microbiome analysis, namely Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt, version: 1.0.0-dev) [41] and Bugbase [42], were selected to pro le the characters of blueberry rhizosphere microbial communities. Speci cally, the functional compositions of blueberry rhizosphere microbial communities were predicted, and the relative abundance of each functional trait that collapsed to levels two and three of the KEGG database (version 66.1, May 1, 2013) was summarized based on the OTU composition. Similarly, the phenotypic compositions of rhizosphere microbial communities, including the content of anaerobic, mobile elements and stress tolerant, were pro led.
Microbial diversity assessment of blueberry rhizosphere microbial communities The number of OTUs, Shannon index and Simpson index of rhizosphere soil samples were selected to evaluate the alpha diversities of rhizosphere microbial communities among three blueberry varieties and bulk soil samples. The alpha diversity was compared using the Kruskal-Wallis test among three blueberry varieties and bulk soil. Bray-Curtis distance and unweighted UniFrac metrics (refers to one of the UniFrac metrics, and it only considers the presence or absence of observed microorganisms) [43] were used to compare the differences of beta diversity among three blueberry varieties and bulk soil samples.
The clustering result of the rhizosphere microbial community was arrayed by principle coordination analysis (PCoA) and visualized using Emperor [44]. Moreover, linear discriminate analysis (LDA) was performed to utilize a linear combination of features to maximize the separation of rhizosphere microbial communities of three blueberry varieties and bulk soil based on taxonomical composition at the phylum, order, and genus levels. On the basis of the Bray-Curtis distance metric of taxonomical composition of the genus level, permutational multivariate analysis of variance (PERMANOVA) [45] was used to evaluate whether the rhizosphere microbial communities are signi cantly different across three blueberry varieties and bulk soil. To determine if other taxa were stable among three blueberry varieties and bulk soil, we identi ed the core microbiome among rhizosphere samples across groups and visualized the results by venn plot and heatmap in R.

Biomarker analysis
Linear discriminate analysis effect size (LEfSe, version 1.0) [46] was applied to select the differentially taxonomical features among rhizosphere microbial communities of three blueberry varieties and bulk soil samples. The p-value for the factorial Kruskal-Wallis test was set at 0.05 to select statistical signi cant taxonomical biomarkers. Biomarker with the logarithmic LDA score higher than 3.5 was de ned as a discriminative biomarker and visualized.
Co-occurrence network in blueberry rhizosphere microbial community The correlations among OTUs of the rhizosphere microbial community of blueberry were calculated by using the SparCC algorithm (https://github.com/hallamlab/utilities/wiki/SparCC), which limits the number of spurious correlation identi ed [47,48]. The threshold of absolute correlations among OTUs was set at 0.8 and the signi cant correlations with p value < 0.05 were visualized in Cytoscape [49] (version 3.7.1). The characteristics of the topological structure of the co-occurrence network were analyzed in igraph package [50] (version 0.7.1) in R.

Results
Differential microbial diversity in blueberry rhizosphere microbial community To pro le the taxonomical compositions of rhizosphere microbial communities of blueberry and compare the taxonomical differences for these three blueberry varieties and bulk soil samples, we sequenced the  Table 1).
The alpha diversities of rhizosphere microbial communities were compared between three blueberry varieties and bulk soil samples using the number of OTUs, Shannon index, and Simpson index (Figs. 1a-1c). We observed that the number of OTUs of microbial communities in bulk soil samples was signi cantly higher than that in three blueberry varieties, and the number of OTUs of rhizosphere microbial communities among three blueberry varieties was also signi cantly different (Kruskal-Wallis test, p < 0.05; Fig. 1a). As for species richness of rhizosphere microbial communities, we found that the Shannon and Simpson indexes of rhizosphere microbial communities of bulk soil samples were signi cantly higher than those of blueberry varieties, except Southern Highbush Blueberry (Figs. 1b-1c).
The similarities of rhizosphere microbial communities were also assessed among three blueberry varieties and bulk soil samples based on Bray-Curtis (Fig. 1d) and unweighted UniFrac distance metrics (Fig. 1e). The results of PCoA based on Bray-Curtis (Fig. 1d) and unweighted UniFrac distance metrics ( Fig. 1e) revealed signi cant differences in taxonomical compositions between microbial communities of blueberry rhizosphere soil and bulk soil (p < 0.001, F = 6.815, One-way PERMANOVA, N = 9,999, Bray-Curtis dissimilarity index). The taxonomical compositions of rhizosphere microbial communities of three blueberry varieties also signi cantly differed (p < 0.001, F = 7.472, One-way PERMANOVA, N = 9,999, Bray-Curtis dissimilarity index).
Differential taxonomical composition in blueberry rhizosphere microbial community To gain insights into the taxonomical compositions of blueberry rhizosphere microbial communities, we strati ed the taxonomical structure of rhizosphere microbial communities at the phylum, order, and genus levels (Fig. 2). We compared the differences in taxonomical compositions not only between rhizosphere microbial communities of blueberry and bulk soil but also among different blueberry varieties.
At the phylum level, we found that Proteobacteria, Actinobacteria, Acidobacteria, Firmicutes, Planctomycetes, and Verrucomicrobia constituted the six most enriched bacterial phyla among rhizosphere microbial community of three blueberry varieties and bulk soil (Fig. 2a). The result of predominant phylum is almost consistent with a previous study, which also reported that Proteobacteria, Actinobacteria, Acidobacteria, Bacteroidetes, Planctomycetes, Chloro exi, and Verrucomicrobia were enriched in the rhizosphere microbiome of blueberry [19]. communities of three blueberry varieties and bulk soil based on the relative abundances of predominant phyla. We observed that rhizosphere microbial communities of three blueberry varieties and bulk soil could be distinctly differentiated by integrating a linear combination of phyla (Fig. 2b). Among the linear combination of phyla, we found that Planctomycetes, Gemmatinonadetes, Chloro exi, and Verrucomicrobia were important for differentiating rhizosphere microbial communities of three blueberry varieties and bulk soil (Fig. 2b).

Differential functional and phenotypic compositions in blueberry rhizosphere microbial community
The functional and phenotypic compositions in blueberry's rhizosphere microbial community were pro led based on their taxonomical compositions (Fig. 3). As to the functional traits that collapsed to level 2 of the KEGG database, we found that the enrichment of enzyme families and environmental adaptation in rhizosphere microbial communities, and the proportion of biosynthesis of other secondary metabolites was higher in Northern Highbush Blueberry (Supplementary Figure 3). The relative abundances of functional traits related to transporters, general function, ABC transporters, DNA repair and recombination proteins, two-component system, and urine metabolism were higher in the rhizosphere microbial community of blueberry varieties and bulk soil (Fig 3a). Moreover, we found that the functional compositions of the rhizosphere microbial communities of Rabbiteye Blueberry signi cantly differed from those of bulk soil (p < 0.05, F = 3.545 One-way PERMANOVA, N = 9,999, Bray-Curtis dissimilarity index) and Southern Highbush Blueberry (p < 0.05, F = 3.3, One-way PERMANOVA, N = 9,999, Bray-Curtis dissimilarity index). The rhizosphere microbial communities of three blueberry varieties and bulk soil could be distinctly distinguished by integrating a linear combination of functional components (Fig. 3b).

Core blueberry rhizosphere microbiome
We extended our analysis to determine which OTUs are stable across in rhizosphere microbial communities of different blueberry varieties and bulk soil. We identi ed 728, 634, 777, and 712 OTUs as the core OTUs in rhizosphere microbial communities of Rabbiteye Blueberry, Northern Highbush Blueberry, Southern Highbush Blueberry and bulk soil (Fig. 4a), respectively. Eventually, 201 OTUs of 1,420 OTUs (14.2%) were identi ed as the core OTUs in rhizosphere microbial communities of blueberry varieties and bulk soil (Fig. 4a, Supplementary Table 2). Many OTU cases mainly a liated with Proteobacteria (78 OTUs), Actinobacteria (41 OTUs), Acidobacteria (34 OTUs), Firmicutes (16 OTUs), Chloro exi (9 OTUs), and Planctomycetes (8 OTUs, Fig. 4b). The distribution of each core OTU in rhizosphere microbial communities of blueberry varieties was different (Fig. 4b), indicating that the relative abundance of core OTUs varied most among different blueberry varieties.

Identi cation of microbial biomarkers for classifying different blueberry varieties
To explore the taxonomical signatures among rhizosphere microbial communities of three blueberry varieties and bulk soil, we conducted LEfSe analysis to identify biomarkers for each blueberry variety based on the taxonomical compositions of rhizosphere microbial communities. Finally, we obtained 28 discriminative biomarkers with logarithmic LDA score > 3.5 (Fig. 5). At the phylum level, we found that Actinobacteria and Planctomycetes were identi ed as the biomarkers for Rabbiteye Blueberry and Southern Highbush Blueberry, respectively, whereas Verrucomincrobia and Chloro exi were detected as the biomarkers for bulk soil (Fig. 5a). At the order level, we observed that Clostridiales, Rhodospirillales, Rhizobiales, Gaiellales, Actinomycetals, Xanthomonadales, and Burkholderiales were identi ed as the biomarkers for three blueberry varieties (Fig. 5).
Patterns of co-occurrence network in blueberry rhizosphere microbial community To gain more insights into the interactions among the microbial members of rhizosphere microbial communities of blueberry varieties, we extended our analysis to explore the patterns of OTUs cooccurrence network from an ecological perspective. The SparCC algorithm was applied to calculate the correlations between OTUs and the signi cant strong correlations (the value of absolute correlations > 0.8 and the p value < 0.05) were chosen to construct the co-occurrence network. The co-occurrence network comprised 198 nodes and 484 edges (Fig. 6). The density and average degree of the co-occurrence network were 0.025 and 4.89, respectively. The clustering coe cient of the co-occurrence network was 0.35 and the co-occurrence network could be clustered into seven clusters. The strong interactions existed between OTUs in the co-occurrence network. The members of co-occurrence network were mainly a liated with Proteobacteria, Actinobacteria, Acidobacteria, Verrucomicrobia, and Firmicutes (Fig. 6). Among the 198 nodes, 74 nodes (37.4%) belonged to core OTUs and these OTUs were mainly a liated with Proteobacteria, Actinobacteria, and Acidobacteria (Fig. 6). The OTUs with the highest average proportions of the co-occurrence network were members of core OTUs of rhizosphere microbial communities of blueberry varieties, which were a liated with Xanthomonadaceae, Koribacteraceae, Gaiekkaceae, and Sinobacteraceae (Fig. 6).

Discussion
To obtain a better understanding of rhizosphere microbiome of blueberry, this pilot study was conducted and mainly focused on the taxonomical, functional, and phenotypic compositions of rhizosphere microbial communities in blueberry. By investigating the compositions of blueberry's rhizosphere microbial communities and comparing the differences in rhizosphere microbial communities among three blueberry varieties, the characterization of blueberry's rhizosphere microbial community and the interactions between rhizosphere microbiota should be understood to provide new opportunities to increase the yield of blueberry [3,15].
Previous studies have reported that plants can shape and recruit protective microorganisms from the soil microbial community to form the rhizosphere microbial community [3,51], leading to a difference between plants' rhizosphere microbial community and bulk soil microbial community. In our study, the alpha diversity and beta diversity of rhizosphere microbial communities of blueberry varieties and bulk soil were signi cantly differed. On the basis of the taxonomical composition, we observed that the microbial diversity of blueberry's rhizosphere microbial communities decreased compared with bulk soil samples. The decrease in the diversity of rhizosphere microbial community was also found in a previous study of blueberry focused on the taxonomical composition of bulk soil and rhizosphere microbial communities [18]. Furthermore, the distribution patterns of three blueberry rhizosphere microbial communities and bulk soil were different at the phylum, order, and genus levels. Phyla Actinobacteria, Firmicutes, and Planctomycetes were dominant in the rhizosphere microbial community of three blueberry varieties. In terms of Firmicutes, previous studies have reported that the members of Firmicutes are identi ed as groups of bacteria that can confer suppressiveness and important in disease suppressiveness in rhizosphere microbiota of plants [3,52]. Similarly, Actinomycetales was enriched in the blueberry rhizosphere microbial community, which was detected as the dominated group in rhizosphere soil alongside with crop growth [53,54]. The differences in rhizosphere microbial communities between three blueberry varieties and bulk soil samples revealed that a series of microbiota were recruited from the soil microbial community to form the rhizosphere microbial community of blueberry. Additionally, there were signi cant differences among the rhizosphere microbial communities of three blueberry varieties by comparing the discrepancy of their rhizosphere microbial communities. These results suggested that blueberry can recruit different microbiota to determine the composition of the rhizosphere microbiome and con rmed that different genotype blueberry varieties recruit various microorganisms to form its speci c rhizosphere microbiome that contributed to its growth and health [55]. These results were consistent with the differences between plant genotypes even in a single gene can contribute a signi cant impact on the rhizosphere microbiome [56].
Moreover, depth functional pro ling analysis revealed that the functional traits were signi cantly different in rhizosphere microbial communities of blueberry varieties and bulk soil. The increase in functional traits a liated with enzyme families, environmental adaptation, and biosynthesis of secondary metabolites were associated with the health of blueberry [57,58]. The phenotypic compositions of different blueberry varieties' rhizosphere microbial communities also exhibited signi cant differences. The proportions of stress tolerant of rhizosphere microbial communities of three blueberry varieties were higher than those of bulk soil, which suggested that the rhizosphere microbial composition contributed to different tolerance to stress tolerant for different blueberry varieties [59]. Overall, the differences in functional and phenotypic compositions of microbial communities between rhizosphere microbial communities of three blueberry varieties and bulk soil also suggested that different genotypes of blueberry hold their own unique microbiome, which contributes to their growth and health. The differences in taxonomical, functional, and phenotypic compositions of microbial communities between rhizosphere of blueberry varieties and bulk soil, even among different blueberry varieties, were determined by blueberry genotypes by actively secreting the compounds that speci cally stimulate or inhibit the members of the microbial community [60].
Besides, there are core microbiota among the rhizosphere microbial communities of blueberry. We identi ed 201 OTUs, which were mainly a liated with Proteobacteria, Actinobacteria, Acidobacteria, Firmicutes, Chloro exi, and Planctomycetes, as the core rhizosphere microbiota for blueberry rhizosphere microbial communities. Previous studies have reported that bene cial rhizosphere microbiota can directly affect the pathogen in the rhizosphere microbial community [60] and produce the antibiotic compounds and lytic enzymes, consumption of pathogen stimulatory compounds and competitions for nutrients for plants [61]. Among the core microbiota of blueberry, we identi ed two OTUs a liated with genus Pseudomonas as bene cial rhizosphere microbiota, because rhizosphere Pseudomonas spp. can produce the antifungal compound 2,4-diacetylphloroglucinol [62]. Moreover, the rhizosphere microbial compositions of three different blueberry varieties could be distinctly separated, and we selected 28 discriminative biomarkers to distinguish these three blueberry varieties.
Finally, the members of the co-occurrence network and their interactions between OTUs provide a deep understanding of the rhizosphere microbiome of blueberry from an ecological perspective. The members of these families of the rhizosphere microbial community contribute to the growth and health of plants.
For example, a previous study reported that the members of Xanthomonadaceae family can be divided into non-pathogenic and pathogenic species that infect humans and plants and these species have diverse effects on plant-related lifestyles [63]. The family Koribacteraceae of the Populus trichocarpa rhizosphere microbiome was reported to be correlated with the production of salicylic acid and populin [64]. Additionally, we observed that Acidbacterium, Salinibacterium, Micrococcus, and Conexibacter were involved in co-occurrence network (Fig. 6). Given the limitation of taxonomical classi cation, the members of these families of rhizosphere microbial communities of blueberry were unclear. Considering the high proportions of these families in co-occurrence network, we need to focus on the functions of these families in future research.

Conclusions
Our ndings highlighted the taxonomical, functional, and phenotypic compositions of the blueberry rhizosphere microbiome and demonstrated the differences of the rhizosphere microbiome in different blueberry varieties. As a result, our study provides an integrative view on the blueberry rhizosphere microbial community and identi es a series of taxa with potential importance from co-occurrence network. The separation of species of core rhizosphere microbiome, especially the bene cial microorganisms, including the non-pathogenic species a liated with genus Pseudomonas and family Xanthomonadaceae, could be used as potential microecologics and microbial fertilizers to maintain the health of blueberry during blueberry production. Given that rhizosphere microbiota harbor fungi and bacteria, and mycorrhizosphere interactions can improve plants tness and soil quality [65], the interactions between bacteria and fungi (especially mycorrhizal fungi) should be emphasized in further study. Our present work allows further investigation into the interactions between bacteria and fungi during blueberry production.