Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Microbiol., 30 July 2025

Sec. Microbiological Chemistry and Geomicrobiology

Volume 16 - 2025 | https://doi.org/10.3389/fmicb.2025.1607082

This article is part of the Research TopicMicrobial Diversity, Function and Process Involved in Redox Reactions of Heavy Metals/Metalloids and Related PollutantsView all 3 articles

Distinct bacterial community structures and arsenic biotransformation gene profiles in dust

Yi YinYi Yin1Yu-Ting LinYu-Ting Lin1Gong-Ren HuGong-Ren Hu1Rui-Lian YuRui-Lian Yu1Xiao-Hui SunXiao-Hui Sun2Yu Yan
Yu Yan1*
  • 1Department of Environmental Science and Engineering, College of Chemical Engineering, Huaqiao University, Xiamen, China
  • 2Department of Bioengineering and Biotechnology, College of Chemical Engineering, Huaqiao University, Xiamen, China

Introduction: Microorganisms, which are ubiquitous in the environment, have evolved a diverse array of arsenic biotransformation genes (ABGs). Dust harbors a wide range of microorganisms. However, the distinct characteristics of bacterial community structures and ABG profiles in dust, compared with those in other environments such as soil and water, remain poorly understood.

Methods: In this study, dust samples were simultaneously collected alongside surrounding soil and seawater samples in Xiamen, a coastal city of China, to investigate the distinct profiles and potential sources of bacterial communities and ABGs in dust using 16S rRNA gene amplicon sequencing and metagenomic sequencing.

Results and discussion: Abundant and diverse bacterial communities and ABGs were detected in dust, revealing significant differences in community structures and ABG profiles compared with those in soil and seawater. Soil was identified as the primary source for both bacterial communities and ABGs in dust through fast expectation–maximization microbial source tracking (FEAST). Acetobacteraceae, which showed significantly greater relative abundance (p < 0.001) in dust than in soil and seawater, was also identified as a keystone taxon in the dust bacterial co-occurrence network. Furthermore, metagenome-assembled genomes (MAGs) affiliated with Acetobacteraceae were effectively recovered from dust via metagenomic binning, and these MAGs harbored an array of ABGs, indicating that Acetobacteraceae could be important hosts for ABGs in dust. Overall, our findings offer new insights into bacterial communities and ABGs in dust, thereby improving our understanding of arsenic biogeochemical cycling.

1 Introduction

Arsenic, a widely distributed toxic metalloid in various environments, originates from both natural processes and human activities (Bowell et al., 2014; Huang et al., 2024). Arsenic contamination poses a considerable risk to both ecosystems and human health, thereby posing a global environmental concern (Chen and Costa, 2021; Tuli et al., 2010). The main concern regarding arsenic contamination has traditionally been associated with its presence in groundwater, surface water, soil, sediment, and biota (Bhattacharya et al., 2002; Raju, 2022). However, arsenic in the air, which exists primarily (> 90%) in particulate form, could also potentially impact the health of the population (Sánchez-Rodas et al., 2015). Dust, a prominent component of atmospheric particulate matter, can facilitate the transport of arsenic (Yin et al., 2013). The total concentration of arsenic in dust and the associated risks have been extensively studied (Balakrishna et al., 2011; Ma et al., 2021). For example, the arsenic concentrations in dust samples collected from kindergartens in Xiamen, China, were found to be approximately 2.2 times greater than the background value of the soil (Ma et al., 2021). Furthermore, the toxicity of arsenic is closely linked to its chemical speciation (Jain and Ali, 2000). The speciation of arsenic in atmospheric particulates, including dust, has also been analyzed recently (Sánchez-Rodas et al., 2015; Xie et al., 2021; Yan et al., 2024; Yin et al., 2017). Generally, the predominant species of arsenic found in dust is inorganic arsenate [As(V)]. In addition, other species, such as inorganic arsenite [As(III)], trimethylarsine oxide [TMAs(V)O], and dimethylarsenate [DMAs(V)], have also been identified in dust.

Microorganisms are ubiquitous in the environment. In addition to soil and seawater, air also serves as a crucial reservoir of microorganisms. Numerous studies have demonstrated that microorganisms are key components of dust (Ding et al., 2020; Leppänen et al., 2018). Moreover, microorganisms can drive the biogeochemical cycling of arsenic, and arsenic biotransformation genes (ABGs) are widespread in microorganisms owing to their longstanding coexistence with arsenic (Ben Fekih et al., 2018; Mukhopadhyay et al., 2002; Oremland and Stolz, 2003). For example, the redox transformation of arsenic can be achieved through As(V) reduction mediated by genes encoding cytoplasmic As(V) reductases (ArsC) or respiratory As(V) reductases (ArrAB) and the oxidation of As(III) mediated by genes encoding As(III) oxidases (AioAB/ArxAB) (Andres and Bertin, 2016; Oremland and Stolz, 2005). Arsenic methylation and demethylation can be facilitated by genes encoding As(III) S-adenosylmethionine methyltransferases (ArsM) and C-As lyases (ArsI), respectively (Qin et al., 2006; Yoshinaga and Rosen, 2014). The distribution and abundance of ABGs have been extensively studied in water, soil, and sediment through the utilization of quantitative PCR (qPCR) (Guo et al., 2019; Qiao et al., 2023; Zecchin et al., 2021; Zhang et al., 2015). In particular, our recent study revealed the presence of abundant ABGs (e.g., the arsM, arsC, and aioA genes) in dust for the first time by using a high-throughput qPCR chip designed for detecting ABGs (Yan et al., 2024; Zhao et al., 2019). Although qPCR method is extensively used for arsenic-related gene detection, this method could only target a limited number of specific ABGs and was prone to non-specific amplification introduced by the primers due to primer design constraints (Song et al., 2022). While qPCR can determine the abundance of ABGs, it cannot provide information on their taxonomic compositions. Metagenomic sequencing technology enables the simultaneous probing of both the abundance and taxonomic profiling of numerous ABGs, thereby characterizing microbial-driven arsenic biotransformation at gene- and species-level resolutions (Li L. et al., 2023). Therefore, metagenomic sequencing can be further employed to elucidate ABGs in the environment more comprehensively. Furthermore, airborne and dust-associated microbial communities can originate from other environments, such as soil and water (Archer et al., 2023; Barberan et al., 2015; Bowers et al., 2011). However, the differences in bacterial community structures and ABG profiles in dust compared with those in other environments, including soil and water, as well as the extent of contribution from these environments to the sources of bacteria and ABGs in dust, still require further investigation.

In this study, we simultaneously collected dust samples and surrounding soil and seawater samples in Xiamen, China, to systematically characterize the distinct profiles and potential sources of bacterial communities and ABGs in dust. This study aimed to (i) characterize the variations in bacterial community structures and ABG profiles across dust, soil, and seawater; (ii) track the sources of bacterial communities and ABGs in dust; and (iii) identify key bacterial taxa and ABGs in dust.

2 Materials and methods

2.1 Sample collection and pretreatment

Xiamen (24°23′N to 24°54′N and 117°52′E to 118°26′E), situated on the southeastern coast of China, is a typical coastal city where the interaction between land and sea creates an environment in which dust, soil, and seawater are closely interconnected. Samples of dust, soil, and seawater were collected from coastal parks situated in Xiamen in September 2023 during a dry period with no rainfall for at least seven consecutive days. The locations of the sampling sites are illustrated in Figure 1. Dust samples were collected by manually sweeping outdoor surfaces with sterilized brushes. Soil samples were collected from a depth of 0–10 cm via a sterile sampling shovel. Seawater samples were collected at depths of 20–30 cm below the sea surface using polyethylene bottles. After thorough mixing, each sample of dust and soil was separated into two parts: one part was air-dried, sieved (<150 μm), and stored at room temperature for determining the total concentrations and speciation of arsenic, while the other portion was frozen at −20°C to facilitate DNA extraction. Each seawater sample was also separated into two parts: one part was acidified with HNO3 for determining the total concentration of arsenic, while the other part was filtered with 0.22 μm mixed cellulose ester (MCE) membranes to collect the microorganisms, after which the membranes were placed in sterile freezing tubes and frozen at −20°C to facilitate DNA extraction.

Figure 1
Map showing Fujian province in southeastern China and Xiamen highlighted with a red circle. The right panel features a detailed map indicating locations of dust, soil, and seawater with symbols: gray circles, brown squares, and blue triangles, respectively. North is indicated. A scale bar denotes distance.

Figure 1. Study area and sampling sites.

2.2 Determination of total arsenic concentrations and arsenic species

The dust and soil samples were digested with a mixture of HNO3 and HCl (3:1, v/v) at 180°C for 24 h (Ma et al., 2021). The digested samples were diluted with 2% (v/v) HNO3, and the acidified seawater samples were then filtered through 0.45 μm filters. The total concentration of arsenic was determined via inductively coupled plasma-mass spectrometry (ICP-MS; 7800, Agilent Technologies, United States) at the Instrumental Analysis Center of Huaqiao University. For quality control, a certified reference material (GSS-5) was used, with recovery rates of arsenic ranging from 97.6% to 111.5%.

