Multikingdom interactions govern the microbiome in subterranean cultural heritage sites.

SignificanceThe conservation of historical relics against microbial biodeterioration is critical to preserving cultural heritages. One major challenge is our limited understanding of microorganisms' dispersal, colonization, and persistence on relics after excavation and opening to external environments. Here, we investigate the ecological and physiological profiles of the microbiome within and outside the Dahuting Han Dynasty Tomb with a 1,800-y history. Actinobacteria dominate the microbiome in this tomb. Via interkingdom signaling mutualism, springtails carry Actinobacteria as one possible source into the tomb from surrounding environments. Subsequently, Actinobacteria produce cellulases combined with antimicrobial substances, which helps them to colonize and thrive in the tomb via intrakingdom competition. Our findings unravel the ecology of the microbiomes colonizing historical relics and provide help for conservation practices.

Temperature thresholds drive the global distribution of soil fungal decomposers.

Unraveling the biogeographic pattern of soil fungal decomposers along temperature gradients-in smooth linearity or an abrupt jump-can help us connect the global carbon cycle to global warming. Through a standardized global field survey, we identify the existence of temperature thresholds that control the global distribution of soil fungal decomposers, leading to abrupt reductions in their proportion (i.e., the relative abundance in the fungal community) immediately after crossing particular air and soil temperature thresholds. For example, small increases over the mean annual temperature threshold of ~9°C result in abrupt reductions in their proportion, paralleling a similar temperature threshold for soil carbon content. We further find that the proportion of soil fungal decomposers is more sensitive to temperature increases under arid conditions. Given the positive correlation between the global distributions of fungal decomposers and soil heterotrophic respiration, the reported temperature-driven abrupt reductions in fungal decomposers could further suppress their driven soil decomposition processes and reduce carbon fluxes from soils to the atmosphere with implications for climate change feedback. This work not only advances the current knowledge on the global distribution of soil fungal decomposers, but also highlights that small changes in temperature around certain thresholds can lead to potential unexpected consequences in global carbon cycling under projected climate change.

Balanced stochastic versus deterministic assembly processes benefit diverse yet uneven ecosystem functions in representative agroecosystems.

Ecological assembly processes, by influencing community composition, determine ecosystem functions of microbiomes. However, debate remains on how stochastic versus deterministic assembly processes influence ecosystem functions such as carbon and nutrient cycling. Towards a better understanding, we investigated three types of agroecosystems (the upland, paddy, and flooded) that represent a gradient of stochastic versus deterministic assembly processes. Carbon and nutrient cycling multifunctionality, characterized by nine enzymes associated with soil carbon, nitrogen, phosphorous and sulfur cycling, was evaluated and then associated with microbial assembly processes and co-occurrence patterns of vital ecological groups. Our results suggest that strong deterministic processes favour microorganisms with convergent functions (as in the upland agroecosystem), while stochasticity-dominated processes lead to divergent functions (as in the flooded agroecosystem). To benefit agroecosystems services, we speculate that it is critical for a system to maintain balance between its stochastic and deterministic assembly processes (as in the paddy agroecosystem). By doing so, the system can preserve a diverse array of functional traits and also allow for particular traits to flourish. To further confirm this speculation, it is necessary to develop a systematic knowledge beyond merely characterizing general patterns towards the associations among community assembly, composition, and ecosystem functions.

Application of bioinformatics to spectral analysis: soil organic carbon structure distinguished by X-ray photoelectron spectroscopy.

Spectroscopy, a powerful tool for analyzing material structure and composition, often encounters difficulties when investigating complex systems, such as soil and water. Here, it is hypothesized that bioinformatic methods based on the definition of the operational taxonomic unit can be applied to spectral analysis and to improve the resolving power of spectroscopic approaches. To test this hypothesis, we investigated SOC structure in response to three fertilization regimes using X-ray photoelectron spectroscopy (XPS). The operational taxonomic unit in spectroscopy (OTUsp) was defined and then the Manhattan plots were performed. Our approach turned out to be successful in determining the discrimination of SOC structure, while the traditional peak fitting method of XPS spectra failed. The results were then validated by chemical extraction analysis. Spectral analysis based on OTUsp can supplement traditional spectral interpretation and enhance its usability for studying complex systems. Graphical abstract.

Characteristics and Complete Genome Analysis of Bacillus asahii OM18, a Bacterium in Relation to Soil Fertility in Alkaline Soils Under Long-Term Organic Manure Amendment.

