Identification of four genes as novel susceptibility loci for early‑onset type 2 diabetes mellitus, metabolic syndrome, or hyperuricemia
- Authors:
- Published online on: May 29, 2018 https://doi.org/10.3892/br.2018.1105
- Pages: 21-36
-
Copyright: © Yamada et al. This is an open access article distributed under the terms of Creative Commons Attribution License.
Abstract
Introduction
Diabetes mellitus (DM) has reached epidemic proportions and affects more than 382 million individuals worldwide (1): The number of patients with DM is expected to increase beyond 592 million individuals by 2035 (1). Approximately 90 to 95% of individuals with DM have type 2 DM (T2DM), the characteristics of which can range from predominant insulin deficiency with relatively minor insulin resistance to predominant insulin resistance with relatively minor insulin deficiency (2). T2DM is a major cause of nephropathy, retinopathy, and neuropathy as well as cardiovascular disease and stroke (3,4). Although obesity resulting from a sedentary lifestyle and overeating is an important risk factor for T2DM, genetic components are involved in the pathogenesis of this condition, given that a positive family history confers a 2.4-fold increased risk for T2DM (5). The heritability of T2DM has been estimated to be 50 to 60% (6).
Genome-wide association studies (GWASs) and meta-analyses have identified >120 susceptibility loci for T2DM (7) in individuals of European (8–13) or African (14) ancestry, in East Asians (15), or in multiple ethnic groups (16,17). Among Japanese, GWASs have identified KCNQ1, UBE2E2, and C2CD4A-B (18–20) as susceptibility genes for T2DM, and a recent meta-analysis identified an additional seven susceptibility loci for this condition (21). Genetic variants that contribute to predisposition to T2DM in Japanese subjects, however, remain to be identified definitively.
Metabolic syndrome (MetS) is a cluster of metabolic traits including abdominal obesity, an increased serum triglycerides, a decreased serum high-density lipoprotein (HDL) cholesterol, high blood pressure (BP), and an increased fasting plasma glucose (FPG) level (22). MetS is a risk factor for atherosclerotic cardiovascular disease, DM (23), and cancer (24). The etiology of MetS is highly complex, with both genetic and environmental factors being thought to play important roles. The heritability of MetS has been estimated to be approximately 50% (25) and traits of this syndrome to be 28 to 48% (26).
GWASs have suggested various loci or genes involved in predisposition to MetS or to traits of this syndrome in individuals of European (27,28) or African (29) ancestry or in Asian Indians (30) or Chinese individuals (31). However, genetic variants that contribute to predisposition to MetS in Japanese individuals remain to be identified definitively.
Circulating uric acid levels are regulated by multiple renal transporters by mediating the excretion or reabsorption of uric acid in the proximal kidney tubules (32). Hyperuricemia is an important risk factor for gout, a common inflammatory arthritis (33), as well as for cardiovascular disease (34) and cancer (35). The heritability of the serum concentration of uric acid has been estimated to be 40% (36), suggesting that genetic variants contribute to regulation of the serum uric acid level by influencing uric acid synthesis, excretion, or reabsorption (36,37).
GWASs have identified single nucleotide polymorphisms (SNPs) significantly associated with the serum uric acid concentration or the prevalence of gout (38–45). A large-scale GWAS in European ancestry populations identified 28 loci that influence the serum concentration of uric acid (46). Although several SNPs have been shown to be associated with gout in Japanese (47,48), genetic variants that contribute to predisposition to hyperuricemia in Japanese remain to be identified definitively.
In a family study of T2DM, a heritability of T2DM was higher in early-onset than late-onset individuals (6). These observations indicate that early-onset T2DM has a strong genetic component (6,49). Similar to T2DM, early-onset forms of MetS (50,51), hyperuricemia, and gout (52,53) have been shown to have strong genetic components. Given that genetic contribution may be greater in early-onset forms of T2DM, MetS, and hyperuricemia than in late-onset forms, statistical power of the genetic association study may be increased by focusing on early-onset subjects with these conditions.
In the present study, we performed exome-wide association studies (EWASs) with the use of human exome array-based genotyping methods to identify genetic variants that confer susceptibility to early-onset T2DM, MetS, or hyperuricemia in Japanese patients. To increase the statistical power of EWASs, early-onset subjects were examined.
Materials and methods
Study subjects
In previous studies, the median age of subjects with T2DM (54), MetS (55), or hyperuricemia (56) was 68, 64, or 62 years, respectively. We thus defined subjects aged ≤65 years as early-onset cases in the present study. A total of 8,102 individuals aged ≤65 years were examined. The subjects were recruited from individuals either who visited outpatient clinics of or were admitted to participating hospitals in Japan (Gifu Prefectural Tajimi Hospital, Tajimi; Gifu Prefectural General Medical Center, Gifu; Japanese Red Cross Nagoya First Hospital, Nagoya; Northern Mie Medical Center Inabe General Hospital, Inabe; Hirosaki University Hospital and Hirosaki Stroke and Rehabilitation Center, Hirosaki) because of various symptoms or for an annual health checkup between October 2002 and March 2014; or who were community-dwelling individuals recruited to a population-based cohort study in Inabe between March 2010 and september 2014 (57).
T2DM was defined according to the criteria of the World Health Organization as described previously (2,58,59). Subjects with T2DM had an FPG level of ≥6.93 mmol/l (126 mg/dl) or a blood hemoglobin A1c content of ≥6.5% or were taking antidiabetes medication. Individuals with T1DM, maturity-onset diabetes of the young, DM associated with mitochondrial diseases or single-gene disorders, pancreatic diseases, or other metabolic or endocrinologic diseases were excluded from the study. Those taking medications that may cause secondary DM were also excluded. The control subjects had an FPG level of <6.05 mmol/l (110 mg/dl), a blood hemoglobin A1c content of <6.2%, and no history of DM or of having taken antidiabetes medication. We thus examined 1,696 subjects with T2DM and 5,711 controls.
Diagnosis of MetS was based on a modified version of the definition proposed by the International Diabetes Federation Task Force on Epidemiology and Prevention; National Heart, Lung, and Blood Institute; American Heart Association; World Heart Federation; International Atherosclerosis Society; and International Association for the Study of Obesity (22). We used cut-off values for waist circumference of ≥90 cm in men or ≥80 cm in women on the basis of a recommendation of the International Diabetes Association (22). A total of 2,296 subjects with MetS thus had three or more of the following five components: i) A waist circumference of ≥90 cm for men or ≥80 cm for women; ii) a serum triglyceride concentration of ≥1.65 mmol/l (150 mg/dl) or drug treatment for elevated triglycerides; iii) a serum HDL-cholesterol concentration of <1.04 mmol/l (40 mg/dl) for men or <1.30 mmol/l (50 mg/dl) for women; iv) a systolic BP of ≥130 mmHg, diastolic BP of ≥85 mmHg, or drug treatment for hypertension; and v) an FPG level of ≥5.50 mmol/l (100 mg/dl) or drug treatment for elevated glucose. History of obesity, dyslipidemia, hypertension, or DM was evaluated with a detailed questionnaire. The control subjects comprised 1,919 individuals who had none of the five traits of MetS.
Hyperuricemia was defined as a serum uric acid concentration of >416 µmol/l (7 mg/dl) or taking of uric acid-lowering medication. Individuals taking drugs that potentially caused secondary hyperuricemia were excluded. The control individuals for the study of hyperuricemia had a serum uric acid concentration of ≤416 µmol/l and had no history of hyperuricemia or gout or of taking uric acid-lowering medication. Thus, we examined 1,365 subjects with hyperuricemia and 6,554 controls.
The 1,040 subjects with both T2DM and MetS as well as 1,884 controls overlapped between the corresponding studies, as did the 375 subjects with both T2DM and hyperuricemia and 4,809 controls as well as the 628 subjects with both MetS and hyperuricemia and 1,771 controls.
EWASs
Venous blood (5–7 ml) was collected into tubes containing 50 mmol/l ethylenediaminetetraacetic acid (disodium salt), peripheral blood leukocytes were isolated, and genomic DNA was extracted from these cells either with the use of a kit (Genomix from Talent Srl, Trieste, Italy or SMITEST EX-R&D from Medical & Biological Laboratories, Co., Ltd., Nagoya, Japan).
The EWAS for T2DM, MetS, or hyperuricemia included 7,407 individuals (1,696 subjects with T2DM, 5,711 controls), 4,215 individuals (2,296 subjects with MetS, 1,919 controls), or 7,919 individuals (1,365 subjects with hyperuricemia, 6,554 controls), respectively. The EWASs were performed with the use of a Human Exome-12 v1.2 DNA Analysis BeadChip or Infinium Exome-24 v1.0 BeadChip (Illumina, San Diego, CA, USA). These exome arrays include putative functional exonic variants selected from ~12,000 individual exome and whole-genome sequences. The exonic content consisted of ~244,000 SNPs from diverse populations, including European, African, Chinese, and Hispanic individuals (60). SNPs contained in only one exome array (~2.6% of all SNPs) were excluded from analysis. We performed quality control (61) as follows: i) Genotyping data with a call rate of <97% were discarded, with the mean call rate for the remaining data being 99.9%. ii) Sex specification was checked for all the samples, and those for which sex phenotype in the clinical records was inconsistent with genetic sex were discarded. iii) Duplicated samples and cryptic relatedness were checked by calculation of identity by descent, with all the pairs of samples showing a value of >0.1875 being inspected and one sample from each pair excluded. iv) Heterozygosity of SNPs was calculated for all the samples, and those with extremely low or high heterozygosity (>3 standard deviations from the mean) were discarded. v) SNPs in sex chromosomes or mitochondrial DNA were excluded from the analysis, as were non-polymorphic SNPs or SNPs with a minor allele frequency of <1.0%. vi) SNPs whose genotype distributions deviated significantly (P<0.01) from the Hardy-Weinberg equilibrium in control individuals were discarded. vii) Genotype data were examined for population stratification by principal components analysis (62), and population outliers were excluded. Totals of 31,210, 31,521, or 31,142 SNPs passed quality control for the T2DM, MetS, and hyperuricemia studies, respectively, and were subjected to analysis.
Statistical analysis
For analysis of characteristics of the study subjects, quantitative or categorical data were compared between individuals with T2DM, MetS, or hyperuricemia and corresponding controls with the unpaired Student's t-test or Pearson's Chi-square test, respectively. Allele frequencies were estimated by the gene counting method, and Fisher's exact test was used to identify departure from Hardy-Weinberg equilibrium. The relationship of allele frequencies of SNPs to T2DM, MetS, or hyperuricemia in the EWAS was examined using the Fisher's exact test. To compensate for multiple comparisons of allele frequencies with T2DM, MetS, or hyperuricemia, we applied Bonferroni's correction for statistical significance of association. The significance level was set at P<1.60×10−6 (0.05/31210), P<1.59×10−6 (0.05/31521), or P<1.61×10−6 (0.05/31142) for the EWAS of T2DM, MetS, or hyperuricemia, respectively. The inflation factor (λ) was 1.04 for T2DM, 1.05 for MetS, or 1.09 for hyperuricemia. Multivariable logistic regression analysis was performed with T2DM, MetS, hyperuricemia as a dependent variable and independent variables including age, sex (0, woman; 1, man), and genotype of each SNP. Genotypes of each SNP were assessed according to dominant [0, AA; 1, AB + BB (A, major allele; B, minor allele)], recessive (0, AA + AB; 1, BB), and additive genetic models, and the P-value, OR, and 95% confidence interval were calculated. Additive models comprised additive 1 (0, AA; 1, AB; 0, BB) and additive 2 (0, AA; 0, AB; 1, BB) scenarios, which were analyzed simultaneously with a single statistical model. The relation of genotypes of identified SNPs to T2DM-, MetS-, or hyperuricemia-related traits was examined by one-way analysis of variance (ANOVA). Bonferroni's correction was also applied to other statistical analysis as indicated. Statistical tests were performed with JMP Genomics version 9.0 software (SAS Institute, Cary, NC, USA).
Results
Characteristics of subjects
The characteristics of the 7,407 subjects enrolled in the T2DM study are shown in Table I. Age, the frequency of men, and the prevalence of obesity, hypertension, dyslipidemia, chronic kidney disease (CKD), and hyperuricemia as well as body mass index (BMI), systolic and diastolic BP, and serum concentrations of triglycerides, creatinine, and uric acid were greater, whereas serum concentration of HDL-cholesterol and estimated glomerular filtration rate (eGFR) were lower, in subjects with T2DM than in controls.
Characteristics of the 4,215 subjects enrolled in the MetS study are shown in Table II. Age, the frequency of men, and the prevalence of smoking, CKD, and hyperuricemia as well as BMI, blood hemoglobin A1c content, and serum concentrations of low-density lipoprotein (LDL) cholesterol, creatinine, and uric acid were greater, whereas eGFR were lower, in subjects with MetS than in controls.
Characteristics of the 7,919 subjects enrolled in the hyperuricemia study are shown in Table III. Age, the frequency of men, and the prevalence of smoking, obesity, hypertension, DM, dyslipidemia, and CKD as well as BMI, systolic and diastolic BP, FPG level, and serum concentrations of triglycerides and creatinine were greater, whereas the serum concentration of HDL-cholesterol and eGFR were lower, in subjects with hyperuricemia than in controls.
EWAS for T2DM, MetS, or hyperuricemia
We examined the relationship of allele frequencies of 31,210 SNPs that passed quality control to T2DM with the use of Fisher's exact test. After Bonferroni's correction, four SNPs were significantly (P<1.60×10−6) associated with T2DM (Table IV). The relationship of allele frequencies of 31,521 SNPs to MetS was examined with Fisher's exact test. After Bonferroni's correction, six SNPs were significantly (P<1.59×10−6) associated with MetS (Table V). The relationship of allele frequencies of 31,142 SNPs to hyperuricemia was also examined with Fisher's exact test. After Bonferroni's correction, nine SNPs were significantly (P<1.61×10−6) associated with hyperuricemia (Table VI).
 | Table IV.The four SNPs significantly (P<1.60×10−6) associated with type 2 diabetes mellitus in the exome-wide association study. |
