DNA was extracted from early enrichments at 37, 50, and 70°C using the MO Bio PowerSoil DNA isolation kit (Qiagen). 16S rRNA gene amplicon libraries were prepared according to the Illumina 16S metagenomic sequencing library preparation protocol (support.illumina.com/documents/documentation/chemistry_documentation/16s/16s-metagenomic-library-prep-guide-15044223-b.pdf). We amplified the V3–V4 region for bacteria and the V4–V6 region for archaea (Supplementary Table S1). The 16S rRNA gene libraries were sequenced at CeBiTec (Bielefeld, Germany) on a MiSeq (Illumina; 2 × 300-bp paired-end run, 100 000 reads per library). Sequences were analyzed in R Statistical Software v 3.5.1 (R-project.org/) with DADA2 v. 1.14 [41 (link)]. DADA2 scripts used for 16S rRNA gene analysis are accessible on GitHub (github.com/dbenitom/Metagenomics_scripts/blob/main/dada2_archaea.R and github.com/dbenitom/Metagenomics_scripts/blob/main/dada2_bacteria.R).
Free full text: Click here