Next Article in Journal
The Molecular Basis of Ubiquitin-Conjugating Enzymes (E2s) as a Potential Target for Cancer Therapy
Next Article in Special Issue
Catalase (CAT) Gene Family in Rapeseed (Brassica napus L.): Genome-Wide Analysis, Identification, and Expression Pattern in Response to Multiple Hormones and Abiotic Stress Conditions
Previous Article in Journal
Genotyping by Sequencing Reveals Genetic Relatedness of Southwestern U.S. Blue Maize Landraces
Previous Article in Special Issue
Genome Wide Analysis of U-Box E3 Ubiquitin Ligases in Wheat (Triticum aestivum L.)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Degeneration Gradient of Poplar Trees Contributes to the Taxonomic, Functional, and Resistome Diversity of Bacterial Communities in Rhizosphere Soils

1
Department of Microbiology, College of Biological Sciences and Biotechnology, Beijing Forestry University, Beijing 100083, China
2
Beijing Advanced Innovation Center for Tree Breeding by Molecular Design, Beijing Forestry University, Beijing 100083, China
*
Author to whom correspondence should be addressed.
J.L. and X.H. contributed equally to this work.
Int. J. Mol. Sci. 2021, 22(7), 3438; https://doi.org/10.3390/ijms22073438
Submission received: 1 March 2021 / Revised: 22 March 2021 / Accepted: 23 March 2021 / Published: 26 March 2021
(This article belongs to the Special Issue Genomic Studies of Plant-Environment Interactions)

Abstract

:
Bacterial communities associated with roots influence the health and nutrition of the host plant. However, the microbiome discrepancy are not well understood under different healthy conditions. Here, we tested the hypothesis that rhizosphere soil microbial diversity and function varies along a degeneration gradient of poplar, with a focus on plant growth promoting bacteria (PGPB) and antibiotic resistance genes. Comprehensive metagenomic analysis including taxonomic investigation, functional detection, and ARG (antibiotics resistance genes) annotation revealed that available potassium (AK) was correlated with microbial diversity and function. We proposed several microbes, Bradyrhizobium, Sphingomonas, Mesorhizobium, Nocardioides, Variovorax, Gemmatimonadetes, Rhizobacter, Pedosphaera, Candidatus Solibacter, Acidobacterium, and Phenylobacterium, as candidates to reflect the soil fertility and the plant health. The highest abundance of multidrug resistance genes and the four mainly microbial resistance mechanisms (antibiotic efflux, antibiotic target protection, antibiotic target alteration, and antibiotic target replacement) in healthy poplar rhizosphere, corroborated the relationship between soil fertility and microbial activity. This result suggested that healthy rhizosphere soil harbored microbes with a higher capacity and had more complex microbial interaction network to promote plant growing and reduce intracellular levels of antibiotics. Our findings suggested a correlation between the plant degeneration gradient and bacterial communities, and provided insight into the role of high-turnover microbial communities as well as potential PGPB as real-time indicators of forestry soil quality, and demonstrated the inner interaction contributed by the bacterial communities.

Graphical Abstract

1. Introduction

The rhizosphere soil influenced by root secretions contains up to one hundred billion microbial cells per gram root [1] and about thirty thousand prokaryotic species [2]. As the second genome of the plant, the complex root-associated microbial community is crucial for nutrition and plant health [3]. A variety of organic compounds including carbohydrate, organic acids, phenol, and other substances secreted by plant roots in the surrounding soil, which are the sources of nutrients and vitamins for various bacteria, stimulating the development of soil microbial communities and changing their activity and ecological distribution [2,4]. Microbes are in competition for nutrients and space with many other microbes in the rhizosphere to affect them to a certain extent through a complex network of interactions. The beneficial microbiome associated with roots, including the so-called plant growth promoting (PGP) bacteria, can contribute to alleviate plant stress by a variety of mechanisms [5]. Among them, PGP bacteria can directly enhance micronutrient uptake and affect phytohormone homeostasis, or indirectly stimulate the plant immune system against phytopathogens [6] and improve soil texture and structure [5,7]. Recent evidence suggests that differences among plant species when grown on the same soil or between plant genotypes have a significant impact on the rhizosphere microbiome [8,9]. However, we were unable to formulate how root microbiota change under different plant phenotype, especially artificial vegetation under natural conditions.
As a woody perennial organism, Populus has become the model for insight into the mechanistic hypotheses related to plant–microbe interactions because of its fast growth rates and full genome sequenced. Several studies have demonstrated that microbial isolates from poplars can enhance the health, growth, development and resistance to biotic and abiotic stresses of their plant hosts [10,11,12]. The core microbiome and niche-level microbiome of poplar were dissected and archaeal/bacterial diversity from rhizosphere soil was richer than that in leaf, stem and root regardless of plant genotype [13]. The composition of poplar rhizosphere microbiome was significantly different across environmental gradients or genotypes [14,15]. However, little is known that the assembly and shift of the microbial consortia within or between Populus phenotypes.
Populus are not only used in paper industry, as a cellulose derived fuel raw materials, but also usually in artificial afforestation to protect soil erosion. Restoration of vegetation is undoubtedly an important measure to reduce air pollution, prevent wind and fix sand, and cope with desertification and global warming [16,17,18]. Many countries and regions have carried out natural forest protection and artificial forest cultivation programs [19], among which the Three North Shelterbelt known as “Green Wall of China” as one of the typical examples. During the decades of the construction of the Three-North Shelterbelt, many studies have made evaluations based on its achievements, and the extensive ecological benefits obtained are beyond doubt [20,21,22]. However, in recent years, with the impact of climate change and human activities, large scale Populus trees degeneration with dieback symptom and partial aggravating trend have been observed in some areas of the artificial forest [23]. Especially in Zhangbei County of Hebei Province, nearly one third of the Populus trees, as the main species in the local artificial forests, exhibit rapid widespread degeneration [24].
Many studies have been put forward to explain the natural factors caused the degradation of the Three North Shelterbelt from the aspects of afforestation methods, soil physical and chemical properties, climate and environmental conditions and tree physiology [19,25,26,27]. Changes in temperature and precipitation exacerbate the frequency and intensity of drought, as well as the timing mismatch between resource demand and local and regional resource supply, increasing plant stress and mortality according to the survey of artificial forests [25,26]. Poplar bacterial canker disease and insufficient soil nutrients were also associated with degeneration and mortality [27,28,29,30]. However, few researches have directly investigated the genetic profile of poplar root microbial community, and there is a knowledge gap in understanding taxonomic, functional, and resistome diversity of rhizosphere microorganisms of poplar under different health conditions. Here, we hypothesized that rhizosphere soil microbial diversity and function varies along a degeneration gradient of poplar. In this study, based on Illumina PE150 metagenome sequencing platform, the microbial community structure, functional, and resistance genes of Populus rhizosphere soils were comparatively analyzed for undegenerated, medium-degenerated, and completely degenerated Populus, respectively (Figure 1a). These results in comparative analysis of the microbiome composition will lay a foundation for further revealing the interaction mechanism between Populus and its rhizosphere microorganisms, and also provide a theoretical basis for rational cultivation of Populus as an ecological tree species.

2. Results and Discussion

2.1. Populus and Rhizosphere Soil Characterization

Three representative plots (C, M and H) were defined based on the different degeneration degrees and the growth characterization of Populus were carried out (Figure 1a). The tree features of H plot were not detected, because the dead trees were recognized to have been cut down and cleared by local authorities, leaving only the stumps. The average diameter at breast height (DBH) and tree height of poplar trees were 15.46 cm and 21.49 m in C plot, but only 11.56 and 18.30 m in M plots, respectively. In M plot, more than 92% of the trees experienced treetop dryness to varying degrees, and the average height of the dry part was 6.82 m. The overall canopy density of C was 0.83, that of M was only 0.54 (Table S1) because of the dieback of the trees. In addition, five-point sampling (1 m × 1 m) was carried out using the diagonal principle in the C, M, and H plots to understand the species richness in each plot (Table S2). The average species coverage of C, M, and H was 17.2%, 71.8%, and 79.4%, respectively. As treetops dry up, a lower canopy density in medium-degraded plot to receive more sunlight, resulting in richer plant species.
Previous research suggested that soil characteristics not only could reflect soil productivity, but also played important roles in shaping the soil bacterial community diversity and structure in the terrestrial ecosystem [31]. Trees healthy had a strong significant effect on all the soil physical and chemical properties. Nine soil chemical characteristics were selected and measured to compare the differences among C, M, and H group (Figure 1b). The content of total nitrogen (TN), total phosphorus (TP), microbial biomass carbon (MBC), microbial biomass nitrogen (MBN), available phosphorus (AP), available potassium (AK), and organic matter (OM) were highest in sample C with 0.166%, 0.143%, 405.7 mg/kg, 77.2 mg/kg, 8.97 mg/kg, 282.94 mg/kg, 36.62 g/kg and markedly (p < 0.05 or p < 0.01) higher than in sample H, which were similar between H and M soils. All the pH were close to neutral, but presented a differing variation: H > M > C samples (rough trend) with 6.72, 6.81, and 7.06, respectively. It is worth noting that the content of total potassium (TK) was highest with 2.21%, but AK was lowest with 119.7 mg/kg in H, respectively. In summary, the soil characteristics with different Populus growth conditions in the same area exhibited obvious differences, especially between C and H plus M. The soil fertility was on full displayed by the physiological and biochemical characteristics, which perhaps lead to the growth distinction of the Populus.