The dust and soil samples were extracted with 0.3 M H3PO4 at a solid–liquid ratio of 1:200 (w/v) for 30 min at 150 rpm and 25°C, followed by an additional sonication step for 30 min (Orero Iserte et al., 2004). Both extracts were filtered through 0.45 μm filters. Arsenic species were determined via high-performance liquid chromatography (HPLC, 1260 II, Agilent Technologies, United States) couple with ICP-MS (7800, Agilent Technologies, USA) at the Instrumental Analysis Center of Huaqiao University. A PRP-X100 anion-exchange column (250 mm × 4.1 mm, 10 μm, Hamilton, United States) was used for analysis. The mobile phase consisted of a mixture of 10 mM NH4NO3 and 10 mM NH4H2PO4 adjusted to pH 6.2 with a flow rate of 1 mL·min−1 (Yin et al., 2011; Yan et al., 2024). The arsenic species in the seawater samples were not determined owing to the low total arsenic concentrations in these samples. For quality control, spike recovery experiments were performed with known amounts of arsenic species spiked to the raw samples (Kara et al., 2021; Pizarro et al., 2003), yielding recovery rates between 86.9% and 97.3% for each individual arsenic species.

2.3 DNA extraction and high-throughput 16S rRNA gene sequencing

Total microbial genomic DNA was extracted from the dust and soil samples as well as the filters that captured microorganisms from the seawater samples via a FastDNA Spin Kit for Soil (MP Biomedicals, United States). The extracted DNA was quantified via a Qubit 3.0 fluorometer (Invitrogen, Belgium). The bacterial 16S rRNA gene (V4-V5 region) was amplified via the 515F/907R primers (Ma et al., 2023). Amplicon sequencing was performed on the Illumina NovaSeq-PE250 platform at Shanghai Personal Biotechnology Co., Ltd., China. High-quality sequences were processed via USEARCH (v.11.0.667) (Edgar, 2010) and subsequently clustered into amplicon sequence variants (ASVs) via UNOISE3 (Edgar, 2016). ASV classification was performed with the Silva database (v. 138) (Quast et al., 2013).

2.4 Metagenomic sequencing

Dust, soil, and seawater DNA were used for library preparation and shotgun metagenomic sequencing. The Illumina NovaSeq platform at Shanghai Personal Biotechnology Co., Ltd., China, was utilized for paired-end sequencing. Approximately 10 Gb of raw sequencing data was generated for each sample and subsequently filtered, subjected to quality control, and trimmed via Fastp (v0.21.0) (Chen et al., 2018b). The clean reads were individually assembled into contigs via Megahit (v.1.2.9) with default parameters, and the quality of the assembly was assessed via Quast (v5.2.0) (Gurevich et al., 2013; Li et al., 2015). Contigs (≥ 200 bp) were translated into protein-coding open reading frames (ORFs) via Prodigal (v.2.6.3) (Hyatt et al., 2010). A non-redundant gene catalog was generated by clustering ORFs at a sequence identity of 90% and a query coverage of 95% via CD-HIT (v4.8.1) (Li and Godzik, 2006). The ABGs were annotated by searching against the prebuilt arsenic biotransformation gene database (AsgeneDB) (Song et al., 2022) using Diamond (Buchfink et al., 2015) with a query coverage ≥ 80% and an identity threshold ≥ 50%. The abundance of ABGs was quantified via Salmon (v0.13.1) (Patro et al., 2017) and normalized to genes per million (GPM) based on the gene length and sequencing depth (Shaffer et al., 2020; Zhang et al., 2023). The taxonomic compositions of ABGs were generated using the “Asgene” R package.1

Furthermore, the clean reads obtained from dust, soil, and seawater samples were co-assembled via Megahit (v.1.2.9) (Li et al., 2015), and the co-assembled contigs were binned into metagenome-assembled genomes (MAGs) via the MetaWRAP pipeline (v1.3.2) (Uritskiy et al., 2018). The resulting MAGs were subjected to further refinement through MetaWRAP and were subsequently dereplicated via dRep (v3.5.0) (Olm et al., 2017). The assessment of MAG quality was conducted using CheckM (v.1.0.12), and only MAGs with high completeness (> 80%) and low contamination (< 10%) were retained for further analysis (Parks et al., 2015). The quantification of the relative abundance and the taxonomic annotation for each MAG were carried out using CoverM (v0.4.0) and the Genome Taxonomy Database Toolkit (GTDB-Tk v2.1.0), respectively (Chaumeil et al., 2019). Both a BLASTP-based search and a hidden Markov model (HMM)-based search were used for ABG annotation in the MAGs (Lavy et al., 2019; Song et al., 2022).

2.5 Statistical analyses

Alpha and beta diversity analyses were performed via the R ‘vegan’ package (v2.5–7) (Oksanen et al., 2011). The sources of the bacterial communities and ABGs in the dust samples were tracked using fast expectation–maximization microbial source tracking (FEAST), a highly efficient community-based microbial source tracking method that employs an expectation–maximization algorithm and demonstrates significant advantages in both computational efficiency and estimation accuracy (Shenhav et al., 2019). The FEAST analysis was performed via the R ‘FEAST’ package (v0.1.0). Key bacterial taxa and ABGs indicative of the different environments were identified by random forest classification models via the R ‘randomForest’ package (v4.6–14) (Liaw and Wiener, 2002). Co-occurrence networks depicting interactions between bacterial communities (ASVs) and ABGs were constructed via Spearman correlation analysis (|r| > 0.50 and p < 0.05) and visualized via Gephi (v0.10.1) (Bastian et al., 2009). Most of the figures were generated via Origin 2024.

3 Results

3.1 Total concentrations and speciation of arsenic

The concentrations of total arsenic in the dust, soil, and seawater samples ranged from 8.49 to 57.18 mg·kg−1, 3.18 to 10.41 mg·kg−1, and 1.50 to 2.49 μg·L−1, respectively (Supplementary Table 1). Additionally, the coefficient of variation (CV) for arsenic in the dust samples was markedly greater at 97.48%, whereas it was 37.99% in the soil samples and 17.43% in the seawater samples. Inorganic arsenic species [As(V) and As(III)] constituted more than 95% of the total arsenic species detected in dust and soil samples (Supplementary Table 2). Furthermore, in dust samples, TMAs(V)O was the main methylated arsenic species, whereas in soil samples, MAs(V) and DMAs(V) were the predominant species identified.

3.2 Variations in bacterial community structures and ABG profiles

Using 16S rRNA gene sequencing, the bacterial community structures in dust, soil, and seawater samples were characterized with an average of 63,533 high-quality sequences per sample. All of the rarefaction curves approached the plateau phase (Supplementary Figure 1), indicating that the diversity was adequately captured by the sequenced reads across different samples. Pseudomonadota (27.2–54.7%) and Actinomycetota (18.2–23.6%) were the dominant phyla across all three sample types (Figure 2A). Notably, Cyanobacteriota and Deinococcota constituted significant proportions of the bacterial communities in dust samples. The bacterial taxonomic structure at the class level varied among the different sample types, with Alphaproteobacteria (40.2%) being predominant in dust, Actinomycetes (14.1%) being predominant in soil, and Gammaproteobacteria (27.0%) being predominant in seawater. At the class level, the abundance of Flavobacteriaceae was significantly (p < 0.05) positively correlated with the proportion of As(III) but negatively correlated with the proportion of As(V) in dust and soil (Supplementary Figure 2). The bacterial alpha diversity index (Shannon index) significantly differed (p < 0.05) among the sample types, with soil exhibiting greater diversity than dust or seawater (Figure 2B). Furthermore, bacterial beta diversity significantly varied among dust, soil, and seawater samples (Bray-Curtis PCoA, Adonis test; p < 0.001; R2 = 0.629; Figure 2C). Specifically, the bacterial communities in the seawater samples were distinct from those in the dust and soil samples along the first principal coordinate axis (PCo1), whereas those in the dust and soil samples were further differentiated along the second principal coordinate axis (PCo2).

Figure 2
A: Sankey diagram showing the distribution of bacterial phyla and classes in dust, soil, and seawater. B: Box plot depicting Shannon diversity indices for dust, soil, and seawater, highlighting significant differences. C: Principal Coordinates Analysis (PCoA) plot displaying beta diversity of bacterial communities from dust, soil, and seawater.

Figure 2. Taxonomic compositions and diversity of bacterial communities in dust, soil, and seawater samples. Taxonomic compositions of the bacterial communities (A). Alpha diversity (Shannon index) of the bacterial communities (B). The Wilcoxon rank-sum test was used to compare the differences in Shannon indices of the bacterial communities among dust, soil, and seawater samples (*p < 0.05, **p < 0.01, ***p < 0.001). Beta diversity (principal coordinate analysis (PCoA) based on the Bray–Curtis distance) of bacterial communities (C).

