Open Access

Identification of 13 novel susceptibility loci for early-onset myocardial infarction, hypertension, or chronic kidney disease

  • Authors:
    • Yoshiji Yamada
    • Kimihiko Kato
    • Mitsutoshi Oguri
    • Hideki Horibe
    • Tetsuo Fujimaki
    • Yoshiki Yasukochi
    • Ichiro Takeuchi
    • Jun Sakuma
  • View Affiliations

  • Published online on: September 4, 2018     https://doi.org/10.3892/ijmm.2018.3852
  • Pages: 2415-2436
  • Copyright: © Yamada et al. This is an open access article distributed under the terms of Creative Commons Attribution License.

Metrics: Total Views: 0 (Spandidos Publications: | PMC Statistics: )
Total PDF Downloads: 0 (Spandidos Publications: | PMC Statistics: )


Abstract

Early‑onset cardiovascular and renal diseases have a strong genetic component. In the present study, exome‑wide association studies (EWASs) were performed to identify genetic variants that confer susceptibility to early‑onset myocardial infarction (MI), hypertension, or chronic kidney disease (CKD) in Japanese individuals. A total of 8,093 individuals aged ≤65 years was enrolled in the study. The EWASs for MI, hypertension, and CKD were performed in 6,926 subjects (1,152 cases, 5,774 controls), 8,080 subjects (3,444 cases, 4,636 controls), and 2,556 subjects (1,051 cases, 1,505 controls), respectively. Genotyping of single nucleotide polymorphisms (SNPs) was performed with Illumina Human Exome‑12 DNA Analysis BeadChip or Infinium Exome‑24 BeadChip arrays. The associations of allele frequencies for 31,245, 31,276, or 31,514 SNPs that passed quality control to MI, hypertension, and CKD, respectively, was examined with Fisher's exact test. Bonferroni's correction for statistical significance of association was applied to compensate for multiple comparisons of genotypes with MI, hypertension, or CKD. The EWASs of allele frequencies revealed that 25, 11, and 11 SNPs were significantly associated with MI (P<1.60x10‑6), hypertension (P<1.60x10‑6), or CKD (P<1.59x10‑6), respectively. Multivariable logistic regression analysis with adjustment for covariates showed that all 25, 11, and 11 SNPs were significantly associated with MI (P<0.0005), hypertension (P<0.0011), or CKD (P<0.0011), respectively. On examination of the results from previous genome‑wide association studies and linkage disequilibrium of the identified SNPs, 11 loci (TMOD4, COL6A3, ADGRL3‑CXCL8‑MARCH1, OR52E4, TCHP‑GIT2, CCDC63, 12q24.1, OAS3, PLCB2‑VPS33B, GOSR2, ZNF77), six loci (MOB3C‑TMOD4, COL6A3, COL6A5, CXCL8‑MARCH1, NFKBIL1‑6p21.3‑NCR3, PLCB2‑VPS33B), and seven loci (MOB3C‑TMOD4, COL6A3, COL6A5, ADGRL3‑CXCL8‑MARCH1, MUC17, PLCB2‑VPS33B, ZNF77) were identified as novel loci significantly associated with MI, hypertension, and CKD, respectively. Furthermore, six genes (TMOD4, COL6A3, CXCL8, MARCH1, PLCB2, VPS33B) were significantly associated with MI, hypertension and CKD; two genes (ADGRL3, ZNF77) with MI and CKD; and two genes (COL6A5, MOB3C) with hypertension and CKD. Therefore, 13 novel loci (MOB3C‑TMOD4, COL6A3, ADGRL3‑CXCL8‑MARCH1, OR52E4, TCHP‑GIT2, CCDC63, 12q24.1, OAS3, PLCB2‑VPS33B, ZNF77, COL6A5, NFKBIL1‑NCR3, MUC17) were identified that confer susceptibility to early‑onset MI, hypertension, or CKD. The determination of genotypes for the SNPs at these loci may provide informative for assessment of the genetic risk for MI, hypertension, or CKD.

Introduction

Coronary artery disease (CAD) and myocardial infarction (MI) are serious clinical conditions that remain the leading causes of mortality in the United States (1). Disease prevention is an important strategy for reducing the overall burden of these conditions, with the identification of biomarkers for disease risk being key for risk prediction and for potential intervention to reduce the chance of future coronary events.

In addition to conventional risk factors for CAD and MI, including hypertension, diabetes mellitus, and dyslipidemia, the importance of genetic factors has been highlighted (2). The heritability of CAD has been estimated to be 40-60% on the basis of family and twin studies (2,3). Genome-wide association studies (GWASs) in European-ancestry populations (4-10), African Americans (11), or Han Chinese (12,13) have identified various genes and loci that confer susceptibility to CAD or MI. A meta-analysis of GWASs for CAD among European-ancestry populations, which included low-frequency variants, identified 202 independent genetic variants at 129 loci with a false discovery rate of <5% (14). These genetic variants together account for ~28% of the heritability of CAD, showing that genetic susceptibility to this condition is largely determined by common variants with small effect sizes (2,14). A more recent meta-analysis for CAD in European-ancestry populations identified 304 independent genetic variants with a false discovery rate of <5%, with these variants explaining 21.2% of CAD heritability (15). Although several single nucleotide polymorphisms (SNPs) have been found to be significantly associated with MI in Japanese individuals (16,17), genetic variants that contribute to susceptibility to MI in the Japanese population remain to be fully elucidated.

Hypertension is a complex multifactorial disorder that results from an interaction between genetic background and both lifestyle and environmental factors (18). The heritability of hypertension has been estimated to be 84% and that of systolic blood pressure (BP) and diastolic BP has been estimated to be 53 and 49%, respectively (19). As hypertension is a major risk factor for CAD and stroke (20), personalized prevention of hypertension is an important public health goal.

GWASs have implicated various loci and genes as influencing BP or conferring susceptibility to hypertension in populations in Europe (6,21-24) and African (25,26) and in individuals from East Asia (27). A trans-ancestry GWAS in individuals of European and East or South Asian ancestry identified 12 loci associated with BP (28). A previous trans-ancestry meta-analysis in individuals of European and South Asian ancestry identified 30 susceptibility loci for BP or hypertension (29). Although an SNP in ADD2 was shown to be a susceptibility locus for hypertension in Japanese individuals (30), genetic variants that confer susceptibility to this condition in Japanese individuals await definitive identification.

Chronic kidney disease (CKD) is also a global public health problem, with affected individuals being at increased risk for end-stage renal disease (ESRD), cardiovascular disease and premature death (31,32). The identification of biomarkers for CKD risk is important to prevent progression to ESRD and to reduce the chance of future cardiovascular adverse events (33). In addition to conventional risk factors for CKD, including diabetes mellitus and hypertension (34), studies have shown the importance of genetic factors in renal function and in the development of CKD (35). The heritability of glomerular filtration rate for creatinine clearance has been estimated to be 52% (36).

GWASs have implicated several genes and loci in renal function or predisposition to CKD or ESRD in populations of European (37-42) or African (43) ancestry, and in renal function-related traits in East Asian populations (44). A meta-analysis of GWASs in European populations identified 53 loci that were significantly associated with estimated glomerular filtration rate (eGFR) (45). Although several SNPs have been shown to be associated with renal function in East Asian populations (44), genetic variants that contribute to predisposition to CKD in the Japanese population remain to be fully elucidated.

A previous study of monozygotic and dizygotic twins found that CAD-associated mortality rate at younger ages was significantly influenced by genetic factors in men and women, whereas the genetic effect was less pronounced at older ages (46). A family history of MI is also more apparent in individuals with early-onset MI compared with those with late-onset MI, suggestive of increased heritability in the former individuals (47,48). Similar to early-onset MI, early-onset forms of hypertension (49) and CKD (35) have been found to have a strong genetic component. As the genetic contribution may be higher in early-onset than in late-onset MI, hypertension and CKD, the statistical power of genetic association studies may be increased by focusing on subjects with early-onset disease (2,35,49).

In the present study, exome-wide association studies (EWASs) for MI, hypertension, and CKD were performed with the use of human exome array-based genotyping methods in order to identify genetic variants that confer susceptibility to these conditions in Japanese individuals. To increase the statistical power of the EWASs, subjects with early-onset disease were examined.

Materials and methods

Study subjects

In our previous EWASs, the median (mean) ages of subjects with MI, hypertension, and CKD were 67 (67.0), 68 (67.2), and 70 (69.9) years, respectively (50-52). Therefore patients with an age of ≤65 years were defined as individuals with early-onset forms of these conditions in the present study. A total of 8,093 Japanese individuals aged ≤65 years (subjects with MI, hypertension, or CKD and corresponding controls) was enrolled in the present study. In the EWASs for MI, hypertension, and CKD, 6,926 individuals (1,152 subjects with MI, 5,774 controls), 8,080 individuals (3,444 subjects with hypertension, 4,636 controls), and 2,556 individuals (1,051 subjects with CKD, 1,505 controls) were examined, respectively. The individuals were recruited as those who visited outpatient clinics 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, Japan; Northern Mie Medical Center Inabe General Hospital, Inabe, Japan; and Hirosaki University Hospital and Hirosaki Stroke and Rehabilitation Center, Hirosaki, Japan) due to the presence of various symptoms or for an annual health checkup between 2002 and 2014, or who were community-dwelling individuals recruited to a population-based cohort study in Inabe between 2010 and 2014 (53).

For the MI study, the diagnosis of MI was based on typical electrocardiographic changes and on increases in the serum activity of creatine kinase (MB isozyme) and in the serum concentration of troponin T. The diagnosis was confirmed by identification of the responsible stenosis in any of the major coronary arteries or in the left main trunk by coronary angiography. The control individuals had no history of MI, CAD, aortic aneurysm, or peripheral artery disease; of ischemic or hemorrhagic stroke; or of other atherosclerotic, thrombotic, embolic, or hemorrhagic disorders. Although certain control individuals had conventional risk factors for MI, including hypertension, diabetes mellitus, dyslipidemia, and CKD, they did not have cardiovascular complications.

In the hypertension study, the subjects with hypertension either had a systolic BP of ≥140 mmHg or diastolic BP of ≥90 mmHg (or both), or had taken antihypertensive medication. Individuals with severe valvular heart disease, congenital malformations of the heart or vessels, renal or endocrinologic diseases that cause secondary hypertension, or drug-induced hypertension were excluded from the study. The control individuals had a systolic BP of <140 mmHg and diastolic BP of <90 mmHg and no history of hypertension or of taking antihypertensive medication. BP was measured at least twice by a skilled physician or nurse, with subjects having first rested in the sitting position for >5 min.

For the CKD study, glomerular filtration rate was estimated with the use of the simplified prediction equation derived from a modified version of that described in the Modification of Diet in Renal Disease Study, as suggested by the Japanese Society of Nephrology (54): eGFR (ml min−1 1.73 m−2) = 194 × [age (years)] −0.287 × [serum creatinine (mg/dl)] −1.094 × [0.739 if female]. The National Kidney Foundation-Kidney Disease Outcomes Quality Initiative guidelines recommend a diagnosis of CKD if eGFR is <60 ml min−1 1.73 m−2 (55). Nonlinear relations between eGFR and the risk of adverse events, including mortality, cardiovascular episodes, and hospitalization, have been demonstrated, with an increased risk being associated with an eGFR of <60 ml min−1 1.73 m−2 and the risk markedly increasing further when values fall <45 ml min−1 1.73 m−2 (56). Therefore, the criterion of an eGFR of <60 ml min−1 1.73 m−2 was used for the diagnosis of CKD. The control individuals for the CKD study had an eGFR of ≥90 ml min−1 1.73 m−2 and had no functional or structural abnormalities of the kidneys or a history of renal disease. Although certain control individuals had hypertension, diabetes mellitus, dyslipidemia, or hyperuricemia, they had no renal complications.

A total of 785 subjects with both MI and hypertension and 3,946 controls overlapped between the corresponding studies, as did 324 subjects with both MI and CKD and 1,127 controls, and 750 subjects with both hypertension and CKD and 955 controls.

EWASs

