Abstract:Biomarkers are often organized into networks, in which the strengths of network connections vary across subjects depending on subject‐specific covariates (e.g., genetic variants). Variation of network connections, as subject‐specific feature variables, has been found to predict disease clinical outcome. In this work, we develop a two‐stage method to estimate biomarker networks that account for heterogeneity among subjects and evaluate network's association with disease clinical outcome. In the first stage, we propose a conditional Gaussian graphical model with mean and precision matrix depending on covariates to obtain covariate‐dependent networks with connection strengths varying across subjects while assuming homogeneous network structure. In the second stage, we evaluate clinical utility of network measures (connection strengths) estimated from the first stage. The second stage analysis provides the relative predictive power of between‐region network measures on clinical impairment in the context of regional biomarkers and existing disease risk factors. We assess the performance of proposed method by extensive simulation studies and application to a Huntington's disease (HD) study to investigate the effect of HD causal gene on the rate of change in motor symptom through affecting brain subcortical and cortical grey matter atrophy connections. We show that cortical network connections and subcortical volumes, but not subcortical connections are identified to be predictive of clinical motor function deterioration. We validate these findings in an independent HD study. Lastly, highly similar patterns seen in the grey matter connections and a previous white matter connectivity study suggest a shared biological mechanism for HD and support the hypothesis that white matter loss is a direct result of neuronal loss as opposed to the loss of myelin or dysmyelination.This article is protected by copyright. All rights reserved

Conditional score matching for high-dimensional partial graphical models

Estimation of High-Dimensional Graphical Models Using Regularized Score Matching

Uncertainty quantification in high-dimensional linear models incorporating graphical structures with applications to gene set analysis

Partial Gromov-Wasserstein Learning for Partial Graph Matching

Estimation of Multiple Networks with Common Structures in Heterogeneous Subgroups.

Identifying disease‐associated biomarker network features through conditional graphical model

Generalized Score Matching

Neural-Kernelized Conditional Density Estimation

Joint estimation of multiple high-dimensional Gaussian copula graphical models

Sufficient dimension reduction for Gaussian copula graphical models

Generalized Precision Matrix for Scalable Estimation of Nonparametric Markov Networks

Estimating Differential Latent Variable Graphical Models with Applications to Brain Connectivity

Neural Score Matching for High-Dimensional Causal Inference

High-Dimensional Differential Parameter Inference in Exponential Family using Time Score Matching

Score-based causal feature selection for cancer risk prediction.

High-dimensional Gaussian Graphical Model for Network-Linked Data.

Integrating Approximate Single Factor Graphical Models.

A Two-Step Method for Estimating High-Dimensional Gaussian Graphical Models

Subgraph Similarity Maximal All-Matching over a Large Uncertain Graph.

A sparse conditional Gaussian graphical model for analysis of genetical genomics data

High-dimensional Mixed Graphical Models