Enrichments along gradients resolve eco-evolutionary forces on subsurface microbiomes

How a single gram of soil harbors billions of microorganisms, each with distinct genomic variants that collectively maintain coherent ecological function(s), is one of microbiology’s grand unsolved problems. A key obstacle is determining which variants contribute to individual- and community-level fitness, in which contexts, and how co-occurring ecotypes interact to divide niche space. Here, using nitrate (NO ₃ ^- )-contaminated subsurface sediment as inoculum, we have performed high throughput enrichments in laboratory media of defined carbon source compositions across ecologically relevant gradients of pH and NO ₃ ^- . Long-read metagenomics and link-community decomposition of co-occurrence networks of taxa across these enrichments has revealed context-specific functional interactions among dominant generalist and lower-abundance specialist denitrifier ecotypes that comprise 53 distinct enriched communities (EnComs) across 288 enrichments derived from a single sediment sample. We identified a single enzymatic difference of alternative NO ₃ ^- reductases (NapAB vs. NarGHI) with differing substrate affinities that provided a mechanistic explanation for competitive niche partitioning between the two dominant taxa, Neorhizobium spp. and Allorhizobium spp., along the NO ₃ ^- gradient. Genome-wide polymorphism ratios (pN/pS) revealed that selective pressures vary systematically with carbon source availability and gradients of pH and NO ₃ ^- , which helps explain the natural biodiversity and functional interactions of ecotypes within denitrifying communities in the subsurface sediment. Our findings show that controlled enrichments along ecological gradients can thus uncover eco-evolutionary forces of selection, drift, and diversification that sculpt the biodiversity of microbial populations in the natural environment.