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.

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

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

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

Genome-Wide Association Studies in Primary Biliary Cirrhosis

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

Clinical Analysis of Primary Biliary Cirrhosis: a Report of 31 Cases

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.

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

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

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

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

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

Association of CCR6 functional polymorphisms with Primary Biliary Cholangitis

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

Clonal Characteristics of Circulating B Lymphocyte Repertoire in Primary Biliary Cholangitis

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

The immunogenetics of primary biliary cirrhosis: A comprehensive review

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

Single-cell RNA sequencing reveals the pro-inflammatory roles of liver-resident Th1-like cells in primary biliary cholangitis

Genome-wide Association Study Identifies 8p21.3 Associated with Persistent Hepatitis B Virus Infection among Chinese.