Diverse Bathyarchaeotal Lineages Dominate Archaeal Communities in the Acidic Dajiuhu Peatland, Central China

Bathyarchaeota are believed to have roles in the carbon cycle in marine systems. However, the ecological knowledge of Bathyarchaeota is limited in peatland ecosystems. Here, we investigated the vertical distribution of Bathyarchaeota community structure using quantitative PCR and high-throughput sequencing technology of ribosomal 16S rRNA gene integrated with detailed chemical profiling in the Dajiuhu Peatland, central China. Eight archaeal phyla were observed in peat samples, which mainly composed of Bathyarchaeota with a mean relative abundance about 88%, followed by Thaumarchaeota (9%). Bathyarchaeota were further split into 17 subgroups, and some subgroups showed habitat specificity to peat horizons with distinct lithological and physicochemical properties, for example, Bathy-6 and Bathy-15 had preference for the acrotelm, Bathy-5b, Bathy-16, and Bathy-19 were enriched in the catotelm, Bathy-5a, Bathy-8, and Bathy-11 were specific for the clay horizon. This spatial distribution pattern of archaeal communities along peat profile was mainly influenced by water content as indicated by RDA ordination and permutational MANOVA, whereas organic matter content exclusively affected Bathyarchaeota distribution along the peat profile significantly. The abundance of archaeal 16S rRNA genes ranged from 105 to 107 copies per gram dry sediment, and the highest archaeal biomass was observed in the periodically oxic mesotelm horizon with more dynamic archaeal interaction relationship as indicated by the network analysis. Bathyarchaeota dominated the archaeal interaction network with 82% nodes, 96% edges, and 71% keystone species. Our results provide an overview of the archaeal population, community structure, and relationship with environmental factors that affect the vertical distribution of archaeal communities and emphasize the ecology of bathyarchaeotal lineages in terrestrial peatland ecosystems.