Abstract

Watersheds are important for supplying fresh water, the quality of which depends on complex interplay involving physical, chemical and biological processes. As water percolates through the soil and underlying weathering rock en route to the river corridor, microorganisms mediate key geochemical transformations, yet the distribution and functional capacities of subsurface microbial communities remain little understood. We have studied metabolic capacities of microbial communities along a meadow to floodplain hillslope transect within the East-River watershed, Colorado, using genome resolved metagenomics and carbon and hydrogen stable isotopes. Very limited strain/species overlap was found at different depths below the ground surface and at different distances along the hillslope, possibly due to restricted hydraulic connectivity after early stages of snowmelt. Functions such as carbon fixation and selenate reduction were prevalent at multiple sites, although the lineages of organisms responsible tend to be location-specific. Based on its abundance, sulfur is significantly more important for microbial metabolism at the floodplain compared to on the hillslope. Nitrification and methylamine oxidation are likely only occurring within the floodplain, with nitrification capacity in shallow soil, and methylamine oxidation in deeper unsaturated sediment. Biogenic methane was detected in deep surface samples, but methanogenic organisms were not identified. Originality-Significance StatementIn a previous study within a hillslope to riparian zone transect of a sub-alpine watershed, the community structure was explored using ribosomal protein S3 genes, and the metabolic potential was hypothesized based on the presence of metabolism related genes. However, tying specific strains and species to metabolic functioning was not discussed as resolved genomes were not available. In the current study, we use genome-resolved metagenomics along with carbon and hydrogen stable isotopes to explore the spatial distribution of biogeochemical processes. By linking taxonomy and function, using multiple functional genes indicative of full metabolic pathways, we detect heterogeneity in the distribution of metabolic potential and the organisms involved with depth and landscape position. Thus, we infer how microbiome genomic variation impacts biogeochemical cycling across the watershed. We found very limited strain/species overlap at different depths below the surface and along the hillslope, possibly due to the restricted site to site hydraulic connectivity, and show that communities are largely distinct in their metabolic capacities. Both proximity to the river and the underlying Mancos shale apparently control species distribution and metabolic potential. Functions such as carbon fixation and selenate reduction were prevalent at multiple sites, although the lineages of organisms responsible tend to be location-specific. Arsenate detoxification was found to be prevalent in the riparian zone whereas selenate reduction was detected within weathered Mancos shale. We conclude that important ecosystem functions are strongly associated with the riparian zone, some of which may have crucial implications as to water quality and human health.