2.2. Sequencing Results

The 18 soil samples subjected to metagenomic sequencing generated more than 101 Gb of raw data (average of 5.6 Gb per sample). After quality control and removal of poplarsequence (host contamination 0.663–0.06%), filtered sequencing data per sample reduced to 97.35% of the average raw data. The assembly of the reads using MEGAHIT software generated an average of 636,794, 573,556 and 504,433 scaftigs in C, M, and H, respectively (Table S3). The lengths of scaftigs N50 were between 646 bp to 850 bp, and those of scaftigs N90 were between 520 bp to 539 bp. These scaftigs (≥500 bp) were clustered into gene catalogs (unigenes) after ORF (Open Reading Frame) prediction and filtering (≥100 nt). 9,660,445 non-redundant ORFs were generated from metagenomic libraries, of which 0.21 million (15.1%) were complete ORFs containing both starting and ending codons. 49.89% of the unigenes can be assigned to taxonomic annotations, which suggesting that there was a huge resource pool for novel functional gene in rhizosphere soil of the shelterbelt poplar. About 99.86% of the total annotated unigenes were prokaryotic (bacteria and archaea). Only 0.06% and 0.01% of the unigenes were annotated as eukaryotic (including fungi, protozoa, algae, and plants) and virus, respectively, indicating that the eukaryotic and viral communities have been seriously underestimated due to potential annotation bias. Several studies found generally high levels of fungal diversity in association with transgenic poplar roots [32,33,34]. However, in this study, eukaryotic abundance was only 0.03%, indicating a unique feature of microbial community in poplar rhizosphere of shelter forest.
The percentage of bacterial was significantly higher in samples from C than those from M and H (p = 0.0012, p = 0.00023), whereas the proportion of archaea increased from sample C to sample H (p = 0.012, p = 0.0011). Three type samples shared 5,530,921 genes, and 161,741 and 108,437 genes were found exclusively in samples H and M, respectively. By contrast, the sample C contained about six times as many peculiar genes as the other two samples (Figure S1a), which suggested that the sample C is annotated to a richer microbial community. The correlation analysis among different samples revealed six parallel samples within each group have higher similarity, and distinct differences between types (Figure S1b).

2.3. Comparative Analysis of Microbial Communities

Bray–Curtis distance matrix was used for sample clustering analysis to construct the total prokaryotic profile among samples at different taxonomic levels from phylum to genus. A total of 160 phyla were recovered from all samples. The most abundant bacteria phyla were Proteobacteria (36.75% ± 6.85%), Acidobacteria (22.2% ± 5.7%), and Actinobacteria (21.9% ± 2.9%) (Figure 2a). The five most abundant prokaryotic families were Sphinogomonadaceae, Bradyrhizobiaceae, Phyllobacteriaceae, Solirubrobacteraceae, and Pyrinomonadaceae, which comprised on average 33.24% of the bacteria communities (Figure 2b). At the genus level, Bradyrhizobium, Sphingomonas, and Solirubrobacter were the most three abundant across all communities (Figure 2c).
By comparing the relative abundance of prokaryotes among the three sampling sites on the above-mentioned taxonomic level, revealed the great difference on microbial communities between C and H groups. Among the three sample groups, Proteobacteria is the most abundant of all the phylum showing significant differences, and the abundance from C to M and H significantly decreased (q = 0.00018, q = 0.0042, Hypothesis Testing), but there was no significant difference between M and H (Figure 2d). The main families, including Bradyrhizobiaceae, Sphingomonadaceae, Phyllobacteriaceae, Nocardioidaceae, Comamonadaceae, Caulobacteraceae, Gemmatimonadaceae, Chitinophagaceae, Burkholderiaceae, and Xanthomonadaceae, and main genera, including Bradyrhizobium, Sphingomonas, Mesorhizobium, Nocardioides, Variovorax, Gemmatimonadetes, Rhizobacter, Pedosphaera, Candidatus Solibacter, Acidobacterium, and Phenylobacterium, showed the similar trend (Figure 2e,f). Moreover, Mesouhizobium, Rhodoplanes, Luteitalea, and Streptomyces were also the microbes with obviously higher abundance in C than in H groups (q < 0.05), which were similar in H and M groups. In addition, the abundance of Thermoleophilum at phylum level, Rubrobacteraceae at family level, Thermoleophilum and Rubrobacter at genus level decreased from H to C (q < 0.01).
In terms of prokaryotes, Proteobacteria, Actinobacteria and Acidobacteria are the dominant phyla, as well as the rhizosphere of many other hosts [35,36]. At the genus level, Bradyrhizobium and Mesorhizobium contains species widely known as nitrogen-fixing bacteria with legumes, and they have an outstanding ability to promote the growth of non-legumes through the production of indole-3-acetic acid (IAA) and siderophores, the solubilization of potassium and phosphate [37]. Members of the Sphingomonas genus increased Arabidopsis thaliana growth rate, alter rhizosphere microbial community structure of plant under drought stress [38], and alleviate Cd stress in oilseed rape through regulation of the GSH-AsA cycle and antioxidative enzymes [39]. Variovorax genus has the potential to promote plant growth, due to its ability to hydrolyze cellulose [40] and to produce plant growth enzymes and substances such as siderophores, IAA and 1-aminocyclopropane-1-carboxylic acid deaminase [41,42,43]. The high relative abundance of Verrucomicrobia was found in grassland soils of a coniferous forest [44], in the rhizosphere of Arabidopsis thaliana [45], poplar [46] and maize [47,48], and in rice roots [49,50,51]. Verrucomicrobia has been evidenced that play critical roles in environmental carbon cycling and (poly) saccharide degradation because of high coding densities for glycoside hydrolase genes in the genome [52,53]. Memebers of Verrucomicrobia colonized in intercellular spaces and inside root cells of rice, and produced IAA and its precursors IAM (indole-acetamide) or IPA (indole-3-pyruvic acid) for the plant growth promotion, showed activity of acid phosphatases for the solubilization of organic phosphates, and possessed also genes for citrate synthase [45]. Members of the Sphingomonas and Variovorax genera can survive in low nutrient environment, often be isolated from oil-contaminated soil due to their unique wide range of xenobiotic-biodegradative abilities [54,55]. Streptomyces bacteria, famous for producing secondary metabolites including antibiotics, may be used as plant probiotics due to the beneficial effects on plants and are not only ubiquitous in soil, but also are isolated from plant roots [56]. In summary, multiple members enriched genera in the healthy poplar rhizosphere are well known as beneficial plant microorganisms, which can maintain the hormonal balance, promote the acquisition of nutrients and development of the root system in its plant host, and thus prevent plant diseases.
It is worth noting that only two genera Thermoleophilum and Rubrobacter have higher content in H than C. Thermoleophilum, Rubrobacter are thermophilic and radiation-resistant strains [57], which may be the reason why they can survive in a relatively poor rhizosphere environment. Studies have shown that the same plant host can recruit similar microorganisms to form a rhizosphere environment on different soils [58]. It can be seen that there is little significant difference in microbial diversity among all the samples, and dominant microbes are widely found in nature. However, compared with H and M samples, C group has a more abundant microbiome. This may indicate a denser material exchange between poplars rhizosphere and host in C group. Importantly, how to link them with drought resistance and disease resistance.

2.4. Correlation Analysis of Soil Characterization and Microbial Communities

