Abstract:Abstract Importance Primary biliary cholangitis (PBC) is a classical autoimmune disease, which is highly influenced by genetic determinants. Many genome-wide association studies (GWAS) have reported that numerous genetic loci were significantly associated with PBC susceptibility. However, the effects of genetic determinants on liver cells and its immune microenvironment for PBC remain unclear. Objective To identify genetics-modulated functional liver cell subsets involved in the pathogenesis of PBC. Design, Setting, and Participants In this present study, 13,239 European participants were collected from IEU open GWAS project on PBC. There were 1,124,241 qualified SNPs used for GWAS analysis. Expression quantitative trait loci (eQTL) data across 49 tissues were downloaded from the GTEx database. Two single cell RNA sequencing (scRNA-seq) profiles and two bulk-based RNA transcriptomes were downloaded from the GEO database. Data collection and analyses were performed from August 2020 to June 2021. Main outcomes and measures We constructed a powerful computational framework to integrate GWAS summary statistics with scRNA-seq data to uncover genetics-modulated liver cell subpopulations. Results Based on our multi-omics integrative analysis, we found that 29 risk genes including ORMDL3, GSNK2B , and DDAH2 were significantly associated with PBC susceptibility. Gene-property analysis revealed that four immune cell types, including Cst3 + dendritic cell, Chil3 + macrophage, Trbc2 + T cell, and Gzma + T cell, were significantly enriched by PBC-risk genes. By combining GWAS summary statistics with scRNA-seq data, we found that cholangiocytes exhibited a notable enrichment by PBC-related genetic association signals (Permuted P < 0.05). The risk gene of ORMDL3 showed the highest expression proportion in cholangiocytes than other liver cells (22.38%). Compared with ORMDL3 + cholangiocytes, there were 71 significantly highly-expressed genes among ORMDL3 - cholangiocytes (FDR < 0.05), such as inflammatory cytokine genes CXCL8, CCL3, IFI16 , and IRF1 . These highly-expressed genes were significantly enriched in numerous biological pathways and functional terms associated with autoimmune diseases (FDR < 0.05). Conclusions and relevance To the best of our knowledge, this is the first study to integrate genetic information with single cell sequencing data for parsing genetics-influenced liver cells for PBC risk. We identified that ORMDL3 - cholangiocytes play important immune-modulatory roles in the etiology of PBC. Key points Question Are genetics factors influenced liver cell subpopulations and its immune microenvironment for PBC? Findings In this comprehensive genomics study based on multi-omics data, genetic determinants were significantly enriched in cholangiocytes and immune cells including subsets of macrophage, dendritic cells, and T cells. ORMDL3 - cholangiocytes have crucial immune-modulatory roles in developing PBC. Meaning Findings suggest that integration of single cell sequencing data with GWAS summary statistics contribute to pinpoint PBC-relevant cell types and risk genes.

Mapping of de novo mutations in primary biliary cholangitis to a disease-specific co-expression network underlying homeostasis and metabolism

Chemokine (C-X-C motif) ligand 13 promotes intrahepatic chemokine (C-X-C motif) receptor 5+ lymphocyte homing and aberrant B-cell immune responses in primary biliary cirrhosis.

A Genome-Wide Association Study Identifies Six Novel Risk Loci for Primary Biliary Cholangitis

Epigenetics of Primary Biliary Cholangitis.

Functional investigation and two-sample Mendelian randomization study of primary biliary cholangitis hub genes

Advances in genetic research on primary biliary cholangitis

Unique DUOX2+ACE2+ small cholangiocytes are pathogenic targets for primary biliary cholangitis

International genome-wide meta-analysis identifies new primary biliary cirrhosis risk loci and targetable pathogenic pathways

Integrated GWAS and mRNA Microarray Analysis Identified IFNG and CD40L as the Central Upstream Regulators in Primary Biliary Cholangitis

Epidemiologic and genetic associations between primary biliary cholangitis and extrahepatic rheumatic diseases

Analysis and functional study on lncRNA expression profile in primary biliary cirrhosis patients

Single cell sequencing analysis identifies genetics-modulated ORMDL3+ cholangiocytes having higher metabolic effects on primary biliary cholangitis

Investigating the causal relationship and potential shared diagnostic genes between primary biliary cholangitis and systemic lupus erythematosus using bidirectional Mendelian randomization and transcriptomic analyses

Clonal Characteristics of Circulating B Lymphocyte Repertoire in Primary Biliary Cholangitis

Dissecting causal relationships between primary biliary cholangitis and extrahepatic autoimmune diseases based on Mendelian randomization

A large-scale plasma proteome Mendelian randomization study identifies novel causal plasma proteins related to primary biliary cholangitis

Regulatory T Cells with a Defect in Inhibition on Co-Stimulation Deteriorated Primary Biliary Cholangitis

Combination of single cell sequencing data and GWAS summary statistics reveals genetically-influenced liver cell types for primary biliary cholangitis

Identification of AKR1B10 As a Key Gene in Primary Biliary Cholangitis by Integrated Bioinformatics Analysis and Experimental Validation

The immunogenetics of primary biliary cirrhosis: A comprehensive review