Water quality data, gene copy number counts, and 16S analyses from Rhode River samples collected at SERC during 2017 and 2018 (NCEI Accession 0288317)

Acquisition Description:
Surface water samples were collected for subsequent microbial community analyses via PCR and shotgun metagenomics. A YSI EXO2 water quality sonde operated by the Smithsonian Environmental Research Center was used to collect physicochemical conditions at time of sample collection.

Shotgun metagenomic DNA and 16S rRNA microbial genes were amplified as described in Sakowski et al. (in revision). Metagenomic DNA was prepped using the Nextera Flex library prep kit. 16S rRNA libraries were prepped by V3-V4 PCR amplification. Both shotgun metagenome and 16S libraries were sequenced on an Illumina MiSeq.

qPCR of 16S rRNA genes and Viral Ribonucleotide Reductase Genes:
Bacterial and viral abundances were estimated by quantitative PCR. For Ribonucleotide Reductase (RNR) quantification, an environmental RNR amplicon was first cloned into chemically competent Escherichia coli cells using the Zero Blunt PCR Cloning Kit (Thermo Scientific) following the manufacturer's protocol. Environmental samples were quantitated for 16S rRNA and RNR gene copy numbers by comparing to a serial dilution of the cloned E. coli cells at known concentrations. All standards and environmental samples were run in triplicate. Three microliters of sample were combined with UltraPure molecular grade water (Thermo, Inc.), SsoAdvanced Universal SYBR Green Supermix (1x final concentration, Bio-Rad Laboratories, Inc.), Forward primer (0.3 mM final concentration), and Reverse primer (0.3 mM final concentration) to a final volume of 25 mL. Samples were amplified on a CFX96 Real-Time PCR Detection System (Bio-Rad Laboratories, Inc.) with the following conditions: denaturing at 98°C for 10 minutes; 45 cycles of denaturing at 98°C for 10 seconds, annealing at 52°C for 30 seconds, and extension at 72°C for 45 seconds; and a final extension of 72°C for 5 minutes.

