Large-scale correlation network construction for unraveling the coordination of complex biological systems.

Advanced measurement and data storage technologies have enabled high-dimensional profiling of complex biological systems. For this, modern multiomics studies regularly produce datasets with hundreds of thousands of measurements per sample, enabling a new era of precision medicine. Correlation analysis is an important first step to gain deeper insights into the coordination and underlying processes of such complex systems. However, the construction of large correlation networks in modern high-dimensional datasets remains a major computational challenge owing to rapidly growing runtime and memory requirements. Here we address this challenge by introducing CorALS (Correlation Analysis of Large-scale (biological) Systems), an open-source framework for the construction and analysis of large-scale parametric as well as non-parametric correlation networks for high-dimensional biological data. It features off-the-shelf algorithms suitable for both personal and high-performance computers, enabling workflows and downstream analysis approaches. We illustrate the broad scope and potential of CorALS by exploring perspectives on complex biological processes in large-scale multiomics and single-cell studies.

Machine learning models of plasma proteomic data predict mood in chronic stroke and tie it to aberrant peripheral immune responses.

Post-stroke depression is common, long-lasting and associated with severe morbidity and death, but mechanisms are not well-understood. We used a broad proteomics panel and developed a machine learning algorithm to determine whether plasma protein data can predict mood in people with chronic stroke, and to identify proteins and pathways associated with mood. We used Olink to measure 1,196 plasma proteins in 85 participants aged 25 and older who were between 5 months and 9 years after ischemic stroke. Mood was assessed with the Stroke Impact Scale mood questionnaire (SIS3). Machine learning multivariable regression models were constructed to estimate SIS3 using proteomics data, age, and time since stroke. We also dichotomized participants into better mood (SIS3 > 63) or worse mood (SIS3 ≤ 63) and analyzed candidate proteins. Machine learning models verified that there is indeed a relationship between plasma proteomic data and mood in chronic stroke, with the most accurate prediction of mood occurring when we add age and time since stroke. At the individual protein level, no single protein or set of proteins predicts mood. But by using univariate analyses of the proteins most highly associated with mood we produced a model of chronic post-stroke depression. We utilized the fact that this list contained many proteins that are also implicated in major depression. Also, over 80% of immune proteins that correlate with mood were higher with worse mood, implicating a broadly overactive immune system in chronic post-stroke depression. Finally, we used a comprehensive literature review of major depression and acute post-stroke depression. We propose that in chronic post-stroke depression there is over-activation of the immune response that then triggers changes in serotonin activity and neuronal plasticity leading to depressed mood.

Julia for biologists.

Major computational challenges exist in relation to the collection, curation, processing and analysis of large genomic and imaging datasets, as well as the simulation of larger and more realistic models in systems biology. Here we discuss how a relative newcomer among programming languages-Julia-is poised to meet the current and emerging demands in the computational biosciences and beyond. Speed, flexibility, a thriving package ecosystem and readability are major factors that make high-performance computing and data analysis available to an unprecedented degree. We highlight how Julia's design is already enabling new ways of analyzing biological data and systems, and we provide a list of resources that can facilitate the transition into Julian computing.

Early prediction and longitudinal modeling of preeclampsia from multiomics.

Preeclampsia is a complex disease of pregnancy whose physiopathology remains unclear. We developed machine-learning models for early prediction of preeclampsia (first 16 weeks of pregnancy) and over gestation by analyzing six omics datasets from a longitudinal cohort of pregnant women. For early pregnancy, a prediction model using nine urine metabolites had the highest accuracy and was validated on an independent cohort (area under the receiver-operating characteristic curve [AUC] = 0.88, 95% confidence interval [CI] [0.76, 0.99] cross-validated; AUC = 0.83, 95% CI [0.62,1] validated). Univariate analysis demonstrated statistical significance of identified metabolites. An integrated multiomics model further improved accuracy (AUC = 0.94). Several biological pathways were identified including tryptophan, caffeine, and arachidonic acid metabolisms. Integration with immune cytometry data suggested novel associations between immune and proteomic dynamics. While further validation in a larger population is necessary, these encouraging results can serve as a basis for a simple, early diagnostic test for preeclampsia.

Integration of mechanistic immunological knowledge into a machine learning pipeline improves predictions.