Using metagenomic sequencing, 28 ABGs were identified across dust, soil, and seawater samples (Figure 3A). The regulatory gene (arsR) presented the highest abundance among the ars genes in dust (3530.02 ± 789.92 GPM), soil (5317.66 ± 459.04 GPM), and seawater (2697.69 ± 193.92 GPM). Significant positive correlations (p < 0.001) were observed among the abundances of different categories of ABGs, except for arsenic respiration (Supplementary Figure 3). The abundance of the arsC gene exhibited a significantly positive correlation (p < 0.01) with the proportion of As(III) in dust and soil, whereas the abundance of the arsM gene showed no significant correlation (p > 0.05) with the proportion of methylated arsenic in dust and soil (Supplementary Figure 4). Overall, the abundance of most ABGs varied across different sample types in the order soil > dust > seawater. However, the Shannon indices of the ABGs were significantly (p < 0.05) lower in the soil samples than in the dust and seawater samples (Figure 3B). The beta diversity of the ABGs also significantly varied among the dust, soil, and seawater samples (Bray-Curtis PCoA, Adonis test; p < 0.001; R2 = 0.713; Figure 3C). Specifically, the soil and seawater samples each produced two clearly separated clusters, whereas the dust samples formed one mildly dispersed cluster. Furthermore, 92.76%, 5.57%, and 1.66% of the ABG sequences were classified as bacteria, archaea, and eukaryota, respectively (Figure 3D). Pseudomonadota, Methanobacteriota, and Ascomycota emerged as the most dominant bacterial, archaeal, and eukaryotic carriers of ABGs, respectively. Notably, Cyanobacteriota were significantly more important contributors to ABGs in dust than in soil and seawater. In terms of the distribution of ABGs in microorganisms, ABGs, including the acr3, arsB, arsC, gstB, arsH, and arsM genes, were widely distributed across bacteria, archaea, and eukaryotes.

Figure 3
Chart consisting of four graphs analyzing arsenic biotransformation genes (ABGs) across different environments. Graph A displays the abundance of ABGs in dust, soil, and seawater. Graph B presents Shannon diversity indices of ABGs for these environments, with significant differences marked. Graph C shows a principal coordinates analysis (PCoA) plot displaying beta diversity of ABGs from dust, soil, and seawater. Graph D shows the count of ABGs across different taxonomic groups. Different colors represent different sample types: dust, soil, and seawater.

Figure 3. Abundance, diversity, and taxonomic compositions of arsenic biotransformation genes (ABGs). The abundances (genes per million, GPM) of ABGs (A). Alpha diversity (Shannon index) of ABGs (B). The Wilcoxon rank-sum test was used to compare the differences in abundances and Shannon indices of ABGs among dust, soil, and seawater samples (*p < 0.05, **p < 0.01, ***p < 0.001). Beta diversity (principal coordinate analysis (PCoA) based on the Bray–Curtis distance) of ABGs (C). Taxonomic compositions of ABGs at the phylum level (D).

3.3 Bacterial and ABG source tracking

Since dust, soil, and seawater may achieve a certain degree of balance after coexisting for long periods at the same site, FEAST analysis was employed to assess their mutual contribution rates for bacteria and ABGs by assuming dust to be a ‘sink’ or a ‘source’ (Figures 4A,B). When dust was assumed to be the ‘sink’, approximately 36.4% of the dust-associated bacteria were sourced from soil, and 2.9% were sourced from seawater. Similarly, 42.6% of the ABGs in dust originated from soil, and 9.1% originated from seawater. When dust was considered the ‘source’, approximately 30.3% and 3.7% of the bacterial communities in soil and seawater, respectively, along with 34.8% and 10.6% of the ABGs in soil and seawater, respectively, originated from dust.

Figure 4
Two diagrams labeled A and B show Fast expectation–maximization microbial source tracking (FEAST) analysis for bacteria and ABGs. In A, contributions are 2.9%, 3.7%, 30.3%, and 36.4%. In B, contributions are 9.1%, 10.6%, 34.8%, and 42.6%. Arrows indicate the relationship of sources and sinks between dust, seawater, and soil.

Figure 4. Fast expectation–maximization microbial source tracking (FEAST) analysis for bacteria (A) and ABGs (B). The direction of the arrows indicates the relationship between sources and sinks, and the percentages reflect the contribution of each source.

3.4 Key bacterial taxa and ABGs

Key bacterial taxa were identified using a random forest classification model and ranked in descending order of their importance (Figure 5A). Among these key bacterial taxa, the Acetobacteraceae (Pseudomonadota), Deinococcaceae (Deinococcota), and Chroococcidiopsidaceae (Cyanobacteriota) presented significantly greater relative abundances (p < 0.001) in dust than in soil and seawater. Furthermore, network analysis was conducted to identify the highly connected taxa within the bacterial community. The bacterial co-occurrence networks in dust, soil, and seawater all displayed scale-free and small-world characteristics (Supplementary Table 3; small-world coefficient (σ) > 1). Compared with those in dust and seawater, the bacterial networks in soil presented a greater number of nodes and a lower proportion of positive correlations (Supplementary Figure 5). A total of one module hub and three connector hubs were identified in the dust network (Supplementary Figure 6). These hubs, which represented keystone taxa, belonged to Acetobacteraceae, Deinococcaceae, and Rhodobacteraceae (Supplementary Table 4).

Figure 5
Panel A shows a dot and bar plot comparing relative abundance and mean decrease accuracy of bacterial families across dust, soil, and seawater, categorized by phylum. Panel B presents shows a dot and bar plot comparing relative proportion and mean decrease accuracy of arsenic biotransformation genes (ABGs) across these environments. Panels C, D, and E display network diagrams, each indicating nodes and edges with varying proportions of positive and negative associations, showing interactions among different bacteria and arsenic biotransformation genes (ABGs).

Figure 5. Key bacterial taxa and arsenic biotransformation genes (ABGs) and the co-occurrence networks of bacterial taxa and ABGs across dust, soil, and seawater. Random forest classification models identifying key bacterial taxa at the family level (A) and ABGs (B). The panels on the left illustrate the relative abundance of these families or the relative proportions of ABGs. The panels on the right show the mean decrease in accuracy for each family or ABG, indicating their importance in the model (*p < 0.05 and **p < 0.01). Network analysis revealing the co-occurrence patterns between ABGs and bacterial communities in dust (C), soil (D), and seawater (E) based on Spearman correlation analysis (|r| > 0.50, p < 0.05).

Fourteen key ABGs were also identified via a random forest classification model (Figure 5B). Among these key ABGs, the relative proportion of the acr3 gene, which is involved in arsenic efflux, was significantly greater (p < 0.001) in dust (4.78%) than in soil (0.89%) or seawater (0.86%). Furthermore, network analysis revealed significant co-occurrence relationships between ABGs and bacterial communities in dust, soil, and seawater, with the arsR, arsC, acr3, and pstSCAB genes emerging as critical nodes (Figures 5CE). The aioB and arsM genes were also identified as key nodes closely associated with bacterial communities in the co-occurrence network of dust. Rubellimicrobium (belonging to Roseobacteraceae), Geodermatophilus (belonging to Geodermatophilaceae), and Craurococcus (belonging to Acetobacteraceae) were closely associated with ABGs in the dust network.

Furthermore, at the genome level, 32, 18, and 35 MAGs with completeness above 80% and contamination below 10% were recovered from dust, soil, and seawater, respectively (Supplementary Figure 7). The predominant phyla for bacterial and archaeal MAGs were Pseudomonadota and Nitrososphaerota, respectively. The relative abundance and presence/absence of ABGs were further analyzed for each MAG in dust (Figure 6; Supplementary Figure 8). The MAG related to Acetobacteraceae (Bin 314) presented the highest relative abundance in dust. Additionally, Chroococcidiopsidaceae MAGs (Bins 467, 541, and 591) presented high relative abundances in dust. Furthermore, although a diverse array of ABGs was observed across different MAGs, a majority of these high-abundance MAGs in dust harbored key ABGs, including the arsR, arsC, aioA, aioB, and arsM genes.

Figure 6
Phylogenetic tree diagram displaying the detected bacterial metagenome-assembled genomes (MAGs) in dust. Branches are color-coded by phylum. The heatmap illustrates the presence (blue) or absence (light blue) of arsenic biotransformation genes (ABGs) for each MAG. Relative abundance of each MAG is indicated with the bars.

Figure 6. Phylogenetic tree of the detected bacterial metagenome-assembled genomes (MAGs) with completeness greater than 80% and contamination less than 10% in dust. The line bar indicates a tree scale of 0.2, and the purple dots represent bootstrap values. The heatmap illustrates the presence (blue) or absence (light blue) of arsenic biotransformation genes (ABGs) for each MAG. The bars enclosed by dashed rings depict the relative abundance of each MAG in dust.

