Novel HLA allele associations with susceptibility, staging, symptomatic state, autoimmune hepatitis and hepatocellular carcinoma events for primary biliary cholangitis in the Japanese population.

Primary biliary cholangitis (PBC) is a rare autoimmune disease with a clear predisposition for human leukocyte antigen (HLA)-DR/DQ-associated loss of immune tolerance for the E2 component of the pyruvate dehydrogenase complex. Three-field-resolution HLA imputation of 1,670 Japanese PBC patients and 2,328 healthy controls was conducted using Japanese population-specific HLA reference panels. Eighteen previously reported Japanese PBC-associated HLA alleles were confirmed and extended to 3-field-resolution, including HLA-DRB1*08:03 to HLA-DRB1*08:03:02, HLA-DQB1*03:01 to HLA-DQB1*03:01:01, HLA-DQB1*04:01 to HLA-DQB1*04:01:01 and HLA-DQB1*06:04 to HLA-DQB1*06:04:01. In addition, additional significant novel HLA alleles were identified, including 3 novel susceptible HLA-DQA1 alleles: HLA-DQA1*03:03:01, HLA-DQA1*04:01:01, HLA-DQA1*01:04:01 and 1 novel protective HLA-DQA1 allele, HLA-DQA1*05:05:01. In addition, PBC patients carrying HLA-DRB1*15:01:01 and HLA-DQA1*03:03:01 would have a higher predisposition toward developing concomitant autoimmune hepatitis (AIH). Further, late-stage and symptomatic PBC shared the same susceptible HLA alleles of HLA-A*26:01:01, HLA-DRB1*09:01:02 and HLA-DQB1*03:03:02. Lastly, HLA-DPB1*05:01:01 was identified as a potential risk HLA allele for development of hepatocellular carcinoma (HCC) in PBC patients. In conclusion, we have extended the current knowledge of HLA allele associations to 3-field resolution and identified novel HLA allele associations with predisposition risk, staging, symptomatic state, and AIH and HCC events for Japanese PBC patients.

rs9459874 and rs1012656 in CCR6/FGFR1OP confer susceptibility to primary biliary cholangitis.

Primary biliary cholangitis (PBC) is a chronic cholestatic autoimmune liver disease that appears to be strongly influenced by genetic factors. Recently, an international meta-analysis of genome-wide association studies (GWAS) identified CC-Motif Chemokine Receptor-6 (CCR6) and FGFR1 Oncogene-Partner (FGFR1OP) as PBC-susceptibility genes. However, the lead single nucleotide polymorphisms (SNPs) of CCR6/FGFR1OP showed low linkage disequilibrium with each other in East Asian and European populations. Additionally, the primary functional variants and the molecular mechanisms responsible for PBC-susceptibility remain unclear. Here, among the PBC-susceptibility SNPs identified by high-density association mapping in our previous meta-GWAS (Patients: n = 10,516; healthy controls: n = 20,772) within the CCR6/FGFR1OP locus, rs9459874 and rs1012656 were identified as primary functional variants. These functional variants accounted for the effects of GWAS-identified lead SNPs in CCR6/FGFR1OP. Additionally, the roles of rs9459874 and rs1012656 in regulating FGFR1OP transcription and CCR6 translation, respectively, were supported by expression quantitative trait loci (eQTL) analysis and gene editing technology using the CRISPR/Cas9 system. Immunohistochemistry showed higher expression of CCR6 protein in the livers of patients with PBC than in those of a non-diseased control. In conclusion, we identified primary functional variants in CCR6/FGFR1OP and revealed the molecular mechanisms by which these variants confer PBC-susceptibility in an eQTL-dependent or -independent manner. The approach in this study is applicable for the elucidation of the pathogenesis of other autoimmune disorders in which CCR6/FGFR1OP is known as a susceptibility locus, as well as 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 × 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.

rs1944919 on chromosome 11q23.1 and its effector genes COLCA1/COLCA2 confer susceptibility to primary biliary cholangitis.

