ferrinetal_nmediatewarmingimpactonsoilhexapods

This study was conducted at the ForHot research site in Iceland (Sigurdsson et al., 2016) between August 2017 and June 2018 (64°0′N, 21°11′W). Soil type was a Brown Andosol (Arnalds, 2015). Mean annual temperature at the site was 5.1 °C. The coldest and warmest temperatures in the neighboring village of Eyrarbakki in 2016 were -12.3 °C and 21.6 °C, respectively. Average annual precipitation for the same year was 1153 mm (Icelandic Meteorological Office, 2016). The vegetation was an unmanaged grassland dominated by <em>Agrostis capillaris</em> L., <em>Galium boreale </em>L. and<em> Anthoxantum odoratum</em> L. Vascular plants cover 46% of the area over a moss mat which covers up to 88% of the ground. Natural N deposition in the area is 1.3 ± 0.1kg N ha-1 y-1 (Leblans et al., 2014). Five transects were established, each one consisting of two 2 x 2 m plots, and each plot with different treatment. Soil cores were collected using an auger to a depth of ~10 cm, excluding the O horizon. Soil cores were sampled seasonally four times: August 2017, corresponding to late growing season; November 2017, at start of winter and initial soil freezing; April 2018, with the first soil thaw in un-warmed soils, and June 2018, in the early part of the growing season. We thus collected a total of 20 core samples for each warming treatment. All samples were immediately sieved to remove roots and stones larger than 2 mm. Fifteen grams of each sample were then frozen in plastic bags in liquid N in the field to immediately stop all biological processes. All frozen samples were freeze-dried in the laboratory. eDNA was extracted from 15 g soil samples belonging to DNA remains (i.e. no alive fauna) as previously described (Taberlet et al., 2012; Zinger et al., 2016). The soil hexapod communities were genetically characterized based on Molecular Operational Taxonomic Units (MOTUs) using the retrieved eDNA and applying a metabarcoding approach. We amplified the 16S mitochondrial rDNA region using the Ins16S_l primer pair (Ins16S_1-F: 5′-TRRGACGAGAAGACCCTATA-3′; Ins16_1-R: 5′-TCTTAATCCAACATCGAGGTC-3′; Clarke et al. 2014). This primer pair, specifically designed for hexapod metabarcoding, introduces a very limited taxonomic bias and performs very well for identifications at the species level throughout the Hexapoda subphylum (e.g. Kocher et al., 2017; Talaga et al., 2017). PCR amplification was performed in triplicate in 20-μL mixtures consisting of 10 μL of AmpliTaq Gold Master Mix (Life Technologies, Carlsbad, USA), 5.84 μL of nuclease-free Ambion water (Thermo Fisher Scientific, Waltham, USA), 0.25 μM each primer, 3.2 μg of bovine serum albumin (Roche Diagnostic, Basel, Switzerland) and 2 μl of DNA template that was diluted 10-fold to reduce PCR inhibition by humic substances. The thermal profile of the PCR amplification was 40 cycles of denaturation at 95 °C (30 s), annealing at 49 °C (30 s) and elongation at 72 °C (60 s), with a final elongation step at 72 °C for 7 min. Tags had at least five differences between them to minimize ambiguities (Coissac et al., 2012). The sequenced multiplexes comprised extractions/PCR blank controls, unused tag combinations and positive controls (Kocher et al., 2017). The PCR products were then sequenced using the MiSeq platform (Illumina Inc., San Diego, USA), with the expected sequencing depth set at 400 000 reads per sample. <br> The sequences were processed using OBITOOLS software (Boyer et al., 2016). Low-quality sequences (containing Ns, alignment scores &lt;50, lengths &lt;140 bp or &gt;320 bp and singletons) were excluded. The remaining sequences were clustered into MOTUs using SUMACLUST (Mercier et al., 2013) at a threshold of sequence similarity of 97%. The hexapod MOTUs were taxonomically assigned using Blast. MOTUs showing &lt;80% similarity with either the local or the EMBL reference databases were removed, leading to 219 MOTUs. These retained MOTUs included taxa from classes Insecta and Entognatha, which both belong to the subphylum Hexapoda. We then applied a post-processing pipeline (Zinger et al., 2021) to minimize PCR and sequencing errors, contaminations and false-positive sequences, and by detailed curation of ecologically incongruent assignments (i.e. taxa with distributions outside the palearctic and neartic ecozones). This conservative approach retained a total of 33 identified species. We then used checklists of Icelandic hexapod species and information from previous studies at the same study site (Fjellberg, 2007; Holmstrup et al., 2018) to assess the performance of our eDNA metabarcoding protocol to properly describe the hexapod communities in the soil.