Classical HLA-DRB1 and DPB1 alleles account for HLA associations with primary biliary cirrhosis

Susceptibility to primary biliary cirrhosis (PBC) is strongly associated with human leukocyte antigen (HLA)-region polymorphisms. To determine if associations can be explained by classical HLA determinants, we studied Italian, 676 cases and 1440 controls, genotyped with dense single-nucleotide polymorphisms (SNPs) for which classical HLA alleles and amino acids were imputed. Although previous genome-wide association studies and our results show stronger SNP associations near DQB1, we demonstrate that the HLA signals can be attributed to classical DRB1 and DPB1 genes. Strong support for the predominant role of DRB1 is provided by our conditional analyses. We also demonstrate an independent association of DPB1. Specific HLA-DRB1 genes (08,11 and14) account for most of the DRB1 association signal. Consistent with previous studies, DRB108 (P=1.59 710-11) was the strongest predisposing allele, whereas DRB111 (P=1.42 710-10) was protective. Additionally, DRB114 and the DPB1 association (DPB103:01; P=9.18 710-7) were predisposing risk alleles. No signal was observed in the HLA class 1 or class 3 regions. These findings better define the association of PBC with HLA and specifically support the role of classical HLA-DRB1 and DPB1 genes and alleles in susceptibility to PBC

LLM-PBC: Logic Learning Machine-Based Explainable Rules Accurately Stratify the Genetic Risk of Primary Biliary Cholangitis

Background: The application of Machine Learning (ML) to genetic individual-level data represents a foreseeable advancement for the field, which is still in its infancy. Here, we aimed to evaluate the feasibility and accuracy of an ML-based model for disease risk prediction applied to Primary Biliary Cholangitis (PBC). Methods: Genome-wide significant variants identified in subjects of European ancestry in the recently released second international meta-analysis of GWAS in PBC were used as input data. Quality-checked, individual genomic data from two Italian cohorts were used. The ML included the following steps: import of genotype and phenotype data, genetic variant selection, supervised classification of PBC by genotype, generation of “if-then” rules for disease prediction by logic learning machine (LLM), and model validation in a different cohort. Results: The training cohort included 1345 individuals: 444 were PBC cases and 901 were healthy controls. After pre-processing, 41,899 variants entered the analysis. Several configurations of parameters related to feature selection were simulated. The best LLM model reached an Accuracy of 71.7%, a Matthews correlation coefficient of 0.29, a Youden’s value of 0.21, a Sensitivity of 0.28, a Specificity of 0.93, a Positive Predictive Value of 0.66, and a Negative Predictive Value of 0.72. Thirty-eight rules were generated. The rule with the highest covering (19.14) included the following genes: RIN3, KANSL1, TIMMDC1, TNPO3. The validation cohort included 834 individuals: 255 cases and 579 controls. By applying the ruleset derived in the training cohort, the Area under the Curve of the model was 0.73. Conclusions: This study represents the first illustration of an ML model applied to common variants associated with PBC. Our approach is computationally feasible, leverages individual-level data to generate intelligible rules, and can be used for disease prediction in at-risk individuals

LLM-PBC: Logic Learning Machine-Based Explainable Rules Accurately Stratify the Genetic Risk of Primary Biliary Cholangitis

Background: The application of Machine Learning (ML) to genetic individual-level data represents a foreseeable advancement for the field, which is still in its infancy. Here, we aimed to evaluate the feasibility and accuracy of an ML-based model for disease risk prediction applied to Primary Biliary Cholangitis (PBC). Methods: Genome-wide significant variants identified in subjects of European ancestry in the recently released second international meta-analysis of GWAS in PBC were used as input data. Quality-checked, individual genomic data from two Italian cohorts were used. The ML included the following steps: import of genotype and phenotype data, genetic variant selection, supervised classification of PBC by genotype, generation of “if-then” rules for disease prediction by logic learning machine (LLM), and model validation in a different cohort. Results: The training cohort included 1345 individuals: 444 were PBC cases and 901 were healthy controls. After pre-processing, 41,899 variants entered the analysis. Several configurations of parameters related to feature selection were simulated. The best LLM model reached an Accuracy of 71.7%, a Matthews correlation coefficient of 0.29, a Youden’s value of 0.21, a Sensitivity of 0.28, a Specificity of 0.93, a Positive Predictive Value of 0.66, and a Negative Predictive Value of 0.72. Thirty-eight rules were generated. The rule with the highest covering (19.14) included the following genes: RIN3, KANSL1, TIMMDC1, TNPO3. The validation cohort included 834 individuals: 255 cases and 579 controls. By applying the ruleset derived in the training cohort, the Area under the Curve of the model was 0.73. Conclusions: This study represents the first illustration of an ML model applied to common variants associated with PBC. Our approach is computationally feasible, leverages individual-level data to generate intelligible rules, and can be used for disease prediction in at-risk individuals

