- Research article
- Open Access
Composition and diversity of bacterial communities in the rhizosphere of the Chinese medicinal herb Dendrobium
BMC Plant Biology volume 21, Article number: 127 (2021)
Dendrobium is a precious herbal that belongs to Orchidaceae and is widely used as health care traditional Chinese medicine in Asia. Although orchids are mycorrhizal plants, most research still focuses on endophytes, and there is still large amount unknown about rhizosphere microorganisms. To investigate the rhizosphere microbial community of different Dendrobium species during the maturity stage, we used high-throughput sequencing to analyze microbial community in rhizosphere soil during the maturity stage of three kinds of Dendrobium species.
In our study, a total of 240,320 sequences and 11,179 OTUs were obtained from these three Dendrobium species. According to the analysis of OTU annotation results, different Dendrobium rhizosphere soil bacteria include 2 kingdoms, 63 phyla, 72 classes, 159 orders, 309 families, 850 genera and 663 species. Among all sequences, the dominant bacterial phyla (relative abundance > 1%) were Proteobacteria, Actinobacteria, Bacteroidetes, Acidobacteria, Firmicutes, Verrucomicrobia, Planctomycetes, Chloroflexi, and Gemmatimonadetes. And through WGCNA analysis, we found the hub flora was also belong to Acidobacteria, Actinobacteria and Proteobacteria.
We found that the rhizosphere bacterial communities of the three kinds of Dendrobium have significant differences, and that the main species of rhizosphere microorganisms of Dendrobium are concentrated in the Proteobacteria, Actinobacteria, and Bacteroidetes. Moreover, the smaller the bacterial level, the greater the difference among Dendrobium species. These results fill knowledge gaps in the rhizosphere microbial community of Dendrobium and provide a theoretical basis for the subsequent mining of microbial functions and the study of biological fertilizers.
Rhizosphere soil refers to the narrow zone of soil affected by root exudations, containing up to 1011 microbial cells and over 30,000 prokaryotic species [1,2,3]. Rhizosphere microbes are significantly different from non-rhizosphere microbes in terms of species, number, and activity . There are abundant microbial resources in rhizosphere soil, which can be 10–1000 times that in bulk soil . The rhizosphere bacterial communities are represented by diverse bacterial taxa, though their abundance varies across the root system under different soil and plant types . In addition, the contents of the Firmicutes and Acidobacteria are also high . Using high-throughput sequencing, Lundberg et al. found that Arabidopsis thaliana had the highest relative abundance of Proteobacteria in the rhizosphere bacterial community, with Pseudomonadaceae being the main population, followed by Bacteroides and Actinomyces . The interaction between rhizosphere microbes and plants can affect material circulation and energy flow. Changes in its community structure and abundance can affect plant growth and development, flowering and fruiting, and plant interaction with phytophagous insects, which is of great significance for plant growth and yield [9, 10]. Dendrobium spp. is the second largest genus in Orchidaceae, with more than 1400 species in the world. There are 74 species and 2 varieties of Dendrobium in Orchidaceae in China, most of which are precious medicinal plants. Dendrobium has many functions, such as benefiting stomach and promoting body fluid, clearing heat and nourishing yin, relieving inflammation and relieving pain, clearing eyesight, and enhancing immunity . Currently, artificial cultivation using non-symbiotic tissue culture does not meet market demand because of slow growth and low survival rates . Therefore, the development of an effective method for propagating these endangered species for both conservation and commercial production is needed. There are three kinds of medicinal Dendrobium in the Ta-pieh Mountain area, which are Dendrobium huoshanense, Dendrobium officinale and Dendrobium moniliforme . All of them are perennial herbs of Dendrobium in Orchidaceae. They contain mainly polysaccharides, alkaloids and bibenzyls. They have antitumor, immunomodulatory, antioxidant, vasodilator and hypoglycemic effects .
At present, the research in Dendrobium mainly focuses on the endophytic bacteria. Indeed, there is a complex relationship between Dendrobium and its endophytes. Previous studies have shown that Sphiugomouas and Mycobacterium bacteria isolated from the roots of D. moschatum (Buch. - ham) SW. could significantly improve the seed germination rate of D. moschatum (Buch. - ham) SW ; and there are many studies using protocorm as material to successfully isolate and obtain effective fungi that promote seed germination [16,17,18]. However, few reports are available regarding Dendrobium and soil microbial communities. There are abundant bacteria, fungi and actinomycetes in the rhizosphere of medicinal plants. Not only can they promote the absorption of soil nutrients by medicinal plants, but they can also improve the yield and quality of medicinal plants; however, they can also cause continuous cropping obstacles for medicinal plants. So, studying Dendrobium medicinal plants and their rhizosphere microbes is of great significance to clarify how rhizosphere microbes can improve the yield and quality of medicinal plants. However, the traditional cultivation methods account for 0.1 to 1% of the environmental microorganisms, which cannot fully reflect the real situation of the environmental microbial community. In recent years, with the rapid development of new generation sequencing technology, based on 16 s rRNA sequence amplification and Illumina Miseq high-throughput sequencing technology, huge amounts of data can be obtained. By means of bioinformatics, the overall microbial community composition can be obtained. With the advantages of high-throughput, low price and short operation cycle, it has been widely used in the study of microbial community structure. Especially in the field of medicinal plants, such as Panax ginseng , Panax notoginseng , Ajuga bracteosa , Origanum vulgare , Lilium davidii  and other plants,these methods are widely used.
They always interact with other microorganisms such as mutualism, parasitism, antagonism. Co-culture is often used to study the interaction between microorganisms , but due to the inculcability of most microorganisms, co-culture has shown great limitations. Therefore, it is necessary to find other reliable ways to study the interaction between microbial communities. The rapid development of high-throughput sequencing has provided a new method for the study of microbial interaction. Based on high-throughput sequencing results, researchers have developed a number of software to predict the interaction among microbial communities, such as CoNet, MENA . These network analyses reflect the co-expression relationship among microbial communities. Weighted gene co-expression network analysis (WGCNA) is a kind of gene expression can be measured with a specific character combined with the analysis of network analysis algorithm . It can divide the gene co-expression network of complex biological processes into several highly correlated characteristic modules, and can establish the correlation between the modules and specific traits, so as to find the hub flora. WGCNA analysis has been widely used in transcriptional data analysis of various plants (such as corn , tomato  and Dioscorea nipponica ). But now, WGCNA analysis also has many applications in microbiome, such as rhizosphere microbes of potatoes .
Therefore, the goal of our work reported here was to characterize the rhizosphere microbial community of different Dendrobium species during the maturity stage. Our specific objective was to describe taxa associated with each Dendrobium species and determine which the environmental factors were related to microbial diversity and community composition. Specifically, we hypothesized that: (1) the dominant genus bacteria of Dendrobium and its comparison with other medicinal plants; (2) rhizosphere community composition will differ between three different Dendrobium; and (3) the hub flora can be obtained through WGCNA analysis.
Materials and methods
Plant material and soil sampling
Dendrobium plants were artificially cultivated in the greenhouse of Anhui Tongjisheng Biotechnology Company, Lu’an, China. The original source was collected by the company from the wild after obtaining local permission. The protocols and the conditions of planting for the growth of protocorm-like bodies were described by our previous study . From the grown plants of the three species of Dendrobium, two-year-old D. huoshanense, D. moniliforme and D. officinale were selected to provide eight replicates of each sample. The voucher specimens were authenticated by Professor Maoyun Yu from West Anhui University and deposited at Jiangsu Key Laboratory of Crop Genetics and Physiology in Yangzhou University, Yangzhou, China (Voucher number: 20A01, 20A02 and 20A03).
To obtain rhizosphere soil, plants were removed from flowerpots and large soil aggregates were removed by hand; soil firmly attached to roots was collected with sterile brushes and regarded as rhizosphere soil. The rhizosphere soil was sampled and sieved to remove plant debris. A part of the soil samples was put into sterile centrifuge tubes, frozen in liquid nitrogen immediately, and then stored at − 80 °C until the soil microbial composition was analyzed.
DNA extraction, PCR amplification and 16S rRNA sequencing
For each flowerpot we obtained a single DNA sample, which was used for 16S rRNA sequencing. Total DNA extraction and concentration were used DNeasy Power Soil Kit (Qiagen, Valencia, CA, USA), Qubit® dsDNA Assay Kit (Life Technologies, CA, USA) and DNA degradation degree and potential contamination were monitored on 1% agarose gels .
For bacterial diversity analysis, V4 hypervariable regions of 16S rRNA genes were amplified with universal primers 515 F (5′-GTGCCAGCMGCCGCGG-3′) and 806 R (5′-GGACTACHVGGGTWTCTAAT-3′) , and then PCR products were sequenced by IonS5™XL sequencing technique platform . Quality filtering on the raw reads were performed under specific filtering conditions to obtain the high-quality clean reads according to the Cutadapt (V1.9.1)  quality-controlled process. Chimeric sequences were removed by comparing with the Silva database  using UCHIME algorithm  to detect chimera sequences, and then the chimera sequences were removed. Operational taxonomic units (OTUs, cutoff 97% sequence identity for 16S rRNA) were clustered by UPARSE software (v7.0.1001) . The phylogenetic taxonomy was assigned according to the Ribosomal Database Project (RDP) classifier at an 80% confidence threshold (Version 2.2)  using the Silva databases for bacteria. Alpha-diversity was described for each sample using the following metrics for the observed species (OTU numbers); Chao1, Shannon and Simpson index, ACE and Good-coverage, which were generated to compare the level of bacterial OTU diversity .
Construction of weighted OTUs co-expression networks and identification of modules
The co-expression network of OTU was constructed using the WGCNA program package in R software . Co-expression module construction was according to Wan et al. methods . We used the function (pickSoftThreshold) to calculate the appropriate weighting coefficient β. The selection of the soft threshold should satisfy that the square of the correlation coefficient is close to 0.8 while ensuring that the connectivity cannot be too low. The WGCNA algorithm was used to construct the modules to extract the OTU information corresponding to each module after the soft threshold is determined.
Interactions analysis of co-expression modules
According to the soft threshold selected above, the adjacency matrix was calculated, and the topological overlap matrix (TOM) was further constructed. The value in the TOM matrix reflects the similarity of the co-expression relationship between two OTUs. The higher the similarity, the closer the value is to 1. Then, we used the function (flashClust) to perform hierarchical clustering and used the dynamic cut tree algorithm (minModuleSize = 25) to cut the gene cluster tree for module division . At the same time, colors were randomly assigned to each module, and principal component analysis was used to calculate the module feature vector of each module. The module feature vector was used to estimate the relationship between the module and the trait, the samples were classified according to the corresponding characteristics, and the module with the correlation coefficient ≥ |0.5| and p value ≤0.01 was selected for further analysis.
Identification of hub OTUs
Hub OTUs are often used as an abbreviation for highly connected OTU, which has a high degree of connectivity in the co-expression module. In our research, according to the size of the module, we classified the top 20 OTUs with the highest correlation as hub OTUs. At the same time, we used Cytoscape (v.3.6.1) to build and visualize the interaction network between OTU and OTU .
The WGCNA software package, stat and ggplot2 packages in R software (Version 2.15.3) were used for PCoA and WGCNA analysis. We conducted the co-occurrence network analysis in different rhizosphere soil compartments. The network analyses were performed using the psych R package (v.1.8.12) . The bipartite network was constructed to feature the shared phyla among different samples . The bipartite network was visualized using Cytoscape (v.3.6.1) .
α- and β-diversity of soil microbial communities
The rhizosphere soil alpha diversity indexes of different Dendrobium species are different (Fig. 1). The Good’s coverage index of Dendrobium library was greater than 98.5%, indicating that the sequencing results reflect the real situation of the bacterial population in the sample. We found that Dendrobium huoshanense (Dh) had the highest Shannon diversity index (10.29), while Dendrobium officinale (Do) had the lowest (9.16), which indicated that Dh rhizosphere soil has the highest diversity of bacterial communities. The index of ACE and Chao1 have the same trend, Dh is the highest, respectively 5225 and 4877, indicating that Dh samples have the highest community richness. Further analysis revealed that the Shannon index and Simpson index of soil bacteria of Dh was significantly different from Do, and Do was significantly different from Dendrobium moniliforme (Dm).
Principal coordinate analysis (PCoA) was performed at the operational taxonomic unit (OTU) level (Fig. 2). The analysis of PCoA using the weighted UniFrac distance, indicated a distinct pattern in the rhizosphere bacterial communities associated with the two axes, explaining 29.58 and 19.62% of the total variation in the Dendrobium rhizosphere soil. Using the unweighted UniFrac distance, the two axes explained 17.64 and 11.79% of the total variation in the Dendrobium rhizosphere soil. Bacterial communities associated with the rhizosphere soil were clustered in three regions according to a PCoA using the unweighted UniFrac distance, corresponding to Dh, Do, and Dm. However, no clustering was detected using the weighted UniFrac distance from different Dendrobium species.
Taxonomic classification and abundance
Rarefaction curves for bacterial communities suggested that changes in OTU density within the different Dendrobium species was sufficiently captured, and the sequencing was relatively comprehensive in covering the microbial communities (Fig. 3a). After quality filtering and processing according to a 97% similarity, a total of 240,320 sequences were obtained from these three Dendrobium species, respectively: 80,101 for Dh, 80,109 for Do and 80,110 for Dm. The OTUs common to the different Dendrobium species are presented in Fig. 3b as a Venn diagram to assess the relationships among the bacterial communities. 5609 OTUs were shared by Dh, Do and Dm; 1041 OTUs were shared by Dm and Do; 926 OTUs were shared by Dh and Do; and 971 OTUs were shared by Dm and Dh. The numbers of OTUs unique to each species were as follows: 1007 for Dm, 909 for Do, and 716 for Dh. A total of 11,179 OTUs were detected in all samples (Fig. 3b).
According to the analysis of OTU annotation results, different Dendrobium rhizosphere soil bacteria include 2 kingdoms, 63 phyla, 72 classes, 159 orders, 309 families, 850 genera and 663 species (Fig. 4). Among all sequences, the dominant bacterial phyla (relative abundance > 1%) were Proteobacteria, Actinobacteria, Bacteroidetes, Acidobacteria, Firmicutes, Verrucomicrobia, Planctomycetes, Chloroflexi, and Gemmatimonadetes, with contributions of 45.41, 17.58, 9.94, 7.77, 3.93, 3.03, 2.88, 1.98 and 1.14%, respectively. The difference at the phylum level of different Dendrobium species is not large, but there were more and more obvious differences in the level of the class afterwards. At the order level, the relative abundance of Alteromonadeles on Dh and Dm was low, while the relative abundance on Do was relatively high. Pseudoalteromonadaceae at the family level and Pseudoalteromonas at the genus level also have the same relative abundance pattern. At the species level, although the top ten species of bacteria in the rhizosphere soils of the three Dendrobium species were the same, the relative expression abundances were very different.
LEfSe analysis and co-occurrence patterns between rhizosphere soil microbiomes
To further elucidate the possible interactions between identified bacterial dependencies in rhizosphere soil samples, linear discriminant analysis (LDA) and effect size (LEfSe) methods were used for quantitative analysis of biomarkers in different species. We detected significant differences in the abundance of bacterial biomarkers from different groups and identified a total of 13 biomarkers from all rhizosphere soil samples as shown in the branching diagram (Fig. 5b). The significant taxa in the Dh were affiliated with diverse phylogenetic groups, including the family Sphingomonadaceae, and order Sphingomoadales. In the Dm, the significantly abundant taxa were the genus Flavobacterium, the family Flavobacteriaceae and the order Flavobacteriales. In the Do, the significant taxa belonged to the phylum Actinobacteria, class unidentified Actinobacteria and Acdobacteriia, orders Pseudonocardiales and Acidobacteriales, family Pseudonocardiaceae and unidentified Acidobacteriales, and genus Crossiella, which were all abundant.
Network analysis is widely used in the exploration of complex microbial communities. In the rhizosphere soil of Dendrobium, the distributions and abundances of microbial genera are significantly different. In order to further explore the interaction between microorganisms in the root micro-ecological environment of Dendrobium, the contributions of different microbial populations to the overall community structure were evaluated using a bipartite association network to visualize the associations between the phylum and rhizospheric compartments (Fig. 7). In the bacterial community, 7 clusters of phyla were generated in the three kinds of Dendrobium species of the bipartite network. In the middle cluster of the bacterial community, 37.5% (3687) of the indicator phyla were associated with all the Dendrobium species, which is a large proportion of the indicator phyla.
Construction of gene co-expression network
We used the WGCNA program package in R software to construct the co-expression network and selected 11,179 OTUs from the sequencing results of the rhizosphere soil microbiome of 24 Dendrobium samples to construct the co-expression network. To begin with, we loaded the R package, read the expression value matrix, calculated and clustered the correlation coefficients of the expression levels of each sample, and drew a heat map of species characteristics. As shown in Fig. 5a, most of the samples are clustered into one group according to species. The relative abundance of OTU is ranked as Dendrobium moniliforme (Dm) > Dendrobium officinale (Do) > Dendrobium huoshanense (Dh). In order to make the co-expression network meet the condition of non-scale distribution, we needed a suitable weighting coefficient β; therefore, the function (pickSoftThreshold) was used to calculate β. While satisfying that the square of the correlation coefficient is close to 0.8, it was ensured that the average connectivity could not be too low. According to these two conditions, we chose β = 26, and the square of the correlation coefficient was 0.85.
The soft threshold obtained above was used to construct a scale-free network. Firstly, the adjacency matrix was calculated based on the expression value matrix, and then the topological overlap matrix reflecting the common expression similarity was deduced. Then, the topological overlap matrix was used for hierarchical clustering to draw a cluster tree that could jointly express the overall distribution characteristics of similarity. Finally, the dynamic shear tree algorithm was used to cut the resulting cluster tree. In this process, OTUs with high similarity in common expression were gathered into the same branch, and different branches of the cluster tree represented different modules, and each module was given a specified color (Fig. 5b). Ten co-expression modules were obtained after constructing the co-expression network (Table 1). Among them, turquoise module had the largest number of OTUs, which is 501, while purple module had the smallest number of OTUs, which was 29.
Interaction analysis of co-expression module
In order to find the interaction relationship between co-expressed network modules, we analyzed the correlation of eigenvectors. Firstly, cluster analysis of eigenvectors was carried out (Fig. 6a). These 10 clusters were divided into two large clusters, including four modules (Brown, Black, turquoise and Red) and the other six modules. Therefore, the correlation between different modules was also significantly different, as shown in Fig. 6b. We analyzed the correlation between module feature vectors and species characteristics. As shown in Fig. 6c, four modules are highly correlated with species characteristics: turquoise module (p-value =7.5× 10− 5, correlation coefficient = − 0.731), Red module (p-value =0.0096, correlation coefficient = − 0.528), yellow module (p-value =0.00028, correlation coefficient = 0.688) and green module (p-value =0.00033, correlation coefficient = 0.682). Therefore, these four modules were further studied as modules of interest below.
Hub OTUs identification
From the above study, we obtained four modules with the most significant differences in species characteristics: turquoise, red, yellow and green. According to OTU notes (Additional file Table 1), we have carried out network analysis on these four modules (Fig. 7). In the turquoise module, there are 501 OTUs distributed in 19 phyla and 73 genera, among which the dominant bacteria phyla are Acidobacteria (13.2%), Actinobacteria (17.4%) and Proteobacteria (31.9%). In the red module, there are 76 OTUs distributed in 11 phyla and 22 genera, among which Proteobacteria (38.2%) and Actinobacteria (23.7%) are the dominant bacteria. In the yellow module, there are 134 OTUs distributed in 10 phyla and 34 genera, among which Proteobacteria (55.2%) is the dominant phylum. In the green module, there are 134 OTUs distributed in 11 phyla and 28 genera, among which the dominant bacteria phyla are Proteobacteria (51.9%) and Acidobacteria (19.8%).
Rhizosphere microorganisms can co-exist with plant roots, colonize and be maintained in roots, and play an important role in promoting plant growth and development . Among these microorganisms, the utilization rate and sensitivity of bacteria to root exudates are far higher than that of fungi, and bacteria are the most active and dominant microorganisms in rhizosphere [47, 48]. In addition, the number and species of rhizosphere microorganisms have a direct impact on soil biochemical activity and nutrient transformation. Under the influence of various complex factors of natural conditions, there are great differences in rhizosphere microbial flora of different plants and even different genotypes of the same plant [49, 50]. So, in this study, bacterial communities associated with the rhizosphere of three different Dendrobium species were characterized by high-throughput sequencing during the maturity stage.
Bacterial communities in rhizosphere soil of Dendrobium
It can be seen from the results that the dominant phyla in the rhizosphere soil of Dendrobium are Proteobacteria, Actinobacteria, Bacteroidetes, Acidobacteria, Firmicutes, Verrucomicrobia, Planctomycetes, Chloroflexi and Gemmatimonadetes. This indicates the richness and diversity of the rhizosphere microbial community of Dendrobium. However, the main dominant bacteria of Dendrobium are still concentrated in Proteobacteria, Actinobacteria and Bacteroidetes. Lundberg et al. found that the relative abundance of Proteobacteria, Bacteroidetes, Actinobacteria, Acidobacteria, Firmicutes and Gemmatimonadetes in Arabidopsis rhizosphere soil were relatively high, and were the dominant bacteria in Arabidopsis rhizosphere soil . Ling et al. studied the rhizosphere microorganisms of watermelon, and found that Acidobacteria, Actinomycetes, Bacteroidetes, Cyanobacteria, Firmicutes and Proteobacteria were the main leading bacteria . Davide et al. found that Actinobacteria, Bacteroidetes and Proteobacteria were dominant in barley rhizosphere soil . Although there are differences in the bacterial communities of different plant rhizosphere soils, the dominant phyla of Actinobacteria, Bacteroidetes and Proteobacteria are the common dominant phyla of Dendrobium and the above-mentioned plants, indicating that these phyla may be the common dominant phyla of all plant rhizosphere bacterial communities. Proteobacteria predominates in all ecosystems, especially in soil systems [53, 54], due to the fact that Proteobacteria contains a large level of physiological, morphological and metabolic diversity, and that Proteobacteria is of great significance to the C and N cycles . Proteobacteria reproduce fast, have good adaptability to unstable carbon sources, and are widely distributed in the global soil environment . In general, the abundance of Proteobacteria or Acidobacteria in soil samples is the greatest. These bacterial groups rich in rhizosphere microorganisms of Dendrobium, are also found to be the dominant communities of other plant rhizosphere microbiomes [57, 58].
Differences of bacterial communities in rhizosphere soils of three different Dendrobium species
According to the diversity indices alpha and beta, there were significant differences among the three Dendrobium species. The relative abundance of Acidobacteriales, Pseudonocardiaceae, Pseudoalteromonas and Pseudomonadales in Dendrobium officinale were higher than those in the other two species. Acidobacteriales is the dominant bacteria in the common plant rhizosphere bacterial community. Because Acidobacteriales can degrade complex root exudates such as cellulose and lignin, it plays a major role in the plant rhizosphere carbon cycle . Pseudoalteromonas sp. secretes a variety of extracellular active substances, including proteins, polysaccharides, brominated compounds, extracellular enzymes, extracellular toxins, antibiotics and so on . These substances have antibacterial, algicidal, bactericidal and cellulose degrading activities [61, 62]. Pseudomonadales is an important group of biocontrol microorganisms, and is also one of the most widely distributed microorganisms in nature. Its rapid reproduction, strong colonization ability and simple nutrition requirements have been widely studied for its inhibition of plant diseases and promotion of plant growth. In the existing studies, Pseudonocardiaceae is mainly associated with cellulose degradation and antibiotic synthesis [63,64,65]. It can be seen from the above that the function of the bacterial community in the rhizosphere of Dendrobium officinale may be stronger than that of the other two kinds of Dendrobium. In Dendrobium moniliforme, the relative abundance of Bacteroidetes is relatively high, and Bacteroidetes is a poor nutrient bacterium, which is suitable for growth in the environment with less absorbable nutrients such as organic matter and available nitrogen .
In bacterial species, we found that Lysobacter soli have a high relative abundance in Dendrobium moniliforme. This fungus has previously been isolated from the soil where ginseng is grown , and has been found in other plants to promote plant activity . Psychrobacter ibarius has been found in Dendrobium officinale. Some researchers have isolated this bacterium from the root plane of Angelica sinensis . This is a fungus related to polysaccharide synthesis. This is consistent with our previous research results . Our previous research also found that among these three kinds of Dendrobium, the polysaccharides of two-year-old Dendrobium are higher than those in the other two kinds of two-year-old Dendrobium.
WGCNA applied to the rhizosphere microbiome in Dendrobium species
In recent years, with the development of high-throughput sequencing technology, network analysis based on high-throughput omics data has become a popular big data processing method due to its high efficiency and convenience . WGCNA, as one of the methods of network analysis, can specifically screen OTUs highly correlated with the target traits, divide them into multiple co-expression modules, and mine the core OTUs, which has been proved to be a fast, accurate and efficient method for data mining in potatoes and other plants [71, 72]. In order to understand the differences in the rhizosphere microbial communities of Dendrobium, this study used the WGCNA package in the R software to construct a co-expression network of the rhizosphere microbial community of Dendrobium, and explored the community structure of the rhizosphere microbes in a new way. A total of 10 modules were identified in the constructed Dendrobium rhizosphere microbial co-expression network, and the 4 modules with the highest correlation with species characteristics were further analyzed. Each module obtained 20 hub OTUs. By annotating these hub OTUs, the dominant bacteria phyla in the rhizosphere microbial community of Proteobacteria, Acidobacteria and Actinobacteria were finally obtained.
Compared with the results of this study, a total of 22 co-expression modules were identified in the rhizospheric co-expression network in the WGCNA analysis of rice root-associated microbiomes. By analyzing the OTUs annotations of the three microbial co-expression networks with the most significant differences among the three niches within the endosphere, rhizoplane and rhizosphere, it is concluded that the hub bacteria phyla in the rice rhizosphere microbial community is Proteobacteria, Actinobacteria and Bacteroidetes, which are highly similar to the dominant bacteria phyla we found in the rhizosphere microbial community of Dendrobium. The dominant microorganisms are relatively abundant in the soil and play an important role in the regulation of ecological functions . Proteobacteria and Actinomycota are enriched in rhizosphere and non-rhizosphere environments, which is related to their strong adaptability. Among them, the Proteobacteria can survive and multiply in most environments and become the dominant bacteria because its outer membrane is mainly composed of lipopolysaccharide, which can protect its internal genetic material from external interference [74, 75]. In addition, the phylum Proteobacteria plays an important role in adapting to environmental changes and resisting adversity stress . Actinomycota plays an important role in the decolorization of the toxic dye-triphenylmethane (TPM) in the soil , and plays a certain role in adaptation to arid environments . Therefore, we speculate that the Proteobacteria and Actinomycota play a pivotal role in the plant rhizosphere microbial community, and they play an important regulatory role in the plant rhizosphere.
Overall, the chemical properties of rhizosphere soil of Dendrobium and the climatic characteristics of Dendrobium cultivation were firstly analyzed. Then, combined with high-throughput sequencing technology, the rhizosphere soil microbial communities of different species of Dendrobium were studied, and the rhizosphere microorganisms and environmental driving factors of different species of Dendrobium were discussed. The main conclusions are as follows: the dominant bacteria in rhizosphere soil of Dendrobium are as follows: Proteobacteria, Bacteroidetes, Actinobacteria, and Acidobacteria. These are the dominant bacteria in the rhizosphere bacterial community. However, there are some differences in the bacterial communities among different Dendrobium species, and the smaller the bacterial level is, the greater the difference is. Through WGCNA analysis, it is found that in the rhizosphere microbial community network of Dendrobium, the hub flora is also in the Bacteroides, Acidobacteria and Proteobacteria.
Availability of data and materials
Raw amplicon sequence data related to this study were deposited in the NCBI Sequence Read Archive (NCBI SRA) under Bioprojects PRJNA638443.
Operational taxonomic units
Ribosomal Database Project
Principal coordinate analysis
Egamberdieva D, Kamilova F, Validov S, Gafurova L, Kucharova Z, Lugtenberg B. High incidence of plant growth-stimulating bacteria associated with the rhizosphere of wheat grown on salinated soil in Uzbekistan. Environ Microbiol. 2008;10(1):1–9.
Mendes R, Kruijt M, De Bruijn I, Dekkers E, van der Voort M, Schneider JH, Piceno YM, DeSantis TZ, Andersen GL, Bakker PA. Deciphering the rhizosphere microbiome for disease-suppressive bacteria. Science. 2011;332(6033):1097–100.
Pangesti N, Pineda A, Pieterse CM, Dicke M, Van Loon JJ. Two-way plant mediated interactions between root-associated microbes and insects: from ecology to mechanisms. Front Plant Sci. 2013;4(414):414.
Hinsinger P, Marschner P. Rhizosphere--perspectives and challenges--a tribute to Lorenz Hiltner 12-17 September 2004—Munich, Germany. Plant Soil. 2006;283(1–2):7–8.
Miransari M. Soil microbes and plant fertilization. Appl Microbiol Biotechnol. 2011;92(5):875–85.
Sharifi R, Ryu C-M. Chatting with a tiny belowground member of the holobiome: communication between plants and growth-promoting rhizobacteria. Adv Bot Res. 2017;82: Elsevier:135–60.
Lareen A, Burton F, Schäfer P. Plant root-microbe communication in shaping root microbiomes. Plant Mol Biol. 2016;90(6):575–87.
Lundberg DS, Lebeis SL, Paredes SH, Yourstone S, Gehring J, Malfatti S, Tremblay J, Engelbrektson A, Kunin V, Del Rio TG. Defining the core Arabidopsis thaliana root microbiome. Nature. 2012;488(7409):86–90.
Schnitzer SA, Klironomos JN, HilleRisLambers J, Kinkel LL, Reich PB, Xiao K, Rillig MC, Sikes BA, Callaway RM, Mangan SA. Soil microbes drive the classic plant diversity–productivity pattern. Ecology. 2011;92(2):296–303.
Wagg C, Jansa J, Schmid B, van der Heijden MG. Belowground biodiversity effects of plant symbionts support aboveground productivity. Ecol Lett. 2011;14(10):1001–9.
Pharmacopoeia CN. Pharmacopoeia of PR China. Beijing: Press of Chemical Industry; 2005.
Chen J, Wang H, Guo S-X. Isolation and identification of endophytic and mycorrhizal fungi from seeds and roots of Dendrobium (Orchidaceae). Mycorrhiza. 2012;22(4):297–307.
Yuan Y, Yu M, Zhang B, Liu X, Zhang J. Comparative nutritional characteristics of the three major Chinese Dendrobium species with different growth years. PLoS One. 2019;14(9):e0222666.
Chen N-D, You T, Li J, Bai L-T, Hao J-W, Xu X-Y. A comparative study of three tissue-cultured Dendrobium species and their wild correspondences by headspace gas chromatography–mass spectrometry combined with chemometric methods. J Food Drug Anal. 2016;24(4):839–47.
Tsavkelova EA, Cherdyntseva TA, Klimova SY, Shestakov AI, Botina SG, Netrusov AI. Orchid-associated bacteria produce indole-3-acetic acid, promote seed germination, and increase their microbial yield in response to exogenous auxin. Arch Microbiol. 2007;188(6):655–64.
Huang H, Zi X-M, Lin H. Gao J-YJjom: host-specificity of symbiotic mycorrhizal fungi for enhancing seed germination, protocorm formation and seedling development of over-collected medicinal orchid, Dendrobium devonianum. J Microbiol. 2018;56(1):42–8.
Meng Y-Y, Shao S-C, Liu S-J, Gao J-Y. Do the fungi associated with roots of adult plants support seed germination? A case study on Dendrobium exile (Orchidaceae). Glob Ecol Conserv. 2019;17:e00582.
Zi X-M, Sheng C-L, Goodale UM, Shao S-C, Gao J-Y. In situ seed baiting to isolate germination-enhancing fungi for an epiphytic orchid, Dendrobium aphyllum (Orchidaceae). Mycorrhiza. 2014;24(7):487–99.
Ying YX, Ding WL, Li Y. Characterization of soil bacterial communities in rhizospheric and nonrhizospheric soil of Panax ginseng. Biochem Genet. 2012;50(11–12):848–59.
Wu Z, Hao Z, Zeng Y, Guo L, Huang L, Chen B. Molecular characterization of microbial communities in the rhizosphere soils and roots of diseased and healthy Panax notoginseng. Antonie Van Leeuwenhoek. 2015;108(5):1059–74.
Kumar G, Kanaujia N, Bafana A. Functional and phylogenetic diversity of root-associated bacteria of Ajuga bracteosa in Kangra valley. Microbiol Res. 2012;167(4):220–5.
Bafana A. Diversity and metabolic potential of culturable root-associated bacteria from Origanum vulgare in sub-Himalayan region. World J Microbiol Biotechnol. 2013;29(1):63–74.
Shang Q, Yang G, Wang Y, Wu X, Zhao X, Hao H, Li Y, Xie Z, Zhang Y, Wang R. Illumina-based analysis of the rhizosphere microbial communities associated with healthy and wilted Lanzhou lily (Lilium davidii var. unicolor) plants grown in the field. World J Microbiol Biotechnol. 2016;32(6):95.
Liu Z, Wei H, Zhang J, Saleem M, He J, Zhong J, Ma RJM. Higher sensitivity of Microbial Network Than Community Structure under Acid Rain. Microorganisms 2021;9(1):118.
Wang H, Wei Z, Mei L, Gu J, Yin S, Faust K, Raes J, Deng Y, Wang Y, Shen QJSB et al: Combined use of network inference tools identifies ecologically meaningful bacterial associations in a paddy soil. 2017.
Weijing W, Wenjie J, Lin H, Haiping D, Yili W, Chunsheng X, Qihua T, Shuxia L, Dongfeng ZJBG. Weighted gene co-expression network analysis of expression data of monozygotic twins identifies specific modules and hub genes related to BMI. BMC Genomics. 2017;18(1):872.
Wang J, Zhang Y, Pan X, Du J, Guo X. Discovery of leaf region and time point related modules and genes in maize (Zea mays L.) leaves by Weighted Gene Co-expression Network analysis (WGCNA) of gene expression profiles of carbon metabolism. J Integr Agric. 2019;18(2):350–60.
DiLeo MV, Strahan GD, den Bakker M, Hoekenga OA. Weighted correlation network analysis (WGCNA) applied to the tomato fruit metabolome. PLoS One. 2011;6(10):e26683.
Sun W, Wang B, Yang J, Wang W, Liu A, Leng L, Xiang L, Song C, Chen S. Weighted gene co-expression network analysis of the Dioscin rich medicinal plant Dioscorea nipponica. Front Plant Sci. 2017;8:789.
Castillo JD, Vivanco JM, DKJME M. Bacterial Microbiome and Nematode Occurrence in Different Potato Agricultural Soils; 2017.
Yuan Y, Zhang J, Kallman J, Liu X, Meng M, Lin J. Polysaccharide biosynthetic pathway profiling and putative gene mining of Dendrobium moniliforme using RNA-Seq in different tissues. BMC Plant Biol. 2019;19(1):1–19.
Chen QL, An XL, Zhu YG, Su JQ, Gillings MR, Ye ZL, Cui L. Application of Struvite alters the antibiotic Resistome in soil, Rhizosphere, and Phyllosphere. Environ Sci Technol. 2017;51(14):8149–57.
Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Huntley J, Fierer N, Owens SM, Betley J, Fraser L, Bauer M. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012;6(8):1621–4.
Li Y, Wang P, Yin J, Jin S, Su W, Tian J, Li T, Yao K. Effects of ornithine alpha-ketoglutarate on growth performance and gut microbiota in a chronic oxidative stress pig model induced by d-galactose. Food Funct. 2020;11(1):472–82.
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17(1):10–2.
Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, Peplies J, Glöckner FO. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2012;41(D1):D590–6.
Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27(16):2194–200.
Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013;10(10):996–8.
Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73(16):5261–7.
Qiao Q, Wang F, Zhang J, Chen Y, Zhang C, Liu G, Zhang H, Ma C, Zhang J. The variation in the Rhizosphere microbiome of cotton with soil type, Genotype and Developmental Stage. Sci Rep. 2017;7(1):3940.
Amrine KCH, Blanco-Ulate B, Cantu D, Alberto DLF. Discovery of Core Biotic Stress Responsive Genes in Arabidopsis by Weighted Gene Co-Expression Network Analysis. J Plos One. 2015;10(3):e0118731.
Wan Q, Tang J, Han Y, Wang D. Co-expression modules construction by WGCNA and identify potential prognostic markers of uveal melanoma. Exp Eye Res. 2018;166:13–20.
Shannon P, Markiel A, Ozier O, et al. Cytoscape: A Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 2003;13(11):2498–504.
Ren Y, Xun W, Yan H, Ma A, Xiong W, Shen Q, Zhang R. Functional compensation dominates the assembly of plant rhizospheric bacterial community. Soil Biol Biochem. 2020;150:107968.
Jiao S, Zhang Z, Yang F, Lin Y, Chen W, Wei G. Temporal dynamics of microbial communities in microcosms in response to pollutants. Mol Ecol. 2017;26(3):923–36.
Prashar P, Kapoor N, Sachdeva S. Rhizosphere: its structure, bacterial diversity and significance. Rev Environ Sci Biotechnol. 2014;13(1):63–77.
Lu Y, Rosencrantz D, Liesack W, Conrad R. Structure and activity of bacterial community inhabiting rice roots and the rhizosphere. Environ Microbiol. 2006;8(8):1351–60.
Singh BK, Munro S, Potts JM, Millard P. Influence of grass species and soil type on rhizosphere microbial community structure in grassland soils. Appl Soil Ecol. 2007;36(2–3):147–55.
Lemanceau P, Corberand T, Gardan L, Latour X, Laguerre G, Boeufgras J, Alabouvette C. Effect of two plant species, flax (Linum usitatissinum L.) and tomato (Lycopersicon esculentum mill.), on the diversity of soilborne populations of fluorescent pseudomonads. Appl Environ Microbiol. 1995;61(3):1004–12.
Wiehe W, Höflich G. Survival of plant growth promoting rhizosphere bacteria in the rhizosphere of different crops and migration to non-inoculated plants under field conditions in north-East Germany. Microbiol Res. 1995;150(2):201–6.
Ling N, Song Y, Raza W, Huang Q, Guo S, Shen Q. The response of root-associated bacterial community to the grafting of watermelon. Plant Soil. 2015;391(1–2):253–64.
Bulgarelli D, Garrido-Oter R, Münch PC, Weiman A, Dröge J, Pan Y, McHardy AC, Schulze-Lefert P. Structure and function of the bacterial root microbiota in wild and domesticated barley. Cell Host Microbe. 2015;17(3):392–403.
Dai Z, Su W, Chen H, Barberán A, Zhao H, Yu M, Yu L, Brookes PC, Schadt CW, Chang SX. Long-term nitrogen fertilization decreases bacterial diversity and favors the growth of Actinobacteria and Proteobacteria in agro-ecosystems across the globe. Glob Chang Biol. 2018;24(8):3452–61.
Lin Y-T, Whitman WB, Coleman DC, Jien S-H, Chiu C-Y. Cedar and bamboo plantations alter structure and diversity of the soil bacterial community from a hardwood forest in subtropical mountain. Appl Soil Ecol. 2017;112:28–33.
Kersters K, De Vos P, Gillis M, Swings J, Vandamme P, Stackebrant E. Introduction to the Proteobacteria. In: Dworkin M, Falkow S, Rosenberg E, Schleifer KH, Stackebrandt E, editors. The Prokaryotes, vol. 5. 3rd ed. Springer: New York; 2006.
Xu J, Zhang Y, Zhang P, Trivedi P, Riera N, Wang Y, Liu X, Fan G, Tang J, Coletta-Filho HD. The structure and function of the global citrus rhizosphere microbiome. Nat Commun. 2018;9(1):1–10.
Mendes LW, Kuramae EE, Navarrete AA, Van Veen JA, Tsai SM. Taxonomical and functional microbial community selection in soybean rhizosphere. ISME J. 2014;8(8):1577–87.
Spain AM, Krumholz LR, Elshahed MS. Abundance, composition, diversity and novelty of soil Proteobacteria. ISME J. 2009;3(8):992–1000.
Ward NL, Challacombe JF, Janssen PH, Henrissat B, Kuske CR. Three genomes from the phylum Acidobacteria provide insight into the lifestyles of these microorganisms in soils. Appl Environ Microbiol. 2009;75(7):2046–56.
Bowman JP. Bioactive compound synthetic capacity and ecological significance of marine bacterial genus Pseudoalteromonas. Marine Drugs. 2007;5(4):220–41.
Tutino ML, Parrilli E, Giaquinto L, Duilio A, Sannia G, Feller G, Marino G. Secretion of α-amylase from Pseudoalteromonas haloplanktis TAB23: two different pathways in different hosts. J Bacteriol. 2002;184(20):5814–7.
Sanchez-Porro C, Mellado E, Bertoldo C, Antranikian G, Ventosa A. Screening and characterization of the protease CP1 produced by the moderately halophilic bacterium Pseudoalteromonas sp. strain CP76. Extremophiles. 2003;7(3):221–8.
Carr G, Derbyshire ER, Caldera E, Currie CR, Clardy J. Antibiotic and antimalarial quinones from fungus-growing ant-associated Pseudonocardia sp. J Nat Prod. 2012;75(10):1806–9.
Cuesta G, Soler A, Alonso J, Ruvira M, Lucena T, Arahal D, Goodfellow M. Pseudonocardia hispaniensis sp. nov., a novel actinomycete isolated from industrial wastewater activated sludge. Antonie Van Leeuwenhoek. 2013;103(1):135–42.
Kaewkla O, Franco CM. Pseudonocardia adelaidensis sp. nov., an endophytic actinobacterium isolated from the surface-sterilized stem of a grey box tree (Eucalyptus microcarpa). Int J Syst Evol Microbiol. 2010;60(12):2818–22.
Wong JTF, Chen X, Deng W, Chai Y, Ng CWW, Wong MH. Effects of biochar on bacterial communities in a newly established landfill cover topsoil. J Environ Manag. 2019;236:667–73.
Srinivasan S, Kim MK, Sathiyaraj G, Kim H-B, Kim Y-J, Yang D-C. Lysobacter soli sp. nov., isolated from soil of a ginseng field. Int J Syst Evol Microbiol. 2010;60(7):1543–7.
Lasudee K, Tokuyama S, Lumyong S, Pathom-Aree W. Mycorrhizal spores associated Lysobacter soli and its plant growth promoting activity. Chiang Mai J Sci. 2017;44:94–101.
Han S-I, Lee H-J, Lee H-R, Kim K-K, Whang K-S. Mucilaginibacter polysacchareus sp. nov., an exopolysaccharide-producing bacterial species isolated from the rhizoplane of the herb Angelica sinensis. Int J Syst Evol Microbiol. 2012;62(3):632–7.
Juan MAY-YC, Li-Feng WANG, Jing-Jing LI, Hao WANG, Yan-Ping FAN, Hui-Yong LI. Identification of gene co-expression modules of maize plant height and ear height by WGCNA. Acta Agron Sin. 2020;46(3):385–94.
Tianyuan Q, Bai JP. Identification of drought-related co-expression modules and hub genes in potato roots based on WGCNA. Acta Agron Sin. 2020;46(7):1033–1051.
Zhao W, Langfelder P, Fuller T, Dong J. Hovarth SJJoBS: weighted gene coexpression network analysis: state of the art. J Biopharm Stat. 2010;20(2):281–300.
Delgado-Baquerizo M, Grinyer J, Reich PB, Singh BK. Relative importance of soil properties and microbial community for soil functionality: insights from a microbial swap experiment. Functional Ecology. 2016;30(11):1862–73.
Song HG, Kim OS, Yoo JJ, Jeon SO, Ahn TS. Monitoring of soil bacterial community and some inoculated Bacteria after prescribed fire in microcosm. J Microbiol. 2004;42(4):285–91.
Whalley WR, Riseley B, Leeds-Harrison PB, Bird NRA. Adderley WPJEJoSS: structural differences between bulk and rhizosphere soil. Eur J Soil Sci. 2010;56(3):353–60.
Campbell B J, Engel A S, Porter M L, et al. The versatile ε-proteobacteria: key players in sulphidic habitats. Nat Rev Microbiol. 2006;4(6):458–468.
Adenan NH, Lim YY, Ting ASY. Discovering Decolorization Potential of Triphenylmethane Dyes by Actinobacteria from Soil. Water Air Soil Pollut. 2020;231(12):560.
Mohammadipanah F, Wink J. Actinobacteria from arid and desert habitats: diversity and biological activity. Front Microbiol. 2015;6:1541.
We would thank Prof. Donald L. DeAngelis, from the Wetland and Aquatic Research Center, U.S. Geological Survey, for his valuable comments and suggestions that have greatly improved the quality of this manuscript.
This project was funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD), the National traditional Chinese medicine standardization project (ZYYS-2014 ). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Ethics approval and consent to participate
Dendrobium huoshanense, Dendrobium officinale and Dendrobium moniliforme used in this study were cultivated by Maoyun Yu’s greenhouse from Anhui Tongjisheng Biotechnology Company, Lu’an, China. No permits were required for the collection of the samples.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zuo, J., Zu, M., Liu, L. et al. Composition and diversity of bacterial communities in the rhizosphere of the Chinese medicinal herb Dendrobium. BMC Plant Biol 21, 127 (2021). https://doi.org/10.1186/s12870-021-02893-y
- Microbial community
- 16 s rRNA
- Rhizosphere soil