| Table V.The six SNPs significantly (P<1.59×10−6) associated with metabolic syndrome in the exome-wide association study. |
| Table VI.The nine SNPs significantly (P<1.61×10−6) associated with hyperuricemia in the exome-wide association study. |
Multivariable logistic regression analysis of the relationship of SNPs to T2DM, MetS, or hyperuricemia
The relationship of the four identified SNPs in the EWAS of T2DM to this condition was further examined by multivariable logistic regression analysis with adjustment for age and sex (Table VII). Three SNPs (rs141569282 of OR4F6, rs140232911 of ZNF860, rs150552771 of LPGAT1) were significantly [P<0.0031 (0.05/16) in at least one genetic model] related to T2DM. The relationship of the six SNPs identified in the EWAS for MetS to this condition was examined by multivariable logistic regression analysis with adjustment for age and sex (Table VIII). All the SNPs were significantly [P<0.0021 (0.05/24)] related to MetS. The relationship of the 9 SNPs identified by the EWAS of hyperuricemia to this condition was also examined by multivariable logistic regression analysis with adjustment for age and sex (Table IX). All SNPs were significantly [P<0.0014 (0.05/36)] related to hyperuricemia.
| Table VII.Association of SNPs to type 2 diabetes mellitus as determined by multivariable logistic regression analysis. |
| Table VIII.Relation of SNPs to metabolic syndrome as determined by multivariable logistic regression analysis. |
| Table IX.Association of SNPs to hyperuricemia as determined by multivariable logistic regression analysis. |
Relationship of SNPs associated with T2DM to FPG level or blood hemoglobin A1c content
We examined the relationship of genotypes for the 3 SNPs associated with T2DM to FPG level or blood hemoglobin A1c content by one-way ANOVA (Table X). A SNP rs141569282 of OR4F6 was significantly [P<0.0083 (0.05/6)] associated with FPG level and blood hemoglobin A1c content, and rs140232911 of ZNF860 with blood hemoglobin A1c content.
| Table X.Relation of SNPs identified in the present study to fasting plasma glucose level and blood hemoglobin A1c content. |
Relationship of SNPs associated with MetS to each trait of MetS
We examined the relationship of genotypes for the six SNPs associated with MetS to waist circumference, serum concentrations of triglycerides and HDL-cholesterol, systolic and diastolic BP, and FPG level by one-way ANOVA (Table XI). Three SNPs (rs1053266 of CCDC6, rs141569282 of OR4F6, rs200716952 of HLA-DQB2) were significantly (P<0.0014) related to all the traits; rs2075291 of APOA5 to waist circumference and serum concentrations of triglycerides and HDL-cholesterol; rs11823828 of OR52E4 to serum concentrations of triglycerides and HDL-cholesterol, systolic BP, and FPG level; and rs17482753 at chromosome 8p21.3 to serum concentrations of triglycerides and HDL-cholesterol.
Relationship of SNPs associated with hyperuricemia to serum concentrations of uric acid
We examined the relationship of the 9 SNPs associated with hyperuricemia to serum concentrations of uric acid by one-way ANOVA (Table XII). Seven SNPs (rs121907892 of SLC22A12, rs3782886 of BRAP, rs11066015 of ACAD10, rs11066280 of HECTD4, rs12229654 at chromosome 12q24.1, rs3775948 of SLC2A9, rs2305335 of HERPUD2) were significantly (P<0.0056) associated with serum concentrations of uric acid.
| Table XII.Relationship of SNPs identified in the present study to the serum concentration of uric acid. |
Linkage disequilibrium (LD) analyses
We examined LD in SNPs associated with MetS or hyperuricemia. For the MetS study, rs11823828 of OR52E4 and rs2075291 of APOA5 were not in LD [square of the correlation coefficient (r2)<0.001]. For the hyperuricemia study, 6 SNPs were located at chromosomal 12q24.11 to 12q24.13. LD plots of these SNPs are shown in Fig. 1. Strong LD was observed between rs10849915 and rs10774610 of CCDC63 (r2, 0.991) and among rs3782886 of BRAP, rs11066015 of ACAD10, and rs11066280 of HECTD4 (r2, 0.876 to 0.941).
Relation of genes, chromosomal loci, and SNPs identified in the present study to phenotypes reported by previous GWASs. We examined the genes, chromosomal loci, and SNPs identified in the present study to phenotypes previously reported by GWASs available in Genome-Wide Repository of Associations Between SNPs and Phenotypes (GRASP) Search database (https://grasp.nhlbi.nih.gov/Search.aspx) developed by the Information Technology and Applications Center, National Center for Biotechnology Information, National Heart, Lung, and Blood Institute, National Institute of Health (Bethesda, MD, USA).
In the T2DM study, none of the three genes or SNPs was shown to be associated with T2DM or diabetes-related traits in previous GWASs (Table XIII). In the MetS study, HLA-DQB2 was shown to be related to T1DM and plasma total cholesterol; rs17482753 at chromosome 8p21.3 to plasma concentrations of triglycerides and HDL-cholesterol; CCDC6 to serum uric acid level; and APOA5 to plasma concentrations of triglycerides, HDL-cholesterol, and LDL-cholesterol (Table XIV). In the hyperuricemia study, SLC2A9, SLC22A12, and BRAP were shown to be related to serum uric acid concentrations; rs12229654 at 12q24.1 to plasma HDL-cholesterol; ACAD10 to plasma LDL-cholesterol and T1DM; and HECTD4 to plasma HDL-cholesterol and LDL-cholesterol and FPG level (Table XV).
| Table XIII.Relationship of genes and SNPs associated with type 2 diabetes mellitus in the present study to previously reported diabetes-related phenotypes. |
| Table XIV.Relationship of genes, chromosomal locus, and SNPs associated with metabolic syndrome in the present study to previously reported metabolic disease-related phenotypes. |
| Table XV.Relationship of genes, chromosomal locus, and SNPs associated with hyperuricemia in the present study to previously reported metabolic disease-related phenotypes. |
Discussion
T2DM, MetS, and hyperuricemia are important public health problems because of the high prevalence of these conditions as well as risk factors for more serious conditions such as cardiovascular disease, cancer, or gout (1,3,4,23,24,33–35). Identification of genetic variants that confer susceptibility to T2DM, MetS, and hyperuricemia are thus clinically important to prevent these conditions. We have now performed EWASs for T2DM, MetS, and hyperuricemia in early-onset subjects with these conditions who likely had greater genetic components compared with late-onset individuals.
In the T2DM study, rs150552771 of LPGAT1, rs140232911 of ZNF860, and rs141569282 of OR4F6 were significantly associated with early-onset T2DM. None of these genes was shown to be associated with T2DM or diabetes-related traits in the previous GWASs. Given that rs150552771 of LPGAT1 was not related to FPG level or blood hemoglobin A1c content, LPGAT1 was removed from new susceptibility locus, even though this discrepancy may be attributable to the effect of medical treatment for T2DM. We have thus newly identified ZNF860 and OR4F6 as susceptibility loci for T2DM. A SNP rs141569282 of OR4F6 was significantly related to FPG level and blood hemoglobin A1c content with the minor A allele being related to decreases in these parameters, while rs140232911 of ZNF860 was related to blood hemoglobin A1c content with the minor T allele being related to an increase in this parameter. Analyses of these traits as well as logistic regression analysis indicate that the A allele of rs141569282 in OR4F6 is protective against T2DM, whereas the T allele of rs140232911 in ZNF860 is a risk factor for this condition.
In the MetS study, six SNPs in five genes and one chromosomal locus were significantly associated with early-onset MetS. Of these genes and locus, HLA-DQB2 (63), rs17482753 at 8p21.3 (63), and APOA5 (63) were previously shown to be related to lipid profiles; and CCDC6 to serum uric acid level (46). OR52E4 or OR4F6 has not been shown to be associated with MetS or MetS-related traits in the previous GWASs. A SNP rs11823828 of OR52E4 was significantly related to serum concentrations of triglycerides and HDL-cholesterol, systolic BP, and FPG level; and rs141569282 of OR4F6 to all traits. We have thus newly identified OR52E4 and OR4F6 as susceptibility loci for MetS. The minor G allele of rs11823828 in OR52E4 was significantly related to increased serum triglycerides, decreased serum HDL-cholesterol, increased systolic BP, and increased FPG level, whereas the minor A allele of rs141569282 in OR4F6 was related to reduced waist circumference, decreased serum triglycerides, increased serum HDL-cholesterol, decreased systolic and diastolic BP, and reduced FPG level. Analyses of these traits and logistic regression analysis indicate that the G allele of rs11823828 in OR52E4 represents a risk factor for MetS, whereas the A allele of rs141569282 in OR4F6 is protective against this condition.
In the hyperuricemia study, 9 SNPs in seven genes and one chromosomal locus were significantly associated with early-onset hyperuricemia. Of these genes and locus, SLC2A9 (46), SLC22A12 (46), and BRAP (46) were previously shown to be related to serum uric acid level; and rs12229654 at 12q24.1 (64), ACAD10 (63), and HECTD4 (63) to lipid profiles. HERPUD2 or CCDC63 has not been shown to be associated with hyperuricemia, gout, or serum uric acid level in the previous GWASs. A SNP rs2305335 of HERPUD2 was significantly related to the serum uric acid level, whereas two SNPs of CCDC63 were not related to this parameter. Therefore CCDC63 was removed from new loci. We have thus newly identified HERPUD2 as a susceptibility locus for hyperuricemia. The minor A allele of rs2305335 in HERPUD2 was significantly related to increased serum uric acid. Examination of this trait and logistic regression analysis indicate that the A allele of rs2305335 in HERPUD2 represents a risk factor or hyperuricemia. Furthermore, OR4F6 was significantly associated with both T2DM and MetS. We thus newly identified four genes (ZNF860, OR4F6, OR52E4, HERPUD2) that confer susceptibility to early-onset T2DM, MetS, or hyperuricemia.
We previously showed that four, five, or three SNPs were related to T2DM (P<1.44×10−4), MetS (P<0.05), or hyperuricemia (P<0.05) determined by multivariable logistic regression analysis with adjustment for age and sex after the initial EWAS screening of allele frequencies in both early-onset and late-onset individuals with these conditions (54–56). The relationship of two of four SNPs [rs138313632 (P=1.11×10−7), rs139012426 (P=4.29×10−5)] to T2DM were replicated (P<0.05) in the present study. The relationship of two of five SNPs [rs1007732 (P=0.0405), rs7350481 (P=3.17×10−5)] to MetS were replicated in the present study. The relationship of two of three SNPs [rs115445569 (P=0.0205), rs60854092 (P=0.0490)] to hyperuricemia were replicated in the present study. These results suggest that genetic variants associated with T2DM, MetS, or hyperuricemia differ, in part, between early-onset and late-onset individuals with these conditions.
There are several limitations to our study: i) Given that the results were not replicated, their validation will be necessary in independent study populations or in other ethnic groups. ii) It is possible that SNPs identified in the present study are in LD with other genetic variants in the same gene or in other nearby genes that are actually responsible for the development of T2DM, MetS, or hyperuricemia. iii) The functional relevance of identified SNPs to the pathogenesis of T2DM, MetS, or hyperuricemia remains to be elucidated.
In conclusion, we have newly identified four genes (ZNF860, OR4F6, OR52E4, HERPUD2) that confer susceptibility to early-onset T2DM, MetS, or hyperuricemia. Determination of genotypes for the SNPs in these genes may prove informative for assessment of the genetic risk for T2DM, MetS, or hyperuricemia in Japanese.
Acknowledgements
Not applicable.
Funding
This study was supported by CREST (JPMJCR1302), Japan Science and Technology Agency (to YYamada, JS and IT).
Availability of data and materials
All data underlying the findings described in the article are available upon request from the corresponding author.
Authors' contributions
YYam contributed to conception and design of the study; to acquisition, analysis, and interpretation of the data; and to drafting of the manuscript. KK, MO, HH and TF each contributed to acquisition of the data and to revision of the manuscript. YYas, IT and JS contributed to analysis and interpretation of the data as well as to revision of the manuscript.
Ethics approval and consent to participate
The study protocol complied with the Declaration of Helsinki and was approved by the Committees on the Ethics of Human Research of Mie University Graduate School of Medicine, Hirosaki University Graduate School of Medicine, and participating hospitals (Gifu Prefectural Tajimi Hospital, Gifu Prefectural General Medical Center, Japanese Red Cross Nagoya First Hospital, Northern Mie Medical Center Inabe General Hospital, and Hirosaki Stroke and Rehabilitation Center). Written informed consent was obtained from all subjects.
Consent for publication
All authors approved submission of the final version of the article for publication.
Competing interests
The authors declare that they have no competing interests.
References
|
Kharroubi AT and Darwish HM: Diabetes mellitus: The epidemic of the century. World J Diabetes. 6:850–867. 2015. View Article : Google Scholar : PubMed/NCBI | |
|
Expert Committee on the Diagnosis and Classification of Diabetes Mellitus: Report of the expert committee on the diagnosis and classification of diabetes mellitus. Diabetes Care. 26(Suppl 1): S5–S20. 2003.PubMed/NCBI | |
|
Ismail-Beigi F: Clinical practice. Glycemic management of type 2 diabetes mellitus. N Engl J Med. 366:1319–1327. 2012. | |
|
Emerging Risk Factors Collaboration1; Sarwar N, Gao P, Seshasai SR, Gobin R, Kaptoge S, Di Angelantonio E, Ingelsson E, Lawlor DA, Selvin E, et al: Diabetes mellitus, fasting blood glucose concentration, and risk of vascular disease: A collaborative meta-analysis of 102 prospective studies. Lancet. 375:2215–2222. 2010. View Article : Google Scholar : PubMed/NCBI | |
|
Stumvoll M, Goldstein BJ and van Haeften TW: Type 2 diabetes: Principles of pathogenesis and therapy. Lancet. 365:1333–1346. 2005. View Article : Google Scholar : PubMed/NCBI | |
|
Almgren P, Lehtovirta M, Isomaa B, Sarelin L, Taskinen MR, Lyssenko V, Tuomi T and Groop L; Botnia Study Group: Heritability and familiality of type 2 diabetes and related quantitative traits in the Botnia Study. Diabetologia. 54:2811–2819. 2011. View Article : Google Scholar : PubMed/NCBI | |
|
Prasad RB and Groop L: Genetics of type 2 diabetes-pitfalls and possibilities. Genes (Basel). 6:87–123. 2015. View Article : Google Scholar : PubMed/NCBI | |
|
Wellcome Trust Case Control Consortium: Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature. 447:661–678. 2007. View Article : Google Scholar : PubMed/NCBI | |
|
Sladek R, Rocheleau G, Rung J, Dina C, Shen L, Serre D, Boutin P, Vincent D, Belisle A, Hadjadj S, et al: A genome-wide association study identifies novel risk loci for type 2 diabetes. Nature. 445:881–885. 2007. View Article : Google Scholar : PubMed/NCBI | |
|
Dupuis JI, Langenberg C, Prokopenko I, Saxena R, Soranzo N, Jackson AU, Wheeler E, Glazer NL, Bouatia-Naji N, Gloyn AL, et al: New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk. Nat Genet. 42:105–116. 2010. View Article : Google Scholar : PubMed/NCBI | |
|
Voight BF, Scott LJ, Steinthorsdottir V, Morris AP, Dina C, Welch RP, Zeggini E, Huth C, Aulchenko YS, Thorleifsson G, et al: MAGIC investigators; GIANT Consortium: Twelve type 2 diabetes susceptibility loci identified through large-scale association analysis. Nat Genet. 42:579–589. 2010. View Article : Google Scholar : PubMed/NCBI | |
|
Morris AP, Voight BF, Teslovich TM, Ferreira T, Segrè AV, Steinthorsdottir V, Strawbridge RJ, Khan H, Grallert H, Mahajan A, et al: DIAbetes Genetics Replication And Meta-analysis (DIAGRAM) Consortium: Large-scale association analysis provides insights into the genetic architecture and pathophysiology of type 2 diabetes. Nat Genet. 44:981–990. 2012. View Article : Google Scholar : PubMed/NCBI | |
|
Scott RA, Lagou V, Welch RP, Wheeler E, Montasser ME, Luan J, Mägi R, Strawbridge RJ, Rehnberg E, Gustafsson S, et al: DIAbetes Genetics Replication and Meta-analysis (DIAGRAM) Consortium: Large-scale association analyses identify new loci influencing glycemic traits and provide insight into the underlying biological pathways. Nat Genet. 44:991–1005. 2012. View Article : Google Scholar : PubMed/NCBI | |
|
Ng MC, Shriner D, Chen BH, Li J, Chen WM, Guo X, Liu J, Bielinski SJ, Yanek LR, Nalls MA, et al: FIND Consortium; eMERGE Consortium; DIAGRAM Consortium; MuTHER Consortium; MEta-analysis of type 2 DIabetes in African Americans Consortium: Meta-analysis of genome-wide association studies in African Americans provides insights into the genetic architecture of type 2 diabetes. PLoS Genet. 10:e10045172014. View Article : Google Scholar : PubMed/NCBI | |
|
Cho YS, Chen CH, Hu C, Long J, Ong RT, Sim X, Takeuchi F, Wu Y, Go MJ, Yamauchi T, et al: MuTHER Consortium: Meta-analysis of genome-wide association studies identifies eight new loci for type 2 diabetes in east Asians. Nat Genet. 44:67–72. 2011. View Article : Google Scholar : PubMed/NCBI | |
|
Mahajan A, Go MJ, Zhang W, Below JE, Gaulton KJ, Ferreira T, Horikoshi M, Johnson AD, Ng MC, Prokopenko I, et al: DIAbetes Genetics Replication And Meta-analysis (DIAGRAM) Consortium; Asian Genetic Epidemiology Network Type 2 Diabetes (AGEN-T2D) Consortium; South Asian Type 2 Diabetes (SAT2D) Consortium; Mexican American Type 2 Diabetes (MAT2D) Consortium; Type 2 Diabetes Genetic Exploration by Nex-generation sequencing in muylti-Ethnic Samples (T2D-GENES) Consortium: Genome-wide trans-ancestry meta-analysis provides insight into the genetic architecture of type 2 diabetes susceptibility. Nat Genet. 46:234–244. 2014. View Article : Google Scholar : PubMed/NCBI | |
|
Zhao W, Rasheed A, Tikkanen E, Lee JJ, Butterworth AS, Howson JMM, Assimes TL, Chowdhury R, Orho-Melander M, Damrauer S, et al: CHD Exome+Consortium; EPIC-CVD Consortium; EPIC-Interact Consortium; Michigan Biobank: Identification of new susceptibility loci for type 2 diabetes and shared etiological pathways with coronary heart disease. Nat Genet. 49:1450–1457. 2017. View Article : Google Scholar : PubMed/NCBI | |
|
Unoki H, Takahashi A, Kawaguchi T, Hara K, Horikoshi M, Andersen G, Ng DP, Holmkvist J, Borch-Johnsen K, Jørgensen T, et al: SNPs in KCNQ1 are associated with susceptibility to type 2 diabetes in East Asian and European populations. Nat Genet. 40:1098–1102. 2008. View Article : Google Scholar : PubMed/NCBI | |
|
Yasuda K, Miyake K, Horikawa Y, Hara K, Osawa H, Furuta H, Hirota Y, Mori H, Jonsson A, Sato Y, et al: Variants in KCNQ1 are associated with susceptibility to type 2 diabetes mellitus. Nat Genet. 40:1092–1097. 2008. View Article : Google Scholar : PubMed/NCBI | |
|
Yamauchi T, Hara K, Maeda S, Yasuda K, Takahashi A, Horikoshi M, Nakamura M, Fujita H, Grarup N, Cauchi S, et al: A genome-wide association study in the Japanese population identifies susceptibility loci for type 2 diabetes at UBE2E2 and C2CD4A-C2CD4B. Nat Genet. 42:864–868. 2010. View Article : Google Scholar : PubMed/NCBI | |
|
Imamura M, Takahashi A, Yamauchi T, Hara K, Yasuda K, Grarup N, Zhao W, Wang X, Huerta-Chagoya A, Hu C, et al: Genome-wide association studies in the Japanese population identify seven novel loci for type 2 diabetes. Nat Commun. 7:105312016. View Article : Google Scholar : PubMed/NCBI | |
|
Alberti KG, Eckel RH, Grundy SM, Zimmet PZ, Cleeman JI, Donato KA, Fruchart JC, James WP, Loria CM and Smith SC Jr; International Diabetes Federation Task Force on Epidemiology and Prevention; Hational Heart, Lung, and Blood Institute; American Heart Association; World Heart Federation; International Atherosclerosis Society; International Association for the Study of Obesity: Harmonizing the metabolic syndrome: A joint interim statement of the International Diabetes Federation Task Force on Epidemiology and Prevention; National Heart, Lung, and Blood Institute; American Heart Association; World Heart Federation; International Atherosclerosis Society; and International Association for the Study of Obesity. Circulation. 120:1640–1645. 2009. View Article : Google Scholar : PubMed/NCBI | |
|
Grundy SM, Cleeman JI, Daniels SR, Donato KA, Eckel RH, Franklin BA, Gordon DJ, Krauss RM, Savage PJ, Smith SC Jr, et al: American Heart Association; National Heart, Lung, and Blood Institute: Diagnosis and management of the metabolic syndrome: An American Heart Association/National Heart, Lung, and Blood Institute Scientific Statement. Circulation. 112:2735–2752. 2005. View Article : Google Scholar : PubMed/NCBI | |
|
Esposito K, Chiodini P, Colao A, Lenzi A and Giugliano D: Metabolic syndrome and risk of cancer: A systematic review and meta-analysis. Diabetes Care. 35:2402–2411. 2012. View Article : Google Scholar : PubMed/NCBI | |
|
Abou Ziki MD and Mani A: Metabolic syndrome: Genetic insights into disease pathogenesis. Curr Opin Lipidol. 27:162–171. 2016. View Article : Google Scholar : PubMed/NCBI | |
|
Vattikuti S, Guo J and Chow CC: Heritability and genetic correlations explained by common SNPs for metabolic syndrome traits. PLoS Genet. 8:e10026372012. View Article : Google Scholar : PubMed/NCBI | |
|
Kraja AT, Vaidya D, Pankow JS, Goodarzi MO, Assimes TL, Kullo IJ, Sovio U, Mathias RA, Sun YV, Franceschini N, et al: A bivariate genome-wide approach to metabolic syndrome: STAMPEED consortium. Diabetes. 60:1329–1339. 2011. View Article : Google Scholar : PubMed/NCBI | |
|
Kristiansson K, Perola M, Tikkanen E, Kettunen J, Surakka I, Havulinna AS, Stancáková A, Barnes C, Widen E, Kajantie E, et al: Genome-wide screen for metabolic syndrome susceptibility Loci reveals strong lipid gene contribution but no evidence for common genetic basis for clustering of metabolic syndrome traits. Circ Cardiovasc Genet. 5:242–249. 2012. View Article : Google Scholar : PubMed/NCBI | |
|
Tekola-Ayele F, Doumatey AP, Shriner D, Bentley AR, Chen G, Zhou J, Fasanmade O, Johnson T, Oli J, Okafor G, et al: Genome-wide association study identifies African-ancestry specific variants for metabolic syndrome. Mol Genet Metab. 116:305–313. 2015. View Article : Google Scholar : PubMed/NCBI | |
|
Zabaneh D and Balding DJ: A genome-wide association study of the metabolic syndrome in Indian Asian men. PLoS One. 5:e119612010. View Article : Google Scholar : PubMed/NCBI | |
|
Zhu Y, Zhang D, Zhou D, Li Z, Li Z, Fang L, Yang M, Shan Z, Li H, Chen J, et al: Susceptibility loci for metabolic syndrome and metabolic components identified in Han Chinese: A multi-stage genome-wide association study. J Cell Mol Med. 21:1106–1116. 2017. View Article : Google Scholar : PubMed/NCBI | |
|
Eraly SA, Vallon V, Rieg T, Gangoiti JA, Wikoff WR, Siuzdak G, Barshop BA and Nigam SK: Multiple organic anion transporters contribute to net renal excretion of uric acid. Physiol Genomics. 33:180–192. 2008. View Article : Google Scholar : PubMed/NCBI | |
|
Choi HK, Mount DB and Reginato AM; American College of Physicians; American Physiological Society: Pathogenesis of gout. Ann Intern Med. 143:499–516. 2005. View Article : Google Scholar : PubMed/NCBI | |
|
Feig DI, Kang DH and Johnson RJ: Uric acid and cardiovascular risk. N Engl J Med. 359:1811–1821. 2008. View Article : Google Scholar : PubMed/NCBI | |
|
Fini MA, Elias A, Johnson RJ and Wright RM: Contribution of uric acid to cancer risk, recurrence, and mortality. Clin Transl Med. 1:162012. View Article : Google Scholar : PubMed/NCBI | |
|
Reginato AM, Mount DB, Yang I and Choi HK: The genetics of hyperuricaemia and gout. Nat Rev Rheumatol. 8:610–621. 2012. View Article : Google Scholar : PubMed/NCBI | |
|
Merriman TR: An update on the genetic architecture of hyperuricemia and gout. Arthritis Res Ther. 17:982015. View Article : Google Scholar : PubMed/NCBI | |
|
Wallace C, Newhouse SJ, Braund P, Zhang F, Tobin M, Falchi M, Ahmadi K, Dobson RJ, Marçano AC, Hajat C, et al: Genome-wide association study identifies genes for biomarkers of cardiovascular disease: Serum urate and dyslipidemia. Am J Hum Genet. 82:139–149. 2008. View Article : Google Scholar : PubMed/NCBI | |
|
Dehghan A, Köttgen A, Yang Q, Hwang SJ, Kao WL, Rivadeneira F, Boerwinkle E, Levy D, Hofman A, Astor BC, et al: Association of three genetic loci with uric acid concentration and risk of gout: A genome-wide association study. Lancet. 372:1953–1961. 2008. View Article : Google Scholar : PubMed/NCBI | |
|
Vitart V, Rudan I, Hayward C, Gray NK, Floyd J, Palmer CN, Knott SA, Kolcic I, Polasek O, Graessler J, et al: SLC2A9 is a newly identified urate transporter influencing serum urate concentration, urate excretion and gout. Nat Genet. 40:437–442. 2008. View Article : Google Scholar : PubMed/NCBI | |
|
Döring A, Gieger C, Mehta D, Gohlke H, Prokisch H, Coassin S, Fischer G, Henke K, Klopp N, Kronenberg F, et al: SLC2A9 influences uric acid concentrations with pronounced sex-specific effects. Nat Genet. 40:430–436. 2008. View Article : Google Scholar : PubMed/NCBI | |
|
Kolz M, Johnson T, Sanna S, Teumer A, Vitart V, Perola M, Mangino M, Albrecht E, Wallace C, Farrall M, et al: EUROSPAN Consortium; ENGAGE Consortium; PROCARDIS Consortium; KORA Study; WTCCC: Meta-analysis of 28,141 individuals identifies common variants within five new loci that influence uric acid concentrations. PLoS Genet. 5:e10005042009. View Article : Google Scholar : PubMed/NCBI | |
|
Yang Q, Köttgen A, Dehghan A, Smith AV, Glazer NL, Chen MH, Chasman DI, Aspelund T, Eiriksdottir G, Harris TB, et al: Multiple genetic loci influence serum urate levels and their relationship with gout and cardiovascular disease risk factors. Circ Cardiovasc Genet. 3:523–530. 2010. View Article : Google Scholar : PubMed/NCBI | |
|
Tin A, Woodward OM, Kao WH, Liu CT, Lu X, Nalls MA, Shriner D, Semmo M, Akylbekova EL, Wyatt SB, et al: CARe and CHARGE Consortia: Genome-wide association study for serum urate concentrations and gout among African Americans identifies genomic risk loci and a novel URAT1 loss-of-function allele. Hum Mol Genet. 20:4056–4068. 2011. View Article : Google Scholar : PubMed/NCBI | |
|
Li C, Li Z, Liu S, Wang C, Han L, Cui L, Zhou J, Zou H, Liu Z, Chen J, et al: Genome-wide association analysis identifies three new risk loci for gout arthritis in Han Chinese. Nat Commun. 6:70412015. View Article : Google Scholar : PubMed/NCBI | |
|
Köttgen A, Albrecht E, Teumer A, Vitart V, Krumsiek J, Hundertmark C, Pistis G, Ruggiero D, O'Seaghdha CM, Haller T, et al: LifeLines Cohort Study; CARDIoGRAM Consortium; DIAGRAM Consortium; ICBP Consortium; MAGIC Consortium: Genome-wide association analyses identify 18 new loci associated with serum urate concentrations. Nat Genet. 45:145–154. 2013. View Article : Google Scholar : PubMed/NCBI | |
|
Matsuo H, Yamamoto K, Nakaoka H, Nakayama A, Sakiyama M, Chiba T, Takahashi A, Nakamura T, Nakashima H, Takada Y, et al: Genome-wide association study of clinically defined gout identifies multiple risk loci and its association with clinical subtypes. Ann Rheum Dis. 75:652–659. 2016. View Article : Google Scholar : PubMed/NCBI | |
|
Nakayama A, Nakaoka H, Yamamoto K, Sakiyama M, Shaukat A, Toyoda Y, Okada Y, Kamatani Y, Nakamura T, Takada T, et al: Eurogout Consortium; Eurogout Consortium: GWAS of clinically defined gout and subtypes identifies multiple susceptibility loci that include urate transporter genes. Ann Rheum Dis. 76:869–877. 2017. View Article : Google Scholar : PubMed/NCBI | |
|
Ng MC, Lee SC, Ko GT, Li JK, So WY, Hashim Y, Barnett AH, Mackay IR, Critchley JA, Cockram CS, et al: Familial early-onset type 2 diabetes in Chinese patients: Obesity and genetics have more significant roles than autoimmunity. Diabetes Care. 24:663–671. 2001. View Article : Google Scholar : PubMed/NCBI | |
|
Bueno AC, Sun K, Martins CS, Elias Junior J, Miranda W, Tao C, Foss-Freitas MC, Barbieri MA, Bettiol H, de Castro M, et al: A novel ADIPOQ mutation (p.M40K) impairs assembly of high-molecular-weight adiponectin and is associated with early-onset obesity and metabolic syndrome. J Clin Endocrinol Metab. 99:E683–E693. 2014. View Article : Google Scholar : PubMed/NCBI | |
|
de Bruin C, Mericq V, Andrew SF, van Duyvenvoorde HA, Verkaik NS, Losekoot M, Porollo A, Garcia H, Kuang Y, Hanson D, et al: An XRCC4 splice mutation associated with severe short stature, gonadal failure, and early-onset metabolic syndrome. J Clin Endocrinol Metab. 100:E789–E798. 2015. View Article : Google Scholar : PubMed/NCBI | |
|
Zivná M, Hůlková H, Matignon M, Hodanová K, Vylet'al P, Kalbácová M, Baresová V, Sikora J, Blazková H, Zivný J, et al: Dominant renin gene mutations associated with early-onset hyperuricemia, anemia, and chronic kidney failure. Am J Hum Genet. 85:204–213. 2009. View Article : Google Scholar : PubMed/NCBI | |
|
Matsuo H, Ichida K, Takada T, Nakayama A, Nakashima H, Nakamura T, Kawamura Y, Takada Y, Yamamoto K, Inoue H, et al: Common dysfunctional variants in ABCG2 are a major cause of early-onset gout. Sci Rep. 3:20142013. View Article : Google Scholar : PubMed/NCBI | |
|
Yamada Y, Sakuma J, Takeuchi I, Yasukochi Y, Kato K, Oguri M, Fujimaki T, Horibe H, Muramatsu M, Sawabe M, et al: Identification of five genetic variants as novel determinants of type 2 diabetes mellitus in Japanese by exome-wide association studies. Oncotarget. 8:80492–80505. 2017. View Article : Google Scholar : PubMed/NCBI | |
|
Yamada Y, Sakuma J, Takeuchi I, Yasukochi Y, Kato K, Oguri M, Fujimaki T, Horibe H, Muramatsu M, Sawabe M, et al: Identification of rs7350481 at chromosome 11q23.3 as a novel susceptibility locus for metabolic syndrome in Japanese individuals by an exome-wide association study. Oncotarget. 8:39296–39308. 2017. View Article : Google Scholar : PubMed/NCBI | |
|
Yamada Y, Sakuma J, Takeuchi I, Yasukochi Y, Kato K, Oguri M, Fujimaki T, Horibe H, Muramatsu M, Sawabe M, et al: Identification of C21orf59 and ATG2A as novel determinants of renal function-related traits in Japanese by exome-wide association studies. Oncotarget. 8:45259–45273. 2017. View Article : Google Scholar : PubMed/NCBI | |
|
Yamada Y, Matsui K, Takeuchi I, Oguri M and Fujimaki T: Association of genetic variants with hypertension in a longitudinal population-based genetic epidemiological study. Int J Mol Med. 35:1189–1198. 2015. View Article : Google Scholar : PubMed/NCBI | |
|
Kuzuya T, Nakagawa S, Satoh J, Kanazawa Y, Iwamoto Y, Kobayashi M, Nanjo K, Sasaki A, Seino Y, Ito C, et al: Committee of the Japan Diabetes Society on the diagnostic criteria of diabetes mellitus: Report of the Committee on the classification and diagnostic criteria of diabetes mellitus. Diabetes Res Clin Pract. 55:65–85. 2002. View Article : Google Scholar : PubMed/NCBI | |
|
World Health Organization and International Diabetes Federation: Definition and diagnosis of diabetes mellitus and intermediate hyperglycemia: Report of a WHO/IDF consultation. World Health Organization; Geneva, Switzerland; pp. 1–46. 2006 | |
|
Grove ML, Yu B, Cochran BJ, Haritunians T, Bis JC, Taylor KD, Hansen M, Borecki IB, Cupples LA, Fornage M, et al: Best practices and joint calling of the HumanExome BeadChip: The CHARGE Consortium. PLoS One. 8:e680952013. View Article : Google Scholar : PubMed/NCBI | |
|
Anderson CA, Pettersson FH, Clarke GM, Cardon LR, Morris AP and Zondervan KT: Data quality control in genetic case-control association studies. Nat Protoc. 5:1564–1573. 2010. View Article : Google Scholar : PubMed/NCBI | |
|
Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA and Reich D: Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 38:904–909. 2006. View Article : Google Scholar : PubMed/NCBI | |
|
Teslovich TM, Musunuru K, Smith AV, Edmondson AC, Stylianou IM, Koseki M, Pirruccello JP, Ripatti S, Chasman DI, Willer CJ, et al: Biological, clinical and population relevance of 95 loci for blood lipids. Nature. 466:707–713. 2010. View Article : Google Scholar : PubMed/NCBI | |
|
Kim YJ, Go MJ, Hu C, Hong CB, Kim YK, Lee JY, Hwang JY, Oh JH, Kim DJ, Kim NH, et al: MAGIC consortium: Large-scale genome-wide association studies in East Asians identify new genetic loci influencing metabolic traits. Nat Genet. 43:990–995. 2011. View Article : Google Scholar : PubMed/NCBI |