Bacillus asahii strain OM18, a bacterium in relation to soil fertility, was isolated from alkaline soils under long-term organic manure application in the North China Plain. B. asahii species play a pivotal role in the promotion of both crop yield and soil fertility via accelerating carbon and phosphorus cycling. However, little is known about the characteristics of B. asahii and its underlying molecular mechanism involved in soil nutrient cycling as well as its potential in promoting crop growth. To this end, we report the characteristics and complete genome analysis of strain OM18, which is relevant to promoting plant growth in phosphorus-deficient alkaline soils. Our results provide a glimpse into the metabolic function of B. asahii OM18.

Two key features influencing community assembly processes at regional scale: Initial state and degree of change in environmental conditions.

Belowground microbial communities strongly influence ecosystem function such that predicting function may rely on understanding ecological processes that assemble communities. Uncertainty remains, however, in what governs the relative contributions of different ecological processes. To help fill this knowledge gap, we test the general hypothesis that both initial state and degree of change in environmental conditions govern the relative contributions of different ecological assembly processes. To do so, we leveraged regional-scale nutrient and organic matter addition experiments and used soil organic matter (SOM) as a proxy of integrated soil environmental conditions. Consistent with our hypothesis, we found that both the initial amount of SOM and the degree of change in SOM-in response to nutrient addition-influenced the relative contributions of different ecological assembly processes. These influences were most clearly observed at the regional scale, suggesting potential scale dependence. More specifically, nutrient additions homogenized bacterial community composition due to enhanced influences of homogenizing dispersal when SOM content was initially high. In contrast, nutrient additions led to divergence in community composition due to variable selection when initial SOM was low and/or when SOM increased significantly in response to nutrient additions. Our findings indicate important connections among initial conditions, degree of change in environmental variables and microbial community assembly processes that may influence ecosystem processes. These conceptual inferences highlight a need to strengthen connections between ecological theory and biogeochemical modelling.

[Effect of X-ray micro-computed tomography on the metabolic activity and diversity of soil microbial communities in two Chinese soils].

OBJECTIVE: X-ray micro-computed tomography (micro-CT) technology, as used in the in situ and nondestructive analysis of soil physical structure, provides the opportunity of associating soil physical and biological assays. Due to the high heterogeneity of the soil matrix, X-ray micro-CT scanning and soil microbial assays should be conducted on the same soil sample. This raises the question whether X-ray micro-CT influences microbial function and diversity of the sample soil to be analyzed. METHODS: To address this question, we used plate counting, microcalorimetry and pyrosequencing approaches to evaluate the effect of X-ray--at doses typically used in micro-CT--on soil microorganisms in a typical soil of North China Plain, Fluvo-aquic soil and in a typical soil of subtropical China, Ultisol soil, respectively. RESULTS: In both soils radiation decreased the number of viable soil bacteria and disturbed their thermogenic profiles. At DNA level, pyrosequencing revealed that alpha diversities of two soils biota were influenced in opposite ways, while beta diversity was not affected although the relative abundances of some guilds were changed. CONCLUSION: These findings indicate that the metabolically active aspects of soil biota are not compatible with X-ray micro-CT; while the beta molecular diversity based on pyrosequencing could be compatible.

Lysinibacillus alkaliphilus sp. nov., an extremely alkaliphilic bacterium, and emended description of genus Lysinibacillus.

A novel aerobic, alkaliphilic, Gram-staining-positive, endospore-forming bacterium, strain OMN17T, was isolated from a typical sandy loam soil under long-term OMN fertilization (half organic manure N plus half mineral N fertilizer) in northern China and was subjected to a polyphasic taxonomic study. The best growth was achieved at 30 °C and pH 8-10 in medium containing 0.5% (w/v) NaCl. The cell-wall peptidoglycan of strain OMN17T was type A4α; (l-Lys-Gly-d-Asp) and the cell-wall sugars were ribose, traces of galactose and arabinose. The only respiratory quinone found in strain OMN17T was menaquinone 7 (MK-7). The major cellular fatty acids were iso-C15 : 0 and anteiso-C15 : 0. The major polar lipids were found to be phosphatidylethanolamine, phosphatidylglycerol and diphosphatidylglycerol. Phylogenetic analysis of strain OMN17T based on the 16S rRNA gene sequence showed that the strain was most closely related to Lysinibacillus halotolerans (97.8%), Lysinibacillus sinduriensis (97.5%), Lysinibacillus chungkukjangi (97.4%) and Lysinibacillus xylanilyticus (97.0%). The DNA-DNA hybridization results indicated that this strain was distinct from other species of the genus Lysinibacillus, the degree of relatedness being 21.8 ± 0.2% with L. halotolerans, 45.6 ± 1.8% with L. sinduriensis, 33.7 ± 1.2% with L. chungkukjangi and 23.7 ± 0.7% with L. xylanilyticus. The DNA G+C content of strain OMN17T was 38.1 mol%. The phenotypic, chemotaxonomic and genetic analyses identified strain OMN17T as a novel species of the genus Lysinibacillus, for which the name Lysinibacillus alkaliphilus sp. nov. is proposed. The type strain is OMN17T ( = DSM 28019T = CCTCC AB 2014073T). An emended description of the genus Lysinibacillus is also provided.

