Metabolite differential connectivity analysis has been successful in investigating potential molecular mechanisms underlying different conditions in biological systems. Correlation and Mutual Information (MI) are two of the most common measures to quantify association and for building metabolite—metabolite association networks and to calculate differential connectivity. In this study, we investigated the performance of correlation and MI to identify significantly differentially connected metabolites. These association measures were compared on (i) 23 publicly available metabolomic data sets and 7 data sets from other fields, (ii) simulated data with known correlation structures, and (iii) data generated using a dynamic metabolic model to simulate real-life observed metabolite concentration profiles. In all cases, we found more differentially connected metabolites when using correlation indices as a measure for association than MI. We also observed that different MI estimation algorithms resulted in difference in performance when applied to data generated using a dynamic model. We concluded that there is no significant benefit in using MI as a replacement for standard Pearson’s or Spearman’s correlation when the application is to quantify and detect differentially connected metabolites.

BACKGROUND. Necrotizing soft-tissue infections (NSTIs) are rapidly progressing infections frequently complicated by septic shock and associated with high mortality. Early diagnosis is critical for patient outcome, but challenging due to vague initial symptoms. Here, we identified predictive biomarkers for NSTI clinical phenotypes and outcomes using a prospective multicenter NSTI patient cohort. METHODS. Luminex multiplex assays were used to assess 36 soluble factors in plasma from NSTI patients with positive microbiological cultures (n = 251 and n = 60 in the discovery and validation cohorts, respectively). Control groups for comparative analyses included surgical controls (n = 20), non-NSTI controls (i.e., suspected NSTI with no necrosis detected upon exploratory surgery, n = 20), and sepsis patients (n = 24). RESULTS. Thrombomodulin was identified as a unique biomarker for detection of NSTI (AUC, 0.95). A distinct profile discriminating mono- (type II) versus polymicrobial (type I) NSTI types was identified based on differential expression of IL-2, IL-10, IL-22, CXCL10, Fas-ligand, and MMP9 (AUC >0.7). While each NSTI type displayed a distinct array of biomarkers predicting septic shock, granulocyte CSF (G-CSF), S100A8, and IL-6 were shared by both types (AUC >0.78). Finally, differential connectivity analysis revealed distinctive networks associated with specific clinical phenotypes. CONCLUSIONS. This study identifies predictive biomarkers for NSTI clinical phenotypes of potential value for diagnostic, prognostic, and therapeutic approaches in NSTIs. TRIAL REGISTRATION. ClinicalTrials.gov NCT01790698. FUNDING. Center for Innovative Medicine (CIMED); Region Stockholm; Swedish Research Council; European Union; Vinnova; Innovation Fund Denmark; Research Council of Norway; Netherlands Organisation for Health Research and Development; DLR Federal Ministry of Education and Research; and Swedish Children’s Cancer Foundation.

Networks and network analyses are fundamental tools of systems biology. Networks are built by inferring pair-wise relationships among biological entities from a large number of samples such that subject-specific information is lost. The possibility of constructing these sample (individual)-specific networks from single molecular profiles might offer new insights in systems and personalized medicine and as a consequence is attracting more and more research interest. In this study, we evaluated and compared LIONESS (Linear Interpolation to Obtain Network Estimates for Single Samples) and ssPCC (single sample network based on Pearson correlation) in the metabolomics context of metabolite–metabolite association networks. We illustrated and explored the characteristics of these two methods on (i) simulated data, (ii) data generated from a dynamic metabolic model to simulate real-life observed metabolite concentration profiles, and (iii) 22 metabolomic data sets and (iv) we applied single sample network inference to a study case pertaining to the investigation of necrotizing soft tissue infections to show how these methods can be applied in metabolomics. We also proposed some adaptations of the methods that can be used for data exploration. Overall, despite some limitations, we found single sample networks to be a promising tool for the analysis of metabolomics data.

 Background: Necrotising soft tissue infections (NSTIs) are rapidly progressing bacterial infections usually caused by either several pathogens in unison (polymicrobial infections) orStreptococcus pyogenes(mono-microbial infection). These infections are rare and are associated with high mortality rates. However, the underlying pathogenic mechanisms in this heterogeneous group remain elusive.Methods: In this study, we built interactomes at both the population and individual levels consisting of host-pathogen interactions inferred from dual RNA-Seq gene transcriptomic profiles of the biopsies from NSTI patients.Results: NSTI type-specific responses in the host were uncovered. TheS. pyogenesmono-microbial subnetwork was enriched with host genes annotated with involved in cytokine production and regulation of response to stress. The polymicrobial network consisted of several significant associations between different species (S. pyogenes,Porphyromonas asaccharolyticaandEscherichia coli) and host genes. The host genes associated withS. pyogenesin this subnetwork were characterised by cellular response to cytokines. We further found several virulence factors including hyaluronan synthase, Sic1, Isp, SagF, SagG, ScfAB-operon, Fba and genes upstream and downstream of EndoS along with bacterial housekeeping genes interacting with the human stress and immune response in various subnetworks between host and pathogen.Conclusions: At the population level, we found aetiology-dependent responses showing the potential modes of entry and immune evasion strategies employed byS. pyogenes, congruent with general cellular processes such as differentiation and proliferation. After stratifying the patients based on the subject-specific networks to study the patient-specific response, we observed different patient groups with different collagens, cytoskeleton and actin monomers in association with virulence factors, immunogenic proteins and housekeeping genes which we utilised to postulate differing modes of entry and immune evasion for different bacteria in relationship to the patients’ phenotype.

During early lactation, dairy cows have a negative energy balance since their energy demands exceed their energy intake: in this study, we aimed to investigate the association between diet and plasma metabolomics profiles and how these relate to energy unbalance of course in the early-lactation stage. Holstein-Friesian cows were randomly assigned to a glucogenic (n = 15) or lipogenic (n = 15) diet in early lactation. Blood was collected in week 2 and week 4 after calving. Plasma metabolite profiles were detected using liquid chromatography–mass spectrometry (LC-MS), and a total of 39 metabolites were identified. Two plasma metabolomic profiles were available every week for each cow. Metabolite abundance and metabolite ratios were used for the analysis using the XGboost algorithm to discriminate between diet treatment and lactation week. Using metabolite ratios resulted in better discrimination performance compared with the metabolite abundances in assigning cows to a lipogenic diet or a glucogenic diet. The quality of the discrimination of performance of lipogenic diet and glucogenic diet effects improved from 0.606 to 0.753 and from 0.696 to 0.842 in week 2 and week 4 (as measured by area under the curve, AUC), when the metabolite abundance ratios were used instead of abundances. The top discriminating ratios for diet were the ratio of arginine to tyrosine and the ratio of aspartic acid to valine in week 2 and week 4, respectively. For cows fed the lipogenic diet, choline and the ratio of creatinine to tryptophan were top features to discriminate cows in week 2 vs. week 4. For cows fed the glucogenic diet, methionine and the ratio of 4-hydroxyproline to choline were top features to discriminate dietary effects in week 2 or week 4. This study shows the added value of using metabolite abundance ratios to discriminate between lipogenic and glucogenic diet and lactation weeks in early-lactation cows when using metabolomics data. The application of this research will help to accurately regulate the nutrition of lactating dairy cows and promote sustainable agricultural development.