Within a host species, habitat and soil type, rather than host genetic background, have larger effects on the overall structure of the microbiome, but the balance of the effects of genetic and soil factors within host habitats on bacteria is less clear [15,47,59,60]. According to principal component analysis (PCA) and the redundancy analysis (RDA) based on the physicochemical properties of the soil and the most abundant prokaryotes among the samples, the three groups were separated obviously at the level of phylum, family and genus. At the level of phylum, Bacteroidetes, Proteobacteria, Verrucomicrobia, Actinobacteria, and Gemmatinonadetes were significantly positively correlated with OM and AK, while negatively correlated with Candidatus Tectomicrobia and Candidatus Rokubacteria (p < 0.05) (Figure 3a). At the family level, Sphingomonadaceae, Comamonadaceae, Nocardioidaceae, Hyphomicrobiaceae, Caulobacteraceae, Bradyrhizobiaceae, Phyllobacteriaceae, and Mycobacteriaceae were significant positive correlation with AK and OM (Figure 3b). At genus level, Nocardioides, Variovorax, Sphingomonas, Phenylobacterium, Rhodoplanes, Mesorhizobium, Bradyrhizobium, Streptomyces, Mycobacterium, and Luteitalea were positive correlated with AK and TN (Figure 3c). Pyrinomonadaceae at family level, and Rubrobacter and Pyrinomonas at genus level were significant positive correlation with pH (Figure 3b,c). Almost all the members positively associated with AK, OM, and TN were enriched in group C, but the members positively associated with pH were enriched in H, which suggested that pH, AK, OM, and TN were the indicators to identify the poplar rhizosphere soil fertility. In particular, AK was the best indicator to distinguish community distribution of group C from the other two sample sat three taxonomic levels (Figure 3).
Potassium is the third essential nutrient factor of plants, which plays a key role during growth, development and metabolism of plants. If potassium supply is insufficient, plants will be more susceptive to diseases and pest with slow-growing, small seeds, poorly developed roots, and lower yields [61,62,63]. Potassium-soluble bacteria in soil can convert insoluble or mineral potassium compounds into soluble forms for use by plants as soil solutions [64]. A variety of organic acids produced by potassium-solubilizing bacteria, such as tartaric acid, lactic, acetic, propionic, citric, oxalic, glycolic, succinic acid, fumaric, malonic and tartaric, mobilized and solubilized the insoluble potassium and structural unavailable forms of potassium compounds [65]. The maintenance of soil quality is a key factor for environmental sustainability, while the increase of microbial diversity with various functions, such as organic matter decomposition, potassium, and phosphorus solubilization, and nitrogen fixation etc., significantly affecting soil quality and sustainability [66]. Our results suggest that changes in bacterial community composition and diversity might be related to the AK influenced by the shifts in biogeochemical processes, which generally limit plants growth and fitness.

2.5. Comparison of Microbial Interaction Networks

To identify potential microbe–microbe interactions in poplar rhizosphere under different health conditions, we construct three bacterial co-occurrence networks (Figure 4). The network topological metrics showed that microbial co-occurrence patterns were markedly different among C, M, and H. The average degree (avgK) and network density values (Table S4) were similarity in C and H, but higher than that in M. Compared with the H group, the co-occurrence network established from C and M samples was more compact and complex, which was intuitively shown in the Figure 4. The network of group C was more homogeneous, while positive links in group M were significantly increased, which may be due to the dense positive correlation of genera belonging to Actinobacteria with each other. These observations suggested that more intimate host-microbe and microbe–microbe interactions occur in the healthy poplar rhizosphere than in the degradated poplar rhizosphere.
Although low abundance genera may contribute significantly to the function of the whole community, genera with high abundance and more connectivity links obviously play a critical role in the structural stability of the entire microbial network, considered as key species. Sphingomonas belong to Proteobacteria, Rubrobacter and Solirubrobacter belong to Actinobacteria, Pyrinomonas belong to Acidobacteria were found to be the key genus of Group C, M, and H, respectively.

2.6. Function Potential Analysis of Microbes

Functional annotations showed approximately 43% of the ORFs were assigned to the KOs by blasting against the KEGG Orthology (KO) database. At the first level of KEGG targeted the metabolism as the main functional pathway (54.24%), followed by genetic information processing (15.36%), environmental information processing (11.78%), cellular processes (9.41%), organismal systems (5.65%), and the human diseases (3.55%) (Figure S2). The 7108 KOs were mainly involved in 405 KEGG pathway at level 3 in total and the top 35 pathways accounted for about 56% of all metabolic pathways.
Plant–microbe and microbe–microbe interactions are very likely to be important factors that would influence the assembly of rhizosphere microbiomes. The samples C and H were enriched and compared in different metabolic pathways (Figure 5a). The KOs involved in known plant–microbe and microbe–microbe interactions, such as quorum sensing and two-component system, were over-represented in C (Figure 5a). These observations suggested that more intimate host-microbe and microbe–microbe interactions occur in the healthy poplar rhizosphere than in the degradated rhizosphere. Tryptophan is an important precursor of auxin biosynthesis. Tryptophan metabolism was enriched in C, indicating that the bacteria in healthy poplar rhizosphere had outstanding ability to promote plant growth. Nutrition in the rhizosphere and surrounding environment of host plants plays a key role in the construction and stabilization of the microbiome. Plant-derived compounds in the rhizosphere soil would be likely nutrient sources for microbes. Consistent with this, ABC transporters responsible for transporting plant-derived nutrients, such as peptide, amino acid, vitamin B12, monosaccharide, oligosaccharide, polysaccharide, lipid, mineral, organic ion, and iron-siderophore, into microbial cells, were over-represented in C group. Degradation of poplar trees may result in reduced secretions and more difficulty for microorganisms to obtain nutrients. Consistent with this, the KOs that include carbon metabolism, synthesis of coenzymes, and nucleic acid metabolism, such as, oxidative phosphorylation, glyoxylate and dicarboxylate metabolism, carbon fixation pathways in prokaryotes, pyruvate metabolism, citrate cycle (TCA cycle), glycolysis/gluconeogenesis, propanoate metabolism, pantothenate and CoA biosynthesis, butanoate metabolism, purine metabolism, and pyrimidine metabolism, were over-represented in H. Such findings were consistent with the fact that healthy poplar rhizosphere microorganisms could acquire a variety of simple nitrogen and carbon sources from root exudation, and therefore would not need to invest in their biosynthesis. The enrichment of pathways involved in the amino acid synthesis and metabolism, such as valine, leucine, isoleucine, arginine, aminoacyl-tRNA biosynthesis, and ribosome, alanine, aspartate, glutamate, glycine, serine, threonine, cysteine, and methionine metabolism, was observed in H. The amino acids recycling depleted in C further suggested that rhizosphere microbes can directly obtain amino acid from root exudates, and that amino acids represent an important nitrogen and/or carbon source. In addition, the KOs involved in adverse variation repair, such as RNA degradation, nucleotide excision repair, homologous recombination and mismatch repair, were enriched in H. These repair pathways and RNA degradation are important ways for cell to respond to environmental changes.
Subsequently, the Spearman correlation coefficient was calculated to reveal the contribution of microbiome to metabolic pathways based on the 17 most abundant genera and pathway profiles with significant differences between groups (Figure 5b). The results show the following multiple important genera, Gemmatirosa, Variovorax, Sphingomonas, Phenylobacterium, Bradyrhizobium, Mesorhizobium, Candidatus Solibacter, Rhizobacter, Rhodoplanes, Nocardioides, Mycobacterium, and Streptomyces, were main driving factors of ABC transporters and two component system. In addition, Acidobacterium and Pedosphaera also drove two component system. Most of above genera, were enriched in large quantities, and ABC transporters and two-component systems were also enriched in C, which indicate that the host–microbe and microbe–microbe interactions were more intimate in healthy poplar rhizosphere than in degradation poplar rhizosphere. The KOs involved in carbon and energy metabolism, amino acid cycling, nucleic acid repair, and homologous recombination were driven by Rubrobacter and Thermoleophilum belong to Actinobacteria. These KOs and the two genera were enriched in the H, which indicated that microbes required stronger basal metabolism to survive with a nutrient-deficient environment in degradation poplar rhizosphere.

2.7. Resistance Genome

