Integrated metagenomics and molecular ecological network analysis of bacterial community composition during the phytoremediation of cadmium-contaminated soils by bioenergy crops

Two energy crops (maize and soybean) were used in the remediation of cadmium-contaminated soils. These crops were used because they are fast growing, have a large biomass and are good sources for bioenergy production. The total accumulation of cadmium in maize and soybean plants was 393.01 and 263.24 μg pot−1, respectively. The rhizosphere bacterial community composition was studied by MiSeq sequencing. Phylogenetic analysis was performed using 16S rRNA gene sequences. The rhizosphere bacteria were divided into 33 major phylogenetic groups according to phyla. The dominant phylogenetic groups included Proteobacteria, Acidobacteria, Actinobacteria, Gemmatimonadetes, and Bacteroidetes. Based on principal component analysis (PCA) and unweighted pair group with arithmetic mean (UPGMA) analysis, we found that the bacterial community was influenced by cadmium addition and bioenergy cropping. Three molecular ecological networks were constructed for the unplanted, soybean- and maize-planted bacterial communities grown in 50 mg kg−1 cadmium-contaminated soils. The results indicated that bioenergy cropping increased the complexity of the bacterial community network as evidenced by a higher total number of nodes, the average geodesic distance (GD), the modularity and a shorter geodesic distance. Proteobacteria and Acidobacteria were the keystone bacteria connecting different co-expressed operational taxonomic units (OTUs). The results showed that bioenergy cropping altered the topological roles of individual OTUs and keystone populations. This is the first study to reveal the effects of bioenergy cropping on microbial interactions in the phytoremediation of cadmium-contaminated soils by network reconstruction. This method can greatly enhance our understanding of the mechanisms of plant-microbe-metal interactions in metal-polluted ecosystems.