Primary biliary cholangitis (PBC) is a chronic, progressive cholestatic liver disease in which intrahepatic bile ducts are destroyed by an autoimmune reaction. Our previous genome-wide association study (GWAS) identified chromosome 11q23.1 as a susceptibility gene locus for PBC in the Japanese population. Here, high-density association mapping based on single nucleotide polymorphism (SNP) imputation and in silico/in vitro functional analyses identified rs1944919 as the primary functional variant. Expression-quantitative trait loci analyses showed that the PBC susceptibility allele of rs1944919 was significantly associated with increased COLCA1/COLCA2 expression levels. Additionally, the effects of rs1944919 on COLCA1/COLCA2 expression levels were confirmed using genotype knock-in versions of cell lines constructed using the CRISPR/Cas9 system and differed between rs1944919-G/G clones and -T/T clones. To our knowledge, this is the first study to demonstrate the contribution of COLCA1/COLCA2 to PBC susceptibility.

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

Genome-wide association studies (GWASs) in European and East Asian populations have identified more than 40 disease-susceptibility genes in primary biliary cholangitis (PBC). The aim of this study is to computationally identify disease pathways, upstream regulators, and therapeutic targets in PBC through integrated GWAS and messenger RNA (mRNA) microarray analysis. Disease pathways and upstream regulators were analyzed with ingenuity pathway analysis in data set 1 for GWASs (1,920 patients with PBC and 1,770 controls), which included 261 annotated genes derived from 6,760 single-nucleotide polymorphisms (P < 0.00001), and data set 2 for mRNA microarray analysis of liver biopsy specimens (36 patients with PBC and 5 normal controls), which included 1,574 genes with fold change >2 versus controls (P < 0.05). Hierarchical cluster analysis and categorization of cell type-specific genes were performed for data set 2. There were 27 genes, 10 pathways, and 149 upstream regulators that overlapped between data sets 1 and 2. All 10 pathways were immune-related. The most significant common upstream regulators associated with PBC disease susceptibility identified were interferon-gamma (IFNG) and CD40 ligand (CD40L). Hierarchical cluster analysis of data set 2 revealed two distinct groups of patients with PBC by disease activity. The most significant upstream regulators associated with disease activity were IFNG and CD40L. Several molecules expressed in B cells, T cells, Kupffer cells, and natural killer-like cells were identified as potential therapeutic targets in PBC with reference to a recently reported list of cell type-specific gene expression in the liver. Conclusion: Our integrated analysis using GWAS and mRNA microarray data sets predicted that IFNG and CD40L are the central upstream regulators in both disease susceptibility and activity of PBC and identified potential downstream therapeutic targets.

NFKB1 and MANBA Confer Disease Susceptibility to Primary Biliary Cholangitis via Independent Putative Primary Functional Variants.

BACKGROUND & AIMS: Primary biliary cholangitis (PBC) is a chronic and cholestatic liver disease that eventually leads to cirrhosis and hepatic failure. We recently identified several susceptibility genes included NFKB1 and MANBA for PBC in the Japanese population by genome-wide association study. However, the primary functional variants in the NFKB1/MANBA region and the molecular mechanism for conferring disease susceptibility to PBC have not yet been clarified. METHODS: We performed high-density association mapping based on a single-nucleotide polymorphism (SNP) imputation analysis, using data from a whole-genome sequence reference panel of 1070 Japanese individuals and the previous genome-wide association study (1389 PBC patients, 1508 healthy controls). Among SNPs (P < 5.0 × 10-7) in the NFKB1/MANBA region, putative primary functional variants and the molecular mechanism for conferring disease susceptibility to PBC were identified by in silico/in vitro functional analysis. RESULTS: Among the SNPs in the NFKB1/MANBA region, rs17032850 and rs227361, which changed the binding of transcription factors lymphoid enhancer-binding factor 1 (LEF-1) and retinoid X receptor α (RXRα), respectively, were identified as putative primary functional variants that regulate gene expression. In addition, expression-quantitative trait locus data and gene editing using a clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9 system supported the potential role of rs17032850 and rs227361 in regulating NFKB1 and MANBA expression, respectively. CONCLUSIONS: We identified independent putative primary functional variants in NFKB1/MANBA and showed the distinct molecular mechanism by which each putative primary functional variant conferred susceptibility to PBC. Our approach was useful to dissect the pathogenesis not only of PBC, but also other digestive diseases in which NFKB1/MANBA has been reported as a susceptibility locus.

Identification of the functional variant driving ORMDL3 and GSDMB expression in human chromosome 17q12-21 in primary biliary cholangitis.

