Illumina bcl2fastq© program was used to demultiplex sequencing data. Fastp and FastQC v0.11.8 (http://www.bioinformatics.babraham.ac.uk, accessed on 4 November 2021) tools allowed checking for quality, adapter trimmed and filtered forward and reversing raw reads.
QIIME software V1.9.1 MiSeq was used to analyze sequencing data [25 (link)], including forward and reverse reads joining, chimera removal, data filtering and taxonomic annotation. To remove chimeric sequences from the reads, the Usearch 6.1 algorithm was used [26 (link)]. Moreover, based on a 97% identity threshold value, reads were clustered into operational taxonomic units (OTUs). PyNAST was used for the alignment of the sequences with reference to the Greengenes core reference database (version 13_8) [27 (link)]. For taxonomic assignment, the UCLUST classifier was used [28 (link)]. The data were expressed as relative abundance.
Free full text: Click here