SRX5576854 station6_rep2 SRP189469 bp0 Total RNA was extracted from filters using a Qiagen RNeasy Mini Kit modifying the lysis step with the addition of Biospec zirconia/silica beads. Synthetic standards (Ambion ERCC RNA Spike-in mix) were then added based upon estimated total RNA (using extraction values from prior samples collected in this region, Alexander et al. 2015, PNAS). The rest of the Qiagen RNeasy Mini Kit protocol was then followed according to the manufacturer's instructions, incorporating the on-column DNase digestion step, using a Qiagen RNase-free DNase kit. Resulting total RNA was eluted with RNase-free water. Extracted samples were sequenced using the polyA pull-down step (Illumina Truseq library preparation kit) to enrich for eukaryotic mRNA. Libraries were sequenced with an Illumina NovaSeq 6000 at the Columbia Genome Center to produce 80 million 100-bp, paired-end reads. SRS4538864 hotlava_stn6_r2 station6_rep2 RNA-Seq METATRANSCRIPTOMIC other Illumina NovaSeq 6000 SRX5576855 station6_rep3 SRP189469 bp0 Total RNA was extracted from filters using a Qiagen RNeasy Mini Kit modifying the lysis step with the addition of Biospec zirconia/silica beads. Synthetic standards (Ambion ERCC RNA Spike-in mix) were then added based upon estimated total RNA (using extraction values from prior samples collected in this region, Alexander et al. 2015, PNAS). The rest of the Qiagen RNeasy Mini Kit protocol was then followed according to the manufacturer's instructions, incorporating the on-column DNase digestion step, using a Qiagen RNase-free DNase kit. Resulting total RNA was eluted with RNase-free water. Extracted samples were sequenced using the polyA pull-down step (Illumina Truseq library preparation kit) to enrich for eukaryotic mRNA. Libraries were sequenced with an Illumina NovaSeq 6000 at the Columbia Genome Center to produce 80 million 100-bp, paired-end reads. SRS4538865 hotlava_stn6_r3 station6_rep3 RNA-Seq METATRANSCRIPTOMIC other Illumina NovaSeq 6000 SRX5576856 station2_rep1 SRP189469 bp0 Total RNA was extracted from filters using a Qiagen RNeasy Mini Kit modifying the lysis step with the addition of Biospec zirconia/silica beads. Synthetic standards (Ambion ERCC RNA Spike-in mix) were then added based upon estimated total RNA (using extraction values from prior samples collected in this region, Alexander et al. 2015, PNAS). The rest of the Qiagen RNeasy Mini Kit protocol was then followed according to the manufacturer's instructions, incorporating the on-column DNase digestion step, using a Qiagen RNase-free DNase kit. Resulting total RNA was eluted with RNase-free water. Extracted samples were sequenced using the polyA pull-down step (Illumina Truseq library preparation kit) to enrich for eukaryotic mRNA. Libraries were sequenced with an Illumina NovaSeq 6000 at the Columbia Genome Center to produce 80 million 100-bp, paired-end reads. SRS4538866 hotlava_stn2_r1 station2_rep1 RNA-Seq METATRANSCRIPTOMIC other Illumina NovaSeq 6000 SRX5576857 station2_rep2 SRP189469 bp0 Total RNA was extracted from filters using a Qiagen RNeasy Mini Kit modifying the lysis step with the addition of Biospec zirconia/silica beads. Synthetic standards (Ambion ERCC RNA Spike-in mix) were then added based upon estimated total RNA (using extraction values from prior samples collected in this region, Alexander et al. 2015, PNAS). The rest of the Qiagen RNeasy Mini Kit protocol was then followed according to the manufacturer's instructions, incorporating the on-column DNase digestion step, using a Qiagen RNase-free DNase kit. Resulting total RNA was eluted with RNase-free water. Extracted samples were sequenced using the polyA pull-down step (Illumina Truseq library preparation kit) to enrich for eukaryotic mRNA. Libraries were sequenced with an Illumina NovaSeq 6000 at the Columbia Genome Center to produce 80 million 100-bp, paired-end reads. SRS4538867 hotlava_stn2_r2 station2_rep2 RNA-Seq METATRANSCRIPTOMIC other Illumina NovaSeq 6000 SRX5576858 station2_rep3 SRP189469 bp0 Total RNA was extracted from filters using a Qiagen RNeasy Mini Kit modifying the lysis step with the addition of Biospec zirconia/silica beads. Synthetic standards (Ambion ERCC RNA Spike-in mix) were then added based upon estimated total RNA (using extraction values from prior samples collected in this region, Alexander et al. 2015, PNAS). The rest of the Qiagen RNeasy Mini Kit protocol was then followed according to the manufacturer's instructions, incorporating the on-column DNase digestion step, using a Qiagen RNase-free DNase kit. Resulting total RNA was eluted with RNase-free water. Extracted samples were sequenced using the polyA pull-down step (Illumina Truseq library preparation kit) to enrich for eukaryotic mRNA. Libraries were sequenced with an Illumina NovaSeq 6000 at the Columbia Genome Center to produce 80 million 100-bp, paired-end reads. SRS4538868 hotlava_stn2_r3 station2_rep3 RNA-Seq METATRANSCRIPTOMIC other Illumina NovaSeq 6000 SRX5576859 station6_rep1 SRP189469 bp0 Total RNA was extracted from filters using a Qiagen RNeasy Mini Kit modifying the lysis step with the addition of Biospec zirconia/silica beads. Synthetic standards (Ambion ERCC RNA Spike-in mix) were then added based upon estimated total RNA (using extraction values from prior samples collected in this region, Alexander et al. 2015, PNAS). The rest of the Qiagen RNeasy Mini Kit protocol was then followed according to the manufacturer's instructions, incorporating the on-column DNase digestion step, using a Qiagen RNase-free DNase kit. Resulting total RNA was eluted with RNase-free water. Extracted samples were sequenced using the polyA pull-down step (Illumina Truseq library preparation kit) to enrich for eukaryotic mRNA. Libraries were sequenced with an Illumina NovaSeq 6000 at the Columbia Genome Center to produce 80 million 100-bp, paired-end reads. SRS4538869 hotlava_stn6_r1 station6_rep1 RNA-Seq METATRANSCRIPTOMIC other Illumina NovaSeq 6000