The abundance features in these datasets were profiled using the standard operating procedure for shotgun metagenomics implemented in the Qiita platform [61]. Specifically, the raw sequencing reads were processed using fastp and Minimap2 to remove low-quality sequences, adapter sequences and sequences...