Pathway-based analysis of primary biliary cirrhosis genome-wide association studies

Genome-wide association studies (GWAS) have successfully identified several loci associated with primary biliary cirrhosis (PBC) risk. Pathway analysis complements conventional GWAS analysis. We applied the recently developed linear combination test for pathways to datasets drawn from independent PBC GWAS in Italian and Canadian subjects. Of the Kyoto Encyclopedia of Genes and Genomes and BioCarta pathways tested, 25 pathways in the Italian dataset (449 cases, 940 controls) and 26 pathways in the Canadian dataset (530 cases, 398 controls) were associated with PBC susceptibility (P<0.05). After correcting for multiple comparisons, only the eight most significant pathways in the Italian dataset had FDR <0.25 with tumor necrosis factor/stress-related signaling emerging as the top pathway (P=7.38 x 10(-)(4), FDR=0.18). Two pathways, phosphatidylinositol signaling and hedgehog signaling, were replicated in both datasets (P<0.05), and subjected to two additional complementary pathway tests. Both pathway signals remained significant in the Italian dataset on modified gene set enrichment analysis (P<0.05). In both GWAS, variants nominally associated with PBC were significantly overrepresented in the phosphatidylinositol pathway (Fisher exact P<0.05). These results point to established and novel pathway-level associations with inherited predisposition to PBC that, on further independent replication and functional validation, may provide fresh insights into PBC etiology

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

Primary biliary cirrhosis (PBC) is a classical autoimmune liver disease for which effective immunomodulatory therapy is lacking. Here we perform meta-analyses of discovery data sets from genome-wide association studies of European subjects (n=2,764 cases and 10,475 controls) followed by validation genotyping in an independent cohort (n=3,716 cases and 4,261 controls). We discover and validate six previously unknown risk loci for PBC (Pcombined<5 × 10(-8)) and used pathway analysis to identify JAK-STAT/IL12/IL27 signalling and cytokine-cytokine pathways, for which relevant therapies exist

X chromosome contribution to the genetic architecture of primary biliary cholangitis.

BACKGROUND & AIMS: Genome-wide association studies in primary biliary cholangitis (PBC) have failed to find X chromosome (chrX) variants associated with the disease. Here, we specifically explore the chrX contribution to PBC, a sexually dimorphic complex autoimmune disease. METHODS: We performed a chrX-wide association study, including genotype data from 5 genome-wide association studies (from Italy, United Kingdom, Canada, China, and Japan; 5244 case patients and 11,875 control individuals). RESULTS: Single-marker association analyses found approximately 100 loci displaying P < 5 × 10-4, with the most significant being a signal within the OTUD5 gene (rs3027490; P = 4.80 × 10-6; odds ratio [OR], 1.39; 95% confidence interval [CI], 1.028-1.88; Japanese cohort). Although the transethnic meta-analysis evidenced only a suggestive signal (rs2239452, mapping within the PIM2 gene; OR, 1.17; 95% CI, 1.09-1.26; P = 9.93 × 10-8), the population-specific meta-analysis showed a genome-wide significant locus in East Asian individuals pointing to the same region (rs7059064, mapping within the GRIPAP1 gene; P = 6.2 × 10-9; OR, 1.33; 95% CI, 1.21-1.46). Indeed, rs7059064 tags a unique linkage disequilibrium block including 7 genes: TIMM17B, PQBP1, PIM2, SLC35A2, OTUD5, KCND1, and GRIPAP1, as well as a superenhancer (GH0XJ048933 within OTUD5) targeting all these genes. GH0XJ048933 is also predicted to target FOXP3, the main T-regulatory cell lineage specification factor. Consistently, OTUD5 and FOXP3 RNA levels were up-regulated in PBC case patients (1.75- and 1.64-fold, respectively). CONCLUSIONS: This work represents the first comprehensive study, to our knowledge, of the chrX contribution to the genetics of an autoimmune liver disease and shows a novel PBC-related genome-wide significant locus

