Taxonomic and functional analyses reveal existence of virulence and antibiotic resistance genes in beach sand bacterial populations

Coastal sands are important natural recreational facilities that have become hotspots for tourism and economic development. However, these sands harbour diverse microbial assemblages that play a critical role in the balance between public health and ecology. In this study, targeted high-throughput sequencing analysis was used to identify sand-borne bacterial populations at four public beaches in Durban. The effect of heavy metal in shaping the distribution of bacterial metacommunities was determined using canonical correspondence analysis (CCA), while the functional gene profiles were predicted using PICRUSt2 analysis. Sequences matching those of the bacterial phylum Proteobacteria were the most abundant in all samples, followed by those of the phyla Firmicutes, Actinobacteria, Bacteroidetes, and Gemmatimonadetes. Genus-level taxonomic analysis showed the presence of 1163 bacterial genera in all samples combined. The distribution of bacterial communities was shaped by heavy metal concentrations, with the distribution of Flavobacteria, Bacteroidia, and Deltaproteobacteria influenced by Pb and Zn, while B and Cr influenced the distribution of Clostridia and Gammaproteobacteria, respectively. Identified antibiotic resistance genes included the peptidoglycan biosynthesis gene II, III, IV, and V, as well as the polymyxin resistance gene, while the virulence genes included the sitA, fimB, aerobactin synthase, and pilL gene. Our findings demonstrate that beach sand-borne bacteria are reservoirs of virulence and antibiotic resistance genes. Contamination of beach sands with heavy metals selects for both heavy metal resistance and antibiotic resistance in beach sand bacterial communities. Children and immunocompromised people engaging in recreational activities on beaches may be exposed to higher risk of infection.