4 Discussion

In this study, dust samples from Xiamen, China, presented significantly higher (p < 0.05) total arsenic concentrations than soil samples did, which is consistent with findings from Nevada, United States (Soukup et al., 2012), and Ahvaz, Iran (Ghanavati et al., 2018). The elevated arsenic levels in dust may be attributed to the fact that dust particles are generally smaller in size and more susceptible to contamination from anthropogenic activities than soil particles are (Soukup et al., 2012). Additionally, the higher CV values for arsenic in the dust samples further suggest that arsenic concentrations in dust are more substantially influenced by variations attributable to diverse anthropogenic activities (Bullard and Baddock, 2019; Chen et al., 2018a). Inorganic arsenic species dominated in both dust and soil samples, but methylated arsenic species were also detected. The detection of MAs(V) and DMAs(V) in soils is possible because of microorganism-mediated biomethylation processes, whereas the occurrence of TMAs(V)O in dust likely results from the atmospheric oxidation of volatile trimethylarsine (TMAs) emitted from terrestrial and marine environments (Qiao et al., 2023; Savage et al., 2019; Yan et al., 2024; Zhang et al., 2015; Zhao et al., 2013).

The diversity of both the bacterial communities and the ABGs significantly differed (p < 0.05) among the dust, soil, and seawater samples. The Shannon indices of bacterial communities were significantly (p < 0.05) greater in soil than in dust and seawater, which is consistent with the recognized role of soil as one of the most extensive reservoirs of biological diversity on Earth (Swift et al., 1998). Furthermore, the majority of ABGs showed significantly (p < 0.05) greater abundances in soil than in dust and seawater, possibly because soil can provide a nutrient-rich medium that supports diverse microbial activities (Yadav et al., 2021). However, the Shannon indices of ABGs were significantly greater (p < 0.05) in dust and seawater than in soil. This observation may be attributed to the harsher environmental conditions in both dust and seawater, with limited nutrient availability likely constraining microbial diversity while promoting gene diversity (Fierer et al., 2012). Compared with those in soil, bacterial networks in dust and seawater presented fewer nodes and higher proportions of positive correlations, further suggesting that bacterial communities exhibit more cooperative relationships when subjected to harsh environmental conditions such as those in dust and seawater (Li et al., 2022). Notably, bacterial beta diversity analysis revealed that bacterial communities in dust are more closely related to those in soil than to those in seawater, which was supported by the greater contribution of soil compared with seawater as a source of bacteria in dust. Additionally, in the PCoA scatter plots visualizing the beta diversity of the ABGs, the dust samples displayed a more dispersed pattern than the soil and seawater samples did, indicating considerable spatial variability in the profiles of the ABGs among the dust samples (Yan et al., 2024).

The source-sink relationships of both bacteria and ABGs were investigated using the FEAST analysis. When dust is considered the ‘sink’, the great contributions of soil to both bacteria and ABGs in dust indicate that soil serves as the dominant source for dust-associated microorganisms. Consistent with these results, previous studies have reported that terrestrial bacterial sources contribute significantly to airborne and dust-associated bacteria (Qi et al., 2023; Salawu-Rotimi et al., 2021). Notably, soil exhibited a greater contribution rate for ABGs than for bacteria in dust. A previous study reported that fungal communities, rather than bacterial communities, in dust were more strongly associated with soil microbiomes since bacteria experience more selective fractionation during the generation of dust (Salawu-Rotimi et al., 2021). Therefore, given that ABGs are associated with both prokaryotic and eukaryotic microorganisms, it is unsurprising that soil significantly contributes to the abundance of ABGs hosted by both bacteria and fungi in dust. Additionally, the arsenic in dust has been reported to originate from terrestrial sources (Ma et al., 2021; Yan et al., 2024). Consequently, soil could not only be a significant source of arsenic but also serve as a considerable contributor to the associated genes in dust. Additionally, although the contribution rate was not high, seawater could also contribute to the presence of some specific bacteria in dust. For instance, Cyanobacteriota was detected as an important bacterial community in dust samples in this study. A previous study reported that the presence of Cyanobacteriota adhering to atmospheric particulate matter is likely attributable to sea and lake spray aerosols, which facilitate the transport of Cyanobacteriota from aquatic environments into the atmosphere (Nie et al., 2023). Both soil and seawater in Tokyo Bay were also identified as consistent and predominant sources of airborne microorganisms in Tokyo, Japan (Uetake et al., 2019). When dust was considered the ‘source’, the contribution of dust to both bacterial communities and ABGs in soil was more significant than that in seawater. Intriguingly, the contribution rates for bacteria and ABGs from dust to seawater were greater than those from seawater to dust, indicating that dust predominantly serves as a primary source of both bacteria and ABGs in seawater. Previous studies have shown that a variety of microorganisms can be transported via dust and aerosol dispersal, influencing marine microbial communities and biogeochemical processes (Berman-Frank and Rahav, 2012; Rahav et al., 2016).

Pseudomonadota were the dominant bacterial phylum according to 16S rRNA gene sequencing and predominated among MAGs recovered through metagenomic binning across dust, soil, and seawater samples, which is consistent with the well-known prevalence of Pseudomonadota in diverse environments (Li et al., 2022; Suh et al., 2015; Wang et al., 2019; Xu et al., 2014). Furthermore, the ubiquitous Pseudomonadota were also the dominant bacterial carriers of ABGs, emphasizing their importance in arsenic biotransformation (Saunders et al., 2019; Zhang et al., 2015). Cyanobacteriota, often identified in bioaerosols (Hu et al., 2020; Liu et al., 2018; Wei et al., 2020), accounted for relatively high proportions of bacterial communities in dust samples. The MAGs affiliated with Cyanobacteriota also presented high relative abundances in dust samples, which is consistent with their large proportions in bacterial communities in dust. Deinococcota also constituted a significant component of bacterial communities in dust, as certain genera within this phylum exhibit remarkable resistance to solar radiation (Xu et al., 2017). Furthermore, according to the random forest analysis, Acetobacteraceae, belonging to Pseudomonadota, Deinococcaceae, belonging to Deinococcota, and Chroococcidiopsidaceae, belonging to Cyanobacteriota, which exhibited high abundances in dust samples, were identified as key taxa. Notably, the keystone taxa within the bacterial co-occurrence network predominantly belonged to Acetobacteraceae and Deinococcaceae. Therefore, Acetobacteraceae and Deinococcaceae not only serve as key taxa characterized by high abundances in dust but also function as keystone taxa with high connectivity within bacterial communities in dust environments. Furthermore, the Acetobacteraceae MAGs (Bin 314) and Chroococcidiopsidaceae MAGs (Bins 467, 541, and 591) presented high relative abundances in dust samples. These two bacterial families are commonly found in extreme environments, such as deserts, underscoring their critical role in adaptation to soil conditions with limited water (Al Ashhab et al., 2021; Caiola and Billi, 2007; Yang et al., 2017). Consequently, they were unsurprisingly key microbial constituents in dust, which can also be considered an arid environment.

Across dust, soil, and seawater samples, highly abundant ABGs involved in arsenic uptake (pstSCAB) (Rosen, 2002), oxidation (moeA and arsH) (Zhu et al., 2017; Chen et al., 2015), reduction (arsC and gstB) (Chrysostomou et al., 2015), As(III) efflux (arsB and acr3) (Garbinski et al., 2019), and methylation (arsM) (Qin et al., 2006) were detected. Additionally, these ABG categories exhibited significant positive correlations with each other (p < 0.001). These findings indicate that multiple arsenic transformation processes commonly coexist within atmospheric, terrestrial, and aquatic environments (Zhu et al., 2017). Additionally, the wide distribution of the acr3, arsB, arsC, gstB, arsH, and arsM genes across bacteria, archaea, and eukaryotes aligns with previous studies reporting their broad presence in prokaryotes and eukaryotes (Chen et al., 2020; Li et al., 2024). Notably, the arrAB gene, which is responsible for As(V) respiration, was detected at extremely low levels in dust, which can be attributed to the prolonged exposure of dust to abundant oxygen in the air (Min et al., 2022; Yan et al., 2024). According to the random forest analysis, the acr3 gene, which presented a significantly greater relative proportion (p < 0.001) in dust samples than in soil and seawater samples, was identified as one of the key ABGs. Therefore, arsenic efflux, widely acknowledged as an effective arsenic detoxification mechanism (Castro-Severyn et al., 2022), plays a pivotal role in arsenic biotransformation and resistance, particularly in dust samples with relatively high arsenic concentrations. Moreover, a significantly positive correlation (p < 0.01) was observed between the abundance of the arsC gene and the proportion of As(III) in dust and soil. This observation is expected since the arsC gene was the key gene involving in As(V) reduction to As(III) (Zhu et al., 2017). Furthermore, the abundance of Flavobacteriaceae exhibited a significantly (p < 0.05) positive correlation with the proportion of As(III) and a negative correlation with the proportion of As(V) in dust and soil. This suggests that Flavobacteriaceae may play an important role in arsenic reduction in the environment. Supporting this, a novel species within the family Flavobacteriaceae was identified as a multi-metal-resistant bacterium harboring the arsC gene (Li S. et al., 2023). However, no significant correlation (p > 0.05) was observed between the abundance of the arsM gene and the proportion of methylated arsenic in dust and soil. Similar findings have been also reported in previous studies, which may be attributed to the influence of other arsenic biotransformation pathways (e.g., arsenic demethylation) and chemical processes on the proportion of methylated arsenic in the environment (Qiao et al., 2023; Yan et al., 2024).