In order to understand whether ARGs in the rhizosphere of poplars affect its adaptation to the environment and the prevention of pathogen infection, we explored the prevalence of ARGs among different samples. Compared with the group H and M, the group C contained a higher number of resistance genes (Figure S3a). More than 90% of the genes in each sample were not annotated in any ARO, indicating that poplar rhizosphere soil was a very huge drug resistance gene pool with great potential exploitation value. By blasting against the CARD database, there were 576 ARGs in common among a total of 633 ARGs annotated in all samples (Figure S3b). The top 30 abundant ARGs accounted for about 64.4% of all the annotated ARGs and the top 10 abundant resistance genes detected were macB, streptomyces_cinnamoneus_EF-Tu, tetA48, mfd, patA, and sav1866, Streptomyces_rishiriensis_parY, TaeA, arlR, and evgs. Macrolides-related ARGs (macB) are the most abundant and common in distribution, and the proportion in each sample exceeded 7% (Figure 6a). The ARGs associated with fluoroquinolone, macrolide, tetracycline, aminocoumarin, elfamycin, and pleuromutilin were also found to be in high abundance (Figure 6b,e). Interesting, the abundance of multi-drug resistance genes in all ARGs accounted for the highest proportion, exceeding 26% in this study among all the samples (Figure 6c,e). Moreover, the multidrug resistance genes, including Evgs, adeL, and emrB in group C significantly increased than those in group H and M (Figure 6b). Currently, multidrug resistance genes have been found in various bacteria, such as Staphylococcus, Bacillus, and Enterococcus. Microorganisms tend to develop multidrug resistance, as an economical powerful generation strategy, to overcome environmental pressures. We learned from the local area that no artificial antibiotics on the poplar samples. As a reservoir of ARGs, soil has obvious overlap with clinically resistant genomes, but the factors that affect the composition of ARGs in soil and their movement between genomes and habitats are largely unknown. All detected ARGs mainly provide microbial resistance through four resistance mechanisms: Antibiotic efflux, antibiotic target alteration, antibiotic target protection, and antibiotic target replacement (Figure 6c). The ARGs were assigned to the three dominant phyla, Proteobacteria, Actinobacteria and Acidobacteria, which from the latter two had the same proportions in the three groups, but which from the former in group C (22%) had higher proportions than in group M (20%) and H (19%) (Figure S4). Moreover, the abundance of ARGs from Sphingomonas, Bradyrhizobium, Mesorhizobium, Phenylobacterium, Variovorax, Rhizobacter, and Sphingopyxis belonging to Proteobacteria, Nocardioides and Kribbella belonging to Actinobacteria, and Luteitalea belonging to Acidobacteria, were significant higher in group C than those in other two groups. However, the abundance of ARGs from Rubrobacter belonging to Actinobacteria were increased from C, M to H (Figure 6d). It was previously reported that Actinobacteria and Proteobacteria were the most abundant potential hosts of multi-drug resistance genes [27], which is also confirmed by our results.
According to statistics based on the abundance of ARGs, antibiotic efflux, antibiotic target alteration and antibiotic target protection are the four main mechanisms that confer host microbial resistance [67]. The C group samples contained more abundant resistance genes, especially more abundant multidrug resistance genes, conferring resistance to more different drug classes, which were directly related to the higher abundance of the microbial community.

3. Materials and Methods

3.1. Sampling Sites and Sample Collection

Soil samples were collected in 2018 during the early August from the poplar-dominated shelterbelts in Zhangbei County, Hebei Province, China, which is a mid-temperate continental monsoon climate with average temperatures and annual precipitation of 12.6 °C and 410 mm, respectively. The three representative plots (20 m × 20 m), which were undegenerated (C, 114.871, 142.00° E, 41.339, 331.99° N), medium-level degenerated (M, 114.872, 047.84° E, 41.335, 539.52° N) and high-level degenerated plot (H, 114.867, 196.47° E, 41.334, 461.14° N), were defined based on the different degeneration degrees (Figure 1). High-level degradation means the trees lost vitality, the leaves fell off and the branches dried entirely, the root systems cannot get any nutrient from the soil (Figure 1a-H). Medium-level degradation means the trees lost a part of vitality, the treetops were dried entirely but the trunk base still alive, there are branches and leaves in the trunk base (Figure 1a-M). Undegradation means the trees live healthy and less dried from the branches or treetops (Figure 1a-C). Six poplar trees were randomly selected from each plot, and the rhizosphere soils were collected at a depth of 10–30 cm using the five-point sampling method. Soils of the same tree were thoroughly mixed as one sample, sieved through 425 μm mesh, placed in a sterilized bag, immediately placed on ice, and transported to the laboratory, where it was stored at −80 °C until physicochemical parameters measurement and DNA extraction.

3.2. Degeneration Data Determination and Soil Characterization

In order to better describe the growth of poplars, diameter at breast height (DBH), tree height, degenerated portion height, and canopy density of all poplars in medium-degenerated (M) and undegenerated (C) quadrat were measured. The physicochemical parameters of the soil were measured as follows. Soil pH was determined using a glass electrode meter (Sartorius, Göttingen, Niedersachsen, Germany) in a suspension of 10 g of soil in 25 mL of distilled water. Soil available phosphorus (AP) was extracted using sodium bicarbonate and then measured by the molybdenum blue method. Soil available potassium (AK) was determined by flame photometry [68,69]. Available nitrogen (AN) was determined by potassium persulfate oxidation. Organic matter (OM) content was determined as described by Walkey and Black [70,71]. Soil total potassium (TK), total phosphorus (TP), total nitrogen (TN), microbiological biomass nitrogen (MBN), and microbiological biomass carbon (MBC) were determined analyzed as previously described [72].

3.3. DNA Extraction, Shotgun Sequencing, and Metagenome Assembly

