Primer Selection and Design for Nitrogen Cycle Genes – ddPCR

SOP-L701 Revision 0.2
Author

Marc Los Huertos

EA Program, Pomona College

Published

2026-06-26

1 Purpose

This procedure describes the selection, design, and validation of primers and probes targeting functional genes of the nitrogen cycle for use with droplet digital PCR (ddPCR). The nitrogen cycle genes covered include those encoding nitrification (amoA for both ammonia-oxidizing bacteria and archaea), denitrification (nirK, nirS, nosZ), nitrogen fixation (nifH), and dissimilatory nitrate reduction to ammonium (nrfA). The SOP covers literature-based primer selection, in silico evaluation, and bench validation on the Bio-Rad QX200 ddPCR system (Bio-Rad Laboratories, 2024; Hindson et al., 2011).

2 Scope

  • Applies to: Design and selection of primers and hydrolysis probes (TaqMan) or EvaGreen-based assays for quantifying nitrogen cycle functional genes from environmental DNA extracts (soil, sediment, water)
  • Does not apply to: Amplicon sequencing library preparation, metagenomics, or RNA-based transcript quantification (RT-ddPCR)
  • Location: Molecular biology bench, Pomona College Biogeochemistry Lab

3 Responsibilities

Role Responsibility
Lab PI Approve primer selections and assay designs
Trained Personnel Conduct in silico analysis, order primers, perform validation runs
New Students Complete training before independent use

