Search

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.

Approximately 10–20% of patients with primary biliary cholangitis (PBC) progress to jaundice stage regardless of treatment with ursodeoxycholic acid and bezafibrate. In this study, we performed a GWAS and a replication study to identify genetic variants associated with jaundice-stage progression in PBC using a total of 1,375 patients (1,202 early-stage and 173 jaundice-stage) in a Japanese population. SNP rs13720, which is located in the 3′UTR of cathepsin Z (CTSZ), showed the strongest association (odds ratio [OR] = 2.15, P = 7.62 × 10−7) with progression to jaundice stage in GWAS. High-density association mapping at the CTSZ and negative elongation factor complex member C/D (NELFCD) loci, which are located within a strong linkage disequilibrium (LD) block, revealed that an intronic SNP of CTSZ, rs163800, was significantly associated with jaundice-stage progression (OR = 2.16, P = 8.57 × 10−8). In addition, eQTL analysis and in silico functional analysis indicated that genotypes of rs163800 or variants in strong LD with rs163800 influence expression levels of both NELFCD and CTSZ mRNA. The present novel findings will contribute to dissect the mechanism of PBC progression and also to facilitate the development of therapies for PBC patients who are resistant to current therapies.

Identification of the primary allele(s) in HLA class II associated diseases remains challenging because of a tight linkage between alleles of HLA-DR and -DQ loci. In the present study, we determined the genotypes of seven HLA loci (HLA-A, -B, -DRB1, -DQA1, -DQB1, -DPA1 and -DPB1) for 1200 Japanese patients with primary biliary cholangitis and 1196 controls. Observation of recombination derivatives facilitated an evaluation of the effects of individual HLA alleles consisting of disease-prone/disease-resistant HLA haplotypes. Consequently, a primary contribution of DQB1*06:04 (odds ratio: 0.19, p = 1.91 × 10−22), DQB1*03:01 (odds ratio: 0.50, p = 6.76 × 10−10), DRB1*08:03 (odds ratio: 1.75, p = 1.01 × 10−7) and DQB1*04:01 (odds ratio: 1.50, p = 9.20 × 10−6) was suggested. Epistasis of the protective DQB1*06:04 to risk conferred by DRB1*08:03 was demonstrated by subpopulation analysis, implicating the presence of an active immunological mechanism that alleviates pathogenic autoimmune reactions. Further, the contribution of the aforementioned HLA alleles as well as an HLA-DP allele, DPB1*02:01 to the association signals of 304 loci among 4103 SNPs in the HLA region at the genome-wide level of significance (p values less than 5 × 10−8) was demonstrated by the stepwise exclusion of the individuals possessing these HLA alleles from the comparison.

Article, Scientific reports.9 :102(2019)

The present study investigated the association between single nucleotide polymorphisms (SNPs) in immune- or DNA repair-related genes and the integration pattern of human papillomavirus (HPV), a promising prognostic marker in cervical cancer. The HPV integration patterns of cervical cancer patients were determined by polymerase chain reaction and in situ hybridization, and categorized as episomal (group A), single-copy or multi-copy tandem repetition integrated (group B), and undetectable HPV types (group C). After sample and SNP quality control, 166,505 SNPs in 161 samples (38, 111, and 12 patients in groups A, B, and C, respectively) were examined. None of the SNPs reached genome-wide significance, and several candidate SNPs for future study were selected, including rs10999435 on chromosome 10q22, rs1322054 on chromosome 9q32-33, and rs10902171 on chromosome 11p15. Luciferase assay identified rs1322054 as the primary functional variant to regulate gene expression in immune cell. Further studies are needed to determine the genetic background of different integration patterns of HPV in cervical cancer patients.

Our recent genome-wide association study found that the NELFCD/CTSZ locus was significantly associated with progression of primary biliary cholangitis (PBC) to jaundice stage in the Japanese population. In this study, we investigated the role of cathepsin Z in the etiology and pathology of PBC. Serum cathepsin Z levels were measured using enzyme-linked immunosorbent assay. The expression and localization of cathepsin Z in liver specimens were analyzed by western blotting and immunohistochemistry. In PBC patients, serum cathepsin Z levels were significantly increased with disease progression. In addition, its levels were positively correlated with alanine transaminase, aspartate transaminase and total bilirubin, and were negatively correlated with platelet count and albumin. Cathepsin Z expression was markedly increased in hepatocytes at later stages of PBC, and its localization was altered from the peri-bile canaliculus to the cytoplasm, where a fraction was no longer colocalized with endosomal/lysosomal vesicles. Similar altered expression of cathepsin Z was observed in end-stage of other cholestatic liver diseases including sepsis, obstructive jaundice, and Alagille syndrome. Our results indicate that altered expression and localization of cathepsin Z in hepatocytes are characteristic features of PBC and other cholestatic liver diseases, and are implicated in the progression of PBC.

The associations of single nucleotide polymorphisms (SNPs) in PLA2R1 and HLA-DQA1, as well as HLA-DRB1*15:01-DQB1*06:02 haplotype with idiopathic membranous nephropathy (IMN) is well known. However, the primary associations of these loci still need to be determined. We used Japanese-specific SNP genotyping array and imputation using 2,048 sequenced Japanese samples to fine-map PLA2R1 region in 98 patients and 413 controls. The most significant SNPs were replicated in a separate sample set of 130 patients and 288 controls. A two-SNP haplotype of intronic and missense SNPs showed the strongest association. The intronic SNP is strongly associated with PLA2R1 expression in the Genotype-Tissue Expression (GTEx) database, and the missense SNP is predicted to alter peptide binding with HLA-DRB1*15:01 by the Immune Epitope Database (IEDB). In HLA region, we performed relative predispositional effect (RPE) tests and identified additional risk alleles in both HLA-DRB1 and HLA-DQB1. We collapsed the risk alleles in each of HLA-DRB1 and HLA-DQB1 into single risk alleles. Reciprocal conditioning of these collapsed risk alleles showed more residual significance for HLA-DRB1 collapsed risk than HLA-DQB1 collapsed risk. These results indicate that changes in the expression levels of structurally different PLA2R protein confer risk for IMN in the presence of risk HLA-DRB1 alleles.

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 −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.

The Jomon and the Yayoi are considered to be the two major ancestral populations of the modern mainland Japanese. The Jomon people, who inhabited mainland Japan, admixed with Yayoi immigrants from the Asian continent. To investigate the population history in the Jomon period (14,500–2,300 years before present [YBP]), we analyzed whole Y-chromosome sequences of 345 Japanese males living in mainland Japan. A phylogenetic analysis of East Asian Y chromosomes identified a major clade (35.4% of mainland Japanese) consisting of only Japanese Y chromosomes, which seem to have originated from indigenous Jomon people. A Monte Carlo simulation indicated that ~70% of Jomon males had Y chromosomes in this clade. The Bayesian skyline plots of 122 Japanese Y chromosomes in the clade detected a marked decrease followed by a subsequent increase in the male population size from around the end of the Jomon period to the beginning of the Yayoi period (2,300 YBP). The colder climate in the Late to Final Jomon period may have resulted in critical shortages of food for the Jomon people, who were hunter-gatherers, and the rice farming introduced by Yayoi immigrants may have helped the population size of the Jomon people to recover.

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 (r2 = 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.