Bacillus fengqiuensis sp. nov., isolated from a typical sandy loam soil under long-term fertilization.

A Gram-staining-positive, endospore-forming, moderately alkaliphilic bacterium, strain NPK15(T), was isolated from a typical sandy loam soil under long-term NPK fertilization in northern China and was subjected to a polyphasic taxonomic study. The diamino acid of the cell-wall peptidoglycan of strain NPK15(T) was found to be meso-diaminopimelic acid and the cell-wall sugars were xylose, glucose and traces of mannose. The only respiratory quinone found in strain NPK15(T) was menaquinone 7 (MK-7). The major cellular fatty acids were iso-C(15 : 0), anteiso-C(15 : 0), C(16 : 0) and C(16 : 1)ω6c/C(16 : 1)ω7c. The major polar lipids were diphosphatidylglycerol, phosphatidylethanolamine and phosphatidylglycerol. Phylogenetic analysis of the strain based on its 16S rRNA gene sequence showed that it was related most closely to 'Bacillus thaonhiensis' KACC 17216 (99.59%), B. songklensis KCTC 13881(T) (99.52%) and B. abyssalis CCTCC AB 2012074(T) (99.00%). DNA-DNA hybridization results indicated that the strain was distinct from other species of the genus Bacillus, the degree of relatedness being 35.4% with B. abyssalis CCTCC AB 2012074(T), 39.7% with B. songklensis KCTC 13881(T) and 51.2% with 'B. thaonhiensis' KACC 17216. The DNA G+C content of strain NPK15(T) was 45.5 mol%. Phenotypic, chemotaxonomic and molecular analyses identified strain NPK15(T) as a member of a novel species of the genus Bacillus, for which the name Bacillus fengqiuensis sp. nov. is proposed. The type strain is NPK15(T) ( = DSM 26745(T) = CCTCC AB 2013156(T)).

[Spatial shifts in methanogenic archaeal community composition in eight Chinese paddy soils].

OBJECTIVE: The methanogenic archaea plays many ecological roles in paddy ecosystems such as the regulation of CH4 emissions. However, knowledge gaps exist about the mechanisms of their spatial shift in population. METHODS: We collected 8 latitudinal paddy soil samples from 20.55 degrees N to 47.43 degrees N in China, and studied their methanogenic archaeal community compositions by PCR-DGGE fingerprinting and phylogenetic analyses. Mechanism of spatial shift in community composition was unraveled by canonical correspondence analysis and Venn diagram. RESULTS: Soil pH was one of the main environmental triggers to the composition of methanogenic archaea community in paddy soils. In addition, community composition shifted along latitudinal gradients. CONCLUSION: It is the first report on biogeography of methanogenic archaeal community in arable soils and its environmental driving factors. The information would contribute to understanding spatial shifts in the transformation of organic matters along Chinese latitudinal gradient.

Urban green space area mitigates the accumulation of heavy metals in urban soils.

