Metabolic niche differentiation and napA evolution stabilize partial denitrification in wastewater ecosystems

Although partial denitrification (PD) is increasingly applied as a nitrite-supplying strategy for anammox-based nitrogen removal, the ecological distribution, metabolic specialization, and genomic determinants of stable nitrite accumulation remain poorly understood at the ecosystem scale. Here, we reconstructed 516 high-quality metagenome-assembled genomes (MAGs) using high-depth metagenomic sequencing of 107 wastewater treatment plants and classified denitrifiers according to their nitrite production or consumption capacities. Of these genomes, 23% (120 MAGs) were classified as partial denitrifiers, 41% (211 MAGs) as complete denitrifiers, and 36% (185 MAGs) as nitrite-reducing denitrifiers, revealing pronounced functional partitioning rather than dominance by complete denitrification pathways. Comparative genomics showed that partial denitrifiers possess metabolic architectures favoring rapid carbon oxidation and NADH generation while exhibiting constrained NADPH production and biosynthetic investment, thereby promoting nitrate-to-nitrite conversion but limiting subsequent nitrite reduction. Nitrite accumulation does not result from incomplete denitrification pathways but from metabolic niche differentiation. These metabolic trade-offs were further associated with the evolutionary divergence of the periplasmic nitrate reductase gene, napA, which displayed distinct sequence characteristics and genomic contexts between partial and complete denitrifiers. Integration of carbohydrate-active enzyme repertoires further revealed metabolic complementarity between partial denitrifiers and anammox bacteria, supporting efficient carbon handoff without direct substrate competition. From an engineering perspective, operating conditions that impose moderate electron limitation, such as low or fluctuating C/N ratios and intermittent carbon feeding, may selectively enrich partial denitrifiers and enhance a stable nitrite supply for PD-anammox systems. Together, these findings identify PD as a predictable ecological state shaped by genome-encoded metabolic specialization and provide a mechanistic basis for designing robust, low-carbon nitrogen-removal processes.