Keystone Species in Anammox Bacterial Communities: Identification, Ecological Impact, and Biotechnological Application

Abstract

This article synthesizes current research on keystone species within anaerobic ammonium oxidation (anammox) bacterial communities, which are crucial for nitrogen removal in both natural and engineered ecosystems. We explore the foundational ecology of these key taxa, from 'Candidatus Scalindua' in marine sediments to 'Candidatus Brocadia' in wastewater treatment systems. Advanced methodological frameworks for identifying keystone species, including top-down network analyses, are detailed alongside their application in optimizing bioreactor performance and community stability. The review further examines troubleshooting for community management and validates findings through cross-ecosystem comparative studies and global pattern analysis. This comprehensive overview provides researchers and biotechnology professionals with actionable insights for harnessing anammox keystone species to enhance nitrogen removal processes in environmental and biomedical applications.

Unveiling the Key Players: An Ecological Introduction to Anammox Keystone Species

The term keystone species, first introduced by Robert Paine in 1969, originally described organisms with disproportionately large ecological impacts relative to their abundance [1]. This concept has since been transferred from macroecology to microecology, where keystone taxa are now recognized as crucial "ecosystem engineers" that drive microbiome structure and functioning [2]. In microbial communities, these taxa exert considerable influence on community assembly, stability, and metabolic functions, with their removal computationally predicted to cause drastic shifts in microbiome structure and performance [3] [2]. Unlike foundational species that contribute significantly to biomass, keystone taxa often operate through specific functional traits rather than numerical dominance, making their identification and characterization both challenging and essential for understanding ecosystem dynamics.

In anaerobic ammonium oxidation (anammox) systems, the identification and understanding of keystone taxa have become particularly valuable for optimizing biological wastewater treatment processes. Anammox consortia represent complex microbial communities where anammox bacteria coexist with ammonia-oxidizing bacteria (AOB), nitrite-oxidizing bacteria (NOB), Chloroflexi bacteria (CFX), and heterotrophic denitrifying bacteria (HDB) in a tightly integrated ecosystem [4]. Within these communities, keystone taxa regulate microbial assemblage patterns and functional traits across different microbial aggregates, guiding community assembly through niche differentiation and environmental filtering [3]. This review explores the conceptual evolution of keystone taxa identification and its critical application in anammox systems for enhanced wastewater treatment performance and stability.

Defining and Identifying Keystone Taxa

Conceptual Definitions and Metrics

Keystone taxa in microbial ecosystems are primarily defined through their community importance, which can be quantified through two distinct approaches [1]. The presence-impact relationship measures how the complete removal or addition of a taxon affects ecosystem traits, while the abundance-impact relationship evaluates how changes in a taxon's abundance influence these traits [1]. In practice, presence-impact measurements align more closely with microbial manipulation techniques, such as targeted removal via antibiotics or addition through bioaugmentation.

A novel top-down identification framework has recently been developed that detects keystones by their total influence on other taxa without requiring detailed reconstruction of interspecific interaction networks [1]. This method uses Empirical Presence-abundance Interrelation (EPI) metrics derived from cross-sectional data, including distance-based measures (D1, D2) and modularity-based approaches (Q) that quantify how strongly a taxon's presence correlates with community-wide abundance profiles [1]. This represents a significant advancement beyond traditional network-based approaches that often rely on pairwise correlation analyses and assume specific functional forms for ecological interactions.

Methodological Approaches for Identification

Table 1: Methodological Approaches for Identifying Keystone Taxa

Method	Underlying Principle	Applications in Anammox Systems	Limitations
Co-occurrence Network Analysis	Identifies taxa with high centrality in microbial correlation networks	Revealed Candidatus Jettenia and Candidatus Kuenenia as keystones in granular sludge systems [3]	Correlation does not imply causation; sensitive to sequencing depth
Top-Down EPI Framework	Quantifies a taxon's influence on community abundance profiles without network reconstruction	Applied to identify keystone candidates in nitrogen-removing bioreactors [1]	Requires sufficient sample size; may detect strongly correlated而非 causative taxa
Metagenome-Assembled Genomes (MAGs)	Recovers genomes from metagenomic data to infer metabolic potential	Revealed Sulfurovum possesses oxidation resistance and electron transport capabilities [2]	Computational complexity; potential for chimeric assemblies
Perturbation Experiments	Directly tests community response to taxon removal or addition	Demonstrated Thauera and Afipia as key denitrifying partners for anammox bacteria [5]	Technically challenging for uncultivated taxa; may disrupt multiple interactions

The 3C-strategy (co-occurrence network analysis, comparative genomics, and co-culture of captured keystone taxa) has emerged as a robust framework for characterizing keystone taxa in complex microbial communities [2]. This integrated approach combines computational identification with experimental validation, enabling researchers to move beyond correlation to causation in establishing keystone functions. In anammox systems, this methodology has revealed how keystone taxa form keystone guilds - functional groups of microorganisms that collectively exert disproportionate influence on community structure and function [2].

Keystone Taxa in Anammox Systems

Anammox Core Consortia and Microbial Interactions

Anammox systems function through a complex microbial community where anammox bacteria serve as the primary nitrogen-removal engines, supported by interacting microbial partners. This anammox core consortium primarily contains anammox bacteria, ammonia-oxidizing bacteria (AOB), nitrite-oxidizing bacteria (NOB), Chloroflexi bacteria (CFX), and heterotrophic denitrifying bacteria (HDB) [4]. Within this consortium, keystone taxa regulate microbial assemblage patterns through both cooperative and competitive interactions that ultimately determine system performance.

Table 2: Key Functional Groups and Their Roles in Anammox Systems

Functional Group	Representative Genera	Primary Ecological Function	Keystone Potential
Anammox Bacteria	Candidatus Brocadia, Candidatus Kuenenia, Candidatus Jettenia [5]	Convert NH₄⁺ and NO₂⁻ to N₂ gas	High - core metabolic engineers
Sulfur-Oxidizing Bacteria	Sulfurovum, Sulfurimonas [2]	Couple sulfur oxidation to nitrate reduction	Context-dependent - keystone under specific conditions
Denitrifying Bacteria	Thauera, Afipia [5]	Reduce nitrate to nitrogen gas using organic carbon	Moderate - functional partners in keystone guilds
Nitrate-Reducing Bacteria	Thioalkalispira [2]	Reduce nitrate while oxidizing sulfur compounds	Moderate - niche-specific keystones

Research has demonstrated that the abundance of anammox bacteria does not always directly correlate with reactor performance, highlighting the importance of microbial interactions mediated by keystone taxa [4]. These interactions include cross-feeding of metabolites, where partner bacteria provide essential growth factors, vitamins, and cofactors to anammox bacteria, while anammox bacteria supply metabolic intermediates to their partners [5]. For instance, metagenome-assembled genomes-based ecological modeling has revealed that dominant denitrifiers like Thauera can provide amino acids, cofactors, and vitamins to anammox bacteria, creating mutualistic relationships that enhance system performance [5].

Documented Keystone Taxa in Anammox Systems

Several specific taxa have been identified as keystones in different anammox configurations. In simultaneous anammox and denitrification (SAD) systems, Candidatus Jettenia demonstrates keystone properties through its unique ability to utilize volatile fatty acids, outcompeting other anammox microorganisms and denitrifiers when organic carbon is present [3]. Similarly, Candidatus Kuenenia has been shown to enhance microbial interactions through metabolic cooperation, with higher abundances leading to increased expression of genes involved in anammox transformation (hzs, nir) and carbon metabolism (fdh, glgA/B/C, acs) [4].

In sulfuretum-anammox coupled systems, sulfur-oxidizing bacteria like Sulfurovum and Sulfurimonas exhibit keystone characteristics through role transitions in response to environmental conditions [2]. These taxa form keystone guilds with other functional microorganisms: Sulfurimonas partners with Thioalkalispira in nitrate-replenished systems, while Sulfurovum mutualizes with PAH-degraders like Novosphingobium and Robiginitalea in benzo[a]pyrene-contaminated systems [2]. These contextual keystone functions demonstrate the environment-dependent nature of microbial keystone roles.

Diagram 1: Role Transitions of Keystone Taxa in Response to Environmental Changes. Keystone taxa such as Sulfurimonas and Sulfurovum transition to keystone roles and form specialized guilds in response to specific environmental factors like nitrate addition or benzo[a]pyrene (BaP) contamination, ultimately reshaping community structure and function [2].

Experimental Approaches for Keystone Taxa Characterization

Enrichment Bioreactor Protocols

Controlled bioreactor systems provide ideal platforms for studying keystone taxa dynamics in anammox communities. The following protocol outlines a standardized approach for anammox enrichment and keystone taxa identification:

• Bioreactor Setup: Establish four anaerobic bioreactors (5L working volume) as replicates, each seeded with 1.0 kg of source sediment (e.g., from eutrophic lakes or wastewater treatment systems) [5]. Maintain anoxic conditions by continuously flushing with argon gas at 0.5 LPM for 30 minutes before each experiment and cover reactors with tin foil to block light.

• Operational Parameters: Maintain temperature at 34±1°C with hydraulic retention time (HRT) of 24-48 hours depending on removal efficiency [5]. Provide continuous clockwise mixing at 60 rpm to enhance microbial contact.

• Feeding Strategy: Use synthetic wastewater containing NH₄⁺ and NO₂⁻ at a molar ratio of 1:1.32 [6]. Essential components include:

  • (NHâ‚„)â‚‚SOâ‚„ and NaNOâ‚‚ as nitrogen sources
  • NaHCO₃ (0.5 g/L) and KHCO₃ (0.5 g/L) as buffering agents
  • KHâ‚‚POâ‚„ (0.027 g/L) as phosphorus source
  • MgSO₄·7Hâ‚‚O (0.02 g/L) and CaCl₂·2Hâ‚‚O (0.136 g/L) as micronutrient sources
  • Trace element solutions I and II (1-1.2 mL/L) [6]

• Monitoring Regimen: Continuously track influent and effluent NH₄⁺-N, NO₂⁻-N, and NO₃⁻-N concentrations. Calculate nitrogen removal efficiency (NRE) weekly using the formula: NRE = [(Influent N - Effluent N) / Influent N] × 100% [4].

• Microbial Community Analysis: Collect biomass samples at regular intervals (e.g., days 0, 30, 60, 90) for 16S rRNA amplicon sequencing, quantitative PCR of functional genes (hzsB, nirS, nirK), and metagenomic analysis [5].

Network Analysis and Keystone Identification Protocol

The computational identification of keystone taxa from sequencing data involves a multi-step process:

• Data Preprocessing: Process raw sequencing data using MOTHUR (version 1.45.3) or QIIME2 pipelines [4]. Filter sequences for quality, remove chimeras, and cluster into operational taxonomic units (OTUs) at 97% similarity threshold.

• Network Construction: Construct co-occurrence networks using SparCC or SPIEC-EASI algorithms to minimize false positives from compositional data [3] [2]. Calculate all pairwise correlations between microbial taxa and create adjacency matrices.

• Topological Analysis: Calculate network topology parameters including:

  • Modularity: Measures degree of compartmentalization in the network (values >0.4 indicate modular structure)
  • Betweenness Centrality: Identifies taxa that connect different modules
  • Zi-Pi Plot Analysis: Classifies nodes into categories (peripherals, connectors, module hubs, and network hubs) based within-module connectivity (Zi) and among-module connectivity (Pi) [6]

• Keystone Identification: Apply the Empirical Presence-abundance Interrelation (EPI) framework using three metrics:

  • D1: Distance-based metric comparing abundance profiles when taxon is present versus absent
  • D2: Alternative distance-based metric with different weighting
  • Q: Modularity-based approach measuring how taxon presence affects community structure [1]

• Validation: Statistically validate identified keystone taxa through permutation tests (n=1000) to determine whether their EPI values are significantly higher than expected by chance [1].

Research Reagent Solutions for Keystone Taxa Studies

Table 3: Essential Research Reagents and Their Applications in Keystone Taxa Studies

Reagent/Category	Specific Examples	Application in Keystone Taxa Research	Technical Notes
Molecular Biology Kits	16S rRNA Amplification Kits (e.g., 515F/806R primers)	Microbial community profiling in anammox systems [5]	Target V4 region for bacterial diversity
qPCR Reagents	Functional gene primers (hzsB, nirS, nirK)	Quantification of anammox and denitrifying populations [5]	hzsB is specific for anammox bacteria
Trace Element Solutions	Solution I: EDTA + FeSO₄·7H₂O; Solution II: EDTA + NaMoO₄·2H₂O + NiCl₂·6H₂O + CuSO₄·5H₂O + CoCl₂·6H₂O + ZnSO₄·7H₂O + MnCl₂·4H₂O [6]	Anammox bioreactor maintenance	Essential for anammox bacterial growth
Metagenomic Sequencing Kits	Illumina NovaSeq, PacBio HiFi	MAG construction and metabolic potential analysis [2]	Long-read technologies improve assembly
Fluorescent Labeling Systems	eGFP-labeling vectors	Tracking keystone taxa in coculture experiments [2]	Visualize microbial interactions
Culture Media Components	Volatile fatty acids (e.g., acetate, propionate)	Testing mixotrophic capabilities of anammox keystones [3]	Candidatus Jettenia can utilize VFAs

Quantitative Effects of Keystone Taxa in Anammox Systems

Performance Metrics and System Stability

Keystone taxa significantly influence nitrogen removal performance and operational stability in anammox systems. Quantitative studies demonstrate that reactors with higher abundances of key keystone taxa achieve superior performance:

Table 4: Quantitative Impacts of Keystone Taxa on Anammox System Performance

Performance Parameter	System with Keystone Taxa	System without Keystone Taxa	Reference
NH₄⁺-N Removal Efficiency	85.92-95.34%	22.96% (initial efficiency)	[5]
NO₂⁻-N Removal Efficiency	Up to 95.34%	Significant accumulation observed	[5]
Start-up Duration	6 days (with 8% anammox inoculum)	85 days (activated sludge only)	[7]
Anammox Bacterial Abundance	Increased from 5.85% to 11.43%	Decreased under stress conditions	[6]
Nitrogen Removal Rate	5.4 gN/L/d in mature granules	<1.0 gN/L/d during lag phase	[7]

The presence of appropriate keystone taxa significantly accelerates anammox system start-up. Research shows that adding just 1-8% of anammox granules to activated sludge inoculum reduces start-up duration from 85 days to 6-0 days, while simultaneously decreasing lag phase duration and cell lysis periods [7]. This demonstrates the critical role of keystone taxa in rapidly establishing functional microbial architecture.

Stress Response and Functional Resilience

Keystone taxa enhance system resilience to environmental fluctuations, particularly nitrogen loading variations commonly encountered in wastewater treatment. Under suboptimal nitrogen loading conditions (<3.68 kg/m³·d), microbial communities with intact keystone taxa employ modular collaboration to counteract loading stress, evidenced by modularity indices of 0.563 and 0.545 during inhibition and starvation phases, respectively [6]. Zi-Pi plot analyses further demonstrate significantly increased inter-module connectivity, indicating reinforced interspecies interactions that help communities resist nitrogen-loading fluctuations [6].

Metagenomic analyses reveal that keystone taxa enhance functional resilience through metabolic versatility. For instance, Sulfurovum exhibits superior oxidation resistance and electron transport capabilities, enabling it to protect anammox guild members from reactive oxygen species generated during stress conditions [2]. Experimental co-culture studies confirm that keystone taxa like Sulfurovum enhance ROS removal, cell growth, and degradation efficiency when partnered with BaP-degrading bacteria, demonstrating their protective role in contaminant-stressed systems [2].

Diagram 2: Keystone Taxa-Mediated Stress Response Mechanisms in Anammox Systems. Keystone taxa coordinate community-wide stress responses through modular collaboration, cross-feeding, and protective mechanisms that maintain system function under fluctuating nitrogen loading rates or antibiotic exposure [6] [2].

The concept of keystone taxa provides a powerful framework for understanding and optimizing anammox systems for wastewater treatment. The transition from classical definitions based on interaction strength to modern top-down identification approaches represents significant methodological progress in microbial ecology [1]. In anammox systems, keystone taxa such as Candidatus Jettenia, Candidatus Kuenenia, Sulfurovum, and Sulfurimonas play disproportionate roles in maintaining system stability, accelerating start-up times, and enhancing functional resilience under stressful conditions [4] [3] [2].

Future research should focus on developing more sophisticated manipulation strategies for keystone taxa, potentially through targeted bioaugmentation or environmental conditioning to promote their establishment and activity. The 3C-strategy (co-occurrence network analysis, comparative genomics, and co-culture) provides a robust framework for identifying and validating keystone taxa in complex microbial communities [2]. Additionally, greater attention to the contextual nature of keystone functions - where a taxon may serve as a keystone under specific environmental conditions but not others - will enhance our ability to predict and manage anammox system performance across varying operational parameters.

As wastewater treatment facilities face increasing challenges from fluctuating loads, inhibitory compounds, and stringent effluent requirements, understanding and harnessing keystone taxa offers promising pathways to more robust, efficient, and resilient nitrogen removal systems. Integrating keystone ecology with process engineering represents the next frontier in optimizing anammox technology for sustainable wastewater management.

Within the broader context of keystone species research in microbial communities, anaerobic ammonium-oxidizing (anammox) bacteria represent a critical functional group in the global nitrogen cycle. These organisms perform the anammox reaction, oxidizing ammonium with nitrite as an electron acceptor under anoxic conditions to produce dinitrogen gas [8]. Since their discovery in the 1990s, anammox bacteria have been identified as keystone species in both natural ecosystems and engineered wastewater treatment systems, where they drive substantial nitrogen loss with lower energy requirements and greenhouse gas emissions compared to conventional nitrogen-removing processes [9] [10]. Their activity can account for 24-67% of N loss in marine sediments and 20-40% in suboxic water columns, fundamentally shaping nitrogen availability in these environments [8]. This whitepaper provides an in-depth examination of four major anammox genera—Scalindua, Brocadia, Kuenenia, and Jettenia—focusing on their distinct ecological niches, physiological adaptations, and roles as keystone species in anammox bacterial communities.

The distribution of anammox genera across environmental gradients demonstrates clear niche partitioning, primarily driven by salinity, temperature, pH, and organic matter content [8] [11]. This specialization positions different genera as keystone species within their respective habitats, where they disproportionately influence nitrogen cycling dynamics.

• Salinity represents the most significant factor governing geographical distribution. Ca. Scalindua dominates saline environments, including marine ecosystems and estuaries, while Ca. Brocadia, Ca. Kuenenia, and Ca. Jettenia primarily inhabit freshwater environments [8] [11] [12]. This divergence is reflected in fundamental genomic and proteomic differences; halophilic Scalindua species share a unique set of genes absent in non-halophilic relatives and exhibit a distinct bias toward acidic amino acids in their proteomes [9].

• Temperature and pH ranges also differentiate anammox genera. While all anammox bacteria are considered mesophilic, their specific temperature optima and tolerance ranges vary, influencing their distribution across climates and bioreactor operations [9]. Similarly, pH tolerance affects community composition, with different genera exhibiting optimal activity at specific pH ranges [8].

• Organic matter content influences niche differentiation through metabolic flexibility. Some anammox bacteria, such as Ca. Anammoxoglobus propionicus, can oxidize short-chain fatty acids like propionate, potentially disguising themselves as denitrifiers and occupying niches with organic carbon present [8]. This metabolic versatility allows certain species to thrive in more complex waste streams.

Table 1: Habitat Preferences of Major Anammox Genera

Genus	Salinity Preference	Common Environments	Relative Abundance	Geographical Distribution
Scalindua	High (Marine)	Oceanic water columns, marine sediments, estuaries [8] [11] [12]	Dominant in marine systems [11]	Global oceans and coastal areas [11]
Brocadia	Low (Freshwater)	Wastewater treatment plants, freshwater sediments, terrestrial ecosystems [8] [13]	Common in engineered systems [13]	Worldwide (in engineered and natural freshwater)
Kuenenia	Low to Moderate	Wastewater treatment plants, can adapt to moderate salinity [8] [13]	Common in engineered systems	Worldwide (primarily in engineered systems)
Jettenia	Low (Freshwater)	Wastewater treatment plants, denitrifying sludge [8] [13]	Common in engineered systems	Worldwide (in engineered and natural freshwater)

The community assembly of anammox bacteria in complex environments like coastal sediments is shaped by a combination of deterministic and stochastic processes. Recent studies indicate that ecological drift predominantly shapes the overall community, while rare species are more susceptible to dispersal limitations and environmental selection [11]. In these networks, Ca. Scalindua often functions as a keystone genus, with rare species playing a crucial role in maintaining the ecological stability of the anammox bacterial community [11].

Physiological and Biochemical Diversity

Underlying the ecological niche differentiation are significant physiological and biochemical variations among the major genera, which affect their growth kinetics, substrate affinities, and metabolic pathways.

Growth Kinetics and Nutrient Uptake

Quantitative physiological characterizations have revealed genus-specific differences in intrinsic microbial growth kinetics, notably the maximum specific growth rate (μmax) and the half-saturation constant (Ks) for substrates like nitrite [8]. These kinetic parameters are critical for predicting population dynamics and competitiveness under given conditions, such as in wastewater treatment reactors where substrate concentrations fluctuate.

Table 2: Physiological Characteristics of Major Anammox Genera

Genus	Nitrite Half-Saturation Constant, Ks (mg-N/L)	Metabolic Flexibility	Salt Tolerance	Key Enzymatic Features
Scalindua	5 - 15 [8]	Primarily chemolithoautotrophic [8]	High (Obligate halophile) [8]	Abundantly expressed cytochrome cd1 NirS [14] [15]
Brocadia	< 5 [8]	Can oxidize some organic acids [8]	Low [13]	Lacks canonical NirS/NirK; suspected HAO-like enzyme reduces NO2- to NH2OH [15]
Kuenenia	~0.4 [8]	Primarily chemolithoautotrophic	Moderate (can adapt to ~30 g NaCl/L) [13]	Possesses NirS and HAO-like enzymes (HAOr, Kuste4574) for NO2- reduction [14] [15]
Jettenia	0.6 - 28 [8]	Can oxidize some organic acids [8]	Moderate (can adapt to ~27.5 g NaCl/L) [13]	Copper-containing NirK [15]

Enzymatic Machinery and Nitrite Reduction Pathways

A key biochemical difference among anammox genera lies in the enzymatic machinery for the first step of the anammox metabolism: the reduction of nitrite. Research indicates significant genus-level redundancy and diversity in nitrite-reducing enzymes, which may enhance adaptability to environmental changes [14] [15].

• Genus-Specific Nitrite Reductases: Ca. Scalindua utilizes a cytochrome cd1-type nitrite reductase (NirS) that is prominently expressed [14]. Ca. Jettenia employs a copper-containing nitrite reductase (NirK) [15]. In contrast, Ca. Brocadia species lack genes for canonical NirS or NirK, suggesting a novel nitrite reductase, potentially a hydroxylamine oxidoreductase (HAO)-like protein that reduces nitrite to hydroxylamine instead of nitric oxide [15].
• Enzymatic Redundancy in Kuenenia: The model organism Ca. Kuenenia stuttgartiensis demonstrates functional redundancy, possessing multiple enzymes capable of nitrite reduction, including NirS and HAO-like proteins (HAOr and the membrane-bound Kuste4574) [14] [15]. Proteomic studies and activity assays suggest that NirS and HAOr are the primary active nitrite reductases, allowing the bacterium to maintain metabolic flexibility under varying nitrite concentrations [14].

The following diagram illustrates the complex experimental workflow used to identify and characterize these nitrite-reducing enzymes in Ca. Kuenenia stuttgartiensis, highlighting the multidisciplinary approach required in this area of research.

Research Reagent Solutions and Methodologies

Studying the physiology and biochemistry of fastidious anammox bacteria requires specialized reagents and methodological approaches. The following table details key research solutions essential for experimental investigations in this field.

Table 3: Essential Research Reagents and Methodologies for Anammox Research

Reagent / Method	Function / Application	Specific Examples / Notes
16S rRNA Gene-Targeted qPCR	Quantification of specific anammox populations in mixed communities [8].	Specific assays developed for Ca. Brocadia, Ca. Kuenenia, Ca. Jettenia, and Ca. Scalindua [8].
Primer Sets Brod541F / Amx820R	Amplification of anammox bacterial 16S rRNA gene for diversity analysis [11].	Used for high-throughput sequencing to analyze community structure in environmental samples [11].
Membrane Inlet Mass Spectrometry (MIMS)	Highly sensitive measurement of gas production/consumption (e.g., Nâ‚‚, NO, Nâ‚‚O) in enzyme activity assays [14].	Used to track NO production from nitrite reductases in protein fractions from Ca. Kuenenia [14].
Fast Protein Liquid Chromatography (FPLC)	Separation and enrichment of active proteins from cell extracts [14].	Size exclusion and anion exchange chromatography used to isolate nitrite-reducing enzymes [14] [15].
Size Exclusion Chromatography	Separates proteins by molecular size; used to fractionate nitrite-reducing activities [15].	Active fractions from Ca. Kuenenia CSTR1 corresponded to 150–200 kDa [15].
Anion Exchange Chromatography	Separates proteins by charge; used for further purification of protein fractions [15].	Resulted in lower activity yields for nitrite reduction, suggesting complex enzyme requirements [15].
Synthetic Wastewater Media	Enrichment and continuous cultivation of anammox bacteria in bioreactors [13].	Typically contains NH₄⁺ (e.g., NH₄Cl), NO₂⁻ (e.g., NaNO₂), bicarbonate buffer, and essential minerals [13].

Key Experimental Protocols

Detailed below are foundational protocols that have enabled the physiological and ecological characterization of anammox bacteria.

• Protocol for Physiological Characterization (Kinetic Parameter Estimation): The maximum specific growth rate (μmax) and half-saturation constant (Ks) are determined using bioreactors operated under substrate-limiting conditions. The specific substrate uptake rate is measured at various substrate concentrations and fitted to the Monod equation [8]. This approach has been successfully applied to compare the kinetics of different species, such as Ca. Brocadia sinica, Ca. Jettenia caeni, and Ca. Kuenenia stuttgartiensis [8] [16].

• Protocol for Investigating Salinity Adaptation: To study long-term adaptation to salinity, continuous stirred-tank reactors (CSTRs) are inoculated with anammox granules. The salinity is gradually increased in a step-wise manner using NaCl. Nitrogen removal performance is monitored via regular measurement of NH₄⁺-N, NO₂⁻-N, and NO₃⁻-N concentrations. Microbial community shifts are tracked using 16S rRNA gene-based qPCR and high-throughput sequencing [13].

The major anammox genera Scalindua, Brocadia, Kuenenia, and Jettenia have evolved distinct genomic, physiological, and biochemical traits that enable them to function as keystone species in their respective habitats. The clear niche differentiation along salinity gradients, with Scalindua dominating marine ecosystems and the other genera prevailing in freshwater and engineered systems, is a fundamental principle governing their distribution. Furthermore, differences in growth kinetics, metabolic flexibility, and the very enzymes responsible for the core anammox reaction underscore a significant evolutionary diversification. Understanding these nuances is not merely an academic exercise; it is critical for optimizing anammox-based wastewater treatment processes and accurately modeling the global nitrogen cycle. Future research, particularly the pursuit of pure cultures and the development of genetic tools, will be essential to unravel the precise mechanistic links between genetic makeup, physiological function, and ecological success in these environmentally crucial bacteria.

Ecological Niche Differentiation Among Anammox Bacteria

Anaerobic ammonium-oxidizing (anammox) bacteria are pivotal players in the global nitrogen cycle, responsible for the conversion of bioavailable nitrogen into dinitrogen gas under anoxic conditions. These planctomycetous bacteria perform the anammox reaction, oxidizing ammonium with nitrite as the electron acceptor, and have transformed our understanding of nitrogen transformations in both natural and engineered ecosystems. Their discovery revealed a previously overlooked pathway accounting for significant nitrogen loss in marine systems—up to 50% of N2 production in some environments—while simultaneously offering revolutionary applications in wastewater treatment.

The ecological success of anammox bacteria across diverse habitats stems from extensive niche differentiation among different species and genera. This differentiation arises from variations in physiological characteristics, metabolic capabilities, and environmental tolerances, enabling distinct anammox taxa to dominate under specific conditions. Understanding these niche specializations is crucial for predicting nitrogen flux in changing environments and optimizing anammox-based technologies. Furthermore, within complex microbial communities, certain anammox species may function as keystone species, disproportionately impacting community structure and function through their metabolic activities and interactions.

This technical review synthesizes current knowledge on ecological niche differentiation among anammox bacteria, with particular emphasis on their roles as keystone species in anaerobic ammonium-oxidizing communities. We examine the physiological basis for niche partitioning, distribution patterns across environmental gradients, experimental approaches for studying anammox ecology, and implications for both natural ecosystems and engineered systems.

Physiological Basis for Niche Differentiation

Growth Kinetics and Substrate Affinity

Fundamental microbial growth kinetics, particularly the maximum specific growth rate (μmax) and half-saturation constant (Ks), create primary differentiation among anammox species. These parameters determine competitive abilities under varying substrate concentrations through the Monod model, establishing a physiological hierarchy in resource acquisition [8].

Species with lower Ks values for ammonium or nitrite possess a competitive advantage in substrate-limited environments. Current data suggest that "Candidatus Brocadia" exhibits relatively high substrate affinity, supporting its prevalence in wastewater treatment systems with fluctuating nitrogen loads. In contrast, "Candidatus Kuenenia" often dominates high-nitrogen environments, while "Candidatus Scalindua" demonstrates adaptations to the low-nitrogen conditions characteristic of marine systems [8]. These kinetic parameters directly influence population dynamics and species succession in both natural and engineered systems.

Salinity Tolerance and Adaptation

Salinity represents one of the most significant environmental filters structuring anammox communities, creating a clear distinction between freshwater and marine taxa. "Candidatus Scalindua" dominates saline environments, with all known enrichments being obligately halophilic, while other genera primarily inhabit freshwater ecosystems [8]. This phylogenetic distinction correlates with fundamental physiological adaptations to osmotic stress.

The dominance of "Candidatus Scalindua" in marine ecosystems extends from coastal sediments to oxygen minimum zones (OMZs), where it can constitute nearly 100% of the anammox community [11]. Meanwhile, "Candidatus Brocadia," "Candidatus Kuenenia," "Candidatus Jettenia," and "Candidatus Anammoxoglobus" predominantly occur in freshwater environments, though some exhibit tolerance to brackish conditions [8]. This salinity-based niche partitioning has profound implications for global nitrogen cycling, with different taxa responsible for nitrogen loss in marine versus terrestrial and freshwater systems.

Oxygen Tolerance and Relationship with Aerobic Microorganisms

Despite being anaerobic processes, anammox bacteria frequently inhabit oxic-anoxic interfaces, necessitating varying degrees of oxygen tolerance. This tolerance enables cooperation with aerobic ammonia-oxidizing archaea (AOA) and bacteria (AOB), which provide the essential nitrite for the anammox reaction [17]. The interaction represents a fascinating syntrophy where aerobic microorganisms create favorable conditions for their anaerobic counterparts.

Different anammox species exhibit distinct oxygen sensitivities and relationships with aerobic nitrifiers. Some "Candidatus Brocadia" genotypes demonstrate remarkable aerotolerance, maintaining metabolic activity and gene expression across a wide range of dissolved oxygen concentrations (0-10 mg/L) in aquifer systems [17]. This tolerance enables them to thrive in habitats with fluctuating oxygen levels, including groundwater, wetlands, and the peripheral regions of wastewater treatment biofilms.

The cooperation between anammox bacteria and aerobic ammonia oxidizers extends beyond mere proximity. Metagenomic analyses reveal coordinated expression patterns between anammox genes (hzsB) and aerobic ammonia oxidation genes (amoA), suggesting tight metabolic coupling in redox transition zones [17]. This relationship exemplifies how niche differentiation in anammox bacteria is influenced not only by their direct environmental tolerances but also by their interactions with other microbial functional groups.

Organic Matter Utilization and Metabolic Versatility

While traditionally considered obligate chemolithoautotrophs, certain anammox bacteria display metabolic flexibility regarding organic compound utilization. This versatility creates another axis for niche differentiation, with implications for their distribution along organic carbon gradients [8].

Some "Candidatus Brocadia" and "Candidatus Anammoxoglobus" strains can oxidize short-chain fatty acids like propionate and acetate, coupling this oxidation to the reduction of nitrate and/or nitrite to ammonium [8]. This metabolic capability allows them to function as facultative denitrifiers under certain conditions, expanding their potential niches beyond strictly autotrophic lifestyles. The capacity for dissimilatory nitrate reduction to ammonium (DNRA) using organic electron donors provides additional ecological flexibility, particularly in carbon-rich environments [18].

Table 1: Physiological Characteristics of Anammox Bacteria Genera

Genus	Salinity Preference	Typical Habitat	Metabolic Flexibility	Notable Characteristics
Candidatus Scalindua	Marine (obligately halophilic)	Marine sediments, oxygen minimum zones	Limited	Dominant in marine systems; phylogenetically distinct from other genera
Candidatus Brocadia	Freshwater	Wastewater treatment systems, freshwater sediments	High (can oxidize volatile fatty acids)	High substrate affinity; common in engineered systems
Candidatus Kuenenia	Freshwater	Wastewater treatment reactors	Moderate	Often dominates high-nitrogen environments
Candidatus Jettenia	Freshwater	Wastewater treatment, freshwater sediments	Moderate	Tolerant to low nitrogen loading rates
Candidatus Anammoxoglobus	Freshwater	Wastewater treatment	High (can oxidize propionate)	Propionate oxidation capability

Environmental Distribution and Niche Specialization

Marine versus Freshwater Ecosystems

The distribution of anammox bacteria across aquatic ecosystems demonstrates clear niche partitioning along salinity gradients. Marine systems, including oxygen minimum zones (OMZs), continental shelf sediments, and deep-sea environments, are overwhelmingly dominated by "Candidatus Scalindua" [11]. This genus exhibits remarkable phylogenetic diversity within marine systems, suggesting further microdiversification and niche specialization among marine populations [11].

In freshwater ecosystems, including lakes, rivers, and groundwater, the anammox community is more diverse, typically comprising "Candidatus Brocadia," "Candidatus Kuenenia," "Candidatus Jettenia," and "Candidatus Anammoxoglobus" [5]. The specific dominance patterns within freshwater systems depend on local conditions such as nitrogen concentration, organic matter content, and dissolved oxygen levels. For instance, "Candidatus Brocadia" frequently dominates eutrophic lakes and wastewater treatment systems, while "Candidatus Jettenia" appears more prevalent in oligotrophic groundwater [17].

Engineered Ecosystems

Engineered ecosystems exhibit distinctive anammox community patterns influenced by operational parameters. Wastewater treatment plants with stable, high-nitrogen loads often favor "Candidatus Kuenenia," while systems with fluctuating loads may maintain "Candidatus Brocadia" or "Candidatus Jettenia" [19] [20]. These distribution patterns reflect differences in growth kinetics, substrate affinity, and resistance to environmental perturbations among genera.

