Demultiplexed paired-end sequences were merged and performed quality control im- plementation (mean sequence quality score ≥ 30) and features table construction (amplicon sequences variants, ASVs) via DADA2 [29] plugin in QIIME2 (version 2019.10) [30]. An even sequence depth of 9452 reads ...