design: 16S rRNA gene amplicon sequencing targeting the V3V4 hypervariable region. Total DNA was extracted from soil samples using the DNeasy PowerSoil kit, PCR-amplified with Illumina overhang-adapted primers, indexed using Illumina UD indexes, and sequenced on an Illumina NextSeq platform using a 2  301 bp paired-end configuration.