Numerous genome-wide association studies (GWAS) have been performed to identify susceptibility genes to various human complex diseases. However, in many cases, neither a functional variant nor a disease susceptibility gene have been clarified. Here, we show an efficient approach for identification of a functional variant in a primary biliary cholangitis (PBC)-susceptible region, chromosome 17q12-21 (ORMDL3-GSDMB-ZPBP2-IKZF3). High-density association mapping was carried out based on SNP imputation analysis by using the whole-genome sequence data from a reference panel of 1,070 Japanese individuals (1KJPN), together with genotype data from our previous GWAS (PBC patients: n = 1,389; healthy controls: n = 1,508). Among 23 single nucleotide polymorphisms (SNPs) with P < 1.0 × 10-8, rs12946510 was identified as the functional variant that influences gene expression via alteration of Forkhead box protein O1 (FOXO1) binding affinity in vitro. Moreover, expression-quantitative trait locus (e-QTL) analyses showed that the PBC susceptibility allele of rs12946510 was significantly associated with lower endogenous expression of ORMDL3 and GSDMB in whole blood and spleen. This study not only identified the functional variant in chr.17q12-21 and its molecular mechanism through which it conferred susceptibility to PBC, but it also illustrated an efficient systematic approach for post-GWAS analysis that is applicable to other complex diseases.

Fine-mapping analysis revealed complex pleiotropic effect and tissue-specific regulatory mechanism of TNFSF15 in primary biliary cholangitis, Crohn's disease and leprosy.

Genetic polymorphism within the 9q32 locus is linked with increased risk of several diseases, including Crohn's disease (CD), primary biliary cholangitis (PBC) and leprosy. The most likely disease-causing gene within 9q32 is TNFSF15, which encodes the pro-inflammatory cytokine TNF super-family member 15, but it was unknown whether these disparate diseases were associated with the same genetic variance in 9q32, and how variance within this locus might contribute to pathology. Using genetic data from published studies on CD, PBC and leprosy we revealed that bearing a T allele at rs6478108/rs6478109 (r(2) = 1) or rs4979462 was significantly associated with increased risk of CD and decreased risk of leprosy, while the T allele at rs4979462 was associated with significantly increased risk of PBC. In vitro analyses showed that the rs6478109 genotype significantly affected TNFSF15 expression in cells from whole blood of controls, while functional annotation using publicly-available data revealed the broad cell type/tissue-specific regulatory potential of variance at rs6478109 or rs4979462. In summary, we provide evidence that variance within TNFSF15 has the potential to affect cytokine expression across a range of tissues and thereby contribute to protection from infectious diseases such as leprosy, while increasing the risk of immune-mediated diseases including CD and PBC.

Disease susceptibility genes shared by primary biliary cirrhosis and Crohn's disease in the Japanese population.

We previously identified TNFSF15 as the most significant susceptibility gene at non-HLA loci for both primary biliary cirrhosis (PBC) and Crohn's diseases (CD) in the Japanese population. The aim of this study is to identify further disease susceptibility genes shared by PBC and CD. We selected 15 and 33 genetic variants that were significantly associated with PBC and CD, respectively, based on previously reported genome-wide association studies of the Japanese population. Next, an association study was independently performed for these genetic variants in CD (1312 CD patients and 3331 healthy controls) and PBC (1279 PBC patients and 1015 healthy controls) cohorts. Two CD susceptibility genes, ICOSLG rs2838519 and IL12B rs6556412, were also nominally associated with susceptibility to PBC (P=3.85 × 10(-2) and P=8.40 × 10(-3), respectively). Three PBC susceptibility genes, CXCR5 rs6421571, STAT4 rs7574865 and NFKB1 rs230534, were nominally associated with susceptibility to CD (P=2.82 × 10(-2), P=3.88 × 10(-2) and P=2.04 × 10(-2), respectively). The effect of ICOSLG and CXCR5 variants were concordant but the effect of STAT4, NFKB1 and IL12B variants were discordant for PBC and CD. TNFSF15 and ICOSLG-CXCR5 might constitute a shared pathogenic pathway in the development of PBC and CD in the Japanese population, whereas IL12B-STAT4-NFKB1 might constitute an opposite pathogenic pathway, reflecting the different balance between Th1 and Th17 in the two diseases.

Human primary biliary cirrhosis-susceptible allele of rs4979462 enhances TNFSF15 expression by binding NF-1.

