Antibiotic resistance genes (ARGs) and metal(loid) resistance genes (MRGs) coexist in organic fertilized agroecosystems based on their correlations in abundance, yet evidence for the genetic linkage of ARG-MRGs co-selected by organic fertilization remains elusive. Here, an analysis of 511 global agricultural soil metagenomes reveals that organic fertilization correlates with a threefold increase in the number of diverse types of ARG-MRG-carrying contigs (AMCCs) in the microbiome (63 types) compared to non-organic fertilized soils (22 types). Metatranscriptomic data indicates increased expression of AMCCs under higher arsenic stress, with co-regulation of the ARG-MRG pairs. Organic fertilization heightens the coexistence of ARG-MRG in genomic elements through impacting soil properties and ARG and MRG abundances. Accordingly, a comprehensive global map was constructed to delineate the distribution of coexistent ARG-MRGs with virulence factors and mobile genes in metagenome-assembled genomes from agricultural lands. The map unveils a heightened relative abundance and potential pathogenicity risks (range of 4-6) for the spread of coexistent ARG-MRGs in Central North America, Eastern Europe, Western Asia, and Northeast China compared to other regions, which acquire a risk range of 1-3. Our findings highlight that organic fertilization co-selects genetically linked ARGs and MRGs in the global soil microbiome, and underscore the need to mitigate the spread of these co-resistant genes to safeguard public health.

Arsenic (As) metabolism pathways and their coupling to nitrogen (N) and carbon (C) cycling contribute to elemental biogeochemical cycling. However, how whole-microbial communities respond to As stress and which taxa are the predominant As-transforming bacteria or archaea in situ remains unclear. Hence, by constructing and applying ROCker profiles to precisely detect and quantify As oxidation (aioA, arxA) and reduction (arrA, arsC1, arsC2) genes in short-read metagenomic and metatranscriptomic datasets, we investigated the dominant microbial communities involved in arsenite (As(III)) oxidation and arsenate (As(V)) reduction and revealed their potential pathways for coupling As with N and C in situ in rice paddies. Five ROCker models were constructed to quantify the abundance and transcriptional activity of short-read sequences encoding As oxidation (aioA and arxA) and reduction (arrA, arsC1, arsC2) genes in paddy soils. Our results revealed that the sub-communities carrying the aioA and arsC2 genes were predominantly responsible for As(III) oxidation and As(V) reduction, respectively. Moreover, a newly identified As(III) oxidation gene, arxA, was detected in genomes assigned to various phyla and showed significantly increased transcriptional activity with increasing soil pH, indicating its important role in As(III) oxidation in alkaline soils. The significant correlation of the transcriptional activities of aioA with the narG and nirK denitrification genes, of arxA with the napA and nirS denitrification genes and of arrA/arsC2 with the pmoA and mcrA genes implied the coupling of As(III) oxidation with denitrification and As(V) reduction with methane oxidation. Various microbial taxa including Burkholderiales, Desulfatiglandales, and Hyphomicrobiales (formerly Rhizobiales) are involved in the coupling of As with N and C metabolism processes. Moreover, these correlated As and N/C genes often co-occur in the same genome and exhibit greater transcriptional activity in paddy soils with As contamination than in those without contamination. Our results revealed the comprehensive detection and typing of short-read sequences associated with As oxidation and reduction genes via custom-built ROCker models, and shed light on the various microbial taxa involved in the coupling of As and N and C metabolism in situ in paddy soils. The contribution of the arxA sub-communities to the coupling of As(III) oxidation with nitrate reduction and the arsC sub-communities to the coupling of As(V) reduction with methane oxidation expands our knowledge of the interrelationships among As, N, and C cycling in paddy soils.

Microbial arsenic (As) methylation in paddy soil produces mainly dimethylarsenate (DMA), which can cause physiological straighthead disease in rice. The disease is often highly patchy in the field, but the reasons remain unknown. We investigated within-field spatial variations in straighthead disease severity, As species in rice husks and in soil porewater, microbial composition and abundance of arsM gene encoding arsenite S-adenosylmethionine methyltransferase in two paddy fields. The spatial pattern of disease severity matched those of soil redox potential, arsM gene abundance, porewater DMA concentration, and husk DMA concentration in both fields. Structural equation modelling identified soil redox potential as the key factor affecting arsM gene abundance, consequently impacting porewater DMA and husk DMA concentrations. Core amplicon variants that correlated positively with husk DMA concentration belonged mainly to the phyla of Chloroflexi, Bacillota, Acidobacteriota, Actinobacteriota, and Myxococcota. Meta-omics analyses of soil samples from the disease and non-disease patches identified 5129 arsM gene sequences, with 71% being transcribed. The arsM-carrying hosts were diverse and dominated by anaerobic bacteria. Between 96 and 115 arsM sequences were significantly more expressed in the soil samples from the disease than from the non-disease patch, which were distributed across 18 phyla, especially Acidobacteriota, Bacteroidota, Verrucomicrobiota, Chloroflexota, Pseudomonadota, and Actinomycetota. This study demonstrates that even a small variation in soil redox potential within the anoxic range can cause a large variation in the abundance of As-methylating microorganisms, thus resulting in within-field variation in rice straighthead disease. Raising soil redox potential could be an effective way to prevent straighthead disease.

Arsenic (As) is a toxic metalloid pervasively present in the environment. Microorganisms have evolved the capacity to metabolize As, and As metabolism genes are ubiquitously present in the environment even in the absence of high concentrations of As.