Article Text

Download PDFPDF

Subtype-specific gout susceptibility loci and enrichment of selection pressure on ABCG2 and ALDH2 identified by subtype genome-wide meta-analyses of clinically defined gout patients
  1. Akiyoshi Nakayama1,2,
  2. Masahiro Nakatochi3,
  3. Yusuke Kawamura1,4,
  4. Ken Yamamoto5,
  5. Hirofumi Nakaoka6,
  6. Seiko Shimizu1,
  7. Toshihide Higashino1,7,
  8. Teruhide Koyama8,
  9. Asahi Hishida9,
  10. Kiyonori Kuriki10,
  11. Miki Watanabe11,
  12. Toru Shimizu12,13,
  13. Keiko Ooyama14,
  14. Hiroshi Ooyama14,
  15. Mitsuo Nagase15,
  16. Yuji Hidaka16,
  17. Daisuke Matsui8,
  18. Takashi Tamura9,
  19. Takeshi Nishiyama11,
  20. Chisato Shimanoe17,18,
  21. Sakurako Katsuura-Kamano19,
  22. Naoyuki Takashima20,21,
  23. Yuya Shirai22,23,
  24. Makoto Kawaguchi1,24,
  25. Mikiya Takao1,25,
  26. Ryo Sugiyama1,
  27. Yuzo Takada26,
  28. Takahiro Nakamura27,
  29. Hiroshi Nakashima28,
  30. Masashi Tsunoda28,
  31. Inaho Danjoh29,
  32. Atsushi Hozawa30,
  33. Kazuyoshi Hosomichi31,
  34. Yu Toyoda32,
  35. Yu Kubota32,
  36. Tappei Takada32,
  37. Hiroshi Suzuki32,
  38. Blanka Stiburkova33,34,
  39. Tanya J. Major35,
  40. Tony R. Merriman35,
  41. Nagato Kuriyama8,
  42. Haruo Mikami36,
  43. Toshiro Takezaki37,
  44. Keitaro Matsuo38,39,
  45. Sadao Suzuki11,
  46. Tatsuo Hosoya40,41,
  47. Yoichiro Kamatani42,43,
  48. Michiaki Kubo44,
  49. Kimiyoshi Ichida40,45,
  50. Kenji Wakai9,
  51. Ituro Inoue6,
  52. Yukinori Okada22,46,
  53. Nariyoshi Shinomiya1,
  54. Hirotaka Matsuo1
  55. on behalf of Japan Gout Genomics Consortium (Japan Gout)
    1. 1 Department of Integrative Physiology and Bio-Nano Medicine, National Defense Medical College, Tokorozawa, Japan
    2. 2 Medical Squadron, Air Base Group, Western Aircraft Control and Warning Wing, Japan Air Self-Defense Force, Kasuga, Japan
    3. 3 Division of Department of Nursing, Nagoya University Graduate School of Medicine, Nagoya, Japan
    4. 4 Department of General Medicine, National Defense Medical College, Tokorozawa, Japan
    5. 5 Department of Medical Biochemistry, Kurume University School of Medicine, Kurume, Japan
    6. 6 Division of Human Genetics, Department of Integrated Genetics, National Institute of Genetics, Mishima, Japan
    7. 7 Graduate School of Information Science and Technology, Hokkaido University, Sapporo, Japan
    8. 8 Department of Epidemiology for Community Health and Medicine, Kyoto Prefectural University of Medicine, Kyoto, Japan
    9. 9 Department of Preventive Medicine, Nagoya University Graduate School of Medicine, Nagoya, Japan
    10. 10 Laboratory of Public Health, School of Food and Nutritional Sciences, University of Shizuoka, Shizuoka, Japan
    11. 11 Department of Public Health, Nagoya City University Graduate School Medical Science, Nagoya, Japan
    12. 12 Midorigaoka Hospital, Takatsuki, Japan
    13. 13 Kyoto Industrial Health Association, Kyoto, Japan
    14. 14 Ryougoku East Gate Clinic, Tokyo, Japan
    15. 15 Nagase Clinic, Tokyo, Japan
    16. 16 Akasaka Central Clinic, Tokyo, Japan
    17. 17 Department of Preventive Medicine, Faculty of Medicine, Saga University, Saga, Japan
    18. 18 Clinical Research Center, Saga University Hospital, Saga, Japan
    19. 19 Department of Preventive Medicine, Institute of Biomedical Sciences, Tokushima University Graduate School, Tokushima, Japan
    20. 20 Department of Health Science, Shiga University of Medical Science, Otsu, Japan
    21. 21 Department of Public Health, Faculty of Medicine, Kindai University, Osaka-Sayama, Japan
    22. 22 Department of Statistical Genetics, Osaka University Graduate School of Medicine, Suita, Japan
    23. 23 Department of Respiratory Medicine and Clinical Immunology, Osaka University Graduate School of Medicine, Suita, Japan
    24. 24 Department of Urology, National Defense Medical College, Tokorozawa, Japan
    25. 25 Department of Surgery, National Defense Medical College, Tokorozawa, Japan
    26. 26 Faculty of Medical Science, Teikyo University of Science, Tokyo, Japan
    27. 27 Laboratory for Mathematics, National Defense Medical College, Tokorozawa, Japan
    28. 28 Department of Preventive Medicine and Public Health, National Defense Medical College, Tokorozawa, Japan
    29. 29 Group of Privacy Controls, Tohoku Medical Megabank Organization, Sendai, Japan
    30. 30 Department of Preventive Medicine and Epidemiology, Tohoku Medical Megabank Organization, Tohoku University, Sendai, Japan
    31. 31 Department of Bioinformatics and Genomics, Graduate School of Advanced Preventive Medical Sciences, Kanazawa University, Kanazawa, Japan
    32. 32 Department of Pharmacy, The University of Tokyo Hospital, Tokyo, Japan
    33. 33 Department of Pediatrics and Adolescent Medicine, First Faculty of Medicine, Charles University and General University Hospital, Prague, Czech Republic
    34. 34 Institute of Rheumatology, Prague, Czech Republic
    35. 35 Department of Biochemisty, University of Otago, Dunedin, New Zealand
    36. 36 Cancer Prevention Center, Chiba Cancer Center Research Institute, Chiba, Japan
    37. 37 Department of International Island and Community Medicine, Kagoshima University Graduate School of Medical and Dental Sciences, Kagoshima, Japan
    38. 38 Division of Cancer Epidemiology and Prevention, Aichi Cancer Center Research Institute, Nagoya, Japan
    39. 39 Department of Epidemiology, Nagoya University Graduate School of Medicine, Nagoya, Japan
    40. 40 Division of Kidney and Hypertension, Department of Internal Medicine, Jikei University School of Medicine, Tokyo, Japan
    41. 41 Department of Pathophysiology and Therapy in Chronic Kidney Disease, Jikei University School of Medicine, Tokyo, Japan
    42. 42 Laboratory for Statistical Analysis, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
    43. 43 Center for Genomic Medicine, Kyoto University Graduate School of Medicine, Kyoto, Japan
    44. 44 RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
    45. 45 Department of Pathophysiology, Tokyo University of Pharmacy and Life Science, Hachioji, Japan
    46. 46 Laboratory of Statistical Immunology, Immunology Frontier Research Center (WPI-IFReC), Osaka University, Suita, Japan
    1. Correspondence to Hirotaka Matsuo, Department of Integrative Physiology and Bio-Nano Medicine, National Defense Medical College, Tokorozawa 359-8513, Japan; hmatsuo{at}ndmc.ac.jp

    Abstract

    Objectives Genome-wide meta-analyses of clinically defined gout were performed to identify subtype-specific susceptibility loci. Evaluation using selection pressure analysis with these loci was also conducted to investigate genetic risks characteristic of the Japanese population over the last 2000–3000 years.

    Methods Two genome-wide association studies (GWASs) of 3053 clinically defined gout cases and 4554 controls from Japanese males were performed using the Japonica Array and Illumina Array platforms. About 7.2 million single-nucleotide polymorphisms were meta-analysed after imputation. Patients were then divided into four clinical subtypes (the renal underexcretion type, renal overload type, combined type and normal type), and meta-analyses were conducted in the same manner. Selection pressure analyses using singleton density score were also performed on each subtype.

    Results In addition to the eight loci we reported previously, two novel loci, PIBF1 and ACSM2B, were identified at a genome-wide significance level (p<5.0×10–8) from a GWAS meta-analysis of all gout patients, and other two novel intergenic loci, CD2-PTGFRN and SLC28A3-NTRK2, from normal type gout patients. Subtype-dependent patterns of Manhattan plots were observed with subtype GWASs of gout patients, indicating that these subtype-specific loci suggest differences in pathophysiology along patients’ gout subtypes. Selection pressure analysis revealed significant enrichment of selection pressure on ABCG2 in addition to ALDH2 loci for all subtypes except for normal type gout.

    Conclusions Our findings on subtype GWAS meta-analyses and selection pressure analysis of gout will assist elucidation of the subtype-dependent molecular targets and evolutionary involvement among genotype, phenotype and subtype-specific tailor-made medicine/prevention of gout and hyperuricaemia.

    • gout/hyperuricaemia
    • genome-wide association study (GWAS)
    • Japanese
    • subtype specific locus
    • selection pressure analysis
    https://creativecommons.org/licenses/by/4.0/

    This is an open access article distributed in accordance with the Creative Commons Attribution 4.0 Unported (CC BY 4.0) license, which permits others to copy, redistribute, remix, transform and build upon this work for any purpose, provided the original work is properly cited, a link to the licence is given, and indication of whether changes were made. See: https://creativecommons.org/licenses/by/4.0/.

    Statistics from Altmetric.com

    Request Permissions

    If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

    Key messages

    What is already known about this subject?

    • Our previous genome-wide association study (GWAS) was performed on broad subtypes of gout with only 945 gout cases. A recent study has revealed genetic adaptive evolution of gout in the Japanese population.

    What does this study add?

    • This is the first GWAS meta-analyses of clinically defined gout with more finely differentiated subtypes using two GWAS platforms with larger samples (3055 cases and 4554 controls). We identified multiple subtype-specific loci including four novel loci such as CD2, which encodes a well-known surface antigen found on all peripheral blood T-cells.

    • The present study showed significant enrichment of selection pressure on two genes, ABCG2 and ALDH2, for gout susceptibility in the Japanese population over the last 2000–3000 years.

    Key messages

    How might this impact on clinical practice or future developments?

    • Our subtype GWASs of gout enabled us to develop subtype-dependent molecular targets that will lead to novel subtype-specific genome tailor-made therapies for gout/hyperuricaemia.

    • The present study also elucidates the Japanese genetic evolution of susceptibility to gout/hyperuricaemia and its subtypes.

    Introduction

    Gout is a well-known disease that manifests as acute and severe non-infectious arthritis.1 2 According to patients’ clinical parameters which reflect its causes,2–5 gout can be classified into four distinct subtypes: the renal underexcretion (RUE) type, renal overload (ROL) type, combined type and normal type, as shown in table 1 and online supplementary figure S1. Because these subtypes reflect causes of gout, genome-wide association studies (GWASs) of these subtypes are also likely to indicate its various genetic and pathophysiological backgrounds. While dividing patients into these subtypes is helpful for understanding patients’ pathophysiology, GWASs of these subtypes have only rarely been conducted, partly because clinical data, including time-consuming urinary collection, are necessary to categorise these subtypes. We previously performed a GWAS with clinically defined gout patients,6 followed by another with broader subtypes:7 RUE gout and ROL gout (table 1 and online supplementary figure S1), that revealed their specific loci. Although we were able to show these associations, this process has its limitations, including the use of a custom chip for replication studies that did not provide comprehensive genetic association searching. We use finely differentiated subtypes in daily clinical settings but there were not sufficient numbers of patients in the previous study7 to enable a GWAS with these finely differentiated subtypes. This prompted us to conduct, for the first time, GWASs with four distinct subtypes using meta-analysis across two GWAS platforms with a larger number of patients. We additionally conducted selection pressure analysis of the Japanese population on gout subtypes with the risk loci identified in the present study in order to investigate the evolutionary selective pressure on the Japanese population over the last 2000–3000 years.

    Table 1

    Subtypes of gout* used in the present study

    Methods

    Study subjects and patients involvement

    We performed subtype genome-wide meta-analyses based on two case–control data sets for gout that included the Japonica Array8 and Illumina Array platforms. Patients with known clinical parameters were recruited from Japanese male outpatients at gout clinics (see online supplementary methods). All 3104 cases were clinically diagnosed as having primary gout according to the criteria established by the American College of Rheumatology,9 and their subtypes were also diagnosed along with their clinical parameters as described previously3 5 6 (table 1 and online supplementary figure S1). As controls, 6081 individuals were assigned from Japanese male participants in the Japan Multi-Institutional Collaborative Cohort Study (J-MICC Study).10 11 This research was done without patient involvement (see online supplementary methods).

    Genotyping and imputation for the Japonica Array data set

    A total of 1048 male clinically defined gout cases and 1179 male controls from the J-MICC Study10 11 were genotyped with the use of a Japonica SNP Array.8 The detail of quality control is described in online supplementary methods. This quality control filtering resulted in the selection of 1028 case subjects and 1167 control subjects as well as 603 009 single-nucleotide polymorphisms (SNPs). Prephasing and imputation were performed using SHAPEIT212 and Minimac3,13 respectively. Postimputation quality control was also performed as described in the online supplementary methods. Ultimately, 1028 case subjects and 952 control subjects as well as 7 529 176 SNPs remained for the GWAS analysis.

    Genotyping and imputation for the Illumina Array data set

    As case data, 2056 male gout cases subjects were genotyped with the use of HumanOmniExpress or HumanOmniExpressExome BeadChip Arrays (Illumina, San Diego, CA, USA). The detail of quality control is described in the online supplementary methods. This quality control filtering resulted in the selection of 2032 case subjects and 4901 control subjects as well as 553 321 SNPs. Postimputation quality control was also performed as described in the online supplementary methods. Ultimately, 2025 case subjects and 3602 control subjects as well as 7 356 207 SNPs remained for the GWAS analysis.

    Association analysis for SNPs and gout

    The association of SNPs with gout was assessed using logistic regression analysis (generalised linear model); the dependent variable was gout label (case=1, control=0), and the independent variables included imputed genotypes of each SNP and covariates. The covariates comprised the first four principal component scores. The effect sizes and standard errors estimated in logistic regression analysis were used in the subsequent meta-analysis. The association analysis was performed with the use of Efficient and Parallelizable Association Container Toolbox (EPACTS). https://genome.sph.umich.edu/wiki/EPACTS).

    Meta-analyses

    The meta-analyses were performed using a total of 3053 cases and 4554 controls from the two data sets (online supplementary table S1–S3). The association results for each SNP across the studies were combined with METAL software14 using the fixed-effects inverse-variance-weighted method. Heterogeneity of effect sizes was assessed via the I 2 index. The meta-analysis included 7 206 774 SNPs and the results from both the Japonica and Illumina Arrays. The genome-wide significance level α was set to a p value of <5×10–8.

    Genetic correlation analysis

    Genetic correlation analysis using linkage disequilibrium score (LDSC) regression analysis15 was conducted to examine the potential genetic overlap between gout subtypes and between each subtype and serum uric acid (SUA) levels. For the regression, we used the 1000 genomes phase_3 East Asian LDSC and summary statistics for high-quality common SNPs present in the HapMap 3 reference panel for each analysis.

    Selection pressure analysis

    The details of genome-wide recent natural selection signature using singleton density score (SDS)16 from high-depth whole genome sequence data of the Japanese population had been described in the previous study.17 Using the same approach,17 we calculated the SDS of gout-risk variants identified in the present study, and evaluated overlaps between enrichment of the natural selection signatures and these variants from each subtype.

    Results

    Subtype GWASs of gout

    Figure 1 displays Manhattan plots of all and four clinical subtypes of gout, and figure 2 shows regional plots of novel loci. Compared with the plotted pattern for all clinically defined gout patients, each gout subtype (RUE type, ROL type, combined type and normal type) shows a subtype-specific plotted pattern and indicates the presence of cause-specific associated genes such as SLC2A9 and ABCG2. Table 2 lists the genome-wide significant loci from subtype GWASs. In total, 10, five, two, seven and three loci were identified at the genome-wide significant level to be associated with all, the RUE type, ROL type, combined type and normal type gout, respectively. Of these, two loci from all gout types, rs76499759 of PIBF1 and rs9926388 of ACSM2B, and another two intergenic loci from normal type gout, rs146978188 of CD2-PTGFRN and rs548944057 of SLC28A3-NTRK2, were detected as novel loci.

    Figure 1

    Manhattan plots of GWASs of subtypes of gout. Clinical subtypes and Manhattan plots of GWASs of all gout types, RUE type gout, combined type gout, ROL type gout and normal type gout are shown. The x-axis represents chromosomal positions and the y-axis shows −log10 p values. The dotted lines indicate the genome-wide significance threshold (p=5.0×10–8). FEUA, fractional excretion of uric acid (%); GWASs, genome-wide association studies; ROL, renal overload; RUE, renal underexcretion; UUE, urinary urate excretion (mg/h/1.73 m2). See table 1 and online supplementary figure S1 for a detailed classification of gout/hyperuricaemia.

    Figure 2

    Regional association plots of novel gout loci. Two loci were revealed to exceed the genome-wide significance level from the meta-analysis of GWASs from all gout patients, and another two loci from normal type gout patients. The highest association signal in each panel is located on (A) PIBF1, (B) ACSM2B, (C) CD2-PTGFRN and (D) SLC28A3-NTRK2. The region within 1 Mb from the single-nucleotide polymorphism (SNP) indicating the lowest p value is shown. (Upper panel) Plots of −log10 p values for the test of SNP association with gout. The SNP showing the lowest p value in the meta-analysis is depicted as a purple diamond. Other SNPs are colour-coded according to the extent of linkage disequilibrium (measured in r2) with the SNP showing the lowest p value. Recombination rates (centimorgans per Mb) estimated from HapMap Phase II data are also plotted. (Lower panel) RefSeq genes. Genomic coordinates are based on NCBI human genome reference sequence build hg19. The r2 data were calculated with 1000 Genomes Project Phase_3 JPT samples.45 GWASs, genome-wide association studies.

    Table 2

    Significant gout loci identified in the present genome-wide meta-analyses

    For all gout cases (table 2), 10 loci showed association at the genome-wide significance level: rs4148155 of ABCG2 (pmeta=1.81×10–101; ORs=2.23), rs671 of ALDH2 (pmeta=3.19×10–54; OR=1.93), rs3775946 of SLC2A9 (pmeta=3.73×10–41; OR=1.63), rs145954970 of SLC22A11 (pmeta=2.25×10–21; OR=12.43), rs3129500 of FAM35A (recently renamed as SHLD2, pmeta=4.34×10–17; OR=1.37), rs1260326 of GCKR (pmeta=2.07×10–13; OR=1.30), rs1010269 of BCAS3 (pmeta=1.81×10–9; OR=1.24), rs2817188 of SLC17A1 (pmeta=1.06×10–8; OR=1.35), rs9926388 of ACSM2B (pmeta=2.30×10–8; OR=1.24) and rs76499759 of PIBF1 (pmeta=2.79×10–8; OR=1.27). Among these 10 loci, PIBF1 and ACSM2B (table 2 and figure 2A,B). were identified for the first time as gout-risk loci at the genome-wide significance level. BCAS3 was identified here for the first time by the GWAS approach with Japanese individuals, while Li et al 18 reported that rs11653176, another SNP of BCAS3, is associated with gout based on a GWAS with a Han Chinese population. We also replicated its association with gout in a Japanese population using a candidate gene approach.19 Other loci have been previously reported to have an association with gout in our previous GWASs6 7 and association studies.20 21 As suggestive loci for gout, seven loci: PDZK1, TACR1-EVA1A, LOC100128993, ARID5B, TOLLIP-AS1-BRSK2, SLC38A1 and MLXIP, were detected (see online supplementary table S4). Of these, PDZK1, a gene encoding a scaffolding protein22 23 such as for urate transporters SLC22A12/URAT1 and ABCG2, was reported to have an association with SUA by a GWAS approach24 25 and with gout as a result of candidate gene approach studies.26–28

    As shown in online supplementary table S2, all the 3053 cases were classified into RUE type gout (654 cases), ROL type gout (486 cases), combined type gout (905 cases) and normal type gout (92 cases) for GWASs of gout subtypes.

    The meta-analysis of a GWAS of the RUE type gout (table 2) showed significant SNPs in the following five loci: rs3775948 of SLC2A9 (pmeta=8.01×10–22; OR=1.89), rs4148155 of ABCG2 (pmeta=2.54×10–16; OR=1.69), rs4646776 of ALDH2 (pmeta=5.80×10–16; OR=1.86), rs9420434 of GLUD1 (pmeta=8.62×10–11; OR=1.54) and rs76741582 of SLC22A11 (pmeta=1.06×10–10; OR=3.93). Since GLUD1 is in LD with SHLD2/FAM35A, for which we previously showed a significant association with gout,7 all of these five loci were previously identified as having an association with gout.6 7 We also detected nine suggestive loci: SMYD3, GCKR-C2orf16, SMARCC1, FRMD4B-MITF, ARL4A-ETV1, C7orf66-EIF3IP1, ASB10, PXDNL and LOC100287896-POLD3, for RUE type gout as shown in online supplementary table S4.

    From ROL type gout (table 2), rs4148155 of ABCG2 (pmeta=9.75×10–46; OR=2.79) and rs11066008 of ACAD10 (ALDH2) (pmeta=4.20×10–12; OR=1.79) were revealed to have an association. We had previously reported both to have a significant association with gout.6 7 20 21 Three suggestive loci, CNPY4, GRID1 and KCNJ2-CASC17, were also identified from ROL type gout (see online supplementary table S4).

    Combined type gout displayed the following seven significant loci (table 2): rs74904971 of ABCG2 (pmeta=1.53×10–56; OR=2.46), rs116873087 of NAA25 (pmeta=1.88×10–27; OR=2.38), rs3775948 of SLC2A9 (pmeta=1.43×10–19; OR=1.67), rs6586063 of SHLD2/FAM35A (pmeta=4.30×10–11; OR=1.49), rs9905274 of BCAS3 (pmeta=5.62×10–9; OR=1.37), rs11231879 of CDC42BPG (pmeta=7.66×10–10; OR=1.44) and rs1260326 of GCKR (pmeta=1.05×10–8; OR=1.37). There are studies on these which report an association between CDC42BPG and hyperuricaemia from the Japanese exome-wide association study,29 and GWAS on SUA from a Korean population.30 The significance around rs11231879 of CDC42BPG was, however, no longer evident when conditioned on rs11231879 itself, nor when conditioned on the secondarily significant SNP, rs56093838 of SLC22A12/URAT1, demonstrating that these signals were from the same locus (see online supplementary figure S2). Because URAT1 is a well-known urate transporter that markedly affects SUA level, these findings indicate that the true associated gene for combined type gout on chromosome 11q13.1 locus is not CDC42BPG, but SLC22A12/URAT1. rs116873087 of NAA25 is in strong LD with rs671 of ALDH2 (r2=0.97 in the 1000 Genomes Project Phase_3: JPT samples). All of the seven loci had therefore been previously identified as having an association with gout.6 7 20 21 Two suggestive loci, NCKAP5-MIR3679 and PRDM8-FGF5, were also identified from combined type gout (see online supplementary table S4).

    Three significant loci were found from normal type gout: rs548944057 of SLC28A3-NTRK2 (pmeta=2.91×10–8; OR=4.04), rs4148155 of ABCG2 (pmeta=3.64×10–8; OR=2.34) and rs146978188 of CD2-PTGFRN (pmeta=4.93×10–8; OR=6.31). Of these, two intergenic loci are novel susceptibility loci for gout (table 2 and figure 2C,D) There were eight suggestive loci: ZNF639-MFN1, RUNX2-CLIC5, DST, HGF, MED27-NTNG2, LINC00944-LINC02372, SV2B and GABPA, for normal type gout as shown in online supplementary table S4.

    The LDSC regression analysis was performed to examine the potential genetic overlap between gout subtypes and between each subtype and SUA levels. Significant positive genetic correlations were observed among these subtypes as well as between these traits and SUA levels (see online supplementary figure S3).

    Selection pressure analysis of gout susceptibility

    We also performed selection pressure analysis of gout on the basis of a previous report on the recent natural selection signature in the Japanese population.17 This analysis enables us to elucidate the genetic risks of gout characterised in the recent evolutionary history (2000–3000 years) of the Japanese population. Because ALDH2 was reported to be subjected to strong selection pressure in the Japanese population,17 and because ABCG2,31 32 as well as ALDH2,20 21 33 is a well-known strong susceptible gene for gout in Japanese, selection pressure analysis was initially performed outside of these two loci. As a result, only combined type gout showed significant enrichment of selection pressure (p=0.026; figure 3 and online supplementary table S5). When the ABCG2 locus was included in the analysis, all of these subtypes except for normal type gout then showed significant enrichment of selection pressure. As expected, analysis including ALDH2 and outside of ABCG2 showed significant enrichment of selection pressure except for normal type gout. This trend also persisted in the analysis with all associated SNPs (figure 3 and online supplementary table S5).

    Figure 3

    Overlap between natural selection signatures and genetic risk of gout and its subtypes in the Japanese population. For each trait, inflation of the selection χ2 value is indicated along the x-axis, and −log10(p) of enrichment is plotted along the y-axis. The horizontal grey line represents significance threshold (p<0.05). Because ABCG2 and ALDH2 are associated with a well-known genetic risk of gout in Japanese individuals, selection pressure analyses without these loci were performed initially (filled circle), and subsequent analyses were conducted with ABCG2 (filled triangle), as well as ABCG2 and ALDH2 (filled square). When calculated with ABCG2 (outside of ALDH2) (filled triangle), all but normal type gout showed significant selective pressure, indicating that ABCG2 is involved in adaptive evolution in Japanese for having higher SUA levels, which can result in gout. Finally, all but normal type gout also showed significant selective pressure with ABCG2 and ALDH2 loci (filled square). ROL, renal overload; RUE, renal underexcretion; SDS, singleton density score; SNPs, single-nucleotide polymorphisms.

    Discussion

    We previously performed GWASs of clinically defined gout cases in the Japanese population and found loci including ABCG2, SLC2A9, ALDH2 (CUX2), GCKR and SHLD2/FAM35A to be associated with gout at a genome-wide significant level.6 7 While our previous study was performed with broader subtypes, it is one of the unique points of this study that the present GWAS is the first to be conducted with four differentiated subtypes: the RUE type, ROL type, combined type and normal type gout (figure 1), which are commonly used in daily clinical settings. Two platforms for GWASs and meta-analyses between them were used to perform a comprehensive genetic association search.

    The results allowed four novel loci to be identified from the present study. The pathophysiological associations of the two novel loci from all gout types, PIBF1 and ACSM2B with gout, are totally unknown. PIBF (progesterone-induced blocking factor) is induced by progesterone and is a mediator that exerts substantial antiabortive activities, including cytokine secretion.34 PIBF1 might be therefore involved in decreased urate production by female hormone and/or decreased inflammatory response to gout attack. ACSM2B (acyl-CoA synthetase medium chain family member 2B) is a predominant transcript in the human liver and an enzyme catalysing the activation of medium-chain length fatty acids.35 Because ASCM2B is involved in the production of ATP, a purine body metabolised to urate, ASCM2B might contribute to gout via that mechanism. While the present study showed significance at SNPs of two genes, it is of course possible that these are mere markers and that the true risk SNPs are present close by. For example, since UMOD, a causative gene of uromodulin-associated kidney disease (previously known as familial juvenile hyperuricemic nephropathy)36 is located 180 kb downstream from rs9926388 of ACSM2B, there might be a relationship between them.

    Another two intergenic loci from normal type gout, rs146978188 of CD2-PTGFRN and rs548944057 of SLC28A3-NTRK2, were detected for the first time to have an association with gout. CD2 is well known as a surface antigen found on all peripheral blood T-cells, and PTGFRN/CD9P-1 encodes prostaglandin F2 receptor inhibitor. Neither of these was previously known to have an association with gout or uric acid. CD101, which is next to PTGFRN (150 kb downstream from rs146978188), is reported to be expressed on macrophages/monocytes and T-cells, to confer a modulatory/coregulatory function, and to be conspicuously downregulated in rheumatoid arthritis patients.37 Because macrophages are a chief contributor to gouty attack, CD101 might be the true susceptible gene for normal type gout. rs146978188 is also in LD with an SNP of SLC22A15/FLIPT1, an orphan transporter gene in the same family as SLC22A12/URAT1, a well-known urate transporter gene that is strongly associated with gout. This transporter gene might have an association with normal type gout. SLC28A3/CNT3 is reported to be an Na+-dependent pyrimidine-selective and purine-selective transporter found predominantly in the intestine and kidney,38 39 which are the main urate excretion pathways. Further analysis is needed to elucidate the relationship between this transporter gene and gout, including normal type gout.

    While the present study revealed only suggestive loci from other subtype GWASs, some of these loci from subtype gout also suggest a relationship with RUE, extra-RUE and/or overproduction of urate, because gout subtypes reflect their clinical parameters (table 1). GWASs with these subtypes should therefore be useful for estimating the expression and function of proteins encoded by identified loci. rs557868370 of SLC38A1, a transporter gene, was detected as a suggestive locus (see online supplementary table S4). Since there is a study reporting the relationship between oxidative stress and SLC38A1/SNAT1,40 it might have a relationship with urate, which also has antioxidative stress effects. Because genetic variants in urate transporter genes such as ABCG2, SLC2A9 and SLC22A12 are well known to cause SUA variation and gout, it is also possible that the SLC38A1/SNAT1 transporter is involved in urate or purine transport.

    Very recently, Tin et al 41 performed transancestry GWAS meta-analysis of SUA and gout, including self-reported gout cases. We compared these SNPs and found five of the 10 loci detected here (GCKR, SLC2A9, ABCG2, SLC17A1 and BCAS3) to be associated with gout (see online supplementary table S6), indicating population differences in the genetic basis of gout.

    Taking into account the evidence of shared genetic background among gout subtypes (see online supplementary figure S3) and the presence of subtype-specific genetic factors of gout (figure 1, table 2), these results will provide helpful information for the development of novel subtype-specific genome tailor-made medicines and/or prevention for gout and hyperuricaemia.

    Adaptive evolution results from adaptation to environmental changes over generations. Selection pressure analysis using SDS in the present study has elucidated which genes have been involved in adaptive evolution over the last 2000–3000 years in the Japanese population. The results revealed that the Japanese population has evolved to have higher SUA levels, which can result in gout, due to the ABCG2 locus in addition to the already-known ALDH2 gene. ABCG2 is now a well-known susceptible gene for hyperuricaemia and gout, especially in the Japanese population.31 42 Few patients had gout before Westernisation of Japan about 150 years ago, which brought more purine-rich foods to Japan. The fact that ABCG2 also has a relationship with SUA levels might thus have caused few problems to Japanese people until recently. The upside of SUA elevation in the Japanese population might include resistance to antioxidative stress, lower cancer risk, neuroprotective effect and longevity.43 44 Selection pressure analysis with other populations will also generate more information on evolutionary association with gout to elucidate this hypothesis.

    In summary, we performed GWASs of all gout as well as of distinct gout subtypes, and identified multiple subtype-specific loci including four novel loci. Selection pressure analysis revealed significant enrichment of selection for the ABCG2 and ALDH2 loci in Japanese gout patients of each subtype. These findings will lead to elucidation of the molecular pathophysiology of each gout/hyperuricaemia subtype and the development of novel subtype-specific genome tailor-made medicine/prevention of gout and hyperuricaemia.

    Acknowledgments

    We would like to thank all the participants for their generous involvement in this study. We are grateful to the staff of the institutions participating in the J-MICC Study for their outstanding assistance in collecting samples and clinical information. We are indebted to M Miyazawa, K Morichika, M Sakiyama, M Nakajima, K Maehara and M Kirihara (National Defense Medical College) for genetic analysis and to N Hamajima at Nagoya University Graduate School of Medicine for sample collection. We are also indebted to H Fujiwara (Midorigaoka Hospital), M Senda and M Hamamura (Ryougoku East Gate Clinic).

    References

    Footnotes

    • Handling editor Josef S Smolen

    • AN, MN, YK, KY and HN contributed equally.

    • Collaborators Members of the Japan Gout Genomics Consortium (Japan Gout) are: Katsuhisa Inoue (Department of Biopharmaceutics, School of Pharmacy, Tokyo University of Pharmacy and Life Sciences, Tokyo), Tomoya Yasujima, Hiroaki Yuasa (Department of Biopharmaceutics, Graduate School of Pharmaceutical Sciences, Nagoya City University, Nagoya, Aichi), Hiroaki Ikezaki, Masayuki Murata (Department of Environmental Medicine and Infectious Disease, Graduate School of Medical Sciences, Kyushu University), Keito Morimoto (Department of Pharmacy, the University of Tokyo Hospital, Tokyo), Mitsuhiro Yokota (Department of Medical Biochemistry, Kurume University School of Medicine, Kurume, Fukuoka), Sahoko Ichihara (Department of Environmental and Preventive Medicine, Jichi Medical University School of Medicine, Shimotsuke, Tochigi), Tatsuaki Matsubara (Department of Internal Medicine, School of Dentistry, Aichi Gakuin University, Nagoya, Aichi), Toshimitsu Ito (Department of Internal Medicine, Self-Defense Forces Central Hospital, Tokyo), Miki Ueno, Kimiko Hayano (Division of Nursing, National Defense Medical College, Tokorozawa, Saitama), Kunio Mizutari, Akihiro Shiotani (Department of Otolaryngology-Head and Neck Surgery, National Defense Medical College, Tokorozawa, Saitama), Yuka Miyoshi, Satoko Suzuki, Satoko Iwasawa (Department of Preventive Medicine and Public Health, National Defense Medical College, Tokorozawa, Saitama), Yuka Aoyagi, Yuka Aoki (Department of Integrative Physiology and Bio-Nano Medicine, National Defense Medical College, Tokorozawa, Saitama), Tatsuya Saitoh (Laboratory of Bioresponse Regulation, Graduate School of Pharmaceutical Sciences, Osaka University, Suita, Osaka).

    • Contributors AN, MNakatochi and HMatsuo conceived and designed this study. KY, HNakaoka, II, YO and NS assisted with research design. TS, KO, HO, MNagase, YH, MKawaguchi, MTakao, RS, THosoya, KI and HMatsuo collected and analysed clinical data of cases. TK, AHishida, KK, MW, DM, TTamura, TNishiyama, CS, SKK, NT, NK, HMikami, TTakezaki, KM, SSuzuki and KW collected and analysed clinical data of controls. AN, YKawamura, KY, SShimizu, THigashino, YTakada, ID, YKamatani, MKubo, NS and HMatsuo performed genetic analysis. AN, MNakatochi, YKawamura, HNakaoka, SShimizu, YS, TNakamura, HNakashima and YO performed statistical analysis. HMatsuo organized this collaborative study. KY, HNakaoka, THigashino, MTsunoda, ID, AHozawa, KH, YToyoda, YKubota, TTakada, HS, BS, TM, TR, THosoya, YKamatani, MKubo, KI, KW, II, YO and NS provided intellectual input and assisted with the preparation of the manuscript. AN, MNakatochi, YKawamura, YO, NS and HMatsuo wrote the manuscript.

    • Funding This study was supported by grants from the Ministry of Education, Culture, Sports, Science and Technology (MEXT) of Japan, including JSPS Kakenhi Grants (Nos. 16H06279, 16H01808, 17H04128, 18KK0247, 19K22786 and 25293145) and a Grant-in-Aid for Scientific Research on Innovative Areas (No. 221S0002), the Ministry of Defense, and the Kawano Masanori Memorial Foundation for the Promotion of Pediatrics and the Gout Research Foundation of Japan. The J-MICC Study was supported by Grants-in-Aid for Scientific Research from MEXT, including those for Priority Areas (No. 17015018) and Innovative Areas (No. 221S0001), as well as by a JSPS Kakenhi grant (No. 16H06277). The BioBank Japan Project was supported by grants from the Japan Agency for Medical Research and Development since April 2015, and the Ministry of Education, Culture, Sports, Science and Technology from April 2003 to March 2015.

    • Competing interests TTakada, NS and HMatsuo have a patent pending based on the work reported in this paper.

    • Patient and public involvement Patients and/or the public were not involved in the design, or conduct, or reporting or dissemination plans of this research.

    • Patient consent for publication Not required.

    • Ethics approval This study was approved by the institutional ethical committees and written consent was obtained from all of its participants. All procedures involved were performed in accordance with the Declaration of Helsinki.

    • Provenance and peer review Not commissioned; externally peer reviewed.

    • Data availability statement No data are available.