Micro- and Macrovascular Complications in Diabetes Mellitus: Preclinical and Clinical StudiesView this Special Issue
Polymorphisms in the Glucagon-Like Peptide 1 Receptor (GLP-1R) Gene Are Associated with the Risk of Coronary Artery Disease in Chinese Han Patients with Type 2 Diabetes Mellitus: A Case-Control Study
Background. Glucagon-like peptide 1 (GLP-1) bestows protective effects upon the cardiovascular system through direct cardiovascular interactions or by improvements to metabolic function. Both these effects are thought to be at least partly mediated by the GLP-1 receptor (GLP-1R). This case-controlled study investigated whether polymorphisms in the GLP-1R gene affect the risk of cardiovascular disease in type 2 diabetic patients in the Chinese Han population. Methods. Eleven haplotype-tagging single nucleotide polymorphisms (SNPs), distributed across 22 kb of the 39 kb GLP-1R gene, were selected and genotyped in diabetic patients from a Chinese Han population. Patients were classified based on the severity of coronary artery stenosis. Coronary artery stenosis was ≥50% in 394 patients (coronary artery disease- (CAD-) positive group), and coronary artery stenosis was <50% in 217 patients (control group). Allele and genotype frequencies were compared between the two groups at all 11 SNPs. Results. When considered in recessive inheritance mode, patients with the GG genotype at rs4714210 had a lower CAD risk than patients with other genotypes ( = 0.442, 95% CI = 0.258–0.757, ), even when other known CAD risk factors were taken into account ( = 0.440, 95% = 0.225–0.863, ). In additive inheritance mode, GG genotype carriers at rs4714210 exhibited a lower risk of CAD than AA carriers ( = 0.475, = 0.232–0.970, ). Conclusion. In type 2 diabetic patients from a Han Chinese population, some variations in the GLP-1R gene were associated with a lower risk of developing CAD.
Coronary artery disease (CAD) is a life-threatening condition that is a frequently occurring complication in patients with type 2 diabetes mellitus (T2DM), with diabetic patients being 2–4 times more likely to develop CAD than nondiabetics . Determination of genetic variants associated with CAD development in T2DM patients may assist in the identification of at-risk individuals and allow targeting of primary prevention and early intervention measures. In recent years, glucagon-like peptide 1 receptor (GLP-1R) agonists such as exenatide and liraglutide have been widely studied because of their glucose-dependent insulinotropic effects  and their other physiological effects such as decrease in fatty acid absorption, increase in satiety, and reduction in body weight . GLP-1Rs are widely expressed in the cardiovascular system , and a number of beneficial effects that protect against coronary heart disease are associated with the GLP-1/GLP-1R signal pathway and its agonist interactions. Previous studies demonstrated that GLP-1 agonists could reduce the rate of the first occurrence of death from cardiovascular causes and nonfatal myocardial infarction among patients with T2DM [5, 6]. GLP-1 agonists were also found to improve heart function, decrease the size of infarct areas in ischemia-reperfusion heart models (pig and mouse) , increase coronary blood flow in isolated mouse heart , and reduce monocyte adhesion and atherosclerotic lesions in apoE−/− mice . Consequently, we reasoned that genetic variation in the GLP-1R gene might affect CAD risk in patients with T2DM. This hypothesis was investigated by examining CAD-positive and CAD-negative patients with T2DM in a Chinese Han population.
2. Materials and Methods
2.1. Ethics Statement
The study protocol and informed consent procedures were approved by the Research Ethics Committees of Peking University First Hospital. Written informed consents were acquired from all subjects participating in this study, in agreement with the 1975 Helsinki Declaration.
Diabetes mellitus was diagnosed according to World Health Organization criteria (1999)  as follows: fasting plasma glucose ≥ 7.0 mmol/l, and/or 2 h plasma glucose ≥ 11.1 mmol/l, or casual plasma glucose (random blood sugar) ≥ 11.1 mmol/l. Patients with type 1 diabetes and subjects with active inflammatory conditions, autoimmune diseases, malignancies, usage of immunosuppressive drugs, and known hematological disorders were excluded. In total, 611 unrelated Chinese Han subjects with T2DM were included in the study: 394 with coronary artery stenosis (CAS) ≥ 50% (CAD-positive group) and 217 with CAS < 50% (control group). Diagnostic procedures were carried out at Peking University First Hospital. CAS ≥ 50% individuals were defined as those who exhibited ≥ 50% stenosis in at least one of the major coronary arteries or their main branches upon cardiac catheterization. Control individuals had <50% coronary stenosis in all main coronary arteries and main branches as determined by cardiac catheterization or high specificity spiral computed tomography (CT) scan . Demographic data and patient cardiovascular risk factor data were collected for all subjects from medical records. These data comprised gender, age, body mass index (BMI), fasting plasma glucose (FPG), history of dyslipidemia, hypertension (blood pressure ≥ 140/90 mmHg or receiving any antihypertensive therapies), and smoking history (“ever” or “never,” with “ever” defined as having smoked more than one cigarette per day for more than 6 months, as per World Health Organization criteria).
2.3. Single Nucleotide Polymorphism Genotyping
Genomic DNA was extracted from peripheral blood using a Whole Blood DNA Extraction Kit (BioTeke).
The GLP-1R gene is located at chromosome 6p21, is 38.9 kb in length, and includes 13 exons. In total, 33 haplotype-tagging single nucleotide polymorphisms (SNPs) were identified at the GLP-1R locus in the CHB (Han Chinese from Beijing) population from the HapMap Phase II database (http://www.hapmap.org) (R#27, , ). Eleven of these SNPs, dispersed across 22,058 bp of the total 38,964 bp of GLP-1R, were selected for further study: rs761387 (T>C), rs2268635 (G>A), rs7769547 (G>A), rs910162 (T>A), rs3765468 (G>A), rs3765467 (G>A), rs3765466 (A>T), rs10305456 (C>T), rs10305518 (T>G), rs1820 (T>A), and rs4714210 (A>G). Target regions were amplified by PCR. Direct DNA sequencing was used for 8 of the 11 selected SNPs, using a MassARRAY system (Sequenom iPLEX assay, San Diego, CA, USA) , and the remaining three SNPs were genotyped using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) analysis. Genotyping success rates were 95–100%, and repeatability rates were 98–100%. To validate PCR-RFLP assays, 5% of amplicons were directly sequenced to confirm the genotypes for each SNP. Concordance rates between RFLP and DNA-sequencing results were 98–100%.
2.4. Data Analysis
Clinical and laboratory data were expressed as means ± SD or percentages. Genotype distributions described departure from Hardy-Weinberg equilibrium at each polymorphic locus. Linkage disequilibrium (LD) and haplotype analysis were performed using Haploview 4.2, with haplotypes estimated using an accelerated expectation-maximization algorithm.
Allele frequencies were determined by gene counting. SNP association with risk of CAD was assessed using the SPSS statistical package (SPSS version 19.0, USA). Qualitative variables were compared using a test, and quantitative variables were compared using an independent samples -test or a Mann–Whitney test. Associations between CAD and genotype were analyzed using multiple logistic regression with adjustment for the following potential confounders: age, gender, BMI, smoking status, positive histories of dyslipidemia and hypertension, and diabetic duration. As a descriptive measure of association between genotypes and outcomes, was considered to be statistically significant and odds ratios (ORs) were calculated with 95% confidence intervals (CIs). Bonferroni correction was used to correct for multiple comparisons. Power and Sample Size Calculation software (version 3.1.2, 2014) was used for power calculations .
3.1. Characteristics of Study Subjects
Males were more likely to be in the CAD-positive T2DM group than in the control T2DM group (); otherwise, no significant differences in phenotypic characteristics were found between groups (Table 1). HbA1c and FPG measurements were acquired after antidiabetic treatment. Genotype distributions at all 11 loci were in agreement with Hardy-Weinberg equilibrium (data not shown). Statistical power was 0.99.
3.2. Allele and Genotype Analysis
The minor allele G at rs4714210, in the 3 untranslated region (UTR), was found more frequently in the control group than in the CAD-positive group, and carriers of the G allele displayed a lower risk of CAD when compared with noncarriers ( = 0.783, 95% CI = 0.613–1.002, ). In dominant inheritance mode, no significant difference in genotype distribution was found between the CAD-positive and control groups (Supplementary Material 1).
The protective effect of the homozygous GG minor allele genotype at rs4714210 was also observed in additive (codominant) inheritance mode using logistic regression analysis (Table 2). GG genotype carriers at rs4714210 exhibited a lower risk of CAD than AA carriers ( = 0.475, = 0.232–0.970, = 0.041, after adjustment for confounders as above), while no such protective effect was observed in heterozygote carriers.
In recessive inheritance mode, the carriers of genotype GG at rs4714210 had a decreased risk of CAD ( = 0.442, 95% CI = 0.258–0.757, ; = 0.440, 95% CIa = 0.225–0.863, ), after adjusting for other known CAD risk factors (gender, age, BMI, smoking status, dyslipidemia history, hypertension history, and diabetic duration) (Table 3).
The other ten SNPs tested in this study displayed similar allele frequencies between the CAD-positive and control groups, and no significant associations were noted between genotype and CAD risk (Tables 2 and 3).
3.3. Haplotype Analysis
Haploview plotting was used to construct haplotypes depending on the physical position and the value of () between each pair of SNPs in one block. Three blocks were delineated as follows: LD block 1 (rs910162, rs3765468, rs3765467, rs765466, and rs10305456), block 2 (rs761387 and rs7769547), and block 3 (rs10305518, rs1820, and rs4714210). The SNP distributions in the three haplotype blocks did not differ significantly between the CAD-positive and control groups (Supplementary Material 2).
GLP-1R is a 463 amino acid member of the class B GPCR secretin family. GLP-1R is a classic seven-transmembrane protein; the C-terminus of GLP-1R interacts with a signaling G protein, and the large N-terminal extracellular domain plays an important role in ligand binding [13, 14]. After binding to the GLP-1 ligand, GLP-1R transmits a signal through a Gas-coupled subunit. This induces an increase in cAMP (cyclic adenosine monophosphate) levels and consequently activates the PKA pathway. The effect of the GLP-1/GLP-1R pathway in the myocardial ischemia-reperfusion model was previously summarized by Ravassa et al. . In particular, GLP-1/GLP-1R can activate the PKA, PI3K, MEK1/2, and eNOS pathways, resulting in cardiovascular protective effects such as reduced apoptosis, improved energy metabolism, reduced inflammation in myocardial cells, and vasodilation in myocardial arteries. Systematically, GLP-1R transmits signals that prompt insulin secretion increases, appetite reduction, metabolism improvement, and lower blood pressure and, as a result, decreases the severity of atherosclerotic lesions . The wide-ranging effects of GLP-1/GLP-1R suggest that variations in the GLP-1R gene may contribute to the risk of CAD.
Here, we found an association between rs4714210 in the GLP-1R gene and CAD risk in T2DM patients in the Chinese Han population. Patients homozygous for the minor allele G at rs4714210 exhibited a 50% lower risk of CAD than other genotype carriers. The mechanisms through which this allele confers protection are unclear. The rs4714210 locus is in the 3 UTR of the GLP-1R gene, and 3 UTRs are thought to play important roles in gene regulation. For example, 3 UTRs can influence chromosome structure, regulate transcription, stabilize mRNA, and modulate translation, thus affecting the stability and transport of the encoded proteins . We therefore speculate that variations at rs4714210 may differentially affect the function of GLP-1R through one or more of these mechanisms, but this remains to be confirmed. It is also likely that the rs4714210 SNP is in strong LD with other SNPs that have biological effects.
GLP-1R SNPs have been confirmed in the association with obesity , pancreatic beta-cell function [18, 19], and T2DM  in different populations. However, there are few studies about the variations of GLP-1R with CAD. In 2016, Scott et al. first observed an effect of GLP-1R genetic variation in Caucasian CAD patients with or without T2DM , identifying an association between GLP-1R rs10305492 and CAD in Caucasians (). Although SNP rs10305492 is not found in the Chinese Han population, SNPs rs10305492 and rs4714210 are in complete LD in Caucasians (), indicating that our results are concordant with Scott’s.
We acknowledge some limitations of this study. Sample size was relatively small, for in the cases were only 26 type 2 diabetes patients with one CAD vessel affected, 83 and 285 for two and three CAD vessels affected, respectively, so we did not stratify the cases and analyze the association with the number of affected vessels. And clinical features were not perfectly matched, and urine albumin creatinine ratio (ACR) was not collected, between the case and control groups. Both of them may introduce bias. Moreover, further functional studies on genetic variations at the GLP-1R locus would be beneficial. If our findings were confirmed through prospective studies, GLP-1R polymorphisms could be used as predictors of CAD risk in patients with T2DM in the Chinese Han population.
In T2DM patients from a Han Chinese population, some variations in the GLP-1R gene were associated with a lower risk of developing CAD.
The study protocol and informed consent procedures were approved by the Research Ethics Committees of Peking University First Hospital (no. 2007-026).
Written informed consents were acquired from all subjects participating in this study, in agreement with the 1975 Helsinki Declaration.
Conflicts of Interest
We declare that we have no conflict of interest.
Ran Lu analyzed data and wrote the manuscript. Xiaowei Ma edited the manuscript. Xiaowei Wei, Ge Bai, Jianwei Zhang, Ruifen Deng, Nan Gu, and Nan Feng analyzed data and contributed to the discussion. Xiaohui Guo reviewed the manuscript and contributed to the discussion. Xiaowei Ma accepts responsibility for the article.
The authors are very grateful to the staff of the Departments of Cardiology, Radiology, Statistics, and others, from Peking University First Hospital, for excellent contributions to data collection. We also wish to thank all the patients who participated in this study. The study was supported by National Natural Science Foundation of China (NSFC 30771033) and National Science and Technology Major Project (2011ZX09307-001-08).
Supplementary Material 1: distribution of SNP genotype frequencies in CAD-positive and control groups in dominant inheritance mode. Supplementary Material 2: association of common haplotypes with CAD risk. (Supplementary Materials)
K. Ban, M. H. Noyan-Ashraf, J. Hoefer, S. S. Bolz, D. J. Drucker, and M. Husain, “Cardioprotective and vasodilatory actions of glucagon-like peptide 1 receptor are mediated through both glucagon-like peptide 1 receptor-dependent and -independent pathways,” Circulation, vol. 117, no. 18, pp. 2340–2350, 2008.View at: Publisher Site | Google Scholar
K. G. Alberti and P. Z. Zimmet, “Definition, diagnosis and classification of diabetes mellitus and its complications. Part 1: diagnosis and classification of diabetes mellitus. Provisional report of a WHO consultation,” Diabetic Medicine, vol. 15, no. 7, pp. 539–553, 1998.View at: Publisher Site | Google Scholar
J. K. Min, A. Dunning, F. Y. Lin et al., “Age- and sex-related differences in all-cause mortality risk based on coronary computed tomography angiography findings results from the International Multicenter CONFIRM (Coronary CT Angiography Evaluation for Clinical Outcomes: An International Multicenter Registry) of 23, 854 patients without known coronary artery disease,” Journal of the American College of Cardiology, vol. 58, no. 8, pp. 849–860, 2011.View at: Publisher Site | Google Scholar
M. P. Millis, “Medium-throughput SNP genotyping using mass spectrometry: multiplex SNP genotyping using the iPLEX® Gold assay,” in Disease Gene Identification, J. DiStefano, Ed., vol. 700 of Methods in Molecular Biology (Methods and Protocols), pp. 61–76, Humana Press, Totowa, NJ, USA, 2011.View at: Publisher Site | Google Scholar
C. H. Tsai and T. W. Dreher, “Second-site suppressor mutations assist in studying the function of the 3 noncoding region of turnip yellow mosaic virus RNA,” Journal of Virology, vol. 66, no. 9, pp. 5190–9, 1992.View at: Google Scholar