Polymetallic contamination drives indigenous microbial community assembly dominated by stochastic processes at Pb-Zn smelting sites.

Indigenous microbial communities in smelting areas are crucial for maintaining fragile ecosystem functions. However, the community assembly process and their responses to polymetallic pollution are poorly understood, especially the taxa in each bin from the amplicons that contributed to the assembly process. Herein, microbial diversity, co-occurrence patterns, assembly process and the intrinsic mechanisms across contamination gradients at a typical PbZn smelting site were systematically unravelled by high-throughput sequencing. The results showed a consistent compositional profile among the indigenous communities across sampling sites, wherein genera KD4-96 from Chloroflexi and Sphingomonas from Proteobacteria emerged as the most abundant taxa. Network modularity of the high- and middle-contaminated communities at Pb and Zn smelting sites was >0.44, indicating that community populations were clustered into modules to resist high heavy metal stress. Stochastic processes dominated the community assembly, with the greatest contribution from drift (DR), which was significantly correlated with Pb, Zn, Cr and Cu contents. What's particular was that the DR-controlled bins were dominated by Proteobacteria (typical r-strategists), while the HoS-controlled bins were by Chloroflexi (typical K-strategists). Furthermore, the proportion of DR in the bins dominated by Sphingomonadaceae (phylum Proteobacteria) increased gradually with the increase of heavy metal contents. These discoveries provide essential insights for community control in restoring and mitigating soil degradation at PbZn smelting sites.