Venous blood (5 or 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 with the use of a Genomix DNA extraction kit (Talent Srl, Trieste, Italy) or SMITEST EX-R&D kit (Medical & Biological Laboratories, Co., Ltd., Nagoya, Japan). The EWASs for MI (1,152 cases, 5,774 controls), hypertension (3,444 cases, 4,636 controls), and CKD (1,051 cases, 1,505 controls) 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 consists of ~244,000 SNPs from European, African, Chinese, and Hispanic individuals (57). SNPs contained in only one of the exome arrays (~2.6% of all SNPs) were excluded from analysis. Quality control was performed (58) 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 each sample, and those for which sex phenotype in the clinical records was inconsistent with genetic sex chromosomes were discarded; iii) duplicate samples and cryptic relatedness were assessed by calculation of identity by descent; all pairs of DNA samples showing an identity by descent of >0.1875 were inspected and one sample from each pair was excluded; iv) the frequency of heterozygosity for SNPs was calculated for all samples, and those with 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 nonpolymorphic SNPs or SNPs with a minor allele frequency of <1.0%; vi) SNPs with genotype distributions deviating significantly (P<0.01) from Hardy-Weinberg equilibrium in control individuals were discarded; vii) genotype data were examined for population stratification by principal components analysis (59), and population outliers were excluded from the analysis. A total of 31,245, 31,276, or 31,514 SNPs passed quality control for the EWASs of MI, hypertension, and CKD, respectively, and these SNPs were subjected to analyses.

Statistical analysis

For analysis of the characteristics of the study subjects, quantitative data were compared between subjects with MI, hypertension, or CKD and corresponding controls with an unpaired t-test. Categorical data were compared between the two groups with Pearson's χ2 test. Allele frequencies were estimated by the gene counting method, and Fisher's exact test was applied to identify departure from Hardy-Weinberg equilibrium. In the EWASs, the relation of allele frequencies of each SNP to MI, hypertension, or CKD was examined with Fisher's exact test. To compensate for multiple comparisons of genotypes with MI, hypertension, or CKD, Bonferroni's correction was applied for statistical significance of association. As 31,245, 31,276, and 31,514 SNPs that passed quality control were examined in the EWASs for MI, hypertension, and CKD, respectively, P<1.60×10−6 (0.05/31,245 and 31,276) for MI and hypertension or P<1.59×10−6 (0.05/31,514) for CKD was considered statistically significant. The inflation factor (λ) was 0.95 for MI, 1.03 for hypertension, and 1.07 for CKD. Multivariable logistic regression analysis was performed with MI as a dependent variable and independent variables included age, sex (0, woman; 1, man), the prevalence of hypertension, diabetes mellitus, and dyslipidemia (0, no history of these conditions; 1, positive history), and the genotype of each SNP. Similar analysis was performed with hypertension as a dependent variable and independent variables of age, sex, and the genotype of each SNP; and with CKD as a dependent variable and independent variables of age, sex, the prevalence of hypertension and diabetes mellitus, and the genotype of each SNP. The genotypes of the SNPs 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, odds ratio, 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 association between the genotypes of SNPs and intermediate phenotypes of MI was examined with Pearson's χ2 test. The association between genotypes of SNPs and systolic or diastolic BP, the serum concentration of creatinine, or eGFR 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 6,926 subjects enrolled in the MI study are shown in Table I. The age, the frequency of men, and the prevalence of obesity, hypertension, diabetes mellitus, dyslipidemia, CKD, and hyperuricemia, in addition to body mass index (BMI), systolic and diastolic BP, fasting plasma glucose (FPG) level, blood glycosylated hemoglobin (hemoglobin A1c) content, and the serum concentrations of triglycerides, creatinine, and uric acid were higher, whereas the serum concentration of high density lipoprotein (HDL)-cholesterol and eGFR were lower, in subjects with MI than in controls.

Table I

Characteristics of subjects with MI and control individuals.

Table I

Characteristics of subjects with MI and control individuals.

CharacteristicControlMyocardial infarctionP-value
Subjects (n)5,7741,152
Age (years)50.6±10.255.6±7.4<0.0001
Sex (men/women, %)52.1/47.985.4/14.6<0.0001
Smoking (%)42.546.70.0105
Obesity (%)31.042.7<0.0001
Body mass index (kg/m2)23.2±3.524.6±3.5<0.0001
Hypertension (%)31.768.9<0.0001
Systolic BP (mmHg)121±18139±26<0.0001
Diastolic BP (mmHg)75±1378±16<0.0001
Diabetes mellitus (%)12.761.3<0.0001
Fasting plasma glucose (mmol/l)5.66±1.787.66±3.39<0.0001
Blood hemoglobin A1c (%)5.72±0.966.88±1.70<0.0001
Dyslipidemia (%)56.984.9<0.0001
Serum triglycerides (mmol/l)1.32±0.981.77±1.25<0.0001
Serum HDL-cholesterol (mmol/l)1.65±0.451.17±0.34<0.0001
Serum LDL-cholesterol (mmol/l)3.18±0.833.15±0.980.8593
Chronic kidney disease (%)10.329.1<0.0001
Serum creatinine (µmol/l)69.8±61.091.9±98.1<0.0001
eGFR (ml min−1 1.73 m−2)78.7±17.171.2±27.8<0.0001
Hyperuricemia (%)15.225.6<0.0001
Serum uric acid (µmol/l)321±89356±105<0.0001

[i] Quantitative data are presented as the mean ± standard deviation and were compared between subjects with MI and controls using the unpaired t-test. Categorical data were compared between the two groups with Pearson's χ2 test. Based on Bonferroni's correction, P<0.0025 (0.05/20) was considered statistically significant. Obesity was defined as a body mass index of ≥25 kg/m2; hypertension as a systolic BP of ≥140 mmHg, diastolic BP of ≥90 mmHg, or the taking of anti-hypertensive medication; diabetes mellitus as a fasting plasma glucose level of ≥6.93 mmol/l, blood hemoglobin A1c content of ≥6.5%, or the taking of antidiabetic medication; dyslipidemia as a serum triglyceride concentration of ≥1.65 mmol/l, serum HDL-cholesterol <1.04 mmol/l, serum LDL-cholesterol ≥3.64 mmol/l, or the taking of anti-dyslipidemic medication; chronic kidney disease as an eGFR of <60 ml min−1 1.73 m−2; and hyperuricemia as a serum uric acid concentration of >416 µmol/l or the taking of uric acid-lowering medication. MI, myocardial infarction; BP, blood pressure; HDL, high density lipoprotein; LDL, low density lipoprotein; eGFR, estimated glomerular filtration rate.

The characteristics of the 8,080 subjects enrolled in the hypertension study are shown in Table II. The age, the frequency of men, and the prevalence of obesity, diabetes mellitus, dyslipidemia, CKD, and hyperuricemia, in addition to the BMI, FPG level, blood hemoglobin A1c content, and the serum concentrations of triglycerides, creatinine, and uric acid were higher, whereas the serum concentration of HDL-cholesterol and eGFR were lower, in subjects with hypertension than in controls.

Table II

Characteristics of subjects with hypertension and control individuals.

Table II

Characteristics of subjects with hypertension and control individuals.

CharacteristicControlHypertensionP-value
Subjects (n)46363444
Age (years)49.2±10.355.8±7.6<0.0001
Sex (men/women, %)52.0/48.067.8/32.2<0.0001
Smoking (%)41.742.90.2829
Obesity (%)24.344.4<0.0001
Body mass index (kg/m2)22.6±3.324.5±3.7<0.0001
Systolic BP (mmHg)114±12144±24<0.0001
Diastolic BP (mmHg)70±1085±14<0.0001
Diabetes mellitus (%)11.139.5<0.0001
Fasting plasma glucose (mmol/l)5.61±1.786.72±2.83<0.0001
Blood hemoglobin A1c (%)5.71±0.966.30±1.48<0.0001
Dyslipidemia (%)51.375.5<0.0001
Serum triglycerides (mmol/l)1.25±0.981.69±1.19<0.0001
Serum HDL-cholesterol (mmol/l)1.63±0.461.41±0.45<0.0001
Serum LDL-cholesterol (mmol/l)3.15±0.853.15±0.880.9206
Chronic kidney disease (%)8.124.0<0.0001
Serum creatinine (µmol/l)65.4±23.989.3±117.6<0.0001
eGFR (ml min−1 1.73 m−2)80.4±18.872.5±21.3<0.0001
Hyperuricemia (%)11.325.2<0.0001
Serum uric acid (µmol/l)311±86349±99<0.0001

[i] Quantitative data are presented as the mean ± standard deviation and were compared between subjects with hypertension and controls with the unpaired t-test. Categorical data were compared between the two groups with Pearson's χ2 test. Based on Bonferroni's correction, P<0.0026 (0.05/19) was considered statistically significant. Obesity was defined as a body mass index of ≥25 kg/m2; hypertension as a systolic BP of ≥140 mmHg, diastolic BP of ≥90 mmHg, or the taking of antihypertensive medication; diabetes mellitus as a fasting plasma glucose level of ≥6.93 mmol/l, blood hemoglobin A1c content of ≥6.5%, or the taking of antidiabetic medication; dyslipidemia as a serum triglyceride concentration of ≥1.65 mmol/l, serum HDL-cholesterol <1.04 mmol/l, serum LDL-cholesterol ≥3.64 mmol/l, or the taking of anti-dyslipidemic medication; chronic kidney disease as an eGFR of <60 ml min−1 1.73 m−2; and hyperuricemia as a serum uric acid concentration of >416 µmol/l or the taking of uric acid-lowering medication. BP, blood pressure; HDL, high density lipoprotein; LDL, low density lipoprotein; eGFR, estimated glomerular filtration rate.

The characteristics of the 2,556 subjects enrolled in the CKD study are shown in Table III. The age, the frequency of men, and the prevalence of obesity, hypertension, diabetes mellitus, dyslipidemia, and hyperuricemia, in addition to the BMI, systolic and diastolic BP, FPG level, blood hemoglobin A1c content, and the serum concentrations of triglycerides and uric acid were higher, whereas the serum concentration of HDL-cholesterol was lower, in subjects with CKD than in controls.

Table III

Characteristics of subjects with CKD and control individuals.

Table III

Characteristics of subjects with CKD and control individuals.

CharacteristicControlCKDP-value
Subjects (n)1,5051,051
Age (years)48.5±10.357.7±6.5<0.0001
Sex (men/women, %)53.8/46.268.9/31.1<0.0001
Smoking (%)43.237.30.0028
Obesity (%)29.138.9<0.0001
Body mass index (kg/m2)23.1±3.824.1±3.5<0.0001
Hypertension (%)36.571.4<0.0001
Systolic BP (mmHg)126±23141±29<0.0001
Diastolic BP (mmHg)76±1480±16<0.0001
Diabetes mellitus (%)24.845.2<0.0001
Fasting plasma glucose (mmol/l)6.27±2.726.88±3.16<0.0001
Blood hemoglobin A1c (%)6.08±1.516.42±1.51<0.0001
Dyslipidemia (%)56.475.3<0.0001
Serum triglycerides (mmol/l)1.35±1.211.69±1.07<0.0001
Serum HDL-cholesterol (mmol/l)1.56±0.471.38±0.48<0.0001
Serum LDL-cholesterol (mmol/l)3.08±0.853.13±0.930.0738
Serum creatinine (µmol/l)51.7±9.3150.3±198.0<0.0001
eGFR (ml min−1 1.73 m−2)103.0±19.146.9±14.6<0.0001
Hyperuricemia (%)8.540.6<0.0001
Serum uric acid (µmol/l)297±80381±106<0.0001

[i] Quantitative data are presented as the mean ± standard deviation and were compared between subjects with CKD and controls with the unpaired t-test. Categorical data were compared between the two groups with Pearson's χ2 test. Based on Bonferroni's correction, P<0.0026 (0.05/19) was considered statistically significant. Obesity was defined as a body mass index of ≥25 kg/m2; hypertension as a systolic BP of ≥140 mmHg, diastolic BP of ≥90 mmHg, or the taking of antihypertensive medication; diabetes mellitus as a fasting plasma glucose level of ≥6.93 mmol/l, blood hemoglobin A1c content of ≥6.5%, or the taking of antidiabetic medication; dyslipidemia as a serum triglyceride concentration of ≥1.65 mmol/l, serum HDL-cholesterol <1.04 mmol/l, serum LDL-cholesterol ≥3.64 mmol/l, or the taking of anti-dyslipidemic medication; chronic kidney disease as an eGFR of <60 ml min−1 1.73 m−2; and hyperuricemia as a serum uric acid concentration of >416 µmol/l or the taking of uric acid-lowering medication. CKD, chronic kidney disease; BP, blood pressure; HDL, high density lipoprotein; LDL, low density lipoprotein; eGFR, estimated glomerular filtration rate.

