Soil Microbial Community Structure and Its Contribution to Carbon Cycling in the Yalu River Estuary Wetland

Abstract

Wetland microbial communities play a vital role in ecosystem functioning, particularly in the intricate processes of carbon cycling. This study employed metagenomic sequencing to investigate the diversity, composition, structural differences, carbon cycling functional gene, and microbial species of soil microbial communities in five distinct soil types of the Yalu River estuary wetland, including shoal soil, bog soil, paddy soil, meadow soil, and brown forest soil. We further explored the influence of environmental factors on both the microbial community structure and carbon cycling functional genes. Our results revealed a bacterial-dominated soil microbial community, constituting about 97.6%. Archaea and fungi represented relatively minor fractions, at 1.9% and 0.4%, respectively. While no significant differences were observed in Chao1 indices between bacterial and fungal communities, the Shannon index revealed notable differences. Both Chao1 and Shannon indices exhibited significant variations within the archaeal communities. The dominant bacterial phyla were Proteobacteria, Actinobacteria, Acidobacteria, Bacteroidetes, and Nitrospirae. Thaumarchaeota, Crenarchaeota, and Euryarchaeota formed the major archaeal phyla, while Ascomycota, Mucoromycota, and Basidiomycota were the dominant fungal phyla. Non-metric multidimensional scaling (NMDS) analysis based on Bray-Curtis distance revealed notable differences in the bacterial, archaeal, and fungal community structures across the samples. Redundancy analysis (RDA) identified key environmental factors for the major phyla. Soil pH, soil organic carbon (SOC), electrical conductivity (EC), and total phosphorus (TP) were the main influencing factors for bacteria, while soil TP, EC, total sulfur (TS), and SOC were the primary drivers for archaeal phyla. Soil total nitrogen (TN) and EC were the main influencing factors for fungal phyla. Analysis of key carbon cycling pathway genes utilizing the Kyoto Encyclopedia of Genes and Genomes (KEGG) database and clustering heatmap revealed some variations in functional gene composition across different soil types. Mantel test indicated that pH, TN, and SOC were the primary environmental factors influencing microbial functional genes associated with soil carbon cycling. Stratified bar chart analysis further demonstrated that the major contributors to carbon cycling originated from corresponding dominatnt phyla and genera of Proteobacteria, Thaumarchaeota, Actinomycetota, Euryarchaeota, and Bacteroidota. The species and relative abundance of microorganisms associated with carbon cycling pathways varied among the samples. These findings provide a crucial reference for informing the conservation and sustainable management of wetland ecosystems in the Yalu River estuary.