Co-occurrence networks between bacterial taxa and ABGs, together with MAG binning, were utilized to explore the relationships between bacteria and ABGs. Craurococcus, a genus within Acetobacteraceae, was closely associated with a variety of ABGs, including acr3, arsA, arsC, and arsM genes in the co-occurrence network between bacterial taxa and ABGs in dust. Furthermore, Acetobacteraceae MAGs were also found to harbor these ABGs. These findings further suggest that Acetobacteraceae are likely important hosts for ABGs within dust environments (Cernava et al., 2018). The arsR, arsC, acr3, and pstSCAB genes represented key nodes in the co-occurrence networks linking ABGs and bacterial communities across dust, soil, and seawater. This finding indicates that after uptake via the PstSCAB system, As(V) can be reduced to As(III) by ArsC and subsequently effluxed by Acr3, thereby forming an essential arsenic detoxification pathway (Garbinski et al., 2019). Furthermore, the arsR, arsC, acr3, aioB, and arsM genes were not only identified as key nodes closely associated with bacterial communities in the co-occurrence network of dust but also harbored by the majority of high-abundance MAGs in dust. These results indicate that dust may possess the capacity for multiple types of arsenic biotransformation, including As(V) reduction, As(III) oxidation, and As(III) methylation (Yan et al., 2024).

5 Conclusion

This study revealed that the bacterial community structures and ABG profiles in dust were significantly different from those in soil and seawater. Soil, rather than seawater, was the main source of both bacteria and ABGs in dust. A variety of bacteria adapted to arid environments, particularly Acetobacteraceae showed high abundance and network connectivity in dust and was identified as a key taxon according to the random forest. The MAGs with high relative abundance in dust were also affiliated with Acetobacteraceae and harbored several key ABGs. Acetobacteraceae could be a potentially important contributor to arsenic biotransformation in dust. Furthermore, ABGs, particularly the arsR, arsC, acr3, aioB, and arsM genes, exhibited close connections with bacterial taxa in the co-occurrence networks between ABGs and bacterial communities in dust and were prevalent within the high-abundance MAGs in dust. Further research is essential to comprehensively investigate the capacity of arsenic biotransformation within dust.

Data availability statement

The 16S rRNA gene sequencing data and metagenomic sequencing data have been deposited in the NCBI SRA under BioProject PRJNA1214267 and PRJNA1214535, respectively.

Author contributions

YiY: Data curation, Formal analysis, Investigation, Methodology, Writing – original draft. Y-TL: Data curation, Writing – review & editing. G-RH: Supervision, Writing – review & editing. R-LY: Writing – review & editing. X-HS: Funding acquisition, Supervision, Writing – review & editing. YuY: Conceptualization, Funding acquisition, Methodology, Project administration, Resources, Supervision, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the National Natural Science Foundation of China (42177103), the Natural Science Foundation of Fujian Province (2023 J05050), the Fundamental Research Funds for the Central Universities (ZQN-1015), Xiamen Key Science and Technology Program (Future Industries Field) (3502Z20234041), and the Youth Science and Technology Innovation Program of Xiamen Ocean and Fisheries Development Special Funds (23YYZP017QCB25).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The authors declare that no Gen AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2025.1607082/full#supplementary-material

Footnotes

References

Al Ashhab, A., Meshner, S., Alexander-Shani, R., Dimerets, H., Brandwein, M., Bar-Lavan, Y., et al. (2021). Temporal and spatial changes in phyllosphere microbiome of acacia trees growing in arid environments. Front. Microbiol. 12:656269. doi: 10.3389/fmicb.2021.656269

PubMed Abstract | Crossref Full Text | Google Scholar

Andres, J., and Bertin, P. N. (2016). The microbial genomics of arsenic. FEMS Microbiol. Rev. 40, 299–322. doi: 10.1093/femsre/fuv050

PubMed Abstract | Crossref Full Text | Google Scholar

Archer, S. D. J., Lee, K. C., Caruso, T., Alcami, A., Araya, J. G., Cary, S. C., et al. (2023). Contribution of soil bacteria to the atmosphere across biomes. Sci. Total Environ. 871:162137. doi: 10.1016/j.scitotenv.2023.162137

PubMed Abstract | Crossref Full Text | Google Scholar

Balakrishna, G., Pervez, S., and Bisht, D. (2011). Source apportionment of arsenic in atmospheric dust fall out in an urban residential area, Raipur, Central India. Atmos. Chem. Phys. 11, 5141–5151. doi: 10.5194/acp-11-5141-2011

Crossref Full Text | Google Scholar

Barberan, A., Ladau, J., Leff, J. W., Pollard, K. S., Menninger, H. L., Dunn, R. R., et al. (2015). Continental-scale distributions of dust-associated bacteria and fungi. Proc. Natl. Acad. Sci. U. S. A. 112, 5756–5761. doi: 10.1073/pnas.1420815112

PubMed Abstract | Crossref Full Text | Google Scholar

Bastian, M., Heymann, S., and Jacomy, M. (2009). Gephi: an open source softwarefor exploring and manipulating networks. Proc. Int. AAAI Conf. Web Soc. Med. 3, 361–362. doi: 10.1609/icwsm.v3i1.13937

Crossref Full Text | Google Scholar

Ben Fekih, I., Zhang, C., Li, Y. P., Zhao, Y., Alwathnani, H. A., Saquib, Q., et al. (2018). Distribution of arsenic resistance genes in prokaryotes. Front. Microbiol. 9:2473. doi: 10.3389/fmicb.2018.02473

PubMed Abstract | Crossref Full Text | Google Scholar

Berman-Frank, I., and Rahav, E. (2012). “Dinitrogen fixation as a source for new production in the Mediterranean Sea: a review” in Life in the Mediterranean Sea: A look at habitat changes. ed. N. E. Stambler (New York: Nova Science Publishers), 199–226.

Google Scholar

Bhattacharya, P., Frisbie, S., Smith, E., Naidu, R., Jacks, G., and Sarkar, B. (2002). “Arsenic in the environment: A global perspective” in Handbook of heavy metals in the environment. ed. P. Bhattacharya (New York: Marcell Dekker, Inc.), 147–215.

Google Scholar

Bowell, R. J., Alpers, C. N., Jamieson, H. E., Nordstrom, D. K., and Majzlan, J. (2014). The environmental geochemistry of arsenic – An overview. Rev. Mineral. Geochem 79, 1–16. doi: 10.2138/rmg.2014.79.1

Crossref Full Text | Google Scholar

Bowers, R. M., Sullivan, A. P., Costello, E. K., Collett, J. L. Jr., Knight, R., and Fierer, N. (2011). Sources of bacteria in outdoor air across cities in the midwestern United States. Appl. Environ. Microbiol. 77, 6350–6356. doi: 10.1128/AEM.05498-11

PubMed Abstract | Crossref Full Text | Google Scholar

Buchfink, B., Xie, C., and Huson, D. H. (2015). Fast and sensitive protein alignment using DIAMOND. Nat. Methods 12, 59–60. doi: 10.1038/nmeth.3176

PubMed Abstract | Crossref Full Text | Google Scholar

Bullard, J., and Baddock, M. (2019). “Dust: sources, entrainment, transport” in Aeolian geomorphology: A new introduction. eds. J. Bullard and M. Baddock (London: Longman), 81–106.

Google Scholar

Caiola, M. G., and Billi, D. (2007). “Chroococcidiopsis from desert to mars” in Algae and cyanobacteria in extreme environments. ed. J. Seckbach, vol. 11 (Berlin: Springer Science and Business Media).

Google Scholar

Castro-Severyn, J., Pardo-Esté, C., Araya-Durán, I., Gariazzo, V., Cabezas, C., Valdés, J., et al. (2022). Biochemical, genomic and structural characteristics of the Acr3 pump in Exiguobacterium strains isolated from arsenic-rich Salar de Huasco sediments. Front. Microbiol. 13:1047283. doi: 10.3389/fmicb.2022.1047283

PubMed Abstract | Crossref Full Text | Google Scholar

Cernava, T., Vasfiu, Q., Erlacher, A., Aschenbrenner, I. A., Francesconi, K., Grube, M., et al. (2018). Adaptions of lichen microbiota functioning under persistent exposure to arsenic contamination. Front. Microbiol. 9:2959. doi: 10.3389/fmicb.2018.02959