EWASs for MI, hypertension, and CKD

The association between allele frequencies for the 31,245, 31,276, and 31,514 SNPs that passed quality control to MI, hypertension, and CKD, respectively, were examined using Fisher's exact test. Following Bonferroni's correction, 25 and 11 SNPs were significantly (P<1.60×10−6) associated with MI (Table IV) or hypertension (Table V), respectively, and 11 SNPs were significantly (P<1.59×10−6) associated with CKD (Table VI).

Table IV

25 SNPs significantly (P<1.60×10−6) associated with myocardial infarction in the exome-wide association study.

Table IV

25 SNPs significantly (P<1.60×10−6) associated with myocardial infarction in the exome-wide association study.

GeneSNPNucleotide substitutionaAmino acid substitutionChromosomePositionMAF (%)Allele ORP-value (allele frequency)
PLCB2rs200787930C/TE1106K15402892981.20.01 8.81×10−26
CXCL8rs188378669G/TE31*4737415681.20.01 1.33×10−25
MARCH1rs61734696G/TQ137K41641973031.20.01 1.50×10−25
VPS33Brs199921354C/TR80Q15910138411.20.01 2.07×10−25
TMOD4rs115287176G/AR277W11511709611.20.01 9.03×10−25
COL6A3rs146092501C/TE1386K22373718611.20.02 1.27×10−24
ZNF77rs146879198G/AR340*1929341091.20.02 1.27×10−24
ADGRL3rs192210727G/TR580I4619096151.30.07 1.12×10−20
ALDH2rs671G/AE504K1211180396227.61.54 2.67×10−19
ACAD10rs11066015G/A1211173020527.50.65 3.35×10−19
BRAPrs3782886A/G1211167268529.31.49 5.52×10−17
HECTD4rs11066280T/A1211237997929.01.49 6.48×10−17
HECTD4rs2074356C/T1211220759725.41.48 4.24×10−15
OR52E4rs11823828T/GF227L11588497336.61.43 1.86×10−12
rs12229654T/G1211097665722.51.39 1.45×10−10
NAA25rs12231744C/TR876K1211203925135.10.74 1.07×10−9
ATXN2rs7969300T/CN248S1211155590838.80.76 8.46×10−9
GOSR2rs1052586T/C174694109748.70.78 1.80×10−7
TCHPrs74416240G/A1210990479313.31.37 2.47×10−7
GIT2rs925368A/GN389S1210995317412.51.37 4.41×10−7
rs2523644A/G6313747078.11.49 4.75×10−7
OAS3rs2072134C/T1211297137117.61.33 5.28×10−7
CCDC63rs10774610T/C1211090243923.71.30 5.42×10−7
rs2596548G/T6313627695.41.60 6.01×10−7
CCDC63rs10849915T/C1211089581823.61.28 7.71×10−7

{ label (or @symbol) needed for fn[@id='tfn4-ijmm-42-05-2415'] } Allele frequencies were analyzed using Fisher's exact test.

a Major allele/minor allele. SNP, single nucleotide polymorphism; MAF, minor allele frequency; OR, odds ratio.

Table V

11 SNPs significantly (P<1.60×10−6) associated with hypertension in the exome-wide association study.

Table V

11 SNPs significantly (P<1.60×10−6) associated with hypertension in the exome-wide association study.

GeneSNPNucleotide substitutionaAmino acid substitutionChromosomePositionMAF (%)Allele OR (allele frequency)P-value
MARCH1rs61734696G/TQ137K41641973031.20.54 1.30×10−7
PLCB2rs200787930C/TE1106K15402892981.20.55 3.33×10−7
COL6A3rs146092501C/TE1386K22373718611.20.55 3.50×10−7
COL6A5rs200982668G/AE2501K31304708941.30.56 3.52×10−7
MOB3Crs139537100C/TR24Q1466150061.20.56 4.76×10−7
VPS33Brs199921354C/TR80Q15910138411.20.55 4.84×10−7
NCR3rs2515920T/A63159483817.30.81 5.02×10−7
TMOD4rs115287176G/AR277W11511709611.20.55 5.49×10−7
CXCL8rs188378669G/TE31*4737415681.20.56 6.93×10−7
rs769177G/A63157983417.20.81 1.35×10−6
NFKBIL1rs2071593C/T63154502218.80.82 1.59×10−6

{ label (or @symbol) needed for fn[@id='tfn6-ijmm-42-05-2415'] } Allele frequencies were analyzed using Fisher's exact test.

a Major allele/minor allele. SNP, single nucleotide polymorphism; MAF, minor allele frequency; OR, odds ratio.

Table VI

11 SNPs significantly (P<1.59×10−6) associated with chronic kidney disease in the exome-wide association study.

Table VI

11 SNPs significantly (P<1.59×10−6) associated with chronic kidney disease in the exome-wide association study.

GeneSNPNucleotide substitutionaAmino acid substitutionChromosomePositionMAF (%)Allele ORP-value (allele frequency)
COL6A5rs200982668G/AE2501K31304708941.30.21 5.08×10−9
MARCH1rs61734696G/TQ137K41641973031.20.22 1.97×10−8
MUC17rs78010183A/TT1305S71010353291.80.29 2.04×10−8
MOB3Crs139537100C/TR24Q1466150061.20.23 4.91×10−8
PLCB2rs200787930C/TE1106K15402892981.20.23 4.94×10−8
CXCL8rs188378669G/TE31*4737415681.20.23 7.83×10−8
VPS33Brs199921354C/TR80Q15910138411.20.23 7.86×10−8
TMOD4rs115287176G/AR277W11511709611.20.23 1.25×10−7
ADGRL3rs192210727G/TR580I4619096151.30.25 1.88×10−7
ZNF77rs146879198G/AR340*1929341091.20.24 2.95×10−7
COL6A3rs146092501C/TE1386K22373718611.20.24 2.99×10−7

{ label (or @symbol) needed for fn[@id='tfn8-ijmm-42-05-2415'] } Allele frequencies were analyzed using Fisher's exact test.

a Major allele/minor allele. SNP, single nucleotide polymorphism; MAF, minor allele frequency; OR, odds ratio.

Multivariable logistic regression analysis of the association between SNPs and MI, hypertension, or CKD

The association of the 25 SNPs identified in the EWAS for MI with this condition was further examined by multivariable logistic regression analysis with adjustment for age, sex, and the prevalence of hypertension, diabetes mellitus, and dyslipidemia (Table VII). All 25 SNPs were significantly [P<0.0005 (0.05/100) in at least one genetic model] associated with MI. The association between the 11 SNPs identified in the EWAS for hypertension and this condition was examined by multivariable logistic regression analysis with adjustment for age and sex (Table VIII). All 11 SNPs were significantly [P<0.0011 (0.05/44)] associated with hypertension. The association between the 11 SNPs identified by the EWAS for CKD and this condition was also examined by multivariable logistic regression analysis with adjustment for age, sex, and the prevalence of hypertension and diabetes mellitus (Table IX). All 11 SNPs were significantly [P<0.0011 (0.05/44)] associated with CKD.

Table VII

Association between SNPs and myocardial infarction as determined by multivariable logistic regression analysis.

Table VII

Association between SNPs and myocardial infarction as determined by multivariable logistic regression analysis.

GeneSNPDominant
Recessive
Additive 1
Additive 2
P-valueOR95% CIP-valueOR95% CIP-valueOR95% CIP-valueOR95% CI
PLCB2rs200787930C/T<0.00010.010.01-0.10<0.00010.010.01-0.10
CXCL8rs188378669G/T<0.00010.010.01-0.10<0.00010.010.01-0.10
MARCH1rs61734696G/T<0.00010.010.01-0.10<0.00010.010.01-0.10
VPS33Brs199921354C/T<0.00010.010.01-0.10<0.00010.010.01-0.10
TMOD4rs115287176G/A<0.00010.010.01-0.11<0.00010.010.01-0.11
COL6A3rs146092501C/T<0.00010.010.01-0.10<0.00010.010.01-0.10
ZNF77rs146879198G/A<0.00010.010.01-0.11<0.00010.010.01-0.11
ADGRL3rs192210727G/T<0.00010.060.02-0.160.1736<0.00010.060.02-0.160.9960
ALDH2rs671G/A<0.00012.001.56-2.55<0.00012.001.70-2.360.00401.461.13-1.90<0.00012.712.08-3.53
ACAD10rs11066015G/A<0.00010.510.39-0.65<0.00010.500.43-0.590.00490.690.53-0.89<0.00010.370.29-0.49
BRAPrs3782886A/G<0.00011.981.67-2.34<0.00011.881.48-2.39<0.00011.841.55-2.20<0.00012.581.99-3.35
HECTD4rs11066280T/A<0.00012.021.71-2.38<0.00011.921.51-2.45<0.00011.881.58-2.24<0.00012.662.05-3.46
HECTD4rs2074356C/T<0.00011.851.57-2.17<0.00011.921.47-2.51<0.00011.731.46-2.06<0.00012.481.87-3.28
OR52E4rs11823828T/G<0.00011.561.30-1.87<0.00012.161.74-2.680.01521.281.05-1.56<0.00012.461.93-3.13
rs12229654T/G<0.00011.641.39-1.92<0.00011.891.41-2.53<0.00011.541.30-1.82<0.00012.271.67-3.07
NAA25rs12231744C/T<0.00010.600.51-0.71<0.00010.490.37-0.64<0.00010.670.56-0.79<0.00010.400.30-0.53
ATXN2rs7969300T/C<0.00010.590.50-0.70<0.00010.550.43-0.71<0.00010.650.55-0.77<0.00010.440.34-0.57
GOSR2rs1052586T/C0.00020.710.59-0.85<0.00010.610.50-0.750.02360.800.66-0.97<0.00010.530.42-0.67
TCHPrs74416240G/A<0.00011.461.23-1.730.2849<0.00011.461.22-1.740.1399
GIT2rs925368A/G<0.00011.481.24-1.760.5500<0.00011.491.25-1.790.3309
rs2523644A/G<0.00011.741.41-2.140.5858<0.00011.761.42-2.170.4407
OAS3rs2072134C/T<0.00011.451.23-1.710.04081.511.02-2.25<0.00011.421.20-1.690.00931.711.14-2.55
CCDC63rs10774610T/C<0.00011.381.18-1.620.00851.491.11-2.010.00091.331.12-1.580.00091.681.24-2.29
rs2596548G/T<0.00011.951.53-2.480.3820<0.00012.041.60-2.610.4739
CCDC63rs10849915T/C<0.00011.381.18-1.630.01841.441.06-1.940.00061.341.14-1.590.00211.631.19-2.22

[i] Multivariable logistic regression analysis was performed with adjustment for age, sex, and the prevalence of hypertension, diabetes mellitus, and dyslipidemia. Based on Bonferroni's correction, P<0.0005 (0.05/100) was considered statistically significant. SNP, single nucleotide polymorphism; OR, odds ratio; CI, confidence interval.

Table VIII

Association between SNPs and hypertension as determined by multivariable logistic regression analysis.

Table VIII

Association between SNPs and hypertension as determined by multivariable logistic regression analysis.

GeneSNPDominant
Recessive
Additive 1
Additive 2
P-valueOR95% CIP-valueOR95% CIP-valueOR95% CIP-valueOR95% CI
MARCH1rs61734696G/T0.00020.620.48-0.790.00020.620.48-0.79
PLCB2rs200787930C/T0.00030.630.49-0.810.00030.630.49-0.81
COL6A3rs146092501C/T0.00020.610.47-0.790.00020.610.47-0.79
COL6A5rs200982668G/A0.00050.650.50-0.830.00050.650.50-0.83
MOB3Crs139537100C/T0.00060.640.50-0.830.00060.640.50-0.83
VPS33Brs199921354C/T0.00030.630.49-0.810.00030.630.49-0.81
NCR3rs2515920T/A<0.00010.740.67-0.820.0930<0.00010.740.66-0.820.02260.720.54-0.95
TMOD4rs115287176G/A0.00040.630.49-0.810.00040.630.49-0.81
CXCL8rs188378669G/T0.00040.630.49-0.810.00040.630.49-0.81
rs769177G/A<0.00010.740.67-0.820.0946<0.00010.740.67-0.830.02380.720.54-0.96
NFKBIL1rs2071593C/T<0.00010.740.67-0.820.2491<0.00010.740.67-0.820.0628