Nitrogen-loading fluctuations significantly impact anammox community structure and function. Studies demonstrate a nonlinear relationship between nitrogen loading and system performance, with both excessive nitrogen loading (>3.68 kg/m³·d) and nitrogen starvation leading to performance deterioration and reduced anammox bacterial abundance [19]. Under such fluctuating conditions, microbial communities adopt modular collaboration strategies, with increased inter-module connectivity indicating reinforced interspecies interactions to resist loading stress [19].

Table 2: Environmental Parameters Governing Anammox Bacteria Distribution

Environmental Factor	Effect on Anammox Community	Dominant Genera Under Specific Conditions
Salinity	Primary determinant separating marine and freshwater communities	Ca. Scalindua (marine); Ca. Brocadia/Kuenenia (freshwater)
Dissolved Oxygen	Shapes community composition in suboxic environments; most activity in anoxic niches	Ca. Brocadia (shows highest aerotolerance in groundwater)
Nitrogen Loading	Influences dominance patterns; different affinities for ammonium/nitrite	Ca. Kuenenia (high nitrogen); Ca. Jettenia (low nitrogen)
Organic Matter	Affects competition with denitrifiers; some can utilize VFAs	Ca. Brocadia and Ca. Anammoxoglobus (can oxidize VFAs)
Temperature	Influences growth rates and activity optima	Mesophilic (most engineered systems) vs. psychrophilic adaptations

Spatial and Temporal Dynamics

Anammox communities demonstrate dynamic shifts in response to changing environmental conditions, both spatially and temporally. In estuary systems, strong spatial heterogeneity occurs across sediment depths and along salinity gradients, with "Candidatus Scalindua" dominating marine stations while "Candidatus Brocadia" and "Candidatus Kuenenia" increase in abundance in upper estuary regions [11].

Temporal shifts in anammox community composition occur in both natural and engineered systems. Laboratory reactors frequently document population successions, such as shifts from "Candidatus Brocadia fulgida"-dominated populations to "Candidatus Brocadia sp.40," or from "Candidatus Brocadia" to "Candidatus Kuenenia stuttgartiensis" [8]. These successions reflect competitive outcomes influenced by changing reactor conditions and operational parameters.

Keystone Role of Anammox Bacteria in Microbial Communities

Anammox Bacteria as Ecosystem Engineers

Anammox bacteria can function as keystone species in nitrogen-cycling microbial communities, disproportionately impacting community structure and function relative to their abundance. Through their metabolic activity, they create ecological niches for other microorganisms while simultaneously removing fixed nitrogen from the system. In coastal sediments, "Candidatus Scalindua" has been identified as a keystone genus, with rare anammox species potentially playing crucial roles in maintaining ecological stability [11].

The keystone role of anammox bacteria extends to their influence on community assembly processes. Analyses of community assembly mechanisms suggest that ecological drift predominantly shapes overall anammox bacterial communities in coastal sediments, while rare species are more susceptible to dispersal limitations and environmental selection [11]. This differential response to assembly processes between abundant and rare taxa contributes to the overall resilience and functional stability of the ecosystem.

Interaction Networks with Co-occurring Microorganisms

Anammox bacteria participate in complex interaction networks with other microbial functional groups, particularly aerobic ammonia oxidizers and denitrifiers. In eutrophic lake sediments, anammox and denitrifying bacteria demonstrate ecological cooperation, with increasing microbial community stability through their coupled metabolisms [5]. Notably, nirS-type denitrifiers show stronger coupling with anammox bacteria than nirK-type denitrifiers, suggesting specialized functional relationships [5].

Metagenome-assembled genomes-based ecological modeling reveals that dominant denitrifiers can provide various materials including amino acids, cofactors, and vitamins for anammox bacteria [5]. This cross-feeding highlights the importance of microbial interactions for nitrogen removal efficiency and demonstrates how anammox bacteria occupy central positions in metabolic interaction networks. The dependence of anammox bacteria on folate produced by Proteobacteria represents another key interaction, with implications for community stability and function [18].

Research Methods for Studying Anammox Niche Differentiation

Molecular Approaches for Community Analysis

Advanced molecular techniques enable comprehensive analysis of anammox community structure and function. Key methodologies include:

• 16S rRNA gene amplicon sequencing: Targets specific hypervariable regions using anammox-specific primers (e.g., Brod541F and Amx820R) to characterize community composition and diversity [11]. This approach revealed "Candidatus Scalindua" as the dominant anammox bacterium in coastal sediments, particularly in the South China Sea [11].

• Quantitative PCR (qPCR): Quantifies functional genes including 16S rRNA genes of anammox bacteria, hydrazine synthase (hzsB), and nitrite reductase (nirS) genes using TaqMan or SYBR Green systems [21]. This method demonstrated significantly greater abundance of anammox bacteria in artificial versus natural water-level fluctuation zones in the Three Gorges Reservoir area [21].

• Metagenomic and metatranscriptomic sequencing: Provides insights into functional potential and gene expression patterns of anammox communities [5] [17]. Metatranscriptomic analysis of aquifer systems revealed active expression of anammox and aerotolerance genes by "Candidatus Brocadiae" genomes under both oxic and dysoxic conditions [17].

Process Rate Measurements

Combining molecular approaches with process rate measurements enables correlation of community structure with function:

• Isotope pairing techniques: Use 15N-labeled substrates (15NH4+ or 15NO2-) to quantify anammox rates in environmental samples and distinguish them from denitrification [11].

• Nitrogen removal efficiency calculations: Monitor influent and effluent nitrogen species in reactor systems to calculate removal rates under different operational conditions [19] [20]. Studies demonstrate nitrogen removal rates up to 1.13 kg N m−3 d−1 in optimized anammox reactors [20].

• Inhibition assays: Employ specific inhibitors like allylthiourea to distinguish anammox from nitrification processes [22].

Physiological Characterization

Determining physiological parameters essential for niche differentiation requires specialized approaches:

• Enrichment cultures: Establish laboratory cultures under controlled conditions to isolate specific anammox bacteria and characterize their physiological traits [8]. Successful enrichment from various inocula including denitrifying sludge ("Candidatus Jettenia caeni"), activated sludge ("Candidatus Brocadia fulgida"), and freshwater sediments ("Candidatus Scalindua") demonstrates differential cultivation requirements [8].

• Microelectrode measurements: Profile oxygen, ammonium, and nitrite gradients at micrometer scales to identify anammox activity zones in biofilms, granules, and sediments [8]. These measurements reveal considerable physicochemical heterogeneity even in small microbial aggregates.

• Kinetic parameter determination: Calculate maximum specific growth rates (μmax) and half-saturation constants (Ks) through controlled batch experiments, providing essential parameters for Monod model predictions of competitive outcomes [8].

Table 3: Essential Research Reagents and Methodologies for Anammox Research

Research Tool	Specific Application	Key Function in Niche Differentiation Studies
Anammox-specific primers (Brod541F/Amx820R)	16S rRNA gene amplification	Target-specific amplification for community analysis
Functional gene primers (hzsB, nirS, nirK)	Quantification of functional genes	Link genetic potential with nitrogen cycling functions
15N-labeled substrates (15NH4+, 15NO2-)	Isotope pairing techniques	Direct measurement of process rates in complex environments
Synthropic wastewater	Enrichment cultures	Selective cultivation of specific anammox taxa
Microsensors (O2, NH4+, NO2-)	Gradient measurements in biofilms	Identify microniches and activity zones
Metagenomic sequencing kits	Community genomic analysis	Reveal metabolic potential and adaptations

Experimental Protocols for Key Investigations

Enrichment and Cultivation Protocol

Successful enrichment of anammox bacteria requires specific conditions favoring their slow growth and specific metabolic requirements:

• Inoculum collection: Collect biomass from anammox-positive environments (wastewater treatment plants, freshwater sediments, or marine sediments) using anaerobic sampling techniques [5].

• Reactor setup: Configure up-flow anaerobic sludge blanket (UASB) reactors or membrane bioreactors (MBR) with temperature control (30-35°C) and complete light exclusion [23] [20].

• Medium composition: Prepare synthetic wastewater containing NH4+ (50-300 mg N/L) and NO2- (50-400 mg N/L) at approximately 1:1.32 molar ratio, essential minerals (KH2PO4, CaCl2·2H2O, MgSO4·7H2O, KHCO3), and trace elements including FeSO4, EDTA, and vitamin solutions [19] [20].

• Operational parameters: Maintain anoxic conditions (DO < 0.05 mg/L), pH 7.0-8.0, and appropriate hydraulic retention time (HRT) based on nitrogen loading rates [20].

• Monitoring: Regularly measure influent and effluent NH4+, NO2-, and NO3- concentrations to calculate nitrogen removal rates and monitor process stability [19].

This protocol successfully enriched anammox bacteria from eutrophic lake sediments within 365 days, achieving maximum NH4+ and NO2- removal efficiencies of 85.92% and 95.34%, respectively [5].

Community Analysis Protocol

Comprehensive analysis of anammox community structure involves molecular biological approaches:

• DNA extraction: Extract total genomic DNA from 0.5 g wet sample using commercial soil DNA extraction kits, with mechanical lysis for efficient cell disruption [11].

• PCR amplification: Amplify anammox-specific 16S rRNA gene fragments using primers Brod541F and Amx820R with the following thermal program: initial denaturation at 95°C for 5 min; 35 cycles of 95°C for 45 s, 56°C for 30 s, and 72°C for 50 s; final extension at 72°C for 10 min [11].

• High-throughput sequencing: Purify amplicons and sequence using Illumina MiSeq or similar platforms with 2×250 bp or 2×300 bp paired-end sequencing [11].

• Bioinformatic analysis: Process raw sequences through quality filtering, chimera removal, and OTU clustering at 97-98% similarity. Classify OTUs using specialized anammox databases [11].

• Statistical analysis: Calculate diversity indices (Shannon, ACE), conduct multivariate analyses (NMDS, ANOSIM), and construct co-occurrence networks to identify ecological patterns [11].

This protocol revealed significant spatial heterogeneity in anammox communities across estuaries, with distinct distribution patterns for rare species [11].

Implications for Environmental Management and Biotechnology

Understanding anammox niche differentiation has significant practical implications for environmental management and biotechnology. In wastewater treatment, selecting appropriate anammox species for specific wastewater characteristics can optimize treatment efficiency. For instance, "Candidatus Brocadia" may be preferable for systems with organic carbon fluctuations due to its metabolic versatility, while "Candidatus Kuenenia" might suit high-strength ammonium wastewater [8] [18].

The knowledge of niche differentiation also informs bioreactor management strategies. Studies demonstrate that nitrogen-loading fluctuations significantly impact anammox community structure and function, with excessive loading (>3.68 kg/m³·d) or nitrogen starvation leading to performance deterioration [19]. By understanding the ecological preferences of different anammox taxa, operators can implement control strategies that maintain optimal conditions for the desired species.

In natural ecosystems, the niche differentiation of anammox bacteria influences nitrogen cycling and ecosystem responses to environmental change. The expansion of oxygen minimum zones in the ocean may favor the low-ammonia concentration (LAC) ecotype of anammox bacteria, potentially altering nitrogen transformation patterns in marine systems [22]. Similarly, anthropogenic nitrogen loading to freshwater systems may shift anammox community composition with consequences for nitrogen removal capacity.

Ecological niche differentiation among anammox bacteria represents a fundamental mechanism structuring nitrogen-cycling communities across diverse ecosystems. This differentiation arises from physiological variations in growth kinetics, salinity tolerance, oxygen sensitivity, and metabolic versatility, creating distinct environmental preferences among different anammox genera. The clear niche partitioning between marine ("Candidatus Scalindua") and freshwater ("Candidatus Brocadia," "Candidatus Kuenenia," "Candidatus Jettenia") taxa highlights the importance of salinity as a primary environmental filter.

Beyond their functional importance in nitrogen cycling, certain anammox bacteria function as keystone species in microbial communities, disproportionately impacting community structure and ecosystem function. Their interactions with aerobic ammonia oxidizers, denitrifiers, and accessory microorganisms create complex ecological networks that influence nitrogen transformation pathways and efficiencies. Understanding these interactions is essential for predicting ecosystem responses to environmental change and optimizing anammox-based biotechnologies.

Future research should focus on further elucidating the physiological basis of niche differentiation, particularly through comparative genomics of closely related species with different environmental distributions. Additionally, investigating the ecological role of rare anammox taxa may reveal their importance in community resilience and functional stability. As molecular techniques continue to advance, particularly in single-cell approaches and meta-omics integration, our understanding of anammox ecology will undoubtedly deepen, revealing new dimensions of their niche specialization and ecological significance.

The Critical Role of Rare Species in Maintaining Community Stability and Function

The anaerobic ammonium oxidation (anammox) process represents one of the most significant discoveries in microbial nitrogen cycling, capable of converting ammonium directly to dinitrogen gas under anoxic conditions. While predominant anammox bacteria such as Candidatus Scalindua, Candidatus Brocadia, and Candidatus Kuenenia have received substantial scientific attention, recent ecological investigations have revealed that rare microbial species play disproportionately critical roles in maintaining community stability and functional integrity. These low-abundance taxa, often constituting less than 0.1% of relative abundance in microbial communities, serve as keystone components that enhance ecosystem resilience through multiple mechanisms including functional redundancy, niche differentiation, and ecological memory [24]. Within anammox systems, these rare species form intricate interaction networks that stabilize community dynamics against environmental fluctuations, thereby ensuring the continuity of nitrogen removal functions essential for both natural biogeochemical cycling and engineered wastewater treatment systems.

The investigation of rare species in anammox environments represents a paradigm shift in microbial ecology, moving beyond the focus on dominant taxa to understand how minority populations contribute to ecosystem services. This whitepaper synthesizes cutting-edge research on the functional significance of rare anammox bacteria, detailing experimental methodologies for their characterization, quantifying their contributions to community stability, and proposing mechanistic frameworks through which these hidden players maintain system functionality under varying environmental conditions.

Ecological Framework of Rare Species in Anammox Communities

Defining Rare Taxa in Microbial Ecology

In anammox bacterial communities, researchers typically classify microbial populations based on their relative abundance and distribution patterns across samples. Based on established ecological frameworks applied to anammox systems, the classification scheme includes:

• Rare Taxa (RT): Populations with an abundance less than 0.1% across all examined samples [24]
• Conditionally Rare Taxa (CRT): Populations that are generally rare (below 1% abundance) but may exceed the 0.1% threshold in specific environmental conditions
• Moderate Taxa (MT): Populations maintaining abundances between 0.1% and 1% across samples
• Abundant Taxa (AT): Populations consistently exceeding 1% relative abundance in the community [24]

This classification is not merely statistical but reflects fundamental ecological strategies. While abundant taxa typically exhibit greater dispersal capabilities and broader environmental tolerance, rare taxa often demonstrate high specialization to specific microenvironmental conditions and exhibit stronger sensitivity to environmental selection pressures [24]. This specialization enables rare species to occupy distinct metabolic niches that complement the activities of dominant community members.

Community Assembly Mechanisms

The formation and maintenance of anammox bacterial communities are governed by complex assembly mechanisms that differentially affect abundant and rare taxa. Comprehensive studies across diverse estuarine and marine environments, including the Changjiang Estuary (CJE), Oujiang Estuary (OJE), Jiulong River Estuary (JLE), and the South China Sea (SCS), have demonstrated that ecological drift predominantly shapes the overall anammox bacterial community structure in coastal sediments [24]. However, the relative influence of various assembly processes differs significantly between abundant and rare microbial fractions:

Table: Community Assembly Mechanisms for Anammox Bacteria

Assembly Mechanism	Impact on Abundant Taxa	Impact on Rare Taxa
Ecological Drift	Primary influence	Moderate influence
Dispersal Limitation	Weaker influence	Stronger influence
Environmental Selection	Moderate influence	Stronger influence
Homogenizing Dispersal	Variable influence	Weaker influence

This differential susceptibility to assembly processes creates a dynamic where rare species are more strongly affected by spatial heterogeneity and local environmental conditions, leading to distinct distribution patterns across geographic gradients [24]. The higher sensitivity of rare taxa to environmental filters makes them particularly responsive to ecosystem changes, positioning them as potential bioindicators of environmental perturbation.

Methodological Approaches for Studying Rare Anammox Species

Sample Collection and Processing

Investigating rare anammox bacteria requires meticulous sampling strategies to ensure adequate representation of low-abundance populations. Core sampling protocols from recent studies involve:

• Sediment Core Collection: Obtain undisturbed sediment cores (typically 15-34 cm depth) using specialized coring equipment aboard research vessels [24]
• High-Resolution Subsampling: Section cores at fine intervals (1-4 cm depending on core depth) to capture vertical stratification of anammox communities [24]
• Chemical Parameter Analysis: Measure sediment characteristics including organic carbon and nitrogen content using elemental analyzers, and quantify porewater nutrients (NO₃⁻, NO₂⁻, NH₄⁺) with nutrient auto-analyzers [24]
• Dissolved Oxygen Profiling: Characterize oxygen gradients at high spatial resolution (0.2 mm intervals) using oxygen microsensors to identify anoxic niches [24]

For DNA-based analyses, preserve 0.5 g of wet sediment samples immediately after collection using appropriate preservation buffers and store at -20°C until extraction to prevent microbial community changes [24].

Molecular Characterization of Anammox Communities

Advanced molecular techniques enable researchers to detect and quantify rare anammox bacteria despite their low abundance:

• DNA Extraction: Use commercial soil DNA extraction kits (e.g., FastDNA SPIN Kit) with mechanical lysis to ensure efficient cell disruption of diverse bacterial groups [24]
• Targeted Amplification: Employ anammox-specific 16S rRNA gene primers (Brod541F and Amx820R) with 35 PCR cycles to ensure detection of rare variants [24]
• High-Throughput Sequencing: Utilize Illumina platforms to generate sufficient sequencing depth (typically >1 million raw sequences) to capture rare members [24]
• Sequence Processing: Denoise raw sequences using Sickle, remove chimeras with QIIME 2, and cluster operational taxonomic units (OTUs) at 98% similarity threshold [24]
• Taxonomic Assignment: Classify sequences against specialized anammox bacterial databases rather than general repositories to improve accuracy for rare lineages [24]

For enhanced functional insights, genome-centric metagenomics can be employed to recover metagenome-assembled genomes (MAGs) of both abundant and rare community members, enabling reconstruction of their metabolic potential [25].

Bioinformatics and Statistical Analyses

Specialized bioinformatic approaches are required to accurately characterize rare members within anammox communities:

• Alpha Diversity Analysis: Calculate Shannon-Wiener diversity and Abundance-based Coverage Estimator (ACE) richness indices using QIIME 2 to assess within-sample diversity [24]
• Beta Diversity Assessment: Perform NMDS (non-metric multidimensional scaling) based on Bray-Curtis dissimilarity and UniFrac distances to visualize community differences [24]
• Network Analysis: Construct co-occurrence networks using correlation approaches to identify keystone taxa and interaction patterns between rare and abundant members [24]
• Phylogenetic Analysis: Reconstruct phylogenetic relationships of anammox lineages to identify evolutionary patterns among rare and dominant taxa [24]
• Differential Abundance Testing: Apply appropriate statistical models (e.g., DESeq2, edgeR) that account for compositionality and sparse distributions in rare taxa analysis

The following experimental workflow illustrates the integrated methodology for investigating rare anammox bacteria:

Quantitative Analysis of Rare Species Contributions

Diversity Patterns Across Environmental Gradients

Comprehensive studies across Chinese coastal systems have revealed distinct distribution patterns of rare anammox bacteria along spatial and environmental gradients. The analysis of three estuaries (Changjiang, Oujiang, and Jiulong River) and the South China Sea demonstrated significant spatial heterogeneity in anammox community composition, characterized by distinct distribution patterns for rare species [24]. Notably, the Jiulong River Estuary (JLE) exhibited the highest Shannon's diversity index, reflecting enhanced species richness and evenness, while the South China Sea (SCS) showed the lowest diversity [24]. Interestingly, the Changjiang Estuary (CJE) demonstrated the highest species richness despite moderate diversity indices, suggesting the presence of numerous rare species contributing to richness metrics.

Table: Anammox Bacterial Diversity Across Coastal Environments

Location	Shannon Diversity Index	Species Richness	Dominant Genera	Rare Taxa Characteristics
Jiulong River Estuary (JLE)	Highest	High	Ca. Brocadia, Ca. Kuenenia	Thrive in native habitats with higher ammonium
Changjiang Estuary (CJE)	Moderate	Highest	Ca. Scalindua	Distinct rare species composition
Oujiang Estuary (OJE)	Moderate	Moderate	Mixed community	Intermediate characteristics
South China Sea (SCS)	Lowest	Lowest	Ca. Scalindua	Limited rare species diversity

The distribution of specific anammox genera further highlights niche differentiation processes. Candidatus Scalindua dominated marine sediments, particularly in the South China Sea, while Candidatus Brocadia and Candidatus Kuenenia were more abundant in estuarine environments, especially the Jiulong River Estuary [24]. Phylogenetic analyses revealed that Candidatus Scalindua exhibited greater diversity compared to other genera, with rare lineages within this genus contributing significantly to this phylogenetic breadth [24].

Co-occurrence Network Topology and Keystone Species

Network analysis of anammox communities in coastal sediments has revealed that rare species play crucial roles in maintaining ecological stability, with Candidatus Scalindua identified as a keystone genus despite variations in its relative abundance across environments [24]. The topological properties of anammox co-occurrence networks demonstrate that:

• Rare species frequently occupy structurally important positions within microbial association networks
• Low-abundance taxa form connective bridges between modules of more abundant organisms
• The removal of certain rare taxa from network models causes disproportionate disintegration of network structure
• Conditionally rare taxa may act as ecological backups that maintain functional stability during environmental perturbations

These network properties suggest that rare anammox bacteria enhance community stability through functional complementarity and by providing ecological resilience to environmental fluctuations. The hidden nature of these contributions means that conventional abundance-based metrics significantly underestimate their ecological importance.

Functional Significance of Rare Species in Anammox Systems

Metabolic Specialization and Niche Differentiation

Rare anammox bacteria contribute to community function through metabolic specialization that complements the activities of dominant taxa. Genomic analyses of anammox granules have revealed substantial functional diversity among community members, with different bacterial groups specializing in particular metabolic transformations [25]. For instance, Chlorobi-affiliated bacteria in anammox systems function as highly active protein degraders, catabolizing extracellular peptides while recycling nitrate to nitrite [25]. This activity supports the anammox process by regenerating essential substrates and maintaining redox balance.

Other heterotrophic bacteria associated with anammox communities contribute to scavenging detritus and peptides produced by anammox bacteria, potentially using alternative electron donors including Hâ‚‚, acetate, and formate to fuel their energy metabolism [25]. This metabolic versatility enables rare species to occupy distinct niches that would otherwise remain unexploited, thereby increasing overall ecosystem efficiency and resource utilization.

Ecological Resilience and Response to Perturbation

Rare species significantly enhance the stability and resilience of anammox communities when facing environmental fluctuations. The "insurance hypothesis" in microbial ecology proposes that rare taxa serve as a reservoir of genetic and functional diversity that can become important under changing conditions. In anammox systems, this manifests through several mechanisms:

• Functional Redundancy: Rare taxa provide backup capacity for critical functions if dominant species are compromised
• Stress Response: Specific rare lineages exhibit enhanced tolerance to environmental stressors like salinity shifts or organic matter fluctuations
• * Metabolic Flexibility*: Rare species often possess specialized metabolic genes that enable survival during nutrient limitation
• Recovery Potential: Following system disturbance, rare taxa can proliferate and maintain ecosystem functions during community reassembly

Experimental evidence from reactor studies demonstrates that anammox communities with higher diversity, including representation of rare species, maintain more stable nitrogen removal performance when confronted with operational perturbations such as load variations or toxic shocks [26]. This functional resilience has practical implications for the design and operation of engineered anammox systems for wastewater treatment.

Research Reagent Solutions for Anammox Studies

Table: Essential Research Reagents and Materials for Anammox Community Analysis

Reagent/Material	Specific Example	Application Purpose	Technical Considerations
DNA Extraction Kit	FastDNA SPIN Kit for Soil	Efficient lysis of diverse bacterial cells	Mechanical beating enhances DNA yield from tough cells
PCR Primers	Brod541F / Amx820R	Specific amplification of anammox bacterial 16S rRNA genes	98% similarity threshold for OTU clustering
Elemental Analyzer	Carlo-Erba EA 2100	Measurement of organic C and N content in sediments	Acidification step removes inorganic carbon
Nutrient Auto-analyzer	AA3 Bran+Luebbe	Quantification of NO₃⁻, NO₂⁻, NH₄⁺ in porewater	High sensitivity for low nutrient concentrations
Oxygen Microsensor	OX50 Unisense	High-resolution DO profiling in sediments	0.2 mm spatial resolution, requires 2-point calibration
Sequencing Platform	Illumina series	High-throughput sequencing of anammox communities	>1 million reads recommended for rare taxa detection
Bioinformatics Package	QIIME 2	Processing and analysis of sequencing data	Specific anammox database improves taxonomic assignment

The critical role of rare species in maintaining community stability and function represents a fundamental paradigm shift in our understanding of anammox ecosystems. Rather than mere ecological passengers, these low-abundance taxa serve as keystone components that enhance functional resilience, provide metabolic versatility, and stabilize community dynamics through intricate interaction networks. The evidence from coastal environments demonstrates that rare anammox bacteria exhibit distinct distribution patterns shaped primarily by dispersal limitation and environmental selection, in contrast to the predominantly stochastic assembly of abundant community members [24].

From a practical perspective, recognizing the importance of rare species has profound implications for biotechnological applications of anammox processes in wastewater treatment. Engineering strategies that promote microbial diversity, including the conservation of rare taxa, may enhance system stability and operational performance under fluctuating conditions [26]. Future research directions should focus on elucidating the specific metabolic contributions of rare anammox bacteria, understanding their dynamics during system perturbations, and developing bioengineering approaches to maintain their beneficial functions in engineered ecosystems. By integrating knowledge of rare species into both ecological theory and biotechnological practice, we can advance toward more sustainable and resilient nitrogen removal systems that harness the full functional potential of microbial communities.

The keystone module concept represents an evolution beyond the classical view of single keystone species, describing a group of multiple keystone species that demonstrate correlated occurrence and collective function within microbial consortia. These modules exert a disproportionately large influence on community structure, stability, and function—greater than would be predicted from the sum of their individual effects [1] [27]. In engineered and natural ecosystems, keystone modules frequently emerge as critical regulatory units that maintain functional robustness against environmental perturbations through synergistic interactions among member taxa.

Within anaerobic ammonium oxidation (anammox) systems, keystone modules play a particularly vital role in maintaining nitrogen removal efficiency under fluctuating operational conditions. The anammox process, mediated by Planctomycetes bacteria, converts ammonium and nitrite directly to dinitrogen gas under anoxic conditions and has revolutionized wastewater treatment over recent decades [28]. However, anammox bacteria are notoriously sensitive to environmental fluctuations, making the stabilizing influence of keystone modules essential for reliable process operation [6] [29]. This technical guide examines the identification, functional mechanisms, and ecological significance of keystone modules within anammox consortia, providing researchers with advanced methodologies for investigating these critical microbial components.

Identification and Characterization of Keystone Modules

Analytical Frameworks for Detection

Identifying keystone modules requires analytical approaches that capture both the individual keystoneness of taxa and their correlated behavior within communities. The Empirical Presence-Abundance Interrelation (EPI) framework offers a top-down methodology that detects keystone taxa by their total influence on community composition without requiring reconstruction of detailed interaction networks [1] [27]. This approach measures how the presence or absence of specific taxa correlates with community-wide abundance profiles, identifying candidate keystones that subsequently can be evaluated for modular co-occurrence patterns.

Figure 1: Analytical workflow for keystone module identification combining cross-sectional data analysis with perturbation experiments.

Network-Based Characterization

Co-occurrence network analysis serves as a powerful complementary approach for detecting keystone modules by revealing interaction patterns that may not be apparent through abundance-based metrics alone. In this methodology, microbial associations are inferred from correlation patterns in abundance data across samples, with nodes representing taxa and edges representing significant positive or negative associations [30] [11]. Keystone modules typically appear as highly interconnected subnetworks with numerous connections to taxa outside their module, positioned as critical hubs within the broader community network.

Application of this approach in anammox reactors has demonstrated that keystone modules frequently include not only anammox bacteria (e.g., Candidatus Jettenia, Candidatus Brocadia) but also cooperative heterotrophs from phyla such as Chloroflexi and Proteobacteria [31] [29]. These associated taxa provide essential functional support to anammox bacteria through metabolic exchanges, with their correlated dynamics serving as a stabilizing mechanism during environmental fluctuations. Research indicates that modules with higher within-module connectivity demonstrate greater functional resilience to nitrogen loading variations in anammox systems [6].

Table 1: Quantitative Metrics for Keystone Module Identification in Microbial Networks

Metric	Calculation Method	Interpretation	Threshold Value
Modularity Index	Quality function optimization of network division into modules	Measures degree of compartmentalization in network; values >0.4 indicate significant modular structure	0.4-0.6 [6]
Within-Module Connectivity (Zi)	Number of connections from a node to other nodes in its own module	Identifies hubs within modules; Zi >2.5 indicates module hubs	>2.5 [30]
Among-Module Connectivity (Pi)	Distribution of a node's connections across different modules	Measures connector nodes; Pi >0.6 indicates connectors between modules	>0.6 [30]
EPI Value	Community importance based on presence-impact measurements	Quantifies keystone influence; higher values indicate greater keystoneness	Species-specific [1]

Keystone Modules in Anammox Systems

Structural and Functional Roles

In anammox bioreactors and natural environments, keystone modules organize around core anammox bacteria that establish the foundational nitrogen-removing capacity, with associated taxa providing critical functional enhancements and stabilization services. The anammox bacteria themselves, including genera such as Candidatus Brocadia, Candidatus Kuenenia, Candidatus Jettenia, and Candidatus Scalindua, typically function as structural anchors within these modules, while cooperating heterotrophs from Proteobacteria, Chloroflexi, and Acidobacteria phyla serve as stabilizing elements [31] [29] [28].

These modules maintain functional stability through different mechanisms depending on environmental conditions. Under optimal nitrogen loading, keystone modules demonstrate balanced internal interactions with both positive and negative associations. However, under suboptimal conditions such as nitrogen starvation or excessive loading, modules undergo structural reorganization characterized by increased modularity and strengthened negative interactions—interpreted as a stress-response mechanism that preserves core functions [6]. This structural flexibility enables anammox systems to maintain nitrogen removal efficiency despite fluctuating influent conditions.

Response to Environmental Stressors

Keystone modules demonstrate remarkable plasticity when confronted with environmental stressors, reorganizing their interaction networks to preserve system functionality. Under temperature stress, cross-kingdom keystone modules (incorporating both bacterial and fungal taxa) display altered interaction strength that correlates directly with stability of microbial carbon metabolic activity [30]. Similarly, under nitrogen loading fluctuations in anammox reactors, keystone modules increase their modular organization (with indices rising to 0.563 during inhibition phases) as a coordinated response to loading stress [6].

Table 2: Keystone Module Responses to Environmental Stressors in Anammox Systems

Stress Type	Module Response	Functional Outcome	Reference
High Nitrogen Loading (>3.68 kg/m³·d)	Increased modularity (index 0.563) and negative interactions	Preservation of anammox activity despite inhibition	[6]
Nitrogen Starvation	Enhanced module connectivity and cooperation	Maintenance of community structure during substrate limitation	[6]
Temperature Stress	Altered cross-kingdom interaction strength	Impacts carbon metabolic stability	[30]
Salinity Perturbation	Closer, more complex connections	Functional adaptation to osmotic stress	[29]
Organic Carbon Addition	Reinforced keystone taxa associations	Enhanced system resilience	[29]

Importantly, keystone modules may incorporate rare taxa that exert influence disproportionate to their abundance. In coastal sediments, Candidatus Scalindua has been identified as a keystone genus despite low relative abundance in some environments, with rare anammox taxa playing crucial roles in maintaining ecological stability of the broader anammox community [11]. This challenges conventional abundance-focused community analyses and highlights the necessity of specialized approaches for detecting these critical but low-abundance components.

Experimental Protocols for Keystone Module Analysis

Reactor Operation and Sampling

Investigating keystone modules in anammox systems requires carefully controlled bioreactor operation with systematic monitoring. The following protocol outlines standard methodology for establishing an anammox enrichment bioreactor suitable for keystone module research:

• Reactor Configuration: Employ an expanded granular sludge bed (EGSB) or up-flow anaerobic sludge bed (UASB) reactor with 5-10 L working volume constructed of plexiglass. Maintain anoxic conditions by continuous argon gas purging (0.5 LPM for 30 minutes before operation) and cover with insulating material to block light and maintain temperature at 34±1°C [6] [31].

• Inoculation: Seed with anammox biomass (approximately 3.0 g/L volatile suspended solids) from established anammox reactors or natural sediments. For studies targeting specific environments, sediment cores (0-20 cm depth) from eutrophic lakes or coastal areas provide appropriate inoculation material [31] [11].

• Feeding Medium Preparation: Prepare synthetic wastewater containing (per liter): (NHâ‚„)â‚‚SOâ‚„ and NaNOâ‚‚ as nitrogen sources (concentration varies by experimental phase), NaHCO₃ (0.5 g), KHCO₃ (0.5 g), KHâ‚‚POâ‚„ (0.027 g), MgSO₄·7Hâ‚‚O (0.02 g), CaCl₂·2Hâ‚‚O (0.136 g), and trace element solutions I and II (1 mL and 1.2 mL respectively) [6]. Maintain NH₄⁺:NO₂⁻ ratio at approximately 1:1.32 based on anammox stoichiometry [28].

• Operational Phases: Structure experiment into distinct operational phases with varying nitrogen loading rates (NLR). A typical design includes: Phase 1 (stable baseline, NLR=1.38±0.01 kg/m³·d), Phase 2 (incremental NLR increase), Phase 3 (starvation and recovery), and Phase 4 (inhibition testing) [6].

• Monitoring: Daily measurement of influent and effluent NH₄⁺-N, NO₂⁻-N, and NO₃⁻-N concentrations to calculate nitrogen removal rates. Regular sampling of biomass for molecular analysis at the end of each operational phase [6] [29].

Molecular Analysis and Network Construction

Comprehensive characterization of keystone modules requires integrated molecular approaches to resolve taxonomic composition, functional potential, and interaction networks:

• DNA Extraction: Extract total genomic DNA from 0.5 g biomass samples using commercial soil DNA extraction kits (e.g., FastDNA SPIN Kit, MP Biomedical) following manufacturer protocols. Quantify and quality-check DNA using fluorometric and spectrophotometric methods [11].

• Amplification and Sequencing:

  • For total bacterial community: Amplify 16S rRNA gene V3-V4 regions using primers 341F/806R [31].
  • For anammox-specific analysis: Amplify anammox 16S rRNA genes using primer set Brod541F/Amx820R with thermal program: 95°C for 5 min; 35 cycles of 95°C for 45s, 56°C for 30s, 72°C for 50s; final extension at 72°C for 10 min [11].
  • For functional genes: Target marker genes such as hydrazine synthase (hzsB) for anammox potential and nirS/nirK for denitrification potential [31].

