, 2008; Lauber et al., 2009). Note that samples were placed in bead tubes containing solution C1 and incubated at 65 °C for 10 min, followed by 2 min of bead beating with the MoBio vortex

adapter; the remaining steps of the extraction were performed as directed by the manufacturer. PCR amplification of bacterial 16S rRNA genes using primers directed at variable regions V1 and V2 (positions 27–338 according to the Escherichia coli 16S rRNA gene numbering scheme) was achieved following the protocol described in our earlier publications (Fierer et al., 2008; Lauber et al., 2009). Briefly, amplicons generated from three PCR reactions per sample were pooled to reduce per-PCR variability and purified using the MoBio Ultra Clean PCR cleanup kit according to the manufacturer’s instructions and quantified (PicoGreen; Invitrogen, Carlsbad, CA). No-template

PCR controls were also performed. GDC0068 PCR products generated from each subsample contained a sub-sample-specific, error-correcting barcode, which allowed us to assemble a UK-371804 single composite sample for pyrosequencing by combining equal amounts of amplicon from each subsample. The composite sample was then gel purified (Qiaquick gel Clean up kit, Qiagen, Valencia, CA) and precipitated with ethanol to remove any remaining contaminants. DNA was sequenced using a Roche 454 FLX pyrosequencer. 16S rRNA gene sequences were processed according to the methods described in our previous publications (Fierer et al., 2008; Hamady et al., 2008). Briefly, sequences

<200 or >300 nt or with average quality scores of <25 were removed from the dataset, as were those with uncorrectable barcodes, ambiguous bases, or if the bacterial 16S rRNA gene-specific primer was absent. Acyl CoA dehydrogenase Sequences were then assigned to the specific subsamples based on their unique 12 nt barcode and then grouped into phylotypes at the 97% level of sequence identity using cd-hit (Li & Godzik, 2006) with a minimum coverage of 97%. We chose to group the phylotypes at 97% identity because this matches the limits of resolution of pyrosequencing (Kunin et al., 2010) and because the branch length so omitted contributes little to the tree and therefore to phylogenetic estimates of β diversity (Hamady et al., 2009). A representative for each phylotype was chosen by selecting the most abundant sequence in the phylotype, with ties being broken by choosing the longest sequence. A phylogenetic tree of the representative sequences was constructed using the Kimura 2-parameter model in Fast Tree (Price et al., 2009) after sequences were aligned with NAST (minimum 150 nt at 75% minimum identity) (DeSantis et al., 2006a) against the GreenGenes database (DeSantis et al., 2006b). Hypervariable regions were screened out of the alignment using PH Lane mask (http://greengenes.lbl.gov/).