Structural and functional responses of bacterial and fungal communities to multiple heavy metal exposure in arid loess

Abstract

Concentration gradients of multiple heavy metals (HMs) in the arid loess region near a smelter were determined. In order to understand the response of soil microbes to multiple HM gradients, bacterial and fungal community structures and functions were analyzed using high-throughput RNA gene sequencing and the PICRUSt method. RDA/PCA analyses revealed that soil pH, HMs, and electrical conductivity (EC) jointly affected the bacterial communities in the soils. The soil microbial community structures responded differently to HMs, EC, and pH. High HMs increased the abundances of the bacterial phyla Actinobacteria, Bacteroidetes, Deinococcus-Thermus, and Chloroflexi, and the genera Blastococcus, Rubrobacter, Quadrisphaera, and Tunicatimonas, whereas they decreased the abundances of the phyla Proteobacteria and Acidobacteria and the genera Streptomyces and Nocardioides. High EC and low pH decreased the abundance of most of the dominant bacterial phyla but increased the abundances of Firmicutes, Deinococcus-Thermus, and Nitrospirae. Furthermore, high HMs and EC reduced the numbers of soil-specific bacterial and fungal groups and drove the succession of certain groups that were highly resistant to increased HMs and EC. In addition, many bacterial and fungal groups exhibited different response patterns to each HM, implying that, in multiple HM-contaminated soils, HMs jointly shaped the microbial communities. PICRUSt analysis suggested that high HMs significantly decreased the total gene abundance and most KEGG modules in the soils. High EC and low pH significantly enhanced the abundances of several two-component system-, electron transfer-, and methanogenesis-related modules. We conclude that excessive multiple HMs and EC principally repressed the microbial activity and severely drove the gradient succession of bacterial and fungal communities in the arid loess region.