• Sequence Processing and Analysis: Process raw sequences through quality filtering, chimera removal, and OTU clustering at 97-98% similarity threshold using QIIME2 or similar pipelines. Classify sequences against specialized databases (e.g., anammox-specific 16S rRNA database) for accurate taxonomic assignment [11].

• Network Construction:

  • Calculate correlation matrices (SparCC, Spearman, or Pearson) from OTU abundance tables, applying appropriate thresholds for statistical significance.
  • Construct co-occurrence networks in R using igraph or similar packages, defining nodes as OTUs and edges as significant correlations.
  • Detect modules using greedy modularity optimization algorithms and calculate topological parameters (degree centrality, betweenness centrality, modularity index) [30] [11].

• Keystone Module Identification: Integrate network topology metrics with cross-sectional EPI analysis to identify candidate keystone modules. Validate through longitudinal sampling or targeted perturbation experiments where feasible [1].

Figure 2: Molecular workflow for keystone module analysis from biomass sampling to network identification.

Research Reagent Solutions for Keystone Module Studies

Table 3: Essential Research Reagents for Keystone Module Investigation in Anammox Systems

Reagent/Category	Specific Examples	Function/Application	Reference
DNA Extraction Kits	FastDNA SPIN Kit for Soil (MP Biomedical)	High-quality metagenomic DNA extraction from complex biomass samples	[11]
PCR Primers	Brod541F/Amx820R (anammox 16S rRNA), 341F/806R (general bacterial), nirS/nirK primers (denitrification)	Target-specific amplification for community composition and functional potential analysis	[31] [11]
Sequencing Kits	Illumina MiSeq Reagent Kits v2/v3	High-throughput amplicon sequencing for community profiling	[31] [11]
Trace Element Solutions	Solution I: EDTA + FeSO₄·7H₂O; Solution II: EDTA + NaMoO₄·2H₂O + NiCl₂·6H₂O + CuSO₄·5H₂O + CoCl₂·6H₂O + ZnSO₄·7H₂O + MnCl₂·4H₂O	Essential micronutrients for anammox and heterotrophic bacterial growth in synthetic media	[6]
Synthetic Wastewater Components	(NH₄)₂SO₄, NaNO₂, NaHCO₃, KHCO₃, KH₂PO₄, MgSO₄·7H₂O, CaCl₂·2H₂O	Controlled substrate and nutrient sources for reactor operation	[6] [31]
Network Analysis Tools	QIIME2, igraph (R), Mothur, PRIMER-e	Bioinformatics processing, statistical analysis, and network construction	[30] [11]

Keystone modules represent a fundamental organizational principle within complex microbial consortia, particularly in engineered anammox systems where functional stability depends on coordinated multispecies interactions. The identification of these modules through integrated approaches combining network analysis with top-down keystone detection frameworks provides unprecedented insights into the ecological mechanisms underlying process stability. Moving forward, several research priorities emerge:

First, developing standardized methodologies for dynamic network analysis will enable researchers to capture temporal patterns in keystone module organization and function, revealing how these modules reassemble in response to environmental fluctuations. Second, integrating metatranscriptomic and metaproteomic approaches with co-occurrence network analysis will elucidate the molecular mechanisms underlying module interactions. Finally, applying keystone module principles to bioprocess engineering offers promising pathways for designing more stable and efficient wastewater treatment systems through targeted management of critical microbial components.

As research in this field advances, the keystone module concept promises to transform our understanding of microbial community ecology while delivering practical strategies for managing complex bioprocess systems. The methodologies and frameworks presented in this technical guide provide researchers with essential tools for investigating these critical ecological units in anammox systems and beyond.

From Detection to Direction: Identifying and Engineering with Anammox Keystone Species

Top-Down Identification Frameworks for Keystone Taxa Without Network Reconstruction

The stability and function of anaerobic ammonium oxidation (anammox) bacterial communities are critical for efficient nitrogen removal in both natural and engineered ecosystems. Within these complex microbial consortia, keystone taxa exert disproportionately large effects on community structure and function relative to their abundance [32]. Identifying these key players is essential for understanding and manipulating these systems for enhanced wastewater treatment and biogeochemical cycling. Traditional methods for identifying keystone species rely heavily on reconstructing detailed correlation networks from compositional data, an approach fraught with challenges due to compositional effects, limited sample sizes, and the inherent difficulty in distinguishing correlation from causation [32].

A paradigm shift is emerging toward top-down identification frameworks that detect keystones by their total influence on other taxa without presupposing pairwise interactions or specific underlying dynamics [32]. When applied to anammox systems, these approaches reveal that keystone species—including various Candidatus genera such as Scalindua, Brocadia, Kuenenia, and Jettenia—are often embedded within keystone modules comprising multiple species with correlated occurrence patterns [32] [11]. These findings fundamentally reshape our understanding of stability and function in anammox communities, suggesting that rare species often play critically important roles in maintaining ecological stability [11] [33].

Theoretical Foundation of Top-Down Identification

Defining Keystone Taxa: From Concept to Metric

The ecological concept of keystone species, first introduced by Paine in 1969, generally refers to native taxa that play an especially important role in ecosystem stability [32]. Power et al. later formalized this concept through "community importance," which evaluates a species' effect on ecosystem traits such as productivity, nutrient cycling, and species richness [32]. Two primary operational definitions emerge from this framework:

• Abundance-impact: Measures how changes in a species' relative abundance affect community traits
• Presence-impact: Measures how the complete removal or introduction of a species affects the community [32]

The presence-impact definition aligns particularly well with microbial manipulation techniques (e.g., antibiotics, probiotics, fecal microbiota transplant) and thus provides a practical foundation for keystone identification in anammox systems [32].

The EPI Framework: Empirical Presence-abundance Interrelation

The top-down framework formalizes keystone identification through the Empirical Presence-abundance Interrelation (EPI), which detects candidate keystones from cross-sectional data by measuring the association between a species' presence/absence pattern and community-wide differences in abundance profiles [32]. The framework implements three complementary measures:

• D₁ and Dâ‚‚: Distance-based measures that quantify differences in abundance profiles between samples where a focal taxon is present versus absent
• Q: A modularity-based measure from network science that identifies taxa whose presence/absence partitions the community into distinct modules [32]

This approach measures a taxon's total influence without calculating pairwise correlation networks and does not assume ecological dynamics are governed primarily by pairwise interactions [32]. The resulting framework is applicable to both perturbation experiments and cross-sectional metagenomic surveys, making it particularly valuable for studying anammox communities in natural environments where controlled manipulations are challenging.

Table 1: Core Definitions in Top-Down Keystone Identification

Term	Definition	Application Context
Keystone Taxa	Native taxa with disproportionately large effects on community structure/stability	Anammox communities in reactors and natural environments
Presence-impact (Iáµ¢)	Effect of removing/introducing species i on abundance profiles of other taxa	Measured directly in perturbation experiments
Empirical Presence-abundance Interrelation (EPI)	Association between species' presence/absence and community abundance profiles	Estimated from cross-sectional data; identifies candidate keystones
D₁ and D₂	Distance-based measures of EPI	Quantifies overall community differences based on taxon presence
Q	Modularity-based measure of EPI	Identifies taxa that partition community into distinct modules

Methodological Implementation in Anammox Systems

Experimental Design and Data Collection

Implementing top-down keystone identification in anammox communities requires careful experimental design spanning both observational and manipulative approaches:

Cross-Sectional Survey Design:

• Collect multiple samples across environmental gradients or temporal sequences
• Ensure sufficient replication within and across conditions (≥3-5 samples per group)
• Record key environmental parameters: NH₄⁺, NO₂⁻, organic carbon, temperature, pH, dissolved oxygen [11] [31]
• Preserve samples appropriately for DNA extraction (e.g., FastDNA SPIN Kit) [11]

Perturbation Experiments:

• Implement species removal via antibiotics, phage targeting, or immunodepletion
• Conduct species addition through bioaugmentation [32]
• Monitor community composition pre- and post-perturbation using 16S rRNA gene sequencing
• Measure functional responses: nitrogen removal rates, substrate utilization, gene expression [33]

Molecular Characterization of Anammox Communities:

• Target 16S rRNA gene with anammox-specific primers (Brod541F/Amx820R) [11]
• Quantify functional genes: hydrazine synthase (hzsB), nitrite reductase (nirS) [31] [34]
• Consider metagenomic sequencing for comprehensive functional profiling [31]
• Analyze sequences using QIIME 2 or similar pipelines with anammox-specific databases [11]

Computational Framework and Analysis Pipeline

The analytical implementation involves calculating EPI measures from abundance tables and presence/absence patterns:

Workflow for top-down keystone identification showing the sequence from data input through preprocessing, calculation of EPI measures, to final candidate keystone identification.

Data Preprocessing Steps:

• Filtering: Remove taxa with prevalence <10% across samples
• Normalization: Convert raw counts to relative abundances or apply CSS normalization
• Presence/Absence Thresholding: Define presence as relative abundance >0.01%

Calculation of EPI Measures: For each taxon i across all samples:

• D₁ⁱ: Mean Bray-Curtis distance between samples where taxon i is present versus absent
• D₂ⁱ: Mean Euclidean distance between abundance vectors of samples where taxon i is present versus absent
• Qⁱ: Modularity of the bipartite network connecting samples to taxa, measuring how well taxon i's presence partitions the community

Statistical Validation:

• Compare observed EPI values to null distributions generated via permutation testing
• Apply false discovery rate correction for multiple comparisons
• Validate candidate keystones through stability analysis (e.g., bootstrap resampling)

Table 2: Key Experimental Parameters for Anammox Keystone Identification

Parameter	Recommended Specification	Rationale
Sequencing Depth	≥10,000 reads/sample	Sufficient coverage for rare taxa detection
Sample Replication	≥5 samples per condition	Enables statistical testing of EPI measures
Taxonomic Resolution	97% 16S rRNA similarity or species-level metagenomic bins	Balance between specificity and computational tractability
Abundance Threshold	>0.01% relative abundance	Filtering of spurious taxa while retaining potentially important rare species
Environmental Metadata	NH₄⁺, NO₂⁻, NO₃⁻, organic C, temperature, pH, DO	Contextualize keystone roles in environmental gradients

Applications in Anammox Community Research

Case Study: Keystone Roles in Estuarine Sediments

Application of top-down frameworks to anammox communities in Chinese coastal sediments revealed Candidatus Scalindua as a keystone genus, with rare species playing crucial roles in maintaining ecological stability [11]. The study analyzed anammox bacterial diversity, community structure, and interspecific relationships across three estuaries and the South China Sea, finding that:

• Rare anammox taxa were more susceptible to dispersal limitations and environmental selection
• Keystone species demonstrated higher connectivity in co-occurrence networks
• Community assembly was predominantly shaped by ecological drift, with deterministic processes stronger for rare species [11]

These findings illustrate how top-down identification reveals organizational principles invisible to network-based methods, particularly the importance of rare species in maintaining the anammox bacterial community in coastal sediments.

Case Study: Carrier-Type Influences on Keystone Interactions

Research on anammox biofilm reactors demonstrated that carrier type significantly affects keystone species interactions by molding microhabitat differences [35]. The study compared three carrier types:

• Sponge carriers: Enabled thicker biofilms (>5mm) with abundant anammox bacteria
• Plastic ring carriers: Resulted in different community composition and nitrogen conversion efficiencies
• Nonwoven fabric carriers: Showed rapid biomass accumulation with distinct keystone species

Carrier type imposed selection pressure on microbial recruitment, altering key functional gene abundances and shaping keystone interactions [35]. This demonstrates the utility of top-down frameworks in optimizing bioreactor design by identifying keystone taxa critical for process stability.

Case Study: Stress Response and Community Resilience

Long-term nitrite stress experiments revealed that rare subcommunities play crucial roles in the counteraction of anammox communities to environmental perturbation [33]. When exposed to 200 mg-N/L nitrite:

• Nitrogen removal rates decreased from 0.51±0.02 to 0.16±0.04 kg-N/(m³·d)
• The community exhibited high resilience after discontinuing nitrite treatment
• Keystone shifts occurred from Candidatus Jettenia to Candidatus Brocadia following stress [33]

Abundant and rare subcommunities employed different reactive strategies, with rare species maintaining functional potential during stress and contributing significantly to community recovery [33]. This highlights the value of top-down approaches in predicting community stability under operational stresses.

The Scientist's Toolkit: Essential Research Reagents

Table 3: Key Research Reagents for Anammox Keystone Identification Studies

Reagent/Category	Specific Examples	Function/Application
DNA Extraction Kits	FastDNA SPIN Kit for Soil	Efficient lysis of anammox bacteria; removes PCR inhibitors common in environmental samples
PCR Primers	Brod541F/Amx820R (16S rRNA); hzsB primers for anammox functional genes	Specific amplification of anammox bacterial sequences; detection and quantification of key functional genes
Sequencing Platforms	Illumina MiSeq/HiSeq for 16S rRNA; NovaSeq for metagenomics	High-throughput sequencing for community profiling; enables sufficient depth for rare taxa detection
Bioinformatics Tools	QIIME 2, Mothur, custom R scripts	Processing of sequencing data; calculation of EPI measures; statistical analysis and visualization
Carrier Materials	Nonwoven fabrics, polyethylene plastics, polyurethane sponges	Biofilm support in reactor studies; creates different microhabitats influencing keystone identities
Isotope Tracers	¹⁵NH₄⁺, ¹⁵NO₂⁻	Measurement of anammox process rates via isotope pairing technique; links taxonomy to function
Activity Assays	Specific anammox activity (SAA) tests, heme c quantification	Assessment of functional capacity; correlates keystone status with process performance
Isotaxiresinol 9,9'-acetonide	Isotaxiresinol 9,9'-acetonide, CAS:252333-72-5, MF:C22H26O6, MW:386.4 g/mol	Chemical Reagent
Fmoc-N-Me-His(Trt)-OH	Fmoc-N-Me-His(Trt)-OH, MF:C41H35N3O4, MW:633.7 g/mol	Chemical Reagent

Integration with Bottom-Up Approaches

While top-down frameworks excel at identifying influential taxa without a priori assumptions, their integration with bottom-up approaches provides a more comprehensive understanding of anammox communities [36]. Bottom-up approaches use prior knowledge of metabolic pathways and potential interactions to rationally design and engineer synthetic microbial consortia [36].

The synergy between these approaches is particularly powerful:

• Top-down identification reveals which species naturally exert disproportionate influence
• Bottom-up engineering tests hypothesized mechanisms through controlled assembly
• Combined approaches accelerate the design of stable, efficient anammox systems for wastewater treatment [36]

This integration is advancing anammox biotechnology toward predictive control of community structure and function, with keystone taxa serving as leverage points for process optimization.

Top-down identification frameworks represent a necessary advancement toward reliable identification of keystone taxa in complex anammox communities. By focusing on a taxon's total influence without requiring difficult network reconstruction, these approaches overcome fundamental limitations of traditional methods and provide robust, biologically interpretable results. The application of these frameworks has revealed several fundamental principles of anammox community ecology, including the importance of rare species, context-dependency of keystone identities, and differential responses to environmental stresses.

As anammox biotechnology continues to develop for wastewater treatment and environmental remediation, top-down keystone identification will play an increasingly important role in diagnosing process stability, predicting responses to perturbation, and designing synthetic communities with enhanced functions. The integration of this approach with bottom-up engineering and metabolic modeling represents the future of predictive microbial community management.

Co-occurrence Network Analysis for Pinpointing Keystone Hubs

The concept of the keystone species, introduced by zoologist Robert T. Paine in 1969, describes a species that exerts a disproportionately large effect on its natural environment relative to its abundance [37]. The removal of such a species triggers dramatic shifts in ecosystem structure and function, even if the species itself constitutes a small portion of the community by biomass [37]. While initially described for macroscopic ecosystems, this concept is fundamental to understanding microbial communities. In the context of anammox (anaerobic ammonium oxidation) bacterial communities, which play a critical role in nitrogen loss from estuarine and marine environments, identifying these keystone taxa is essential for predicting ecosystem function and stability [24].

Modern studies reveal that rare species can fulfill these keystone roles in microbial systems. Research on anammox bacteria in Chinese estuaries found that rare species were more susceptible to dispersal limitations and environmental selection, yet Candidatus Scalindua was identified as a keystone genus through co-occurrence network analysis [24]. Similarly, other soil microbiome studies have demonstrated that keystone taxa carry out specialized metabolic functions, such as nitrogen and phosphonate metabolism, which are crucial for maintaining community stability [38]. Co-occurrence network analysis has thus emerged as a powerful statistical tool for pinpointing these critical hubs within complex microbial communities, including those driving the anammox process.

Theoretical Foundations of Co-occurrence Network Analysis

What is a Co-occurrence Network?

A co-occurrence network is a graphical representation of statistically significant relationships between operational taxonomic units (OTUs), amplicon sequence variants (ASVs), or other microbial taxa across multiple samples. In these networks, nodes represent individual microbial taxa, and edges (connecting lines) represent strong positive (e.g., mutualism, synergism) or negative (e.g., competition, antagonism) correlations between them. The structure of these networks reveals the organizational principles of a microbial community and helps identify which members play critically important roles.

Network analysis allows researchers to move beyond simple diversity metrics to understand the complex web of interactions. This is particularly valuable in anammox communities, where the synergistic relationships between different bacterial groups determine the overall rate of nitrogen loss [24]. By analyzing the topology (connection pattern) of these networks, scientists can identify which taxa act as central connectors, and are therefore potential keystone hubs.

Key Network Metrics for Identifying Keystone Hubs

Keystone hubs are typically identified by their unique positions and high connectivity within the network. The table below summarizes the key metrics used to quantify a node's importance.

Table 1: Key Network Metrics for Identifying Keystone Hubs

Metric	Definition	Interpretation for Keystone Status
Within-Module Connectivity (Zi)	Measures how well a node is connected to other nodes within its own module (a group of highly interconnected nodes).	A node with Zi > 2.5 is considered a module hub, indicating it is a central connector within a specific functional group.
Among-Module Connectivity (Pi)	Measures how well a node is connected to nodes in different modules.	A node with Pi > 0.62 is considered a connector, indicating it links different modules and facilitates cross-module communication.
Degree	The total number of connections (edges) a node has to other nodes.	A high degree indicates a highly connected taxon, which may be critical for network stability.
Closeness Centrality	Measures how close a node is to all other nodes in the network.	A high closeness centrality suggests the taxon is efficient in influencing or communicating with the entire network.
Betweenness Centrality	Measures the number of shortest paths that pass through a node.	A high betweenness centrality identifies a taxon that acts as a critical "bridge" between different parts of the network.

Keystone hubs are often defined as nodes that are both module hubs (high Zi) and connectors (high Pi), making them critical for integrating functions within and between different ecological niches in the community [39] [38]. For example, in a study of forested wetland ecotones, researchers identified nine such network hubs as the most important nodes in the microbial co-occurrence networks [39].

Application to Anammox Bacterial Communities

Anammox Community Assembly and Network Patterns

Anammox bacteria are slow-growing chemolithoautotrophs that belong to the phylum Planctomycetes. Key genera include Candidatus Scalindua, Brocadia, Kuenenia, and Jettenia. Understanding the co-occurrence patterns of these bacteria with other community members is vital because their activity is heavily influenced by synergistic and competitive interactions.

A 2025 study analyzed anammox bacterial communities in three estuaries and the South China Sea, providing a clear example of this application [24]. The research revealed significant spatial heterogeneity in the anammox bacteria across the four regions, with distinct distribution patterns for rare species. The analysis of the community assembly mechanism suggested that ecological drift (a stochastic process) predominantly shaped the overall anammox bacterial community in the coastal sediments. However, the study also found that rare species were more susceptible to deterministic processes like dispersal limitations and environmental selection.

Crucially, the co-occurrence network analysis identified Candidatus Scalindua as a keystone genus. Furthermore, it highlighted that rare species may play a indispensable role in maintaining the ecological stability of the anammox bacterial community [24]. This finding challenges the simple view that abundant taxa are always the most important, underscoring the power of network analysis to reveal the hidden roles of low-abundance organisms.

Table 2: Quantitative Data from an Anammox Bacterial Community Study [24]

Parameter	Changjiang Estuary (CJE)	Oujiang Estuary (OJE)	Jiulong River Estuary (JLE)	South China Sea (SCS)
Shannon's Diversity Index	Information Missing	Information Missing	Highest	Lowest
Species Richness (ACE)	Greatest	Information Missing	Information Missing	Information Missing
Predominant Anammox Genus	Candidatus Scalindua	Candidatus Scalindua	Candidatus Brocadia & Kuenenia	Candidatus Scalindua
Community Assembly Driver	Ecological Drift	Ecological Drift	Ecological Drift	Ecological Drift
Role of Rare Species	Susceptible to dispersal limitation & environmental selection, crucial for network stability.	Susceptible to dispersal limitation & environmental selection, crucial for network stability.	Susceptible to dispersal limitation & environmental selection, crucial for network stability.	Susceptible to dispersal limitation & environmental selection, crucial for network stability.

Workflow for Constructing and Analyzing Co-occurrence Networks

The following diagram illustrates the end-to-end workflow for performing a co-occurrence network analysis on amplicon sequencing data, from initial bioinformatics processing to the biological interpretation of keystone hubs.

Detailed Experimental Protocols

Sample Collection and DNA Sequencing for Anammox Studies

The foundational step for a robust network analysis is proper sample collection and sequencing. The protocol from the anammox study in Chinese coastal sediments is exemplary [24].

• Sample Collection: Sediment cores were collected from various estuaries (e.g., Changjiang Estuary, Jiulong River Estuary) and the South China Sea. Cores were subsampled at high-resolution intervals (e.g., 1-2 cm for CJE and SCS, 4 cm for JLE), resulting in 15 to 17 subsamples per core. This high-resolution sampling is critical for capturing fine-scale variations.
• DNA Extraction: Total DNA was extracted from 0.5 g of each wet sediment sample using a commercial kit like the FastDNA SPIN Kit for Soil. The DNA's quantity and purity were assessed using a Qubit fluorometer and NanoDrop spectrophotometer.
• PCR Amplification: The anammox bacterial 16S rRNA gene was specifically amplified using the primer pair Brod541F and Amx820R. The PCR reaction mix included Premix Ex Taq, primers, DNA template, bovine serum albumin (BSA) to reduce inhibition, and RNase-free water.
• Thermal Cycler Program:
  • Initial Denaturation: 95°C for 5 min.
  • 35 Cycles of:
    • Denaturation: 95°C for 45 s.
    • Annealing: 56°C for 30 s.
    • Extension: 72°C for 50 s.
  • Final Extension: 72°C for 10 min.
• Sequencing: The amplicons were verified via gel electrophoresis and then subjected to high-throughput sequencing on an Illumina platform.

Bioinformatics and Statistical Pipeline

After sequencing, the data must be processed to construct the network.

• Sequence Processing: Raw sequences were denoised using programs like Sickle. Operational Taxonomic Units (OTUs) were clustered at a 98% similarity threshold in QIIME 2, and chimeras were removed. The OTUs were aligned and annotated using a specialized anammox bacterial gene database.
• Correlation Inference: For microbial data, use a robust correlation method that accounts for compositionality, such as SparCC (Sparse Correlations for Compositional Data). Alternatively, Pearson or Spearman correlations can be used on properly transformed data. Calculate all pairwise correlations between OTUs/ASVs.
• Significance Filtering: Apply a p-value correction (e.g., Benjamini-Hochberg False Discovery Rate) to control for multiple comparisons. Retain only correlations that are statistically significant (e.g., p-adjusted < 0.01) and exceed a minimum correlation strength threshold (e.g., |r| > 0.6). This filtered set of correlations forms the adjacency matrix for the network.
• Network Construction and Analysis: Import the adjacency matrix into a network analysis tool like igraph (in R) or NetworkX (in Python). Calculate all network-level (e.g., modularity) and node-level metrics (e.g., degree, betweenness centrality, Zi, Pi).

The following diagram outlines the logical relationships and decision process for classifying node roles based on the within-module connectivity (Zi) and among-module connectivity (Pi) metrics, which is the final, crucial step for pinpointing keystone hubs.

The Scientist's Toolkit: Essential Reagents and Materials

Table 3: Research Reagent Solutions for Anammox Network Analysis

Item	Function / Application	Example Product / Kit
Soil DNA Extraction Kit	To extract high-quality, inhibitor-free microbial DNA from complex sediment samples for downstream PCR.	FastDNA SPIN Kit for Soil (MP Biomedical) [24]; PowerSoil DNA Isolation Kit (Mo Bio) [39] [38]
Anammox-Specific PCR Primers	To specifically amplify the 16S rRNA gene fragment of anammox bacteria, reducing background noise from other community members.	Brod541F & Amx820R [24]
PCR Reagents	To amplify the target gene region for sequencing. Includes polymerase, buffers, dNTPs, and BSA to improve amplification from environmental samples.	Premix Ex Taq (Takara Bio) [24]
High-Through Sequencing Platform	To generate the raw sequence data required for community profiling and subsequent network analysis.	Illumina MiSeq [39]
Bioinformatics Software Suite	For processing raw sequences, including quality control, OTU/ASV picking, taxonomy assignment, and diversity analysis.	QIIME 2 [24]
Network Analysis Tools	To construct, visualize, and analyze co-occurrence networks and calculate key topological metrics.	igraph (R), NetworkX (Python), Gephi, Cytoscape
Statistical Programming Environment	To perform correlation calculations, statistical filtering, and integrate all analytical steps.	R [24], Python
Stearoyl coenzyme A lithium	Stearoyl coenzyme A lithium, CAS:193402-48-1, MF:C39H66Li4N7O17P3S, MW:1057.73	Chemical Reagent
Dabcyl-RGVVNASSRLA-Edans	Dabcyl-RGVVNASSRLA-Edans, CAS:163265-38-1, MF:C73H109N23O18S, MW:1628.9 g/mol	Chemical Reagent

Co-occurrence network analysis provides a powerful computational framework to move beyond microbial taxonomy and understand the functional ecology of complex communities. By applying this methodology, researchers can pinpoint keystone hubs—highly connected taxa that are critical for maintaining community structure and function. In anammox bacterial communities, this approach has revealed that genera like Candidatus Scalindua and various rare species act as such keystones, orchestrating the community's response to environmental gradients and sustaining the critical process of nitrogen loss [24]. As these methods become more sophisticated and integrated with meta-omics data, they will profoundly enhance our ability to predict the behavior of microbial ecosystems and harness their functions in biotechnology and environmental management.

Linking Keystone Species to Nitrogen Conversion Rates and Functional Gene Abundance

In the engineered ecosystems of anaerobic ammonium oxidation (anammox) reactors, keystone species are not merely participants but fundamental architects of community stability and function. These highly connected taxa orchestrate microbial interactions, driving the assembly of the community and exerting disproportionate influence on ecosystem processes, including the critical function of nitrogen removal [40]. The stability of the entire anammox microbiome is intrinsically linked to its biodiversity, with higher phylogenetic diversity fostering greater resistance to perturbation [40]. Within this complex consortium, which includes anammox bacteria, ammonia-oxidizing bacteria (AOB), nitrite-oxidizing bacteria (NOB), Chloroflexi (CFX), and heterotrophic denitrification bacteria (HDB), keystone taxa embed specialized metabolic functions that are essential for maintaining system performance [4] [40]. This technical guide explores the definitive links between these keystone species, quantitative nitrogen conversion rates, and the abundance of key functional genes, providing a framework for optimizing anammox processes for wastewater treatment.

Keystone Taxa and Their Identified Roles in Anammox Consortia

The anammox core is a complex microbial community where keystone taxa sustain system function through specialized metabolic roles and interspecies interactions. The removal of these keystone taxa can trigger dramatic shifts in microbiome composition and function [40]. Molecular techniques such as co-occurrence network analysis and machine learning classification are powerful tools for identifying these pivotal organisms [40].

Table 1: Identified Keystone Taxa and Their Functional Roles in Anammox Systems

Keystone Taxon/Group	Proposed Functional Role	Impact on Community/Process
Chloroflexi & Proteobacteria	Mitigation of apoptotic cell tissue effects; Provision of folate and molybdopterin cofactor for AnAOB central carbon metabolism [29].	Alternate in importance over time to maintain functional stability; reinforce interspecies interactions to resist nitrogen-loading fluctuations [29] [6].
Ca. Kuenenia (AnAOB)	Directly catalyzes the anammox reaction; potential production of acetate and glycogen to enhance microbial interactions and biofilm formation [4].	Higher relative abundance correlates with increased biomass of biofilm and flocs, and enhanced nitrogen removal efficiency (NRE) [4].
Nitrospira & Gemmatimonas	Carry out specialized "keystone functions" such as nitrogen metabolism and phosphonate/phosphinate metabolism [40].	Their specialized metabolic functions are essential for soil microbiome stability, with implications for analogous functions in engineered anammox systems.
Heterotrophic Denitrification Bacteria (HDB)	Perform the "nitrite loop" – reducing nitrate (a by-product of the anammox reaction) back to nitrite [4].	Increases the overall Nitrogen Removal Efficiency (NRE) by recycling a limiting substrate, as the theoretical nitrate production by anammox is 11% [4].

The dynamics of these keystone species are not static. Research shows that in a stable anammox reactor, the subnetworks of anammox bacteria and the overall microbial network can vary significantly over time, with Chloroflexi and Proteobacteria alternately playing important roles in maintaining stability [29]. Furthermore, microbial communities adopt modular collaboration to counteract environmental stress, such as fluctuations in nitrogen loading, with increased inter-module connectivity reinforcing interspecies interactions [6].

Quantitative Links: Nitrogen Conversion and Functional Genes

The influence of keystone taxa extends directly to measurable process performance and the genetic potential of the community. The abundance of anammox bacteria, often acting as keystone organisms, does not always directly reflect reactor performance; the microbial interactions and competition within the consortium are equally critical [4].

Nitrogen Loading and Conversion Rates

The Nitrogen Loading Rate (NLR) has a demonstrable, nonlinear relationship with system performance and anammox bacterial abundance. A study investigating anammox reactors treating rare earth tailings leachate revealed a clear optimal range for NLR [6].

Table 2: Response of Anammox System to Nitrogen-Loading Rate (NLR) Fluctuations [6]

Operational Phase	Nitrogen-Loading Rate (NLR)	Anammox Bacterial Abundance	System Performance Observation
Phase P1 (Baseline)	1.38 ± 0.01 kg/m³·d	5.85%	Baseline stable performance.
Phase P3 (Synchronous Enhancement)	Increasing (up to 3.68 kg/m³·d)	Increased to 11.43%	Denitrification efficiency and anammox bacterial abundance showed synchronous enhancement.
Phase P2/P4 (Inhibition/Starvation)	>3.68 kg/m³·d (Inhibition) or Nitrogen Starvation	Reduced	Performance deterioration and reduced anammox bacterial abundance.

Functional Gene Abundance

The presence and activity of keystone taxa are reflected in the abundance of functional genes responsible for nitrogen transformation and associated metabolic pathways. Metagenomic and proteomic analyses demonstrate that keystone taxa drive the abundance of critical genes.

Table 3: Key Functional Genes and Metabolic Pathways in Anammox Consortia

Functional Gene / Pathway	Function / Encoded Enzyme	Association with Keystone Taxa/Consortia
hzs	Hydrazine synthase; catalyzes the combination of nitric oxide and ammonium to form hydrazine, a key anammox intermediate [9]	Higher abundance in reactors with higher Ca. Kuenenia [4].
nir	Nitrite reductase; reduces nitrite to nitric oxide [9]	Higher abundance in reactors with higher Ca. Kuenenia [4].
Nitrogen Metabolism	Specialized metabolic pathway for nitrogen transformation	Identified as a keystone function carried out by specific taxa like Nitrospira [40].
Carbon Metabolism (fdh, glgA/B/C, acs)	Formate dehydrogenase, glycogen synthesis, Acetyl-CoA synthetase; related to carbon utilization and storage	Higher in consortia with robust keystone interactions, suggesting Ca. Kuenenia may produce acetate and glycogen to enhance microbial interactions [4].

Experimental Protocols for Investigating Keystone Species

Reactor Operation and Performance Monitoring

A standard methodology for investigating anammox keystone species involves operating laboratory-scale bioreactors, such as Up-flow Anaerobic Sludge Bed (UASB) or Expanded Granular Sludge Bed (EGSB) reactors, under controlled conditions [4] [6].

Detailed Protocol:

• Reactor Setup: Use a UASB reactor with a working volume of 1.0-10 L, filled with sponge or other biomass carriers. Maintain temperature at room temperature (~20°C) or 35±1°C using a water jacket or insulating foam [4] [29] [6].
• Inoculation: Seed the reactor with mature anammox sludge. The initial inoculated abundance of anammox bacteria can be a key variable, for example, comparing high (>20%) and low (<8%) relative abundances of Candidatus Kuenenia [4].
• Feeding: Feed the reactor with synthetic wastewater. A typical composition includes [6]:
  • (NHâ‚„)â‚‚SOâ‚„ and NaNOâ‚‚ as nitrogen sources (controlled at a stoichiometric ratio of NH₄⁺:NO₂⁻ ≈ 1:1.32)
  • NaHCO₃ and KHCO₃ as buffer and inorganic carbon source
  • Essential minerals: KHâ‚‚POâ‚„, MgSO₄·7Hâ‚‚O, CaCl₂·2Hâ‚‚O
  • Trace element solutions I and II (containing EDTA, FeSO₄·7Hâ‚‚O, and other micronutrients like Mo, Ni, Cu, Co, Zn, Mn)
• Operation: Maintain hydraulic retention time (HRT) and nitrogen loading rate (NLR) as per experimental design. Strictly control dissolved oxygen (DO) below 0.1 mg·L⁻¹ to preserve anoxic conditions [4].
• Performance Monitoring: Regularly monitor influent and effluent concentrations of NH₄⁺-N, NO₂⁻-N, and NO₃⁻-N. Calculate performance metrics: Nitrogen Removal Efficiency (NRE), NLR, and Nitrogen Removal Rate (NRR) [6].

Microbial Community and Network Analysis

This protocol outlines the steps for resolving the microbial community structure and identifying keystone taxa.

Detailed Protocol:

• DNA Extraction: Collect sludge samples (flocs and biofilm) from the reactor at different time points. Extract total genomic DNA using a commercial kit (e.g., PowerSoil DNA Isolation Kit) [40].
• Sequencing: Perform high-throughput sequencing of the 16S rRNA gene to profile the microbial community. For functional insights, shotgun metagenomic sequencing can be employed [40].
• Bioinformatic Processing: Process raw sequencing data using platforms like MOTHUR or QIIME2 to obtain operational taxonomic units (OTUs) or amplicon sequence variants (ASVs) [40].
• Co-occurrence Network Construction: Construct microbial association networks using statistical tools (e.g., in R). Calculate correlation matrices (e.g., SparCC) and build networks in programs like Gephi or using R packages. Define connections based on strong correlations and statistical significance [29] [40].
• Keystone Taxon Identification: Analyze network topology to identify keystone taxa based on metrics such as [29] [40]:
  • Within-Module Connectivity (Zi): Measures how well a node connects to others within its own module.
  • Among-Module Connectivity (Pi): Measures how well a node connects to different modules. Nodes with high Zi and Pi are considered network hubs, while those with low Zi but high Pi are connectors; both can be classified as keystone taxa.
