β-Glucuronidase Pattern Predicted From Gut Metagenomes Indicates Potentially Diversified Pharmacomicrobiomics

β-glucuronidases (GUS) of intestinal bacteria remove glucuronic acid from glucoronides, reversing phase II metabolism of the liver and affecting the level of active deconjugated metabolites deriving from drugs or xenobiotics. Two hundred seventy-nine non-redundant GUS sequences are known in the gut microbiota, classified in seven structural categories (NL, L1, L2, mL1, mL2, mL1,2, and NC) with different biocatalytic properties. In the present study, the intestinal metagenome of 60 healthy subjects from five geographically different cohorts was assembled, binned, and mined to determine qualitative and quantitative differences in GUS profile, potentially affecting response to drugs and xenobiotics. Each metagenome harbored 4-70 different GUS, altogether accounting for 218. The amount of intestinal bacteria with at least one GUS gene was highly variable, from 0.7 to 82.2%, 25.7% on average. No significant difference among cohorts could be identified, except for the Ethiopia (ETH) cohort where GUS-encoding bacteria were significantly less abundant. The structural categories were differently distributed among the metagenomes, but without any statistical significance related to the cohorts. GUS profiles were generally dominated by the category NL, followed by mL1, L2, and L1. The GUS categories most involved in the hydrolysis of small molecules, including drugs, are L1 and mL1. Bacteria contributing to these categories belonged to , , , , , , , and . Bacteria harboring L1 GUS were generally scarcely abundant (