[i] Multivariable logistic regression analysis was performed with adjustment for age and sex. Based on Bonferroni's correction, P<0.0011 (0.05/44) was considered statistically significant. SNP, single nucleotide polymorphism; OR, odds ratio; CI, confidence interval.

Table IX

Association between SNPs and chronic kidney disease as determined by multivariable logistic regression analysis.

Table IX

Association between SNPs and chronic kidney disease as determined by multivariable logistic regression analysis.

GeneSNPDominant
Recessive
Additive 1
Additive 2
P-valueOR95% CIP-valueOR95% CIP-valueOR95% CIP-valueOR95% CI
COL6A5rs200982668G/A0.00020.250.12-0.520.00020.250.12-0.52
MARCH1rs61734696G/T0.00040.260.12-0.550.00040.260.12-0.55
MUC17rs78010183A/T0.00080.350.19-0.650.00080.350.19-0.65
MOB3Crs139537100C/T0.00060.270.13-0.570.00060.270.13-0.57
PLCB2rs200787930C/T0.00070.270.13-0.570.00070.270.13-0.57
CXCL8rs188378669G/T0.00040.260.12-0.550.00040.260.12-0.55
VPS33Brs199921354C/T0.00070.270.13-0.580.00070.270.13-0.58
TMOD4rs115287176G/A0.00070.270.13-0.580.00070.270.13-0.58
ADGRL3rs192210727G/T0.00060.270.12-0.570.48760.00070.270.13-0.570.9964
ZNF77rs146879198G/A0.00090.280.13-0.590.00090.280.13-0.59
COL6A3rs146092501C/T0.00100.280.13-0.600.00100.280.13-0.60

[i] Multivariable logistic regression analysis was performed with adjustment for age, sex, and the prevalence of hypertension and diabetes mellitus. Based on Bonferroni's correction, P<0.0011 (0.05/44) was considered statistically significant. SNP, single nucleotide polymorphism; OR, odds ratio; CI, confidence interval.

Association between SNPs associated with MI and intermediate phenotypes

The present study examined the associations between the 25 SNPs associated with MI and intermediate phenotypes of this condition, including hypertension, diabetes mellitus, hypertriglyceridemia, hypo-HDL-cholesterolemia, hyper-low density lipoprotein (LDL)-cholesterolemia, CKD, obesity, and hyperuricemia, using Pearson's χ2 test. The SNPs rs200787930 of PLCB2, rs188378669 of CXCL8, rs61734696 of MARCH1, rs199921354 of VPS33B, rs115287176 of TMOD4, rs146092501 of COL6A3, rs146879198 of ZNF77, and rs192210727 of ADGRL3 were significantly [P<0.0003 (0.05/200)] associated with hypo-HDL-cholesterolemia, hyper-LDL-cholesterolemia, and CKD, with rs61734696 also being associated with hypertension. The SNPs rs671 of ALDH2, rs11066015 of ACAD10, and rs11066280 and rs2074356 of HECTD4 were significantly associated with hypo-HDL-cholesterolemia, hyper-LDL-cholesterolemia, and hyperuricemia. The SNPs rs3782886 of BRAP and rs12229654 at 12q24.1 were significantly associated with hypo-HDL-cholesterolemia and hyperuricemia. The SNP rs11823828 of OR52E4 was significantly associated with diabetes mellitus, rs12231744 of NAA25 with hypo-HDL-cholesterolemia, rs7969300 of ATXN2 with hyper-triglyceridemia, and rs10774610 and rs10849915 of CCDC63 with hyperuricemia (Table X).

Table X

Association between SNPs associated with myocardial infarction and intermediate phenotypes.

Table X

Association between SNPs associated with myocardial infarction and intermediate phenotypes.

GeneSNPHypertensionDMHyper-TGHypo-HDLHyper-LDLCKDObesityHyperuricemia
PLCB2rs200787930C/T0.00040.00260.5013 <0.0001 <0.0001 <0.00010.04820.7764
CXCL8rs188378669G/T0.00070.00390.8922 <0.0001 <0.0001 <0.00010.02590.9412
MARCH1rs61734696G/T0.00020.00130.5876 <0.0001 <0.0001 <0.00010.03810.7835
VPS33Brs199921354C/T0.00060.00260.5939 <0.0001 <0.0001 <0.00010.03980.9038
TMOD4rs115287176G/A0.00060.00730.9087 <0.0001 <0.0001 <0.00010.05780.8457
COL6A3rs146092501 C/T0.00040.00580.8385 <0.0001 <0.0001 <0.00010.04590.8305
ZNF77rs146879198G/A0.00140.00530.8625 <0.0001 <0.0001 <0.00010.02450.9643
ADGRL3rs192210727G/T0.01060.02620.6918 <0.0001 <0.0001 <0.00010.05010.5669
ALDH2rs671G/A0.00170.02310.0253 <0.00010.00020.02970.0479 <0.0001
ACAD10rs11066015G/A0.00140.03000.0261 <0.00010.00020.03760.0578 <0.0001
BRAPrs3782886A/G0.00120.06750.0233 <0.00010.00080.03290.0548 <0.0001
HECTD4rs11066280T/A0.00520.21610.0093 <0.00010.00020.01720.1293 <0.0001
HECTD4rs2074356C/T0.02450.21480.0042 <0.00010.00030.02180.2238 <0.0001
OR52E4rs11823828T/G0.0429 <0.00010.09170.23410.83120.01610.07680.1053
rs12229654T/G0.02000.43090.0996 <0.00010.02900.06330.2134 <0.0001
NAA25rs12231744C/T0.09130.52140.01510.00020.21750.85860.35430.0055
ATXN2rs7969300T/C0.03680.63410.00010.00050.71890.36310.31550.0047
GOSR2rs1052586T/C0.51750.08430.31620.07180.17790.65250.43070.5766
TCHPrs74416240G/A0.42440.02790.06300.05940.01420.00370.34090.1148
GIT2rs925368A/G0.34920.05150.03370.08730.01050.00200.43990.0998
rs2523644A/G0.44640.25810.56510.07970.06740.03970.42030.3215
OAS3rs2072134C/T0.64770.09360.19020.00040.00300.18620.25800.0185
CCDC63rs10774610T/C0.35730.21330.26230.00300.01020.00550.5909 <0.0001
rs2596548G/T0.27490.02680.26150.21380.32230.07510.90650.2475
CCDC63rs10849915T/C0.33860.22560.28910.00250.01350.00600.5869 <0.0001

[i] Data are P-values. The association between genotypes of each SNP and intermediate phenotypes was examined with Pearson's χ2 test. Based on Bonferroni's correction, P<0.0003 (0.05/200) were considered statistically significant and are shown in bold. SNP, single nucleotide polymorphism; DM, diabetes mellitus; hyper-TG, hypertriglyceridemia; hypo-HDL, hypo-high density lipoprotein-cholesterolemia; hyper-LDL, hyper-low density lipoprotein-cholesterolemia; CKD, chronic kidney disease.

Association between SNPs associated with hypertension and systolic or diastolic BP

The present study examined the associations between genotypes for the 11 SNPs associated with hypertension and systolic or diastolic BP by one-way ANOVA (Table XI). All 11 SNPs were significantly [P<0.0023 (0.05/22)] associated with systolic BP, whereas only rs2071593 of NFKBIL1 was associated with diastolic BP.

Table XI

Association between SNPs associated with hypertension and systolic or diastolic BP.

Table XI

Association between SNPs associated with hypertension and systolic or diastolic BP.

GeneSNPSystolic BP (mmHg)P-valueDiastolic BP (mmHg)P-value
MARCH1rs61734696G/TGGGT <0.0001GGGT0.0239
128±24119±1677±1475±12
PLCB2rs200787930 C/TCCCT <0.0001CCCT0.0482
128±24119±1677±1475±12
COL6A3rs146092501 C/TCCCT <0.0001CCCT0.0584
128±24119±1577±1475±12
COL6A5rs200982668G/AGGGA <0.0001GGGA0.0341
128±24119±1677±1475±12
MOB3Crs139537100 C/TCCCT <0.0001CCCT0.0359
128±24119±1677±1475±12
VPS33Brs199921354 C/TCCCT <0.0001CCCT0.0348
128±24119±1677±1475±12
NCR3rs2515920T/ATTTAAA <0.0001TTTAAA0.0062
128±24125±22126±2377±1476±1376±13
TMOD4rs115287176G/AGGGA <0.0001GGGA0.0793
128±24119±1677±1475±12
CXCL8rs188378669G/TGGGT <0.0001GGGT0.0703
128±24119±1677±1475±12
rs769177G/AGGGAAA0.0001GGGAAA0.0094
128±24125±22126±2377±1476±1376±13
NFKBIL1rs2071593C/TCCCTTT <0.0001CCCTTT0.0012
128±24125±22127±2377±1476±1376±13

[i] Data are presented as the mean ± standard deviation and were compared among genotypes by one-way analysis of variance. Based on Bonferroni's correction, P<0.0023 (0.05/22) was considered statistically significant and are shown in bold. SNP, single nucleotide polymorphism; BP, blood pressure.

Association between SNPs associated with CKD and the serum concentration of creatinine or eGFR

The present study examined the associations between the 11 SNPs associated with CKD and the serum concentration of creatinine or eGFR by one-way ANOVA (Table XII). All 11 SNPs were significantly [P<0.0023 (0.05/22)] associated with eGFR, however, none were associated with the serum concentration of creatinine.

Table XII

Association between SNPs associated with chronic kidney disease and the serum concentration of creatinine and eGFR.

Table XII

Association between SNPs associated with chronic kidney disease and the serum concentration of creatinine and eGFR.

GeneSNPSerum creatinine (µmol/l)P-valueeGFR (ml min−1 1.73 m−2)P-value
COL6A5rs200982668G/AGGGA0.0281GGGA <0.0001
93.6±138.561.8±34.479.3±32.995.2±20.2
MARCH1rs61734696G/TGGGT0.0352GGGT <0.0001
93.5±138.362.4±35.079.3±32.994.7±20.3
MUC17rs78010183A/TAAAT0.0068AAAT <0.0001
93.8±139.059.2±18.879.4±33.093.0±19.9
MOB3Crs139537100 C/TCCCT0.0374CCCT <0.0001
93.5±138.462.4±35.579.3±32.994.9±20.8
PLCB2rs200787930 C/TCCCT0.0380CCCT <0.0001
93.5±138.362.5±35.479.3±32.994.8±20.8
CXCL8rs188378669G/TGGGT0.0397GGGT <0.0001
93.5±138.362.6±35.679.3±32.995.0±20.9
VPS33Brs199921354 C/TCCCT0.0406CCCT <0.0001
93.5±138.462.6±35.879.3±32.994.8±21.1
TMOD4rs115287176G/AGGGA0.0426GGGA <0.0001
93.4±138.262.8±35.879.4±32.994.6±20.9
ADGRL3rs192210727G/TGGGTTT0.1294GGGTTT0.0002
93.2±137.563.0±35.864.579.3±32.994.4±21.494.3
ZNF77rs146879198G/AGGGA0.0491GGGA <0.0001
93.4±138.163.1±36.479.4±32.994.3±21.1
COL6A3rs146092501 C/TCCCT0.0466CCCT <0.0001
93.4±138.262.8±36.479.4±32.994.6±21.2

[i] Data are presented as the mean ± standard deviation and were compared among genotypes by one-way analysis of variance. Based on Bonferroni's correction, P<0.0023 (0.05/22) was considered statistically significant and are shown in bold. SNP, single nucleotide polymorphism; eGFR, estimated glomerular filtration rate.

Linkage disequilibrium analyses

The linkage disequilibrium (LD) among SNPs associated with MI, hypertension, or CKD was examined. For the MI study, rs192210727 of ADGRL3, rs188378669 of CXCL8, and rs61734696 of MARCH1 were in LD [square of the correlation coefficient (r2), 0.907-0.972]. Two SNPs (rs2596548, rs2523644) at chromosome 6p21.3 were also in LD (r2=0.624). LD plots of the 13 SNPs located at chromosome 12q24.11 to 12q24.13 are shown in Fig. 1. There was significant LD between rs74416240 of TCHP and rs925368 of GIT2 (r2=0.937); between rs10849915 and rs10774610 of CCDC63 (r2=0.992); and among rs3782886 of BRAP, rs11066015 of ACAD10, rs671 of ALDH2, and rs2074356 and rs11066280 of HECTD4 (r2=0.814-0.994). A significant LD (r2=0.994) was also apparent between rs200787930 of PLCB2 and rs199921354 of VPS33B.

