A new framework for exploratory network mediator analysis in omics data.

Omics methods are widely used in basic biology and translational medicine research. More and more omics data are collected to explain the impact of certain risk factors on clinical outcomes. To explain the mechanism of the risk factors, a core question is how to find the genes/proteins/metabolites that mediate their effects on the clinical outcome. Mediation analysis is a modeling framework to study the relationship between risk factors and pathological outcomes, via mediator variables. However, high-dimensional omics data are far more challenging than traditional data: (1) From tens of thousands of genes, can we overcome the curse of dimensionality to reliably select a set of mediators? (2) How do we ensure that the selected mediators are functionally consistent? (3) Many biological mechanisms contain nonlinear effects. How do we include nonlinear effects in the high-dimensional mediation analysis? (4) How do we consider multiple risk factors at the same time? To meet these challenges, we propose a new exploratory mediation analysis framework, medNet, which focuses on finding mediators through predictive modeling. We propose new definitions for predictive exposure, predictive mediator, and predictive network mediator, using a statistical hypothesis testing framework to identify predictive exposures and mediators. Additionally, two heuristic search algorithms are proposed to identify network mediators, essentially subnetworks in the genome-scale biological network that mediate the effects of single or multiple exposures. We applied medNet on a breast cancer data set and a metabolomics data set combined with food intake questionnaire data. It identified functionally consistent network mediators for the exposures' impact on the outcome, facilitating data interpretation.

Bayesian functional analysis for untargeted metabolomics data with matching uncertainty and small sample sizes.

Untargeted metabolomics based on liquid chromatography-mass spectrometry technology is quickly gaining widespread application, given its ability to depict the global metabolic pattern in biological samples. However, the data are noisy and plagued by the lack of clear identity of data features measured from samples. Multiple potential matchings exist between data features and known metabolites, while the truth can only be one-to-one matches. Some existing methods attempt to reduce the matching uncertainty, but are far from being able to remove the uncertainty for most features. The existence of the uncertainty causes major difficulty in downstream functional analysis. To address these issues, we develop a novel approach for Bayesian Analysis of Untargeted Metabolomics data (BAUM) to integrate previously separate tasks into a single framework, including matching uncertainty inference, metabolite selection and functional analysis. By incorporating the knowledge graph between variables and using relatively simple assumptions, BAUM can analyze datasets with small sample sizes. By allowing different confidence levels of feature-metabolite matching, the method is applicable to datasets in which feature identities are partially known. Simulation studies demonstrate that, compared with other existing methods, BAUM achieves better accuracy in selecting important metabolites that tend to be functionally consistent and assigning confidence scores to feature-metabolite matches. We analyze a COVID-19 metabolomics dataset and a mouse brain metabolomics dataset using BAUM. Even with a very small sample size of 16 mice per group, BAUM is robust and stable. It finds pathways that conform to existing knowledge, as well as novel pathways that are biologically plausible.

Severe disease during both primary and secondary dengue virus infections in pediatric populations.

Dengue is a global epidemic causing over 100 million cases annually. The clinical symptoms range from mild fever to severe hemorrhage and shock, including some fatalities. The current paradigm is that these severe dengue cases occur mostly during secondary infections due to antibody-dependent enhancement after infection with a different dengue virus serotype. India has the highest dengue burden worldwide, but little is known about disease severity and its association with primary and secondary dengue infections. To address this issue, we examined 619 children with febrile dengue-confirmed infection from three hospitals in different regions of India. We classified primary and secondary infections based on IgM:IgG ratios using a dengue-specific enzyme-linked immunosorbent assay according to the World Health Organization guidelines. We found that primary dengue infections accounted for more than half of total clinical cases (344 of 619), severe dengue cases (112 of 202) and fatalities (5 of 7). Consistent with the classification based on binding antibody data, dengue neutralizing antibody titers were also significantly lower in primary infections compared to secondary infections (P ≤ 0.0001). Our findings question the currently widely held belief that severe dengue is associated predominantly with secondary infections and emphasizes the importance of developing vaccines or treatments to protect dengue-naive populations.