4 Definitions

  • Amplicon length: The size in base pairs (bp) of the PCR product generated by a primer pair; for ddPCR, optimal amplicon length is 60–200 bp; shorter amplicons partition more efficiently into droplets and produce cleaner positive/negative separation (Bio-Rad Laboratories, 2024; Hindson et al., 2011).
  • amoA (ammonia monooxygenase alpha subunit): The gene encoding the alpha subunit of ammonia monooxygenase; used as a functional marker for ammonia-oxidizing bacteria (AOB; Nitrosomonas-like) and ammonia-oxidizing archaea (AOA; Thaumarchaeota); quantified separately because AOB and AOA primers differ (Francis et al., 2005; Rotthauwe et al., 1997).
  • Annealing temperature (\(T_a\)): The temperature at which primers bind to their complementary template sequences during PCR cycling; determined by gradient PCR optimization; typically 55–65 °C for nitrogen cycle gene assays; primer pairs should be within 2 °C of each other (Bio-Rad Laboratories, 2024).
  • ddPCR (droplet digital PCR): A third-generation PCR technology that partitions a reaction into approximately 20,000 nanoliter droplets and amplifies target DNA to endpoint within each droplet; the fraction of positive droplets is fit to a Poisson distribution to calculate absolute target concentration without a standard curve (Hindson et al., 2011).
  • DNRA (dissimilatory nitrate reduction to ammonium): A microbial process that reduces NO3- to NH4+, retaining nitrogen in the ecosystem rather than losing it as gas; marker gene: nrfA (Levy-Booth et al., 2014).
  • EvaGreen: A fluorescent intercalating dye used as a non-specific reporter in ddPCR and qPCR; binds double-stranded DNA; less specific than hydrolysis probes but useful for assay development and for targets without validated probe sequences (Bio-Rad Laboratories, 2024).
  • FunGene: A curated functional gene database (http://fungene.cme.msu.edu/) that provides reference sequences, primer coverage analysis, and taxonomy for nitrogen cycle genes including amoA, nirK, nirS, nosZ, and nifH (Fish et al., 2013).
  • gBlock (synthetic DNA fragment): A double-stranded DNA fragment (IDT) with a sequence matching the target amplicon; used as a positive control for ddPCR assay validation and as a template for limit-of-detection determination (Hindson et al., 2011).
  • Hydrolysis probe (TaqMan probe): A short oligonucleotide labelled with a 5’ fluorophore (e.g., FAM) and a 3’ quencher (e.g., Iowa Black FQ); cleaved by Taq polymerase during extension, releasing the fluorophore and generating a signal; provides higher specificity than intercalating dyes and is required for multiplexed ddPCR (Bio-Rad Laboratories, 2024).
  • In silico primer evaluation: Computational assessment of primer specificity, coverage, and secondary structure using tools such as NCBI Primer-BLAST, IDT OligoAnalyzer, and FunGene; a required step before ordering primers for bench validation (Fish et al., 2013; Ye et al., 2012).
  • nifH (nitrogenase reductase): The gene encoding nitrogenase iron protein; the universal marker for biological nitrogen fixation (diazotrophy); highly conserved across diverse diazotrophic bacteria and archaea (Levy-Booth et al., 2014; Poly et al., 2001).
  • nirK / nirS: Genes encoding the two structurally distinct nitrite reductases in denitrifiers: nirK encodes a copper-containing (Cu) enzyme; nirS encodes a cytochrome cd1 enzyme; the two forms serve the same function and are often both targeted to capture the full denitrifying community (Braker et al., 2000; Henry et al., 2006).
  • nosZ: The gene encoding nitrous oxide reductase, which reduces N2O to N2 in the final step of denitrification; divided into Clade I (classical denitrifiers; Henry et al. (2006)) and Clade II (diverse non-denitrifiers; Jones et al. (2013)); the sum of both clades estimates the community capacity to consume the greenhouse gas N2O.
  • nrfA: The gene encoding cytochrome c nitrite reductase, the enzyme catalyzing DNRA; a functional marker for organisms that retain nitrogen in the ecosystem as NH4+ rather than releasing it as N2O or N2 (Welsh et al., 2014).
  • Rain (ddPCR): Droplets with intermediate fluorescence amplitude that fall between the positive and negative clusters in a ddPCR amplitude plot; indicates suboptimal primer binding, non-specific amplification, or partial inhibition; reduced by increasing annealing temperature or decreasing primer concentration (Lievens et al., 2016).

5 Safety and Hazards

WarningHazards
  • Chemical: ddPCR droplet generation oil contains proprietary surfactants; avoid skin and eye contact. EvaGreen and other intercalating dyes are potential mutagens at high concentrations; handle with gloves.
  • Physical: The QX200 Droplet Generator uses pressurized cartridges. Follow manufacturer instructions for cartridge handling.
  • Biological: Environmental DNA extracts may originate from samples containing pathogens. Handle with standard precautions.

5.1 Required PPE

5.2 Emergency Procedures

For chemical splash, flush with water for 15 minutes. Refer to SOP-L003 Hazardous Materials Handling.

6 Target Genes Overview

The following table summarizes the nitrogen cycle gene targets, the enzymatic steps they represent, and key references for established primer sets (Braker et al., 2000; Henry et al., 2006; Poly et al., 2001; R"osch et al., 2002; Rotthauwe et al., 1997; Welsh et al., 2014):

Gene Enzyme Process Organism Group
amoA Ammonia monooxygenase Ammonia oxidation (nitrification) AOB, AOA
nirK Copper-containing nitrite reductase Denitrification (NO2- to NO) Denitrifiers
nirS Cytochrome cd1 nitrite reductase Denitrification (NO2- to NO) Denitrifiers
nosZ Nitrous oxide reductase Denitrification (N2O to N2) Denitrifiers (Clade I and II)
nifH Nitrogenase reductase Nitrogen fixation (N2 to NH3) Diazotrophs
nrfA Cytochrome c nitrite reductase DNRA (NO2- to NH4+) DNRA organisms

7 Procedure

7.1 Literature-Based Primer Selection

  1. Search the published literature for validated primer sets targeting each gene of interest. Prioritized references include Rotthauwe et al. (1997), Francis et al. (2005), Braker et al. (2000), Henry et al. (2006), Poly et al. (2001), Jones et al. (2013), and Welsh et al. (2014).
  2. For each candidate primer pair, record: primer names, sequences (5’ to 3’), amplicon length, annealing temperature, and the original publication.
  3. Evaluate published primer sets against the following criteria:
    • Amplicon length between 60 and 200 bp (Bio-Rad Laboratories, 2024; Hindson et al., 2011)
    • Demonstrated use in environmental samples similar to your matrix
    • Broad coverage of the target gene clade verified by in silico analysis (see §6.3)
    • Minimal off-target amplification in non-target organisms
Note

Many classic nitrogen cycle primers (e.g., nosZ primers from Henry et al. (2006)) were designed for conventional qPCR with amplicons of 250–400 bp. These may need to be redesigned or replaced with shorter alternatives optimized for ddPCR (Bio-Rad Laboratories, 2024).

7.3 In Silico Primer Evaluation

  1. Check primer specificity using NCBI Primer-BLAST against the nr/nt database, restricting the search to bacteria and/or archaea as appropriate (Ye et al., 2012).
  2. Verify primer coverage against curated functional gene databases:
  3. Check for primer dimers and secondary structures using IDT OligoAnalyzer (https://www.idtdna.com/calc/analyzer) or Primer3 (https://primer3.ut.ee/).
  4. Confirm that the melting temperature (\(T_m\)) of each primer falls between 55 and 65 °C, with no more than 2 °C difference between forward and reverse primers.
  5. For hydrolysis probe (TaqMan) assays, design the probe to have a \(T_m\) 5–10 °C higher than the primers, with no G at the 5’ end (quenches FAM fluorescence), and a length of 15–30 bases (Bio-Rad Laboratories, 2024).

7.4 Primer and Probe Ordering

  1. Order desalted primers at 25 nmol scale from IDT (Integrated DNA Technologies) or equivalent supplier.
  2. Order hydrolysis probes HPLC-purified with 5’-FAM (or HEX for multiplexing) and 3’-Iowa Black FQ quencher.
  3. For positive controls, order synthetic double-stranded DNA fragments (gBlocks) matching the target amplicon sequence. Include at least 50 bp of flanking sequence on each side of the primer binding sites.
  4. Resuspend lyophilized primers in nuclease-free water to a stock concentration of 100 µM. Store at –20 °C.
  5. Prepare 10 µM working stocks by diluting 1:10 in nuclease-free water. Store at –20 °C.

7.5 ddPCR Thermal Cycling Optimization

  1. Set up a temperature gradient experiment on the Bio-Rad C1000 Touch to determine the optimal annealing temperature for each primer set.
  2. Use a thermal gradient spanning \(T_m\) ± 5 °C (typically 55–65 °C).
  3. Use positive control gBlock DNA diluted to approximately 1,000 copies per µL as template.
  4. Recommended cycling protocol (Bio-Rad Laboratories, 2024):
    • Enzyme activation: 95 °C for 10 min
    • 40 cycles: 94 °C for 30 s, then annealing/extension at gradient temperature for 60 s
    • Signal stabilization: 98 °C for 10 min
    • Hold at 4 °C (indefinite)
    • Ramp rate: 2 °C/s for all steps
  5. Read droplets on the QX200 Droplet Reader.
  6. Select the annealing temperature that produces the clearest separation between positive and negative droplet clusters in the QuantaSoft amplitude plot.
Tip

If positive and negative droplet populations overlap (rain), try increasing the annealing temperature by 1–2 °C or reducing the primer concentration from 900 nM to 450 nM (Lievens et al., 2016).

7.6 Assay Validation

  1. Specificity: Run no-template controls (NTC) and extraction blanks from SOP-L700 to confirm the absence of false positives.
  2. Dynamic range: Prepare a serial dilution of gBlock positive control (10 to 100,000 copies per reaction) and verify linearity of the ddPCR response across the range (Huggett et al., 2013).
  3. Limit of detection (LOD): Determine the lowest copy number that produces a statistically significant positive call above NTC background; typically 1–10 copies per 20 µL reaction (Hindson et al., 2011).
  4. Matrix effects: Run environmental DNA extracts at multiple dilutions (undiluted, 1:5, 1:10) to check for PCR inhibition. If the 1:10 dilution value (back-calculated) disagrees by more than 20% from the 1:5 value, inhibition is likely present (Ra vcki et al., 2014).
  5. Document all validation results in the primer/probe registry (Appendix A) before approving the assay for routine use.

8 Primer and Probe Registry

Maintain a lab registry of all validated and candidate assays. The registry must include at minimum: gene target, primer names and sequences, probe sequence (if applicable), amplicon length, validated annealing temperature, LOD, matrix tested, validation date, and analyst initials. Store the registry as a version-controlled file in the lab GitHub repository.

9 References

Bio-Rad Laboratories. (2024). QX200 droplet digital PCR system instruction manual. https://www.bio-rad.com/en-us/product/qx200-droplet-digital-pcr-system
Braker, G., Zhou, J., Wu, L., Devol, A. H., & Tiedje, J. M. (2000). Nitrite reductase genes (
textitnirK and
textitnirS) as functional markers to investigate diversity of denitrifying bacteria in Pacific Northwest marine sediment communities. Applied and Environmental Microbiology, 66(5), 2096–2104. https://doi.org/10.1128/AEM.66.5.2096-2104.2000
Fish, J. A., Chai, B., Wang, Q., Sun, Y., Brown, C. T., Tiedje, J. M., & Cole, J. R. (2013). FunGene: The functional gene pipeline and repository. Frontiers in Microbiology, 4, 291. https://doi.org/10.3389/fmicb.2013.00291
Francis, C. A., Roberts, K. J., Beman, J. M., Santoro, A. E., & Oakley, B. B. (2005). Ubiquity and diversity of ammonia-oxidizing archaea in water columns and sediments of the ocean. Proceedings of the National Academy of Sciences, 102(41), 14683–14688. https://doi.org/10.1073/pnas.0506625102
Henry, S., Bru, D., Stres, B., Hallet, S., & Philippot, L. (2006). Quantitative detection of the
textitnosZ gene, encoding nitrous oxide reductase, and comparison of the abundances of 16S rRNA,
textitnarG,
textitnirK, and
textitnosZ genes in soils. Applied and Environmental Microbiology, 72(8), 5181–5189. https://doi.org/10.1128/AEM.00231-06
Hindson, B. J., Ness, K. D., Masquelier, D. A., Belgrader, P., Heredia, N. J., Makarewicz, A. J., Bright, I. J., Lucero, M. Y., Hiddessen, A. L., Legler, T. C., Kitano, T. K., Hodel, M. R., Petersen, J. F., Wyatt, P. W., Steenblock, E. R., Shah, P. H., Bousse, L. J., Troup, C. B., Mellen, J. C., … Colston, B. W. (2011). High-throughput droplet digital PCR system for absolute quantitation of DNA copy number. Analytical Chemistry, 83(22), 8604–8610. https://doi.org/10.1021/ac202028g
Huggett, J. F., Foy, C. A., Benes, V., Emslie, K., Garson, J. A., Haynes, R., Hellemans, J., Kubista, M., Mueller, R. D., Nolan, T., Pfaffl, M. W., Shipley, G. L., Vandesompele, J., Wittwer, C. T., & Bustin, S. A. (2013). Digital PCR as a novel technology and its potential implications for molecular diagnostics. Clinical Chemistry, 59(12), 1691–1702. https://doi.org/10.1373/clinchem.2013.214742
Jones, C. M., Graf, D. R. H., Bru, D., Philippot, L., & Hallin, S. (2013). A comprehensive survey of soil
textitnosZ genes and transcripts reveals that nosZ clade II are important drivers of N
textsubscript2O reduction. The ISME Journal, 7(11), 2106–2119. https://doi.org/10.1038/ismej.2013.100
Levy-Booth, D. J., Prescott, C. E., & Grayston, S. J. (2014). Microbial functional genes involved in nitrogen fixation, nitrification and denitrification in forest ecosystems. Soil Biology and Biochemistry, 75, 11–25. https://doi.org/10.1016/j.soilbio.2014.03.021
Lievens, A., Jacchia, S., Kagkli, D., Savini, C., & Querci, M. (2016). Enhanced analysis of real-time PCR data by using a variable efficiency model: FPK-PCR. Nucleic Acids Research, 44(10), e95. https://doi.org/10.1093/nar/gkw007
Poly, F., Monrozier, L. J., & Bally, R. (2001). Comparison of
textitnifH gene pools in soils and soil microenvironments with contrasting properties. Applied and Environmental Microbiology, 67(5), 2255–2262. https://doi.org/10.1128/AEM.67.5.2255-2262.2001
Ra vcki, N., Dreo, T., Gutierrez-Aguirre, I., Blejec, A., & Ravnikar, M. (2014). Reverse transcriptase droplet digital PCR shows high resilience to PCR inhibitors from plant, soil and water samples. Plant Methods, 10, 42. https://doi.org/10.1186/s13007-014-0042-6
R"osch, C., Mergel, A., & Bothe, H. (2002). Biodiversity of denitrifying and dinitrogen-fixing bacteria in an acid forest soil. Applied and Environmental Microbiology, 68(8), 3818–3829. https://doi.org/10.1128/AEM.68.8.3818-3829.2002
Rotthauwe, J.-H., Witzel, K.-P., & Liesack, W. (1997). The ammonia monooxygenase structural gene
textitamoA as a functional marker: Molecular fine-scale analysis of natural ammonia-oxidizing populations. Applied and Environmental Microbiology, 63(12), 4704–4712. https://doi.org/10.1128/aem.63.12.4704-4712.1997
Welsh, A., Chee-Sanford, J. C., Connor, L. M., L"offler, F. E., & Sanford, R. A. (2014). Refined
textitnrfA phylogeny improves PCR-based
textitnrfA gene detection. Applied and Environmental Microbiology, 80(7), 2110–2119. https://doi.org/10.1128/AEM.03443-13
Ye, J., Coulouris, G., Zaretskaya, I., Cutcutache, I., Rozen, S., & Madden, T. L. (2012). Primer-BLAST: A tool to design target-specific primers for polymerase chain reaction. BMC Bioinformatics, 13, 134. https://doi.org/10.1186/1471-2105-13-134

11 Training Record

Before designing or ordering primers independently, personnel must:

  1. Read this SOP in full
  2. Complete the Bio-Rad ddPCR Foundations eLearning module
  3. Demonstrate competent use of Primer-BLAST and IDT OligoAnalyzer
  4. Perform a supervised gradient optimization experiment and select the annealing temperature with guidance from the PI
  5. Be signed off by the PI