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

Abstract

Genetical genomics experiments have now been routinely conducted to measure both the genetic markers and gene expression data on the same subjects. The gene expression levels are often treated as quantitative traits and are subject to standard genetic analysis in order to identify the gene expression quantitative loci (eQTL). However, the genetic architecture for many gene expressions may be complex, and poorly estimated genetic architecture may compromise the inferences of the dependency structures of the genes at the transcriptional level. In this paper we introduce a sparse conditional Gaussian graphical model for studying the conditional independent relationships among a set of gene expressions adjusting for possible genetic effects where the gene expressions are modeled with seemingly unrelated regressions. We present an efficient coordinate descent algorithm to obtain the penalized estimation of both the regression coefficients and the sparse concentration matrix. The corresponding graph can be used to determine the conditional independence among a group of genes while adjusting for shared genetic effects. Simulation experiments and asymptotic convergence rates and sparsistency are used to justify our proposed methods. By sparsistency, we mean the property that all parameters that are zero are actually estimated as zero with probability tending to one. We apply our methods to the analysis of a yeast eQTL data set and demonstrate that the conditional Gaussian graphical model leads to a more interpretable gene network than a standard Gaussian graphical model based on gene expression data alone.

Published in at http://dx.doi.org/10.1214/11-AOAS494 the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org)

References in corpus (6)

• High-dimensional graphs and variable selection with the Lasso
• Sparsistency and rates of convergence in large covariance matrix estimation
• Network exploration via the adaptive LASSO and SCAD penalties
• A sparse conditional Gaussian graphical model for analysis of genetical genomics data
• Large-Sample Learning of Bayesian Networks is NP-Hard
• Causal graphical models in systems genetics: A unified framework for joint inference of causal network and genetic architecture for correlated phenotypes

Cited by in corpus (29)

• Generalized Network Psychometrics: Combining Network and Latent Variable Models
• A sparse conditional Gaussian graphical model for analysis of genetical genomics data
• Simultaneous Variable and Covariance Selection with the Multivariate Spike-and-Slab Lasso
• Penalized Maximum Likelihood Estimation of Multi-layered Gaussian Graphical Models
• Estimating causal structure using conditional DAG models
• Inference for Sparse Conditional Precision Matrices
• Sparse Ising Models with Covariates
• Joint estimation of sparse multivariate regression and conditional graphical models
• Joint Association Graph Screening and Decomposition for Large-scale Linear Dynamical Systems
• A Projection Based Conditional Dependence Measure with Applications to High-dimensional Undirected Graphical Models
• Graph quilting: graphical model selection from partially observed covariances
• Joint Mean-Covariance Estimation via the Horseshoe with an Application in Genomic Data Analysis
• The conditional censored graphical lasso estimator
• An explicit mean-covariance parameterization for multivariate response linear regression
• Covariate-Assisted Bayesian Graph Learning for Heterogeneous Data
• Structured Regularization for conditional Gaussian Graphical Models
• Neuroimaging Feature Extraction using a Neural Network Classifier for Imaging Genetics
• Partial Gaussian Graphical Model Estimation
• A covariance-enhanced approach to multi-tissue joint eQTL mapping with application to transcriptome-wide association studies
• Change point detection in dynamic Gaussian graphical models: the impact of COVID-19 pandemic on the US stock market
• Bayesian Variable Selection in Multivariate Nonlinear Regression with Graph Structures
• Bayesian Analysis of Multiway Tables in Association Studies: A Model Comparison Approach
• Interaction Pursuit Biconvex Optimization
• New insights for the multivariate square-root lasso
• SOFAR: large-scale association network learning
• High-dimensional Gaussian graphical model for network-linked data
• Gene regulatory network in single cells based on the Poisson log-normal model
• Combining Smoothing Spline with Conditional Gaussian Graphical Model for Density and Graph Estimation
• Objective Bayes Covariate-Adjusted Sparse Graphical Model Selection