• Functional Prediction: Use tools like PICRUSt2 to predict the functional composition of the metagenome from 16S rRNA data and annotate pathways against databases like KEGG [6] [40].

The Scientist's Toolkit: Essential Research Reagents and Materials

Successful investigation into anammox keystone species requires a suite of specific reagents and laboratory materials.

Table 4: Essential Research Reagents and Materials for Anammox Keystone Species Studies

Category	Item	Function / Application
Bioreactor Components	UASB/EGSB Reactor vessel, peristaltic pumps, sponge/plastic carriers, temperature control system	Creates the engineered ecosystem for anammox biomass cultivation and study under controlled conditions [4] [6].
Culture Media Components	(NH₄)₂SO₄, NaNO₂, NaHCO₃, KH₂PO₄, MgSO₄·7H₂O, CaCl₂, Trace element solutions (EDTA, FeSO₄, ZnSO₄, etc.)	Constitutes the synthetic wastewater, providing essential substrates, buffers, minerals, and micronutrients for anammox consortia [6].
DNA Analysis Kits	PowerSoil DNA Isolation Kit (or equivalent)	Standardized and efficient extraction of high-quality genomic DNA from complex sludge samples [40].
Sequencing & Bioinformatics	16S rRNA gene primers, Shotgun metagenomic library prep kits, MOTHUR/QIIME2 software, PICRUSt2 tool, R with relevant packages (e.g., igraph, phyloseq)	For profiling microbial community structure, predicting functional potential, and constructing co-occurrence networks [40].
Analytical Instruments	Spectrophotometer/Flow Analyzer, pH meter, Scale	For routine monitoring of nitrogen species concentrations (NH₄⁺-N, NO₂⁻-N, NO₃⁻-N) and other water quality parameters [4] [6].
Diphenylsulfane-d1	Diphenylsulfane-d1, CAS:180802-01-1, MF:C12H10S, MW:196.34 g/mol	Chemical Reagent
D-Alanine-d3	D-Alanine-d3, MF:C3H7NO2, MW:92.11 g/mol	Chemical Reagent

The intricate relationship between keystone species, nitrogen conversion rates, and functional gene abundance forms the bedrock of a stable and efficient anammox ecosystem. Keystone taxa such as Ca. Kuenenia, Chloroflexi, and Proteobacteria do not operate in isolation but are embedded within a complex network where their specialized metabolic functions—from direct nitrogen removal to cross-feeding and structural support—sustain the entire community's function [4] [29] [40]. The quantitative links, evidenced by the synchronized enhancement of anammox abundance and denitrification efficiency with increasing NLR (up to a threshold) and the upregulation of genes like hzs and nir in robust consortia, provide a clear roadmap for diagnostics and optimization [4] [6]. A deep understanding of these relationships, enabled by the detailed experimental and analytical protocols outlined in this guide, empowers researchers to steer anammox communities toward greater resilience and performance, advancing the application of this critical technology in sustainable wastewater treatment.

Carrier Type Selection for Enrichment and Retention of Keystone Anammox Bacteria

The anaerobic ammonium oxidation (anammox) process is universally regarded as a robust biological nitrogen removal approach with significant energy-saving potential [35]. The anammox bacterial community, or the "anammox core," is a complex microbial consortium primarily featuring anammox bacteria but also including ammonia-oxidizing bacteria (AOB), nitrite-oxidizing bacteria (NOB), Chloroflexi, and heterotrophic denitrifying bacteria (HDB) [4]. Among these, keystone anammox species play a disproportionately large role in maintaining community structure and function. However, their slow growth rates, high sensitivity to environmental fluctuations, and propensity for loss in effluent present major challenges for widespread application [35] [41].

Biofilm-based reactors provide a promising solution by offering a favorable microenvironment for microbial growth, safeguarding against extreme conditions, and ensuring longer sludge retention times [35]. The selection of an appropriate carrier type is critically important for the robust and efficient nitrogen removal performance of anammox systems in engineering applications, as it directly influences microbial recruitment, keystone species interactions, and metabolic functionality [35]. This technical guide synthesizes current research to provide a comprehensive framework for selecting carriers that optimize the enrichment and retention of keystone anammox bacteria, framed within the context of keystone species research in anammox bacterial communities.

Carrier Type Induces Differential Biofilm Morphology and Activity

The physical and chemical characteristics of carriers significantly influence biofilm development patterns, which in turn affect microbial activity and reactor performance. Different carrier types impose distinct selection pressures on microbial recruitment by potentially molding microhabitat differences [35].

Biofilm Growth Patterns and Mass Transfer Characteristics

Carrier morphology determines the development patterns of biofilm. Research has demonstrated that anammox biomass exhibits different growth states on various carriers:

• Planar Growth Pattern: Observed on carriers like vinylon fabrics of cord-like polyvinyl formal (CP) and polyethylene ring plastic (PHC), where biomass forms clusters of different sizes [35].
• Three-Dimensional Growth State: Found in elastic cosine sponge (ECS) carriers, creating a more complex biofilm architecture [35].

The structural differences between macroporous and microporous carriers lead to significant variations in mature biofilm characteristics. Macroporous carriers (e.g., plastic rings) typically develop thin biofilms with planar growth, while microporous carriers (e.g., polyurethane sponges) form thick biofilms (>5 mm) with three-dimensional growth patterns [42]. These structural differences profoundly impact mass transfer and microbial distribution.

Table 1: Biofilm Characteristics and Nitrogen Removal Performance on Different Carrier Types

Carrier Type	Physical Structure	Biofilm Morphology	Dominant Anammox Genera	Nitrogen Removal Rate	Key Advantages	Implementation Challenges
Elastic Cosine Sponge (ECS)	Microporous, 3D structure	Thick biofilm (>5 mm), 3D growth	Candidatus Brocadia, Candidatus Jettenia	0.88 kg N/m³·d [35]	Rapid early-phase AnAOB enrichment, high biomass retention	Potential clogging, substrate diffusion limitation at high biomass
Polyurethane Sponge	Microporous, reticulated structure	Thick biofilm, suspended biomass	Candidatus Brocadia	Higher specific activity [42]	Synergy between anammox and denitrifying bacteria	Requires controlled biomass growth to maintain efficiency
Plastic Ring (PHC)	Macroporous, planar surface	Thin biofilm, planar growth	Candidatus Kuenenia	Efficient at lower loads [35]	Reduced clogging risk, stable performance	Lower initial enrichment rate
Cord-like Polyvinyl (CP)	Hybrid structure	Planar growth, clustered biomass	Varies with operational conditions	Moderate [35]	Balanced performance	Complex biofilm control
Nonwoven Fibrous Material (ETEX)	Fibrous, high surface area	Rapid biomass accumulation	Community-dependent	1.5-3x higher per chip [43]	Rapid start-up, efficient retention	Buoyancy management required

Biomass Retention and Carrier Clogging Considerations

Carrier selection must balance biomass retention capacity with long-term operational stability. Microporous carriers like sponges demonstrate superior biomass retention capabilities but face potential clogging issues with excess biomass, leading to poor substrate diffusion and decreased removal efficiency [35]. Macroporous carriers generally exhibit better long-term stability with reduced clogging potential but may have lower initial enrichment rates [42].

Experimental evidence suggests that ECS carriers, despite clogging risks, are easier for enriching and retaining anaerobic ammonium oxidation bacteria (AnAOB) compared to CP and PHC carriers during the early phase of anammox biofilm formation, making them preferable for reactor start-up [35]. The nonwoven fibrous material used in ETEK biochips showed a five-fold more rapid accumulation of activated sludge biomass upon reactor launching compared to foamed polyethylene carriers [43].

Carrier Effects on Microbial Community Structure and Keystone Species

Carrier type exerts significant influence on anammox community assembly by shaping the relative abundance of keystone species and their symbiotic partners. Understanding these ecological dynamics is essential for targeted enrichment of desired anammox consortia.

Selection Pressure on Microbial Recruitment

Different carrier surfaces create distinct microniches that impose selection pressure on microbial recruitment. Plastic surfaces can function as either net autotrophic "hot spots" or organic carbon pools, directly influencing the preference of microbes toward carriers [35]. This selection pressure significantly affects the relative abundance of keystone anammox genera:

• Carrier-Specific Genus Enrichment: Research has demonstrated that carrier type can significantly alter the dominant anammox species. For instance, the elastic cosine sponge (ECS) carrier promoted the enrichment of Candidatus Brocadia, while Candidatus Kuenenia became more dominant in suspended sludge within the same reactor [35].
• Rare Species Importance: Rare anammox species play critical roles in maintaining ecological stability in coastal sediments, with distinct distribution patterns observed across different environments [11]. Carrier selection should consider the preservation of rare taxa that may contribute to community resilience.

Microbial Interactions and Cross-Feeding Mechanisms

Carrier type affects keystone species interactions by potentially molding microhabitat differences [35]. Complex microbial trophic networks develop within biofilms, with different functional microbes establishing close relationships:

• Syntrophic Partnerships: Mature anammox biofilms typically feature anammox bacteria aggregated in the inner layer, while AOB predominantly occupy the outer layer, creating favorable conditions for stable partial nitritation/anammox (PN/A) processes [42].
• Cross-Feeding Relationships: Metagenome-assembled genomes have revealed that dominant denitrifiers can provide various essential materials such as amino acids, cofactors, and vitamins for anammox bacteria [5] [31]. This cross-feeding highlights the importance of microbial interactions for enhancing anammox nitrogen removal.
• Ecological Cooperation: The ecological cooperation between anammox and denitrifying bacteria tends to increase microbial community stability, indicating potential coupling between these functional groups [5] [31]. The nirS-type denitrifiers show stronger coupling with anammox bacteria than nirK-type denitrifiers during enrichment [5].

Experimental Protocols for Evaluating Carrier Performance

Rigorous experimental methodologies are essential for quantitatively assessing carrier performance in anammox systems. The following protocols provide standardized approaches for comparing carrier efficacy.

Carrier Selection and Reactor Setup

Materials and Reagents:

• Candidate carriers (e.g., ECS, polyurethane sponge, plastic rings, nonwoven materials)
• Inoculum: Mature anammox sludge (VSS concentration ~4920 ± 148 mg/L)
• Synthetic wastewater components: NH₄⁺-N and NO₂⁻-N sources, mineral supplements
• Bioreactor systems with temperature and mixing control
• Analytical instruments: HPLC, IC, microelectrodes, molecular biology tools

Procedure:

• Carrier Preparation: Cut carriers to uniform sizes (e.g., 20mm cubes for sponges). Pre-treat if necessary to modify surface properties.
• Reactor Inoculation: Seed reactors with mature anammox sludge at appropriate volatile suspended solid (VSS) concentrations [35].
• Carrier Addition: Add carriers to reactors at specified fill ratios (typically 10-15% of reactor working volume) [42] [43].
• System Operation: Operate reactors in up-flow mode with controlled hydraulic retention time (HRT). Maintain temperature at 32-34°C and dissolved oxygen below 0.1 mg·L⁻¹ [35] [4].
• Performance Monitoring: Regularly monitor nitrogen concentrations in influent and effluent, biomass accumulation, and microbial community dynamics.

Biofilm Characterization and Microbial Community Analysis

Protocol for Biofilm Morphology Analysis:

• Sampling: Extract carrier samples at regular intervals throughout operation.
• Microscopy: Examine biofilm structure using scanning electron microscopy (SEM) and optical microscopy [42].
• Biomass Quantification: Measure volatile suspended solids (VSS) attached to carriers versus in suspension [35].
• Activity Assessment: Conduct specific anammox activity (SAA) tests using batch assays [35].
• EPS Analysis: Quantify extracellular polymeric substance production, including polysaccharides and proteins [35].

Molecular Analysis Workflow:

• DNA Extraction: Extract total DNA from biofilm samples using commercial kits.
• Quantitative PCR: Target functional genes (hzsB, nirS, nirK) to quantify anammox and denitrifying bacteria [5].
• High-Throughput Sequencing: Perform 16S rRNA gene amplicon sequencing to characterize microbial community structure [11].
• Metagenomic Analysis: Employ shotgun metagenomics for functional potential assessment [5] [31].
• Network Analysis: Construct co-occurrence networks to identify keystone species and microbial interactions [11].

The Scientist's Toolkit: Essential Research Reagents and Materials

Table 2: Essential Research Materials for Anammox Carrier Studies

Category	Specific Items	Function/Application	Example Sources/References
Carrier Types	Elastic cosine sponge (ECS), Polyurethane sponge, Polyethylene rings, Nonwoven fibrous materials	Provide attachment surfaces for biofilm formation; create distinct microhabitats	[35] [42] [43]
Analytical Instruments	Scanning Electron Microscope (SEM), Nutrient Auto-analyzer, Microelectrodes (O₂, NO₂⁻), Qubit Fluorometer, NanoDrop Spectrophotometer	Biofilm visualization; nitrogen species quantification; microgradient measurement; DNA quantification	[42] [11]
Molecular Biology Reagents	FastDNA SPIN Kit for soil, Premix Ex Taq, Brod541F/Amx820R primers, hzsB gene primers, nirS/nirK gene primers	DNA extraction; PCR amplification; anammox bacteria quantification; denitrifier detection	[5] [11]
Synthetic Wastewater Components	(NH₄)₂SO₄, NaNO₂, NaHCO₃, KH₂PO₄, CaCl₂·2H₂O, MgCl₂, FeSO₄·7H₂O, Trace element solutions	Provide controlled nutrient environment for anammox growth; maintain optimal pH and mineral balance	[35] [43]
Bioreactor Systems	Up-flow anaerobic sludge blanket (UASB) reactors, Moving bed biofilm reactors (MBBR), Sequencing batch reactors	Maintain controlled hydraulic retention time; temperature regulation; mixing provision	[4] [43]
Pristane-d40	Pristane-d40, CAS:16416-35-6, MF:C19H40, MW:308.8 g/mol	Chemical Reagent	Bench Chemicals
8'-Oxo-6-hydroxydihydrophaseic acid	8'-Oxo-6-hydroxydihydrophaseic acid, MF:C15H20O7, MW:312.31 g/mol	Chemical Reagent	Bench Chemicals

Operational Parameters and Implementation Strategies

Successful implementation of carrier-based anammox systems requires careful consideration of multiple operational parameters that influence keystone species enrichment and retention.

Optimal Operational Conditions for Different Carrier Types

Table 3: Operational Parameters for Carrier-Based Anammox Systems

Parameter	Optimal Range	Impact on Keystone Species	Carrier-Specific Considerations
Temperature	32-34°C [5] [43]	Influences growth rates of different anammox species; affects metabolic activity	Microporous carriers provide better insulation against temperature fluctuations
Hydraulic Retention Time (HRT)	5.56 ± 0.2 h to 24 h [4] [43]	Determines contact time between substrates and biomass; affects shear forces	Thick biofilms on sponges tolerate shorter HRTs; thin biofilms require longer contact times
Dissolved Oxygen (DO)	<0.1 mg·L⁻¹ [4]	Critical for anammox bacteria protection; regulates AOB activity for NO₂⁻ supply	Structured carriers create oxygen gradients enabling aerobic/anoxic zone coexistence
pH	7.5-8.5 [43]	Affects enzyme activity and substrate availability (free ammonia)	Biofilms buffer against pH fluctuations better than suspended systems
Nitrogen Loading Rate	0.15-2.7 kg N/(m³·d) [42]	High loads may favor specific species (e.g., Ca. Brocadia sinica at >270 mg NH₄⁺-N/L)	Microporous carriers generally support higher loading rates due to greater biomass retention
Carrier Fill Ratio	10-15% of reactor volume [42] [43]	Determines total available surface area for biofilm development	Optimal ratio balances attachment surface with mixing efficiency and reactor volume

Species Selection and Niche Differentiation

Operational parameters can be manipulated to selectively enrich specific anammox species with desirable characteristics. Different anammox bacteria exhibit distinct ecological niche preferences:

• Substrate Concentration Selection: High ammonium and nitrite concentrations (>270 mg NH₄⁺-N/L and >340 mg NO₂⁻-N/L) can shift the dominant population from Candidatus Brocadia caroliniensis to Candidatus Brocadia sinica [44].
• Organic Carbon Effects: The presence of complex organic carbon can lead to the emergence of heterotrophic communities that coevolve with Candidatus B. caroliniensis, reaching up to 50% relative abundance [44].
• Nitric Oxide Influence: Supplementation with nitric oxide can inhibit heterotrophic growth and promote the coexistence of Candidatus Jettenia with Candidatus B. caroliniensis [44].
• Morphological Preferences: Different anammox species show distinct attachment preferences—Candidatus B. caroliniensis and Candidatus Jettenia preferentially form biofilms on surfaces, whereas Candidatus Brocadia sinica tends to form granules in suspension [44].

Carrier type selection fundamentally influences the enrichment and retention of keystone anammox bacteria through multiple mechanisms: by shaping biofilm morphology and mass transfer characteristics, imposing selection pressure on microbial community assembly, modulating keystone species interactions, and determining functional gene expression related to nitrogen transformation. Microporous carriers like sponges and nonwoven materials typically support thicker biofilms with higher biomass retention and nitrogen removal rates, making them ideal for rapid start-up phases, though they may require management of clogging risks. Macroporous carriers offer more stable long-term operation with reduced clogging potential but may exhibit slower initial enrichment.

The optimal carrier choice must align with specific operational objectives, whether for mainstream or sidestream applications, and consider the targeted anammox species based on their morphological preferences and metabolic capabilities. Future research should focus on developing novel carrier materials with tailored surface properties and spatial architectures that selectively enhance the retention of keystone anammox species while promoting beneficial microbial partnerships. Such advanced carrier designs, combined with optimized operational strategies, will accelerate the implementation of robust anammox processes for sustainable wastewater treatment.

In the complex ecosystem of anammox (anaerobic ammonium oxidation) communities, certain microbial species play disproportionately large roles in maintaining community structure and function. These keystone species engage in intricate networks of metabolic cross-feeding, exchanging essential metabolites and antioxidants that determine the stability and efficiency of the entire consortium. The anammox process represents a resource-efficient biological wastewater treatment technology that converts ammonium (NH₄⁺-N) and nitrite (NO₂⁻-N) into nitrogen gas (N₂) [45]. However, anammox bacteria rarely function in isolation; they exist within complex microbial networks where cross-feeding interactions significantly enhance community resistance to ecological disturbances and improve overall community productivity [45] [31].

Understanding these cooperative relationships is crucial for optimizing anammox processes in both engineered and natural systems. Metabolic cross-feeding refers to the process where bacteria exchange metabolites with other microorganisms, which are essential for growth and metabolism while simultaneously expanding the ecological niche of the cross-feeders [45]. In anammox consortia, these interactions have been observed through the exchange of amino acids, cofactors, vitamins, and even protective antioxidants that enable community survival under stressful conditions [45] [46]. This guide explores the mechanistic basis of these interactions, with particular emphasis on how keystone species orchestrate community function through targeted metabolite exchange.

Metabolic Cross-Feeding Mechanisms in Anammox Communities

Vitamin and Antioxidant Cross-Feeding

Under high ammonium stress, symbiotic bacteria within anammox consortia provide essential antioxidants to anammox bacteria, notably vitamin B6 and methionine, which enhance the antioxidant capacity of anammox bacteria [45]. This cross-feeding of highly effective antioxidants represents a crucial cooperation mechanism that benefits anammox bacteria resisting high NH₄⁺-N concentrations and variable dissolved oxygen (DO) conditions [45]. The provision of these compounds enables anammox bacteria to mitigate oxidative stress and maintain metabolic functionality under adverse environmental conditions.

Parallel mechanisms have been observed in gut microbiome systems, where dietary antioxidants like ergothioneine are metabolized through cross-feeding pathways. In these systems, Clostridium symbiosum transforms ergothioneine into thiourocanic acid (TUA), which is subsequently utilized by Bacteroides xylanisolvens for enhanced energy production under anaerobic conditions [47] [48]. This conserved cross-feeding mechanism across diverse microbial ecosystems highlights the fundamental importance of antioxidant exchange in maintaining microbial community stability.

Amino Acid and Cofactor Exchange

Anammox bacteria and their symbiotic partners engage in reciprocal exchange of essential nutrients. While symbionts provide critical vitamins and antioxidants, anammox bacteria reciprocate by synthesizing costly amino acids for their partners [45]. However, under high NH₄⁺-N conditions (1785.46 ± 228.5 mg/L), anammox bacteria strategically reduce amino acid supply to symbiotic bacteria to conserve metabolic costs, demonstrating dynamic regulation of cross-feeding in response to environmental stress [45].

Additionally, symbiotic bacteria including Armatimonadetes and Proteobacteria provide essential secondary metabolites such as molybdopterin cofactor (MOCO) and folate to anammox bacteria [46]. These metabolites are indispensable for anammox bacterial growth and activity as they significantly influence carbon fixation and acetyl-CoA production [45] [46]. Acidobacteriota-affiliated bacteria further contribute to anammox activity by synthesizing exopolysaccharides that facilitate consortium aggregation [46].

Nitric Oxide (NO) as a Key Cross-Fed Metabolite

Nitric oxide serves as a crucial interactive mediator among anammox bacteria (AnAOB), ammonia-oxidizing bacteria (AOB), and denitrifying bacteria (DNB) [46]. For typical AOB like Nitrosomonas, NO is released either during hydroxylamine (NH₂OH) oxidation at atmospheric O₂ levels or via NO₂⁻ reduction under O₂-limited conditions [46]. Heterotrophic denitrifying bacteria, such as Chloroflexi and Anaerolineae, also encode nitrite reductase (Nir) that produces NO [46].

Table 1: Nitrogen Removal Performance in Anammox Biofilter

Parameter	Value	Context
Anammox Contribution to N Removal	91.3%	With only 14.4% AnAOB abundance [46]
Nitrogen Removal Rate (NRR)	1432.8 ± 298.69 mg N/L/d	In high NH₄⁺-N reactor (SNAD1) [45]
NRR	214.52 ± 88.55 mg N/L/d	In lower NH₄⁺-N reactor (SNAD2) [45]
NH₄⁺-N Removal Efficiency	Up to 85.92%	During anammox enrichment [31]
NO₂⁻-N Removal Efficiency	Up to 95.34%	During anammox enrichment [31]

As a lipophilic gas molecule, NO can freely diffuse among different bacterial communities [46]. AnAOB can avoid NO₂⁻-N reduction and convert NO directly to dinitrogen gas with ammonium as an electron donor, resulting in no NO₃⁻-N production [46]. This NO-dependent anammox pathway represents a more efficient nitrogen removal route that minimizes nitrate residue.

Quantitative Analysis of Microbial Community Adaptations

Community Structure Shifts Under Ammonium Stress

High ammonium concentrations trigger significant restructuring of anammox microbial communities. Under high NH₄⁺-N conditions (297.95 ± 54.84 mg/L in SNAD1 reactor versus 76.03 ± 34.01 mg/L in SNAD2 reactor), approximately 26.1% of bacterial generalists switch to specialists to increase community stability and functional heterogeneity [45]. This metabolic specialization enhances the community's ability to cope with high NH₄⁺-N conditions while maintaining nitrogen removal performance.

The species richness of microbial communities is higher under high NH₄⁺-N conditions compared to lower concentration environments [45]. Anammox bacteria are predominantly distributed in biofilms (8.42-11.17% abundance) rather than suspended sludges (0.46-0.60%) [45], highlighting the importance of biofilm formation for community aggregation and metabolite exchange under stressful conditions.

Table 2: Microbial Community Response to High NH₄⁺-N Conditions

Adaptation Mechanism	Observed Change	Functional Significance
Generalist-to-Specialist Shift	26.1% of bacteria	Increases community stability and functional heterogeneity [45]
Vitamin B6 Uptake	Enhanced by anammox bacteria	Improves antioxidant capacity against high NH₄⁺-N and variable DO [45]
Amino Acid Supply Reduction	By anammox bacteria	Saves metabolic cost under high NH₄⁺-N stress [45]
V/A-type ATPase Upregulation	In anammox bacteria	Counters cellular alkalization caused by free ammonia [45]
Biofilm Association	8.42-11.17% anammox abundance	Facilitates metabolite exchange and stress resistance [45]

Functional Gene Expression and Metabolic Responses

Metagenomic analyses reveal that anammox bacteria upregulate V/A-type ATPase to strive against cellular alkalization caused by free ammonia [45]. This adaptation is crucial for maintaining intracellular pH homeostasis under high ammonium conditions. Additionally, the upregulation of ABC transporter proteins enhances bacterial resistance to high strength NH₄⁺-N by accelerating membrane transport and promoting metabolite exchange between microorganisms [45].

The expression of genes related to NO generation in AOB and DNB, as well as anammox activation on the NH₄⁺+NO→N₂ pathway, demonstrates how functional gene regulation facilitates cross-feeding interactions [46]. Meta-omics analysis confirms that NO serves as an interactive medium that significantly influences nitrogen transformation pathways in anammox systems.

Experimental Protocols for Studying Cross-Feeding

Anammox Enrichment Bioreactor Setup

To investigate ecological interactions between anammox and denitrifying bacteria, set up anaerobic bioreactors with the following specifications [31]:

• Reactor Configuration: Use 5L anaerobic bioreactors modified with inlet and outlet ports for continuous flow operation
• Microbial Carriers: Implement polyurethane sponge fillers to enhance microbial attachment
• Inoculum: Seed with 1.0 kg of freshwater lake sediments (e.g., from eutrophic Lake Donghu) containing approximately 3.0 g/L VSS
• Temperature Control: Maintain at 34 ± 1°C using temperature controllers
• Anoxic Conditions: Flush continuously with argon gas for 30 minutes at 0.5 LPM flow rate before experiments and cover with tin foil to block light
• Hydraulic Retention Time: Set initially at 48 hours, potentially reducing to 24 hours based on removal efficiency
• Mixing: Provide clockwise water circulation at 60 rpm to enhance microbial contact

Synthetic Wastewater Composition

Prepare synthetic wastewater with the following composition [46] [31]:

• Primary Substrates: NHâ‚„Cl (0.76 g/L providing NH₄⁺-N) and NaNOâ‚‚ (1.38 g/L providing NO₂⁻-N)
• Essential Minerals:
  • CaCl₂·2Hâ‚‚O: 0.135 g/L
  • KHâ‚‚POâ‚„: 0.027 g/L
  • FeSO₄·7Hâ‚‚O: 9.0 mg/L
  • MgClâ‚‚: 0.26 g/L
• Trace Elements: Add according to standard anammox cultivation recipes

Metagenomic Analysis and Metagenome-Assembled Genomes (MAGs)

For comprehensive analysis of cross-feeding interactions [45] [31]:

• Sequencing: Generate approximately 345 million reads using high-throughput metagenomic sequencing
• Quality Control: Process reads through quality filtering (approximately 340.5 million reads retained after QC)
• MAG Construction: Assemble high-quality MAGs with completeness >90% and contamination <10%
• Functional Annotation: Annotate genes for nitrogen cycling, metabolite transport, and stress response
• Metabolic Reconstruction: Reconstruct metabolic networks to identify potential cross-feeding pathways

Research Reagent Solutions for Cross-Feeding Studies

Table 3: Essential Research Reagents for Microbial Cross-Feeding Studies

Reagent/Category	Specific Examples	Function/Application
Molecular Biology Kits	DNA/RNA extraction kits	Nucleic acid isolation for metagenomic and transcriptomic analysis [45] [31]
Sequencing Reagents	Illumina sequencing chemistry	High-throughput metagenomic sequencing [45]
Bioinformatics Tools	MAG reconstruction software	Genome binning and functional annotation [45] [31]
Culture Media Components	NHâ‚„Cl, NaNOâ‚‚, mineral supplements	Synthetic wastewater for reactor operation [46] [31]
Inhibition Assay Reagents	Specific metabolic pathway inhibitors	Validation of cross-feeding pathways [46]
Metabolite Standards	Vitamin B6, ergothioneine, NO donors	Metabolite tracing and quantification [45] [47]
Antibody-Based Detection	Enzyme-specific antibodies	Protein expression validation [46]

Visualization of Cross-Feeding Pathways and Interactions

The strategic manipulation of microbial cross-feeding relationships presents promising opportunities for enhancing the performance and stability of anammox systems in wastewater treatment and environmental bioremediation. By identifying and nurturing keystone species within these communities, and understanding their metabolic dependencies, we can design more resilient and efficient microbial consortia. The dynamic nature of these interactions—where microbes adjust their cooperation strategies in response to environmental conditions—provides a blueprint for engineering robust microbial ecosystems capable of withstanding operational stresses and maintaining high nitrogen removal efficiency.

Future research should focus on quantitative modeling of these cross-feeding networks, developing techniques for real-time monitoring of metabolite exchange, and exploring the potential of targeted probiotic amendments to enhance specific cross-feeding pathways. As we deepen our understanding of these complex microbial relationships, we move closer to harnessing the full potential of engineered microbial communities for sustainable environmental biotechnology applications.

Sustaining the Core: Managing Anammox Keystone Species for System Stability

The integration of anaerobic ammonium oxidation (anammox) with denitrification processes represents a paradigm shift in sustainable wastewater treatment, leveraging microbial ecology to achieve efficient nitrogen removal. While the anammox process, mediated by planctomycetal bacteria, provides an autotrophic pathway for converting ammonium and nitrite directly to nitrogen gas, its practical application faces significant challenges including slow microbial growth rates, sensitivity to environmental perturbations, and the accumulation of nitrate byproducts [49]. The coupling of anammox with denitrification processes addresses these limitations by establishing a synergistic relationship between the autotrophic anammox bacteria (AnAOB) and heterotrophic denitrifying bacteria (DNB) [50].

Within these complex microbial consortia, keystone species play a disproportionately large role in maintaining ecosystem stability and function. These critical organisms mediate interactions between AnAOB and DNB, transforming potential competition into cooperation through metabolic cross-feeding and niche differentiation [5] [4]. Understanding the identity, functional roles, and dynamics of these keystone taxa is essential for optimizing system performance, particularly under the low carbon-to-nitrogen (C/N) conditions characteristic of mainstream wastewater treatment [51]. This technical guide examines the keystone dynamics that mitigate competition and enhance stability in anammox-denitrification coupled systems, providing researchers with both theoretical frameworks and practical methodologies for harnessing these ecological interactions.

Keystone Species and Their Metabolic Roles in System Stability

Keystone species in anammox-denitrification systems function as ecological engineers that create favorable conditions for AnAOB persistence and activity. Through specific metabolic functions and interaction networks, these taxa help resolve the inherent competition between AnAOB and DNB for common substrates like nitrite, instead establishing cooperative cross-feeding relationships that enhance overall system performance.

Identification and Functional Characterization of Keystone Taxa

Molecular analyses and network correlation studies have identified several consistent keystone taxa in high-performing anammox-denitrification systems:

• Chloroflexi and Proteobacteria: These phyla demonstrate dynamic and often alternating dominance in stable anammox systems, where they mitigate the effects of apoptotic cell tissue and contribute to sludge granulation [29]. Metagenome-assembled genomes indicate these organisms provide critical metabolic intermediates including amino acids, cofactors, and vitamins that support AnAOB growth and metabolism [5]. Their presence correlates strongly with system stability, with more negative interactions between AnAOB and heterotrophs associated with enhanced functional stability [29].

• Thauera and Afipia: As dominant denitrifying genera, these bacteria exhibit strong cooperative coupling with AnAOB, particularly in systems enriched from natural environments like lake sediments [5]. These nirS-type denitrifiers show stronger ecological coupling with anammox bacteria than nirK-type denitrifiers, facilitating efficient nitrite supply for the anammox process while minimizing carbon competition [5].

• Candidatus Jettenia and Candidatus Brocadia: These AnAOB genera frequently emerge as keystone taxa in coupled systems, with their relative dominance shifting in response to environmental conditions [5] [52]. Under seasonal cooling, Candidatus Brocadia demonstrates remarkable resilience, with both absolute and relative abundances increasing significantly (by 429.1% and 343.5%, respectively) as temperature decreased from 27.8°C to 7.5°C in pilot-scale implementations [52].

Table 1: Key Functional Roles of Keystone Microbial Groups in Anammox-Denitrification Systems

Keystone Group	Phylogenetic Affiliation	Primary Ecological Function	Impact on System Performance
Chloroflexi	Chloroflexi	EPS biosynthesis, sludge aggregation, metabolite provision	Enhances granular stability and AnAOB retention
Proteobacteria	Proteobacteria	Organic matter degradation, vitamin and cofactor synthesis	Supports AnAOB metabolism and reduces inhibition
Thauera	Proteobacteria (Beta)	Partial denitrification (NO₃⁻ to NO₂⁻), SMP utilization	Stable nitrite supply for anammox without carbon competition
Candidatus Brocadia	Planctomycetes	Core anammox metabolism (NH₄⁺ + NO₂⁻ → N₂)	Direct nitrogen removal under mainstream conditions

Metabolic Cross-Feeding and the "Nitrite Loop"

The conceptual "nitrite loop" represents a critical metabolic architecture facilitated by keystone species in coupled systems [4]. In this framework, keystone denitrifiers like Thauera perform partial denitrification, reducing nitrate to nitrite rather than completing the pathway to nitrogen gas. This generated nitrite then becomes available for AnAOB metabolism, creating a sustainable cycle that minimizes the accumulation of nitrate, which represents approximately 11% of the nitrogen transformed by the anammox reaction [4].

Simultaneously, AnAOB and associated keystone species generate soluble microbial products (SMPs) and extracellular polymeric substances (EPS) through metabolism and biomass decay [50]. These microbial byproducts, including tryptophan and humic acid-like compounds, serve as effective carbon sources for keystone denitrifiers, creating a reciprocal cross-feeding relationship that reduces dependence on external organic carbon [50]. This metabolic exchange transforms the potentially competitive relationship between AnAOB and DNB into a cooperative one, enhancing system stability and nitrogen removal efficiency.

Quantitative Performance Metrics of Keystone-Mediated Systems

The impact of robust keystone dynamics manifests in measurable performance enhancements across various system configurations and operating conditions. The following quantitative data illustrate the performance benefits achieved when keystone-mediated interactions are established in anammox-denitrification systems.

Table 2: Nitrogen Removal Performance Across Different Anammox-Denitrification Configurations

System Configuration	Operating Conditions	Nitrogen Removal Efficiency	Key Contributing Factors	Reference
PD/A-FeS@ZVI	Low C/N wastewater	97.06 ± 1.41% TNRE	Fe-S cycling enhanced nitrate reduction	[51]
Pilot-scale PDA	27.8-20.0°C	75.0 ± 4.6% NRE	Partial denitrification supplying NO₂⁻ to anammox	[52]
Pilot-scale PDA	10.0-7.5°C	70.4 ± 4.5% NRE	Anammox contribution to N-remaintained at 39.7 ± 6.7%	[52]
Immobilized filler SAD	Using endogenous SMPs	85.6% TNRE	SMP-mediated cooperation between AnAOB and DNB	[50]
UASB with keystone enrichment	Room temperature, NLR = 0.96 kg-N·m⁻³·d⁻¹	91.21% ARE, 80.13% NRE	Microbial interactions in biofilm and flocs	[4]

