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.

NELFCD and CTSZ loci are associated with jaundice-stage progression in primary biliary cholangitis in the Japanese population

Genome-wide association study identifies TNFSF15 and POU2AF1 as susceptibility loci for primary biliary cirrhosis in the Japanese population

Multiple Genetic Variants Associated with Primary Biliary Cirrhosis in a Han Chinese Population

A genome-wide association study identified PTPN2 as a population-specific susceptibility gene locus for primary biliary cholangitis

[Analysis of disease-pathway by identifying susceptible genes to primary biliary cirrhosis]

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

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

Genome Wide Association Studies of Specific Antinuclear Autoantibody Sub-phenotypes in Primary Biliary Cholangitis.

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

Fine mapping of the MHC region identifies major independent variants associated with Han Chinese primary biliary cholangitis

IL-17A produced by invariant natural killer T cells and CD3(+)CD56(+)alpha Galcer-CD1d tetramer(-) T cells promote liver fibrosis in patients with primary biliary cholangitis

The greater impact of PNPLA3 polymorphism on liver-related events in Japanese non-alcoholic fatty liver disease patients: A multicentre cohort study

MGAT5/TMEM163 Variant is Associated with Prognosis in Ursodeoxycholic Acid-Treated Patients with Primary Biliary Cholangitis

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

Association of CCR6 functional polymorphisms with Primary Biliary Cholangitis

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

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

Keratin Variants Are Overrepresented in Primary Biliary Cirrhosis and Associate with Disease Severity

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

Novel HLA-DRB1 alleles contribute risk for disease susceptibility in primary biliary cholangitis

Corrigendum To: "an International Genome-Wide Meta-Analysis of Primary Biliary Cholangitis: Novel Risk Loci and Candidate Drugs" [J Hepatol 75 (2021) 572-581].