Despite that the heavy metals in urban soils pose a threat to public health, the critical factors that influence their concentrations in urban soils are not well understood. In this study, we conducted a survey of surface soil samples from urban green spaces in Shanghai, to analyze the concentrations of the key heavy metals. The results showed that Zn was the most abundant metal with an average concentration of 122.99 mg kg-1, followed by Pb (32.72 mg kg-1) and Cd (0.23 mg kg-1). All concentrations were found to be below the risk screening values defined by the National Environmental Quality Standards for soils of development land in China (GB36600-2018), indicating no current risk in Shanghai. However, there was a clear accumulation of heavy metals, as the mean concentrations were significantly higher than the background values. Furthermore, we explored the relationships between key heavy metals with population density, GDP and green space area. Both Spearman correlation and Random Forest analysis indicated that per capita green space area (pGSA) and population density were the most crucial factors influencing the status of heavy metals in urban soils, unlike edaphic factors e.g. SOM content in farmland soils. Specifically, there was a significantly positive linear correlation between heavy metal concentrations and population density, with correlation coefficients ranging from 0.3 to 0.4. However, the correlation with pGSA was found to be non-linear. The nonlinear regression analysis revealed threshold values between heavy metals concentrations and pGSA (e.g Zn 22.22 m2, Pb 24.92 m2, and Cd 25.92 m2), with a sharp reduction in heavy metal concentrations below the threshold and a slow reduction above the threshold. It suggests that an increase in per capita green space area can mitigate the accumulation of heavy metals caused by growing population density, but the effect is limited after the threshold. Our findings not only provide insights into the distribution patterns of heavy metals in the urban soils at the local scale, but also contribute to the urban greening at the global scale and offer guidance for city planning in the face of increasing population densities over the coming decades.

Eco-evolutionary strategies for relieving carbon limitation under salt stress differ across microbial clades.

With the continuous expansion of saline soils under climate change, understanding the eco-evolutionary tradeoff between the microbial mitigation of carbon limitation and the maintenance of functional traits in saline soils represents a significant knowledge gap in predicting future soil health and ecological function. Through shotgun metagenomic sequencing of coastal soils along a salinity gradient, we show contrasting eco-evolutionary directions of soil bacteria and archaea that manifest in changes to genome size and the functional potential of the soil microbiome. In salt environments with high carbon requirements, bacteria exhibit reduced genome sizes associated with a depletion of metabolic genes, while archaea display larger genomes and enrichment of salt-resistance, metabolic, and carbon-acquisition genes. This suggests that bacteria conserve energy through genome streamlining when facing salt stress, while archaea invest in carbon-acquisition pathways to broaden their resource usage. These findings suggest divergent directions in eco-evolutionary adaptations to soil saline stress amongst microbial clades and serve as a foundation for understanding the response of soil microbiomes to escalating climate change.

The role of metal nanoparticles in influencing arbuscular mycorrhizal fungi effects on plant growth.

A knowledge gap still remains concerning the in situ influences of nanoparticles on plant systems, partly due to the absence of soil microorganisms. Arbuscular mycorrhizal fungi (AMF) can form a mutualistic symbiosis with the roots of over 90% of land plants. This investigation sought to reveal the responses of mycorrhizal clover (Trifolium repens) to silver nanoparticles (AgNPs) and iron oxide nanoparticles (FeONPs) along a concentration gradient of each. FeONPs at 3.2 mg/kg significantly reduced mycorrhizal clover biomass by 34% by significantly reducing the glomalin content and root nutrient acquisition of AMF. In contrast, no negative effects of AgNPs at concentrations over 0.1 mg/kg were observed; however, AgNPs at 0.01 mg/kg inhibited mycorrhizal clover growth. In response to the elevated AgNPs content, the ability of AMF to alleviate AgNPs stress (via increased growth and ecological behaviors) was enhanced, which decreased Ag content and the activities of antioxidant enzymes in plants. These results were further supported by X-ray microcomputed tomography. Our findings suggest that in soil ecosystem, the influence of nanometals on plant systems would be more complicated than expected, and more attention should be focused on plant responses in combination with those of soil microorganisms.

The impact of Ag nanoparticles on methane emission in two typical paddy soils.

Numerous applications of Ag nanoparticles (AgNPs) have increased the likelihood of their release and accumulation in agroecosystem. Thus far, few studies have evaluated the impacts of AgNPs to soil methane emissions and the microbial dynamics. In this study, microcosmic experiments were conducted to investigate the responses of methanogenic processes from two paddy soils (Cambisols and Ultisols) subjected to four AgNPs doses (0.1, 1, 10 and 50 mg/kg). The results showed that 0.1 and 1 mg/kg AgNPs had no significant effects on CH4 emissions, but 50 mg/kg AgNPs increased soil CH4 emissions in both paddy soils. The aggravation effect of AgNPs on CH4 emissions was more apparent in Ultisols compared to Cambisols paddy soils. Real-time PCR suggested that 50 mg/kg AgNPs significantly increased the ratio of methanogenic to bacterial gene for both paddy soils. Amplicon sequencing indicated that methanogenic community was clustered into a separate group after 50 mg/kg AgNPs exposure. Structural equation model illustrated that Methanosarcinales was both significantly responded to AgNPs in Cambisols and Ultisols soils; however, Methanocellales significantly responded to AgNPs only in Cambisols soils. Subsequently, uncontrolled use of AgNPs may account as an environmental risk due to the potentially increased soil CH4 emissions in paddy ecosystems.

