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.

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

Disease predisposition of human leukocyte antigen class II genes influences the gut microbiota composition in patients with primary biliary cholangitis

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

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

Genome-Wide Association Studies in Primary Biliary Cirrhosis

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

The immunogenetics of primary biliary cirrhosis: A comprehensive review

Genetic link between primary biliary cholangitis and connective tissue diseases in European populations: A two-sample Mendelian randomization study

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

HLA class II alleles differing by a single amino acid associate with clinical phenotype and outcome in patients with primary membranous nephropathy.

Association of CCR6 functional polymorphisms with Primary Biliary Cholangitis

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

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

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

Fine Mapping Analysis of HLA-DP/DQ Gene Clusters on Chromosome 6 Reveals Multiple Susceptibility Loci for HBV Infection

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

Fine Mapping of the MHC Region Identifies Novel Variants Associated with HBV-Related Hepatocellular Carcinoma in Han Chinese

Replication Study and Meta-Analysis Indicate a Suggestive Association of RUNX3 Locus with Primary Biliary Cholangitis

Causal association between systemic lupus erythematosus and primary biliary cholangitis: A bidirectional Mendelian randomization study

Genome Wide Association Studies of Specific Antinuclear Autoantibody Sub-phenotypes in 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