PubMed Abstract | Crossref Full Text | Google Scholar

Chaumeil, P. A., Mussig, A. J., Hugenholtz, P., and Parks, D. H. (2019). GTDB-Tk: a toolkit to classify genomes with the genome taxonomy database. Bioinformatics 36, 1925–1927. doi: 10.1093/bioinformatics/btz848

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, J., Bhattacharjee, H., and Rosen, B. P. (2015). ArsH is an organoarsenical oxidase that confers resistance to trivalent forms of the herbicide monosodium methylarsenate and the poultry growth promoter roxarsone. Mol. Microbiol. 96, 1042–1052. doi: 10.1111/mmi.12988

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, Q. Y., and Costa, M. (2021). Arsenic: a global environmental challenge. Annu. Rev. Pharmacol. Toxicol. 61, 47–63. doi: 10.1146/annurev-pharmtox-030220-013418

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, S., Jiang, N., Huang, J., Xu, X., Zhang, H., Zang, Z., et al. (2018a). Quantifying contributions of natural and anthropogenic dust emission from different climatic regions. Atmos. Environ. 191, 94–104. doi: 10.1016/j.atmosenv.2018.07.043

Crossref Full Text | Google Scholar

Chen, S. C., Sun, G. X., Yan, Y., Konstantinidis, K. T., Zhang, S. Y., Deng, Y., et al. (2020). The great oxidation event expanded the genetic repertoire of arsenic metabolism and cycling. Proc. Natl. Acad. Sci. 117, 10414–10421. doi: 10.1073/pnas.2001063117

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, S., Zhou, Y., Chen, Y., and Gu, J. (2018b). fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34, i884–i890. doi: 10.1093/bioinformatics/bty560

PubMed Abstract | Crossref Full Text | Google Scholar

Chrysostomou, C., Quandt, E. M., Marshall, N. M., Stone, E., and Georgiou, G. (2015). An alternate pathway of arsenate resistance in E. coli mediated by the glutathione S-transferase Gst B. ACS Chem. Biol. 10, 875–882. doi: 10.1021/cb500755j

PubMed Abstract | Crossref Full Text | Google Scholar

Ding, L. J., Zhou, X. Y., and Zhu, Y. G. (2020). Microbiome and antibiotic resistome in household dust from Beijing. China. Environ. Int. 139:105702. doi: 10.1016/j.envint.2020.105702

PubMed Abstract | Crossref Full Text | Google Scholar

Edgar, R. C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461. doi: 10.1093/bioinformatics/btq461

PubMed Abstract | Crossref Full Text | Google Scholar

Edgar, R. C. (2016). UNOISE2: improved error-correction for Illumina 16S and ITS amplicon sequencing. bio Rxiv 2016:081257. doi: 10.1101/081257

Crossref Full Text | Google Scholar

Fierer, N., Ferrenberg, S., Flores, G. E., González, A., Kueneman, J., Legg, T., et al. (2012). From animalcules to an ecosystem: application of ecological concepts to the human microbiome. Annu. Rev. Ecol. Evol. Syst. 43, 137–155. doi: 10.1146/annurev-ecolsys-110411-160307

Crossref Full Text | Google Scholar

Garbinski, L. D., Rosen, B. P., and Chen, J. (2019). Pathways of arsenic uptake and efflux. Environ. Int. 126, 585–597. doi: 10.1016/j.envint.2019.02.058

PubMed Abstract | Crossref Full Text | Google Scholar

Ghanavati, N., Nazarpour, A., and De Vivo, B. (2018). Ecological and human health risk assessment of toxic metals in street dusts and surface soils in Ahvaz, Iran. Environ. Geochem. Health 41, 875–891. doi: 10.1007/s10653-018-0184-y

PubMed Abstract | Crossref Full Text | Google Scholar

Guo, T., Li, L., Zhai, W., Xu, B., Yin, X., He, Y., et al. (2019). Distribution of arsenic and its biotransformation genes in sediments from the East China Sea. Environ. Pollut. 253, 949–958. doi: 10.1016/j.envpol.2019.07.091

PubMed Abstract | Crossref Full Text | Google Scholar

Gurevich, A., Saveliev, V., Vyahhi, N., and Tesler, G. (2013). QUAST: quality assessment tool for genome assemblies. Bioinformatics 29, 1072–1075. doi: 10.1093/bioinformatics/btt086

PubMed Abstract | Crossref Full Text | Google Scholar

Hu, Z., Liu, H., Zhang, H., Zhang, X., Zhou, M., Lou, L., et al. (2020). Temporal discrepancy of airborne total bacteria and pathogenic bacteria between day and night. Environ. Res. 186:109540. doi: 10.1016/j.envres.2020.109540

PubMed Abstract | Crossref Full Text | Google Scholar

Huang, Y., Miu, Q., Kwong, R. W., Zhang, D., Fan, Y., Zhou, M., et al. (2024). Leveraging the one health concept for arsenic sustainability. Ecoenviron. Health 3, 392–405. doi: 10.1016/j.eehl.2024.02.006

Crossref Full Text | Google Scholar

Hyatt, D., Chen, G.-L., LoCascio, P. F., Land, M. L., Larimer, F. W., and Hauser, L. J. (2010). Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics 11:119. doi: 10.1186/1471-2105-11-119

PubMed Abstract | Crossref Full Text | Google Scholar

Jain, C. K., and Ali, I. (2000). Arsenic: occurrence, toxicity and speciation techniques. Water Res. 34, 4304–4312. doi: 10.1016/S0043-1354(00)00182-2

Crossref Full Text | Google Scholar

Kara, S., Chormey, D.S., Saygılar, A., and Bakırdere, S. (2021). Arsenic speciation in rice samples for trace level determination by high performance liquid chromatography-inductively coupled plasma-mass spectrometry. Food Chem. 356:129706. doi: 10.1016/j.foodchem.2021.129706

Crossref Full Text | Google Scholar

Lavy, A., Matheus Carnevali, P. B., Keren, R., Bill, M., Wan, J., Tokunaga, T. K., et al. (2019). Taxonomically and metabolically distinct microbial communities with depth and across a hillslope to riparian zone transect. bioRxiv 2019:768572. doi: 10.1101/768572

Crossref Full Text | Google Scholar

Leppänen, H. K., Täubel, M., Jayaprakash, B., Vepsäläinen, A., Pasanen, P., and Hyvärinen, A. (2018). Quantitative assessment of microbes from samples of indoor air and dust. J. Expo. Sci. Environ. Epidemiol. 28, 231–241. doi: 10.1038/jes.2017.24

PubMed Abstract | Crossref Full Text | Google Scholar

Li, L., Chen, S., Xue, X., Chen, J., Tian, J., Huo, L., et al. (2024). Purifying selection drives distinctive arsenic metabolism pathways in prokaryotic and eukaryotic microbes. ISME Commun. 4:106. doi: 10.1093/ismeco/ycae106

PubMed Abstract | Crossref Full Text | Google Scholar

Li, W., and Godzik, A. (2006). Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 22, 1658–1659. doi: 10.1093/bioinformatics/btl158

PubMed Abstract | Crossref Full Text | Google Scholar

Li, H., Liu, P. Q., Luo, Q. P., Ma, J. J., Yang, X. R., Yan, Y., et al. (2022). Spatiotemporal variations of microbial assembly, interaction, and potential risk in urban dust. Environ. Int. 170:107577. doi: 10.1016/j.envint.2022.107577

PubMed Abstract | Crossref Full Text | Google Scholar

Li, D., Liu, C. M., Luo, R., Sadakane, K., and Lam, T. W. (2015). MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics 31, 1674–1676. doi: 10.1093/bioinformatics/btv033

PubMed Abstract | Crossref Full Text | Google Scholar

Li, S., Mao, S., Mu, W., Guo, B., Li, C., Zhu, Q., et al. (2023). Flavobacterium potami sp. nov., a multi-metal resistance genes harbouring bacterium isolated from shallow river silt. Antonie Van Leeuwenhoek 116, 265–280. doi: 10.1007/s10482-022-01802-z

PubMed Abstract | Crossref Full Text | Google Scholar

Li, L., Meng, D., Yin, H., Zhang, T., and Liu, Y. (2023). Genome-resolved metagenomics provides insights into the ecological roles of the keystone taxa in heavy-metal-contaminated soils. Front. Microbiol. 14:1203164. doi: 10.3389/fmicb.2023.1203164

PubMed Abstract | Crossref Full Text | Google Scholar

Liaw, A., and Wiener, A. (2002). Classification and regression by randomForest. R News 2, 18–22.

Google Scholar

Liu, H., Zhang, X., Zhang, H., Yao, X., Zhou, M., Wang, J., et al. (2018). Effect of air pollution on the total bacteria and pathogenic bacteria in different sizes of particulate matter. Environ. Pollut. 233, 483–493. doi: 10.1016/j.envpol.2017.10.070

