Inferring Correlation Networks from Genomic Survey Data

TLDR

High‑throughput sequencing such as 16S rRNA profiling reveals microbial community structure, but correlation analysis of the resulting compositional data can produce unreliable, artifactual dependencies. The authors show that community diversity modulates compositional effects and introduce SparCC, a method that estimates correlations from compositional data. Simulations and Human Microbiome Project data reveal that compositional effects generate widespread artifactual correlations, and that SparCC uncovers a dense ecological network while the standard approach produces three spurious interactions per true one and misses 60 % of true links, especially in low‑diversity samples.

Abstract

High-throughput sequencing based techniques, such as 16S rRNA gene profiling, have the potential to elucidate the complex inner workings of natural microbial communities - be they from the world's oceans or the human gut. A key step in exploring such data is the identification of dependencies between members of these communities, which is commonly achieved by correlation analysis. However, it has been known since the days of Karl Pearson that the analysis of the type of data generated by such techniques (referred to as compositional data) can produce unreliable results since the observed data take the form of relative fractions of genes or species, rather than their absolute abundances. Using simulated and real data from the Human Microbiome Project, we show that such compositional effects can be widespread and severe: in some real data sets many of the correlations among taxa can be artifactual, and true correlations may even appear with opposite sign. Additionally, we show that community diversity is the key factor that modulates the acuteness of such compositional effects, and develop a new approach, called SparCC (available at https://bitbucket.org/yonatanf/sparcc), which is capable of estimating correlation values from compositional data. To illustrate a potential application of SparCC, we infer a rich ecological network connecting hundreds of interacting species across 18 sites on the human body. Using the SparCC network as a reference, we estimated that the standard approach yields 3 spurious species-species interactions for each true interaction and misses 60% of the true interactions in the human microbiome data, and, as predicted, most of the erroneous links are found in the samples with the lowest diversity.