Innovative strategy for the conservation of a millennial mausoleum from biodeterioration through artificial light management.

Artificial lights can cause critical microbial biodeterioration of heritage monuments by promoting the outbreak of phototrophic microbiomes when they are used for touristic viewing. Here, with the ultimate aim of providing innovative solutions for the conservation and visiting of such monuments, we conducted a pioneering two-year in situ manipulative experiment to evaluate the impacts of different artificial light wavelengths (i.e., blue, green and red lights compared to white light) on the phototrophic microbiome of a millennial Chinese imperial mausoleum. Our results show that artificial light can shape the ecophysiological features of the phototrophic bacteriome in this monument and reduce its potential for further biodeterioration. In general, Cyanobacteria dominated (42.0% of the total relative abundance) the phototrophic bacteriome of this cultural relic; however, they were also very sensitive to the choice of artificial light. Compared to white light, monochromatic light, especially green light, reduced Cyanobacteria abundances (18.6%) by decreasing photosynthetic pigment abundances (42.9%); decreased the abundances of heterotrophic species belonging to Proteobacteria (4.5%) and the proportion of genes (6.1%) associated with carbon (i.e., carbon fixation), nitrogen (i.e., denitrification), and sulfur (i.e., dissimilatory sulfate reduction) cycling; and further decreased organic acid (10.1-14.1%) production of the phototrophic bacteriome, which is known to be involved in biodeterioration. Taken together, our findings constitute a major advancement in understanding how light wavelengths influence the phototrophic microbiome in cultural relics, and we found that artificial lights with certain wavelengths (e.g., green light) can help long-term conservation while allowing tourism activities.

Water availability creates global thresholds in multidimensional soil biodiversity and functions.

Soils support an immense portion of Earth's biodiversity and maintain multiple ecosystem functions which are essential for human well-being. Environmental thresholds are known to govern global vegetation patterns, but it is still unknown whether they can be used to predict the distribution of soil organisms and functions across global biomes. Using a global field survey of 383 sites across contrasting climatic and vegetation conditions, here we showed that soil biodiversity and functions exhibited pervasive nonlinear patterns worldwide and are mainly governed by water availability (precipitation and potential evapotranspiration). Changes in water availability resulted in drastic shifts in soil biodiversity (bacteria, fungi, protists and invertebrates) and soil functions including plant-microbe interactions, plant productivity, soil biogeochemical cycles and soil carbon sequestration. Our findings highlight that crossing specific water availability thresholds can have critical consequences for the provision of essential ecosystem services needed to sustain our planet.

Long-term balanced fertilization decreases arbuscular mycorrhizal fungal diversity in an arable soil in North China revealed by 454 pyrosequencing.

A balanced fertilization can increase crop yields partly due to stimulated microbial activities and growths. In this study, we investigated arbuscular mycorrhizal fungi (AMF) in arable soils to determine the optimal practices for an effective fertilization. We used pyrosequencing-based approach to study AMF diversity, as well as their responses to different long-term (>20 years) fertilizations, including OM (organic manure) and mix chemical fertilizers of NP (nitrogen-phosphorus), NK (nitrogen-potassium), and NPK (nitrogen-phosphorus-potassium). Results revealed that 124,998 of 18S rRNA gene fragments were dominated by Glomeromycota with 59,611 sequences, generating 70 operational taxonomic units (OTUs), of which the three largest families were Glomeraceae, Gigasporaceae and Acaulosporaceae. In Control and NK plots, AMF diversity and richness significantly decreased under long-term P fertilizations, such as NP, NPK, and OM. Concomitantly, the AMF community structure shifted. Supported by canonical correspondence analysis, we hereby propose that long-term balanced fertilization, especially P fertilizer with additional N fertilizer, helps the build-up of soil nutrients. Consequently, some AMF community constituents are sacrificed, propelled by the self-regulation of plant-AMF-microbes system, resulting in an agro-ecosystem with a better sustainability. This knowledge would be valuable toward better understandings of AMF community in agro-ecosystem, and long-term ecosystem benefits of the balanced fertilization.

