Materials and Reagents

Pipette tips (20 μl, 200 μl, 1,000 μl) (Thermo Fisher Scientific, Thermo Scientific^TM, catalog numbers: 3521-HR , 3551-HR , 3101-05-HR )
1.5 ml Eppendorf tubes (Eppendorf, catalog number: T9661-1000EA )
Strip PCR tubes and caps (Roche Diagnostics, catalog number: 1667009001 )
PowerSoil^TM DNA Isolation Kit (Mo Bio Laboratories, catalog number: 12888-100 )
Bovine serum albumin (BSA) (Sigma-Aldrich, catalog number: A7906 )
Premix Taq (Takara Bio, catalog number: R004A )
dNTP mixture
Tris-HCl (pH 8.3)
Potassium chloride (KCl)
Magnesium chloride (MgCl₂)
AxyPrep^TM PCR Cleanup Kit (Corning, Axygen^®, catalog number: AP-PCR-250 )
Primers

Equipment

Pipettes (0-10 μl, 10-100 μl, 100-1,000 μl) (Eppendorf, catalog numbers: k03030 , k03031 , k03032 )
Vortex-genie^TM 2 (Mo Bio Laboratories, model: Vortex-Genie^® 2 )
PICO 17 centrifuge (Thermo Fisher Scientific, Thermo Scientific^TM, model: Heraeus^TM Pico^TM 17 )
Automated thermal cycler TP600 (Takara Bio, model: TP600 )

Software

Mothur (www.mothur.org)
BLAST+ (ftp://ftp.ncbi.nlm.nih.gov/blast/executables/blast+/LATEST/)
Megan5 (http://ab.inf.uni-tuebingen.de/software/megan5/)

Procedure

Sample collection and nucleic acid extraction
Extract genomic DNA from soil or sediment samples (0.25 g wet weight) using a PowerSoil^TM DNA Isolation Kit (Mo Bio Laboratories, USA) according to the manufacturer’s instructions.
THDP-PCR method (Figure 1)

Figure 1. Principle of THDP-PCR
1. First step of THDP-PCR
2. PCR clean-up
  Purify the entire volume (50 μl) of the first-step PCR product using an AxyPrep^TM PCR Cleanup Kit (Axygen, China) according to the manufacturer’s instructions to remove unbound degenerate primers.
  Note: Unbound long tagged 3’-end unmatched degenerate primers will also bind to templates but not trigger the polymerization and thereby prevent the binding of the single tag primer to the first-step PCR products, so PCR clean-up between the two steps is necessary.
3. Second step of THDP-PCR
  1. PCR reaction recipe. Final volume: 50 μl
    H₂O                                                  17-19 μl
    Tag-barcode (10 μM)                       4 μl
    PCR clean-up product                     2-4 μl
    Note: Including the negative-control sample.
    Premix Taq (see Table 1)                 25 μl
  2. Automated thermal cycler program
    1) Denaturation at 94 °C for 5 min
    2) 35-40 cycles at:
        94 °C for 1 min
        59 °C for 1 min
        72 °C for 1 min
    3) 72 °C for 10 min
  Note: The used clean-up product volume and cycle numbers in the second step PCR can be adjusted to suit different kinds of samples. For some low-abundance samples, try to increase the clean-up product volume to 4 μl and cycle numbers to 40.
Pyrosequencing
Mix together THDP-PCR products from different samples and sequence via Roche_454 pyrosequencing (GS FLX Titanium System) according to the manufacturer’s instructions.

Data analysis

Trim pyrosequencing reads using the command trim.seqs in mothur with length > 300 and q > 25, and no ambiguous sequences allowed.
Translate the trimmed sequences into amino acid sequences using FrameBot tools in Functional Gene Pipeline/Repository (http://fungene.cme.msu.edu/) to remove potential non-specific sequences which cannot be translated into CuMMO-related amino acid sequences.
Note: The nucleotide sequences which can pass the translation are retrieved using the command get.seqs in mothur from trimmed sequences, based on this step’s results. The following analysis steps are based on the retrieved nucleotide sequences, not polypeptides.
Identify and remove chimeras using the commands chimera.uchime and remove.seqs in mothur.
Split the sequences into different samples by the barcode sequence using the command split.seqs in mothur.
Classify the split sequences with BLAST and the lowest common ancestor method in MEGAN using the expanded CuMMO database.
Note: The original database and parameters used in this step are interpreted by Dumont et al. (2014) in the article and supplementary files. The reference database used in the classification is expanded by Wang et al. (2017) from the original database to suit the broad coverage of the highly degenerate primers, and the detailed information is listed in the article and supplementary files (Wang et al., 2017).

Acknowledgments

This protocol was adapted from our previous studies (Wang et al., 2017). This work was supported by the National Natural Science Foundation of China (NSFC) [grant numbers 31470222, 31170114].

References

Dumont, M. G., Luke, C., Deng, Y. C. and Frenzel, P. (2014). Classification of pmoA amplicon pyrosequences using BLAST and the lowest common ancestor method in MEGAN. Front Microbiol 5: 34.
Ledeker, B. M. and De Long, S. K. (2013). The effect of multiple primer-template mismatches on quantitative PCR accuracy and development of a multi-primer set assay for accurate quantification of pcrA gene sequence variants. J Microbiol Methods 94(3): 224-231.
Tuomivirta, T. T., Yrjala, K. and Fritze, H. (2009). Quantitative PCR of pmoA using a novel reverse primer correlates with potential methane oxidation in Finnish fen. Res Microbiol 160(10): 751-756.
Wang, J. G., Xia, F., Zeleke, J., Zou, B., Rhee, S. K. and Quan, Z. X. (2017). An improved protocol with a highly degenerate primer targeting copper-containing membrane-bound monooxygenase genes for community analysis of methane- and ammonia-oxidizing bacteria. FEMS Microbiol Ecol 93(3): fiw244.