The data demonstrate that systems with established keystone dynamics maintain robust nitrogen removal efficiency even under challenging conditions, such as low temperatures and low C/N ratios. The pilot-scale implementation of partial denitrification-anammox (PDA) maintained approximately 70% nitrogen removal efficiency even at temperatures as low as 7.5-10°C, highlighting the remarkable resilience conferred by proper microbial community structure [52]. Notably, the anammox process contributed 39.7 ± 6.7% to the total nitrogen removal even at these low temperatures, underscoring the importance of maintaining viable AnAOB populations through keystone partnerships [52].

Methodological Framework for Analyzing Keystone Dynamics

Establishing and monitoring keystone dynamics in anammox-denitrification systems requires specialized methodological approaches spanning reactor operation, molecular biology, and metabolic analysis. The following protocols provide a standardized framework for investigating these complex ecological interactions.

Reactor Operation and Performance Monitoring

System Configuration: Utilize up-flow anaerobic sludge blanket (UASB) reactors or sequencing batch reactors (SBRs) with integrated biofilm carriers [5] [4]. For immobilized cell systems, employ polyvinyl alcohol (PVA) gel carriers to separately entrap anammox and denitrifying biomass, creating defined interaction zones [50].

Operating Parameters:

• Maintain hydraulic retention time (HRT) at 5.56 ± 0.2 hours for enriched systems [4]
• Control dissolved oxygen (DO) below 0.1 mg·L⁻¹ to support anammox activity [4]
• For temperature stress tests, operate across a range from 7.5°C to 35°C to assess keystone-mediated resilience [52]
• Maintain a nitrogen loading rate (NLR) of 0.96-2.96 kg-N·m⁻³·d⁻¹ depending on system maturity [29] [4]

Performance Monitoring:

• Daily measurement of nitrogen species (NH₄⁺, NO₂⁻, NO₃⁻) using standard methods [50]
• Regular analysis of extracellular polymeric substances (EPS) composition, particularly the protein/polysaccharide (PN/PS) ratio, which significantly influences sludge characteristics and microbial interactions [53]
• Quantification of soluble microbial products (SMPs), especially tryptophan-like compounds, which serve as electron donors for denitrifiers [50]

Molecular Characterization of Keystone Taxa

DNA Extraction and Amplification:

• Extract genomic DNA from both biofilm and suspended biomass fractions using commercial kits [5]
• Amplify target genes using specific primer sets:
  • 16S rRNA genes for total bacterial community analysis [5]
  • hzsB genes for specific quantification of anammox bacteria [52]
  • nirS and nirK genes for denitrifier community profiling [5]

Quantitative and Sequencing Approaches:

• Apply quantitative PCR (qPCR) to track absolute abundances of functional genes over time [52]
• Conduct 16S rRNA amplicon sequencing to determine relative microbial abundances and identify potential keystone taxa through network correlation analysis [29]
• Perform metagenomic sequencing for deeper functional insights, followed by metagenome-assembled genome (MAG) analysis to elucidate metabolic pathways [5]

Network Analysis:

• Construct microbial co-occurrence networks using correlation metrics from abundance data [29]
• Identify keystone taxa based on high betweenness centrality and connectivity measures within the network [4]
• Analyze network topological features (modularity, connectivity) across different operational phases or conditions [29]

Figure 1: Experimental workflow for analyzing keystone dynamics in anammox-denitrification systems, spanning from sample collection through to metabolic modeling.

Metabolic Interaction Analysis

Isotope Tracing Experiments:

• Apply ¹⁵N isotope pairing techniques to differentiate anammox and denitrification contributions to Nâ‚‚ production [54]
• Use ¹⁵NO₃⁻ and ample ¹⁴NH₄⁺ under anoxic incubation conditions [54]
• Calculate anammox contribution based on ²⁹N2 production while accounting for codenitrification effects [54]

Metabolite Profiling:

• Analyze SMP composition using liquid chromatography-mass spectrometry (LC-MS) [50]
• Identify specific utilizable carbon sources like tryptophan that mediate cross-feeding between AnAOB and DNB [50]
• Track metabolite exchanges in defined co-culture systems to validate putative cross-feeding relationships [5]

Visualization of Keystone-Mediated Metabolic Pathways

The stability of anammox-denitrification systems depends heavily on the metabolic interactions facilitated by keystone species. The following diagram illustrates the key metabolic exchanges and transformations that define these cooperative relationships.

Figure 2: Metabolic cross-feeding network in anammox-denitrification systems showing how keystone species mediate resource exchange to mitigate competition.

Essential Research Reagents and Methodological Tools

Investigating keystone dynamics in anammox-denitrification systems requires specialized reagents and methodological tools. The following table summarizes critical resources for studying these complex microbial communities.

Table 3: Essential Research Reagents and Methodological Tools for Keystone Dynamics Investigation

Category	Specific Reagents/Tools	Application Purpose	Technical Notes
Molecular Biology	hzsB gene primers (hzsB1597F/hzsB1857R)	Anammox bacteria quantification	Target hydrazine synthase beta-subunit gene [52]
	nirS/nirK gene primers	Denitrifier community analysis	nirS-type shows stronger anammox coupling [5]
	16S rRNA gene primers (e.g., 515F/806R)	Total community profiling	Network construction and keystone identification [29]
Isotope Tracers	¹⁵N-labeled nitrate (¹⁵NO₃⁻)	Nitrogen pathway differentiation	Calculate process contributions to N₂ production [54]
	¹⁵N-labeled ammonium (¹⁵NH₄⁺)	Anammox activity verification	Confirm anammox pathway in complex systems [54]
Analytical Standards	SMP reference compounds (tryptophan, humic acids)	Metabolite identification and quantification	HPLC/LC-MS calibration for cross-feeding studies [50]
Process Enhancers	FeS/Zero-Valent Iron (ZVI) composites	Electron donor supplementation	Enhance nitrate reduction in low C/N systems [51]
	PVA-based immobilization matrices	Cell entrapment and segregation	Study defined microbial interactions [50]

Keystone species play an indispensable role in mitigating competition and fostering stability in anammox-denitrification coupled systems. Through metabolic cross-feeding, niche differentiation, and the establishment of efficient nutrient loops, these critical organisms transform potentially competitive microbial relationships into cooperative networks that enhance system functionality and resilience. The strategic management of keystone dynamics enables maintained nitrogen removal efficiency under challenging conditions, including low temperatures and low carbon availability, which represents a significant advancement for sustainable wastewater treatment.

Future research should focus on several promising directions. First, developing targeted enrichment strategies for specific keystone taxa could accelerate system startup and enhance stability. Second, exploring the molecular communication mechanisms underlying keystone functions, particularly quorum sensing and response systems, may reveal new opportunities for process control. Third, engineering synthetic microbial consortia that incorporate defined keystone functions could provide more predictable and robust system performance. Finally, integration of multi-omics data through machine learning approaches promises to uncover deeper insights into the complex interaction networks governing these systems. By advancing our understanding and application of keystone dynamics, researchers can further optimize these ecologically-engineered systems for sustainable nitrogen management across diverse operational scenarios.

Adjusting Bacterial Cooperation Under High Ammonium Stress

In the ecological framework of anaerobic ammonium oxidation (anammox) systems, the concept of keystone species extends beyond mere abundance to encompass functional criticality within the microbial network. These consortia represent remarkable examples of mutualistic interdependence, where anammox bacteria (belonging to the Brocadiales order) rely on synergistic relationships with specialized symbiotic partners to maintain ecosystem functioning under environmental stress [55]. Under optimal conditions, these communities maintain a delicate balance of metabolic exchange, but this equilibrium faces significant challenges when ammonium concentrations rise to inhibitory levels commonly encountered in full-scale wastewater treatment plants (WWTPs) treating sludge digester liquor, where levels can exceed 1700 mg/L NH₄⁺-N [55].

The keystone role of anammox bacteria becomes particularly evident under high ammonium stress, where their metabolic capabilities determine the system's functional stability. Recent research reveals that microbial communities employ sophisticated socio-microbiological strategies to counteract ammonium inhibition, including dynamic adjustments to metabolic cross-feeding, modification of community architecture, and activation of specific cellular detoxification mechanisms [55]. Understanding these adaptive responses provides crucial insights for optimizing anammox processes under the fluctuating ammonium loads characteristic of real-world wastewater treatment scenarios, while advancing fundamental knowledge of microbial ecology in engineered systems.

Metabolic Cross-Feeding Adjustments Under Ammonium Stress

The stability of anammox consortia under high ammonium stress depends significantly on the rewiring of metabolic interactions between anammox bacteria and their symbiotic partners. Through metagenomic analyses of full-scale reactors, researchers have documented profound adjustments in cross-feeding relationships that directly impact community resilience.

Table 1: Metabolic Exchange Alterations Under High Ammonium Conditions

Exchanged Metabolite	Function	Change Under High NH₄⁺	Impact on Anammox Bacteria
Vitamin B₆	Antioxidant protection	Increased supply from symbiotic bacteria	Enhanced resistance to oxidative stress from high NH₄⁺ and variable DO
Methionine	Essential amino acid	Increased supply from symbiotic bacteria	Improved cellular maintenance and stress tolerance
Other amino acids	Building blocks for proteins	Reduced supply from anammox bacteria	Metabolic cost-saving strategy for anammox bacteria
Public goods (via DSF signaling)	Quorum-sensing molecules	Significant upregulation (p < 0.05)	Enhanced community coordination and stress response

A particularly crucial protective mechanism involves the enhanced provision of vitamin B₆ by symbiotic bacteria to anammox bacteria under high ammonium conditions [55]. Vitamin B₆ serves as a potent antioxidant, helping anammox bacteria counter the oxidative stress associated with elevated ammonium concentrations and fluctuating dissolved oxygen (DO) levels. This cross-feeding represents a fundamental keystone interaction, where the metabolic output of symbiotic partners directly determines the stress resilience of the functionally critical anammox bacteria.

Concurrently, anammox bacteria implement energy conservation strategies by reducing their metabolic investment in amino acid production for the community [55]. This reallocation of resources enables them to maintain essential anaerobic ammonium oxidation functions while under physiological duress. The net effect is a shift in the economic balance of the microbial ecosystem, where anammox bacteria prioritize self-preservation over community provision when facing ammonium inhibition.

Community Structure Reorganization

Under high ammonium stress, anammox consortia undergo substantial structural reorganization at both taxonomic and functional levels. Research demonstrates that approximately 26.1% of bacterial generalists switch to specialist niches under high NH₄⁺-N conditions (297.95 ± 54.84 mg/L versus 76.03 ± 34.01 mg/L in series-connected full-scale reactors) [55]. This niche differentiation enhances functional heterogeneity and increases community stability when confronting environmental stress.

The competitive dynamics among anammox species also shift under adverse conditions. Different anammox bacteria exhibit genus-specific stress response strategies influenced by variations in their metabolic versatility and intrinsic growth kinetics [8]. For instance, some species demonstrate superior adaptation capabilities through their ability to utilize alternative substrates or engage in dissimilatory nitrate reduction to ammonium (DNRA) when ammonium becomes inhibitory [18]. This functional diversity within the anammox guild provides functional redundancy, a key ecosystem property that enhances reactor resilience to fluctuating ammonium loads.

Rare subcommunity taxa play disproportionately important roles during stress periods, despite their low abundance under stable conditions [33]. Under nitrite stress (a correlated stressor often accompanying high ammonium conditions), rare species can become dominant contributors to community function, exhibiting differential reactive strategies compared to abundant taxa [33]. This demonstrates how cryptic keystone species emerge during perturbation events to maintain ecosystem functioning.

Cellular Defense Mechanisms

At the cellular level, anammox bacteria activate specific molecular defense systems to counteract ammonium toxicity. A critical mechanism involves the upregulation of V/A-type ATPase, which helps anammox bacteria maintain intracellular pH homeostasis against cellular alkalization caused by free ammonia (the neutral form that increases with higher ammonium concentrations and pH) [55]. This ATP-dependent proton pumping is essential for combating the chemiosmotic challenges associated with high ammonium environments.

Additionally, anammox consortia enhance their collective antioxidant defenses through the coordinated action of superoxide dismutase (Sod) and catalase (Cat) activities [56]. While freshwater anammox species typically exhibit limited Sod activity, marine "Ca. Scalindua sp." demonstrates remarkably high oxygen tolerance (IC₅₀ = 18.0 µM, DOmax = 51.6 µM) attributed to its elevated Sod activity (22.6 ± 1.9 U/mg-protein) coupled with moderate Cat activity (1.6 ± 0.7 U/mg-protein) [56]. Although studied in the context of oxygen stress, these enzymatic systems likely contribute to handling reactive oxygen species generated under various stress conditions, including high ammonium.

The expression of genes related to carbon metabolism, two-component systems, and quorum sensing shows significant upregulation (p < 0.05) by 0.4-6.6 folds under organic pollutant stress [57], demonstrating the comprehensive nature of the cellular response to environmental challenges. These regulatory systems enable anammox bacteria to sense environmental changes and coordinate appropriate physiological adjustments at the population level.

Experimental Methodologies for Investigating Bacterial Cooperation

Continuous-Flow Reactor Systems

Long-term stress response studies typically employ continuous-flow reactors with controlled feeding regimes. The following protocol outlines the establishment of such systems for investigating ammonium stress adaptation:

• Reactor Configuration: Use laboratory-scale upflow anaerobic sludge blanket (UASB) reactors or sequencing batch reactors (SBRs) with working volumes typically ranging from 1-8 L [58] [33]. Maintain complete biomass retention through carrier attachment or granular sludge formation.

• Operating Conditions: Maintain temperature at 32-37°C using water jackets or external temperature control systems [18] [43]. Control pH at 7.3-8.3 using NaHCO₃ buffering [18] [43]. Keep dissolved oxygen concentrations below inhibitory thresholds (typically <0.1 mg/L for freshwater species) through inert gas bubbling (e.g., 95% Ar/5% COâ‚‚) [56].

• Feeding Medium Composition: Prepare synthetic wastewater containing (per liter): (NHâ‚„)â‚‚SOâ‚„ (0.094-0.942 g, providing 20-200 mg N-NH₄⁺/L), NaNOâ‚‚ (providing stoichiometric amounts based on ammonium), KHâ‚‚POâ‚„ (44-54 mg), NaHCO₃ (840-2100 mg), and essential trace elements including FeSO₄·7Hâ‚‚O (9.0 mg), EDTA·4Na (5.0 mg), and other micronutrients [56] [43].

• Stress Application: Gradually increase ammonium concentrations in the feed from non-inhibitory levels (≤100 mg N/L) to stress conditions (≥300 mg N/L) over several weeks to monitor adaptive responses [55].

Diagram 1: Experimental stress application workflow

Metagenomic and Metatranscriptomic Analysis

Omics approaches provide powerful tools for deciphering microbial interactions at the molecular level:

• Biomass Sampling and Preservation: Collect biomass samples (10-50 mL) from reactors at different stress phases. Centrifuge at 10,000 × g for 10 minutes at 4°C. Immediately freeze pellets in liquid nitrogen and store at -80°C until DNA/RNA extraction [55].

• Nucleic Acid Extraction: Use commercial kits with bead-beating step optimized for hard-to-lyse anammox bacteria. Validate DNA/RNA quality using spectrophotometry (A260/A280 ≈ 1.8-2.0) and electrophoresis [55].

• Sequencing and Bioinformatics: Perform shotgun metagenomic sequencing on Illumina platforms (≥10 Gb data per sample). For metatranscriptomics, include mRNA enrichment steps. Process raw sequences through quality filtering, assembly, gene prediction, and functional annotation using appropriate databases (KEGG, COG, NCBI NR) [55].

• Metabolic Reconstruction: Identify potential cross-feeding relationships by mapping annotated genes to metabolic pathways. Quantify expression levels of key stress response genes (e.g., V/A-type ATPase, antioxidant enzymes, vitamin biosynthesis pathways) [55].

Metabolic Interaction Validation

Putative interactions identified through omics analyses require experimental validation:

• Stable Isotope Probing: Use ¹³C-labeled compounds (e.g., ¹³COâ‚‚, ¹³C-acetate) to track carbon flow through microbial networks. Identify cross-fed metabolites using nanoSIMS or GC-MS analysis [55].

• Co-culture Experiments: Islete key symbiotic bacteria using dilution-to-extinction methods in defined media. Establish co-cultures with anammox bacteria under controlled conditions with and without stress factors [55].

• Metabolite Supplementation: Test the protective effects of putative cross-fed metabolites (e.g., vitamin B₆, methionine) by supplementing reactors under stress conditions and monitoring performance improvements and stress marker expression [55].

The Scientist's Toolkit: Essential Research Reagents and Materials

Table 2: Key Research Reagents for Investigating Anammox Bacterial Cooperation

Reagent/Material	Function	Application Examples	Technical Considerations
Percoll density gradient medium	Separation of planktonic anammox cells	Obtaining highly enriched cultures (>99.8% purity) for physiological studies [56]	Enables study of intrinsic species characteristics without interference from consortium members
Hydrocyclone separator	Density-based biomass fractionation	Granule size classification; enrichment of anammox biomass [59]	Critical for studying population heterogeneity and functional specialization
Non-woven fibrous carriers	Biomass attachment and retention	Biofilm studies; community assembly under stress [43] [33]	High surface area promotes rapid biomass accumulation (5× faster than conventional carriers)
Formate and folate supplements	Metabolic effectors	Studying metabolic flexibility; stress response modulation [18]	Formate stimulates DNRA; folate supports one-carbon metabolism in anammox bacteria
Inhibitor stocks (NHâ‚„Cl, NaNOâ‚‚)	Stress application	Establishing inhibition kinetics; resilience testing [33]	Prepare fresh solutions; use gradual concentration increases to study adaptation
RNA stabilization buffers	Nucleic acid preservation	Metatranscriptomic time-series during stress response [55]	Critical for capturing instantaneous transcriptional responses to ammonium shock
Heme quantification reagents	Biomass activity assessment	Spectrophotometric measurement of anammox biomass [59]	Correlation between heme content and anammox activity (R² = 0.91)
Methyl heptadecanoate-d33	Methyl heptadecanoate-d33, CAS:1219804-81-5, MF:C18H36O2, MW:317.7 g/mol	Chemical Reagent	Bench Chemicals
CYM 9484	CYM 9484, CAS:1383478-94-1, MF:C27H31N3O3S2, MW:509.7 g/mol	Chemical Reagent	Bench Chemicals

The sophisticated adjustments in bacterial cooperation under high ammonium stress illustrate the remarkable plasticity of anammox consortia, highlighting the dynamic nature of keystone interactions in engineered ecosystems. The observed shifts—from modified metabolic cross-feeding to community restructuring and cellular defense activation—collectively enable functional persistence under inhibitory conditions. From an ecological perspective, these adaptations demonstrate how microbial societies redistribute resources and reorganize to maintain critical ecosystem functions under stress.

For engineering applications, these insights suggest novel bioaugmentation strategies wherein specific symbiotic partners could be harnessed to enhance anammox reactor resilience. The identification of vitamin B₆ and methionine as critical stress-alleviating metabolites points to potential microbial relationship management approaches through targeted metabolite supplementation. Furthermore, understanding the conditions that promote the emergence of stress-resistant rare species could inform reactor operation strategies that maintain functional diversity as a buffer against process upsets.

Future research should focus on precisely quantifying the metabolic costs of these adaptive strategies and their long-term stability, potentially through integrated multi-omics time-series analyses across gradual stress induction scenarios. Such studies would further illuminate the fundamental principles governing keystone interactions in engineered microbial ecosystems while advancing the practical implementation of robust anammox technologies for sustainable wastewater treatment.

The Impact of Carrier Materials on Keystone Species Interactions and Biofilm Development

In the pursuit of sustainable wastewater treatment, the anaerobic ammonium oxidation (anammox) process has emerged as a revolutionary biotechnology for energy-efficient nitrogen removal. The efficacy of this process hinges upon the functionality of slow-growing keystone species, particularly anammox bacteria (AnAOB) such as Candidatus Brocadia, Candidatus Kuenenia, and Candidatus Jettenia, which operate in syntrophy with ammonia-oxidizing bacteria (AOB) and heterotrophic communities. The inherent challenge of retaining these slow-growing microbes within continuous-flow systems has propelled the development and utilization of specialized carrier materials. These materials are not merely passive scaffolds; they actively shape the microbial ecosystem by modulating the physical and chemical microenvironments, thereby influencing keystone species selection, biofilm architecture, and ultimately, process stability and efficiency. This technical guide examines the profound impact of carrier materials on these critical biological interactions, framing the discussion within the broader context of keystone species management in anammox research. It synthesizes current knowledge on carrier design principles, their influence on microbial niche partitioning, and provides detailed methodologies for studying these complex systems.

Theoretical Foundations: Carrier Mechanisms in Microbial Ecology

Dynamics of Stratified Microbial Niche Formation on Carrier Surfaces

Carrier materials facilitate the initial attachment and subsequent development of stratified biofilms, a process critical for the coexistence of aerobes and anaerobes in single-stage systems. The process initiates with the landing campaign, where microorganisms in the wastewater flow penetrate the hydrodynamic boundary layer and approach the carrier surface. The subsequent adhesion and colonization phase is governed by carrier surface properties, including hydrophobicity, zeta potential, and roughness, which dictate the thermodynamic feasibility of bacterial attachment [60]. This initial conditioning film determines the trajectory of biofilm succession.

The maturation of the biofilm culminates in the formation of a stratified microbial structure, which is paramount for the one-stage partial nitritation-anammox (PN/A) process. AOB, such as Nitrosomonas, preferentially colonize the outer, oxygen-rich layers of the biofilm where they catalyze the oxidation of ammonium to nitrite. This consumption of oxygen creates anoxic microniches in the deeper biofilm strata, providing an ideal habitat for the oxygen-sensitive AnAOB [60]. This spatial organization, facilitated by the carrier's three-dimensional structure, allows for the direct transfer of nitrite from AOB to AnAOB, minimizing competition from nitrite-oxidizing bacteria (NOB) and enabling efficient nitrogen removal.

Material Properties Governing Microbial Adhesion and Biofilm Development

The capacity of a carrier to foster biofilm development is intrinsically linked to its physicochemical properties. Key characteristics include:

• Surface Free Energy (SFE): The thermodynamic compatibility between bacterial cells and the carrier surface is a primary determinant of adhesion. A smaller difference in SFE between the bacterium (γbv) and the carrier (γsv) results in a more negative adhesion energy (ΔFadh), thereby promoting attachment [61]. For instance, the SFE of Nitrosomonas europaea is approximately 48.3 mJ m⁻². Modifying polyvinyl alcohol (PVA) biocarriers to match this SFE significantly enhanced AOB adhesion force to 6.3 nN, compared to a mere 1.2 nN on conventional high-density polyethylene (HDPE) [61].
• Surface Morphology and Porosity: Macroporous carriers typically support two-dimensional biofilm growth (0.4–2.3 mm thick), whereas microporous architectures like polyurethane sponges facilitate three-dimensional growth, yielding biofilms exceeding 5 mm in thickness [60]. This internal porosity is crucial for protecting AnAOB from hydraulic shear and environmental perturbations.
• Chemical Composition and Functional Groups: Materials like melamine and nylon, which possess amine groups and higher surface energy (45.7–48.9 mJ m⁻²), demonstrate superior nitrifier enrichment compared to low-energy materials like HDPE and acetal [61]. The presence of essential trace elements in mineral carriers (e.g., quartz, bentonite, zeolite) can further stimulate the activity of specific autotrophic communities [62] [63].

Table 1: Key Carrier Properties and Their Influence on Biofilm Formation

Property	Impact on Biofilm	Experimental Measurement
Surface Free Energy (SFE)	Determines thermodynamic adhesion energy; smaller SFE difference enhances bacterial attachment.	Contact angle goniometry followed by Neumann's equation of state [61].
Surface Roughness	Increases effective surface area and provides sheltered sites for initial cell attachment.	Atomic Force Microscopy (AFM); profilometry [60].
Porosity & Pore Size	Governs biofilm thickness, density, and protection; creates diffusion-limited anoxic microniches.	Mercury intrusion porosimetry; micro-CT scanning [60].
Surface Charge (Zeta Potential)	Influences initial electrostatic interactions between cells and the carrier surface.	Zeta potential analyzers [60].
Hydrophobicity	Mediates the deposition and adhesion of microbial cells, which often exhibit hydrophobic surfaces.	Water contact angle measurement [60].

Quantitative Performance of Selected Carrier Materials

The selection of an appropriate carrier material has a direct and quantifiable impact on the startup time and nitrogen removal performance of anammox systems. Different materials offer distinct advantages.

