Enhancing inference of differential gene expression in metatranscriptomes from human microbial communities

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 essential for RNA-sequencing analysis, current metatranscriptomic approaches have only been benchmarked on simulated data, resulting in a lack of standard practices for analysis of real datasets. Here, we use mock communities (defined mixtures of microbial cells with known properties) to quantitatively assess robustness and susceptibility of current approaches to various confounders including organisms' low relative abundance, differential abundance, low prevalence, global transcriptional output changes, and compositional effects. We show that no current method is robust to all confounders and method performance on simulated data does not generalize to real datasets. We then apply the same approaches to MTX datasets generated from gnotobiotic mice colonized with defined consortia of human bacterial strains and show that the method nominated by the mock community comparisons successfully inferred cross-feeding dynamics that were subsequently validated in vitro. Finally, using metagenome-assembled genomes from a human clinical study, we leverage genome-level sequencing depth and detection of genes to exclude low information samples on a per-organism basis to overcome confounding low prevalence and enhance differential expression inference. We conclude that MTX benchmarking on real, non-simulated datasets can and should guide choice of methods and their implementation, enabling inference and validation of microbial metabolic strategies and interactions in vivo.