The dense network of interconnected cellular signalling responses that are quantifiable in peripheral immune cells provides a wealth of actionable immunological insights. Although high-throughput single-cell profiling techniques, including polychromatic flow and mass cytometry, have matured to a point that enables detailed immune profiling of patients in numerous clinical settings, the limited cohort size and high dimensionality of data increase the possibility of false-positive discoveries and model overfitting. We introduce a generalizable machine learning platform, the immunological Elastic-Net (iEN), which incorporates immunological knowledge directly into the predictive models. Importantly, the algorithm maintains the exploratory nature of the high-dimensional dataset, allowing for the inclusion of immune features with strong predictive capabilities even if not consistent with prior knowledge. In three independent studies our method demonstrates improved predictions for clinically relevant outcomes from mass cytometry data generated from whole blood, as well as a large simulated dataset. The iEN is available under an open-source licence.

A longitudinal study of the post-stroke immune response and cognitive functioning: the StrokeCog study protocol.

BACKGROUND: Stroke increases the risk of cognitive impairment even several years after the stroke event. The exact mechanisms of post-stroke cognitive decline are unclear, but the immunological response to stroke might play a role. The aims of the StrokeCog study are to examine the associations between immunological responses and long-term post-stroke cognitive trajectories in individuals with ischemic stroke. METHODS: StrokeCog is a single-center, prospective, observational, cohort study. Starting 6-12 months after stroke, comprehensive neuropsychological assessment, plasma and serum, and psychosocial variables will be collected at up to 4 annual visits. Single cell sequencing of peripheral blood monocytes and plasma proteomics will be conducted. The primary outcome will be the change in global and domain-specific neuropsychological performance across annual evaluations. To explain the differences in cognitive change amongst participants, we will examine the relationships between comprehensive immunological measures and these cognitive trajectories. It is anticipated that 210 participants will be enrolled during the first 3 years of this 4-year study. Accounting for attrition, an anticipated final sample size of 158 participants with an average of 3 annual study visits will be available at the completion of the study. Power analyses indicate that this sample size will provide 90% power to detect an average cognitive change of at least 0.23 standard deviations in either direction. DISCUSSION: StrokeCog will provide novel insight into the relationships between immune events and cognitive change late after stroke.

VoPo leverages cellular heterogeneity for predictive modeling of single-cell data.

High-throughput single-cell analysis technologies produce an abundance of data that is critical for profiling the heterogeneity of cellular systems. We introduce VoPo (https://github.com/stanleyn/VoPo), a machine learning algorithm for predictive modeling and comprehensive visualization of the heterogeneity captured in large single-cell datasets. In three mass cytometry datasets, with the largest measuring hundreds of millions of cells over hundreds of samples, VoPo defines phenotypically and functionally homogeneous cell populations. VoPo further outperforms state-of-the-art machine learning algorithms in classification tasks, and identified immune-correlates of clinically-relevant parameters.

Multimodal network diffusion predicts future disease-gene-chemical associations.

MOTIVATION: Precision medicine is an emerging field with hopes to improve patient treatment and reduce morbidity and mortality. To these ends, computational approaches have predicted associations among genes, chemicals and diseases. Such efforts, however, were often limited to using just some available association types. This lowers prediction coverage and, since prior evidence shows that integrating heterogeneous data is likely beneficial, it may limit accuracy. Therefore, we systematically tested whether using more association types improves prediction. RESULTS: We study multimodal networks linking diseases, genes and chemicals (drugs) by applying three diffusion algorithms and varying information content. Ten-fold cross-validation shows that these networks are internally consistent, both within and across association types. Also, diffusion methods recovered missing edges, even if all the edges from an entire mode of association were removed. This suggests that information is transferable between these association types. As a realistic validation, time-stamped experiments simulated the predictions of future associations based solely on information known prior to a given date. The results show that many future published results are predictable from current associations. Moreover, in most cases, using more association types increases prediction coverage without significantly decreasing sensitivity and specificity. In case studies, literature-supported validation shows that these predictions mimic human-formulated hypotheses. Overall, this study suggests that diffusion over a more comprehensive multimodal network will generate more useful hypotheses of associations among diseases, genes and chemicals, which may guide the development of precision therapies. AVAILABILITY AND IMPLEMENTATION: Code and data are available at https://github.com/LichtargeLab/multimodal-network-diffusion. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.