Classical HLA-DRB1 and DPB1 alleles account for HLA associations with primary biliary cirrhosis.

Susceptibility to primary biliary cirrhosis (PBC) is strongly associated with human leukocyte antigen (HLA)-region polymorphisms. To determine if associations can be explained by classical HLA determinants, we studied Italian, 676 cases and 1440 controls, genotyped with dense single-nucleotide polymorphisms (SNPs) for which classical HLA alleles and amino acids were imputed. Although previous genome-wide association studies and our results show stronger SNP associations near DQB1, we demonstrate that the HLA signals can be attributed to classical DRB1 and DPB1 genes. Strong support for the predominant role of DRB1 is provided by our conditional analyses. We also demonstrate an independent association of DPB1. Specific HLA-DRB1 genes (*08, *11 and *14) account for most of the DRB1 association signal. Consistent with previous studies, DRB1*08 (P=1.59 7 10(-11)) was the strongest predisposing allele, whereas DRB1*11 (P=1.42 7 10(-10)) was protective. Additionally, DRB1*14 and the DPB1 association (DPB1*03:01; P=9.18 7 10(-7)) were predisposing risk alleles. No signal was observed in the HLA class 1 or class 3 regions. These findings better define the association of PBC with HLA and specifically support the role of classical HLA-DRB1 and DPB1 genes and alleles in susceptibility to PBC

X Chromosome Contribution to the Genetic Architecture of Primary Biliary Cholangitis

BACKGROUND & AIMS: Genome-wide association studies in primary biliary cholangitis (PBC) have failed to find X chromosome (chrX) variants associated with the disease. Here, we specifically explore the chrX contribution to PBC, a sexually dimorphic complex autoimmune disease. METHODS: We performed a chrX-wide association study, including genotype data from 5 genome-wide association studies (from Italy, United Kingdom, Canada, China, and Japan; 5244 case patients and 11,875 control individuals). RESULTS: Single-marker association analyses found approximately 100 loci displaying P < 5 x 10(-4), with the most significant being a signal within the OTUD5 gene (rs3027490; P = 4.80 x 10(-6); odds ratio [OR], 1.39; 95% confidence interval [CI], 1.028-1.88; Japanese cohort). Although the transethnic meta-analysis evidenced only a suggestive signal (rs2239452, mapping within the PIM2 gene; OR, 1.17; 95% CI, 1.09-1.26; P = 9.93 x 10(-8)), the population-specific meta-analysis showed a genome-wide significant locus in East Asian individuals pointing to the same region (rs7059064, mapping within the GRIPAP1 gene; P = 6.2 x 10(-9); OR, 1.33; 95% CI, 1.21-1.46). Indeed, rs7059064 tags a unique linkage disequilibrium block including 7 genes: TIMM17B, PQBP1, PIM2, SLC35A2, OTUD5, KCND1, and GRIPAP1, as well as a superenhancer (GH0XJ048933 within OTUD5) targeting all these genes. GH0XJ048933 is also predicted to target FOXP3, the main T-regulatory cell lineage specification factor. Consistently, OTUD5 and FOXP3 RNA levels were up-regulated in PBC case patients (1.75- and 1.64-fold, respectively). CONCLUSIONS: This work represents the first comprehensive study, to our knowledge, of the chrX contribution to the genetics of an autoimmune liver disease and shows a novel PBC-related genome-wide significant locus