“…Diamond was used to align all contigs to the GenBank nr and KEGG databases for taxonomic assignment and KEGG ID annotation, respectively (Buchfink et al, 2015; Kanehisa & Goto, 2000). We performed differential gene expression analyses as described above for Daphnia , after we mapped, using Kallisto (Bray et al, 2016), eRNA reads to the reference genomes and transcriptomes of species known a priori to persist in the communities ( A. falcatus, S. quadricauda , and P. subcapitata ) and highly abundant species identified by SqueezeMeta ( Brachionus plicatilis , Volvox cateri , Stylonychia lemnae , Chlamydomonas eustigma ) (Aeschlimann et al, 2014; Han et al, 2019; Hirooka et al, 2017; Nag Dasgupta et al, 2018; Prochnik et al, 2010; Schomaker & Dudycha, 2021; Suzuki et al, 2018). We used the SQMtools R package to further analyze the SqueezeMeta data and created stacked-bar plots to represent the most relatively abundant taxa (Puente-Sánchez et al, 2020).…”