PubMed Abstract | Crossref Full Text | Google Scholar

Ma, J. W., Wu, Y. Q., Xu, C. L., Luo, Z. X., Yu, R. L., Hu, G. R., et al. (2023). Inhibitory effect of polyethylene microplastics on roxarsone degradation in soils. J. Hazard. Mater. 454:131483. doi: 10.1016/j.jhazmat.2023.131483

PubMed Abstract | Crossref Full Text | Google Scholar

Ma, J. J., Yan, Y., Chen, X. J., Niu, Z. R., Yu, R. L., and Hu, G. R. (2021). Incorporating bioaccessibility and source apportionment into human health risk assessment of heavy metals in urban dust of Xiamen, China. Ecotoxicol. Environ. Saf. 228:112985. doi: 10.1016/j.ecoenv.2021.112985

PubMed Abstract | Crossref Full Text | Google Scholar

Min, D., Cheng, L., Liu, D. F., Liu, J. Q., Li, W. W., and Yu, H. Q. (2022). Single strain-triggered biogeochemical cycle of arsenic. Environ. Sci. Technol. 56, 16410–16418. doi: 10.1021/acs.est.2c02015

PubMed Abstract | Crossref Full Text | Google Scholar

Mukhopadhyay, R., Rosen, B. P., Phung, L. T., and Silver, S. (2002). Microbial arsenic: from geocycles to genes and enzymes. FEMS Microbiol. Rev. 26, 311–325. doi: 10.1111/j.1574-6976.2002.tb00617.x

PubMed Abstract | Crossref Full Text | Google Scholar

Nie, C., Geng, X., Zhang, R., Wang, L., Li, L., and Chen, J. (2023). Abundant cyanobacteria in autumn adhering to the heating, ventilation, and air-conditioning (HVAC) in Shanghai. Microorganisms 11:1835. doi: 10.3390/microorganisms11071835

PubMed Abstract | Crossref Full Text | Google Scholar

Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., Minchin, P. R., O’Hara, R. B., et al. (2011). Vegan: Comunity ecology package. Avaliable online at: http://CRAN.R-project.org/package=vegan.

Google Scholar

Olm, M. R., Brown, C. T., Brooks, B., and Banfield, J. F. (2017). dRep: a tool for fast and accurate genomic comparisons that enables improved genome recovery from metagenomes through de-replication. ISME J. 11, 2864–2868. doi: 10.1038/ismej.2017.126

PubMed Abstract | Crossref Full Text | Google Scholar

Oremland, R. S., and Stolz, J. F. (2003). The ecology of arsenic. Science 300, 939–944. doi: 10.1126/science.1081903

PubMed Abstract | Crossref Full Text | Google Scholar

Oremland, R. S., and Stolz, J. F. (2005). Arsenic, microbes and contaminated aquifers. Trends Microbiol. 13, 45–49. doi: 10.1016/j.tim.2004.12.002

PubMed Abstract | Crossref Full Text | Google Scholar

Orero Iserte, L., Roig-Navarro, A. F., and Hernández, F. (2004). Simultaneous determination of arsenic and selenium species in phosphoric acid extracts of sediment samples by HPLC-ICP-MS. Anal. Chim. Acta 527, 97–104. doi: 10.1016/j.aca.2004.08.001

Crossref Full Text | Google Scholar

Parks, D. H., Imelfort, M., Skennerton, C. T., Hugenholtz, P., and Tyson, G. W. (2015). Check M: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 25, 1043–1055. doi: 10.1101/gr.186072.114

PubMed Abstract | Crossref Full Text | Google Scholar

Patro, R., Duggal, G., Love, M. I., Irizarry, R. A., and Kingsford, C. (2017). Salmon provides fast and bias-aware quantification of transcript expression. Nat. Methods 14, 417–419. doi: 10.1038/nmeth.4197

PubMed Abstract | Crossref Full Text | Google Scholar

Pizarro, I., Gómez, M., Cámara, C., and Palacios, M.A. (2003). Arsenic speciation in environmental and biological samples: extraction and stability studies. Anal. Chim. Acta. 495, 85–98. doi: 10.1016/j.aca.2003.08.009

Crossref Full Text | Google Scholar

Qi, J., Huang, Z., Xue, F., Gao, Z., Maki, T., Zhang, Z., et al. (2023). Aridification alters the diversity of airborne bacteria in drylands of China. Atmos. Environ. 315:120135. doi: 10.1016/j.atmosenv.2023.120135

Crossref Full Text | Google Scholar

Qiao, J., Liu, J., Palomo, A., Bostick, B. C., Phan, K., Zheng, Y., et al. (2023). Prevalence of methylated arsenic and microbial arsenic methylation genes in paddy soils of the Mekong Delta. Environ. Sci. Technol. 57, 9675–9682. doi: 10.1021/acs.est.3c00210

PubMed Abstract | Crossref Full Text | Google Scholar

Qin, J., Rosen, B. P., Zhang, Y., Wang, G., Franke, S., and Rensing, C. (2006). Arsenic detoxification and evolution of trimethylarsine gas by a microbial arsenite S-adenosylmethionine methyltransferase. Proc. Natl. Acad. Sci. U. S. A. 103, 2075–2080. doi: 10.1073/pnas.0506836103

PubMed Abstract | Crossref Full Text | Google Scholar

Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219

PubMed Abstract | Crossref Full Text | Google Scholar

Rahav, E., Ovadia, G., Paytan, A., and Herut, B. (2016). Contribution of airborne microbes to bacterial production and N2 fixation in seawater upon aerosol deposition. Geophys. Res. Lett. 43, 719–727. doi: 10.1002/2015GL066898

Crossref Full Text | Google Scholar

Raju, N. J. (2022). Arsenic in the geo-environment: a review of sources, geochemical processes, toxicity and removal technologies. Environ. Res. 203:111782. doi: 10.1016/j.envres.2021.111782

PubMed Abstract | Crossref Full Text | Google Scholar

Rosen, B. P. (2002). Biochemistry of arsenic detoxification. FEBS Lett. 529, 86–92. doi: 10.1016/s0014-5793(02)03186-1

PubMed Abstract | Crossref Full Text | Google Scholar

Salawu-Rotimi, A., Lebre, P. H., Vos, H. C., Fister, W., Kuhn, N., Eckardt, F. D., et al. (2021). Gone with the wind: microbial communities associated with dust from emissive farmlands. Microb. Ecol. 82, 859–869. doi: 10.1007/s00248-021-01717-8

PubMed Abstract | Crossref Full Text | Google Scholar

Sánchez-Rodas, D., de la Campa, A. M. S., and Alsioufi, L. (2015). Analytical approaches for arsenic determination in air: a critical review. Anal. Chim. Acta 898, 1–18. doi: 10.1016/j.aca.2015.09.043

PubMed Abstract | Crossref Full Text | Google Scholar

Saunders, J. K., Fuchsman, C. A., McKay, C., and Rocap, G. (2019). Complete arsenic-based respiratory cycle in the marine microbial communities of pelagic oxygen-deficient zones. PNAS 116, 9925–9930. doi: 10.1073/pnas.1818349116

PubMed Abstract | Crossref Full Text | Google Scholar

Savage, L., Carey, M., Williams, P. N., and Meharg, A. A. (2019). Maritime deposition of organic and inorganic arsenic. Environ. Sci. Technol. 53, 7288–7295. doi: 10.1021/acs.est.8b06335

PubMed Abstract | Crossref Full Text | Google Scholar

Shaffer, M., Borton, M. A., McGivern, B. B., Zayed, A. A., La Rosa, S. L., Solden, L. M., et al. (2020). DRAM for distilling microbial metabolism to automate the curation of microbiome function. Nucleic Acids Res. 48, 8883–8900. doi: 10.1093/nar/gkaa621

PubMed Abstract | Crossref Full Text | Google Scholar

Shenhav, L., Thompson, M., Joseph, T. A., Briscoe, L., Furman, O., Bogumil, D., et al. (2019). FEAST: fast expectation-maximization for microbial source tracking. Nat. Methods 16, 627–632. doi: 10.1038/s41592-019-0431-x

PubMed Abstract | Crossref Full Text | Google Scholar

Song, X., Li, Y., Stirling, E., Zhao, K., Wang, B., Zhu, Y., et al. (2022). Asgene DB: a curated orthology arsenic metabolism gene database and computational tool for metagenome annotation. NAR Genom. Bioinform. 4:lqac080. doi: 10.1093/nargab/lqac080

PubMed Abstract | Crossref Full Text | Google Scholar

Soukup, D., Buck, B., Goossens, D., Ulery, A., McLaurin, B. T., Baron, D., et al. (2012). Arsenic concentrations in dust emissions from wind erosion and off-road vehicles in the Nellis dunes recreational area, Nevada, USA. Aeolian Res. 5, 77–89. doi: 10.1016/j.aeolia.2011.11.001

Crossref Full Text | Google Scholar