16S rRNA amplification and sequencing:
16S rRNA genes were amplified from surface water samples in a 25 mL PCR reaction with the following conditions: three microliters of column-purified DNA were combined with UltraPure molecular grade water (Thermo, Inc.), 10X buffer (1x final concentration), dNTPs (0.1mM each final concentration), 16S forward primer 27F (0.3 mM final concentration), 16S reverse primer PE_16S_V4_E786_R (0.3 mM final concentration), bovine serum albumin (0.02 mg/mL final concentration), and Phusion High-Fidelity DNA Polymerase (0.5U; New England BioLabs, Inc. PCR reactions were combined with 150 mL of 4% UMIL EM90 oil (4% UMIL EM90 oil, 0.05% TritonX-100 v/v in mineral oil; Universal Preserv-A-Chem, Inc.) and emulsified by vortexing at max speed (~2,700 rpm) for one minute on a Vortex Genie 2 (MoBio). Emulsions were loaded as 50 mL aliquots and amplified with the following conditions: denaturation at 94°C for 3 minutes; 33 cycles of denaturation at 94°C for 10 seconds, annealing at 54°C for 30 seconds, and extension at 72°C for 45 seconds; and a final extension of 72°C for 5 minutes (C1000, BioRad Labs., Inc.).

PCR oil emulsions were broken with isobutanol. Briefly, PCR aliquots were pooled in a 1.5mL microcentrifuge tube and combined with 100 uL of sterile 5M NaCl solution and 1 mL of isobutanol. Samples were vortexed briefly to mix and centrifuged at 16,000 x g for 1 minute. The bottom aqueous layer was retained, and DNA was purified by spin column purification (Zymo, Inc.). DNA was eluted in 20 uL of Tris-HCl.

Purified DNA was run on a 1.5% agarose gel (UltraPure Agarose, ThermoFisher Scientific). 16S rRNA gene bands were visualized under blue light excitation, extracted, and gel purified (Zymo, Inc.) Purified DNA was eluted into 20 mL of Tris-HCl. Barcodes and Illumina adapters were added to 16S rRNA gene amplicon products in two subsequent limited PCR steps. Barcodes were added as follows: two microliters of purified DNA were combined with UltraPure molecular grade water (Thermo, Inc.), 10X buffer (1x final concentration), dNTPs (0.1mM each final concentration), 16S forward primer PE_16S_V4_U515F (0.3 mM final concentration), 16S rRNA gene reverse primer with 8-mer barcodes PE_IV_XXX (0.3 mM final concentration), and Phusion High-Fidelity DNA Polymerase (0.5U; New England BioLabs, Inc.). Samples were amplified with the following conditions: denaturing at 98°C for 30 seconds; 8 cycles of denaturing at 98°C for 10 seconds, annealing at 54°C for 30 seconds, and extension at 72°C for 45 seconds; and a final extension of 72°C for 5 minutes. DNA was purified by spin column purification (Zymo, Inc.) and eluted into 20 mL Tris-HCl. Illumina adapters were then added as above with the following primers: Illumina adapter forward primer PE-III-PCR-F (0.3 mM final concentration) and Illumina adapter reverse primer Barcode_Rev (0.3 mM final concentration). Samples were amplified with the following conditions: denaturing at 98°C for 30 seconds; 5 cycles of denaturing at 98°C for 10 seconds, annealing at 54°C for 30 seconds, and extension at 72°C for 45 seconds; and a final extension of 72°C for 5 minutes. DNA was purified by spin column purification (Zymo, Inc.) and eluted into 20 mL Tris-HCl. 16S rRNA gene amplicon libraries were sequenced on an Illumina MiSeq (2 x 300 bp) at the Genetic Core Research Facility at Johns Hopkins University.

Shotgun metagenomic library preparation and sequencing:
Surface water samples were filtered through a 0.2µm PES filter. DNA was extracted from filters for shotgun sequencing from water samples with the DNeasy PowerWater kit (Qiagen) following the manufacturer’s protocol with the following amendment: 20 µL of proteinase K was combined with 1 mL of solution PW1 in the bead tube. The bead tube was incubated at 65°C for ten minutes prior to bead beating. Libraries were prepared with the Nextera DNA Flex Library Prep kit (Illumina, Inc.) following the manufacture's protocol and sequenced on an Illumina MiSeq (2 x 300 bp) at the Genetic Core Research Facility at Johns Hopkins University.

epicPCR of environmental samples:
Surface water glycerol samples (25% v/v) were thawed on ice and one mL was added to three replicate 1.5 mL microcentrifuge tubes per sample. Samples were centrifuged at 25,000 x g for 10 minutes and resuspended after supernatant removal to reduce free viral particles. Thirty microliters of each sample was combined with UltraPure molecular grade water (Thermo, Inc.), 10X buffer (1x final concentration), dNTPs (0.1mM each final concentration), Viral Forward primer (1.0 mM final concentration), Viral Reverse_519R fusion primers (R1 and R2 combined, 0.01 mM each final concentration), 16S Reverseeverse primer (1.0 mM final concentration), bovine serum albumin (0.02 mg/mL final concentration), Tween-20 (0.8% v/v final concentration), and Phusion High-Fidelity DNA Polymerase (1.5U; New England BioLabs, Inc.) to a final volume of 75 mL. PCR reactions were combined with 450 mL of 4% UMIL EM90 oil (4% UMIL EM90 oil, 0.05% TritonX-100 v/v in mineral oil; Universal Preserv-A-Chem, Inc.) and emulsified by vortexing at max speed (~2,700 rpm) for one minute. Emulsions were loaded as 50 mL aliquots and amplified with the following conditions: denaturation at 94°C for 3 minutes; 33 cycles of denaturation at 94°C for 10 seconds, annealing at 54°C for 30 seconds, and extension at 72°C for 45 seconds; and a final extension of 72°C for 5 minutes (C1000, BioRad Labs., Inc.).

PCR oil emulsions were broken with isobutanol. Briefly, PCR aliquots were pooled in a 1.5mL microcentrifuge tube and combined with 100 uL of sterile 5M NaCl solution and 1 mL of isobutanol. Samples were vortexed briefly to mix and centrifuged at 16,000 x g for 1 minute. The bottom aqueous layer was retained, and DNA was purified by spin column purification (Zymo, Inc.). DNA was eluted in 20 uL of Tris-HCl.

Amplicons were enriched by nested PCR and subsequently sequenced on a PacBio Sequel with Sequel v3 chemistry (University of Maryland Institute for Genome Sciences). Circular consensus sequences were obtained from raw reads with the following parameters: minimum signal-to-noise ratio (SNR): 3, minimum length: 500bp, minimum passes: 10, minimum read score: 0.75, minimum predicted accuracy: 0.90.