Abstract:Background: Understanding genetic underpinnings of immune-mediated inflammatory diseases is crucial to improve treatments. Single-cell RNA sequencing (scRNA-seq) identifies cell states expanded in disease, but often overlooks genetic causality due to cost and small genotyping cohorts. Conversely, large genome-wide association studies (GWAS) are commonly accessible. Methods: We present a 3-step robust benchmarking analysis of integrating GWAS and scRNA-seq to identify genetically relevant cell states and genes in inflammatory diseases. First, we applied and compared the results of two recent algorithms based on networks (scGWAS) or single-cell disease scores (scDRS), according to accuracy/sensitivity and interpretability (M. J. Zhang et al. 2022; Jia et al. 2022). While previous studies focused on coarse cell types, we used disease-specific, fine-grained single-cell atlases (183,742 and 228,211 cells) and GWAS data (Ns of 97,1,73 and 45,975) for rheumatoid arthritis (RA) and ulcerative colitis (UC) (F. Zhang et al. 2023; Ishigaki et al. 2022; Smillie et al. 2019; de Lange et al. 2017). Second, given the lack of scRNA-seq for many diseases with GWAS, we further tested the tools' resolution limits by differentiating between similar diseases with only one fine-grained scRNA-seq atlas. Lastly, we provide a novel evaluation of noncoding SNP incorporation methods by testing which enabled the highest sensitivity/accuracy of known cell-state calls. Results: We first found that single-cell based tool scDRS called superior numbers of supported cell states, like MERTK+ myeloid cells in RA, which were overlooked by network-based scGWAS. While scGWAS was advantageous for gene exploration, scDRS captured cellular heterogeneity of disease-relevance without single-cell genotyping. For noncoding SNP integration, we found a key trade-off between statistical power and confidence with positional (e.g. MAGMA) and non-positional approaches (e.g. chromatin-interaction, eQTL). Even when directly incorporating noncoding SNPs through 5' scRNA-seq measures of regulatory elements, non disease-specific atlases gave misleading results by not containing disease-tissue specific transcriptomic patterns. Despite this criticality of tissue-specific scRNA-seq, we showed that scDRS enabled deconvolution of two similar diseases with a single fine-grained scRNA-seq atlas and separate GWAS. Indeed, we identified supported and novel genetic-phenotype linkages separating RA and ankylosing spondylitis, and UC and crohn's disease. Overall, while noting evolving single-cell technologies, our study provides key findings for integrating expanding fine-grained scRNA-seq, GWAS, and noncoding SNP resources to unravel the complexities of inflammatory diseases.

Integrative analysis of 3604 GWAS reveals multiple novel cell type-specific regulatory associations

FORGEdb: a tool for identifying candidate functional variants and uncovering target genes and mechanisms for complex diseases

Joint analysis of functional genomic data and genome-wide association studies of 18 human traits

DeepWAS: Multivariate Genotype-Phenotype Associations by Directly Integrating Regulatory Information Using Deep Learning

A New Approach for Discovering Functional Links Connecting Non-Coding Regulatory Variants to Gene Targets

Single-cell based integrative analysis of transcriptomics and genetics reveals robust associations and complexities for inflammatory diseases

Integrative modeling of eQTLs and cis-regulatory elements suggest mechanisms underlying cell type specificity of eQTLs

Multitissue H3K27ac profiling of GTEx samples links epigenomic variation to disease

GWAS4D: multidimensional analysis of context-specific regulatory variant for human complex diseases and traits

Genetic Analyses of Blood Cell Structure for Biological and Pharmacological Inference

Discovery of target genes and pathways at GWAS loci by pooled single-cell CRISPR screens

sc2GWAS: a comprehensive platform linking single cell and GWAS traits of human

Inferring causal cell types of human diseases and risk variants from candidate regulatory elements

Heritability enrichment of specifically expressed genes identifies disease-relevant tissues and cell types

Integrative analysis of omics summary data reveals putative mechanisms underlying complex traits

Supplementary Material For: Inferring Gene-Disease Association by an Integrative Analysis of Eqtl Genome-Wide Association Study and Protein-Protein Interaction Data

Inferring Gene-Disease Association by an Integrative Analysis of eQTL Genome-Wide Association Study and Protein-Protein Interaction Data.

Variant-to-function dissection of rare non-coding GWAS loci with high impact on blood traits

Gene and Pathway-Based Analysis: Second Wave of Genome-wide Association Studies

DNA-binding factor footprints and enhancer RNAs identify functional non-coding genetic variants

Powerful mapping of cis-genetic effects on gene expression across diverse populations reveals novel disease-critical genes