Fig 2 . Taxonomic composition at family and phylum level (A ), and taxonomic (B ) and functional (C ) richness per dataset. In A , unclassified families and those out of the most five abundant families within their phylum were grouped under “Others”. Similarly, Phyla out of the top six most abundant were grouped under “Others”. In B and C , barplots represent the total richness per dataset while boxplots indicate richness per sample and dataset. Taxonomic richness comparison is made at family level. In C , Humann2 columns represent the metagenomics data. EC= Enzyme commission numbers, and KO = KEGG ortholog numbers.