A genome-wide association study (GWAS) identified tumor necrosis factor superfamily member 15 (TNFSF15) as the strongest associated gene with susceptibility to primary biliary cirrhosis (PBC) outside the HLA loci in the Japanese population. However, causal functional variants of the TNFSF15 locus and the molecular mechanism underlying disease susceptibility have not been clarified. Here, to identify the functional causal variants of the TNFSF15 locus, integrated analysis comprising in silico analysis, a case-control association study and in vitro functional analysis was performed. Initially, 32 functional candidate single-nucleotide polymorphisms (SNPs) in the expression regulatory motifs, the coding region, or the untranslated regions (UTRs) of the TNFSF15 locus were selected by in silico analysis. By the case-control association studies using PBC patients (n = 1279) and healthy controls (n = 1091) in the Japanese population, rs4979462 [P = 1.85 × 10(-14) (our previous study)], rs56211063 (P = 2.21 × 10(-14)), and rs55768522 (r(2) = 1 with rs4979462) were likely candidates for causal variants. Among these SNPs, rs4979462 was identified as the causal variant by in vitro functional analysis using luciferase assay and electrophoretic mobility shift assay (EMSA). Super-shift assay clarified that PBC-susceptible allele of rs4979462 generated a novel NF-1 binding site. Moreover, higher endogenous TNFSF15 protein and mRNA expression levels were observed in individuals with the PBC-susceptible allele of rs4979462. This study identified the causal variant for PBC susceptibility in the TNFSF15 locus and clarified its underlying molecular mechanism. TNFSF15 and NF-1 are considered to be potential targets for the treatment of PBC.

New susceptibility and resistance HLA-DP alleles to HBV-related diseases identified by a trans-ethnic association study in Asia.

Previous studies have revealed the association between SNPs located on human leukocyte antigen (HLA) class II genes, including HLA-DP and HLA-DQ, and chronic hepatitis B virus (HBV) infection, mainly in Asian populations. HLA-DP alleles or haplotypes associated with chronic HBV infection or disease progression have not been fully identified in Asian populations. We performed trans-ethnic association analyses of HLA-DPA1, HLA-DPB1 alleles and haplotypes with hepatitis B virus infection and disease progression among Asian populations comprising Japanese, Korean, Hong Kong, and Thai subjects. To assess the association between HLA-DP and chronic HBV infection and disease progression, we conducted high-resolution (4-digit) HLA-DPA1 and HLA-DPB1 genotyping in a total of 3,167 samples, including HBV patients, HBV-resolved individuals and healthy controls. Trans-ethnic association analyses among Asian populations identified a new risk allele HLA-DPB1*09 ∶ 01 (P = 1.36 × 10(-6); OR= 1.97; 95% CI, 1.50-2.59) and a new protective allele DPB1*02 ∶ 01 (P = 5.22 × 10(-6); OR = 0.68; 95% CI, 0.58-0.81) to chronic HBV infection, in addition to the previously reported alleles. Moreover, DPB1*02 ∶ 01 was also associated with a decreased risk of disease progression in chronic HBV patients among Asian populations (P = 1.55 × 10(-7); OR = 0.50; 95% CI, 0.39-0.65). Trans-ethnic association analyses identified Asian-specific associations of HLA-DP alleles and haplotypes with HBV infection or disease progression. The present findings will serve as a base for future functional studies of HLA-DP molecules in order to understand the pathogenesis of HBV infection and the development of hepatocellular carcinoma.

Systemic and local expression levels of TNF-like ligand 1A and its decoy receptor 3 are increased in primary biliary cirrhosis.

BACKGROUND & AIMS: Through a genome-wide association study of a Japanese population, we recently identified TNFSF15, a gene encoding TNF-like ligand 1A (TL1A), as a susceptibility gene for primary biliary cirrhosis (PBC). We investigated the clinical significance of TL1A and one of its receptors, decoy receptor 3 (DcR3), in PBC. METHODS: We analysed the systemic and local expression of TL1A and DcR3 in 110 PBC patients and 46 healthy controls using enzyme-linked immunosorbent assay, quantitative polymerase chain reaction and immunohistochemical staining. RESULTS: Serum TL1A levels were significantly increased in PBC patients at both early and late stages as compared with healthy controls, and its levels were significantly decreased in early-stage PBC patients after ursodeoxycholic acid (UDCA) treatment. TL1A was immunohistochemically localized to biliary epithelial cells, Kupffer cells, blood vessels and infiltrating mononuclear cells in the PBC liver. In addition, TL1A messenger RNA expression was increased in the PBC liver as compared with the non-diseased liver. Serum DcR3 levels were also significantly increased in PBC patients, and were significantly decreased after UDCA treatment in early-stage PBC patients. CONCLUSIONS: These results indicate that TL1A and DcR3 may play an important role in the pathogenesis of PBC.

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

