Distribution and variation of MGEs and ARGs
The concentration of MGEs and ARGs in summer and in winter was quantified via qPCR (Fig. 1). In the air samples, intl1 was the main MGE and belongs to the integron. The intl1 concentration (6.9 ± 0.3 lg copies/m3) in the WA samples was the highest (Fig. S3; Table S3). In the WA samples, the high abundance of intl1 might enhance the transfer of airborne ARGs. The abundance of airborne ARG and MGE subtypes (except qnrA and intl1) in winter was obviously higher than that in summer (p < 0.01; Fig. 1a). The average concentration of the ARG subtypes (except strA and ermB) in the fecal samples taken in summer is higher than those taken in winter (Fig. 1b).
Among all the resistance subtypes, tetM is the dominant ARG in the chicken farm. The tetM concentration was 7.9 ± 1.3 lg copies/m3 in air and 11.8 ± 0.7 lg copies/g in the feces (Fig. S3). In the fecal samples, the tet genes concentration in a previous study (~ 9 lg copies/g [16]) is lower than that in this study. And the abundance of airborne tet genes is much higher than the reported values (~ 2 lg copies/m3 [17]), indicating the high health risk of the chicken farm’s indoor air. In addition, this study conducted a correlation analysis between the individual ARGs and the total concentration of all the ARGs in air during winter and summer. tetM presented the strongest correlation with the total concentration of all the ARG subtypes (Fig. S4, r > 0.95, p < 0.01). Therefore, tetM could be a potential indicator in the regular monitoring of airborne ARGs in the chicken farm.
Resistance mechanisms of ARGs
The specificity of antibiotics to bacteria is crucial for its usage in clinic for curing and preventing bacterial infections. However, the wide spread of ARB is raising a serious public health concern. Four main resistance mechanisms are responsible for the antibiotics in bacteria (Fig. 2a). The enzymatic function of acquired genetic traits that could cause antibiotic inactivation is a common mechanism in various pathogenic bacteria (such as Vibrio cholera) [18]. In this study, the percentage of this mechanism in all the samples is greater than 13.0% (Figs. 2b–2e, Table S4). The reduction permeability, an important defense mechanism, can prevent the bacterial cytoplasm of several antibiotics [19]. Many antibiotics have a high specific affinity for bacterial target protein, so they can inhibit indigenous cellular functions [20]. The antibiotic targets include (i) target alteration, (ii) target replacement, and (iii) target protection (Fig. 2a). Target alteration and replacement change the composition of the target molecule [21]. Target protection only inhibits the physical association between the target protein and antibiotics, which does change the nature of target protein [20]. According to the literature, 13 distinct tetracycline ribosomal protection protein classes mediate target protection [22], of which tetM and tetO are the best characterized [23, 24]. In Figs. 2b–2e, the percentage of antibiotic target protection is the highest (> 31.0%) in resistance mechanisms in all the samples. The ARG subtypes in this mechanism include tetM, tetO, and qnrA. The active efflux of antibiotics is encoded by the genes that are commonly present in the genome of bacteria. Some efflux proteins could recognize and pump out several antibiotics, such as tetracyclines, aminoglycosides, β-lactams, nalidixic acid, and macrolides (Das et al., 2020). In this study, tetG belongs to the antibiotic efflux, and its percentage is above 12.4% in all the samples.
Distribution of bacterial community
The bacterial communities in both samples were determined. The Shannon and Chao1 indexes were respectively used to evaluate the diversity and richness of bacteria [25] (Figs. 3a & 3b). In contrast to those in feces, the bacterial communities in air could be easily mixed together from different areas due to ventilation systems [26]. Meanwhile, the bacteria richness of the SA samples is higher than that of the WA samples (p < 0.01, Figs. 3a & 3b), which was also determined using the ACE indexes (Fig. S5). This result suggests that the high temperatures in summer (24.9°C on average) are conducive to the growth of bacterial communities compared with the low temperatures in winter (13.0°C on average, Table S5).
Firmicutes, Actinobacteria, Bacteroidetes, and Proteobacteria are the predominant bacteria in all the samples (Figs. S6–S7; at the phylum level), and the same result was obtained in pig manure [27]. The relative concentration of Firmicutes increased by approximately 20.9% in the feces compared with that in the air samples (Figs. S6–S7). However, Actinobacteria decreased by approximately 30.5% in the feces compared with that in the air samples (Figs. S6–S7). At the genus level, Lactobacillus (WA: 6.6%±0.3%, WF: 1.7%±1.0%, SA: 5.7%±3.0%, and SF: 12.0%±6.0%) was found in all the samples (Fig. 3c). Acinetobacter was detected in the SA and SF samples and caused various diseases (e.g., sepsis and pneumonia) [28]. Meanwhile, Ochrobactrum (8.6%±1.3%), Rhodococcus (4.2%±3.1%), Clostridium (2.7%±1.5%), and Terrisporobacter (2.1%±1.2%) were only detected in the SA samples. The concentration of pathogens (e.g., Corynebacterium) in the SA samples is higher than that in the other samples [29], Meanwhile, some strains that belong to the Rhodococcus genus could cause broncho-pneumonia [30]. Therefore, more attention should be paid to the atmosphere of chicken farms, especially in summer.
Variation in bacterial community in both air and feces
Figure 4 shows the biomarkers among the microbes in the fecal and air samples taken in summer and winter revealed by the LEfSe algorithm (p < 0.05, LDA scores < 3.5). The bacteria with no significant difference are in yellow, and the significantly different biomarkers are colored with the group (Fig. 4). Taking the fecal and air samples obtained in winter as bacterial examples (Fig. 4a), the number of airborne biomarkers (53) is lower than that of bacteria with no significant difference (249; Table S6). Figure 4b shows that the airborne bacteria contained 272 species of bacteria (at the genus level), which were significantly different from those in the fecal samples taken in summer (Fig. 4b, Table S6). Airborne bacteria may migrate from outside the chicken farm during summer. The fecal biomarkers in winter included Pseudomonadaceae and Erysipelothrix, which could be reproduced and grown at 4°C and 15°C, respectively (Fig. S8). The numbers of fecal biomarkers in summer and in winter were 10 and 27, respectively (Fig. 4c, Table S6). No obvious difference in the bacterial composition of the fecal samples taken in summer and in winter was observed. Consistent with the highest diversity of bacterial communities in the air samples indicated by the Shannon index, the number of airborne biomarkers was significantly higher in summer (271) than in winter (0; Fig. 4d, Table S6).
Influence of bacteria, environmental factors, and MGEs on composition and abundance of ARGs
In the environment, the abundance and dissemination of ARGs are associated with the bacteria [31]. Various environmental factors (for example wind speed, RH, and temperature) play significant roles in the spread and occurrence of bioaerosols [32, 33]. We checked for significant (p < 0.05) and strong (R > 0.5) positive correlations among the bacterial taxa (at the genus level, Top 20), the environmental factors, and the ARGs/MGEs (Fig. 5). Such correlations can imply the hosts of airborne ARGs [34]. The possible host of airborne ARGs/MGEs is Acinetobater (Fig. 5). In addition, the possible hosts of ARGs/MGEs is fewer in WA than in the SA samples. Eight of the top 20 genera (such as Corynebacterium and Ruminoccaceae) and the ARGs/MGEs had a significant correlation in summer (p < 0.05), indicating that the profiles of airborne ARGs might be significantly shaped by the contribution of these bacteria (Fig. 5b). However, Rothia and Kocuria were negatively correlated with the ARGs, indicating that the appearance of these bacteria would negatively affect the dissemination and occurrence of ARGs in bioaerosols (Figs. 5a & S9).
The chicken farm had no ventilation in winter, whereas the mechanical ventilation system was operated in summer. The wind speed in the chicken farm was 0.00 ± 0.0 m/s in winter (Table S5). Thus, the PM2.5 concentration during winter (472.6 µg/m3) was higher than that during summer (93.1 µg/m3; Table S5). PM2.5 and Jeotgalicoccus had a positive correlation in winter (Fig. 5a), indicating that particles might contribute to the growth of bacteria. At high RH (> 80%), the concentration of the total airborne aerobic bacteria (such as Corynebacterium) can decrease [35, 36]. The RH had a negative influence on most of the ARGs/MGEs subtypes, which may be attributed to the decrease of the possible host bacteria (Corynebacterium) (Fig. 5b). Temperature had a positive effect on some ARGs/MGEs. The increase in temperature (Table S5) facilitated the proliferation and spread of the ARGs.
ARGs transmission route in the atmospheric environment
ARGs could transfer in various bacteria through either HGT or VGT. The co-exclusion and co-occurrence relationships in MGEs, bacterial communities, and ARGs in the chicken farm were through the network to study the transmission route of ARGs in bioaerosols [37] (Fig. 5). On the basis of the number of co-occurrence edges in MGEs, bacterial communities, and ARGs, the transmission route of ARGs in bioaerosols was quantified (Table S7). In winter, most of the subtypes of ARGs and MGEs had more co-occurrence edges than the bacterial communities in the WA samples (Fig. 5a). Both Tp614 and intl1 (Plasmid and Integron, respectively) showed obvious correlations with tetG, tetO, and tetM among the visualized linkages (p > 0.05). The airborne ARGs in winter exhibited a higher possibility of HGT (PH = 71.4%) than VGT (28.6%; Eqs. 1 & 2), indicating that HGT might play a key role in airborne ARGs transfer during winter. By comparison, the subtypes of ARGs with bacterial communities had more co-occurrence edges than the MGEs in the SA samples (Fig. 5b). Therefore, the airborne ARGs in summer had a higher possibility of VGT (Pv = 57.1%) than HGT (42.9%; Eqs. 1 & 3). This finding suggests that the VGT of ARGs may play an important role in the spread of airborne ARGs in summer.
Correlations of ARGs and MGEs among all samples
The succession of ARG compositions in all the samples was determined by principal component analysis (PCA; Fig. 6). PCA illustrates that the profiles of the ARGs in the fecal samples in summer and in winter had high similarity. Nevertheless, the ARG profiles in the summer fecal samples significantly separated along PC1, indicating that the SF samples had different resistome profiles compared with the WF samples (Fig. 6). For example, the concentration of ermB in the SF samples is significantly lower than that in the WF sample (p < 0.01; Fig. 1b). The profiles of the airborne ARGs collected around the feces (from the SA samples) and feces-borne ARGs (from the SF samples) in summer show distinct distribution patterns, which are similar to the conclusion on the bacterial community (Fig. 4). Furthermore, among the atmospheric environmental samples taken in summer and in winter (SA and WA samples), the ARGs showed an obvious distinction (Fig. 6). The concentration of ARGs in the SA samples was nearly three or four orders of magnitude lower than that in WA samples (Fig. 1a). The samples were divided into three distinct clusters on the basis of the sample types: the cluster constituted by SA samples, the cluster constituted by WA samples, and the cluster constituted by the WF and SF samples.