Empirical Bayes Estimation of Sample-Specific Correlation Networks for Cross-Sectional and Longitudinal Omics Data
Saha, E.; Chen, H.; Ghosal, R.
Show abstract
Sample-specific molecular networks, which capture how an individual's gene regulatory structure may deviate from the population as a whole, are important tools for understanding heterogeneity in health and disease. However, most existing methods either do not accommodate uncertainty quantification or cannot handle missing data, non-Gaussian modalities, and longitudinal structure common in modern multi-omic studies. For example, copy number alterations in cancer genomic studies are non-Gaussian and often only partially observed, and their regulatory effects on gene expression may vary substantially across patients, variation that a single, population-averaged network cannot capture. To address these challenges, we propose PRISM (Personalized Relational Inference via Sample-specific Modeling), an empirical Bayes framework for estimating sample-specific covariance matrices from cross-sectional and longitudinal multi-omic data. PRISM models each individual sample's covariance matrix as drawn from an inverse-Wishart prior whose parameters are estimated from the data and allowed to vary across subpopulations or evolve over time, yielding a closed-form posterior and edge-level uncertainty quantification, ensuring scalability to high-dimensional omics data. In simulation experiments, PRISM achieves point-estimate accuracy comparable to or better than existing methods, substantially improved uncertainty calibration under non-Gaussian marginals, and lower computational cost, while maintaining robust coverage under moderate missingness. We apply PRISM to a yeast cell-cycle time course dataset, recovering time-varying biologically meaningful regulatory relationships independently documented in literature, and to paired copy number and gene expression data from ovarian cancer patients, identifying copy-number-expression interactions associated with survival, including TFs HNF4G and PAX8 which have independent evidence of clinical relevance in this cancer.
Matching journals
The top 4 journals account for 50% of the predicted probability mass.
Similar papers in this journal
- Trajectory-guided dimensionality reduction for multi-sample single-cell RNA-seq data reveals biologically relevant sample-level heterogeneity 95%
- GRNFomer: Accurate Gene Regulatory Network Inference Using Graph Transformer 95%
- bayNorm: Bayesian gene expression recovery, imputation and normalisation for single cell RNA-sequencing data 95%
Similar papers in this journal
Similar papers in this journal
"Similar papers" are the closest papers from that journal in the model's embedding space. They show what the match is built on, but the ranking comes mostly from a classifier over the whole training set, not from these examples alone.