For the identification of susceptibility loci for primary biliary cirrhosis (PBC), a genome-wide association study (GWAS) was performed in 963 Japanese individuals (487 PBC cases and 476 healthy controls) and in a subsequent replication study that included 1,402 other Japanese individuals (787 cases and 615 controls). In addition to the most significant susceptibility region, human leukocyte antigen (HLA), we identified two significant susceptibility loci, TNFSF15 (rs4979462) and POU2AF1 (rs4938534) (combined odds ratio [OR] = 1.56, p = 2.84 × 10(-14) for rs4979462, and combined OR = 1.39, p = 2.38 × 10(-8) for rs4938534). Among 21 non-HLA susceptibility loci for PBC identified in GWASs of individuals of European descent, three loci (IL7R, IKZF3, and CD80) showed significant associations (combined p = 3.66 × 10(-8), 3.66 × 10(-9), and 3.04 × 10(-9), respectively) and STAT4 and NFKB1 loci showed suggestive association with PBC (combined p = 1.11 × 10(-6) and 1.42 × 10(-7), respectively) in the Japanese population. These observations indicated the existence of ethnic differences in genetic susceptibility loci to PBC and the importance of TNF signaling and B cell differentiation for the development of PBC in individuals of European descent and Japanese individuals.

Genome-wide association study confirming association of HLA-DP with protection against chronic hepatitis B and viral clearance in Japanese and Korean.

Hepatitis B virus (HBV) infection can lead to serious liver diseases, including liver cirrhosis (LC) and hepatocellular carcinoma (HCC); however, about 85-90% of infected individuals become inactive carriers with sustained biochemical remission and very low risk of LC or HCC. To identify host genetic factors contributing to HBV clearance, we conducted genome-wide association studies (GWAS) and replication analysis using samples from HBV carriers and spontaneously HBV-resolved Japanese and Korean individuals. Association analysis in the Japanese and Korean data identified the HLA-DPA1 and HLA-DPB1 genes with P(meta) = 1.89×10⁻¹² for rs3077 and P(meta) = 9.69×10⁻¹° for rs9277542. We also found that the HLA-DPA1 and HLA-DPB1 genes were significantly associated with protective effects against chronic hepatitis B (CHB) in Japanese, Korean and other Asian populations, including Chinese and Thai individuals (P(meta) = 4.40×10⁻¹9 for rs3077 and P(meta) = 1.28×10⁻¹5 for rs9277542). These results suggest that the associations between the HLA-DP locus and the protective effects against persistent HBV infection and with clearance of HBV were replicated widely in East Asian populations; however, there are no reports of GWAS in Caucasian or African populations. Based on the GWAS in this study, there were no significant SNPs associated with HCC development. To clarify the pathogenesis of CHB and the mechanisms of HBV clearance, further studies are necessary, including functional analyses of the HLA-DP molecule.

Establishment of human cell lines showing circadian rhythms of bioluminescence.

We have established human retinal pigment epithelial cell lines stably expressing the luciferase gene, driven by the human Bmal1 promoter, to obtain human-derived cells that show circadian rhythms of bioluminescence after dexamethasone treatment. The average circadian period of bioluminescence for the obtained clones was 24.07+/-0.48 h. Lithium (10 mM) in the medium significantly lengthened the circadian period of bioluminescence, which is consistent with previous reports, while 2 mM or 5 mM lithium had no effect. This is the first report on the establishment of human-derived cell lines that proliferate infinitely and show circadian rhythms of bioluminescence, and also the first to investigate the effects of low-dose lithium on the circadian rhythms of human-derived cells in vitro. The established cells will be useful for various in vitro studies of human circadian rhythms and for the development of new therapies for human disorders related to circadian rhythm disturbances.