Spatial and temporal variations of the community structure and abundance of Candidatus Methanoperedens nitroreducens-like archaea in paddy soils

Candidatus Methanoperedens nitroreducens (M. nitroreducens)-like archaea, catalyzing anaerobic methane oxidation coupled with nitrate reduction, could play an important role in reducing methane emissions from paddy ecosystems. Currently, the spatial and temporal variations of the communities of M. nitroreducens-like archaea and the factors regulating their distribution in paddy soils are poorly known. Here, we examined the diversity, community composition and abundance of these anaerobic methanotrophs at four representative depths (0–10, 10–20, 20–30 and 30–40 cm) of paddy soils across four rice growth stages (tillering, jointing, flowering and milky stages). High-throughput sequencing of M. nitroreducens-like archaeal mcrA genes showed that their community structure varied among different depths. In contrast, the community structure of these archaea remained relatively stable across rice growth stages. Further, both soil depth and rice growth stage had a significant impact on the abundance of mcrA genes. The gene abundance varied from 7.6 × 105 to 5.0 × 106 copies g−1 dry soil with higher values in 20–30 cm soils, and the abundance was found to be decreased from tillering to milky stages. The soil organic carbon content and ammonium content were the key factors affecting the community structure of these archaea. Our results indicated that both soil depth and crop growth stage greatly affected the abundance of M. nitroreducens-like archaea, while their community structure seemed to be more responsive to soil depth. Taken together, both space and time scales should be considered for better understanding the role of M. nitroreducens-like archaea reducing methane emissions from paddy ecosystems.