“…This assumption, which is made by most metabarcoding studies (Bik et al, 2012;Blaxter et al, 2005;Taberlet, Bonin, Zinger, & Coissac, 2018), allowed us to infer relative occurrences of each protist taxonomic unit among all samples of each datasets and allowed to infer the co-occurrence and co-exclusion networks. However, all clustering programs used to construct OTUs make assumptions about the best ways to handle the environmental sequencing data (Callahan et al, 2016;Caron & Hu, 2018;Mahé et al, 2015;Nebel, Pfabel, Stock, Dunthorn, & Stoeck, 2011;Rognes, Flouri, Nichols, Quince, & Mahé, 2016;Zhang, Kapli, Pavlidis, & Stamatakis, 2013) and these assumptions, along with the choice of molecular markers, may or may not lead to under-or over-estimations of species diversity. Here the reads were clustered into OTU with the program Swarm Mahé, Rognes, Quince, Vargas, & Dunthorn, 2014), which uses local clustering thresholds and a breaking phase to construct the OTUs.…”