Abstract Multi-omics analyses are increasingly employed in microbiome studies to obtain a holistic view of molecular changes occurring within microbial communities exposed to different conditions. However, it is not always clear to what extent each omics data type contributes to our understanding of the community dynamics and whether they are concordant with each other. Here we map the molecular response of a synthetic community of 32 human gut bacteria to three non-antibiotic drugs by using five omics layers, namely 16S rRNA gene profiling, metagenomics, metatranscriptomics, metaproteomics, and metabolomics. Using this controlled setting, we find that all omics methods with species resolution in their readouts are highly consistent in estimating relative species abundances across conditions. Furthermore, different omics methods complement each other in their ability to capture functional changes in response to the drug perturbations. For example, while nearly all omics data types captured that the antipsychotic drug chlorpromazine selectively inhibits Bacteroidota representatives in the community, the metatranscriptome and metaproteome suggested that the drug induces stress responses related to protein quality control and metabolomics revealed a decrease in polysaccharide uptake, likely caused by Bacteroidota depletion. Taken together, our study provides insights into how multi-omics datasets can be utilised to reveal complex molecular responses to external perturbations in microbial communities.

Consistency across multi-omics layers in a drug-perturbed gut microbial community