Patterns of microbial arsenic detoxification genes in low-arsenic continental paddy soils

Microbes mediate the arsenic detoxification in paddy soils, determining the fate of arsenic in soils and its availability to rice plants, yet little is known about the structures and abundances of functional genes as well as the driving forces in low-arsenic paddy fields. To depict the arsenic detoxification functional gene patterns, 429 soil samples were collected from 39 paddy fields across four climatic zones in China, with the arsenic contents ranged from 9.76 to 19.74 mg kg^-1. GeoChip, a microarray-based metagenomic technique, was used to analyze the functional genes involved in arsenic detoxification. A total of three arsenic detoxification gene families were detected, aoxB, arxA (arsenite oxidase), and arsM (methyltransferase). Both the diversity and abundance of functional genes varied significantly among sampling sites (p < 0.05) and decreased along the arsenic gradient. Arsenic detoxification genes were carried by bacteria, archaea, and eukaryotes. Redundancy analysis showed that soil samples were grouped according to both climatic zones they located in and arsenic gradients at the continental scale. Soil pH, average annual temperature (AAT), arsenic, annual average precipitation (AAP), and CEC were the most important factors in shaping the functional structure. Structural equation modeling showed that AAT (r = 0.21), pH (r = -0.20), and arsenic contents (r = -0.11) directly affected the arsenic detoxification gene abundances. These findings provide an overall picture of microbial communities involved in arsenic detoxification in paddy soils and reveal the importance of climatic factors in shaping functional genes across a large spatial scale.