Suh, S. S., Park, M., Hwang, J., Kil, E. J., Jung, S. W., Lee, S., et al. (2015). Seasonal dynamics of marine microbial community in the South Sea of Korea. PLoS One 10:e0131633. doi: 10.1371/journal.pone.0131633

PubMed Abstract | Crossref Full Text | Google Scholar

Swift, M., Andren, O., Brussaard, L., Briones, M., Couteaux, M. M., Ekschmitt, K., et al. (1998). Global change, soil biodiversity, and nitrogen cycling in terrestrial ecosystems: three case studies. Glob. Change Biol. 4, 729–743. doi: 10.1046/j.1365-2486.1998.00207.x

Crossref Full Text | Google Scholar

Tuli, R., Chakrabarty, D., Trivedi, P. K., and Tripathi, R. D. (2010). Recent advances in arsenic accumulation and metabolism in rice. Mol. Breed. 26, 307–323. doi: 10.1007/s11032-010-9412-6

Crossref Full Text | Google Scholar

Uetake, J., Tobo, Y., Uji, Y., Hill, T. C., DeMott, P. J., Kreidenweis, S. M., et al. (2019). Seasonal changes of airborne bacterial communities over Tokyo and influence of local meteorology. Front. Microbiol. 10:1572. doi: 10.3389/fmicb.2019.01572

PubMed Abstract | Crossref Full Text | Google Scholar

Uritskiy, G. V., DiRuggiero, J., and Taylor, J. (2018). MetaWRAP-a flexible pipeline for genome-resolved metagenomic data analysis. Microbiome 6:158. doi: 10.1186/s40168-018-0541-1

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, Y., Liu, Y., Wang, J., Luo, T., Zhang, R., Sun, J., et al. (2019). Seasonal dynamics of bacterial communities in the surface seawater around subtropical Xiamen Island, China, as determined by 16S rRNA gene profiling. Mar. Pollut. Bull. 142, 135–144. doi: 10.1016/j.marpolbul.2019.03.035

PubMed Abstract | Crossref Full Text | Google Scholar

Wei, M., Liu, H., Chen, J., Xu, C., Li, J., Xu, P., et al. (2020). Effects of aerosol pollution on PM2.5-associated bacteria in typical inland and coastal cities of northern China during the winter heating season. Environ. Pollut. 262:114188. doi: 10.1016/j.envpol.2020.114188

PubMed Abstract | Crossref Full Text | Google Scholar

Xie, J., Niu, X. D., Xie, J. J., He, K. Q., Shi, M. D., Yu, S. J., et al. (2021). Distribution and chemical speciation of arsenic in different sized atmospheric particulate matters. J. Environ. Sci. 108, 1–7. doi: 10.1016/j.jes.2021.02.010

PubMed Abstract | Crossref Full Text | Google Scholar

Xu, H. J., Li, S., Su, J. Q., Nie, S. A., Gibson, V., Li, H., et al. (2014). Does urbanization shape bacterial community composition in urban park soils? A case study in 16 representative Chinese cities based on the pyrosequencing method. FEMS Microbiol. Ecol. 87, 182–192. doi: 10.1111/1574-6941.12215

PubMed Abstract | Crossref Full Text | Google Scholar

Xu, A. L., Song, Z. W., Lang, X. L., Chen, X., and Xia, Y. (2017). Seasonal variability in bacterial and fungal diversity and community composition of the near-surface atmosphere in coastal megacity. Aerobiologia 33, 555–575. doi: 10.1007/s10453-017-9489-z

Crossref Full Text | Google Scholar

Yadav, A. N., Kour, D., Kaur, T., Devi, R., Yadav, A., Dikilitas, M., et al. (2021). Biodiversity, and biotechnological contribution of beneficial soil microbiomes for nutrient cycling, plant growth improvement and nutrient uptake. Biocatal. Agric. Biotechnol. 33:102009. doi: 10.1016/j.bcab.2021.102009

Crossref Full Text | Google Scholar

Yan, Y., Ma, J. J., Liang, X. P., Yin, Y., Wu, Y. Q., Yu, R. L., et al. (2024). Occurrence and spatiotemporal distribution of arsenic biotransformation genes in urban dust. Environ. Int. 190:108823. doi: 10.1016/j.envint.2024.108823

PubMed Abstract | Crossref Full Text | Google Scholar

Yang, Z. W., Salam, N., Hua, Z. S., Liu, B. B., Han, M. X., Fang, B. Z., et al. (2017). Siccirubricoccus deserti gen. Nov., sp. nov., a proteobacterium isolated from a desert sample. Int. J. Syst. Evol. Microbiol. 67, 4862–4867. doi: 10.1099/ijsem.0.002397

PubMed Abstract | Crossref Full Text | Google Scholar

Yin, X., Wang, L., Duan, G., and Sun, G. (2011). Characterization of arsenate transformation and identification of arsenate reductase in a green alga Chlamydomonas reinhardtii. J. Environ. Sci. 23, 1186–1193. doi: 10.1016/S1001-0742(10)60492-5

Crossref Full Text | Google Scholar

Yin, R., Wang, D., Deng, H., Shi, R., and Chen, Z. (2013). Heavy metal contamination and assessment of roadside and foliar dust along the outer-ring highway of Shanghai, China. J. Environ. Qual. 42, 1724–1732. doi: 10.2134/jeq2013.04.0150

PubMed Abstract | Crossref Full Text | Google Scholar

Yin, X., Wang, L., Yu, X., Du, S., Zhang, H., and Zhang, Z. (2017). Arsenic accumulation and speciation of PM2.5 and relevant health risk assessment in Jinan, China. Pol. J. Environ. Stud. 26, 949–954. doi: 10.15244/pjoes/66714

Crossref Full Text | Google Scholar

Yoshinaga, M., and Rosen, B. P. (2014). A C. As lyase for degradation of environmental organoarsenical herbicides and animal husbandry growth promoters. Proc. Natl. Acad. Sci. U. S. A. 111, 7701–7706. doi: 10.1073/pnas.1403057111

PubMed Abstract | Crossref Full Text | Google Scholar

Zecchin, S., Crognale, S., Zaccheo, P., Fazi, S., Amalfitano, S., Casentini, B., et al. (2021). Adaptation of microbial communities to environmental arsenic and selection of arsenite-oxidizing bacteria from contaminated groundwaters. Front. Microbiol. 12:634025. doi: 10.3389/fmicb.2021.634025

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, C., Fang, Y. X., Yin, X., Lai, H., Kuang, Z., Zhang, T., et al. (2023). The majority of microorganisms in gas hydrate-bearing subseafloor sediments ferment macromolecules. Microbiome 11:37. doi: 10.1186/s40168-023-01482-5

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, S. Y., Zhao, F. J., Sun, G. X., Su, J. Q., Yang, X. R., Li, H., et al. (2015). Diversity and abundance of arsenic biotransformation genes in paddy soils from southern China. Environ. Sci. Technol. 49, 4138–4146. doi: 10.1021/acs.est.5b00028

PubMed Abstract | Crossref Full Text | Google Scholar

Zhao, F. J., Harris, E., Yan, J., Ma, J., Wu, L., Liu, W., et al. (2013). Arsenic methylation in soils and its relationship with microbial ars M abundance and diversity, and as speciation in rice. Environ. Sci. Technol. 47, 7147–7154. doi: 10.1021/es304977m

PubMed Abstract | Crossref Full Text | Google Scholar

Zhao, Y., Su, J. Q., Ye, J., Rensing, C., Tardif, S., Zhu, Y. G., et al. (2019). AsChip: a high-throughput qPCR chip for comprehensive profiling of genes linked to microbial cycling of arsenic. Environ. Sci. Technol. 53, 798–807. doi: 10.1021/acs.est.8b03798

PubMed Abstract | Crossref Full Text | Google Scholar

Zhu, Y. G., Xue, X. M., Kappler, A., Rosen, B. P., and Meharg, A. A. (2017). Linking genes to microbial biogeochemical cycling: lessons from arsenic. Environ. Sci. Technol. 51, 7326–7339. doi: 10.1021/acs.est.7b00689

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: arsenic, arsenic biotransformation genes, bacterial communities, dust, metagenomes

Citation: Yin Y, Lin Y-T, Hu G-R, Yu R-L, Sun X-H and Yan Y (2025) Distinct bacterial community structures and arsenic biotransformation gene profiles in dust. Front. Microbiol. 16:1607082. doi: 10.3389/fmicb.2025.1607082

Received: 07 April 2025; Accepted: 15 July 2025;
Published: 30 July 2025.

Edited by:

Xian-Chun Zeng, China University of Geosciences Wuhan, China

Reviewed by:

Xiaomin Li, South China Normal University, China
Andrés Esteban Marcoleta, University of Chile, Chile

Copyright © 2025 Yin, Lin, Hu, Yu, Sun and Yan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yu Yan, eXlhbkBocXUuZWR1LmNu

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.