SSU raw reads were quality filtered and processed using LotuS2 V2.34 (zkurt et al., 2022) into amplicon sequence variants (ASVs) using DADA2. To ensure that only AM fungi ASVs were retained, an edited version of SILVA v138.2 (Quast et al., 2013) with Glomeromycota sequences removed was used to identify and remove non-AM fungal ASVs. Taxonomy of the retained ASVs was assigned using MaarjAM virtual taxa (VT) with adjusted nomenclature (pik et al., 2010), using 98% coverage and 97% similarity thresholds. To assign mycorrhizal status, AM fungi were identified from the two SSU datasets (processed with  MaarjAM) by sub-setting ASVs with taxonomic assignment to the classes Glomeromycetes, Archaeosporomycetes, and Paraglomeromycetes.