For the hypertension study, there was significant LD between rs139537100 of MOB3C and rs115287176 of TMOD4 (r2=0.984); between rs188378669 of CXCL8 and rs61734696 of MARCH1 (r2=0.972); among rs2071593 of NFKBIL1, rs769177 at 6p21.3, and rs2515920 of NCR3 (r2=0.883-0.980); and between rs200787930 of PLCB2 and rs199921354 of VPS33B (r2=0.994).

For the CKD study, complete LD (r2=1.000) was apparent between rs139537100 of MOB3C and rs115287176 of TMOD4, and between rs200787930 of PLCB2 and rs199921354 of VPS33B. Significant LD (r2=0.905-0.952) was also observed among rs192210727 of ADGRL3, rs188378669 of CXCL8, and rs61734696 of MARCH1.

Associations between the genes, chromosomal loci and SNPs identified in the present study and the phenotypes examined in previous GWASs

The present study examined the association of the genes, chromosomal loci and SNPs identified in the present study with the phenotypes previously investigated by GWASs available in the Genome-Wide Repository of Associations Between SNPs and Phenotypes Search database (https://grasp.nhlbi.nih.gov/Search.aspx) developed by the Information Technology and Applications Center at the National Center for Biotechnology Information (National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD, USA) (60,61).

For the MI study, chromosomal region 6p21.3, ATXN2, BRAP, ACAD10, ALDH2, NAA25, and HECTD4 were previously shown to be associated with MI or CAD (Table XIII). The remaining 15 genes or loci identified in the present study have not been previously found to be associated with MI or CAD, although GOSR2 was previously found to be associated with systolic BP. For the hypertension study, none of the genes, chromosomal loci, or SNPs identified were found to be associated with systolic or diastolic BP or with hypertension in previous GWASs, although 6p21.3 was previously shown to be associated with CAD (Table XIV). For the CKD study, none of the genes or SNPs identified were shown to be associated with renal function-related traits or CKD in previous GWASs (Table XV).

Table XIII

Association between genes, chromosomal loci, and SNPs associated with MI in the present study and previously examined cardiovascular disease-related phenotypes.

Table XIII

Association between genes, chromosomal loci, and SNPs associated with MI in the present study and previously examined cardiovascular disease-related phenotypes.

Gene/chr. locusSNPChr.PositionPreviously examined phenotypes
TMOD4rs1152871761151170961None
COL6A3rs1460925012237371861None
ADGRL3rs192210727461909615None
CXCL8rs188378669473741568Total cholesterol (23063622)
MARCH1rs617346964164197303Adiponectin concentrations (20887962)
6p21.3rs2596548631362769MI (21971053)
rs252364431374707
OR52E4rs11823828115884973None
TCHPrs7441624012109904793None
GIT2rs92536812109953174None
CCDC63rs1084991512110895818None
rs10774610110902439
12q24.1rs1222965412110976657HDL-cholesterol (21909109)
ATXN2rs796930012111555908MI (19820697, 19198610), CAD (19820697, 23202125, 21060863)
BRAPrs378288612111672685MI (21971053, 19820697, 21107343), CAD (21971053, 21572416, 23202125, 19820697, 23364394, 21060863)
ACAD10rs1106601512111730205CAD (23364394, 23202125)
ALDH2rs67112111803962MI (21971053), CAD (21971053, 21572416, 23202125)
NAA25rs1223174412112039251MI (19820697, 23202125), CAD (21971053, 19820697)
HECTD4rs207435612112207597MI (19820697), CAD (21971053, 21572416, 22751097, 19820697, 23364394, 23202125)
rs11066280112379979
OAS3rs207213412112971371HDL-cholesterol (21909109)
PLCB2rs2007879301540289298Triglycerides (23063622)
VPS33Brs1999213541591013841Type 2 diabetes (22885922)
GOSR2rs10525861746941097Systolic BP (21909110, 21909115)
ZNF77rs146879198192934109None

[i] Data were obtained from the Genome-Wide Repository of Associations Between SNPs and Phenotypes Search database (https://grasp.nhlbi.nih.gov/Search.aspx) with P<1.0×10−6. Numbers in parentheses are PubMed IDs. SNP, single nucleotide polymorphism; MI, myocardial infarction; Chr., chromosome; HDL, high density lipoprotein; CAD, coronary artery disease; BP, blood pressure.

Table XIV

Association between genes, chromosomal loci, and SNPs associated with hypertension in the present study and previously examined cardiovascular disease-related phenotypes.

Table XIV

Association between genes, chromosomal loci, and SNPs associated with hypertension in the present study and previously examined cardiovascular disease-related phenotypes.

Gene/chr. locusSNPChr.PositionPreviously examined phenotypes
MOB3Crs139537100146615006None
TMOD4rs1152871761151170961None
COL6A3rs1460925012237371861None
COL6A5rs2009826683130470894None
CXCL8rs188378669473741568Total cholesterol (23063622)
MARCH1rs617346964164197303Adiponectin concentrations (20887962)
NFKBIL1rs2071593631545022Total cholesterol (20686565)
6p21.3rs769177631579834Total cholesterol (20686565), triglycerides (20686565), CAD (21971053)
NCR3rs2515920631594838Total cholesterol (20686565)
PLCB2rs2007879301540289298Triglycerides (23063622)
VPS33Brs1999213541591013841Type 2 diabetes (22885922)

[i] Data were obtained from the Genome-Wide Repository of Associations Between SNPs and Phenotypes Search database (https://grasp.nhlbi.nih.gov/Search.aspx) with a P<1.0×10−6. Numbers in parentheses are PubMed IDs. SNP, single nucleotide polymorphism; Chr., chromosome; CAD, coronary artery disease.

Table XV

Association between genes and SNPs associated with chronic kidney disease in the present study and previously examined cardiovascular and renal phenotypes.

Table XV

Association between genes and SNPs associated with chronic kidney disease in the present study and previously examined cardiovascular and renal phenotypes.

GeneSNPChr.PositionPreviously examined phenotypes
MOB3Crs139537100146615006None
TMOD4rs1152871761151170961None
COL6A3rs1460925012237371861None
COL6A5rs2009826683130470894None
ADGRL3rs192210727461909615None
CXCL8rs188378669473741568Total cholesterol (23063622)
MARCH1rs617346964164197303Adiponectin concentrations (20887962)
MUC17rs780101837101035329None
PLCB2rs2007879301540289298Triglycerides (23063622)
VPS33Brs1999213541591013841Type 2 diabetes (22885922)
ZNF77rs146879198192934109None

[i] Data were obtained from the Genome-Wide Repository of Associations Between SNPs and Phenotypes Search database (https://grasp.nhlbi.nih.gov/Search.aspx) with a P<1.0×10−6. Numbers in parentheses are PubMed IDs. SNP, single nucleotide polymorphism; Chr., chromosome.

Gene Ontology (GO) analysis of genes identified in the present study

The potential biological functions of the 18 genes (MOB3C, TMOD4, COL6A3, ADGRL3, CXCL8, MARCH1, OR52E4, TCHP, GIT2, CCDC63, OAS3, PLCB2, VPS33B, ZNF77, COL6A5, NFKBIL1, NCR3, MUC17) identified in the present study were examined with the use of the QuickGO database of Gene Ontology and GO Annotations (https://www.ebi.ac.uk/QuickGO; European Bioinformatics Institute, European Molecular Biology Laboratory, Hinxton, UK) (62,63). The GO analysis predicted various biological functions for the 18 genes, including those associated with muscle contraction (TMOD4), cell adhesion (COL6A3, COL6A5), G protein-coupled receptor signaling (ADGRL3, OR52E4, PLCB2), the inflammatory response (CXCL8), the immune response (MARCH1, OAS3 NCR3), apoptosis (TCHP), GTPase activity (GIT2), protein transport (VPS33B), NF-κB signaling (NFKBIL1), and lectin receptor signaling (MUC17) (Table XVI).

Table XVI

Gene Ontology analysis of the 18 genes identified in the present study.

Table XVI

Gene Ontology analysis of the 18 genes identified in the present study.

GeneFunctionBiological process
MOB3CProtein binding, metal ion binding Uncharacterized
TMOD4Actin binding, tropomyosin bindingMuscle contraction, actin filament organization, myofibril assembly, pointed-end actin filament capping
COL6A3Serine-type endopeptidase inhibitor activity, extracellular matrix structural constituent conferring tensile strength, collagen-containing extracellular matrixCell adhesion, muscle organ development, negative regulation of endopeptidase activity, extracellular matrix organization
ADGRL3G protein-coupled receptor activity, calcium ion binding, protein binding, carbohydrate bindingNeuron migration, signal transduction, G protein-coupled receptor signaling pathway, positive regulation of synapse assembly, cell-cell adhesion via plasma-membrane adhesion molecules
CXCL8Cytokine activity, interleukin-8 receptor binding, protein bindingPositive regulation of angiogenesis, induction of positive chemotaxis, defense response, inflammatory response, immune response, cell cycle arrest, negative regulation of G protein-coupled receptor protein signaling pathway, negative regulation of cell proliferation, cytokine-mediated signaling pathway, calcium-mediated signaling, regulation of cell adhesion, positive regulation of neutrophil chemotaxis, receptor internalization, response to endoplasmic reticulum stress, intracellular signal transduction, PKR-like endoplasmic reticulum kinase-mediated unfolded protein response, neutrophil activation, cellular responses to fibroblast growth factor stimulus, lipopolysaccharide, interleukin-1, and tumor necrosis factor
MARCH1Ubiquitin-protein transferase activity, zinc ion binding, MHC protein binding, ubiquitin protein ligase activityProtein polyubiquitination, immune response, antigen processing and presentation of peptide antigen via MHC class II
OR52E4G protein-coupled receptor activity, olfactory receptor activitySignal transduction, G protein-coupled receptor signaling pathway, sensory perception of smell, response to stimulus, detection of chemical stimulus involved in sensory perception of smell
TCHPProtein bindingApoptotic process, negative regulation of cell growth, negative regulation of cilium assembly
GIT2GTPase activator, kinase activity, metal ion bindingPhosphorylation, positive regulation of GTPase activity
CCDC63Cell differentiationSpermatogenesis, spermatid development, cilium movement, outer dynein arm assembly
OAS3Nucleotide binding, 2′-5′-oligoadenylate synthetase activity, RNA binding, protein binding, ATP binding, nucleotidyltransferase activity, metal ion bindingImmune response, nucleobase-containing compound metabolic process, defense response to virus, negative regulation of viral genome replication, interferon-g-mediated signaling pathway, regulation of ribonuclease activity
PLCB2 Phosphatidylinositol phospholipase C activity, calcium ion binding, phosphoric diester hydrolase activityPhospholipid metabolic process, intracellular signal transduction, G protein-coupled receptor signaling pathway, activation of phospholipase C activity, Wnt signaling pathway, calcium-modulating pathway, inositol phosphate metabolic process
VPS33BProtein binding, AP-3 adaptor complex, clathrin complexVesicle docking involved in exocytosis, endosome organization, protein transport, vesicle-mediated transport, peptidyllysine hydroxylation, melanosome localization, lysosome localization, collagen metabolic process, membrane fusion, platelet a granule organization
ZNF77RNA polymerase II transcription factor activity, sequence-specific DNA binding, metal ion bindingRegulation of DNA-templated transcription, regulation of transcription by RNA polymerase II
COL6A5Collagen-containing extracellular matrix, protein binding, extracellular matrix structural constituent conferring tensile strengthCell adhesion
NFKBIL1Protein bindingInhibitor of NF-κB kinase/NF-κB signaling, cytoplasmic sequestering of transcription factor, negative regulation of lipopolysaccharide-mediated signaling pathway, NF-κB transcription factor activity, tumor necrosis factor production, Toll-like receptor signaling pathway
NCR3Signaling receptor binding, protein binding, protein homodimerization activity, cell adhesion molecule bindingRegulation of immune response, homophilic and heterophilic cell-cell adhesion via plasma-membrane adhesion molecules, cell recognition, positive regulation of natural killer cell-mediated cytotoxicity, susceptibility to natural killer cell-mediated cytotoxicity and T cell-mediated cytotoxicity
MUC17Protein binding, PDZ domain binding, extracellular matrix constituent, lubricant activityStimulatory C-type lectin receptor signaling pathway, O-glycan processing, cellular homeostasis

[i] Data for predicted functions and biological processes for the genes were obtained from the database of Gene Ontology and GO Annotations (QuickGO; https://www.ebi.ac.uk/QuickGO). GO, Gene Ontology; MHC, major histocompatibility complex; NF-κB, nuclear factor-κB.

Network analysis of genes identified in the present study

Network analysis of the 18 genes identified in the present study was performed to predict functional gene-gene interactions with the use of the GeneMANIA Cytoscape plugin (http://apps.cytoscape.org/apps/genemania; Donnelly Centre for Cellular and Biomolecular Research, University of Toronto, Toronto, ON, Canada) (64-66) and Cytoscape v3.4.0 software (http://www.cytoscape.org; The Cytoscape Consortium, San Diego, CA, USA) (67). Three sets of 50 genes were selected from the DisGeNET database (http://www.disgenet.org/web/DisGeNET; Integrative Biomedical Informatics Group, Research Programme on Biomedical Informatics, Barcelona Biomedical Research Park, Barcelona, Spain) (68,69) according to the rank order (high to low) of scores for association with MI, hypertension, or CKD (Fig. 2). The network analysis revealed that the 13 (Fig. 2A), 10 (Fig. 2B), and 11 (Fig. 2C) genes found to be associated with MI, hypertension, and CKD, respectively, in the present study have potential direct or indirect interactions with the corresponding sets of 50 genes previously shown to be associated with these conditions.

Figure 2

(A) Network analysis of the 13 genes (TMOD4, COL6A3, ADGRL3, CXCL8, MARCH1, OR52E4, TCHP, GIT2, CCDC63, OAS3, PLCB2, VPS33B, ZNF77) associated with MI in the present study (red circles) was performed to predict functional gene-gene interactions with the use of the GeneMANIA Cytoscape plugin (http://apps.cytoscape.org/apps/genemania) and Cytoscape v3.4.0 software (http://www.cytoscape.org). A total of 50 genes (green circles) selected from the DisGeNET database (http://www.disgenet.org/web/DisGeNET) according to the rank order (high to low) of scores for association with MI were applied to the analysis. Interactions between red circles or between red and green circles are indicated by bold lines. Molecules represented by gray circles are putative mediators of interactions between genes. (B) Network analysis of the 10 genes (MOB3C, TMOD4, COL6A3, COL6A5, CXCL8, MARCH1, NFKBIL1, NCR3, PLCB2, VPS33B) associated with hypertension in the present study (red circles) was performed in the same manner; 50 genes (green circles) selected from the DisGeNET database according to the rank order of scores for association with hypertension were applied to the analysis. MI, myocardial infarction. (C) Network analysis of the 11 genes (MOB3C, TMOD4, COL6A3, COL6A5, ADGRL3, CXCL8, MARCH1, MUC17, PLCB2, VPS33B, ZNF77) associated with CKD in the present study (red circles) was performed in the same manner; 50 genes (green circles) selected from the DisGeNET database according to the rank order of scores for association with CKD were applied to the analysis. CKD, chronic kidney disease.

Discussion

As the prevalence of cardiovascular and renal diseases is increasing and is therefore a key public health concern (1,20,31,32), the identification of genetic variants that confer susceptibility to MI, hypertension, and CKD is important to prevent these conditions. In the present study, EWASs for MI, hypertension, and CKD were performed in subjects with early-onset forms of these conditions, with genetic factors likely to be more important in these individuals compared with those with late-onset forms.

In the MI study, it was found that 25 SNPs of 20 genes and two chromosomal loci were significantly associated with early-onset MI. Among these genes and loci, 6p21.3 (16), ATXN2 (9), BRAP (16), ACAD10 (9), ALDH2 (9), NAA25 (9), and HECTD4 (9) have previously been shown to be associated with MI or CAD. Therefore, 16 SNPs of 14 genes and one locus were newly identified. However, significant LD was detected among three SNPs in ADGRL3, CXCL8, and MARCH1; between two SNPs in TCHP and GIT2; between two SNPs in CCDC63; and between two SNPs in PLCB2 and VPS33B. Therefore, 11 novel loci (TMOD4, COL6A3, ADGRL3-CXCL8-MARCH1, OR52E4, TCHP-GIT2, CCDC63, 12q24.1, OAS3, PLCB2-VPS33B, GOSR2, ZNF77) were identified that confer susceptibility to MI.

In the hypertension study, 11 SNPs of 10 genes and one chromosomal locus were significantly associated with early-onset hypertension. None of these genes or loci have been found to be associated with systolic or diastolic BP or with hypertension in previous GWASs. However, there was significant LD between two SNPs in MOB3C and TMOD4; between two SNPs in CXCL8 and MARCH1; among three SNPs of NFKBIL1, 6p21.3, and NCR3; and between two SNPs in PLCB2 and VPS33B. Therefore six novel loci (MOB3C-TMOD4, COL6A3, COL6A5, CXCL8-MARCH1, NFKBIL1-6p21.3-NCR3, PLCB2-VPS33B) were identified that confer susceptibility to hypertension.

In the CKD study, 11 SNPs in 11 genes were significantly associated with early-onset CKD. None of these genes have been shown to be associated with renal function-related traits or CKD in previous GWASs. However, there was significant LD between two SNPs in MOB3C and TMOD4; between two SNPs in PLCB2 and VPS33B; and among three SNPs in ADGRL3, CXCL8, and MARCH1. Therefore, seven new loci (MOB3C-TMOD4, COL6A3, COL6A5, ADGRL3-CXCL8-MARCH1, MUC17, PLCB2-VPS33B, ZNF77) were identified that confer susceptibility to CKD.

Furthermore, rs115287176 of TMOD4, rs146092501 of COL6A3, rs188378669 of CXCL8, rs61734696 of MARCH1, rs200787930 of PLCB2, and rs199921354 of VSP33B were significantly associated with all three diseases (MI, hypertension, and CKD); rs192210727 of ADGRL3 and rs146879198 of ZNF77 were associated with both MI and CKD; and rs200982668 of COL6A5 and rs139537100 of MOB3C were associated with both hypertension and CKD. In addition, GOSR2, which was associated with MI, has previously been shown to be associated with systolic BP (23); and 6p21.3, which was associated with hypertension, has previously been shown to be associated with CAD (16). Therefore, 13 loci (MOB3C-TMOD4, COL6A3, ADGRL3-CXCL8-MARCH1, OR52E4, TCHP- GIT2, CCDC63, 12q24.1, OAS3, PLCB2-VPS33B, ZNF77, COL6A5, NFKBIL1-NCR3, MUC17) were newly identified that confer susceptibility to MI, hypertension, or CKD. Genes, chromosomal loci, and SNPs identified in the present study, in particular, those significantly associated with two or three diseases, may prove to be informative for clinical practice.

In the MI study, associations of CXCL8 to plasma total cholesterol (70), of MARCH1 to serum adiponectin concentration (71), of 12q24.1 and OAS3 to plasma HDL-cholesterol (72), of PLCB2 to plasma triglycerides (70), of VPS33B to type 2 diabetes mellitus (73), and of GOSR2 to systolic BP (23) have been shown in previous GWASs. In the present study, eight loci (TMOD4, COL6A3, ADGRL3-CXCL8-MARCH1, OR52E4, CCDC63, 12q24.1, PLCB2-VPS33B, ZNF77) newly associated with MI were significantly associated with one or more intermediate phenotypes; in particular, five loci (TMOD4, COL6A3, ADGRL3-CXCL8-MARCH1, PLCB2-VPS33B, ZNF77) associated with two or three diseases were associated with three or four intermediate phenotypes. As these intermediate phenotypes are important risk factors for MI, the association between these loci and MI may be attributable, at least in part, to their effects on intermediate phenotypes. By contrast, three loci (TCHP-GIT2, OAS3, and GOSR2) associated with MI were not associated with intermediate phenotypes. The molecular mechanisms underlying the association between these loci and MI remain to be elucidated.

In the hypertension study, although the six loci (MOB3C-TMOD4, COL6A3, COL6A5, CXCL8-MARCH1, NFKBIL1-6p21.3-NCR3, PLCB2-VPS33B) identified in the present study were all significantly associated with systolic BP, they have not been shown to be associated with BP or hypertension by previous GWASs. The associations of NFKBIL1 and NCR3 to plasma total cholesterol (74), of 6p21.3 to plasma total cholesterol and triglycerides (74), and of CXCL8 (70), MARCH1 (71), PLCB2 (70), and VPS33B (73) to lipid or glucose metabolism have been shown in previous GWASs. The effects of these loci on lipid or glucose metabolism may influence BP or hypertension, although the functional relevance of these loci to the pathogenesis of hypertension remains to be fully elucidated.

In the CKD study, although the seven loci (MOB3C-TMOD4, COL6A3, COL6A5, ADGRL3-CXCL8-MARCH1, MUC17, PLCB2-VPS33B, ZNF77) identified in the present study were all significantly associated with eGFR, they have not been shown to be associated with renal function-related traits or CKD in previous GWASs. The associations of CXCL8 (70), MARCH1 (71), PLCB2 (70), and VPS33B (73) to lipid or glucose metabolism have been demonstrated in previous GWASs. The effects of these loci on lipid or glucose metabolism may contribute, at least in part, to renal function or the development of CKD, although the molecular mechanisms underlying the association between these loci and CKD remain to be elucidated.

Previous GWASs have identified potential biological pathways underlying the associations between genetic loci and CAD (75,76). Network analysis of functional gene-gene interactions may be informative with regard to the clarification of biological processes underlying CAD and to the identification of therapeutic targets for this condition (77). Therefore the present study performed GO and network analyses to predict biological processes associated with the identified genes, and the interactions between these genes and those previously shown to be associated with MI, hypertension, or CKD.

GO analysis of the 13 genes associated with MI suggested that the functions of TMOD4 (actin filament organization); COL6A3 (cell adhesion and extracellular matrix organization); ADGRL3, OR52E4, and PLCB2 (G protein-coupled receptor signaling); CXCL8 (angiogenesis and inflammatory response); MARCH1 and OAS3 (immune response); TCHP (apoptotic process); GIT2 (regulation of GTPase activity); and VPS33B (vesicle-mediated transport and platelet α-granule organization) may be associated with the pathogenesis of MI. For the 10 genes associated with hypertension, the functions of TMOD4, COL6A3, COL6A5 (cell adhesion), CXCL8, MARCH1, NFKBIL1 (NF-κB signaling), NCR3 (immune response), PLCB2, and VPS33B may be associated with the pathogenesis of hypertension. Among the 11 genes associated with CKD, the functions of COL6A3, COL6A5, ADGRL3, CXCL8, MARCH1, MUC17 (cellular homeostasis), PLCB2, and VPS33B may be associated with the pathogenesis of CKD. Network analysis showed that the 13, 10, and 11 genes found to be associated with MI, hypertension, and CKD, respectively, in the present study had direct or indirect interactions with the corresponding sets of 50 genes selected from the DisGeNET database (68,69). However, the underlying molecular mechanisms of these interactions remain to be elucidated.

In our previous studies, it was shown that nine, nine, and four SNPs were associated with MI (P<0.01), hypertension (P<0.01), and CKD (P<0.05), respectively, as determined by multivariable logistic regression analysis with adjustment for covariates following an initial EWAS screening of allele frequencies among subjects with early-onset and late-onset forms of these conditions (50-52). The associations between five of the nine SNPs [rs202103723 (P=0.0107), rs188212047 (P=0.0039), rs1265110 (P=0.0004), rs9258102 (P=0.0373), rs439121 (P=0.0063)] and MI were replicated (P<0.05) in the present study. The associations between six of the nine SNPs [rs150854849 (P=0.0123), rs2867125 (P=0.0111), rs202069030 (P=6.40×10−10), rs139012426 (P=5.50×10−6), rs12229654 (P=0.0003), rs201633733 (P=0.0356)] and hypertension were replicated (P<0.05) in the present study. By contrast, none of the associations between the four SNPs and CKD were replicated in the present study, indicative of the importance of age in the development of CKD. These results suggested that genetic variants associated with MI, hypertension, and CKD differ, in part, between individuals with early-onset and late-onset disease.

There were several limitations to the present study: i) As that the results were not replicated, their validation is necessary in independent study populations or in other ethnicities; 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 MI, hypertension, or CKD; iii) the functional relevance of identified SNPs to the pathogenesis of MI, hypertension, and CKD remains to be elucidated. Due to the lack of experiments for functional analyses, the association of the SNPs identified in the present study with MI, hypertension, and CKD requires careful interpretation.

In conclusion, the present study identified 13 loci (MOB3C-TMOD4, COL6A3, ADGRL3-CXCL8-MARCH1, OR52E4, TCHP- GIT2, CCDC63, 12q24.1, OAS3, PLCB2-VPS33B, ZNF77, COL6A5, NFKBIL1-NCR3, MUC17) that confer susceptibility to MI, hypertension, or CKD. Determination of the genotypes for the SNPs at these loci may provide informative for assessment of the genetic risk for MI, hypertension, and CKD.

Funding

The present study was supported by the CREST (grant no. JPMJCR1302), Japan Science and Technology Agency to YYamada, JS, and IT).

Availability of data and materials

All data underlying the findings described are available upon request from the corresponding author.

Authors' contributions

YYamada 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. YYasukochi, IT, and JS contributed to analysis and interpretation of the data and to revision of the manuscript. All authors read and approved the final 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 participants.

Patient consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Acknowledgments

Not applicable.

References

1 

Benjamin EJ, Virani SS, Callaway CW, Chamberlain AM, Chang AR, Cheng S, Chiuve SE, Cushman M, Delling FN, Deo R, et al: Heart disease and stroke statistics-2018 update: A report from the american heart association. Circulation. 137:e67–e492. 2018. View Article : Google Scholar : PubMed/NCBI

2 

McPherson R and Tybjaerg-Hansen A: Genetics of coronary artery disease. Circ Res. 118:564–578. 2016. View Article : Google Scholar : PubMed/NCBI

3 

Peden JF and Farrall M: Thirty-five common variants for coronary artery disease: The fruits of much collaborative labour. Hum Mol Genet. 20:R198–R205. 2011. View Article : Google Scholar : PubMed/NCBI

4 

McPherson R, Pertsemlidis A, Kavaslar N, Stewart A, Roberts R, Cox DR, Hinds DA, Pennacchio LA, Tybjaerg-Hansen A, Folsom AR, et al: A common allele on chromosome 9 associated with coronary heart disease. Science. 316:1488–1491. 2007. View Article : Google Scholar : PubMed/NCBI

5 

Helgadottir A, Thorleifsson G, Manolescu A, Gretarsdottir S, Blondal T, Jonasdottir A, Jonasdottir A, Sigurdsson A, Baker A, Palsson A, et al: A common variant on chromosome 9p21 affects the risk of myocardial infarction. Science. 316:1491–1493. 2007. View Article : Google Scholar : PubMed/NCBI

6 

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

7 

Samani NJ, Erdmann J, Hall AS, Hengstenberg C, Mangino M, Mayer B, Dixon RJ, Meitinger T, Braund P, Wichmann HE, et al: Genomewide association analysis of coronary artery disease. N Engl J Med. 357:443–453. 2007. View Article : Google Scholar : PubMed/NCBI

8 

Schunkert H, König IR, Kathiresan S, Reilly MP, Assimes TL, Holm H, Preuss M, Stewart AF, Barbalic M, Gieger C, et al: Large-scale association analysis identifies 13 new susceptibility loci for coronary artery disease. Nat Genet. 43:333–338. 2011. View Article : Google Scholar : PubMed/NCBI

9 

CARDIoGRAMplusC4D Consortium; Deloukas P, Kanoni S, Willenborg C, Farrall M, Assimes TL, Thompson JR, Ingelsson E, Saleheen D, Erdmann J, et al: Large-scale association analysis identifies new risk loci for coronary artery disease. Nat Genet. 45:25–33. 2013. View Article : Google Scholar

10 

van der Harst P and Verweij N: Identification of 64 novel genetic loci provides an expanded view on the genetic architecture of coronary artery disease. Circ Res. 122:433–443. 2018. View Article : Google Scholar :

11 

Lettre G, Palmer CD, Young T, Ejebe KG, Allayee H, Benjamin EJ, Bennett F, Bowden DW, Chakravarti A, Dreisbach A, et al: Genome-wide association study of coronary heart disease and its risk factors in 8,090 african americans: The NHLBI CARe Project. PLoS Genet. 7:e10013002011. View Article : Google Scholar : PubMed/NCBI

12 

Wang F, Xu CQ, He Q, Cai JP, Li XC, Wang D, Xiong X, Liao YH, Zeng QT, Yang YZ, et al: Genome-wide association identifies a susceptibility locus for coronary artery disease in the chinese han population. Nat Genet. 43:345–349. 2011. View Article : Google Scholar : PubMed/NCBI

13 

Lu X, Wang L, Chen S, He L, Yang X, Shi Y, Cheng J, Zhang L, Gu CC, Huang J, et al: Genome-wide association study in han chinese identifies four new susceptibility loci for coronary artery disease. Nat Genet. 44:890–894. 2012. View Article : Google Scholar : PubMed/NCBI

14 

Nikpay M, Goel A, Won HH, Hall LM, Willenborg C, Kanoni S, Saleheen D, Kyriakou T, Nelson CP, Hopewell JC, et al: A comprehensive 1000 Genomes-based genome-wide association meta-analysis of coronary artery disease. Nat Genet. 47:1121–1130. 2015. View Article : Google Scholar : PubMed/NCBI

15 

Nelson CP, Goel A, Butterworth AS, Kanoni S, Webb TR, Marouli E, Zeng L, Ntalla I, Lai FY, Hopewell JC, et al: Association analyses based on false discovery rate implicate new loci for coronary artery disease. Nat Genet. 49:1385–1391. 2017. View Article : Google Scholar : PubMed/NCBI

16 

Takeuchi F, Yokota M, Yamamoto K, Nakashima E, Katsuya T, Asano H, Isono M, Nabika T, Sugiyama T, Fujioka A, et al: Genome-wide association study of coronary artery disease in the japanese. Eur J Hum Genet. 20:333–340. 2012. View Article : Google Scholar :

17 

Hirokawa M, Morita H, Tajima T, Takahashi A, Ashikawa K, Miya F, Shigemizu D, Ozaki K, Sakata Y, Nakatani D, et al: A genome-wide association study identifies PLCL2 and AP3D1-DOT1L-SF3A2 as new susceptibility loci for myocardial infarction in japanese. Eur J Hum Genet. 23:374–380. 2015. View Article : Google Scholar :

18 

Padmanabhan S, Caulfield M and Dominiczak AF: Genetic and molecular aspects of hypertension. Circ Res. 116:937–959. 2015. View Article : Google Scholar : PubMed/NCBI

19 

Byars SG, Ewbank D, Govindaraju DR and Stearns SC: Natural selection in a contemporary human population. Proc Natl Acad Sci USA. 107(Suppl 1): S1787–S1792. 2010. View Article : Google Scholar

20 

Rapsomaniki E, Timmis A, George J, Pujades-Rodriguez M, Shah AD, Denaxas S, White IR, Caulfield MJ, Deanfield JE, Smeeth L, et al: Blood pressure and incidence of twelve cardiovascular diseases: Lifetime risks, healthy life-years lost, and age-specific associations in 1.25 million people. Lancet. 383:1899–1911. 2014. View Article : Google Scholar : PubMed/NCBI

21 

Levy D, Ehret GB, Rice K, Verwoert GC, Launer LJ, Dehghan A, Glazer NL, Morrison AC, Johnson AD, Aspelund T, et al: Genome-wide association study of blood pressure and hypertension. Nat Genet. 41:677–687. 2009. View Article : Google Scholar : PubMed/NCBI

22 

Newton-Cheh C, Johnson T, Gateva V, Tobin MD, Bochud M, Coin L, Najjar SS, Zhao JH, Heath SC, Eyheramendy S, et al: Genome-wide association study identifies eight loci associated with blood pressure. Nat Genet. 41:666–676. 2009. View Article : Google Scholar : PubMed/NCBI

23 

International Consortium for Blood Pressure Genome-Wide Association Studies; Ehret GB, Munroe PB, Rice KM, Bochud M, Johnson AD, Chasman DI, Smith AV, Tobin MD, Verwoert GC, et al: Genetic variants in novel pathways influence blood pressure and cardiovascular disease risk. Nature. 478:103–109. 2011. View Article : Google Scholar : PubMed/NCBI

24 

Warren HR, Evangelou E, Cabrera CP, Gao H, Ren M, Mifsud B, Ntalla I, Surendran P, Liu C, Cook JP, et al: Genome-wide association analysis identifies novel blood pressure loci and offers biological insights into cardiovascular risk. Nat Genet. 49:403–415. 2017. View Article : Google Scholar : PubMed/NCBI

25 

Adeyemo A, Gerry N, Chen G, Herbert A, Doumatey A, Huang H, Zhou J, Lashley K, Chen Y, Christman M and Rotimi C: A genome-wide association study of hypertension and blood pressure in african americans. PLoS Genet. 5:e10005642009. View Article : Google Scholar : PubMed/NCBI

26 

Liang J, Le TH Edwards, DRV Tayo, BO Gaulton, KJ Smith, JA Lu Y, Jensen RA, Chen G, Yanek LR, et al: Single-trait and multi-trait genome-wide association analyses identify novel loci for blood pressure in african-ancestry populations. PLoS Genet. 13:e10067282017. View Article : Google Scholar : PubMed/NCBI

27 

Kato N, Takeuchi F, Tabara Y, Kelly TN, Go MJ, Sim X, Tay WT, Chen CH, Zhang Y, Yamamoto K, et al: Meta-analysis of genome-wide association studies identifies common variants associated with blood pressure variation in east asians. Nat Genet. 43:531–538. 2011. View Article : Google Scholar : PubMed/NCBI

28 

Kato N, Loh M, Takeuchi F, Verweij N, Wang X, Zhang W, Kelly TN, Saleheen D, Lehne B, Mateo Leach I, et al: Trans-ancestry genome-wide association study identifies 12 genetic loci influencing blood pressure and implicates a role for DNA methylation. Nat Genet. 47:1282–1293. 2015. View Article : Google Scholar : PubMed/NCBI

29 

Surendran P, Drenos F, Young R, Warren H, Cook JP, Manning AK, Grarup N, Sim X, Barnes DR, Witkowska K, et al: Trans-ancestry meta-analyses identify rare and common variants associated with blood pressure and hypertension. Nat Genet. 48:1151–1161. 2016. View Article : Google Scholar : PubMed/NCBI

30 

Kato N, Miyata T, Tabara Y, Katsuya T, Yanai K, Hanada H, Kamide K, Nakura J, Kohara K, Takeuchi F, et al: High-density association study and nomination of susceptibility genes for hypertension in the japanese national project. Hum Mol Genet. 17:617–627. 2008. View Article : Google Scholar

31 

Webster AC, Nagler EV, Morton RL and Masson P: Chronic kidney disease. Lancet. 389:1238–1252. 2017. View Article : Google Scholar

32 

Weiner DE, Tighiouart H, Amin MG, Stark PC, MacLeod B, Griffith JL, Salem DN, Levey AS and Sarnak MJ: Chronic kidney disease as a risk factor for cardiovascular disease and all-cause mortality: A pooled analysis of community-based studies. J Am Soc Nephrol. 15:1307–1315. 2004. View Article : Google Scholar : PubMed/NCBI

33 

James MT, Hemmelgarn BR and Tonelli M: Early recognition and prevention of chronic kidney disease. Lancet. 375:1296–1309. 2010. View Article : Google Scholar : PubMed/NCBI

34 

Yamagata K, Ishida K, Sairenchi T, Takahashi H, Ohba S, Shiigai T, Narita M and Koyama A: Risk factors for chronic kidney disease in a community-based population: A 10-year follow-up study. Kidney Int. 71:159–166. 2007. View Article : Google Scholar

35 

Vivante A and Hildebrandt F: Exploring the genetic basis of early-onset chronic kidney disease. Nat Rev Nephrol. 12:133–146. 2016. View Article : Google Scholar : PubMed/NCBI

36 

Bochud M, Elston RC, Maillard M, Bovet P, Schild L, Shamlaye C and Burnier M: Heritability of renal function in hypertensive families of African descent in the Seychelles (Indian Ocean). Kidney Int. 67:61–69. 2005. View Article : Google Scholar

37 

Köttgen A, Glazer NL, Dehghan A, Hwang SJ, Katz R, Li M, Yang Q, Gudnason V, Launer LJ, Harris TB, et al: Multiple loci associated with indices of renal function and chronic kidney disease. Nat Genet. 41:712–717. 2009. View Article : Google Scholar : PubMed/NCBI

38 

Chambers JC, Zhang W, Lord GM, van der Harst P, Lawlor DA, Sehmi JS, Gale DP, Wass MN, Ahmadi KR, Bakker SJ, et al: Genetic loci influencing kidney function and chronic kidney disease. Nat Genet. 42:373–375. 2010. View Article : Google Scholar : PubMed/NCBI

39 

Köttgen A, Pattaro C, Böger CA, Fuchsberger C, Olden M, Glazer NL, Parsa A, Gao X, Yang Q, Smith AV, et al: New loci associated with kidney function and chronic kidney disease. Nat Genet. 42:376–384. 2010. View Article : Google Scholar : PubMed/NCBI

40 

Böger CA, Gorski M, Li M, Hoffmann MM, Huang C, Yang Q, Teumer A, Krane V, O'Seaghdha CM, Kutalik Z, et al: Association of eGFR-related loci identified by GWAS with incident CKD and ESRD. PLoS Genet. 7:e10022922011. View Article : Google Scholar : PubMed/NCBI

41 

Pattaro C, Köttgen A, Teumer A, Garnaas M, Böger CA, Fuchsberger C, Olden M, Chen MH, Tin A, Taliun D, et al: Genome-wide association and functional follow-up reveals new loci for kidney function. PLoS Genet. 8:e10025842012. View Article : Google Scholar : PubMed/NCBI

42 

Gorski M, van der Most PJ, Teumer A, Chu AY, Li M, Mijatovic V, Nolte IM, Cocca M, Taliun D, Gomez F, et al: 1000 Genomes-based meta-analysis identifies 10 novel loci for kidney function. Sci Rep. 7:468352017. View Article : Google Scholar

43 

Liu CT, Garnaas MK, Tin A, Kottgen A, Franceschini N, Peralta CA, de Boer IH, Lu X, Atkinson E, Ding J, et al: Genetic association for renal traits among participants of African ancestry reveals new loci for renal function. PLoS Genet. 7:e10022642011. View Article : Google Scholar : PubMed/NCBI

44 

Okada Y, Sim X, Go MJ, Wu JY, Gu D, Takeuchi F, Takahashi A, Maeda S, Tsunoda T, Chen P, et al: Meta-analysis identifies multiple loci associated with kidney function-related traits in east asian populations. Nat Genet. 44:904–909. 2012. View Article : Google Scholar : PubMed/NCBI

45 

Pattaro C, Teumer A, Gorski M, Chu AY, Li M, Mijatovic V, Garnaas M, Tin A, Sorice R, Li Y, et al: Genetic associations at 53 loci highlight cell types and biological pathways relevant for kidney function. Nat Commun. 7:100232016. View Article : Google Scholar : PubMed/NCBI

46 

Marenberg ME, Risch N, Berkman LF, Floderus B and de Faire U: Genetic susceptibility to death from coronary heart disease in a study of twins. N Engl J Med. 330:1041–1046. 1994. View Article : Google Scholar : PubMed/NCBI

47 

Nora JJ, Lortscher RH, Spangler RD, Nora AH and Kimberling WJ: Genetic-epidemiologic study of early-onset ischemic heart disease. Circulation. 61:503–508. 1980. View Article : Google Scholar : PubMed/NCBI

48 

Roncaglioni MC, Santoro L, D'Avanzo B, Negri E, Nobili A, Ledda A, Pietropaolo F, Franzosi MG, La Vecchia C, Feruglio GA, et al: Role of family history in patients with myocardial infarction. An italian case-control study. GISSI-EFRIM investigators. Circulation. 85:2065–2072. 1992. View Article : Google Scholar : PubMed/NCBI

49 

Niiranen TJ, McCabe EL, Larson MG, Henglin M, Lakdawala NK, Vasan RS and Cheng S: Heritability and risks associated with early onset hypertension: Multigenerational, prospective analysis in the framingham heart study. BMJ. 357:j19492017. View Article : Google Scholar : PubMed/NCBI

50 

Yamada Y, Sakuma J, Takeuchi I, Yasukochi Y, Kato K, Oguri M, Fujimaki T, Horibe H, Muramatsu M, Sawabe M, et al: Identification of STXBP2 as a novel susceptibility locus for myocardial infarction in Japanese individuals by an exome-wide association study. Oncotarget. 8:33527–33535. 2017.PubMed/NCBI

51 

Yamada Y, Sakuma J, Takeuchi I, Yasukochi Y, Kato K, Oguri M, Fujimaki T, Horibe H, Muramatsu M, Sawabe M, et al: Identification of polymorphisms in 12q24.1, ACAD10, and BRAP as novel genetic determinants of blood pressure in Japanese by exome-wide association studies. Oncotarget. 8:43068–43079. 2017.PubMed/NCBI

52 

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.PubMed/NCBI

53 

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

54 

Matsuo S, Imai E, Horio M, Yasuda Y, Tomita K, Nitta K, Yamagata K, Tomino Y, Yokoyama H and Hishida A; Collaborators developing the Japanese equation for estimated GFR: Revised equations for estimated GFR from serum creatinine in japan. Am J Kidney Dis. 53:982–992. 2009. View Article : Google Scholar : PubMed/NCBI

55 

National Kidney Foundation: K/DOQI clinical practice guidelines for chronic kidney disease: Evaluation, classification, and stratification. Am J Kidney Dis. 39(Suppl 1): S1–S266. 2002.PubMed/NCBI

56 

Go AS, Chertow GM, Fan D, McCulloch CE and Hsu CY: Chronic kidney disease and the risks of death, cardiovascular events, and hospitalization. N Engl J Med. 351:1296–1305. 2004. View Article : Google Scholar

57 

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

58 

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

59 

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

60 

Leslie R, O'Donnell CJ and Johnson AD: GRASP: Analysis of genotype-phenotype results from 1390 genome-wide association studies and corresponding open access database. Bioinformatics. 30:i185–i194. 2014. View Article : Google Scholar : PubMed/NCBI

61 

Eicher JD, Landowski C, Stackhouse B, Sloan A, Chen W, Jensen N, Lien JP, Leslie R and Johnson AD: GRASP v2.0: An update on the genome-wide repository of associations between SNPs and phenotypes. Nucleic Acids Res. 43:D799–D804. 2015. View Article : Google Scholar :

62 

Binns D, Dimmer E, Huntley R, Barrell D, O'Donovan C and Apweiler R: QuickGO: A web-based tool for gene ontology searching. Bioinformatics. 25:3045–3046. 2009. View Article : Google Scholar : PubMed/NCBI

63 

Huntley RP, Binns D, Dimmer E, Barrell D, O'Donovan C and Apweiler R: QuickGO: A user tutorial for the web-based gene ontology browser. Database (Oxford). 2009:bap0102009. View Article : Google Scholar

64 

Warde-Farley D, Donaldson SL, Comes O, Zuberi K, Badrawi R, Chao P, Franz M, Grouios C, Kazi F, Lopes CT, et al: The GeneMANIA prediction server: Biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 38:W214–W220. 2010. View Article : Google Scholar : PubMed/NCBI

65 

Montojo J, Zuberi K, Rodriguez H, Kazi F, Wright G, Donaldson SL, Morris Q and Bader GD: GeneMANIA cytoscape plugin: Fast gene function predictions on the desktop. Bioinformatics. 26:2927–2928. 2010. View Article : Google Scholar : PubMed/NCBI

66 

Montojo J, Zuberi K, Rodriguez H, Bader GD and Morris Q: GeneMANIA: Fast gene network construction and function prediction for Cytoscape. F1000Res. 3:1532014.PubMed/NCBI

67 

Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B and Ideker T: Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 13:2498–2504. 2003. View Article : Google Scholar : PubMed/NCBI

68 

Piñero J, Queralt-Rosinach N, Bravo À, Deu-Pons J, Bauer-Mehren A, Baron M, Sanz F and Furlong LI: DisGeNET: A discovery platform for the dynamical exploration of human diseases and their genes. Database (Oxford). 2015. pp. bav0282015, View Article : Google Scholar

69 

Piñero J, Bravo À, Queralt-Rosinach N, Gutiérrez-Sacristán A, Deu-Pons J, Centeno E, García-García J, Sanz F and Furlong LI: DisGeNET: A comprehensive platform integrating information on human disease-associated genes and variants. Nucleic Acids Res. 45:D833–D839. 2017. View Article : Google Scholar :

70 

Asselbergs FW, Guo Y, van Iperen EP, Sivapalaratnam S, Tragante V, Lanktree MB, Lange LA, Almoguera B, Appelman YE, Barnard J, et al: Large-scale gene-centric meta-analysis across 32 studies identifies multiple lipid loci. Am J Hum Genet. 91:823–838. 2012. View Article : Google Scholar : PubMed/NCBI

71 

Jee SH, Sull JW, Lee JE, Shin C, Park J, Kimm H, Cho EY, Shin ES, Yun JE, Park JW, et al: Adiponectin concentrations: A genome-wide association study. Am J Hum Genet. 87:545–552. 2010. View Article : Google Scholar : PubMed/NCBI

72 

Kim YJ, Go MJ, Hu C, Hong CB, Kim YK, Lee JY, Hwang JY, Oh JH, Kim DJ, Kim NH, et al: 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

73 

Morris AP, Voight BF, Teslovich TM, Ferreira T, Segrè AV, Steinthorsdottir V, Strawbridge RJ, Khan H, Grallert H, Mahajan A, et al: 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

74 

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

75 

Klarin D, Zhu QM, Emdin CA, Chaffin M, Horner S, McMillan BJ, Leed A, Weale ME, Spencer CCA, Aguet F, et al: Genetic analysis in UK Biobank links insulin resistance and transendothelial migration pathways to coronary artery disease. Nat Genet. 49:1392–1397. 2017. View Article : Google Scholar : PubMed/NCBI

76 

Howson JMM, Zhao W, Barnes DR, Ho WK, Young R, Paul DS, Waite LL, Freitag DF, Fauman EB, Salfati EL, et al: Fifteen new risk loci for coronary artery disease highlight arterial-wall-specific mechanisms. Nat Genet. 49:1113–1119. 2017. View Article : Google Scholar : PubMed/NCBI

77 

Lempiäinen H, Brænne I, Michoel T, Tragante V, Vilne B, Webb TR, Kyriakou T, Eichner J, Zeng L, Willenborg C, et al: Network analysis of coronary artery disease risk genes elucidates disease mechanisms and druggable targets. Sci Rep. 8:34342018. View Article : Google Scholar : PubMed/NCBI

Related Articles

Journal Cover

November-2018
Volume 42 Issue 5

Print ISSN: 1107-3756
Online ISSN:1791-244X

Sign up for eToc alerts

Recommend to Library

Copy and paste a formatted citation
x
Spandidos Publications style
Yamada Y, Kato K, Oguri M, Horibe H, Fujimaki T, Yasukochi Y, Takeuchi I and Sakuma J: Identification of 13 novel susceptibility loci for early-onset myocardial infarction, hypertension, or chronic kidney disease. Int J Mol Med 42: 2415-2436, 2018
APA
Yamada, Y., Kato, K., Oguri, M., Horibe, H., Fujimaki, T., Yasukochi, Y. ... Sakuma, J. (2018). Identification of 13 novel susceptibility loci for early-onset myocardial infarction, hypertension, or chronic kidney disease. International Journal of Molecular Medicine, 42, 2415-2436. https://doi.org/10.3892/ijmm.2018.3852
MLA
Yamada, Y., Kato, K., Oguri, M., Horibe, H., Fujimaki, T., Yasukochi, Y., Takeuchi, I., Sakuma, J."Identification of 13 novel susceptibility loci for early-onset myocardial infarction, hypertension, or chronic kidney disease". International Journal of Molecular Medicine 42.5 (2018): 2415-2436.
Chicago
Yamada, Y., Kato, K., Oguri, M., Horibe, H., Fujimaki, T., Yasukochi, Y., Takeuchi, I., Sakuma, J."Identification of 13 novel susceptibility loci for early-onset myocardial infarction, hypertension, or chronic kidney disease". International Journal of Molecular Medicine 42, no. 5 (2018): 2415-2436. https://doi.org/10.3892/ijmm.2018.3852