design: The PCR products were sequenced on the Illumina MiSeq platform, using MiSeq Reagent Kit V3 600-cycles (Illumina, San Diego, USA), according to the manufacturers specifications.DNA sequence analysis was performed in the R Project environment, using the DADA2 (Callahan et al., 2016) and phyloseq (McMurdie and Holmes, 2013) packages. Briefly, multiplexed sequences were assigned to a sample metadata table. Then, FASTQ files were filtered to remove primers, adapters, barcodes, and chimeras. Phred score was settled to Q > 30 to remove low-quality sequences, resulting in 2 x 250 bp paired-end reads. Finally, sequences were assigned to Amplicon Sequencing Variants (ASVs) and classified according to the SILVA 132 database (Quast et al., 2012).