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.

ChIP-seq analysis found IL21R, a target gene of GTF2I–the susceptibility gene for primary biliary cholangitis in Chinese Han

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

Association of CCR6 functional polymorphisms with Primary Biliary Cholangitis

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

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

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

Primary Application of Genechip for Detecting Interleukin-1B Susceptibility Gene of Gastric Cancer

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

IL-17 and IL-21 polymorphisms in relation to HBV related hepatocellular carcinoma in Chinese Han population

Association Analysis Between Snps In Il-28b Gene And The Progress Of Hepatitis B Infection In Han Chinese

Genetic variants in IL33 and IL1RL1 genes confer susceptibility to HBV-related liver cirrhosis in Chinese Han population

Genetic variants in five novel loci including CFB and CD40 predispose to chronic hepatitis B.

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

GWAS identifies novel susceptibility loci on 6p21.32 and 21q21.3 for hepatocellular carcinoma in chronic hepatitis B virus carriers.

A regulatory variant at 19p13.3 is associated with primary biliary cholangitis risk and ARID3A expression

Analyzing the gene regulatory network in hepatitis B patients by single-cell ATAC sequencing

[Bioinformatic analysis of differentially expressed genes and immune cell infiltration in primary biliary cholangitis patients with poor response to ursodeoxycholic acid therapy]

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

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

Genome-wide association study identifies new loci associated with risk of HBV infection and disease progression.

ANXA1 is identified as a key gene associated with high risk and T cell infiltration in primary sclerosing cholangitis