Abstract
Background:
This study examined the association in a South Indian population with gestational diabetes mellitus (GDM) of type 2 diabetes risk variants that have previously conferred susceptibility to GDM in other populations.
Subjects and Methods:
The study groups comprised 518 women with GDM and 910 pregnant women with normal glucose tolerance (NGT). Women with GDM were recruited from a tertiary diabetes center in Chennai, in south India, and NGT women were selected from antenatal clinics also in Chennai. Genomic DNA was isolated from whole blood using the phenol chloroform method. Twelve previously reported GDM-associated single nucleotide polymorphisms (SNPs) in or near nine loci were genotyped using the MassARRAY™ system (Sequenom, San Diego, CA).
Results:
Among the 12 SNPs genotyped, 11 SNPs were in Hardy–Weinberg equilibrium and had a call rate of >95%. Of the 11 SNPs previously associated with GDM in other populations, significant association was observed only with the rs7754840 and rs7756992 SNPs of the CDK5 regulatory subunit associated protein 1-like 1 (CDKAL1) gene in this population. The minor alleles of the SNPs rs7754840 and rs7756992 showed significant susceptibility to GDM with an odds ratio of 1.34 (95% confidence interval, 1.12–1.60; P=0.0013) and 1.45 (95% confidence interval, 1.21–1.72; P=0.00004), respectively.
Conclusions:
The rs7754840 and rs7756992 SNPs of the CDKAL1 gene were found to be associated with GDM in this south Indian population. This is the first study describing genetic susceptibility of GDM in Asian Indians.
Introduction
G
Evidence is also accumulating that susceptibility to GDM has a genetic component. 9 Studies have shown clustering of type 2 diabetes mellitus (T2DM) and impaired glucose tolerance in families with GDM and evidence for higher prevalence of T2DM in mothers of women with GDM. 10 Earlier studies have indicated that the main pathophysiological defects that occur in GDM may be the same as those observed with T2DM, such as marked insulin resistance, impairment of insulin secretion, and abnormal utilization of glucose. 9,11 Several studies, including genome-wide association studies (GWASs) and meta-analysis reports, have provided evidence that GDM and T2DM may share similar risk factors and genetic background. 12 –17
Recently, replication of T2DM-associated genes identified in various GWASs has been carried out in GDM patients in different Asian populations from Korea, China, and Malaysia, 4,18 –20 but there are no data for Asian Indians, to the best of our knowledge. Although the genetics of T2DM have been extensively studied in Asian Indians using different approaches such as candidate gene and GWAS, 21 –29 GDM, which is probably one of the early indicators of future occurrence of T2DM in women, has not been well studied in India. We have previously reported that associations of certain genetic variants were different in Asian Indians compared with other ethnic groups. 21 –23,27,30 The replication of GWASs in our south Indian population showed that out of 45 single nucleotide polymorphisms (SNPs) that showed an association in white T2DM populations, only six—namely, CDK5 regulatory subunit associated protein 1-like 1 (CDKAL1) (rs7754840, rs7756992, and rs6931514), HHEX (rs7923837), and BAZ1B (rs12056034)—showed association with T2DM in our population. 26 Moreover, a GWAS on south Asians with T2DM has shown the association of six novel genes (AP3S2 rs2028299, GRB14 rs3923113, HMG20A rs7178572, HNF4A rs4812829, ST6GAL1 rs16861329, and VPS26A rs1802295) that were not found in Europeans. 28 Therefore, it is likely that GDM might also show differences in genetic association in our population. Hence, it is relevant to study the role of genetic factors in GDM in Asian Indians.
The present study was designed to invesigate whether the SNPs that have previously shown association with GDM in other populations are also associated with GDM in this south Indian population. To our knowledge, this is the first major study of genetic susceptibility of GDM in Asian Indians.
Subjects and Methods
Study subjects
Pregnant women without previous diagnosis of glucose intolerance were selected for the present study. The study groups comprised 518 women with GDM and 910 pregnant women with normal glucose tolerance (NGT). Women with GDM were recruited from a tertiary diabetes center in Chennai. Diagnosis of GDM is done using International Association of Diabetes and Pregnancy Study Groups criteria: one or more of fasting plasma glucose level of ≥5.1 mmol/L, 1-h plasma glucose level of ≥10.0 mmol/L, and 2-h plasma glucose level of ≥8.5 mmol/L following a 75-g oral glucose tolerance test. Women with NGT (fasting plasma glucose level of <5.6 mmol/L [100 mg/dL] and 2-h postprandial blood glucose level of <7.8 mmol/L [140 mg/dL]) were selected from antenatal clinics in Chennai. Women with advanced maternal age (>40 years), prepregnant body mass index of ≥30 kg/m2, previous history of still birth, a macrosomic baby, or a baby with congenital abnormality, pregnancy-induced or preexisting hypertension, and history of other complicated pregnancies were excluded from the study. The subjects included in the study were all unrelated individuals of south Indian origin and Dravidian ethnicity as per the self-reported information. The study has been approved by the Madras Diabetes Research Foundation Ethics Committee, and written informed consent was obtained from all study subjects included in the study.
SNP selection and genotyping
Genomic DNA was isolated from whole blood using the phenol chloroform method. The SNPs that were selected for the present study include 12 previously reported GDM-associated SNPs from nine loci: CDKAL1 (rs7754840 and rs7756992), glucokinase (hexokinase 4) (GCK) (rs1799884 and rs4607517), insulin-like growth factor 2 mRNA binding protein 2 (IGF2BP2) (rs4402960), potassium inwardly-rectifying channel, subfamily J, member 11 (KCNJ11) (rs5219), potassium voltage-gated channel, KQT-like subfamily Q, member 1 (KCNQ1) (rs2237892), melatonin receptor 1A (MTNR1A) (rs2119882), melatonin receptor 1B (MTNR1B) (rs1387153 and rs10830963), serine racemase (SRR) (rs391300), and transcription factor 7-like 2 (TCF7L2) (rs12255372). All the SNPs were genotyped using the MassARRAY™ system (Sequenom, San Diego, CA) following the manufacturer's instructions. SpectroTYPER™ software (Sequenom) automatically called genotypes, and only conservative and moderate calls were accepted for this study. Ten percent of the samples were replicated, and the discordance rate was less than 0.1%. The call rate required for each assay was ≥90%.
Statistical analysis
Hardy–Weinberg equilibrium was tested using Pearson χ2 statistics with a threshold of P<0.05 in cases and controls separately for each SNP. Statistical analyses were performed in SPSS version 15.0 software (SPSS, Inc., Chicago, IL) and PLINK, an open-source tool (
Results
The clinical and biochemical parameters of the groups of GDM and pregnant NGT women are presented in Table 1. The mean age, fasting plasma glucose, 2-h postload plasma glucose, and glycosylated hemoglobin were significantly higher in the GDM group than in controls (P<0.001).
Data are mean±SD values.
Indicates a statistically significant value (P value<0.05).
BMI, body mass index; GDM, gestational diabetes mellitus; NGT, normal glucose tolerance.
In total, 12 SNPs (in or near nine loci) were genotyped in 518 women with GDM and 910 pregnant women with NGT, of which 11 SNPs were in Hardy–Weinberg equilibrium and had a call rate of >95%: rs7754840, 97.2%; rs7756992, 96.7%; rs1799884, 97.4%; rs4607517, 99.1%; rs4402960, 97.1%; rs5219, 97.8%; rs2119882, 97.4%; rs1387153, 97.4%; rs10830963, 97.9%; rs391300, 95.4%; and rs12255372, 97.3%. The KCNQ1 (rs2237892) SNP with a call rate of 98.2% was excluded from further analysis as it deviated from Hardy–Weinberg equilibrium. The P value of Hardy–Weinberg equilibrium for the failed SNP rs2237892 remained significant even after multiple corrections using Bonferroni's correction (P>0.05/12=0.004). The distributions of the allelic frequencies of all the GDM-associated SNPs were compared between the GDM and pregnant NGT women (Table 2). Of the 11 SNPs in or near eight loci (CDKAL1, GCK, IGF2BP2, KCNJ11, MTNR1A, MTNR1B, SRR, and TCF7L2) that were tested, significant association with GDM was seen only with two SNPs of the CDKAL1 gene: rs7754840 (P= 0.0013) and rs7756992 (P=0.00004). The association of rs7754840 and rs7756992 with GDM remained significant even after multiple corrections using Bonferroni's correction (P>0.05/11=0.0045). Association with GDM was not observed in the other nine SNPs. However, the SNP of the TCF7L2 gene (rs12255372) showed a moderately significant association (P=0.054) with GDM.
Indicates a statistically significant value (P value<0.05). They are less than the significance threshold after Bonferroni's correction (0.05/11=0.0045).
GDM, gestational diabetes mellitus; MAF, minor allele frequency; NGT, normal glucose tolerance; SNP, single nucleotide polymorphism.
Logistic regression analysis was performed for the CDKAL1 gene SNPs rs7754840 and rs7756992 under multiple genetic models, and the findings are presented in Table 3. The minor alleles of the SNPs rs7754840 and rs7756992 showed significant susceptibility to GDM with ORs of 1.34 (95% CI, 1.12–1.60; P=0.0013) and 1.45 (95% CI, 1.21–1.72; P=0.00004) respectively. Under the dominant model, the SNP rs7756992 showed significant association with GDM (adjusted OR=1.37; 95% CI, 1.07–1.74; P=0.012), whereas SNP rs7754840 did not show any association (adjusted OR=1.20; 95% CI, 0.94–1.50; P=0.138). Adjusting this P value by the number of competing models (three) for two SNPs using Bonferroni's correction, we found that the association finding did not remain significant (Bonferroni threshold P=0.05/6=0.008); therefore the association of both of the SNPs did not follow the dominant mode of inheritance. Although the heterozygous variant genotype showed significant association with GDM in both the SNPs, the homozygous variant failed to show any association. This showed that the genotypic distribution of both the SNPs was not compatible under an additive model. However, the SNPs rs7754840 (adjusted OR=1.93; 95% CI, 1.20–3.00; P=0.005) and rs7756992 (adjusted OR=2.12; 95% CI, 1.35–3.30; P=0.001) of the CDKAL1 gene showed association with GDM under the recessive model. Therefore the association of the CDKAL1 gene SNPs with GDM fits well into a recessive model of inheritance.
Data are number of subjects (%).
Data are adjusted (for age) odds ratio (95% confidence interval); P value.
Indicates a statistically significant difference (P value<0.05). They are less than the significance threshold after Bonferroni's correction.
Indicates a statistically significant difference (P value<0.05).
GDM, gestational diabetes mellitus; NGT, normal glucose tolerance; SNP, single nucleotide polymorphism.
Analysis of linkage disequilibrium found the SNPs of CDKAL1 (rs7754840 and rs7756992; D′=0.91), GCK (rs1799884 and rs4607517; D′=1.0), and MTNR1B (rs1387153 and rs10830963; D′=0.80) were in a strong linkage disequilibrium. Table 4 shows the comparision of diplotype frequencies of CDKAL1 (rs7754840 and rs7756992), GCK (rs1799884 and rs4607517), and MTNR1B (rs10830963 and rs1387153) gene variants in the study subjects. The CC/GG diplotype CDKAL1 gene with a frequency of 4.1% in pregnant women with NGT and 8.8% in GDM was significantly different, and this diplotype showed association with a 2.1 times increased risk of GDM (adjusted P=0.019; OR=2.08; 95% CI, 1.13–3.84).
Indicates a statistically significant difference (P value<0.05).
Odds ratio (OR) adjusted for age.
CI, confidence interval; GDM, gestational diabetes mellitus; NGT, normal glucose tolerance; SNP, single nucleotide polymorphism.
Discussion
Our study is the first major genetic study of GDM from India. Recent studies have shown that the T2DM-associated gene variants were also found to be associated with GDM susceptibility. 4,16 –20 Among the genes that were investigated in the present study, CDKAL1 is the only gene that showed association with GDM in our south Indian population. Both SNPs of the CDKAL1 gene—rs7754840 and rs7756992—showed association with GDM even after Bonferroni's correction with a power of 61% and 82%, respectively.
The impact of CDKAL1 gene variants in our population can be noticed in that the minor alleles of rs7754840 and rs7756992 are shown to be the risk alleles and are significantly associated with GDM with an increased risk of 1.34 and 1.45 times, respectively. Logistic regression analysis performed for the CDKAL1 gene SNPs rs7754840 and rs7756992 showed that the association of the CDKAL1 gene SNPs with GDM fits well into a recessive model of inheritance.
The study of these CDKAL1 gene SNPs with GDM in the Malaysian population showed its strongest evidence of association under the additive model. 4 Analysis of different genetic models showed that the association of these SNPs with T2DM in this south Indian population best fits the additive and dominant models. 26 The minor alleles of both the associated CDKAL1 gene variants were shown to be the risk alleles. Previous studies among Asian populations have reported the association between SNPs rs7754840 and rs7756992 and GDM risk. 15,16,19 In the Korean population, the minor alleles of the SNPs rs7754840 and rs7756992 showed significant susceptibility to GDM with ORs of 1.55 (95% CI, 1.34–1.79; P=4.17×10–9) and 1.39 (95% CI, 1.20–1.61; P=9.14×10–9), respectively. 19 However, no significant association was found in the Chinese population. 20 The meta-analysis among Asian populations indicated that the C allele of rs7754840 was significantly associated with risk of GDM with a pooled OR of 1.40 (95% CI, 1.13–1.72; P=0.002). 17 In the present study, haplotype analysis of the CDKAL1 gene SNPs showed that both SNPs of CDKAL1 (rs7754840 and rs7756992) were in a strong linkage disequilibrium, and the risk allele haplotype CG was found to be significantly associated with GDM (P=0.0216). The risk of rs7754840 and rs7756992 SNPs is strengthened by the diplotypes containing the risk alleles of both these SNPs. The CC/GG diplotype of the CDKAL1 gene showed a stronger association with GDM (OR=2.1) compared with the single risk allele.
Although the function of the CDKAL1 gene is still unclear, earlier studies have shown that this gene may play a role in β-cell function by inhibiting activity of cyclin-dependent kinase 5, which prevents a decrease in insulin secretion due to glucotoxicity. 31,32 The CDKAL1 protein has been reported to function as a tRNA modification enzyme, and CDKAL1 risk variants affect the protein translation, thereby impairing the conversion of proinsulin to insulin in pancreatic β-cells. 33,34
The CDKAL1 gene has been previously reported to be associated with T2DM in this south Indian population. 26 The association of the SNP rs7756992 with GDM risk has been reported in the Malaysian population. 4 The rs7754840 SNP has shown association with GDM in earlier studies, mostly among Asian populations. 4,15,19,20 The directionality of the association of the CDKAL1 gene SNPs between cases and controls was found to be similar to that observed in other studies in GDM women. 4,15,19,20 The frequency of the C allele of the rs7754840 variant of the CDKAL1 gene in our study is lower when compared with that among other Asian ethnic groups like Korean (approximately 56%), Chinese (approximately 47%), and Malaysian (approximately 38%). 4,15,19,20 Similarly, the frequency of the G allele of the rs7756992 variant of the CDKAL1 gene is also lower compared with the Malaysian study (approximately 44%). 4
In the present study, we observed that most of the common T2DM-associated SNPs that showed susceptibility to GDM in other populations did not show any association with GDM in this south Indian population. Earlier studies have reported the significant association of the SNPs GCK (rs1799884 and rs4607517), IGF2BP2 (rs4402960), KCNJ11 (rs5219), MTNR1A (rs2119882), MTNR1B (rs1387153 and rs10830963), SRR (rs391300), and TCF7L2 (rs12255372) with GDM in whites and Asians, mainly Koreans, Chinese, and Malayasians. 4,16 –20 The TCF7L2 gene SNP rs12255372 showed only a moderate significant association (P=0.054) with GDM in the present study, although this SNP has previously shown association with T2DM in Asian Indians. 30 The association of rs12255372 with GDM has been reported in other populations like Mexican Americans, whites, and Koreans. 14,19,35,36 The frequency of the T allele was found to be lower (approximately 24%) in our GDM group compared with other populations such as Mexican Americans (approximately 39%), 14 whites (approximately 31%), 35,36 and Asians (Koreans; approximately 40%). 19 Analysis of the present study was limited to candidate SNPs that have been validated in multiple large studies mainly with the concern of multiple testing. This is the limitation of our study.
The most significant finding of the present study is the identification of SNPs of the CDKAL1 gene, which is found to be strongly associated with GDM in this south Indian population. Testing these SNPs with a larger sample size and in a prospective study would help to validate our findings and also help in evaluating these SNPs for their ability to predict development of future diabetes mellitus in this population.
Footnotes
Acknowledgments
This study was supported by funding from the MDRF Innovative Research Fund of the Madras Diabetes Research Foundation, Chennai, India.
Author Disclosure Statement
No competing financial interests exist.
