Abstract Many microbes associate with higher eukaryotes and impact their vitality. To engineer microbiomes for host benefit, we must understand the rules of community assembly and maintenance that, in large part, demand an understanding of the direct interactions among community members. Toward this end, we have developed a Poisson-multivariate normal hierarchical model to learn direct interactions from the count-based output of standard metagenomics sequencing experiments. Our model controls for confounding predictors at the Poisson layer and captures direct taxon–taxon interactions at the multivariate normal layer using an ℓ1 penalized precision matrix. We show in a synthetic experiment that our method handily outperforms state-of-the-art methods such as SparCC and the graphical lasso (glasso). In a real in planta perturbation experiment of a nine-member bacterial community, we show our model, but not SparCC or glasso, correctly resolves a direct interaction structure among three community members that ..
