Community structure explains antibiotic resistance gene dynamics over a temperature gradient in soil

Soils are reservoirs of antibiotic resistance genes (ARGs), but environmental dynamics of ARGs are largely unknown. Long-term disturbances offer opportunities to examine microbiome responses at scales relevant for both ecological and evolutionary processes and can be insightful for studying ARGs. We examined ARGs in soils overlying the underground coal seam fire in Centralia, PA, which has been burning since 1962. As the fire progresses, previously hot soils can recover to ambient temperatures, which creates a gradient of fire impact. We examined metagenomes from surface soils along this gradient to examine ARGs using a gene-targeted assembler. We targeted 35 clinically relevant ARGs and two horizontal gene transfer-related genes (intI and repA). We detected 17 ARGs in Centralia: AAC6-Ia, adeB, bla_A, bla_B, bla_C, cmlA, dfra12, intI, sul2, tetA, tetW, tetX, tolC, vanA, vanH, vanX and vanZ. The diversity and abundance of bla_A, bla_B, dfra12 and tolC decreased with soil temperature, and changes in ARGs were largely explained by changes in community structure. We observed sequence-specific biogeography along the temperature gradient and observed compositional shifts in bla_A, dfra12 and intI. These results suggest that increased temperatures can reduce soil ARGs but that this is largely due to a concomitant reduction in community-level diversity.