Beyond the snapshot: identification of the timeless, enduring indicator microbiome informing soil fertility and crop production in alkaline soils.

BACKGROUND: Microorganisms are known to be important drivers of biogeochemical cycling in soil and hence could act as a proxy informing on soil conditions in ecosystems. Identifying microbiomes indicative for soil fertility and crop production is important for the development of the next generation of sustainable agriculture. Earlier researches based on one-time sampling have revealed various indicator microbiomes for distinct agroecosystems and agricultural practices as well as their importance in supporting sustainable productivity. However, these microbiomes were based on a mere snapshot of a dynamic microbial community which is subject to significant changes over time. Currently true indicator microbiomes based on long-term, multi-annual monitoring are not available. RESULTS: Here, using samples from a continuous 20-year field study encompassing seven fertilization strategies, we identified the indicator microbiomes ecophysiologically informing on soil fertility and crop production in the main agricultural production base in China. Among a total of 29,184 phylotypes in 588 samples, we retrieved a streamlined consortium including 2% of phylotypes that were ubiquitously present in alkaline soils while contributing up to half of the whole community; many of them were associated with carbon and nutrient cycling. Furthermore, these phylotypes formed two opposite microbiomes. One indicator microbiome dominated by Bacillus asahii, characterized by specific functional traits related to organic matter decomposition, was mainly observed in organic farming and closely associated with higher soil fertility and crop production. The counter microbiome, characterized by known nitrifiers (e.g., Nitrosospira multiformis) as well as plant pathogens (e.g., Bacillus anthracis) was observed in nutrient-deficit chemical fertilizations. Both microbiomes are expected to be valuable indictors in informing crop yield and soil fertility, regulated by agricultural management. CONCLUSIONS: Our findings based on this more than 2-decade long field study demonstrate the exciting potential of employing microorganisms and maximizing their functions in future agroecosystems. Our results report a "most-wanted" or "most-unwanted" list of microbial phylotypes that are ready candidates to guide the development of sustainable agriculture in alkaline soils.

Microbial Community Composition and Activity in Saline Soils of Coastal Agro-Ecosystems.

Soil salinity is a serious problem for agriculture in coastal regions. Nevertheless, the effects of soil salinity on microbial community composition and their metabolic activities are far from clear. To improve such understanding, we studied microbial diversity, community composition, and potential metabolic activity of agricultural soils covering non-, mild-, and severe-salinity. The results showed that salinity had no significant effect on bacterial richness; however, it was the major driver of a shift in bacterial community composition and it significantly reduced microbial activity. Abundant and diverse of microbial communities were detected in the severe-salinity soils with an enriched population of salt-tolerant species. Co-occurrence network analysis revealed stronger dependencies between species associated with severe salinity soils. Results of microcalorimetric technology indicated that, after glucose amendment, there was no significant difference in microbial potential activity among soils with the three salinity levels. Although the salt prolonged the lag time of microbial communities, the activated microorganisms had a higher growth rate. In conclusion, salinity shapes soil microbial community composition and reduces microbial activity. An addition of labile organic amendments can greatly alleviate salt restrictions on microbial activity, which provides new insight for enhancing microbial ecological functions in salt-affected soils.

Organic fertilizer potentiates the transfer of typical antibiotic resistance gene among special bacterial species.

The propagation of antibiotic resistance genes (ARGs) in environments has evoked many attentions, however, how to identify their host pathogenic bacteria in situ remains a great challenge. Here we explored the bacterial host distribution and dissemination of a typical ARG, sul1 gene, in agricultural soils through the simultaneous detection of sul1 and its host 16S rRNA gene by emulsion paired isolation and concatenation PCR (epicPCR). Compared to chemical fertilizer, organic fertilizer (chicken manure) led to a higher prevalence of sul1 gene in the soil, and dominant bacterial hosts of sul1 gene were classified into Proteobacteria and Bacteroidetes phyla. Additionally, significant higher diversity of antibiotic resistance bacteria (ARB), higher rate of horizontal gene transfer (HGT), higher rate of mobile genetic elements (MGE) and higher proportion of pathogens were all observed in the treatment of organic fertilizer. This study alerts potential health risks of manure applications in agricultural soils.