Systematic evaluation of metatranscriptomic differential gene expression in silico, in vitro, and in vivo enables elucidation of inter-species cross-feeding.
Metatranscriptomic (MTX) sequencing quantifies gene expression from the collective genomes of microbial communities (microbiomes), enabling assessment of functional activity rather than functional potential. While differential expression testing is instrumental to RNA-sequencing analysis, current metatranscriptomic approaches have been benchmarked only on simulated data and not under real operating conditions, resulting in a lack of standard practices. Here, we evaluate the performance of statistical differential expression methods on both simulated datasets and data collected from real bacterial 'mock communities' designed for this purpose. We assess the robustness of individual methods to organisms' low relative abundance, differential abundance, low prevalence, and transcription rate changes, showing that no existing methods perform adequately across all confounding conditions. We then apply the same approaches to metatranscriptomic datasets generated from gnotobiotic mice colonized with defined consortia of human bacterial strains and show that the method nominated by our mock community comparisons successfully inferred cross-feeding dynamics which were validated in vitro. We conclude that MTX method benchmarking on real, not simulated, datasets can and should optimize model implementation, enabling inference and validation of cross-feeding and other inter-species and host-microbe dynamics from in vivo studies.