• Polyurethane (PU) Porous Material: In an up-flow anammox reactor, PU carriers enabled the successful start-up and enrichment of anammox bacteria within 73 days. The system achieved remarkable removal rates of 97.87% for ammonia and 99.96% for nitrite, with the dominant keystone species being Candidatus Brocadia (20.44% relative abundance) [64]. The high porosity and specific surface area of PU provide an ideal habitat for biofilm development and biomass retention.
• Graphene Oxide-Composite (GO@PU): The integration of GO with PU to form a composite carrier further enhances performance. In an A2/O process treating mainstream municipal wastewater, the GO@PU carrier increased the nitrogen removal efficiency (NRE) from 60% to 84% and the nitrogen removal rate (NRR) from 110 to 151 g N/(m³·d), while also shortening the system start-up period by 10 days [65]. GO's large surface area and biocompatibility are believed to promote microbial metabolism and enrichment.
• Mineral Carriers: Materials such as quartz, kaolin clay, bentonite clay, and zeolite have proven effective for in situ permeable barriers in cold (10°C), nitrogen-polluted aquifers. These carriers supported the indigenous anammox community (dominated by Candidatus Scalindua) and, when bioaugmented with a "warm" community (dominated by *Candidatus Kuenenia"), further improved nitrogen removal rates [62] [63]. Their efficacy is linked to surface morphology and the presence of essential trace elements.

Table 2: Performance Comparison of Different Carrier Materials in Anammox Systems

Carrier Material	Reactor Configuration	Key Keystone Species	Key Performance Metric	Start-up/Enhancement Effect
Polyurethane (PU)	Up-flow biofilm reactor	Candidatus Brocadia	NH₄⁺-N removal: 97.87%; NO₂⁻-N removal: 99.96% [64]	Successful start-up and enrichment in 73 d [64]
GO@PU Composite	A2/O process	AnAOB retained	NRE increased from 60% to 84%; NRR from 110 to 151 g N/(m³·d) [65]	Start-up time reduced by 10 d [65]
Modified PVA	MBBR	Nitrosomonas europaea (AOB)	NH₄⁺-N removal: 91% (vs. 76% on HDPE) [61]	60-fold higher AOB adhesion than HDPE [61]
Mineral Carriers	In-situ permeable barrier	Candidatus Scalindua, C. Kuenenia	Enhanced N-removal in cold groundwater (10°C) [62] [63]	Supported indigenous and bioaugmented anammox communities [62]

Experimental Protocols for Carrier Evaluation

Protocol 1: Evaluating Microbial Adhesion Force via Atomic Force Microscopy (AFM)

Objective: To quantitatively measure the adhesion force between target keystone species (e.g., AOB) and a candidate carrier material.

Materials and Reagents:

• Probe Functionalization: Cell culture of target microorganism (e.g., Nitrosomonas europaea), glutaraldehyde (2.5% v/v) or polyethyleneimine (PEI) solution, phosphate-buffered saline (PBS).
• Substrate Preparation: Flat substrates of the candidate carrier material.
• Equipment: Atomic Force Microscope with a liquid probe holder, tipless cantilevers, V-shaped cantilevers for colloid probe method.

Procedure:

• Probe Functionalization: Tipless cantilevers are sterilized and coated with a cell-adhesive polymer like PEI. A concentrated suspension of the target microorganism is deposited onto the cantilever and fixed with a gentle glutaraldehyde solution (or using PEI adsorption) to immobilize cells without compromising viability [61].
• Force Measurement: The cell-functionalized probe is approached onto a flat, sterile sample of the carrier material submerged in a relevant liquid medium (e.g., synthetic wastewater). Upon contact, a dwell time (e.g., 1-5 seconds) is applied to allow for cell-surface interaction. The probe is then retracted.
• Data Acquisition: The force-distance curve during retraction is recorded. The adhesion force is determined from the "pull-off" force observed in the retraction curve. A minimum of 100 force curves should be collected from different locations on the sample to ensure statistical significance.
• Data Analysis: The average adhesion force and standard deviation are calculated. A higher mean adhesion force indicates superior compatibility between the microorganism and the carrier material.

Protocol 2: Assessing Biofilm Formation and Nitrogen Removal in Laboratory Bioreactors

Objective: To evaluate the long-term performance of a carrier in enriching anammox keystone species and facilitating nitrogen removal.

Materials and Reagents:

• Bioreactor: Laboratory-scale reactor (e.g., MBBR, up-flow anaerobic biofilm reactor) with temperature control (e.g., 30-35°C).
• Carrier: Test carrier material, loaded at 40-50% of reactor volume.
• Inoculum and Feed: Anammox-active sludge (e.g., MLSS ~9500 mg/L). Synthetic wastewater containing NHâ‚„Cl (N source, 30-100 mg/L NH₄⁺-N), NaNOâ‚‚ (N source, 40-130 mg/L NO₂⁻-N), and micronutrient solutions (CaClâ‚‚, KHâ‚‚POâ‚„, MgSO₄·7Hâ‚‚O, trace elements) [64].

Procedure:

• Reactor Setup: The reactor is filled with the carrier and inoculated with anammox sludge. It is initially left in batch mode for 24 hours to facilitate initial attachment.
• Continuous Operation: Continuous feeding is initiated with a long hydraulic retention time (HRT, e.g., 8-48 h). The influent substrate concentrations are maintained at a low level initially (e.g., 30 mg/L NH₄⁺-N, 40 mg/L NO₂⁻-N) [64].
• Performance Monitoring: Effluent is sampled regularly (e.g., 2-3 times per week) and analyzed for NH₄⁺-N, NO₂⁻-N, and NO₃⁻-N concentrations using standard methods (e.g., colorimetric methods). pH and dissolved oxygen (DO) are also monitored.
• Enrichment Phase: Once stable removal is observed, the nitrogen loading rate is increased by either raising influent concentration or decreasing HRT, pushing the community toward higher activity.
• Microbial Community Analysis: At the end of the experiment or key phases, biofilm samples are scraped from the carriers. DNA is extracted and analyzed via 16S rRNA gene amplicon sequencing (e.g., targeting V3-V4 region with primers 338F/806R) and/or qPCR to quantify the abundance and diversity of anammox keystone species and associated bacteria [64].

Visualization of Microbial Succession and Interactions

The following diagram illustrates the ecological succession and metabolic interactions between keystone species within a carrier-based partial nitritation-anammox (PN/A) biofilm, as revealed by longitudinal metagenomic studies [66] [25].

Diagram 1: Keystone species succession and stratified interactions in a PN/A biofilm on a carrier surface. The process begins with the attachment of early colonizers, whose metabolic activity facilitates the subsequent arrival of anammox bacteria (AnAOB) and the formation of a metabolically stratified structure. Operational changes, such as increased nitrogen loading, can trigger a shift in the dominant AOB population, highlighting the dynamic nature of these engineered ecosystems [66] [25].

The Scientist's Toolkit: Essential Research Reagents and Materials

Table 3: Key Research Reagents and Materials for Investigating Carriers and Biofilms

Reagent/Material	Function/Application	Specific Examples / Notes
Polyurethane (PU) Sponge	High-porosity carrier for 3D biofilm growth and biomass retention.	Apparent density ~25.81 mg/cm³; bulk density ~14.52 mg/cm³ [64].
Polyvinyl Alcohol (PVA)	Hydrophilic, modifiable biocarrier polymer.	Can be esterified with ε-caprolactone to tailor Surface Free Energy [61].
Graphene Oxide (GO)	Nanomaterial enhancer for microbial metabolism and activity.	Often used as a composite (e.g., GO@PU) at loadings of ~10 mg/L [65].
Mineral Carriers	Components for in-situ permeable barriers in groundwater bioremediation.	Includes quartz, kaolin clay, bentonite clay, and natural zeolites [62] [63].
Trace Element Solutions	Essential micronutrients for autotrophic bacteria (AOB, AnAOB).	Typically contains Fe, Zn, Co, Cu, Ni, Mo, Mn, Se, B, and chelators like EDTA [64].
DNA Extraction Kit	Isolation of metagenomic DNA from carrier biofilms for community analysis.	e.g., FastDNA SPIN kit for soil [66].
16S rRNA Primers	Amplicon sequencing of prokaryotic communities from biofilm DNA.	Universal primers, e.g., 338F / 806R for the V3-V4 hypervariable region [64].
Atomic Force Microscope	Quantitative measurement of single-cell adhesion forces to carrier materials.	Requires functionalized cantilevers with microbial cells or carrier substrates [61].

Carrier materials are pivotal engineering tools that directly manipulate the ecology of anammox systems by shaping the habitat for keystone species. The strategic design of carriers—optimizing surface energy, porosity, and chemical composition—governs the initial adhesion, spatial organization, and metabolic interplay between AOB and AnAOB. This direct influence on biofilm architecture and community succession enables faster process start-up, robust operational stability, and enhanced nitrogen removal efficiency. Future research will likely focus on "smart" carriers with tailored surface chemistries, controlled-release properties for micronutrients, and perhaps even biosensing capabilities. A deeper, systems-level understanding of the metabolic networks between anammox bacteria and their synergistic heterotrophic partners, as facilitated by these advanced materials, will be crucial for pushing the boundaries of mainstream anammox application and achieving new levels of sustainability in wastewater treatment.

Strategies for Preventing Keystone Species Loss and Community Collapse

The anaerobic ammonium oxidation (anammox) process represents a critically important microbial pathway for nitrogen removal in both natural and engineered ecosystems. Within the complex microbial consortia facilitating this process, certain taxa function as keystone species whose presence and ecological activity disproportionately impact community structure, functional stability, and ecosystem performance. The loss of these keystone species can trigger community collapse, leading to impaired nitrogen cycling and system failure. In estuarine and coastal sediments, anammox bacteria play an essential role in nitrogen loss, with Candidatus Scalindua identified as a predominant and often keystone genus [11]. Recent research has demonstrated that rare species in these communities may also play crucial roles in maintaining ecological stability [11]. This technical guide synthesizes cutting-edge research to present evidence-based strategies for identifying, preserving, and protecting keystone species within anammox systems, with particular emphasis on conservation approaches that enhance community resilience to environmental perturbation. The principles outlined herein are essential for maintaining the stability and functionality of both natural and engineered ecosystems dependent on anammox processes for nitrogen removal.

Identification and Quantification of Keystone Species

Accurate identification of keystone species is the foundational step in developing effective conservation strategies. Traditional approaches that relied solely on binary network descriptions have proven insufficient, as they ignore both the strength of trophic links and critical indirect effects [67]. Modern identification methodologies employ integrated approaches combining high-throughput sequencing, co-occurrence network analysis, and quantitative flux modeling to pinpoint species whose impact on community stability exceeds their relative abundance.

Advanced Network Analysis Techniques

The application of weighted network analysis represents a significant advancement over simple binary connectivity metrics. In quantitative food webs, keystone species identification should incorporate:

• Weighted Betweenness Centrality (wBC): Measures the fraction of shortest paths passing through a node, with link weights based on quantitative fluxes
• Weighted Closeness Centrality (wCC): Assesses how quickly a node can interact with others in the network via weighted connections
• Topological Importance: Evaluates a node's impact on network connectivity and function at multiple scales [68]

Research comparing these indices has demonstrated that link weight and indirect effects have stronger influences on food-web stability than the simple removal of highly connected 'hubs' [67]. This paradigm shift emphasizes that both quantitative fluxes and species dissipating their effects across many links should be of primary concern in conservation planning. In anammox systems, co-occurrence network analysis has identified Candidatus Scalindua as a keystone genus, with rare species playing unexpectedly important roles in maintaining network structure [11].

Experimental Protocols for Keystone Species Identification

Protocol 1: Quantitative Network Construction for Anammox Communities

• Sample Collection: Collect sediment cores or bioreactor samples at appropriate depth intervals (e.g., 1-4 cm depending on system)
• DNA Extraction: Extract total DNA from 0.5 g wet sediment/sludge using FastDNA SPIN Kit for soil
• Gene Amplification: Amplify anammox bacterial 16S rRNA gene with primers Brod541F and Amx820R
• High-Throughput Sequencing: Perform sequencing on Illumina platform; denoise sequences using Sickle program
• OTU Assignment: Cluster sequences into operational taxonomic units (OTUs) at 98% similarity using QIIME 2
• Network Construction: Build co-occurrence networks using SparCC or similar method; calculate network metrics
• Keystone Identification: Identify keystone taxa based on high wBC and wCC values, and phylogenetic significance [11]

Protocol 2: Stable Isotope-Based Interaction Strength Quantification

• Sample Preparation: Collect primary producers and consumers from ecosystem
• Stable Isotope Analysis: Measure δ15N and δ13C values for all specimens
• Mixing Models: Use SIAR or MixSIAR to quantify resource contributions to consumers
• Link Weight Assignment: Calculate interaction strengths based on proportional dietary contributions
• Weighted Network Metrics: Compute wBC and wCC using customized scripts [68]

Table 1: Key Metrics for Keystone Species Identification in Anammox Systems

Metric	Calculation Method	Interpretation	Application in Anammox Systems
Degree Centrality	Number of direct connections	Identifies highly connected species	Limited utility; may overestimate hub importance
Weighted Betweenness Centrality (wBC)	Proportion of shortest paths passing through a node, weighted by interaction strength	Identifies species critical for information/material flow	Candidatus Scalindua shows high wBC in coastal sediments [11]
Weighted Closeness Centrality (wCC)	Inverse of sum of shortest paths to all other nodes, weighted by interaction strength	Identifies species that can quickly interact with others	Rare taxa may have high wCC despite low abundance [11]
Nearest Taxon Index (NTI)	Measures phylogenetic clustering	Indicates environmental filtering vs. competitive exclusion	Anammox communities show distinct phylogenetic differentiation [11]

Critical Threats to Keystone Species in Anammox Systems

Understanding the specific threats facing keystone species is essential for developing targeted conservation strategies. Research across diverse ecosystems has identified several primary drivers of keystone species loss with particular relevance to anammox bacterial communities.

Environmental Pressures and Habitat Degradation

Anammox bacteria, particularly keystone taxa like Candidatus Scalindua, demonstrate heightened sensitivity to specific environmental parameters:

• Dissolved Oxygen (DO) Fluctuations: Anammox bacteria are strictly anaerobic and experience inhibition at DO concentrations >0.1 mg·L⁻¹ [4]
• Temperature Variability: Optimal temperature ranges for anammox bacteria are narrow (30-37°C); deviations reduce activity and growth rates [69]
• Organic Matter Loading: Elevated organic carbon can shift competitive advantages toward heterotrophic denitrifiers, suppressing anammox processes [11]
• Chemical Inhibitors: Specific wastewater constituents (e.g., pharmaceuticals, heavy metals) disproportionately affect anammox bacteria [69]

Spatial heterogeneity in anammox bacterial communities across estuaries like the Changjiang Estuary (CJE), Oujiang Estuary (OJE), Jiulong River Estuary (JLE), and the South China Sea (SCS) demonstrates the niche specialization of different taxa [11]. This specialization makes keystone species particularly vulnerable to habitat changes that exceed their narrow tolerance ranges.

Community Assembly Disruption

The mechanisms governing anammox community assembly play a crucial role in maintaining keystone species. Analysis of coastal sediments indicates that ecological drift predominantly shapes the overall anammox bacterial community, while rare species (including some keystone taxa) are more susceptible to dispersal limitations and environmental selection [11]. Environmental changes that disrupt these natural assembly processes can lead to disproportionate losses of keystone species through:

• Barriers to Dispersal: Physical barriers or habitat fragmentation that prevent natural colonization processes
• Homogenization of Environmental Conditions: Reduction in microenvironment diversity that supports specialized taxa
• Stochastic Extinction Events: Random population fluctuations that disproportionately affect low-abundance keystone species

Preservation Strategies for Anammox Keystone Species

Biomass Preservation Protocols

Effective preservation of anammox biomass is crucial for maintaining keystone species, enabling reactivation of failed systems, and ensuring conservation of microbial diversity. Multiple preservation methods have been experimentally validated:

Protocol 3: Cryopreservation of Anammox Biomass

• Pre-treatment: Clean biomass and add cryoprotectant agents (CPAs) including trehalose, sodium alginate, skim milk, dimethyl sulfoxide (DMSO), or tryptone soy broth (TSB)
• Packaging: Transfer treated biomass to cryovials under anaerobic conditions
• Freezing: Implement controlled-rate freezing (1°C/min) to -80°C
• Storage: Maintain at -80°C or in liquid nitrogen for long-term preservation
• Reactivation: Thaw rapidly at 37°C and inoculate into pre-conditioned reactor [69]

Protocol 4: Gel Encasement Preservation

• Biomass Preparation: Concentrate anammox biomass via centrifugation
• Gel Formation: Immobilize biomass in polyvinyl alcohol (PVA) or sodium alginate gels
• Dehydration: Partially dehydrate gels to stabilize metabolic activity
• Storage: Store at 4°C under anaerobic conditions
• Reactivation: Rehydrate and directly transfer to reactor systems [69]

Table 2: Comparison of Anammox Biomass Preservation Techniques

Method	Preservation Duration	Recovery Efficiency	Key Advantages	Limitations
Refrigeration (4°C)	Short-term (weeks)	60-80% activity retention	Simple protocol, low cost	Rapid viability decline
Cryopreservation (-80°C)	Long-term (years)	70-90% activity retention	Maximum preservation duration	Requires specialized equipment
Gel Encasement	Medium-term (months)	65-85% activity retention	Protects from environmental stress	Complex preparation
Lyophilization	Long-term (years)	50-70% activity retention	Easy storage and transport	Significant initial activity loss

Community-Based Conservation Approaches

Given the importance of microbial interactions in anammox systems, conservation strategies must extend beyond individual species to encompass the entire community:

• Maintenance of Cross-Feeding Relationships: Preserve metabolic interdependencies where anammox bacteria (e.g., Ca. Kuenenia) produce acetate and glycogen that enhance microbial interactions and support biofilm formation [4]
• Functional Group Retention: Protect non-anammox community members (e.g., Chloroflexi, Proteobacteria) that constitute the "anammox core" and support keystone species function [4]
• Biofilm Integrity Conservation: Prioritize preservation of intact biofilm structures that provide optimal microenvironments for keystone species

Recent research demonstrates that microbial interactions play a keystone role in rapid anammox sludge proliferation and biofilm formation, with specific taxa like Ca. Kuenenia enhancing metabolic cooperation through production of microbial substrates [4].

Reactivation and Restoration Protocols

Following preservation or system disturbance, effective reactivation protocols are essential for restoring keystone species and community functionality.

System Reactivation Strategies

Protocol 5: Reactivation of Preserved Anammox Biomass

• Reactor Preparation: Establish anaerobic conditions with DO <0.1 mg·L⁻¹ in UASB or similar reactor configuration
• Nutrient Medium: Provide synthetic wastewater with NH₄⁺ (70-100 mg·L⁻¹) and NO₂⁻ (70-100 mg·L⁻¹)
• Inoculation: Introduce preserved biomass at volatile suspended solids (VSS) of 4150 ± 200 mg·L⁻¹
• Gradual Loading: Initiate at low nitrogen loading rate (0.5 kg-N·m⁻³·d⁻¹) and increase gradually
• Performance Monitoring: Track ammonium removal efficiency (ARE), nitrogen removal efficiency (NRE), and specific anammox activity (SAA) [69]

Protocol 6: Biofilm Enhancement for Keystone Species Recovery

• Carrier Selection: Provide appropriate attachment surfaces (sponges, polyurethane carriers)
• Inoculation Strategy: Combine floccular and biofilm-based inoculation
• Extracellular Polymeric Substances (EPS) Enhancement: Optimize conditions for EPS production (C/N ratio, controlled starvation periods)
• Quorum Sensing Activation: Utilize N-acyl-homoserine lactone-mediated signaling to enhance biofilm formation [69]
• Community Analysis: Monitor keystone species abundance via qPCR and network analysis

Performance Monitoring and Stability Assessment

Successful reactivation requires comprehensive monitoring to ensure keystone species recovery and community stability:

• Molecular Monitoring: Track keystone species abundance using qPCR targeting specific 16S rRNA genes
• Network Analysis: Regularly assess co-occurrence networks to verify restoration of keystone species centrality
• Functional Metrics: Monitor nitrogen transformation rates, heme c content, and EPS composition
• Stability Indicators: Evaluate resistance to perturbation and recovery trajectory

Research indicates that reactivation performance varies significantly with preservation method, with cryopreservation typically yielding the highest recovery rates of anammox activity (70-90% retention) [69].

The Scientist's Toolkit: Research Reagent Solutions

Table 3: Essential Research Reagents for Keystone Species Studies

Reagent/Kit	Application	Function	Example Use in Anammox Research
FastDNA SPIN Kit for Soil (MP Biomedical)	DNA extraction	Efficient lysis and purification of microbial DNA from sediment/sludge	DNA extraction for 16S rRNA sequencing of anammox communities [11]
Brod541F and Amx820R Primers	PCR amplification	Specific amplification of anammox bacterial 16S rRNA gene	Target amplification for high-throughput sequencing [11]
QIIME 2 Platform	Bioinformatic analysis	Processing and analysis of high-throughput sequencing data	OTU assignment and diversity calculations [11]
Cryoprotectant Agents (Trehalose, DMSO)	Biomass preservation	Protect cellular integrity during freezing	Cryopreservation of anammox biomass for long-term storage [69]
Stable Isotopes (¹⁵N, ¹³C)	Metabolic tracing	Quantification of material flows through food webs	Construction of quantitative food webs for keystone identification [68]
Extracellular Polymeric Substances (EPS) Extraction Kits	Biofilm characterization	Isolation and quantification of EPS components	Monitoring biofilm development during reactor reactivation [69]

Visualizing Keystone Species Protection Strategies

The following diagram illustrates the integrated approach to keystone species conservation in anammox systems, highlighting the interconnected nature of identification, preservation, and reactivation strategies:

Keystone Species Conservation Workflow

This integrated approach emphasizes the cyclical nature of effective keystone species conservation, where continuous monitoring informs refinement of identification criteria and preservation protocols.

Protecting keystone species in anammox bacterial communities requires a multifaceted approach that integrates advanced identification methods, effective preservation techniques, and strategic reactivation protocols. The strategies outlined in this technical guide emphasize the critical importance of:

• Moving beyond binary network analysis to incorporate quantitative flux measurements and indirect effects in keystone species identification
• Implementing appropriate preservation protocols that maintain both keystone species viability and their essential ecological interactions
• Adopting community-based conservation approaches that recognize the embeddedness of keystone species within complex microbial networks

As research continues to reveal the intricate relationships between anammox keystone species and overall ecosystem function, these strategies will become increasingly vital for maintaining the stability and functionality of both natural and engineered nitrogen cycle systems. Future directions should focus on refining quantitative network approaches, developing more effective preservation methodologies for complex microbial communities, and establishing standardized monitoring frameworks for detecting early warning signs of keystone species decline. Through implementation of these evidence-based strategies, researchers and engineers can significantly enhance the resilience of anammox systems to environmental perturbation and prevent the community collapse that typically follows keystone species loss.

Optimizing Environmental Conditions to Support Keystone Metabolic Functions

Within the complex microbial consortia of anaerobic ammonium oxidation (anammox) systems, keystone species perform specialized metabolic functions that are critical for community stability and ecosystem functioning. These taxa, though often not the most abundant, exert a disproportionate influence on community structure and biogeochemical cycling [40]. In anammox reactors and natural environments, the metabolic activities of these keystone organisms—including nitrogen transformation, carbon metabolism, and the production of essential metabolites—are highly sensitive to environmental conditions [70] [4]. This technical guide synthesizes current research on optimizing key operational parameters to support these keystone functions, framed within the broader context of keystone species research in anammox bacterial communities.

The ecological importance of keystone taxa extends beyond their functional roles to their position within microbial interaction networks. Studies of anammox granular sludge have revealed that rare microbial sub-communities often exhibit contrasting assemblage patterns and metabolic functions compared to abundant taxa, yet play indispensable roles in maintaining ecological stability [71] [11]. Understanding how environmental factors shape these keystone functions is essential for advancing anammox applications in wastewater treatment and environmental biotechnology.

Key Environmental Factors Affecting Keystone Metabolic Functions

Nitrogen Loading and Substrate Fluctuations

The nitrogen-loading rate (NLR) represents a critical control parameter influencing anammox system performance and keystone taxon stability. Research demonstrates a nonlinear relationship between NLR and system performance, with both excessive loading and nitrogen starvation leading to functional deterioration [6].

Table 1: Response of Anammox Systems to Nitrogen Loading Variations

Nitrogen Loading Condition	Anammox Bacterial Abundance	Denitrification Efficiency	Microbial Community Response
Optimal NLR (~1.38-3.68 kg/m³·d)	Increases (from 5.85% to 11.43%)	Enhanced	Stable community structure with balanced interactions
Excessive NLR (>3.68 kg/m³·d)	Reduced	Deterioration	Increased modularity (0.563 index); reinforced interspecies interactions
Nitrogen Starvation	Reduced	Deterioration	Increased modularity (0.545 index); stress response mechanisms

The mechanisms underlying these responses involve both direct physiological inhibition and community-level adaptations. Under high substrate conditions, residual ammonia nitrogen and nitrite can be converted to free ammonia and free nitrous acid, which diffuse into cells, alter intracellular pH, and cause cellular damage [6]. Microbial communities counteract loading stress through modular collaboration, with increased inter-module connectivity indicating reinforced interspecies interactions [6].

Inorganic Carbon Availability

Inorganic carbon serves as an essential assimilation carbon source for chemoautotrophic anammox bacteria and significantly influences community assembly patterns and metabolic functions. Different IC/TN ratios select for distinct functional groups within anammox consortia [71]:

Table 2: Optimal IC/TN Ratios for Key Functional Groups in Anammox Systems

Functional Group	Optimal IC/TN Ratio	Key Metabolic Functions	Representative Genera
Anammox Bacteria	0.62-1.24	Nitrogen removal via anammox pathway	Candidatus Brocadia, Candidatus Kuenenia
Ammonia-Oxidizing Bacteria	<0.31	Partial nitritation, ammonia oxidation	Nitrosomonas
Nitrite-Oxidizing Bacteria	>1.24	Nitrite oxidation, nitrate production	Nitrospira
Heterotrophic Denitrifiers	Variable	Carbon metabolism, nitrate reduction	Thauera, Burkholderiales

Lower IC/TN ratios (<0.31) benefit ammonia-oxidizing bacteria and Candidatus Brocadia, while moderate ratios (0.62-1.24) favor most anammox bacteria [71]. This specialization reflects niche partitioning and differential carbon utilization strategies among keystone taxa.

Dissolved Oxygen and Redox Conditions

Oxygen concentration represents a decisive factor shaping anammox community structure and function. Anammox bacteria are obligate anaerobes, but micro-aeration can create favorable conditions for complementary nitrogen-cycling microorganisms [72]. In full-scale wastewater treatment plants, dissolved oxygen below 0.5 mg/L combined with high ammonia-nitrogen loading stimulates anammox bacterial growth, though DO control alone is insufficient for stable process establishment [72].

The interplay between oxygen sensitivity and microbial cooperation creates ecological niches for keystone taxa with different oxygen affinities. The presence of ammonia-oxidizing bacteria under micro-aerobic conditions can generate the necessary nitrite for anammox metabolism, creating a syntrophic relationship that enhances overall nitrogen removal [4].

Organic Carbon and C/N Ratio

The influent C/N ratio significantly influences the balance between autotrophic and heterotrophic processes in anammox systems. When the C/N ratio exceeds 3, heterotrophic denitrifying bacteria proliferate rapidly and compete with anammox bacteria for nitrite, slowing anammox growth [72]. However, moderate organic carbon availability supports cross-feeding interactions where anammox bacteria provide metabolites for heterotrophic partners who in turn remove potential inhibitors.

Recent research indicates that keystone anammox bacteria like Candidatus Kuenenia might produce acetate and glycogen to enhance microbial interactions and facilitate biofilm formation [4]. This metabolic versatility underscores the importance of understanding anammox keystone functions beyond nitrogen transformation alone.

Experimental Protocols for Investigating Keystone Functions

Reactor Configuration and Operation for Nitrogen Loading Studies

Objective: To investigate the response of anammox keystone taxa and their metabolic functions to fluctuating nitrogen loads.

Methodology:

• Reactor Setup: Employ an expanded granular sludge bed reactor with 10 L working volume constructed of Plexiglas and insulated for temperature control [6].
• Inoculation: Use mature anammox sludge as seed material, with initial nitrogen loading rate of 1.38 ± 0.01 kg/m³·d and nitrogen removal rate of 1.15 ± 0.02 kg/m³·d after acclimation [6].
• Operational Phases:
  • Phase 1 (Days 1-15): Maintain stable NLR at 1.38 ± 0.01 kg/m³·d as control
  • Phase 2 (Days 16-65): Increase NH₄⁺-N concentration in 30 mg/L increments
  • Phase 3 (Days 66-105): Rapidly decrease then gradually increase substrate concentration
  • Phase 4 (Days 106-120): Implement nitrogen starvation conditions [6]
• Analytical Monitoring: Daily measurement of nitrogen species in effluent; periodic sampling for microbial community analysis.

This protocol enables systematic investigation of NLR effects on anammox performance, keystone taxon abundance, and metabolic pathway expression.

Assessing Microbial Community Assembly Under Inorganic Carbon Stress

Objective: To elucidate the assemblage patterns and functional traits of abundant versus rare microbial sub-communities under different IC/TN ratios.

Methodology:

• Reactor Configuration: Operate an expanded granular sludge blanket reactor with 2.0 L working volume [71].
• IC/TN Manipulation:
  • Phase I: Baseline operation with IC/TN ~1.24
  • Phase II: Reduce IC/TN to 0.31
  • Phase III: Restore IC/TN to 1.24
  • Phase IV: Increase IC/TN to 1.80 [71]
• Microbial Community Analysis:
  • DNA extraction and 16S rRNA amplicon sequencing
  • Classification of taxa into abundant (>1% in all samples) and rare (<0.1% in all samples) sub-communities [71]
  • Co-occurrence network analysis using SparCC method with 100 iterations
  • Functional prediction with PICRUSt2 and KEGG pathway annotation [71]
• Statistical Analysis: Calculate network topological features (modularity, connectivity, centrality) to identify keystone taxa.

This approach reveals how IC stress reshapes microbial interactions and selects for specialized keystone functions.

Metabolic Pathways and Microbial Interactions

Key Metabolic Functions of Anammox Keystone Taxa

Anammox keystone taxa perform specialized metabolic functions that sustain the broader microbial community and ecosystem processes:

• Nitrogen Metabolism: The core anammox metabolism converts ammonium and nitrite to dinitrogen gas, with hydrazine synthase (hzs) and nitrite reductase (nir) as key enzymes [4]. This pathway is particularly important in estuarine and coastal sediments where anammox accounts for substantial Nâ‚‚ production [11].

• Carbon Fixation and Metabolism: Anammox bacteria utilize the Wood-Ljungdahl pathway for carbon fixation, with formate dehydrogenase (fdh) and carbon monoxide dehydrogenase (cdh) as critical enzymes [4]. These pathways provide organic carbon for heterotrophic community members.

• Phosphonate and Phosphinate Metabolism: Specialized phosphorus metabolism pathways have been identified as keystone functions in complex microbial communities, carried out by specific bacterial taxa including Nitrospira and Gemmatimonas [40].

The diagram below illustrates the key metabolic pathways and microbial interactions centered around anammox keystone taxa:

Diagram Title: Metabolic Networks in Anammox Keystone Consortia

Microbial Cross-Feeding and Cooperative Interactions

Keystone functions in anammox systems are sustained through complex interaction networks involving multiple microbial groups:

• Nitrite Looping: Heterotrophic denitrifying bacteria reduce nitrate produced by the anammox reaction back to nitrite, creating a metabolic loop that improves nitrogen removal efficiency [4]. This cross-feeding demonstrates how keystone functions emerge from multi-species cooperation.

• Metabolite Exchange: Anammox bacteria produce organic compounds including acetate and glycogen that support heterotrophic partners, who in turn help maintain redox balance and remove inhibitory substances [4].

• Biofilm Formation: Specific microbial groups initiate biofilm formation by attaching to carriers, followed by growth of anammox bacteria and symbiotic species that rely on extracellular polymeric substances secreted by pioneer colonizers [4].

The stability of these interaction networks depends critically on environmental conditions, particularly nitrogen loading, inorganic carbon availability, and oxygenation.

Research Reagent Solutions and Essential Materials

Table 3: Essential Research Reagents and Materials for Anammox Keystone Function Studies

Reagent/Material	Specification/Composition	Function/Application	Key Considerations
Synthetic Wastewater Base	NaHCO₃ (0.5 g/L), KHCO₃ (0.5 g/L), KH₂PO₄ (0.027 g/L), MgSO₄·7H₂O (0.02 g/L), CaCl₂·2H₂O (0.136 g/L	Provides essential ions and buffer capacity	Maintain alkalinity and pH stability; phosphate for microbial growth
Nitrogen Sources	(NH₄)₂SO₄ and NaNO₂ as requested	Ammonium and nitrite substrates for anammox metabolism	Maintain NH₄⁺:NO₂⁻ ratio at ~1:1.32 [6]
Trace Element Solution I	EDTA (5 g/L), FeSO₄·7H₂O (5 g/L)	Iron chelation and delivery	Prevents iron precipitation and improves bioavailability
Trace Element Solution II	EDTA (5 g/L), NaMoO₄·2H₂O (0.22 g/L), NiCl₂·6H₂O (0.19 g/L), CuSO₄·5H₂O (0.25 g/L), CoCl₂·6H₂O (0.24 g/L), ZnSO₄·7H₂O (0.43 g/L), MnCl₂·4H₂O (0.99 g/L)	Essential micronutrients for metalloenzymes	Critical for anammox metabolism including hydrazine synthase
DNA Extraction Kit	PowerSoil DNA Isolation Kit or equivalent	Microbial community DNA extraction	Standardized for reproducibility in molecular analyses
PCR Reagents	Premix Ex Taq, specific primers (Brod541F, Amx820R) [11]	Amplification of anammox bacterial 16S rRNA genes	Enables community composition and diversity analysis

Optimizing environmental conditions to support keystone metabolic functions requires integrated management of multiple parameters, including nitrogen loading rate, inorganic carbon availability, dissolved oxygen, and organic carbon inputs. The specialized metabolic functions embedded in keystone taxa—particularly nitrogen metabolism, carbon fixation, and phosphorus transformation—are highly sensitive to these environmental conditions [40]. Maintaining operational parameters within optimal ranges supports not only the abundance of keystone taxa but, more importantly, their functional expression and ecological interactions.

Future research should focus on elucidating the specific metabolic pathways and regulatory mechanisms that enable keystone taxa to maintain community stability under fluctuating environmental conditions. Advanced approaches combining metagenomics, metatranscriptomics, and metabolic network analysis will further illuminate the complex interplay between environmental factors, keystone functions, and ecosystem processes in anammox systems.

Cross-System Verification: Validating Keystone Roles from Bioreactors to Global Ecosystems

Within the broader thesis on keystone species in anammox bacterial communities, this guide addresses the critical challenge of tracking the stability of these pivotal organisms over time. Keystone species are not necessarily the most abundant taxa but are those whose impact on community structure and function is disproportionately large relative to their abundance. In engineered anammox systems, these species exert fundamental influence on nitrogen removal performance, community assembly, and functional resilience [29]. Longitudinal studies that monitor these species throughout reactor enrichment and operation are therefore essential for diagnosing stability thresholds, predicting process failure, and designing robust bioaugmentation strategies. This technical guide synthesizes current methodologies and insights for effectively tracking keystone species dynamics, with a specific focus on anammox bioreactors as model systems.

Keystone Species in Anammox Communities: Ecological Foundations

In anammox bioreactors, keystone species often include specific anammox bacteria such as Candidatus Brocadia, Candidatus Jettenia, and Candidatus Kuenenia, but also extend to synergistic heterotrophic partners from phyla like Chloroflexi and Proteobacteria [29] [31]. These associated heterotrophs can form stable, cooperative relationships with anammox bacteria, providing essential vitamins, cofactors (e.g., molybdopterin), and amino acids, while consuming potentially inhibitory organic carbon [31] [29]. The stability of these keystone partnerships is not static but varies with operational conditions. For instance, the dominance of specific anammox species can shift reversibly in response to factors like nitrogen load, with Candidatus Brocadia caroliniensis giving way to Candidatus Brocadia sinica at high nitrite concentrations (>340 mg N/L) [73]. Furthermore, microbial networks respond to internal and external pressures by altering interaction patterns; increased negative interactions between anammox bacteria and heterotrophs have been observed as a mechanism to maintain functional stability over time [29].

Methodological Framework for Longitudinal Tracking

Tracking keystone species requires a multi-faceted approach that links community structure with function over time. The following workflow outlines the integrated methodology.

Experimental Workflow for Longitudinal Tracking

The following diagram illustrates the comprehensive workflow for designing and executing a longitudinal study of keystone species in bioreactors.

Core Analytical Techniques

3.2.1 High-Throughput Sequencing

• 16S rRNA Gene Amplicon Sequencing: Tracks general community composition and diversity shifts over time. The V3-V4 hypervariable region is commonly targeted using 515F/806R primers [74] [75]. This method provides a cost-effective way to monitor broad community changes but offers limited functional insights.
• Metagenomic Sequencing: Provides deep functional insights by sequencing all genetic material in a sample. This allows for the reconstruction of Metagenome-Assembled Genomes (MAGs) and analysis of functional genes. Quality trimming (e.g., using Sickle at Phred quality <20) and co-assembly with tools like MEGAHIT are standard pre-processing steps [74].
• Metatranscriptomic Sequencing: Reveals actively expressed genes, providing a dynamic view of microbial activity and responses to environmental changes that cannot be deduced from DNA-based metagenomics alone.

3.2.2 Network Analysis Co-occurrence networks are constructed from abundance data to infer potential microbial interactions. Keystone species are identified through topological features using:

• Betweenness Centrality: Identifies nodes that act as bridges between different parts of the network.
• Zi-Pi Plot Analysis: Classifies nodes into categories (e.g., network hubs, connectors) based on their within-module connectivity (Zi) and among-module connectivity (Pi) [6]. Modules are subgroups of highly interconnected taxa, and their dynamics (e.g., modularity indices) reflect community stability in response to perturbations [6] [29].

3.2.3 Functional Profiling

• PICRUSt2: Predicts functional potential from 16S rRNA data by mapping to databases like KEGG [6].
• MAG Analysis: Provides high-resolution genomic information for specific taxa, enabling the reconstruction of metabolic pathways and potential interactions. Quality thresholds (e.g., >50% completion and <10% contamination) are critical for reliable analysis [74].

Table 1: Key Bioinformatic Tools for Keystone Species Analysis

Tool Name	Primary Function	Application in Keystone Species Tracking
QIIME 2/Mothur	16S rRNA amplicon data processing	Processes raw sequencing data into Amplicon Sequence Variants (ASVs) for diversity and composition analysis
MetaWRAP	Metagenomic binning and analysis	Integrates multiple binning algorithms (MetaBAT2, MaxBin2, CONCOCT) to generate high-quality MAGs
CheckM	MAG quality assessment	Evaluates completion and contamination of MAGs using lineage-specific marker genes
Cytoscape/Gephi	Network visualization and analysis	Visualizes and calculates topological properties of microbial co-occurrence networks
PICRUSt2	Functional prediction	Infers functional potential of communities from 16S rRNA gene data

Quantitative Assessment of Keystone Species Dynamics

Longitudinal tracking in anammox systems has revealed several key patterns in keystone species behavior under different operational regimes.

Table 2: Keystone Species Responses to Operational Parameters in Anammox Bioreactors

Operational Parameter	Impact on Keystone Species	Functional Outcome	Reference
Nitrogen Loading Rate (NLR) Fluctuations	Nonlinear response: abundance increases with moderate NLR but decreases sharply beyond threshold (>3.68 kg/m³·d)	Synchronous enhancement then deterioration of nitrogen removal efficiency; increased network modularity (0.545-0.563)	[6]
Carrier Type (Biofilm vs. Granular)	Ca. B. caroliniensis and Ca. Jettenia prefer biofilm attachment; Ca. B. sinica forms granules	Impacts biomass retention strategy and reactor configuration selection	[35] [73]
Historical Environmental Instability	Proliferation of specialist taxa; reduced resistance to subsequent disturbances	Increased susceptibility to deterministic assembly; lower functional resistance to organic loading shocks	[74]
Organic Carbon Supplementation	Enables coexistence of Ca. B. caroliniensis with heterotrophic communities; inhibits some anammox species with NO present	Emergence of heterotrophic networks that coevolve with anammox bacteria; differential species selection	[73]
Long-term Stable Operation	Dynamic subnetworks with alternating importance of Chloroflexi and Proteobacteria; shift in keystone species	Maintenance of functional stability despite changing species interactions	[29]

Essential Research Reagent Solutions

The following table details key reagents and materials required for longitudinal tracking studies in anammox bioreactors.

Table 3: Essential Research Reagents for Keystone Species Tracking

Reagent/Material	Specification/Function	Application Example
DNA Extraction Kit	Phenol:chloroform-based or commercial kits (e.g., CTAB method); must effectively lyse tough bacterial cell walls	Total community DNA extraction for subsequent sequencing [74]
PCR Primers	515F/806R for 16S rRNA; hzsB for anammox-specific detection; nirS/nirK for denitrifier community	Amplification of target genes for community composition and functional group analysis [31] [74]
Synthetic Wastewater Components	(NH₄)₂SO₄ (NH₄⁺-N source); NaNO₂ (NO₂⁻-N source); NaHCO₃/KHCO₃ (buffer & inorganic carbon)	Maintains defined substrate conditions for anammox enrichment and function [31] [6]
Trace Element Solutions	Contains EDTA, FeSO₄·7H₂O, and micronutrients (Mo, Ni, Cu, Co, Zn, Mn)	Essential for anammox bacterial growth and metabolism [31] [6]
Biofilm Carriers	Polyurethane sponge fillers, polyethylene hollow fiber membranes, plastic hollow cylinders	Provides attachment surface for biofilm-forming anammox species; enhances biomass retention [35] [31]
Sequencing Standards	Mock communities with known composition; internal standard genes	Validates sequencing accuracy and quantifies absolute abundances in community data

Signaling and Metabolic Interactions in Anammox Consortia

Understanding the metabolic dependencies between anammox bacteria and their keystone partners is crucial for interpreting stability data. The following diagram illustrates these complex interactions.

As illustrated, anammox bacteria (keystone species) convert ammonium (NH₄⁺) and nitrite (NO₂⁻) to nitrogen gas (N₂) while secreting extracellular polymeric substances (EPS) [29]. Associated heterotrophic bacteria (keystone partners) consume this EPS and organic carbon, while in return providing essential nutrients including amino acids, cofactors, and vitamins that support anammox metabolism [31] [29]. This cross-feeding creates interdependent relationships that stabilize the community. Longitudinal disruptions to these metabolic exchanges often precede keystone species collapse and system failure.

Longitudinal tracking of keystone species in enrichment bioreactors requires integrating advanced molecular tools with ecological theory to decipher the complex dynamics that underpin functional stability. The methodologies outlined in this guide—from multi-omics integration to network analysis and metabolic reconstruction—provide a robust framework for identifying, monitoring, and interpreting keystone species behavior over time. Within the broader thesis context, this approach reveals that keystone stability in anammox systems is not merely a function of species abundance but emerges from the dynamic interplay between environmental conditions, network architecture, and metabolic interdependencies. Future research should focus on establishing quantitative thresholds for keystone species stability and developing early-warning indicators that can predict community collapse before functional failure occurs.

Metagenome-Assembled Genomes (MAGs) for Unraveling Keystone Metabolism

Metagenome-Assembled Genomes (MAGs) have emerged as a transformative tool in microbial ecology, enabling the decoding of metabolic blueprints and ecological roles of uncultured microorganisms. Within complex bacterial communities, keystone taxa exert disproportionate influence on community structure and function despite their low abundance. This technical guide explores the integrated application of MAGs and multi-omics approaches to unravel keystone metabolism in anaerobic ammonium oxidation (anammox) systems. We detail methodologies for MAG reconstruction from complex microbial aggregates, computational pipelines for predicting metabolic interactions, and experimental frameworks for validating keystone functions. The synthesis of genomic data from anammox consortia reveals that keystone taxa sustain system stability through specialized metabolic functions including nitrogen transformation, cofactor synthesis, and metabolite exchange. This whitepaper provides researchers with advanced protocols for extracting meaningful biological insights from MAGs, thereby accelerating the discovery of keystone organisms and their metabolic networks in engineered ecosystems.

The study of keystone species—organisms with disproportionate ecological impact—has been revolutionized by culture-independent genomic techniques. Metagenome-Assembled Genomes (MAGs) represent chromosome-level reconstructions of microbial genomes directly from environmental sequences, bypassing cultivation requirements. In anammox systems, where keystone taxa maintain ecosystem stability through specialized metabolic functions [38], MAGs provide unprecedented access to genetic determinants of their keystone roles.

Taxonomic novelty discovery through MAGs has revealed extensive uncharacterized diversity in anammox systems. For instance, in Lake Barkol sediment samples, approximately 97% of 309 reconstructed MAGs could not be classified at the species level, indicating substantial taxonomic novelty [76]. Similarly, analysis of subsea tunnel biofilms yielded three MAGs representing a novel family of anammox bacteria with average amino acid identity (AAI) of only 48-50% compared to known genera [77]. This expanded taxonomic resolution enables precise identification of keystone populations that drive community assembly despite their potentially low abundance.

The functional profiling capabilities of MAGs allow researchers to move beyond taxonomic inventories to mechanistic understanding of keystone functions. By annotating metabolic pathways, stress response systems, and biosynthetic capabilities, MAGs reveal how keystone taxa influence community stability. In anammox consortia, MAG-based analyses have identified keystone metabolic functions including "nitrogen metabolism" and "phosphonate and phosphinate metabolism" carried out by specific bacterial taxa such as Nitrospira and Gemmatimonas [38]. These specialized metabolic pathways, embedded within keystone taxa, create dependency networks that sustain microbiome stability under environmental perturbation.

Methodological Framework for MAGs Analysis

Experimental Design and Sample Preparation

Robust MAG reconstruction begins with strategic experimental design that accounts for microbial community heterogeneity across different habitat microenvironments. In anammox systems, distinct microbial assemblages occur in suspended flocs versus granular biofilms [4] [3], requiring separate processing for comprehensive keystone taxon recovery.

Table 1: Sample Collection and Processing Strategies for Anammox Systems

Sample Type	Collection Method	Biomass Processing	Preservation	Considerations
Suspended Flocs	Continuous filtration through 10-μm, 3-μm, and 0.22-μm polycarbonate membranes [76]	Sequential filtration with DNA extraction from 3-μm and 0.22-μm filters [76]	Immediate freezing at -80°C on dry ice [76]	Retains both free-living and particle-associated communities
Granular Sludge	Gravity settling followed by brief centrifugation [4]	Homogenization or direct lysis of intact granules	Flash freezing in liquid Nâ‚‚	Preserves spatial structure; may require disruption for DNA yield
Biofilms	Scraping from carrier surfaces [18]	Direct lysis or mechanical disruption	Storage in RNAlater or rapid freezing	Matrix-rich samples require enhanced lysis protocols

DNA extraction represents a critical step influencing MAG quality. For anammox sludge samples, the FastDNA SPIN Kit for Soil (MP Biomedicals) has been successfully employed [78], while water samples may require specialized kits like the ALFA-SEQ Advanced Water DNA Kit [76]. Extraction efficiency should be verified through fluorometric quantification, with recommended thresholds of OD₂₆₀/OD₂₈₀ ratios between 1.8-2.0 [76]. For multi-omics integration, parallel RNA extraction should utilize compatible methods that preserve transcript integrity.

Sequencing Strategies and MAG Reconstruction

Library preparation and sequencing parameters directly impact assembly completeness. For MAG reconstruction from anammox communities, Illumina HiSeq 4000 platforms generating 150-250 bp paired-end reads with 350 bp insert sizes have proven effective [78]. Higher-throughput platforms like NovaSeq or long-read technologies (PacBio, Nanopore) can enhance contiguity through hybrid assembly approaches.

The MAG reconstruction workflow involves multiple computational stages:

• Quality Control: Adapter trimming and read filtering using tools like fastp [79] or Trimmomatic
• Metagenome Assembly: Multiple k-mer approach with assemblers such as MEGAHIT or metaSPAdes
• Binning: Grouping contigs into putative genomes using composition and abundance metrics with tools like MetaBAT2, MaxBin2, or CONCOCT
• Bin Refinement: Dereplication and quality improvement using DAS Tool and MetaWRAP
• Quality Assessment: Evaluation of completeness and contamination with CheckM or similar tools

Table 2: Minimum Information Standards for MAG Quality Assessment

Quality Parameter	Threshold (High Quality)	Threshold (Medium Quality)	Assessment Tool
Completeness	>90%	>70%	CheckM [77]
Contamination	<5%	<10%	CheckM [77]
Strain Heterogeneity	<5%	<10%	CheckM
Presence of Marker Genes	16S rRNA, ≥1 tRNA	tRNA genes sufficient	Barrnap, tRNAscan-SE
N50	>10 kbp	>5 kbp	Assembly metrics

For taxonomic classification of anammox MAGs, the GTDB-Tk tool kit with the Genome Taxonomy Database provides standardized phylogenetic placement beyond traditional 16S rRNA classification [77]. This is particularly valuable for identifying novel anammox lineages, as demonstrated by the discovery of a new anammox family in subsea tunnel biofilms [77].

Metabolic Reconstruction and Interaction Networks

Functional annotation forms the foundation for identifying keystone metabolic capabilities. The recommended pipeline includes:

• Gene Prediction: Prodigal for protein-coding sequences
• Functional Assignment: Diamond BLASTP against KEGG, COG, and UniRef databases [79]
• Pathway Analysis: Manual curation of key metabolic pathways using KEGG Mapper or MetaCyc
• Comparative Genomics: Average Amino Acid Identity (AAI) calculations for phylogenetic relatedness [77]

For keystone function identification, several analytical approaches have proven effective:

• Co-occurrence networks: Statistical correlations across samples to identify tightly associated taxa
• Machine learning classification: Random forest models to identify taxa predictive of community stability [38]
• Metabolic complementarity analysis: Gap filling in metabolic pathways to identify cross-feeding dependencies
• Multi-omics integration: Correlation of metagenomic potential with metatranscriptomic and metaproteomic activity [78]

Figure 1: Computational Workflow for MAG Reconstruction and Keystone Taxon Identification

Keystone Metabolism in Anammox Systems

Nitrogen Cycling Pathways

Anammox consortia represent model systems for studying keystone metabolism due to their well-defined central function—anaerobic ammonium oxidation—and complex interdependencies. MAG-based analyses have revealed specialized nitrogen transformation capabilities distributed across community members.

The core anammox metabolism, carried out by Planctomycetes members like Candidatus Brocadia, Candidatus Kuenenia, and Candidatus Jettenia, involves the coupling of ammonium oxidation to nitrite reduction with hydrazine as an intermediate [18]. Key diagnostic genes include hydrazine synthase (hzs) and hydrazine dehydrogenase (hdh), which show distinct distribution patterns among anammox genera. For instance, hzs abundance in biofilm sludge was approximately 486 times higher than in granules in one comparative study [80].

Beyond the core pathway, MAGs reveal nitrogen cycling complementarity between anammox bacteria and coexisting populations:

• Nitrite looping: Heterotrophic denitrifiers like Thauera and Afipia reduce nitrate produced by anammox bacteria back to nitrite, creating a sustainable nitrite supply [31] [3]
• Ammonium regeneration: Dissimilatory nitrate reduction to ammonium (DNRA) by organisms like Ignavibacterium recycles nitrogen back to anammox substrates [3]
• Partial nitritation: Ammonia-oxidizing bacteria (AOB) like Nitrosomonas provide the necessary nitrite from ammonium under oxygen-limited conditions [78]

Table 3: Key Nitrogen Cycling Genes in Anammox Consortia

Gene	Function	Keystone Taxa	Application in MAG Analysis
hzsA	Hydrazine synthesis, alpha subunit	Ca. Brocadia, Ca. Jettenia [80]	Primary marker for anammox potential; quantifies functional abundance
hdh	Hydrazine dehydrogenase	Ca. Kuenenia, Ca. Scalindua	Confirms complete anammox pathway in MAGs
nirS	Cytochrome cd1 nitrite reductase	Thauera, Denitratisoma [31]	Distinguishes denitrifier-associated nitrite reduction
narG	Nitrate reductase	Ca. Brocadia, Ca. Kuenenia [18]	Identifies capacity for dissimilatory nitrate reduction
amoA	Ammonia monooxygenase	Nitrosomonas	Marks ammonia oxidizers in partial nitritation-anammox systems

Carbon and Energy Metabolism

Keystone metabolism in anammox systems extends beyond nitrogen transformation to encompass sophisticated carbon fixation and energy conservation strategies. MAGs from anammox reactors have revealed diverse carbon fixation pathways, including the Wood-Ljungdahl (reductive acetyl-CoA) pathway, which is folate-dependent in anammox bacteria [18].

Metabolic interactions around carbon sources create important dependencies:

• Volatile fatty acid (VFA) utilization: Certain anammox bacteria (e.g., Candidatus Brocadia fulgida) can oxidize formate, acetate, or propionate coupled to nitrate reduction [18]
• Glycogen metabolism: Candidatus Kuenenia may produce glycogen and acetate to enhance microbial interactions in biofilms [4]
• Cross-feeding: Heterotrophic bacteria provide essential metabolites including amino acids, cofactors, and vitamins to anammox bacteria [31]

The folate dependency of anammox bacteria represents a key metabolic interaction. Anammox bacteria lack complete folate biosynthesis pathways and rely on Proteobacteria and other community members to provide this essential cofactor for the Wood-Ljungdahl pathway [18]. This creates a obligate cross-feeding relationship that stabilizes the community structure.

Stress Response and Environmental Adaptation

Keystone taxa maintain system functionality under environmental fluctuations through specialized stress response mechanisms identified via MAG analysis:

• Osmoadaptation: Distinct "salt-in" (ion transport systems) and "salt-out" (compatible solute accumulation) strategies observed in hypersaline systems [76]
• Cold adaptation: EPS secretion and aggregation promoted at 25-15°C to protect anammox consortia [79]
• Oxygen tolerance: Detoxification systems including superoxide dismutase and catalase genes in anammox bacteria and their partners

Figure 2: Metabolic Interactions in Anammox Consortia

Research Reagent Solutions and Tools

Table 4: Essential Research Reagents and Computational Tools for MAG Analysis

Category	Product/Tool	Specific Application	Considerations
DNA Extraction	FastDNA SPIN Kit for Soil (MP Biomedicals) [78]	Biomass from sludge and granular systems	Effective for difficult-to-lyse environmental samples
DNA Extraction	ALFA-SEQ Advanced Water DNA Kit [76]	Low-biomass water samples	Optimized for filtration concentrates
Library Prep	Illumina DNA Prep Kit	Metagenomic sequencing	Compatible with various Illumina platforms
Sequencing	Illumina HiSeq 4000 [78]	High-throughput metagenomes	150-250 bp paired-end reads recommended
Quality Control	fastp [79]	Adapter trimming and quality filtering	Rapid processing of large datasets
Assembly	MEGAHIT, metaSPAdes	Metagenome assembly	Multiple k-mer strategies improve contiguity
Binning	MetaBAT2, MaxBin2	MAG reconstruction from assemblies	Ensemble approaches improve completeness
Taxonomy	GTDB-Tk [77]	Standardized taxonomic classification	Essential for novel organism identification
Functional Analysis	DIAMOND [79]	Fast sequence similarity searches	BLAST-compatible with improved speed
Pathway Analysis	KEGG, MetaCyc	Metabolic pathway reconstruction	Manual curation required for novel pathways
Network Analysis	Co-occurrence networks	Keystone taxon identification	Integrated with machine learning approaches [38]

Metagenome-Assembled Genomes have fundamentally transformed our capacity to identify and characterize keystone taxa within complex microbial consortia. The technical framework presented herein enables researchers to reconstruct metabolic networks and interaction patterns that sustain ecosystem functionality. In anammox systems, MAG applications have revealed that keystone metabolism centers on specialized nitrogen transformation, obligate cross-feeding of essential metabolites, and stress response coordination.

Future methodological advances will likely focus on long-read sequencing integration to improve MAG contiguity, single-cell genomics to resolve strain-level variation, and machine learning algorithms to predict keystone functions from genomic features. The integration of metabolic modeling with MAG data will further enhance our ability to predict community dynamics and engineer consortia with enhanced functions. As these tools mature, MAG-based keystone analysis will expand beyond anammox systems to diverse microbial habitats, advancing our fundamental understanding of microbiome assembly, stability, and function across engineered and natural ecosystems.

Global Patterns of Anammox and Denitrification Coupling Across Aquatic Ecosystems

Nitrogen overload poses a significant threat to aquatic ecosystems worldwide, driving eutrophication and associated environmental degradation. Within the nitrogen cycle, two microbial processes—anaerobic ammonium oxidation (anammox) and denitrification—collectively drive nitrogen loss from aquatic environments. While denitrification has long been recognized as the dominant pathway, the discovery and characterization of anammox have revolutionized our understanding of aquatic nitrogen cycling. These processes exhibit complex interactions, competing for substrates while potentially engaging in metabolic cooperation through microbial consortia.

This technical guide examines the global patterns of coupling between anammox and denitrification across diverse aquatic ecosystems, with particular emphasis on the role of keystone species in maintaining the stability and function of anammox bacterial communities. Understanding these patterns and mechanisms is crucial for accurately modeling global nitrogen budgets and developing innovative wastewater treatment technologies.

Global Patterns and Rates Across Aquatic Ecosystems

Comprehensive analysis of global datasets reveals distinct patterns in the distribution and activity of anammox and denitrification processes across aquatic ecosystems. A synthesis of 2539 observations from 136 peer-reviewed studies demonstrates that although these processes co-occur across diverse environments, their relative contributions to nitrogen loss vary substantially.

Table 1: Global Rates of Anammox and Denitrification in Aquatic Ecosystems

Ecosystem Type	Median Anammox Rate (nmol-N g⁻¹ day⁻¹)	Median Denitrification Rate (nmol-N g⁻¹ day⁻¹)	Typical Rana/den Ratio	Dominant Anammox Genera
Rivers	1471.38 ± 1366.09	968.67 ± 419.42	>0.5	Candidatus Brocadia, Candidatus Kuenenia
Lakes & Reservoirs	89.94	284.93	0.129 (global median)	Candidatus Jettenia, Candidatus Brocadia
Wetlands	21.55 (95% CI: 8.21–58.90)	171.76 (95% CI: 65.40–519.25)	0.129 (global median)	Candidatus Brocadia
Estuaries	1.92–264	171.76 (95% CI: 65.40–519.25)	Variable	Candidatus Scalindua
Marine Sediments	1.92–264	171.76 (95% CI: 65.40–519.25)	<0.1	Candidatus Scalindua

[10] [81] [82]

Globally, denitrification dominates nitrogen loss across most aquatic ecosystems, with a median rate of 171.76 nmol-N g⁻¹ day⁻¹ compared to 21.55 nmol-N g⁻¹ day⁻¹ for anammox. The median ratio of anammox to denitrification (Rana/den) stands at 0.129 globally, though notable exceptions exist where anammox exceeds denitrification, particularly in river ecosystems [10]. Inland aquatic ecosystems demonstrate significantly higher anammox rates (89.94–1471.38 nmol-N g⁻¹ day⁻¹) compared to terrestrial soils (38.4 nmol-N g⁻¹ day⁻¹) and marine sediments (1.92–264 nmol-N g⁻¹ day⁻¹), establishing them as critical hotspots for anammox-mediated nitrogen loss [10].

Spatial heterogeneity in process rates reflects ecosystem-specific conditions. River systems exhibit the highest anammox rates, potentially exceeding denitrification in some cases, while estuarine and marine environments show lower but still significant activity [10]. This spatial variation underscores the importance of local environmental conditions in regulating these microbial processes.

Table 2: Ecosystem-Specific Anammox Community Composition

Ecosystem Type	Dominant Anammox Genera	Relative Abundance	Key Environmental Drivers
Groundwater Aquifers	Candidatus Brocadia	80-99.9%	Stable anaerobic conditions, low disturbance
Rivers & Lakes	Candidatus Brocadia, Candidatus Jettenia	60-95%	Organic carbon, ammonium availability
Estuaries	Candidatus Scalindua	70-90%	Salinity gradient, nitrite availability
Marine Sediments	Candidatus Scalindua	>95%	Salinity, organic matter content
Wastewater Systems	Candidatus Brocadia, Candidatus Kuenenia	Variable	Nitrogen loading, temperature, reactor configuration

[24] [31] [82]

Keystone Species in Anammox Microbial Communities

Keystone species play disproportionately large roles in maintaining the structure and function of anammox bacterial communities, despite their often low relative abundance. These critical taxa include certain anammox bacteria themselves as well as associated microorganisms that facilitate the anammox process through metabolic interactions.

Among anammox bacteria, Candidatus Scalindua functions as a keystone genus in coastal and marine sediments, where it demonstrates extensive diversification and forms the backbone of microbial co-occurrence networks [24]. In terrestrial aquatic systems and engineered ecosystems, Candidatus Brocadia and Candidatus Kuenenia frequently occupy keystone roles, with the latter showing particular importance in biofilm formation and stability [4] [31].

Rare taxa contribute significantly to community stability despite their low abundance. These conditionally rare taxa exhibit high habitat specificity and become active under particular environmental conditions, providing functional resilience to the microbial community [24]. Network analysis has revealed that rare species frequently occupy central positions in anammox co-occurrence networks, potentially acting as stabilizers that maintain ecological structure under fluctuating conditions [24].

The interaction between anammox bacteria and denitrifying bacteria represents another dimension of keystone relationships. Certain denitrifying bacteria, including Thauera and Afipia, provide essential metabolites such as amino acids, cofactors, and vitamins to anammox bacteria, establishing cross-feeding relationships that enhance nitrogen removal efficiency [31]. This metabolic cooperation illustrates how keystone functions can emerge from interactions between different microbial groups.

Keystone Roles in Anammox Communities

Microbial Interactions and Coupling Mechanisms

The relationship between anammox and denitrification processes extends beyond simple competition to include complex cooperative interactions mediated by microbial consortia. Understanding these coupling mechanisms is essential for predicting nitrogen fluxes in natural ecosystems and optimizing engineered systems.

Competitive Interactions

Anammox and denitrifying bacteria compete for common substrates, particularly nitrite (NO₂⁻), which serves as a key electron acceptor for both processes [10]. This competition creates a dynamic balance where the dominance of either process depends on environmental conditions. Denitrifying bacteria typically outperform anammox bacteria in carbon-rich environments, as they can utilize organic matter more efficiently [10]. However, anammox bacteria gain a competitive advantage in carbon-limited conditions due to their autotrophic metabolism [31].

The competition extends beyond nitrite to include organic carbon sources. While anammox bacteria are primarily autotrophic, some species can oxidize short-chain fatty acids and may compete with denitrifiers for these electron donors [31]. This versatile metabolism allows certain anammox bacteria to function as facultative organotrophs under specific conditions, further blurring the functional boundaries between these processes.

Cooperative Relationships

Despite their competitive aspects, anammox and denitrification frequently exhibit cooperative relationships in aquatic ecosystems. Denitrifying bacteria can enhance anammox performance through the "nitrite loop," whereby they reduce nitrate (NO₃⁻) produced by the anammox reaction back to nitrite, thus providing additional substrate for anammox bacteria [4]. This cooperation is particularly important because the anammox reaction theoretically produces approximately 11% nitrate, which could otherwise accumulate and potentially inhibit the process [4].

Metabolic cross-feeding represents another form of cooperation, where denitrifying bacteria provide essential metabolites to anammox bacteria. Genomic analyses have revealed that dominant denitrifiers can supply amino acids, cofactors, and vitamins to anammox bacteria, creating symbiotic relationships that enhance the overall efficiency of nitrogen removal [31]. This cross-feeding is particularly important in biofilm systems, where physical proximity facilitates metabolite exchange [4].

Community Assembly and Network Interactions

The assembly of anammox bacterial communities is governed by both deterministic and stochastic processes, with deterministic factors (environmental selection) playing a stronger role in shaping community composition [82]. Ecological drift predominantly structures the overall anammox bacterial community in coastal sediments, while rare species are more susceptible to dispersal limitations and environmental selection [24].

Co-occurrence network analyses reveal minimal competition among anammox bacterial species, suggesting that environmental factors such as anaerobic conditions and ecosystem stability, rather than substrate limitation, are primary determinants of community structure [82]. These networks also identify anammox bacteria as keystone taxa, with Candidatus Scalindua frequently occupying central positions in marine sediments [24].

Research Methods and Experimental Protocols

Field Sampling and Sediment Core Collection

Standardized sampling protocols are essential for comparative studies of anammox and denitrification across aquatic ecosystems. Sediment cores should be collected using appropriate coring devices (e.g., gravity corers, box corers) based on water depth and sediment characteristics [24] [81]. For comprehensive analysis, cores should be sectioned at high-resolution intervals (1-4 cm depending on research objectives) to capture vertical stratification of microbial processes [24].

Pore water extraction should be performed immediately after sectioning using rhizons or centrifugation methods. Samples for molecular analysis should be preserved at -80°C, while those for process rate measurements should be processed fresh or stored under appropriate conditions [24] [81]. Parallel collection of water column samples provides context for understanding sediment-water exchanges.

Quantification of Process Rates

Slurry Assays for Denitrification and Anammox Potential:

• Prepare sediment slurries by homogenizing fresh sediment with degassed, artificial seawater or site water in gas-tight vials [81].
• Pre-incubate slurries overnight under anoxic conditions to consume residual oxygen, nitrate, and nitrite.
• Amend slurries with isotopically labeled substrates (¹⁵N-NH₄⁺, ¹⁵N-NO₃⁻, or ¹⁵N-NH₄⁺ + ¹⁴N-NO₃⁻) at environmentally relevant concentrations (typically 100 μM) [81].
• Sacrifice replicate vials at multiple time points (0, 0.5, 17, and 24 hours) by adding ZnClâ‚‚ to terminate biological activity.
• Analyze Nâ‚‚ production using gas chromatography-mass spectrometry (GC-MS) to distinguish between ²⁸Nâ‚‚, ²⁹Nâ‚‚, and ³⁰Nâ‚‚ [81].
• Calculate potential rates of denitrification and anammox based on the production of labeled Nâ‚‚ species over time.

Batch Tests for Partial Denitrification Kinetics:

• Set up batch tests with activated sludge at varying suspended solids concentrations (e.g., 40 mg-VSS/L) and carbon-to-nitrogen (C/N) ratios (e.g., 2-4) [83].
• Monitor nitrate, nitrite, and COD concentrations over time (typically 4-8 hours) to determine nitrite transformation rates (NTR) and nitrate removal efficiency (NRE).
• Determine kinetic parameters (ηNO₃, ηNOâ‚‚, KS1, KS2) by fitting data to the Activated Sludge Model No. 3 (ASM3) or modified versions [83].

Experimental Workflow for Coupling Studies

Molecular Analysis of Microbial Communities

DNA Extraction and Quantification:

• Extract total genomic DNA from 0.5 g of wet sediment using commercial kits (e.g., FastDNA SPIN Kit for Soil) following manufacturer's protocols [24].
• Assess DNA quality and quantity using fluorometric methods (e.g., Qubit fluorometer) and spectrophotometry (e.g., NanoDrop) [24].

Target Gene Amplification and Sequencing:

• Amplify anammox-specific 16S rRNA genes or functional genes (hzsB) using primers such as Brod541F/Amx820R [24].
• Amplify denitrification genes (nirS, nirK, nosZ) using appropriate primer sets [31].
• Perform high-throughput sequencing on amplified products using Illumina or similar platforms [24] [31].
• Process raw sequences through quality filtering, chimera removal, and OTU clustering at 97-98% similarity [24].

Metagenomic and Metatranscriptomic Analysis:

• Prepare sequencing libraries from extracted DNA/RNA using standard protocols.
• Perform shotgun sequencing to capture comprehensive functional potential.
• Reconstruct metagenome-assembled genomes (MAGs) to elucidate metabolic capabilities of key taxa [31].
• Annotate genes involved in nitrogen cycling and cross-feeding pathways.

Kinetic Modeling Approaches

Various kinetic models have been applied to describe anammox and denitrification processes:

Grau Second-Order Substrate Removal Model:

• Applied to both batch and continuous systems
• Effective for describing substrate removal at varying loading rates [84]

Modified Stover-Kincannon Model:

• Originally developed for attached-growth systems
• Adapted for anammox processes by replacing surface area with reactor volume [84]

Activated Sludge Model Framework:

• ASM1, ASM2d, and ASM3 have been modified to incorporate anammox processes
• Requires accurate determination of heterotrophic biomass proportions (e.g., 66.4% based on metagenomic sequencing) [83]

Process Kinetics Validation:

• Model selection should be based on statistical criteria (R², MSE)
• The Grau second-order and Stover-Kincannon models typically show best performance for anammox systems [84]

The Scientist's Toolkit: Essential Research Reagents and Materials

Table 3: Essential Research Reagents for Anammox and Denitrification Studies

Category	Specific Reagents/Materials	Application Purpose	Key Considerations
Molecular Biology	FastDNA SPIN Kit for Soil	DNA extraction from sediment samples	Effective for difficult environmental matrices
	Brod541F/Amx820R primers	Anammox 16S rRNA gene amplification	Specificity for anammox bacterial detection
	nirS/nirK primers	Denitrifier community analysis	Targets different nitrite reductase genes
Isotope Tracers	¹⁵N-NH₄⁺ (as chloride or sulfate salts)	Anammox process tracing	≥98% isotopic purity required
	¹⁵N-NO₃⁻ (as potassium or sodium salts)	Denitrification process tracing	Distinguish N₂ sources
Process Rate Measurements	ZnClâ‚‚ (50% w/v)	Biological reaction termination	Preserves Nâ‚‚ composition for GC-MS
	Degassed artificial seawater	Slurry assay medium	Controls salinity while removing interferents
	Exetainer vials (Labco Ltd)	Anaerobic incubations	Gas-tight for reliable Nâ‚‚ measurement
Analytical Standards	Certified N₂ gas mixtures (²⁸N₂, ²⁹N₂, ³⁰N₂)	GC-MS calibration	Essential for isotope pairing calculations
	Nutrient standards (NH₄⁺, NO₂⁻, NO₃⁻)	Ion chromatography calibration	Accurate concentration determination
Bioreactor Components	Polyurethane sponge carriers	Microbial biofilm support	High surface area for biomass retention
	Argon gas	Anoxic condition maintenance	Oxygen-free environment for anammox

[4] [24] [81]

Emerging Research Technologies and Approaches

Recent technological advances are transforming our ability to study anammox and denitrification coupling in aquatic ecosystems:

Machine Learning and Computer Vision:

• Deep learning models (e.g., ResNet50d neural network) can predict anammox sludge activity based on color features with high accuracy (R² = 0.984) [85].
• This non-destructive approach allows real-time monitoring of anammox process performance based on sludge coloration, which correlates with cytochrome c content and bacterial activity [85].

Metagenome-Assembled Genomes (MAGs):

• MAG reconstruction from metagenomic data enables metabolic modeling of uncultured anammox and denitrifying bacteria [31].
• These models reveal cross-feeding interactions and metabolic dependencies within microbial consortia [31].

Multi-Omics Integration:

• Combined analysis of metagenomics, metatranscriptomics, and metaproteomics provides comprehensive insights into functional potential, gene expression, and protein synthesis in complex communities [31].
• This approach helps identify active metabolic pathways and regulatory responses to environmental changes.

Enhanced Kinetic Modeling:

• Integration of microbial community data from metagenomic sequencing improves parameter estimation in kinetic models [83].
• Hybrid models incorporating both process kinetics and microbial dynamics offer more accurate predictions of system performance under varying conditions [83] [84].

The coupling of anammox and denitrification processes across aquatic ecosystems represents a complex interplay of competition and cooperation mediated by diverse microbial communities. Keystone species, including certain anammox bacteria and their associated denitrifiers, play critical roles in maintaining ecosystem function and stability. Understanding the global patterns and underlying mechanisms of this coupling is essential for accurate nitrogen budgeting and developing sustainable approaches to mitigate nitrogen pollution in an era of global environmental change.

Comparative Analysis of Keystone Species in Natural vs. Engineered Ecosystems

The concept of the keystone species, formally introduced by zoologist Robert T. Paine in 1969, describes an organism that exerts a disproportionately large influence on its natural environment relative to its abundance [37]. The removal of such a species triggers dramatic shifts in ecosystem structure and function, often leading to radically altered states or a collapse in biodiversity [86]. While this paradigm was founded on observations in natural ecosystems, the critical role of strongly interacting species extends into engineered environments designed for specific functions, such as wastewater treatment.

This review frames its analysis within the context of a broader thesis on keystone species in anaerobic ammonium oxidation (anammox) bacterial communities. It provides a comparative examination of the identity, ecological roles, and functional consequences of keystone species across natural ecosystems and engineered anammox systems, synthesizing foundational ecology with cutting-edge environmental biotechnology.

Foundational Concepts and Definitions

Keystone Species and Related Concepts

A keystone species is defined by its outsized ecological impact, which is not reflected by its biomass or abundance [86] [37]. The classic example is Paine's Pisaster ochraceus sea star, whose removal from a tidal plain led to a mussel monopoly and a halving of species biodiversity within a year [86]. Keystone species have low functional redundancy, meaning no other species can fill their ecological niche if they are lost [86].

Several related concepts are often discussed alongside keystone species:

• Ecosystem Engineers: Species that create, modify, maintain, or destroy habitats. Beavers are a quintessential example, building dams that transform streams into wetlands and create habitats for numerous other species [86] [87].
• Foundation Species: Species that play a major role in creating or maintaining a habitat itself, such as corals that build the physical structure of reefs [86].
• Umbrella Species: Species with large habitat needs, where protecting them consequently protects many other co-occurring species [86].
• Indicator Species: Organisms sensitive to environmental changes, providing early warning signs of ecosystem degradation [86].

While these roles can overlap, the defining feature of a keystone species is its disproportionate influence in maintaining ecosystem structure and biodiversity.

Keystone Species in Natural Ecosystems

Classic Paradigms and Mechanisms

In natural ecosystems, keystone species typically operate through top-down regulation, mutualistic interactions, or physical ecosystem engineering.

Table 1: Classic Keystone Species in Natural Ecosystems

Keystone Species	Ecosystem	Mechanism of Action	Impact of Removal
Gray Wolf (Canis lupus) [86]	Greater Yellowstone Ecosystem (Temperate Forest)	Top-down predation pressure regulating elk populations.	Elk overbrowsing; reduced willow and aspen; erosion of stream banks; decline of beaver and songbird populations.
Sea Otter (Enhydra lutris) [37]	North Pacific Kelp Forests	Predation on herbivorous sea urchins.	Urchin population explosion; overgrazing of kelp forests; collapse of kelp forest ecosystem.
Ochre Sea Star (Pisaster ochraceus) [86] [37]	Intertidal Zones	Predation on competitive dominant mussels.	Mussel monopoly; crowding out of other invertebrate and algal species; biodiversity cut in half.
African Elephant (Loxodonta spp.) [86] [37]	Savanna Grasslands	Herbivory and physical destruction of trees and shrubs.	Conversion of grassland to woodland or forest; loss of grazing areas for other herbivores.
Beaver (Castor canadensis) [86] [87]	Riparian Zones	Ecosystem engineering via dam building.	Loss of wetland habitats; reduced habitat complexity; changes in hydrology and nutrient cycling.
Green-backed Firecrown (Sephanoides sephanoides) [86]	Patagonian Woodlands	Keystone mutualist: pollinates 20% of local plant species.	Reproductive failure of dependent plants; collapse of pollinator-dependent habitat pockets.

Trophic Cascades and Ecosystem Stability

The case of the gray wolf in Yellowstone National Park provides a seminal example of a trophic cascade initiated by the loss of a keystone predator [86]. After wolves were extirpated in the early 20th century, elk populations exploded and overgrazed vegetation, which in turn affected physical geography, leading to eroded stream banks and increased water temperatures [86]. The reintroduction of wolves in the 1990s demonstrated the restorative power of a keystone species, leading to reduced elk pressure, recovery of willow and aspen, and the return of beaver and songbird populations [86]. This case underscores the role of keystone species in maintaining not only biodiversity but also fundamental abiotic processes.

Keystone Species in Engineered Anammox Ecosystems

The Anammox Process and its Microbial Community

Engineered anammox ecosystems are designed for the autotrophic removal of nitrogen from wastewater. The process is mediated by anammox bacteria, which are slow-growing, anaerobic chemolithoautotrophs belonging to the phylum Planctomycetes [4] [5]. They convert ammonium (NH₄⁺) and nitrite (NO₂⁻) directly into dinitrogen gas (N₂), offering a more energy-efficient alternative to conventional nitrification-denitrification [4] [88].

Anammox bioreactors are not monocultures but complex microbial consortia, termed the "anammox core" [4]. This community includes:

• Anammox bacteria (e.g., Candidatus Brocadia, Kuenenia, Jettenia): The primary drivers of nitrogen removal [4] [5].
• Ammonia-Oxidizing Bacteria (AOB) and Nitrite-Oxidizing Bacteria (NOB): Contribute to the substrate pool.
• Heterotrophic Denitrifying Bacteria (HDB) (e.g., Thauera, Afipia): Compete for substrates but also perform critical cross-feeding [4] [5].
• Other symbionts (e.g., Chloroflexi): With roles in aggregate structure and metabolic byproduct consumption [4].

Keystone Roles in the Anammox Consortium

In these engineered systems, the keystone function is not held by a single species but is an emergent property of specific, strongly interacting microbial guilds. The anammox bacteria themselves are foundational, but the stability and performance of the entire ecosystem rely on keystone interactions.

1. Anammox Bacteria as Keystone Mutualists: Anammox bacteria engage in cross-feeding and mutualistic interactions with denitrifying bacteria. Metagenome-assembled genomes have revealed that dominant denitrifiers can provide essential materials like amino acids, cofactors, and vitamins to anammox bacteria, supporting their growth and activity [5]. Conversely, some anammox bacteria, such as Ca. Kuenenia, might produce metabolites like acetate and glycogen that benefit other community members [4]. This cross-feeding creates a cooperative network that enhances the system's nitrogen removal efficiency.

2. Keystone Role of Microbial Interactions in Biofilm Formation: Microbial interactions are "keystone" to rapid sludge proliferation and biofilm formation [4]. The initial attachment of specific microbes to carriers and the subsequent production of Extracellular Polymeric Substances (EPS) create a scaffold for anammox bacteria and their symbionts to form robust biofilms [4]. These biofilms are crucial for protecting the slow-growing anammox bacteria from environmental perturbations.

3. Functional Keystone of the "Nitrite Loop": Heterotrophic denitrifiers play a keystone role by performing a "nitrite loop" [4]. Anammox metabolism produces nitrate (NO₃⁻), which can accumulate and inhibit the process. Denitrifiers reduce this nitrate back to nitrite, the essential substrate for anammox. This functional group maintains the stoichiometric balance of the ecosystem, preventing product inhibition and enhancing the overall Nitrogen Removal Efficiency (NRE).

Table 2: Keystone Functional Groups in Engineered Anammox Ecosystems

Functional Group	Role in the Consortium	Keystone Mechanism	Impact of Disruption
Anammox Bacteria (e.g., Ca. Brocadia, Ca. Kuenenia) [4] [5]	Primary N₂ production from NH₄⁺ and NO₂⁻.	Central metabolism driving the ecosystem's core function.	Collapse of nitrogen removal; process failure.
Heterotrophic Denitrifying Bacteria (e.g., Thauera, Afipia) [4] [5]	Reduces NO₃⁻ to NO₂⁻; consumes organic carbon.	"Nitrite Loop": Recycles inhibitory nitrate into anammox substrate.	Nitrate accumulation; reduced NRE; potential process imbalance.
Biofilm-Forming Microbes [4]	Secrete EPS and initiate biofilm formation on carriers.	Ecosystem Engineering: Creates a protected, structured habitat for the consortium.	Poor biomass retention; washout of slow-growing anammox bacteria; system instability.
nirs-type Denitrifiers [5]	Possess the nirS gene for nitrite reduction.	Strong Cooperative Coupling: Shows stronger metabolic coupling with anammox bacteria than nirK-type.	Weaker network stability; potential decrease in functional redundancy and resilience.

Comparative Analysis: Natural vs. Engineered Keystones

The following diagram illustrates the fundamental differences in the organization and keystone relationships within these two ecosystem types.

Key Differences and Similarities

Aspect	Natural Ecosystems	Engineered Anammox Ecosystems
Keystone Identity	Often a single, macroscopic animal species [86].	A functional guild of microscopic organisms; a network property [4] [5].
Primary Mechanism	Top-down predation, physical ecosystem engineering, mutualism [86] [37].	Metabolic cross-feeding, niche differentiation, and biofilm-mediated ecosystem engineering [4] [5].
Ecosystem Goal	Biodiversity, stability, energy flow [86].	Optimized, robust, and efficient nitrogen removal [4] [88].
Response to Stress	Trophic cascade, regime shift, biodiversity loss [86] [89].	Process inhibition (e.g., by Cr(VI)), microbial succession, network reorganization [88].
Redundancy	Low functional redundancy for keystone species [86].	Partial redundancy within functional guilds (e.g., multiple denitrifiers), but core anammox function is irreplaceable [4].
Commonality	Disproportionate impact on ecosystem structure and function relative to abundance.	Disproportionate impact on system performance and stability relative to abundance.

Experimental Approaches and Research Toolkit

The study of keystone species in both natural and engineered environments relies on a combination of observational, experimental, and molecular techniques.

Foundational Ecological Methods

• Manipulation Experiments: The most direct method, pioneered by Robert Paine, involves the removal or addition of a putative keystone species and monitoring the ecosystem-level response [86] [89]. In anammox systems, this can be simulated by inhibiting specific microbial groups.
• Long-Term Monitoring: Tracking population dynamics and ecosystem properties over time, as demonstrated in the 40-year study of the Yellowstone wolf reintroduction [86].
• Network Analysis: Using co-occurrence patterns to infer microbial interactions and identify keystone taxa based on network topology metrics like centrality [90] [89].

Molecular and "Omics" Tools for Microbial Consortia

Engineered ecosystem research heavily depends on advanced molecular tools to deconstruct the microbial black box.

Table 3: Research Reagent Solutions for Anammox Keystone Species Research

Reagent / Tool Category	Specific Examples & Targets	Function in Research
Molecular Probes & Primers	16S rRNA gene primers (general community); hzsB gene (anammox); nirS/K genes (denitrifiers) [5].	Quantifying and tracking the abundance and diversity of specific functional groups in the consortium.
Metagenomic Sequencing	Shotgun metagenomics of total DNA from biofilm or flocs [4] [5].	Revealing the total genetic potential and reconstructing Metagenome-Assembled Genomes (MAGs) to infer metabolic pathways.
Activity Assays	Hydrazine Dehydrogenase (HDH) activity; Heme c content; Specific Anammox Activity (SAA) [88].	Directly measuring the metabolic activity and health of the anammox bacteria under different conditions.
Stable Isotope Probing	¹⁵N-labeled ammonium (¹⁵NH₄⁺) or nitrite (¹⁵NO₂⁻) [5].	Tracing the flow of nitrogen through the different microbial groups, confirming metabolic pathways and cross-feeding.
Chemical Inhibitors	Allylthiourea (inhibits AOB); Cr(VI) stressor [88].	Experimentally manipulating specific functional groups to infer their role and identify keystone interactions.

The following workflow outlines a typical integrated protocol for analyzing keystone interactions in an anammox system.

Detailed Experimental Protocol for Anammox Reactor Operation:

• Reactor Setup: Use an Up-flow Anaerobic Sludge Blanket (UASB) reactor with a working volume of 1.0-1.5 L, filled with sponge or other porous carriers to support biofilm formation [4] [88].
• Inoculation: Seed with mature anammox sludge. The initial abundance of anammox bacteria (e.g., Ca. Kuenenia) in the seed sludge is a critical variable affecting subsequent community structure and function [4].
• Feed Composition: Use synthetic wastewater containing ammonium (NH₄⁺, e.g., 100 mg/L as N) and nitrite (NO₂⁻, e.g., 130 mg/L as N), along with essential minerals and bicarbonate as an inorganic carbon source [88].
• Operational Control: Maintain strict anoxic conditions (Dissolved Oxygen < 0.1 mg/L), temperature at 34 ± 1°C or room temperature, and a controlled Hydraulic Retention Time (HRT), typically between 5-48 hours depending on the study phase [4] [5].
• Monitoring & Sampling: Continuously monitor influent and effluent concentrations of NH₄⁺, NO₂⁻, and NO₃⁻ to calculate Nitrogen Removal Efficiency (NRE) and Rate (NRR). Periodically sample biomass (both biofilm and suspended flocs) for molecular and metabolic activity analysis [4] [88].

This comparative analysis reveals that while the scale and constituents differ profoundly, the fundamental principle of the keystone species—a disproportionate impact on ecosystem stability and function—holds true in both natural and engineered worlds. The gray wolf exerts top-down control to maintain biodiversity in Yellowstone, while in an anammox reactor, a metabolic guild of denitrifiers performs a critical "nitrite loop" to maintain functional efficiency. The shift from a single macroscopic predator to a network of microscopic mutualists and engineers reflects the adaptation of ecological principles to human-designed systems.

Understanding engineered ecosystems like anammox reactors through the lens of keystone species and ecological networks provides a powerful framework for optimizing their design and operation. It moves the focus from a single bacterium to the synergistic interactions that underpin system resilience and performance. This ecological perspective is indispensable for advancing the next generation of robust, self-regulating environmental biotechnologies.

Physiological and Genomic Signatures Validating Keystone Status

The identification of keystone species within anaerobic ammonium oxidation (anammox) bacterial communities represents a critical frontier in microbial ecology and environmental biotechnology. This whitepaper synthesizes contemporary research to delineate how integrated physiological and genomic signatures validate the keystone status of specific anammox taxa. Through multi-omics approaches and advanced community analysis, we demonstrate that keystone species—including both abundant core organisms and critical rare taxa—disproportionately influence community stability, metabolic networking, and ecosystem functionality. These findings provide a mechanistic framework for predicting community assembly and designing targeted interventions in both engineered and natural ecosystems.

The conceptual framework of keystone species, originally developed in macro-ecology, has profound implications for understanding the structure and function of anammox consortia. In the context of anammox systems, keystone species are defined as those whose impact on community structure, stability, and metabolic function is disproportionately large relative to their abundance [11]. These species operate through multiple mechanisms: they may act as metabolic pioneers that establish critical biogeochemical pathways, form architectural scaffolds for biofilm development, or generate essential metabolites that cross-feed auxiliary community members. The remarkable functional resilience observed in anammox systems facing environmental perturbations—such as nitrite inhibition, salinity fluctuations, and temperature variations—is increasingly attributed to these keystone taxa rather than merely to overall community diversity [33] [4].

Anammox-based nitrogen removal technologies have emerged as sustainable alternatives to energy-intensive conventional processes, with full-scale installations now exceeding 100 globally [91] [9]. However, the practical implementation of anammox processes continues to face challenges related to slow startup times, sensitivity to environmental perturbations, and unpredictable community dynamics. Resolving the keystone status of specific community members provides a mechanistic basis for addressing these limitations through targeted bioaugmentation strategies, optimized reactor inoculations, and predictive modeling of community assembly. This whitepaper integrates findings from genomic, physiological, and ecological studies to establish a definitive signature-based framework for identifying and validating keystone species in anammox communities.

Physiological Signatures of Keystone Status

Metabolic Activity and Environmental Plasticity

Keystone anammox species exhibit distinctive physiological signatures that enable them to exert disproportionate influence on community function. Among the most significant is their ability to maintain metabolic activity across fluctuating environmental conditions while supporting the functional stability of the broader community. Experimental evidence demonstrates that specific anammox taxa display remarkable plasticity in their responses to environmental stressors such as nitrite inhibition, with certain species capable of maintaining metabolic activity at nitrite concentrations as high as 280 mg-N/L while others are completely inhibited [33]. This differential tolerance creates ecological opportunities for keystone species to stabilize community function during perturbation events.

The anammox bacterium Candidatus Brocadia sinica exemplifies this signature, maintaining optimal anaerobic ammonium oxidation activity (SAA) within pH ranges of 7.0–8.0 and temperature ranges of 30–40°C, while other species like Ca. Brocadia sapporiensis exhibit adaptation to lower pH ranges [91]. This niche specialization enables complementary functioning under variable reactor conditions. Furthermore, keystone species demonstrate superior resilience mechanisms, with documented recovery of nitrogen removal rates from 0.16 ± 0.04 kg-N/(m³·d) to 0.29 ± 0.06 kg-N/(m³·d) following the cessation of nitrite stress at 200 mg-N/L [33]. This recovery capacity often precedes and facilitates the restoration of broader community function, highlighting the disproportionate role of these taxa in maintaining ecosystem stability.

Structural Specialization and Biofilm Formation

Distinctive cellular architectures represent another key physiological signature of anammox keystone species. The unique crateriform structures observed on anammox cell surfaces—uniformly distributed with diameters of approximately 5 nm—facilitate enhanced material exchange and communication between microorganisms [92]. These specialized structures are genetically stable features that potentially originate from flagellar degeneration and function as conduits connecting the cytoplasmic membrane, outer membrane, and anammoxosome membrane. This structural adaptation enables efficient substrate transport and metabolic coordination within the community.

Keystone species further demonstrate pronounced capabilities in biofilm formation and architectural scaffolding. Research indicates that anammox bacteria preferentially enrich in biofilm regions, where they achieve higher relative abundances compared to suspended flocs due to reduced mass transfer resistance and enhanced protection from environmental stressors [4]. The biofilm formation process follows a specific ecological succession wherein pioneer species initially colonize carrier surfaces, followed by the growth of anammox bacteria and their symbiotic partners supported by extracellular polymeric substances (EPS) secreted by these foundational colonizers [4]. This structured development creates heterogeneous microenvironments that support metabolic division of labor and cross-feeding relationships essential for community stability.

Table 1: Physiological Signatures of Keystone Anammox Species

Signature Type	Specific Characteristics	Functional Significance	Representative Taxa
Metabolic Plasticity	Maintains activity at NO₂⁻ concentrations up to 280 mg-N/L; Recovers from stress (0.16 to 0.29 kg-N/(m³·d) NRR)	Stabilizes community function during perturbation	Candidatus Jettenia, Candidatus Brocadia
Niche Specialization	pH optimum 7.0-8.0 (Ca. B. sinica); Lower pH adaptation (Ca. B. sapporiensis)	Enables functional complementarity	Candidatus Brocadia sinica, Candidatus Brocadia sapporiensis
Structural Adaptation	Crateriform structures (5 nm diameter); Anammoxosome compartment (38.9%-66% cell volume)	Enhances substrate transport and metabolic efficiency	All anammox bacteria
Biofilm Architecture	Preferential enrichment in biofilm; EPS-mediated colonization	Provides environmental resistance and facilitates symbiosis	Candidatus Kuenenia

Genomic and Molecular Signatures

Genomic Features and Adaptive Gene Sets

The keystone status of anammox bacteria is strongly supported by distinctive genomic signatures that enable niche specialization and functional dominance. Comparative genomic analyses reveal that salt-adapted species like Candidatus Scalindua possess unique gene sets absent in non-halophilic relatives, including specialized sodium ion transporters (Na+/H+ antiporters) and distinct amino acid composition biases with significant enrichment of acidic residues [91] [9]. These genomic adaptations represent evolutionary responses to environmental pressures and contribute to the functional specialization that underpins keystone roles in specific habitats.

Pan-genome analyses further illuminate the genetic foundations of keystone functionality, demonstrating that only 8.1% of genes constitute the core genome shared across anammox species, while a substantial proportion of genes are accessory and potentially contribute to niche specialization [91]. This genetic diversity facilitates the establishment of keystone taxa through mechanisms including horizontal gene transfer and environmental selection. Phylogenetic evidence indicates that the last common ancestor of contemporary anammox bacteria existed approximately 2.1 billion years ago, providing an extensive evolutionary timeline for the development of these specialized genomic signatures [91]. The relatively small core genome highlights the importance of accessory genes in establishing the metabolic versatility that characterizes keystone taxa.

Metabolic Network Topology and Keystone Genes

Beyond individual genomic features, the topological properties of metabolic networks provide compelling signatures of keystone status. Integrated multi-omic analyses—incorporating metagenomics, metatranscriptomics, and metaproteomics—enable the reconstruction of community-wide metabolic networks where specific genes occupy strategically important positions [93]. Genes exhibiting both high relative expression (transcript and protein abundance relative to gene copy number) and high betweenness centrality within metabolic networks are operationally defined as "keystone genes" [93]. These genes encode enzymes that occupy critical choke points in community metabolism and disproportionately influence overall ecosystem function.

Empirical studies have identified several keystone gene categories within anammox consortia, including those involved in glycerolipid metabolism (particularly triacylglycerol lipase), nitrogen metabolism, and fatty acid biosynthesis [93]. The expression of these genes is often linked to known lipid-accumulating populations such as Candidatus Microthrix parvicella and nitrogen-transforming organisms like Nitrosomonas spp. and Rhodococcus spp. [93]. The identification of these genetically encoded critical functions enables researchers to move beyond taxonomy-based community profiling toward a mechanistic understanding of how specific metabolic capabilities influence community stability and ecosystem function.

Table 2: Genomic and Molecular Signatures of Keystone Anammox Species

Signature Category	Specific Features	Analytical Method	Functional Interpretation
Adaptive Gene Sets	Na+/H+ antiporters; Acidic amino acid bias; Unique gene complement	Comparative genomics	Environmental specialization (e.g., halotolerance)
Genomic Architecture	8.1% core genome; Extensive accessory genome	Pan-genome analysis	Metabolic versatility and niche adaptation
Metabolic Network Features	High betweenness centrality; High relative expression (transcript/protein to gene copy ratio)	Integrated multi-omics	Identification of "keystone genes" critical to community function
Key Gene Categories	Glycerolipid metabolism (triacylglycerol lipase); Nitrogen metabolism; Fatty acid biosynthesis	Metatranscriptomics & metaproteomics	Rate-limiting steps in community metabolism

Experimental Methodologies for Keystone Validation

Community Manipulation and Perturbation Studies

Experimental validation of keystone status requires carefully designed community manipulation approaches that assess functional responses to the presence, absence, or reduction of putative keystone taxa. Long-term nitrite stress experiments—conducted over 588 days in parallel biofilm reactors—provide a robust methodological framework for identifying keystone species through their differential responses to environmental perturbation [33]. These studies employ reactor systems with identical seeding sludge under controlled conditions, with systematic monitoring of nitrogen removal efficiency, microbial community dynamics, and functional resilience.

The experimental workflow involves phased operational strategies, beginning with baseline performance establishment without stress, followed by progressive increases in nitrite loading (e.g., 0-200 mg-N/L), and concluding with a recovery phase after stress termination [33]. Throughout these phases, integrated sampling for chemical analysis (NH₄⁺, NO₂⁻, NO₃⁻ concentrations), extracellular polymeric substances (EPS) characterization, and biomolecular extraction (DNA, RNA) enables correlation of process performance with community structural and functional changes. This approach successfully identified a shift from Candidatus Jettenia to Candidatus Brocadia as a key indicator of nitrite stress response, with the latter demonstrating keystone resilience properties [33].

Multi-Omic Integration and Network Analysis

The comprehensive identification of keystone species necessitates integrated multi-omic approaches that resolve community members alongside their functional contributions. A validated methodological framework begins with coordinated biomolecular extraction from community samples, sequentially recovering DNA, RNA, and proteins from a single sample to maintain functional linkage [93]. High-throughput sequencing of DNA (metagenomics) and RNA (metatranscriptomics), coupled with mass spectrometry-based proteomics, generates complementary data layers that collectively illuminate the functional capacity, expression patterns, and catalytic activity of community members.

Bioinformatic processing involves quality-controlled assembly of sequencing reads, prediction of protein-coding genes, and non-redundant cataloging of gene products [93]. Functional annotation against specialized databases (e.g., Kyoto Encyclopedia of Genes and Genomes) supports reconstruction of community-wide metabolic networks, which are subsequently analyzed to identify genes with disproportionately important roles through their network topology (high betweenness centrality) and expression patterns (high relative expression) [93]. This integrated methodology enables researchers to move beyond correlation-based co-occurrence networks toward mechanistic identification of keystone functions and the taxa that encode them.

Co-occurrence Network Analysis and Community Assembly Assessment

Network-based analytical approaches provide powerful methodological tools for identifying keystone species through their patterns of interaction within microbial communities. Co-occurrence network analysis begins with high-throughput 16S rRNA gene amplicon sequencing of anammox communities across environmental gradients or temporal series, followed by rigorous data preprocessing to account for sequencing artifacts and compositional effects [11]. Statistical correlation measures (e.g., SparCC, SPIEC-EASI) that address compositional data challenges are applied to derive robust co-occurrence patterns, which are then visualized as network graphs where nodes represent operational taxonomic units (OTUs) and edges represent significant positive or negative associations.

Within these co-occurrence networks, keystone species are identified through topological metrics including betweenness centrality (measuring a node's role as a connector), degree (number of connections), and closeness centrality (proximity to other nodes) [11]. Complementary analyses assess community assembly mechanisms by quantifying the relative contributions of deterministic processes (like environmental selection) and stochastic processes (like ecological drift) through null modeling approaches [11]. This methodological framework has revealed that rare anammox taxa are particularly susceptible to dispersal limitations and environmental selection, while abundant taxa exhibit greater dispersal capability, with both potentially serving complementary keystone roles in different contexts [11].

The Scientist's Toolkit: Research Reagent Solutions

Table 3: Essential Research Reagents and Methodologies for Keystone Species Identification

Reagent/Method	Specific Application	Function in Keystone Species Research
Primer Set Brod541F/Amx820R	16S rRNA gene amplification	Specific detection and quantification of anammox bacteria in complex communities [11]
FastDNA SPIN Kit for Soil	DNA extraction from complex samples	Efficient biomolecular recovery from biofilm and sediment matrices [11]
Qiagen AllPrep DNA/RNA/Protein Mini Kit	Coordinated multi-omic extraction	Simultaneous recovery of DNA, RNA, and proteins from single samples for functional correlation [93]
Kyoto Encyclopedia of Genes and Genomes (KEGG)	Functional annotation	Metabolic pathway mapping and network reconstruction [93]
QIIME 2 Pipeline	Community sequence analysis	Processing of amplicon sequencing data for diversity and network analysis [11]
Non-metric Multidimensional Scaling (NMDS)	Beta-diversity visualization	Spatial heterogeneity analysis of anammox communities across environments [11]
Metabolic Network Betweenness Centrality	Topological network analysis	Identification of critically positioned genes in community metabolic networks [93]

Ecological Implications and Applications

Keystone Roles in Natural Ecosystems

In natural environments, keystone anammox species perform critical ecosystem functions that extend beyond their nitrogen-cycling roles. Research across estuarine and coastal systems—including the Changjiang Estuary, Oujiang Estuary, Jiulong River Estuary, and South China Sea—has demonstrated that Candidatus Scalindua functions as a keystone genus in marine sediments, where it dominates the anammox community and influences nitrogen loss patterns [11]. The ecological dominance of this genus in marine environments highlights how keystone taxa can shape biogeochemical cycling across ecosystem boundaries.

The stability of natural anammox communities is maintained through a complex interplay between abundant and rare keystone species. While abundant taxa like Candidatus Scalindua contribute significantly to nitrogen flux, rare taxa provide critical functional resilience through mechanisms that include rapid response to changing conditions and maintenance of functional diversity [11]. Network analyses have revealed that rare species frequently occupy positions that enhance network connectivity and stability, potentially serving as "insurance" populations that can expand under favorable conditions to maintain ecosystem function [11]. This dual-keystone strategy, incorporating both high-abundance core species and low-abundance but functionally critical rare species, represents a sophisticated ecological adaptation that ensures functional stability across fluctuating environmental conditions.

Biotechnological Applications and Engineering Design

The identification and validation of keystone species in anammox communities enables targeted bioengineering strategies for wastewater treatment applications. Inoculation strategies that leverage keystone principles have demonstrated significantly enhanced startup times and process stability. Experimental evidence shows that reactors inoculated with communities containing higher relative abundances of keystone Candidatus Kuenenia achieved superior biomass accumulation (biofilm and flocs) and nitrogen removal performance compared to those with lower initial keystone abundance [4]. This effect stems from the enhanced expression of critical functional genes including hydrazine synthase (hzs) and nitrite reductase (nir), alongside key carbon metabolism genes (fdh, glgA/B/C, acs) that support the production of acetate and glycogen for cross-feeding within the community [4].

Carrier-based biofilm systems represent another application of keystone principles, leveraging the natural propensity of specific anammox taxa to form architectural scaffolds that support broader community development [4]. The strategic design of carrier surfaces to promote the initial attachment and growth of keystone pioneer species can accelerate the establishment of functional anammox biofilms. Furthermore, operational strategies that maintain environmental conditions favorable to keystone taxa—including controlled nitrite levels, optimal pH ranges, and temperature stability—can promote the functional dominance of these critical community members [33]. These engineered applications demonstrate how mechanistic understanding of keystone species can translate to improved biotechnological performance.

The validation of keystone status in anammox communities through integrated physiological and genomic signatures represents a paradigm shift in microbial ecology and environmental biotechnology. The signatures detailed in this whitepaper—including metabolic plasticity, structural specialization, genomic adaptation, and network topological properties—provide a robust framework for identifying taxa that disproportionately influence community structure and function. The methodological advances in multi-omic integration, community perturbation studies, and network analysis now enable researchers to move beyond correlative relationships toward mechanistic understanding of keystone functions.

Future research directions should prioritize the development of dynamic models that incorporate keystone concepts to predict community responses to environmental change. Additionally, systematic exploration of rare biosphere members promises to reveal novel keystone taxa with potentially unique functional capabilities. Translation of these ecological insights into engineering design principles will accelerate the implementation of robust anammox-based technologies for sustainable wastewater treatment. As the resolution of multi-omic methodologies continues to improve, so too will our capacity to identify and leverage keystone species for enhanced ecosystem function and biotechnological innovation.

Conclusion

The identification and understanding of keystone species within anammox communities represent a paradigm shift in managing microbial ecosystems for enhanced nitrogen removal. These pivotal taxa, often rare yet disproportionately influential, dictate community structure, functional stability, and resilience to environmental stressors. The integration of top-down identification frameworks with network analysis and genomic validation provides a powerful toolkit for pinpointing these critical organisms across diverse ecosystems. Future research must focus on elucidating the specific metabolic interactions that underpin keystone functionality, developing precise manipulation strategies for these taxa in biotechnological applications, and scaling insights from bioreactors to predict ecosystem responses to global change. Harnessing the power of anammox keystone species promises more stable, efficient, and predictable nitrogen removal systems, with significant implications for wastewater treatment, environmental remediation, and our fundamental understanding of microbial ecology.

References

• 1. Top-down identification of keystone taxa in the microbiome [https://pmc.ncbi.nlm.nih.gov/articles/PMC10319726/]
• 2. Characterization of two keystone taxa, sulfur-oxidizing, and ... [https://microbiomejournal.biomedcentral.com/articles/10.1186/s40168-023-01583-1]
• 3. Keystone taxa regulate microbial assemblage patterns and ... [https://www.sciencedirect.com/science/article/abs/pii/S0960852419310089]
• 4. Microbial interactions play a keystone role in rapid ... [https://www.sciencedirect.com/science/article/abs/pii/S0960852423010404]
• 5. Ecological interactions and the underlying mechanism of ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC10116762/]
• 6. Dynamic Response Mechanisms of Anammox Reactors ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12029708/]
• 7. Quantitative effect of adding percentages of anammox ... [https://www.sciencedirect.com/science/article/abs/pii/S0301479723021497]
• 8. Ecological niche differentiation among anammox bacteria [https://www.sciencedirect.com/science/article/abs/pii/S004313542030004X]
• 9. Genomic and Physiological Signatures of Evolution in ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12445922/]
• 10. Global patterns and drivers of coupling between anammox ... [https://www.nature.com/articles/s43247-024-01980-w]
• 11. Critical roles of rare species in the anaerobic ammonium ... [https://link.springer.com/article/10.1007/s42995-025-00315-8]
• 12. A comprehensive insight into the abundance and ... [https://www.sciencedirect.com/science/article/abs/pii/S0025326X23013504]
• 13. Long-term adaptation of two anammox granules with ... [https://www.sciencedirect.com/science/article/abs/pii/S0048969722075969]
• 14. Nitric oxide‐forming nitrite reductases in the anaerobic ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12485887/]
• 15. Fractionating proteins with nitrite-reducing activity in “ ... [https://www.frontiersin.org/journals/microbiology/articles/10.3389/fmicb.2025.1483703/full]
• 16. Experimental validation of the Monod growth model with ... [https://www.sciencedirect.com/science/article/pii/S0043135424017834]
• 17. Metabolic Diversity and Aero-Tolerance in Anammox ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC8862662/]
• 18. New Insight Into the Interspecies Shift of Anammox ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC8851195/]
• 19. Dynamic Response Mechanisms of Anammox Reactors ... [https://www.mdpi.com/2076-2607/13/4/899]
• 20. Characterization of a Microbial Community in an Anammox ... [https://www.mdpi.com/2073-4441/9/11/829]
• 21. Niche differentiation of anaerobic ammonium-oxidizing ... [https://pubmed.ncbi.nlm.nih.gov/41059904/]
• 22. Niche Differentiation of Aerobic and Anaerobic Ammonia ... [https://www.frontiersin.org/journals/microbiology/articles/10.3389/fmicb.2019.02141/full]
• 23. Bacterial community structure of a lab-scale anammox ... [https://pubmed.ncbi.nlm.nih.gov/25270790/]
• 24. Critical roles of rare species in the anaerobic ammonium ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12413392/]
• 25. Metabolic network analysis reveals microbial community ... [https://www.nature.com/articles/ncomms15416]
• 26. Establishing and maintaining stable anammox pathway for ... [https://www.sciencedirect.com/science/article/abs/pii/S1385894724105633]
• 27. Top-down identification of keystone taxa in the microbiome [https://pubmed.ncbi.nlm.nih.gov/37402745/]
• 28. Anammox [https://en.wikipedia.org/wiki/Anammox]
• 29. Functional stability correlates with dynamic microbial ... [https://www.sciencedirect.com/science/article/abs/pii/S0960852422018909]
• 30. The interaction strength of keystone module in cross ... [https://www.sciencedirect.com/science/article/abs/pii/S0929139325000447]
• 31. Ecological interactions and the underlying mechanism of ... [https://microbiomejournal.biomedcentral.com/articles/10.1186/s40168-023-01532-y]
• 32. Top-down identification of keystone taxa in the microbiome [https://www.nature.com/articles/s41467-023-39459-5]
• 33. The counteraction of anammox community to long-term ... [https://www.sciencedirect.com/science/article/abs/pii/S0048969722001528]
• 34. Advances in methods for detection of anaerobic ammonium ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC3082692/]
• 35. Carrier type affects anammox community assembly ... [https://www.sciencedirect.com/science/article/abs/pii/S0960852422017552]
• 36. Top-down and bottom-up microbiome engineering ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC11287213/]
• 37. Keystone species [https://en.wikipedia.org/wiki/Keystone_species]
• 38. Specialized metabolic functions of keystone taxa sustain soil ... [https://microbiomejournal.biomedcentral.com/articles/10.1186/s40168-020-00985-9]
• 39. The characteristics of soil microbial co-occurrence ... [https://www.frontiersin.org/journals/microbiology/articles/10.3389/fmicb.2023.1160683/full]
• 40. Specialized metabolic functions of keystone taxa sustain soil ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC7849160/]
• 41. Anammox-based technologies: A review of recent ... [https://www.sciencedirect.com/science/article/abs/pii/S1001074224000214]
• 42. Carrier type induces anammox biofilm structure and the ... [https://www.sciencedirect.com/science/article/abs/pii/S0960852421005885]
• 43. Comparative Investigation of the Anammox Process Using ... [https://www.mdpi.com/2073-4441/16/23/3363]
• 44. Controlling anammox speciation and biofilm attachment ... [https://www.nature.com/articles/s41598-022-26069-2]
• 45. Adjusted bacterial cooperation in anammox ... - PubMed Central [https://pmc.ncbi.nlm.nih.gov/articles/PMC11460484/]
• 46. Anammox activity improved significantly by the cross-fed ... [https://www.sciencedirect.com/science/article/abs/pii/S0043135423014264]
• 47. Metabolic cross-feeding of a dietary antioxidant enhances ... [https://pubmed.ncbi.nlm.nih.gov/40763732/]
• 48. Gut bacteria cooperate to convert dietary compound into ... [https://microbiomepost.com/gut-bacteria-cooperate-to-convert-dietary-compound-into-energy-boosting-microbial-health/]
• 49. Combined partial denitrification/anammox process for ... [https://www.sciencedirect.com/science/article/abs/pii/S2213343722018516]
• 50. Nitrogen removal performance and influencing factors of in ... [https://www.sciencedirect.com/science/article/abs/pii/S2214714422005943]
• 51. Synergistic enhancement of nitrogen removal in partial ... [https://www.sciencedirect.com/science/article/abs/pii/S2213343725048961]
• 52. Pilot-scale implementation of mainstream anammox for ... [https://www.nature.com/articles/s41467-024-54805-x]
• 53. Reason and control strategy for denitrification ... [https://www.sciencedirect.com/science/article/abs/pii/S0013935124013616]
• 54. Co-Occurring Anammox, Denitrification, and ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC3536082/]
• 55. Adjusted bacterial cooperation in anammox community to ... [https://www.sciencedirect.com/science/article/pii/S2589914724000483]
• 56. Oxygen tolerance and detoxification mechanisms of highly ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC10156729/]
• 57. Interspecies cooperation, degradation and self-adaptation ... [https://pubmed.ncbi.nlm.nih.gov/41192179/]
• 58. Regulation mechanism and resistant strain based ... [https://www.sciencedirect.com/science/article/abs/pii/S1385894725085857]
• 59. Comparative evaluation of multiple methods to quantify ... [https://www.sciencedirect.com/science/article/abs/pii/S0043135414007003]
• 60. Advancements and Challenges in Carrier-Mediated ... [https://www.sciencedirect.com/science/article/pii/S0043135425017762]
• 61. Enhance adhesion of ammonia-oxidizing bacteria on ... [https://sustainenvironres.biomedcentral.com/articles/10.1186/s42834-025-00263-z]
• 62. Effect of Mineral Carriers on Biofilm Formation and ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC9598201/]
• 63. Effect of Mineral Carriers on Biofilm Formation and ... [https://www.mdpi.com/2079-7737/11/10/1421]
• 64. Start-Up and Bacterial Enrichment of an Anammox Reactor ... [https://www.mdpi.com/2073-4441/16/15/2116]
• 65. Enhancement of mainstream Anammox by the composite ... [https://www.sciencedirect.com/science/article/abs/pii/S1385894725077605]
• 66. Biofilm colonization and succession in a full-scale partial ... [https://microbiomejournal.biomedcentral.com/articles/10.1186/s40168-024-01762-8]
• 67. Weighting and indirect effects identify keystone species in ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC5008267/]
• 68. Analysis of keystone species in a quantitative network ... [https://www.sciencedirect.com/science/article/abs/pii/S1476945X24000205]
• 69. Preservation and reactivation of anammox biomass: A mini ... [https://www.sciencedirect.com/science/article/abs/pii/S2213343724001945]
• 70. A review of anammox metabolic response to environmental ... [https://pubmed.ncbi.nlm.nih.gov/36773633/]
• 71. Abundant and rare microbial sub-communities in anammox ... [https://www.sciencedirect.com/science/article/abs/pii/S0960852418308010]
• 72. Micro-aeration and low influent C/N are key environmental ... [https://www.sciencedirect.com/science/article/abs/pii/S0043135423005778]
• 73. Controlling anammox speciation and biofilm attachment strategy using... [https://www.nature.com/articles/s41598-022-26069-2?error=cookies_not_supported&code=07a85bc5-d706-45e9-a3b7-5a96070d23ac]
• 74. Environmental instability reduces shock resistance by ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC11958701/]
• 75. Characterization of the Spatial Variation of Microbial ... [https://www.mdpi.com/2073-4441/13/9/1157]
• 76. Metagenomics-assembled genomes reveal microbial ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12174138/]
• 77. Metagenomic evidence of a novel family of anammox ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC9325076/]
• 78. Exploring the effects of operational mode and microbial ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC6714388/]
• 79. Metagenomic insights into the symbiotic relationship in ... [https://www.sciencedirect.com/science/article/abs/pii/S0043135422011290]
• 80. Metagenomic analysis of anammox communities in three ... [https://pubmed.ncbi.nlm.nih.gov/26568531/]
• 81. Denitrification and Anammox in Tropical Aquaculture ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC3433476/]
• 82. Global Distribution of Anaerobic Ammonia Oxidation ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC6861858/]
• 83. Determination of partial denitrification kinetic model ... [https://www.sciencedirect.com/science/article/abs/pii/S0960852423004030]
• 84. Investigation of the process stability of different anammox ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC9057419/]
• 85. Research on the Denitrification Efficiency of Anammox ... [https://www.mdpi.com/2073-4441/17/14/2084]
• 86. Role of Keystone Species in an Ecosystem [https://education.nationalgeographic.org/resource/role-keystone-species-ecosystem/]
• 87. Ecosystem Engineers: Creating Habitat for Others [https://www.oriannesociety.org/science-of-scales/ecosystem-engineers-creating-habitat-for-others/]
• 88. Deciphering the microbial patterns of anammox process ... [https://www.sciencedirect.com/science/article/abs/pii/S0304389421008141]
• 89. 5.4 Keystone Species and Ecosystem Engineers [https://fiveable.me/conservation-biology/unit-5/keystone-species-ecosystem-engineers/study-guide/uWWU8SHowQgpHPsF]
• 90. Changes in keystone species attenuate the complexity and ... [https://www.sciencedirect.com/science/article/pii/S1470160X25012245]
• 91. 综述：厌氧氨氧化细菌的基因组与生理特征演化印记 [https://www.ebiotrade.com/newsf/2025-9/20250920011800361.htm]
• 92. 厌氧氨氧化菌火山口结构 [https://actamicro.ijournals.cn/html/actamicrocn/2023/1/20230106.htm]
• 93. Comparative integrated omics: identification of key ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC5515219/]