DNA from the soil samples was extracted using the PowerSoil DNA Isolation kit (Cellgeno, Beijing, China) following the manufacturer’s instructions. DNA quality was monitored on 1% agarose gels and concentration was measured using Qubit® dsDNA Assay Kit in Qubit® 2.0 Flurometer (Life Technologies, Carlsbad, CA, USA). Sequencing libraries were generated using NEBNext® Ultra™ DNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA) according to the manufacturer’s recommendations and index codes were added to attribute sequences to each sample [73]. The library preparation was then sequenced on IlluminaHiSeq PE150 platform at Novogene Company (Beijing, China) with paired-end reads generated. Raw sequencing data was preprocessed to remove low-quality and the host-originated reads to acquire the clean data for subsequent analysis. The clean data has been uploaded into the Sequence Read Archive of NCBI under BioProject PRJNA563959 with the accession number SRP241077.
Clean data of 18 samples were assembled by MEGAHIT (EBI-metagenomic, https://github.com/voutcn/megahit, accessed on 20 October 2018) with parameters“-presets meta-large (—min-count 2—k-min 27—k-max 87—k-step 10)” [74] to generate scaffolds. Assembled scaffolds were broken from N connection and obtained the scaftigs which were then filtered the fragment shorter than 500 bp for statistical analysis [75].

3.4. Taxonomic Profiling

Scaftigs (≥500 bp) from each sample were used for open reading frame (ORF) prediction [76,77] through MetaGeneMark (http://exon.gatech.edu/meta_gmhmmp.cgi, accessed on 5 November 2018). For ORF (≤100 nt) predicted, CD-HIT software (https://github.com/weizhongli/cdhit, accessed on 10 November 2018) was adopted to remove redundancy and obtain the unique initial gene catalogue. The clean reads from each sample successfully mapped to the gene catalogue, and the genes that contain reads ≤2 in all 18 soil samples were filtered to obtain non-redundant gene catalogue (unigenes) for further analysis. Based on the number of mapped reads and the length of gene, the abundance information of genes in each sample were counted.
To generate the taxonomic information of the unigenes, DIAMOND software (downloaded from https://github.com/bbuchfink/diamond, accessed on 15 December 2018) was taken to blast the unigenes to the sequences of bacteria, archaea, fungi and virus extracted from the NCBI microNR database (blastp, e-value ≤ 1 × 10−5) [78]. Based on the lowest common ancestor (LCA) algorithm implemented in MEGAN software (downloaded from https://www.wsi.uni-tuebingen.de/lehrstuehle/algorithms-in-bioinformatics/software/megan6/, accessed on 25 December 2018), the taxonomic annotation of the unigenes was assigned and the relative abundances of each taxa was calculated.

3.5. Network Construction

Networks were constructed for C, M, and H based on the relative abundance of the top 100 genera, account for about 80% of all annotated genera, respectively. The adjacency matrixes (threshold filtering: Correlation coefficient r > 0.7; significant p < 0.05) were achieved by using rcorr function (R: Hmisc) based on Spearman correlation [79] and visualized by using GePHI (version 0.9.2, downloaded from https://gephi.org/, accessed on 10 February 2019).

3.6. Function Annotations

Functional annotations of unigenes were performed through against the Kyoto Encyclopedia of Genes and Genomes (KEGG) database by DIAMOND (blastp, e-value ≤ 1 × 10−5). In brief, according to predefined sets in the KEGG database, the identified KO genes were further annotated into different pathways. The relative abundance of different functional hierarchy and the gene number table of each sample in each taxonomy hierarchy were obtained.

3.7. Resistance Gene Annotation

The unigenes were aligned to CARD database (https://card.mcmaster.ca/home, accessed on 30 January 2019) by Resistance Gene Identifier (RGI) built in CARD with the parameter setting blastp, e-value ≤ 1 × 10−30 [80]. Based on the aligned result, abundance distribution of the resistance genes in each samples, the species attribution analysis and the resistance mechanism of resistance genes analysis were also conducted.

3.8. Data Analysis

Six replicate samples each group were analyzed. Unweighted pair-group method with arithmetic mean (UPGMA) clustering based on Bray-Curtis distances were performed to investigate beta-diversity patterns. Metastats method [81] was applied to perform hypothesis testing on species abundance and p-values were adjusted to q-values by false discovery rate (FDR) correction to screen species with significant differences. Canoco (version 5.1, Wageningen University & Research, Wageningen, The Netherlands) was used to analyze the effects of soil physicochemical factors on species with Pearson correlation. Metabolic pathways and resistance gene heat maps were converted to Z values based on relative abundance, and the maps were presented using Origin (version 9.1, OriginLab, Northampton, MA, USA).

4. Conclusions

Through our analysis in this paper, the microbiome composition and functionality in rhizosphere of poplar are likely influenced by the health condition of poplar and available potassium (AK) of soil. We identified several microbes, Bradyrhizobium, Sphingomonas, Mesorhizobium, Nocardioides, Variovorax, Gemmatimonadetes, Rhizobacter, Pedosphaera, Candidatus Solibacter, Acidobacterium, and Phenylobacterium, may serve as useful bio-indicators of soil fertility. The highest abundance of multidrug resistance genes and the four mainly microbial resistance mechanisms (antibiotic efflux, antibiotic target protection, antibiotic target alteration and antibiotic target replacement) in healthy poplar rhizosphere corroborated the relationship between soil fertility and microbial activity. The healthy rhizosphere soil harbored microbes with a higher capacity and had more complex microbial interaction network to promote plant growing and reduce intracellular levels of antibiotics. The relationships among rhizosphere soil microbial communities, soil fertility and the host poplar presented here set the stage for a valuable reference database for protection and management policies of forestry. Therefore, we are striving for the continuous sampling of these as well as additional sites in the Three North Shelterbelt to provide new insight into the interactions between microbial communities and the health conditions of poplar.

Supplementary Materials

Supplementary Materials are available online at https://www.mdpi.com/article/10.3390/ijms22073438/s1.

Author Contributions

Y.M. designed and revised the manuscript. J.L. and X.H. performed the bioinformatics analysis and wrote the paper, with assistance from Y.M. Sampling work was carried out by J.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by Fundamental Research Funds for the Central Universities (2019ZY31); the Fundamental Research Funds for the Central Universities (2017ZY14) and the Forestry Industry Research Special Funds for Public Welfare Projects (201304409); National Natural Science Foundation of China (31770110).

Data Availability Statement

The sequencing data has been uploaded into the Sequence Read Archive of NCBI under BioProject PRJNA563959 with the accession number SRP241077 (https://www.ncbi.nlm.nih.gov/sra/?term=SRP241077, accessed on 20 February 2021).

Acknowledgments

We thank Jianbo Xie for the guidance of sequence analysis.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. 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–9. [Google Scholar] [CrossRef]
  2. Mendes, R.; Kruijt, M.; de Bruijn, I.; Dekkers, E.; van der Voort, M.; Schneider, J.H.M.; Piceno, Y.M.; DeSantis, T.Z.; Andersen, G.L.; Bakker, P.A.H.M.; et al. Deciphering the Rhizosphere Microbiome for Disease-Suppressive Bacteria. Science 2011, 332, 1097–1100. [Google Scholar] [CrossRef] [PubMed]
  3. Edwards, J.A.; Santos-Medellin, C.M.; Liechty, Z.S.; Bao, N.; Lurie, E.; Eason, S.; Phillips, G.; Sundaresan, V. Compositional shifts in root-associated bacterial and archaeal microbiota track the plant life cycle in field-grown rice. PLoS Biol. 2018, 16, e2003862. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Bakker, M.G.; Schlatter, D.C.; Otto-Hanson, L.; Kinkel, L.L. Diffuse symbioses: Roles of plant-plant, plant-microbe and microbe-microbe interactions in structuring the soil microbiome. Mol. Ecol. 2014, 23, 1571–1583. [Google Scholar] [CrossRef]
  5. Rolli, E.; Marasco, R.; Vigani, G.; Ettoumi, B.; Mapelli, F.; Deangelis, M.L.; Gandolfi, C.; Casati, E.; Previtali, F.; Gerbino, R.; et al. Improved plant resistance to drought is promoted by the root-associated microbiome as a water stress-dependent trait. Environ. Microbiol. 2015, 17, 316–331. [Google Scholar] [CrossRef]
  6. Castrillo, G.; Teixeira, P.J.P.L.; Paredes, S.H.; Law, T.F.; de Lorenzo, L.; Feltcher, M.E.; Finkel, O.M.; Breakfield, N.W.; Mieczkowski, P.; Jones, C.D.; et al. Root microbiota drive direct integration of phosphate stress and immunity. Nature 2017, 543, 513–518. [Google Scholar] [CrossRef]
  7. Toju, H.; Peay, K.G.; Yamamichi, M.; Narisawa, K.; Hiruma, K.; Naito, K.; Fukuda, S.; Ushio, M.; Nakaoka, S.; Onoda, Y.; et al. Core microbiomes for sustainable agroecosystems. Nat. Plants 2018, 4, 733. [Google Scholar] [CrossRef]
  8. Wei, Z.; Gu, Y.; Friman, V.-P.; Kowalchuk, G.A.; Xu, Y.; Shen, Q.; Jousset, A. Initial soil microbiome composition and functioning predetermine future plant health. Sci. Adv. 2019, 5, eaaw0759. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Chen, T.; Nomura, K.; Wang, X.; Sohrabi, R.; Xu, J.; Yao, L.; Paasch, B.C.; Ma, L.; Kremer, J.; Cheng, Y.; et al. A plant genetic network for preventing dysbiosis in the phyllosphere. Nature 2020, 580, 653–657. [Google Scholar] [CrossRef]
  10. Timm, C.M.; Pelletier, D.A.; Jawdy, S.S.; Gunter, L.E.; Henning, J.A.; Engle, N.; Aufrecht, J.; Gee, E.; Nookaew, I.; Yang, Z.; et al. Two Poplar-Associated Bacterial Isolates Induce Additive Favorable Responses in a Constructed Plant-Microbiome System. Front. Plant Sci. 2016, 7, 497. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Henning, J.A.; Weston, D.J.; Pelletier, D.A.; Timm, C.M.; Jawdy, S.S.; Classen, A.T. Root bacterial endophytes alter plant phenotype, but not physiology. Peerj 2016, 4, e2606. [Google Scholar] [CrossRef] [PubMed]
  12. Taghavi, S.; Garafola, C.; Monchy, S.; Newman, L.; Hoffman, A.; Weyens, N.; Barac, T.; Vangronsveld, J.; van der Lelie, D. Genome Survey and Characterization of Endophytic Bacteria Exhibiting a Beneficial Effect on Growth and Development of Poplar Trees. Appl. Environ. Microbiol. 2009, 75, 748–757. [Google Scholar] [CrossRef] [Green Version]
  13. Cregger, M.A.; Veach, A.M.; Yang, Z.K.; Crouch, M.J.; Vilgalys, R.; Tuskan, G.A.; Schadt, C.W. The Populus holobiont: Dissecting the effects of plant niches and genotype on the microbiome. Microbiome 2018, 6, 31. [Google Scholar] [CrossRef] [PubMed]
  14. Gottel, N.R.; Castro, H.F.; Kerley, M.; Yang, Z.M.; Pelletier, D.A.; Podar, M.; Karpinets, T.; Uberbacher, E.; Tuskan, G.A.; Vilgalys, R.; et al. Distinct Microbial Communities within the Endosphere and Rhizosphere of Populus deltoides Roots across Contrasting Soil Types. Appl. Environ. Microbiol. 2011, 77, 5934–5944. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Shakya, M.; Gottel, N.; Castro, H.; Yang, Z.K.; Gunter, L.; Labbe, J.; Muchero, W.; Bonito, G.; Vilgalys, R.; Tuskan, G.; et al. A Multifactor Analysis of Fungal and Bacterial Community Structure in the Root Microbiome of Mature Populus deltoides Trees. PLoS ONE 2013, 8, e76382. [Google Scholar] [CrossRef]
  16. Wei, W.; Wang, B.; Niu, X. Soil Erosion Reduction by Grain for Green Project in Desertification Areas of Northern China. Forests 2020, 11, 473. [Google Scholar] [CrossRef]
  17. Deng, L.; Kim, D.-G.; Li, M.; Huang, C.; Liu, Q.; Cheng, M.; Shangguan, Z.; Peng, C. Land-use changes driven by ‘Grain for Green’ program reduced carbon loss induced by soil erosion on the Loess Plateau of China. Glob. Planet. Chang. 2019, 177, 101–115. [Google Scholar] [CrossRef]
  18. Gong, G.; Liu, J.; Shao, Q.; Zhai, J. Sand-Fixing Function under the Change of Vegetation Coverage in a Wind Erosion Area in Northern China. J. Resour. Ecol. 2014, 5, 105–114. [Google Scholar] [CrossRef]
  19. Huang, L.; Wang, B.; Niu, X.; Gao, P.; Song, Q. Changes in ecosystem services and an analysis of driving factors for China’s Natural Forest Conservation Program. Ecol. Evol. 2019, 9, 3700–3716. [Google Scholar] [CrossRef]
  20. Cao, S.; Suo, X.; Xia, C. Payoff from afforestation under the Three-North Shelter Forest Program. J. Clean. Prod. 2020, 256. [Google Scholar] [CrossRef]
  21. Qiu, B.; Chen, G.; Tang, Z.; Lu, D.; Wang, Z.; Chen, C. Assessing the Three-North Shelter Forest Program in China by a novel framework for characterizing vegetation changes. ISPRS J. Photogramm. Remote Sens. 2017, 133, 75–88. [Google Scholar] [CrossRef]
  22. Wang, C.; Chu, X.; Zhan, J.; Wang, P.; Zhang, F.; Xin, Z. Factors Contributing to Efficient Forest Production in the Region of the Three-North Shelter Forest Program, China. Sustainability 2020, 12, 302. [Google Scholar] [CrossRef] [Green Version]
  23. Chunya, Z.; Zhongqi, X.; Changming, M.; Shoujia, S.; Tengfei, Y. The Factors Influencing the Poplar Shelterbelt Degradation in the Bashang Plateau of Northwest Hebei Province. For. Resour. Manag. 2018, 1, 9–16. (In Chinese) [Google Scholar]
  24. Huan, Z.; Jun, C.; Hua-bing, W.; Bo, S.; Guo-dong, J.; Zi-qiang, L.; Xin-xiao, Y.; Jia, Z. Water utilization characteristics of the degraded poplar shelterbelts in Zhangbei, Hebei, China. Chin. J. Appl. Ecol. 2018, 29, 1381–1388. (In Chinese) [Google Scholar]
  25. Jump, A.S.; Ruiz-Benito, P.; Greenwood, S.; Allen, C.D.; Kitzberger, T.; Fensham, R.; Martinez-Vilalta, J.; Lloret, F. Structural overshoot of tree growth with climate variability and the global spectrum of drought-induced forest dieback. Glob. Chang. Biol. 2017, 23, 3742–3757. [Google Scholar] [CrossRef]
  26. Sanchez-Salguero, R.; Julio Camarero, J. Greater sensitivity to hotter droughts underlies juniper dieback and mortality in Mediterranean shrublands. Sci. Total Environ. 2020, 721, 137599. [Google Scholar] [CrossRef] [PubMed]
  27. Ji, Y.; Zhou, G.; Li, Z.; Wang, S.; Zhou, H.; Song, X. Triggers of widespread dieback and mortality of poplar (Populus spp.) plantations across northern China. J. Arid Environ. 2020, 174, 104076. [Google Scholar] [CrossRef]
  28. Shimizu, H. An overview of the “Three-North” Shelterbelt project in China. J. For. Stud. China 2012, 14, 70–79. [Google Scholar] [CrossRef]
  29. Zhang, F.; Xing, Z.; Rees, H.W.; Dong, Y.; Li, S.; Meng, F. Assessment of effects of two runoff control engineering practices on soil water and plant growth for afforestation in a semi-arid area after 10 years. Ecol. Eng. 2014, 64, 430–442. [Google Scholar] [CrossRef]
  30. Zhang, Y.; Peng, C.; Li, W.; Tian, L.; Zhu, Q.; Chen, H.; Fang, X.; Zhang, G.; Liu, G.; Mu, X.; et al. Multiple afforestation programs accelerate the greenness in the ‘Three North’ region of China from 1982 to 2013. Ecol. Indic. 2016, 61, 404–412. [Google Scholar] [CrossRef]
  31. Liu, J.; Dang, P.; Gao, Y.; Zhu, H.; Zhu, H.; Zhao, F.; Zhao, Z. Effects of tree species and soil properties on the composition and diversity of the soil bacterial community following afforestation. For. Ecol. Manag. 2018, 427, 342–349. [Google Scholar] [CrossRef]
  32. Kaldorf, M.; Fladung, M.; Muhs, H.J.; Buscot, F. Mycorrhizal colonization of transgenic aspen in a field trial. Planta 2002, 214, 653–660. [Google Scholar] [CrossRef] [PubMed]
  33. Stefani, F.O.P.; Moncalvo, J.-M.; Seguin, A.; Berube, J.A.; Hamelin, R.C. Impact of an 8-Year-Old Transgenic Poplar Plantation on the Ectomycorrhizal Fungal Community. Appl. Environ. Microbiol. 2009, 75, 7527–7536. [Google Scholar] [CrossRef] [Green Version]
  34. Danielsen, L.; Thuermer, A.; Meinicke, P.; Buee, M.; Morin, E.; Martin, F.; Pilate, G.; Daniel, R.; Polle, A.; Reich, M. Fungal soil communities in a young transgenic poplar plantation form a rich reservoir for fungal root communities. Ecol. Evol. 2012, 2, 1935–1948. [Google Scholar] [CrossRef] [PubMed]
  35. Hao, D.-C.; Song, S.-M.; Mu, J.; Hu, W.-L.; Xiao, P.-G. Unearthing microbial diversity of Taxus rhizosphere via MiSeq high-throughput amplicon sequencing and isolate characterization. Sci. Rep. 2016, 6, 22006. [Google Scholar] [CrossRef] [Green Version]
  36. Dai, Z.; Su, W.; Chen, H.; Barberan, A.; Zhao, H.; Yu, M.; Yu, L.; Brookes, P.C.; Schadt, C.W.; Chang, S.X.; et al. 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, 3452–3461. [Google Scholar] [CrossRef] [PubMed]
  37. Menendez, E.; Perez-Yepez, J.; Hernandez, M.; Rodriguez-Perez, A.; Velazquez, E.; Leon-Barrios, M. Plant Growth Promotion Abilities of Phylogenetically Diverse Mesorhizobium Strains: Effect in the Root Colonization and Development of Tomato Seedlings. Microorganisms 2020, 8, 412. [Google Scholar] [CrossRef] [Green Version]
  38. Luo, Y.; Wang, F.; Huang, Y.; Zhou, M.; Gao, J.; Yan, T.; Sheng, H.; An, L. Sphingomonas sp. Cra20 Increases Plant Growth Rate and Alters Rhizosphere Microbial Community Structure of Arabidopsis thaliana Under Drought Stress. Front. Microbiol. 2019, 10, 1221. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Wang, Q.; Ge, C.; Xu, S.A.; Wu, Y.; Sahito, Z.A.; Ma, L.; Pan, F.; Zhou, Q.; Huang, L.; Feng, Y.; et al. The endophytic bacterium Sphingomonas SaMR12 alleviates Cd stress in oilseed rape through regulation of the GSH-AsA cycle and antioxidative enzymes. BMC Plant Biol. 2020, 20, 63. [Google Scholar] [CrossRef]
  40. Talia, P.; Sede, S.M.; Campos, E.; Rorig, M.; Principi, D.; Tosto, D.; Hopp, H.E.; Grasso, D.; Cataldi, A. Biodiversity characterization of cellulolytic bacteria present on native Chaco soil by comparison of ribosomal RNA genes. Res. Microbiol. 2012, 163, 221–232. [Google Scholar] [CrossRef] [PubMed]
  41. Natsagdorj, O.; Sakamoto, H.; Santiago, D.M.O.; Santiago, C.D.; Orikasa, Y.; Okazaki, K.; Ikeda, S.; Ohwada, T. Variovorax sp. Has an Optimum Cell Density to Fully Function as a Plant Growth Promoter. Microorganisms 2019, 7, 82. [Google Scholar] [CrossRef] [Green Version]
  42. Sun, S.-L.; Yang, W.-L.; Fang, W.-W.; Zhao, Y.-X.; Guo, L.; Dai, Y.-J. The Plant Growth-Promoting Rhizobacterium Variovorax boronicumulans CGMCC 4969 Regulates the Level of Indole-3-Acetic Acid Synthesized from Indole-3-Acetonitrile. Appl. Environ. Microbiol. 2018, 84, e00298-18. [Google Scholar] [CrossRef] [Green Version]
  43. Buenger, W.; Jiang, X.; Mueller, J.; Hurek, T.; Reinhold-Hurek, B. Novel cultivated endophytic Verrucomicrobia reveal deep-rooting traits of bacteria to associate with plants. Sci. Rep. 2020, 10, 8692. [Google Scholar] [CrossRef]
  44. Bergmann, G.T.; Bates, S.T.; Eilers, K.G.; Lauber, C.L.; Caporaso, J.G.; Walters, W.A.; Knight, R.; Fierer, N. The under-recognized dominance of Verrucomicrobia in soil bacterial communities. Soil Biol. Biochem. 2011, 43, 1450–1455. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Lebeis, S.L.; Paredes, S.H.; Lundberg, D.S.; Breakfield, N.; Gehring, J.; McDonald, M.; Malfatti, S.; del Rio, T.G.; Jones, C.D.; Tringe, S.G.; et al. Salicylic acid modulates colonization of the root microbiome by specific bacterial taxa. Science 2015, 349, 860–864. [Google Scholar] [CrossRef] [Green Version]
  46. Utturkar, S.M.; Cude, W.N.; Robeson, M.S., Jr.; Yang, Z.K.; Klingeman, D.M.; Land, M.L.; Allman, S.L.; Lu, T.-Y.S.; Brown, S.D.; Schadt, C.W.; et al. Enrichment of Root Endophytic Bacteria from Populus deltoides and Single-Cell-Genomics Analysis. Appl. Environ. Microbiol. 2016, 82, 5698–5708. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Peiffer, J.A.; Spor, A.; Koren, O.; Jin, Z.; Tringe, S.G.; Dangl, J.L.; Buckler, E.S.; Ley, R.E. Diversity and heritability of the maize rhizosphere microbiome under field conditions. Proc. Natl. Acad. Sci. USA 2013, 110, 6548–6553. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Aguirre-von-Wobeser, E.; Rocha-Estrada, J.; Shapiro, L.R.; de la Torre, M. Enrichment of Verrucomicrobia, Actinobacteria and Burkholderiales drives selection of bacterial community from soil by maize roots in a traditional milpa agroecosystem. PLoS ONE 2018, 13, e0208852. [Google Scholar] [CrossRef] [PubMed]
  49. Sessitsch, A.; Hardoim, P.; Doering, J.; Weilharter, A.; Krause, A.; Woyke, T.; Mitter, B.; Hauberg-Lotte, L.; Friedrich, F.; Rahalkar, M.; et al. Functional Characteristics of an Endophyte Community Colonizing Rice Roots as Revealed by Metagenomic Analysis. Mol. Plant Microbe Interact. 2012, 25, 28–36. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Edwards, J.; Johnson, C.; Santos-Medellin, C.; Lurie, E.; Podishetty, N.K.; Bhatnagar, S.; Eisen, J.A.; Sundaresan, V. Structure, variation, and assembly of the root-associated microbiomes of rice. Proc. Natl. Acad. Sci. USA 2015, 112, E911–E920. [Google Scholar] [CrossRef] [Green Version]
  51. Krings, M.; Hass, H.; Kerp, H.; Taylor, T.N.; Agerer, R.; Dotzler, N. Endophytic cyanobacteria in a 400-million-yr-old land plant: A scenario for the origin of a symbiosis? Rev. Palaeobot. Palynol. 2009, 153, 62–69. [Google Scholar] [CrossRef]
  52. He, S.; Stevens, S.L.R.; Chan, L.-K.; Bertilsson, S.; del Rio, T.G.; Tringe, S.G.; Malmstrom, R.R.; McMahon, K.D. Ecophysiology of Freshwater Verrucomicrobia Inferred from Metagenome-Assembled Genomes. mSphere 2017, 2, e00277-17. [Google Scholar] [CrossRef] [Green Version]
  53. Pold, G.; Conlon, E.M.; Huntemann, M.; Pillay, M.; Mikhailova, N.; Stamatis, D.; Reddy, T.B.K.; Daum, C.; Shapiro, N.; Kyrpides, N.; et al. Genome Sequence of Verrucomicrobium sp. Strain GAS474, a Novel Bacterium Isolated from Soil. Genome Announc. 2018, 6, e01451-17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Noraini, C.H.C.; Morad, N.; Norli, I.; Teng, T.T.; Ogugbue, C.J. Methylene Blue Degradation by Sphingomonas paucimobilis under Aerobic Conditions. Water Air Soil Pollut. 2012, 223, 5131–5142. [Google Scholar] [CrossRef]
  55. Ravintheran, S.K.; Sivaprakasam, S.; Loke, S.; Lee, S.Y.; Manickam, R.; Yahya, A.; Croft, L.; Millard, A.; Parimannan, S.; Rajandas, H. Complete genome sequence of Sphingomonas paucimobilis AIMST S2, a xenobiotic-degrading bacterium. Sci. Data 2019, 6, 280. [Google Scholar] [CrossRef]
  56. Worsley, S.F.; Newitt, J.; Rassbach, J.; Batey, S.F.D.; Holmes, N.A.; Murrell, J.C.; Wilkinson, B.; Hutchings, M.I. Streptomyces Endophytes Promote Host Health and Enhance Growth across Plant Species. Appl. Environ. Microbiol. 2020, 86, e01053-20. [Google Scholar] [CrossRef]
  57. Ferreira, A.C.; Nobre, M.F.; Moore, E.; Rainey, F.A.; Battista, J.R.; da Costa, M.S. Characterization and radiation resistance of new isolates of Rubrobacter radiotolerans and Rubrobacter xylanophilus. Extrem. Life Extrem. Cond. 1999, 3, 235–238. [Google Scholar] [CrossRef] [PubMed]
  58. Berendsen, R.L.; Pieterse, C.M.; Bakker, P.A. The rhizosphere microbiome and plant health. Trends Plant Sci. 2012, 17, 478–486. [Google Scholar] [CrossRef]
  59. Bulgarelli, D.; Rott, M.; Schlaeppi, K.; van Themaat, E.V.L.; Ahmadinejad, N.; Assenza, F.; Rauf, P.; Huettel, B.; Reinhardt, R.; Schmelzer, E.; et al. Revealing structure and assembly cues for Arabidopsis root-inhabiting bacterial microbiota. Nature 2012, 488, 91–95. [Google Scholar] [CrossRef]
  60. Lundberg, D.S.; Lebeis, S.L.; Paredes, S.H.; Yourstone, S.; Gehring, J.; Malfatti, S.; Tremblay, J.; Engelbrektson, A.; Kunin, V.; del Rio, T.G.; et al. Defining the core Arabidopsis thaliana root microbiome. Nature 2012, 488, 86–94. [Google Scholar] [CrossRef] [Green Version]
  61. White, P.J.; George, T.S.; Dupuy, L.X.; Karley, A.J.; Valentine, T.A.; Wiesel, L.; Wishart, J. Root traits for infertile soils. Front. Plant Sci. 2013, 4, 193. [Google Scholar] [CrossRef] [Green Version]
  62. Armengaud, P.; Breitling, R.; Amtmann, A. Coronatine-Insensitive 1 (COI1) Mediates Transcriptional Responses of Arabidopsis thaliana to External Potassium Supply. Mol. Plant 2010, 3, 309–405. [Google Scholar] [CrossRef]
  63. Troufflard, S.; Mullen, W.; Larson, T.R.; Graham, I.A.; Crozier, A.; Amtmann, A.; Armengaud, P. Potassium deficiency induces the biosynthesis of oxylipins and glucosinolates in Arabidopsis thaliana. BMC Plant Biol. 2010, 10, 172. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Meena, V.S.; Maurya, B.R.; Verma, J.P. Does a rhizospheric microorganism enhance K+ availability in agricultural soils? Microbiol. Res. 2014, 169, 337–347. [Google Scholar] [CrossRef] [PubMed]
  65. Dong, X.; Lv, L.; Wang, W.; Liu, Y.; Yin, C.; Xu, Q.; Yan, H.; Fu, J.; Liu, X. Differences in Distribution of Potassium-Solubilizing Bacteria in Forest and Plantation Soils in Myanmar. Int. J. Environ. Res. Public Health 2019, 16, 700. [Google Scholar] [CrossRef] [Green Version]
  66. Mark Ibekwe, A.; Ors, S.; Ferreira, J.F.S.; Liu, X.; Suarez, D.L. Seasonal induced changes in spinach rhizosphere microbial community structure with varying salinity and drought. Sci. Total. Environ. 2017, 579, 1485–1495. [Google Scholar] [CrossRef] [Green Version]
  67. Sun, S.; Gao, H.; Liu, Y.; Jin, L.; Wang, R.; Wang, X.; Wang, Q.; Yin, Y.; Zhang, Y.; Wang, H. Co-existence of a novel plasmid-mediated efflux pump with colistin resistance gene mcr in one plasmid confers transferable multidrug resistance in Klebsiella pneumoniae. Emerg. Microbes Infect. 2020, 9, 1102–1113. [Google Scholar] [CrossRef] [PubMed]
  68. Wu, L.; Wang, J.; Huang, W.; Wu, H.; Chen, J.; Yang, Y.; Zhang, Z.; Lin, W. Plant-microbe rhizosphere interactions mediated by Rehmannia glutinosa root exudates under consecutive monoculture. Sci. Rep. 2015, 5, 15871. [Google Scholar] [CrossRef] [Green Version]
  69. Xia, Z.; Bai, E.; Wang, Q.; Gao, D.; Zhou, J.; Jiang, P.; Wu, J. Biogeographic Distribution Patterns of Bacteria in Typical Chinese Forest Soils. Front. Microbiol. 2016, 7, 1106. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Delgado Restrepo, O.M.; Menjivar Flores, J.C.; Muñoz Arboleda, F. Influence of management systems on the nitrogen mineralization and fertilization of sugarcane. Rev. Fac. Nac. Agron. Medellín 2016, 69, 7755–7762. [Google Scholar] [CrossRef]
  71. Inclan, R.; De la Torre, D.; Benito, M.; Rubio, A. Soil CO2 efflux in a mixed pine-oak forest in Valsain (central Spain). Sci. World J. 2007, 7 (Suppl. 1), 166–174. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Bao, S.D. Soil and Agricultural Chemistry Analysis; Agricultural Press: Beijing, China, 2000; pp. 355–356. [Google Scholar]
  73. Wang, J.-H.; Lu, J.; Zhang, Y.-X.; Wu, J.; Luo, Y.; Liu, H. Metagenomic analysis of antibiotic resistance genes in coastal industrial mariculture systems. Bioresour. Technol. 2018, 253, 235–243. [Google Scholar] [CrossRef]
  74. Li, D.; Liu, C.-M.; Luo, R.; Sadakane, K.; Lam, T.-W. MEGAHIT: An ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics 2015, 31, 1674–1676. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Nielsen, H.B.; Almeida, M.; Juncker, A.S.; Rasmussen, S.; Li, J.; Sunagawa, S.; Plichta, D.R.; Gautier, L.; Pedersen, A.G.; Le Chatelier, E.; et al. Identification and assembly of genomes and genetic elements in complex metagenomic samples without using reference genomes. Nat. Biotechnol. 2014, 32, 822–828. [Google Scholar] [CrossRef]
  76. Oh, J.; Byrd, A.L.; Deming, C.; Conlan, S.; Kong, H.H.; Segre, J.A.; Progra, N.C.S. Biogeography and individuality shape function in the human skin metagenome. Nature 2014, 514, 59–64. [Google Scholar] [CrossRef] [Green Version]
  77. Zhu, W.; Lomsadze, A.; Borodovsky, M. Ab initio gene identification in metagenomic sequences. Nucleic Acids Res. 2010, 38, e132. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Buchfink, B.; Xie, C.; Huson, D.H. Fast and sensitive protein alignment using DIAMOND. Nat. Methods 2015, 12, 59–60. [Google Scholar] [CrossRef]
  79. Shi, Y.; Fan, K.K.; Li, Y.T.; Yang, T.; He, J.S.; Chu, H.Y. Archaea Enhance the Robustness of Microbial Co-occurrence Networks in Tibetan Plateau Soils. Soil Sci. Soc. Am. J. 2019, 83, 1093–1099. [Google Scholar] [CrossRef]
  80. Jia, B.; Raphenya, A.R.; Alcock, B.; Waglechner, N.; Guo, P.; Tsang, K.K.; Lago, B.A.; Dave, B.M.; Pereira, S.; Sharma, A.N.; et al. CARD 2017: Expansion and model-centric curation of the comprehensive antibiotic resistance database. Nucleic Acids Res. 2017, 45, D566–D573. [Google Scholar] [CrossRef]
  81. White, J.R.; Nagarajan, N.; Pop, M. Statistical methods for detecting differentially abundant features in clinical metagenomic samples. PLoS Comput. Biol. 2009, 5, e1000352. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Poplar tree growth conditions at the sampling site and comparison of physical and chemical factors between groups. (a) Actual diagram of different growth states of poplar. Poplars in undegraded condition (C); Medium-degraded poplars from the treetop (M); Completely degraded poplars (H); (b) variation in soil properties between different groups, with different parallel samples in the group set as the black dots. ** p < 0.01, * p < 0.05 (After Shapiro-Wilk and Levene’s test, Anova with LSD method or Kruskal–Wallis tests and multiple comparison with Bonferroni method were carried out.).
Figure 1. Poplar tree growth conditions at the sampling site and comparison of physical and chemical factors between groups. (a) Actual diagram of different growth states of poplar. Poplars in undegraded condition (C); Medium-degraded poplars from the treetop (M); Completely degraded poplars (H); (b) variation in soil properties between different groups, with different parallel samples in the group set as the black dots. ** p < 0.01, * p < 0.05 (After Shapiro-Wilk and Levene’s test, Anova with LSD method or Kruskal–Wallis tests and multiple comparison with Bonferroni method were carried out.).
Ijms 22 03438 g001
Figure 2. Comparative analysis of bacterial communities at different taxonomic levels. (ac) Unweighted pair-group method with arithmetic mean (UPGMA) clustering based on Bray–Curtis distances and the relative abundance of the most abundant bacteria at the phylum, family, and genus level; (df) Indicator bacteria with significant differences between C, M, and H groups. “*” represents adjust q value < 0.05, “**” represents adjust q value < 0.01.
Figure 2. Comparative analysis of bacterial communities at different taxonomic levels. (ac) Unweighted pair-group method with arithmetic mean (UPGMA) clustering based on Bray–Curtis distances and the relative abundance of the most abundant bacteria at the phylum, family, and genus level; (df) Indicator bacteria with significant differences between C, M, and H groups. “*” represents adjust q value < 0.05, “**” represents adjust q value < 0.01.
Ijms 22 03438 g002
Figure 3. Redundancy analysis for the microbes among samples at (a) phylum, (b) family, and (c) genus levels with edaphic physicochemical factors, respectively.
Figure 3. Redundancy analysis for the microbes among samples at (a) phylum, (b) family, and (c) genus levels with edaphic physicochemical factors, respectively.
Ijms 22 03438 g003
Figure 4. Community networks of C, M and H set up co-occurrence based on top 100 genera respectively. Each node represents a genus, and the size and color represent its abundance and the annotation color of the phylum which it belongs, respectively. Green line—positive correlation, blue line—negative correlation.
Figure 4. Community networks of C, M and H set up co-occurrence based on top 100 genera respectively. Each node represents a genus, and the size and color represent its abundance and the annotation color of the phylum which it belongs, respectively. Green line—positive correlation, blue line—negative correlation.
Ijms 22 03438 g004
Figure 5. (a) Heat map analysis of the top 35 most abundant Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways in at least one soil sample using normalized abundance. Heat map is color-coded based on row z-scores. (b) Spearman correlation analysis between genus with significant difference between C and H and metabolic pathway on level 3. “*” represents q value < 0.05; “**” represents q value < 0.01.
Figure 5. (a) Heat map analysis of the top 35 most abundant Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways in at least one soil sample using normalized abundance. Heat map is color-coded based on row z-scores. (b) Spearman correlation analysis between genus with significant difference between C and H and metabolic pathway on level 3. “*” represents q value < 0.05; “**” represents q value < 0.01.
Ijms 22 03438 g005
Figure 6. (a) The bar width of the bars between ARG types and samples correlates to the percentages of respective ARGs in these samples (in ppm). The different colors in the circle represent different samples and ARGs. (b) Heat map analysis of the top 30 most abundant ARGs in at least one group sample using normalized abundance and color-coded based on row z-scores. (c) Histogram of four resistance mechanisms of ARGs between groups. (d) Histogram of ARGs species attribution analysis. Different colors represent corresponding genera under three phyla. (e) Abundance statistics of the top 30 abundant ARGs-related antibiotics.
Figure 6. (a) The bar width of the bars between ARG types and samples correlates to the percentages of respective ARGs in these samples (in ppm). The different colors in the circle represent different samples and ARGs. (b) Heat map analysis of the top 30 most abundant ARGs in at least one group sample using normalized abundance and color-coded based on row z-scores. (c) Histogram of four resistance mechanisms of ARGs between groups. (d) Histogram of ARGs species attribution analysis. Different colors represent corresponding genera under three phyla. (e) Abundance statistics of the top 30 abundant ARGs-related antibiotics.
Ijms 22 03438 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, J.; He, X.; Sun, J.; Ma, Y. A Degeneration Gradient of Poplar Trees Contributes to the Taxonomic, Functional, and Resistome Diversity of Bacterial Communities in Rhizosphere Soils. Int. J. Mol. Sci. 2021, 22, 3438. https://doi.org/10.3390/ijms22073438

AMA Style

Liu J, He X, Sun J, Ma Y. A Degeneration Gradient of Poplar Trees Contributes to the Taxonomic, Functional, and Resistome Diversity of Bacterial Communities in Rhizosphere Soils. International Journal of Molecular Sciences. 2021; 22(7):3438. https://doi.org/10.3390/ijms22073438

Chicago/Turabian Style

Liu, Juan, Xiangwei He, Jingya Sun, and Yuchao Ma. 2021. "A Degeneration Gradient of Poplar Trees Contributes to the Taxonomic, Functional, and Resistome Diversity of Bacterial Communities in Rhizosphere Soils" International Journal of Molecular Sciences 22, no. 7: 3438. https://doi.org/10.3390/ijms22073438

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop