Skip to main content

Genetic variants of the GLP-1R gene affect the susceptibility and glucose metabolism of gestational diabetes mellitus: a two-center nested case‒control study

Abstract

Background

Gestational diabetes mellitus (GDM) is the most common complication during pregnancy, occurring under the combined action of environmental and genetic factors. Genetic variants of glucagon-like peptide-1 receptor (GLP-1R) have been reported to affect insulin secretion and susceptibility to type 2 diabetes. This study aimed to explore the role of GLP-1R polymorphisms in GDM and glucose metabolism.

Methods

A two-center nested case‒control study was designed, including 200 pregnant women with GDM and 200 pregnant women without GDM genotyped for five tag SNPs of GLP-1R using Sanger sequencing. Logistic regression was used to evaluate the relationship between GLP-1R polymorphisms and GDM risk. Glucose and insulin concentrations were measured based upon the 75 g oral glucose tolerance test (OGTT). Beta cell function of different genotypes was estimated with the 60 min insulinogenic index (IGI60) and OGTT-derived disposition index (DI).

Results

Mutant genotype AG + GG of tag SNP rs6458093 nominally increased GDM risk (p = 0.049), especially among subjects younger than 35 years (p = 0.024) and with BMI no less than 24 (p = 0.041), after adjusting for confounders. Meanwhile, compared with subjects with wild genotype AA, subjects with genotype AG + GG of rs6458093 also showed nominally significantly lower IGI60 (p = 0.032) and DI (p = 0.029), as well as significantly higher 75 g OGTT-based 1 h glucose load plasma glucose levels (p = 0.045). Moreover, the mutant heterozygous genotype GA of tag SNP rs3765467 nominally decreased GDM risk among subjects older than 35 years (p = 0.037) but showed no association with insulin secretion and glucose homeostasis.

Conclusions

Tag SNP rs6458093 of GLP-1R was nominally associated with increased GDM risk and affected beta cell function and postprandial glucose metabolism, while tag SNP rs3765467 of GLP-1R was nominally associated with decreased GDM risk, providing evidence for molecular markers and etiological study of GDM.

Introduction

Gestational diabetes mellitus (GDM) is one of the most common complications during pregnancy, affecting 14% of pregnant women worldwide [1]. The prevalence of GDM among pregnant Chinese women increased remarkably from 4% in 2010 to 21% in 2020 [2]. Macrosomia and neonatal hypoglycemia are the most common complications of GDM [3]. In addition, women with a history of GDM are more likely to develop type 2 diabetes mellitus (T2DM) [4] and cardiovascular disease [5], while their offspring are more likely to become obese [6] and develop T2DM [7] in the future.

Insulin resistance and pancreatic beta cell (β cell) dysfunction are considered to be the main pathogenesis of GDM. To meet the energy needs of the growing fetus by limiting the mother's consumption of glucose, insulin sensitivity decreases due to hormones, obesity, the placenta, cytokines, adipokines, decreased physical activity and so on [8]. Pregnant women with normal glucose tolerance (NGT) can compensate for increasing insulin resistance by increasing the secretion of insulin, while women with GDM have insufficient insulin secretion [9]. In other words, excellent β cell function secreting insulin plays a crucial role in glucose homeostasis during pregnancy.

It is well known that food can stimulate intestinal mucosa L cells to release incretin hormones, which stimulate insulin secretion in a glucose-dependent manner, a phenomenon known as the “incretin” effect [10, 11]. Insulin secretion stimulated by a meal is partly due to direct actions of blood glucose on β cells, while up to two-thirds of postprandial insulin is due to the incretin effect [11, 12]. Glucagon-like peptide-1 (GLP-1) is one of the most important incretin hormones [13], acting through glucagon-like peptide-1 receptor (GLP-1R) signaling. It is reasonable to assume that dysfunction of GLP-1R may affect the function of β cells secreting insulin by inhibiting the interaction between GLP-1 and GLP-1R. Although the development of GDM is closely associated with dysfunction of β cells, the relationship between GDM susceptibility and GLP-1R genetic variation that directly influences the quality or quantity of GLP-1R and β cell function remains unknown.

A previous study found that two SNPs (rs6923761 and rs3765467) in GLP-1R could alter insulin secretion in response to exogenous GLP-1 [14]. In recent years, rs3765467 and rs10305492 of the GLP-1R gene were reported to affect β cell function and promote apoptosis of β cells [15]. Meanwhile, rs3765467 was identified as a susceptibility locus for T2DM in the Japanese population in a meta-analysis of 4 genome-wide association studies (GWAS) [16], and rs10305492 was reported to be linked to lower fasting plasma glucose (FPG) and T2DM risk [17, 18]. Moreover, rs1042044 and rs367543060 of the GLP-1R gene were reported to be associated with susceptibility to T2DM [19, 20]. It is well known that GDM and T2DM share the same pathogenesis (insulin resistance and pancreatic β cell dysfunction) [21] and share a common genetic background [22]. The percentage of GDM women diagnosed with T2DM was 12% higher for each additional year after pregnancy [23]. The estimated risks for T2DM were 19.72% at 10 years, 29.36% at 20 years, 39.00% at 30 years, 48.64% at 40 years, and 58.27% at 50 years [24]. Therefore, GDM is believed to be the early stage of T2DM [25]. The above studies provided a rationale for the hypothesis that GLP-1R polymorphisms might be involved in the pathogenesis of GDM. Herein, we validated the correlation of GLP-1R gene polymorphisms with GDM risk and evaluated the effects of GLP-1R polymorphisms on β cell function and glucose homeostasis to provide preliminary evidence for the genetic and etiological study of GDM.

Materials and methods

Study design and participants

It is a two-center nested case‒control study. All subjects were from the Jingzhou maternal and child health cohort. Pregnant women who registered in Jingzhou Hospital Affiliated to Yangtze University and Gongan County Maternal and Child Health Care Hospital in Hubei Province of China from April 2021 to December 2021, and were willing to undergo routine pregnancy check-ups in these two hospitals and provide informed consent for participation could be recruited into the cohort. All pregnant women were recruited into the cohort on the day of the first pregnancy check-up. The diagnosis of GDM was based upon the Health Industry Standard of the People's Republic of China of “Diagnosis criteria for gestational diabetes mellitus” (WS 331-2011) which adopts recommendations of the diagnostic criteria of the International Association of Diabetes and Pregnancy Study Groups (IADPSG). All pregnant women in the cohort were routinely tested for FPG, glycosylated hemoglobin a1c (GHbA1c), or random plasma glucose at the first prenatal visit to identify prepregnancy diabetes (overt diabetes). Pregnant women were diagnosed as having overt diabetes and excluded from our study when (1) FPG ≥ 7.0 mmol/L (126 mg/dL); (2) GHbA1c ≥ 6.5% (DCCT/UKPDS standardized); or (3) having typical symptoms of hyperglycemia or hyperglycemic crisis and random plasma glucose ≥ 11.1 mmol/L (200 mg/dL) (the values of FPG and GHbA1c need to be rechecked on another day when there are no obvious symptoms of hyperglycemia). All the pregnant women without prepregnancy diabetes underwent a 75 g OGTT at 24–28 gestational weeks. Subjects were diagnosed as GDM when at least one value from the 75 g oral glucose tolerance test (OGTT) was no less than the following levels: FPG 5.1 mmol/L, 1-h plasma glucose (1hPG) 10.0 mmol/L and 2-h plasma glucose (2hPG) 8.5 mmol/L. Pregnant women in the cohort who underwent the 75 g OGTT at 24–28 gestational weeks and met the above diagnostic criteria were included in the case group (GDM group). Age-, basal BMI- and hospital visit date of OGTT-matched normal glucose tolerance (NGT) individuals in the cohort were selected as controls (NGT group) at a frequency of 1:1. When multiple matched controls were available, the subject who matched best was selected. Subjects who met the following criteria were excluded: (1) had other types of diabetes; (2) had heart, liver or kidney disease, hypertension, thyroid dysfunction and polycystic ovary syndrome (PCOS) identified before the diagnosis of GDM (before 24–28 gestational weeks); (3) had twin or multiple pregnancies; (4) conceived by assisted reproduction technology; (5) had taken drugs affecting glucose metabolism during pregnancy; and (6) had incomplete pregnancy check-up information. Finally, a total of 400 pregnant women (200 cases and 200 controls) were included in our study. The flow chart for the selection of subjects is shown in Fig. 1.

Fig. 1
figure 1

Flow chart for the selection of subjects. *The controls were matched with the cases according to age, basal BMI, hospital and hospital visit date of OGTT

Two milliliters of EDTA anticoagulant blood samples and serum samples based upon a 75 g OGTT were collected and then stored at − 80 °C until DNA extraction and insulin measurement. Single nucleotide polymorphism (SNP) genotyping was performed on all 400 subjects, among which insulin levels were measured only in 218 subjects (109 patients and 109 controls).

Selection of tag SNPs

We did not plan to select tag SNPs with data from public databases, such as the 1000 genome, because we would perform SNP genotyping and identify SNPs that influence GDM susceptibility using direct sequencing. Selecting tag SNPs based on the sequencing and genotyping data from our 400 study subjects is more representative of our study population. Two PCR fragments covering two exons and four introns of the GLP-1R gene, which contain high-density SNPs with a minor allele frequency (MAF) beyond 0.10 in the East Asian population and contain at least one of the few missense variant SNPs reported to be related to the susceptibility or pathogenesis of diabetes were analyzed by Sanger sequencing. Finally, 12 SNPs with MAFs beyond 0.10 in our study population were identified, including 2 missense variant SNPs, 1 synonymous change and 9 intronic variant SNPs. There were 4 SNPs (rs1042044, rs10305478, rs10305479, rs6458093) in fragment 1 (located in 39073424–39074250 Homo sapiens chromosome 6, GRCh38), including the missense variant SNP rs1042044 reported to be associated with susceptibility to T2DM [19], and 8 SNPs (rs910168, rs910167, rs910165, rs910163, rs910162, rs3765468, rs3765467, rs3765466) in fragment 2 (located in 39065201–39066245 Homo sapiens chromosome 6, GRCh38), including the missense variant SNP rs3765467 reported to be associated with insulin secretion, β cell apoptosis and T2DM [15]. The genotypes of the 400 subjects and chromosomal position information of the 12 SNPs were input into Haploview 4.2 software for linkage disequilibrium (LD) analysis and selection of tag SNPs [26]. Finally, Haploview identified one LD block and five tag SNPs (rs10305478 C > T, rs6458093 A > G, rs1042044 A > C, rs3765468 G > A, rs3765467 G > A) that could capture the 12 SNPs. Information for 12 SNPs captured by the 5 tag SNPs is shown in Table 1. Given the strong linkage of the 12 SNPs, 5 tag SNPs rather than 12 SNPs were used in the following analyses.

Table 1 Information on the tag SNPs

Laboratory measurements

Genomic DNA was isolated following the guidelines of the TIANamp Blood DNA Kit DP348 (Tiangen, China). PCR amplicons were designed to be approximately 1000 bp for optimal Sanger sequencing analysis and included as many SNPs as possible. Two primer pairs were designed for amplification of two selected fragments synthesized by Sangon Biotech (Wuhan) Co., Ltd. The primer sequences of fragment 1 (827 bp) were as follows: forward—5′-CAGGATAGTGAGGATTCAATG-3′, reverse—5′-CACTGTACTTGGCTCTCC-3′. The primer sequences of fragment 2 (1045 bp) were as follows: forward—5′-ACTCAGTGCCAACCTTGTTC-3′, reverse—5′-CCACCGTGTAGATGATGTAGAG-3′. Standard PCRs for the two fragments were performed in a final volume of 25 μL, containing 2 μL (50 ng) of genomic DNA, 2.5 μL 10 × Buffer (Mg2+ plus) (Takara), 2 μL dNTPs (Takara), 0.125 μL Taq DNA polymerase (5 U/μL) (Takara), and 1.25 μL forward and reverse primers (20 μM) (Sangon Biotech). The amplification conditions of fragments 1/2 were initial denaturation at 95 °C for 3 min, followed by 35 cycles of 98 °C for 10 s, 64. 4 °C (fragment 1)/63. 3 °C (fragment 2) for 30 s and 72 °C for 1 min. The specificity of the PCR products was verified by electrophoresis in a 2% agarose gel stained with GelRed nucleic acid dye (Biotium) before sequencing. Sanger sequencing was used for SNP genotyping, completed by Sangon Biotech (Wuhan) Co., Ltd.

The levels of glucose were measured by the glucose oxidase method on biochemical analyzers of AU680 and AU5800 (Beckman Coulter, America). The 75 g OGTT-based fasting insulin (FINS), 1-h insulin (1hINS) and 2-h insulin (2hINS) were measured by enzyme-linked immunosorbent assay (ELISA) using Insulin ELISA BS-0783 (Bensheng (Tianjin) Health Technology Co., Ltd., China). The composite (Matsuda) insulin sensitivity index (ISIMatsuda) was calculated to assess insulin sensitivity. The 60 min insulinogenic index (IGI60) and OGTT-derived disposition index (DI) which have been reported to be more sensitive than the conventional index—homoeostasis model assessment of β cell function (HOMA-β) were calculated to dynamically estimate β cell function [27]. ISIMatsuda was calculated as 10,000/[(FPG × FINS)1/2 (GLUmean × INSmean)1/2] using fasting, 1-h, and 2-h values of OGTT. IGI60 was calculated as (1hINS-FINS)/(1hPG-FPG) and we estimated DI by multiplying IGI60 with ISIMatsuda to reflect β cell function adjusting for insulin sensitivity.

Ethical approval

All subjects provided informed consent for participation, and the study was approved by the ethics committees of Jingzhou Hospital Affiliated to Yangtze University and Gongan County Maternal and Child Health Care Hospital. This study followed the principles of the Declaration of Helsinki.

Statistical analysis

Normally distributed data are presented as the mean ± SD. Skewed distribution data are displayed as the median (interquartile range). P–P plots were used for normality tests, and the Levene test was used for homogeneity of variance tests. Two-sample t test for continuous variables conforming to normal distribution, Mann‒Whitney U test for continuous variables with nonnormal distribution and Chi-square test or Mann‒Whitney U test for categorical variables were used to evaluate the difference in means and proportions between case and control group. ANOVA and Kruskal‒Wallis tests were used to assess the difference in means of metabolic parameters between different genotypes. Multiple comparisons were conducted by the Bonferroni test. Logistic regression was performed to evaluate the association of GLP-1R gene polymorphisms and GDM risk. p < 0.05 was accepted as statistically significant. Analyses were performed with SPSS 25.0 software. Hardy–Weinberg equilibrium (HWE) was evaluated using the SNPStats program (http://bioinfo.iconcologia.net/SNPstats). LD analysis and selection of tag SNPs were calculated using HaploView 4.2 software.

Results

Baseline characteristics of participants

To eliminate the effects of confounding factors and keep the baseline characteristics of the GDM and NGT groups as consistent as possible, we matched the two groups according to age, basal BMI, hospital and hospital visiting date of OGTT. Moreover, we collected as much information about the subjects as possible, such as sociodemographic and clinical characteristics, to compare and identify potential confounders. As shown in Table 2, the GDM and NGT groups had comparable age, height, weight, first-trimester BMI, diastolic blood pressure and gestational age at enrollment. There was also no significant difference between the two groups in the component ratios of occupation, education level, history of pregnancy and parturition, domicile address and nature. However, systolic blood pressure in the GDM group was significantly higher than that in the NGT group (p = 0.006), which would be considered a confounder to adjust for covariates in subsequent analyses. Moreover, although age and first-trimester BMI showed no difference between the two groups, they were also used as confounders to adjust for covariates to eliminate their weak influence.

Table 2 Baseline characteristics of the study population

Relationship between GLP-1R polymorphisms and GDM risk

The genotype distributions of the 5 tag SNPs showed no deviation from HWE in cases and controls. Table 3 shows the genotype frequency distribution of the 5 tag SNPs in the GDM and NGT groups. We analyzed five genetic models for each SNP. Logistic regression in Table 3 showed that rs6458093 was associated with an approximately 1.54-fold increased risk of GDM under the dominant model (p = 0.044), and this effect remained nominally significant after adjusting for age, first-trimester BMI and systolic blood pressure (p = 0.049).

Table 3 Association between GLP-1R gene polymorphisms and gestational diabetes mellitus susceptibility

Age and BMI are recognized as the most important environmental risk factors for GDM, so we performed analyses stratified by age and BMI. As shown in Table 4, genotype GA of rs3765467 decreased GDM risk among subjects older than 35 years old after adjusting for age, first-trimester BMI and systolic blood pressure (p = 0.037). Genotype AG of rs6458093 was associated with increased GDM risk only among pregnant women younger than 35 years old after adjusting for age, first-trimester BMI and systolic blood pressure (p = 0.029). Moreover, rs6458093 was associated with increased GDM risk under the dominant model (AG + GG vs. AA) among pregnant women younger than 35 years old (p = 0.024) and with BMI no less than 24 (p = 0.041) after adjusting for confounders.

Table 4 Association between GLP-1R polymorphisms and risk of GDM after stratification by age and BMI

Effects of GLP-1R polymorphisms on glucose metabolism and beta cell function

We compared and estimated metabolic parameters of different genotypes, including FPG, 1hPG, 2hPG, FINS, 1hINS, 2hINS, ISImatsuda, IGI60 and DI. As shown in Fig. 2, genotype AG + GG of rs6458093 (A > G) showed significantly lower IGI60 and DI than genotype AA, as well as significantly higher 1hPG. No significant difference in metabolic parameters was found among different genotypes of rs3765467 (Fig. 3).

Fig. 2
figure 2

Effects of rs6458093 on glucose metabolism. A Effects of rs6458093 on glucose parameters. B Effects of rs6458093 on insulin parameters. C Effects of rs6458093 on insulin sensitivity and beta cell function. FPG: fasting plasma glucose; 1hPG: 1-h plasma glucose; 2hPG: 2-h plasma glucose; FINS: fasting insulin; 1hINS: 1-h insulin; 2hINS: 2-h insulin; ISImatsuda: Matsuda insulin sensitivity index; DI: disposition index; IGI60: 60 min insulinogenic index

Fig. 3
figure 3

Effects of rs3765467 on glucose metabolism. A Effects of rs3765467 on glucose parameters. B Effects of rs3765467 on insulin parameters. C Effects of rs3765467 on insulin sensitivity and beta cell function. FPG: fasting plasma glucose; 1hPG: 1-h plasma glucose; 2hPG: 2-h plasma glucose; FINS: fasting insulin; 1hINS: 1-h insulin; 2hINS: 2-h insulin; ISImatsuda: Matsuda insulin sensitivity index; DI: disposition index; IGI60: 60 min insulinogenic index

Discussion

The prevalence of GDM continues to increase steadily [28]. However, the etiology and pathogenesis of GDM remain vague due to its complexity under the combined action of environmental and genetic factors. Currently, environmental risk factors for GDM, such as age, BMI [29], sleep disorders [30], and polycystic ovary syndrome [31], have been well explained, but genetic risk factors still have not been well defined.

In this study, we conducted a two-center nested case‒control study based on the important role of GLP-1R in glucose homeostasis [32, 33] aimed to explore genetic risk factors for GDM. GLP-1R is a kind of G protein-coupled receptors [34]. When GLP-1 binds to GLP-1R, adenylate cyclase (AC) is activated by G protein to increase the intracellular concentration of cAMP. On the one hand, increased cAMP leads to the closure of K+ channels on the cytomembrane, depolarization of cell, opening of voltage-dependent Ca2+ channels, influx of extracellular Ca2+, increase of intracellular Ca2+ concentration, and synthesis and release of insulin finally. On the other hand, the increased cAMP activates protein kinase A (PKA) to phosphorylate related proteins, which further stimulates the transcription and translation of insulin gene finally [32, 34,35,36]. It was reported that very little incretin effect remained and glucose tolerance was pathologically impaired in double incretin receptors knockout mice [37, 38], illustrating that GLP-1R dysfunction affected the incretin effect producing insulin. Besides stimulating insulin synthesis and secretion, GLP-1R also mediates GLP-1 to stimulate β cell proliferation and neogenesis and inhibit β cell apoptosis [39]. Furthermore, abundant evidence suggests that SNPs of GLP-1R, such as rs10305420, rs3765467 and rs6923761, lead to GLP-1R dysfunction, subsequently affecting the binding of GLP-1 or GLP-1R agonists (GLP-1RAs) to GLP-1R, and ultimately affecting the hypoglycemic effects of GLP-1RAs and inhibitor of dipepdityl peptidase-4 (DPP-4) (the enzyme responsible for incretin degradation) [40,41,42]. In the same way, SNPs of GLP-1R impaired receptor function, subsequently affecting β cell insulin secretion function and ultimately affecting the efficacy of Roux-en-Y gastric bypass [15]. The above studies indicate that GLP-1R gene polymorphism is related to β cell insulin secretion function, which may affect glucose metabolism and susceptibility of GDM.

In our study, we identified two GLP-1R tag SNPs, rs6458093 and rs3765467, associated with GDM susceptibility, and rs6458093 also affected β cell function and postprandial glucose metabolism. The tag SNP rs6458093 was located in the intron region that may influence gene function by changing the level, location, stability, or timing of gene expression [43], directly influencing the quality or quantity of the gene product GLP-1R. Tag SNP rs3765467 is a missense mutation resulting in substitution of glutamine for arginine at position 131, directly influencing the function of GLP-1R [14]. As mentioned above, rs6458093 and rs3765467 might affect the interaction of GLP-1 and GLP-1R by affecting the quality or quantity of GLP-1R, then affecting β cell insulin secretion function, and ultimately altering glucose metabolism and susceptibility of GDM.

As rare missense variant SNPs with a high frequency of MAF in GLP-1R, rs1042044 and rs3765467 are often the focus of research. Although rs1042044 had nothing to do with insulin secretion [15], it was reported to be associated with susceptibility of T2DM in a Chinese literature [19] and the risk of papillary thyroid cancer among the Egyptian population in a recent study [44]. Moreover, rs1042044 was also differentially associated with brain functional connectivity in individuals with low versus high severity of alcohol use [45]. Yapici-Eser et al. [46] suggested a possible association of rs1042044 with anhedonia but no association with depression diagnosis. However, a correlation between rs1042044 and GDM was not found in our study. For rs3765467, Korean patients with T2DM carrying mutant allele A (GA/AA) showed a significantly better hypoglycemic effect to DPP-4 inhibitors than those with genotype GG [47], indicating that mutant genotype GA/AA could improve the function of GLP-1R mediating insulin secretion, while the major genotype GG might be a risk factor for disease. Research from the Mayo Clinic also suggested that subjects with genotype GA of rs3765467 exhibited a 100% increase in the total insulin secretion response to infused GLP-1 in the presence of hyperglycemia compared with subjects carrying the GG genotype [14]. Nishiya et al. [48] also found that the GG genotype was a significant risk factor for decreased insulin secretion. Moreover, rs3765467 was reported to be associated with T2DM risk [16]. Although no evidence in our study showed that rs3765467 affected β cell function and glucose metabolism, it was found that the GA genotype of rs3765467 was related to a decreased risk of GDM among pregnant women older than 35 years. This inconsistent result might be explained by the fact that rs3765467 was mainly related to a decreased risk of GDM among pregnant women older than 35 years. Regrettably, in our study, it was not feasible to compare metabolic parameters of different genotypes stratified by age due to the small sample size of subjects older than 35 years.

To our knowledge, rs6458093, rs10305478 and rs3765468 have not been reported to be related to disease or insulin secretion, so we know little about their role. However, our study found that rs6458093 was not only associated with GDM susceptibility but also affected postprandial glucose homeostasis and β cell function. The findings again implied that genetic variations in GLP-1R may affect susceptibility to GDM by affecting β cell function, especially postprandial insulin secretion and glucose homeostasis, which could be explained by the fact that incretin hormones mainly mediate postprandial insulin secretion and glucose homeostasis [49].

With a high MAF in the European population, the missense variant SNP rs6923761 is the most popular study locus. It was reported to be associated with a glucose response to a mixed meal [50], a smaller glycemic response to treatment with DPP-4 inhibitors [42, 51], delayed gastric-emptying in obese patients treated with liraglutide (GLP-1RAs) or exenatide (GLP-1 analog) [52], weight loss and metabolic change in diabetic patients treated with liraglutide [53] and so on. Despite many studies supporting the hypothesis that rs6923761 may be associated with GDM susceptibility and glucose metabolism, rs6923761 was not included in our study due to the low MAF of < 0.01 in the East Asian population. Regrettably, the missense variant SNP rs10305420 was not included in our study because it was localized in a fragment with a low SNP density. Future studies exploring the relationship between rs10305420 and GDM susceptibility are reasonable, as rs10305420 was reported to be associated with reduced glycemic response and less weight loss response to GLP-1R agonists [54].

Furthermore, our findings also have some enlightenment for research on T2DM. Kwak et al. even suggested the identification of genetic risk loci in GDM as a complementary approach to improve our insight into the genetics of T2DM given the similar genetic basis between GDM and T2DM [55]. It remains to be determined whether rs6458093 and rs3765467 also play roles in the development of T2DM. Moreover, GLP-1-based therapeutics have been rapidly developed and used clinically for the treatment of T2DM. Despite many advantages, individual differences in treatment outcomes are obvious. Some SNPs in GLP-1R, such as rs6923761 [42, 50], rs10305420 [41], and rs3765467 [40], have been reported to lead to individual differences in these therapeutics. The strong association between GLP-1R gene polymorphisms and the hypoglycemic effects of DPP-4 inhibitors and GLP-1RAs makes the GLP-1R gene a candidate gene for precision medicine in diabetes [56]. Future studies are needed to determine whether rs6458093 accounts for interindividual differences in response to GLP-1-based therapy.

Our study also has some limitations. First, although there is a large population in our country and there are 1018 pregnant women in our cohort, the sample size in our study was limited to some extent due to our study design (strict inclusion and exclusion criteria, and strict 1:1 frequency matching between the case group and control group based on a nested case‒control study aimed to minimize the influence of confounding factors as much as possible), so subgroup analyses and stratification analyses of metabolic parameters for rs3765467 could not be performed. Second, lifestyle was not considered in our study since GDM is a lifestyle-related disease. Third, none of the reported p values in this study were corrected for multiple testing. The p values of < 0.01 were significant after Bonferroni correction for 5 SNPs. Fourth, we can only use some indices calculated to estimate β cell function due to the particular period of pregnancy, while the gold standard for assessing β cell function remains the glucose clamp study. Each calculated index has some limitations.

In conclusion, our study identified two tag SNPs of GLP-1R nominally associated with GDM susceptibility, and one of them affected beta cell function and postprandial glucose metabolism, providing evidence for the etiology and genetic study of GDM and some inspirations for T2DM.

Availability of data and materials

The datasets used in this study are available from the corresponding author on reasonable request.

Abbreviations

GDM:

Gestational diabetes mellitus

T2DM:

Type 2 diabetes mellitus

β cell:

Beta cell

NGT:

Normal glucose tolerance

GLP-1:

Glucagon-like peptide-1

GLP-1R:

Glucagon-like peptide-1 receptor

SNP:

Single nucleotide polymorphism

MAF:

Minor allele frequency

LD:

Linkage disequilibrium

OGTT:

Oral glucose tolerance test

FPG:

Fasting plasma glucose

1hPG:

75 g OGTT-based 1-h plasma glucose

2hPG:

75 g OGTT-based 2-h plasma glucose

FINS:

Fasting insulin

1hINS:

75 g OGTT-based 1-h insulin

2hINS:

75 g OGTT-based 2-h insulin

ELISA:

Enzyme-linked immunosorbent assay

ISIMatsuda :

Matsuda insulin sensitivity index

IGI60 :

Insulinogenic index at 60 min

DI:

75 g OGTT-derived disposition index

HWE:

Hardy–Weinberg equilibrium

BMI:

Body mass index

GLP-1RAs:

Glucagon-like peptide-1 receptor agonists

DPP-4:

Dipepdityl peptidase-4

GHbA1c:

Glycosylated hemoglobin a1c

References

  1. Wang H, Li N, Chivese T, et al. IDF diabetes atlas: estimation of global and regional gestational diabetes mellitus prevalence for 2021 by international association of diabetes in pregnancy study group’s criteria. Diabetes Res Clin Pract. 2022;183:109050. https://doi.org/10.1016/j.diabres.2021.109050.

    Article  Google Scholar 

  2. Zhu H, Zhao Z, Xu J, et al. The prevalence of gestational diabetes mellitus before and after the implementation of the universal two-child policy in China. Front Endocrinol (Lausanne). 2022;13:960877. https://doi.org/10.3389/fendo.2022.960877.

  3. Gorban de Lapertosa S, Sucani S, Salzberg S, et al. Prevalence of gestational diabetes mellitus in Argentina according to the Latin American Diabetes Association (ALAD) and International Association of Diabetes and Pregnancy Study Groups (IADPSG) diagnostic criteria and the associated maternal-neonatal complications. Health Care Women Int. 2021;42(4–6):636–656. https://doi.org/10.1080/07399332.2020.1800012.

  4. Vounzoulaki E, Khunti K, Abner SC, Tan BK, Davies MJ, Gillies CL. Progression to type 2 diabetes in women with a known history of gestational diabetes: systematic review and meta-analysis. BMJ. 2020;369:m1361. Published 2020 May 13. https://doi.org/10.1136/bmj.m1361.

  5. Li J, Song C, Li C, Liu P, Sun Z, Yang X. Increased risk of cardiovascular disease in women with prior gestational diabetes: a systematic review and meta-analysis. Diabetes Res Clin Pract. 2018;140:324–38. https://doi.org/10.1016/j.diabres.2018.03.054.

    Article  Google Scholar 

  6. Faith MS, Hittner JB, Hurston SR, et al. Association of infant temperament with subsequent obesity in young children of mothers with gestational diabetes mellitus. JAMA Pediatr. 2019;173(5):424–33. https://doi.org/10.1001/jamapediatrics.2018.5199.

    Article  Google Scholar 

  7. Clausen TD, Mathiesen ER, Hansen T, et al. High prevalence of type 2 diabetes and pre-diabetes in adult offspring of women with gestational diabetes mellitus or type 1 diabetes: the role of intrauterine hyperglycemia. Diabetes Care. 2008;31(2):340–6. https://doi.org/10.2337/dc07-1596.

    Article  Google Scholar 

  8. McIntyre HD, Catalano P, Zhang C, Desoye G, Mathiesen ER, Damm P. Gestational diabetes mellitus. Nat Rev Dis Primers. 2019;5(1):47. https://doi.org/10.1038/s41572-019-0098-8.

    Article  Google Scholar 

  9. Kampmann U, Knorr S, Fuglsang J, Ovesen P. Determinants of maternal insulin resistance during pregnancy: an updated overview. J Diabetes Res. 2019;2019:5320156. https://doi.org/10.1155/2019/5320156.

    Article  CAS  Google Scholar 

  10. Elrick H, Stimmler L, Hlad CJ Jr, Arai Y. Plasma insulin response to oral and intravenous glucose administration. J Clin Endocrinol Metab. 1964;24:1076–1082. https://doi.org/10.1210/jcem-24-10-1076.

  11. Hampton SM, Morgan LM, Tredger JA, Cramb R, Marks V. Insulin and C-peptide levels after oral and intravenous glucose. Contribution of enteroinsular axis to insulin secretion. Diabetes. https://doi.org/10.2337/diab.35.5.612.

  12. Nauck MA, Homberger E, Siegel EG, et al. Incretin effects of increasing glucose loads in man calculated from venous insulin and C-peptide responses. J Clin Endocrinol Metab. 1986;63(2):492–8. https://doi.org/10.1210/jcem-63-2-492.

    Article  CAS  Google Scholar 

  13. Kreymann B, Williams G, Ghatei MA, Bloom SR. Glucagon-like peptide-1 7–36: a physiological incretin in man. Lancet. 1987;2(8571):1300–4. https://doi.org/10.1016/s0140-6736(87)91194-9.

    Article  CAS  Google Scholar 

  14. Sathananthan A, Man CD, Micheletto F, et al. Common genetic variation in GLP1R and insulin secretion in response to exogenous GLP-1 in nondiabetic subjects: a pilot study. Diabetes Care. 2010;33(9):2074–6. https://doi.org/10.2337/dc10-0200.

    Article  CAS  Google Scholar 

  15. Li W, Li P, Li R, et al. GLP1R single-nucleotide polymorphisms rs3765467 and rs10305492 affect β cell insulin secretory capacity and apoptosis through GLP-1. DNA Cell Biol. 2020;39(9):1700–10. https://doi.org/10.1089/dna.2020.5424.

    Article  CAS  Google Scholar 

  16. Suzuki K, Akiyama M, Ishigaki K, et al. Identification of 28 new susceptibility loci for type 2 diabetes in the Japanese population. Nat Genet. 2019;51(3):379–86. https://doi.org/10.1038/s41588-018-0332-4.

    Article  CAS  Google Scholar 

  17. Wessel J, Chu AY, Willems SM, et al. Low-frequency and rare exome chip variants associate with fasting glucose and type 2 diabetes susceptibility. Nat Commun. 2015;6:5897. https://doi.org/10.1038/ncomms6897.

    Article  CAS  Google Scholar 

  18. Scott RA, Freitag DF, Li L, et al. A genomic approach to therapeutic target validation identifies a glucose-lowering GLP1R variant protective for coronary heart disease. Sci Transl Med. 2016;8(341):341ra76. https://doi.org/10.1126/scitranslmed.aad3744.

  19. Fang B, Luo J, Li Y, Li X, Yang Y. Association between GLP1R rs1042044 gene polymorphism and type 2 diabetes mellitus. J Xinjiang Med Univ. 2021;44(01):8–11.

    Google Scholar 

  20. Tokuyama Y, Matsui K, Egashira T, Nozaki O, Ishizuka T, Kanatsuka A. Five missense mutations in glucagon-like peptide 1 receptor gene in Japanese population. Diabetes Res Clin Pract. 2004;66(1):63–9. https://doi.org/10.1016/j.diabres.2004.02.004.

    Article  CAS  Google Scholar 

  21. Plows JF, Stanley JL, Baker PN, Reynolds CM, Vickers MH. The pathophysiology of gestational diabetes mellitus. Int J Mol Sci. 2018;19(11):3342. https://doi.org/10.3390/ijms19113342.

    Article  CAS  Google Scholar 

  22. Pervjakova N, Moen GH, Borges MC, et al. Multi-ancestry genome-wide association study of gestational diabetes mellitus highlights genetic links with type 2 diabetes [published online ahead of print, 2022 Feb 26]. Hum Mol Genet. 2022;ddac050. https://doi.org/10.1093/hmg/ddac050.

  23. Dennison RA, Chen ES, Green ME, et al. The absolute and relative risk of type 2 diabetes after gestational diabetes: a systematic review and meta-analysis of 129 studies. Diabetes Res Clin Pract. 2021;171:108625. https://doi.org/10.1016/j.diabres.2020.108625.

    Article  Google Scholar 

  24. Li Z, Cheng Y, Wang D, et al. Incidence Rate of Type 2 Diabetes Mellitus after Gestational Diabetes Mellitus: A Systematic Review and Meta-Analysis of 170,139 Women. J Diabetes Res. 2020;2020:3076463. Published 2020 Apr 27. https://doi.org/10.1155/2020/3076463.

  25. Ma Y, Wang N, Gu L, et al. Postpartum assessment of the beta cell function and insulin resistance for Chinese women with previous gestational diabetes mellitus. Gynecol Endocrinol. 2019;35(2):174–8. https://doi.org/10.1080/09513590.2018.1512094.

    Article  CAS  Google Scholar 

  26. Coussens AK, van Daal A. Linkage disequilibrium analysis identifies an FGFR1 haplotype-tag SNP associated with normal variation in craniofacial shape. Genomics. 2005;85(5):563–73. https://doi.org/10.1016/j.ygeno.2005.02.002.

    Article  CAS  Google Scholar 

  27. Ohn JH, Kwak SH, Cho YM, et al. 10-year trajectory of β-cell function and insulin sensitivity in the development of type 2 diabetes: a community-based prospective cohort study. Lancet Diabetes Endocrinol. 2016;4(1):27–34. https://doi.org/10.1016/S2213-8587(15)00336-8.

    Article  CAS  Google Scholar 

  28. Yan B, Yu Y, Lin M, et al. High, but stable, trend in the prevalence of gestational diabetes mellitus: a population-based study in Xiamen, China. J Diabetes Investig. 2019;10(5):1358–64. https://doi.org/10.1111/jdi.13039.

    Article  CAS  Google Scholar 

  29. Leng J, Shao P, Zhang C, et al. Prevalence of gestational diabetes mellitus and its risk factors in Chinese pregnant women: a prospective population-based study in Tianjin, China. PLoS ONE. 2015;10(3):e0121029. Published 2015 Mar 23. https://doi.org/10.1371/journal.pone.0121029.

  30. Zhang X, Zhang R, Cheng L, et al. The effect of sleep impairment on gestational diabetes mellitus: a systematic review and meta-analysis of cohort studies. Sleep Med. 2020;74:267–77. https://doi.org/10.1016/j.sleep.2020.05.014.

    Article  Google Scholar 

  31. Qiu Y, Zhang X, Ni Y. Association between polycystic ovarian syndrome and risk of gestational diabetes mellitus: a meta-analysis. Gynecol Obstet Invest. 2022;87(2):150–8. https://doi.org/10.1159/000521728.

    Article  CAS  Google Scholar 

  32. Perfetti R, Merkel P. Glucagon-like peptide-1: a major regulator of pancreatic beta-cell function. Eur J Endocrinol. 2000;143(6):717–725. https://doi.org/10.1530/eje.0.1430717.

  33. Hayes MR, De Jonghe BC, Kanoski SE. Role of the glucagon-like-peptide-1 receptor in the control of energy balance. Physiol Behav. 2010;100(5):503–510. https://doi.org/10.1016/j.physbeh.2010.02.029.

  34. Fletcher MM, Halls ML, Christopoulos A, Sexton PM, Wootten D. The complexity of signalling mediated by the glucagon-like peptide-1 receptor. Biochem Soc Trans. 2016;44(2):582–588. https://doi.org/10.1042/BST20150244.

  35. Gomez E, Pritchard C, Herbert TP. cAMP-dependent protein kinase and Ca2+ influx through L-type voltage-gated calcium channels mediate Raf-independent activation of extracellular regulated kinase in response to glucagon-like peptide-1 in pancreatic beta-cells. J Biol Chem. 2002;277(50):48146–48151. https://doi.org/10.1074/jbc.M209165200.

  36. Kang G, Chepurny OG, Rindler MJ, et al. A cAMP and Ca2+ coincidence detector in support of Ca2+-induced Ca2+ release in mouse pancreatic beta cells. J Physiol. 2005;566(Pt 1):173–188. https://doi.org/10.1113/jphysiol.2005.087510.

  37. Ahrén B, Yamada Y, Seino Y. The Incretin Effect in Female Mice With Double Deletion of GLP-1 and GIP Receptors. J Endocr Soc. 2019;4(2):bvz036. Published 2019 Dec 23. https://doi.org/10.1210/jendso/bvz036.

  38. Hansotia T, Baggio LL, Delmeire D, et al. Double incretin receptor knockout (DIRKO) mice reveal an essential role for the enteroinsular axis in transducing the glucoregulatory actions of DPP-IV inhibitors. Diabetes. 2004;53(5):1326–35. https://doi.org/10.2337/diabetes.53.5.1326.

    Article  CAS  Google Scholar 

  39. Campbell JE, Drucker DJ. Pharmacology, physiology, and mechanisms of incretin hormone action. Cell Metab. 2013;17(6):819–37. https://doi.org/10.1016/j.cmet.2013.04.008.

    Article  CAS  Google Scholar 

  40. Guan Z, Du Y, Li R, et al. Association between glucagon-like peptide-1 receptor gene polymorphism and treatment response to GLP1R agonists in Chinese patients with type 2 diabetes: a prospective cohort study. Eur J Clin Pharmacol. 2022;78(5):793–9. https://doi.org/10.1007/s00228-021-03249-z.

    Article  CAS  Google Scholar 

  41. Jensterle M, Pirš B, Goričar K, Dolžan V, Janež A. Genetic variability in GLP-1 receptor is associated with inter-individual differences in weight lowering potential of liraglutide in obese women with PCOS: a pilot study. Eur J Clin Pharmacol. 2015;71(7):817–24. https://doi.org/10.1007/s00228-015-1868-1.

    Article  CAS  Google Scholar 

  42. Űrgeová A, Javorský M, Klimčáková L, et al. Genetic variants associated with glycemic response to treatment with dipeptidylpeptidase 4 inhibitors. Pharmacogenomics. 2020;21(5):317–23. https://doi.org/10.2217/pgs-2019-0147.

    Article  CAS  Google Scholar 

  43. Liu K, Huo H, Jia W, et al. RAB40C gene polymorphisms rs62030917 and rs2269556 are associated with an increased risk of lumbar disc herniation development in the Chinese Han population. J Gene Med. 2021;23(4):e3252. https://doi.org/10.1002/jgm.3252.

    Article  CAS  Google Scholar 

  44. Abdul-Maksoud RS, Elsayed WSH, Rashad NM, Elsayed RS, Elshorbagy S, Hamed MG. GLP-1R polymorphism (rs1042044) and expression are associated with the risk of papillary thyroid cancer among the Egyptian population. Gene. 2022;834:146597. https://doi.org/10.1016/j.gene.2022.146597.

    Article  CAS  Google Scholar 

  45. Farokhnia M, Fede SJ, Grodin EN, et al. Differential association between the GLP1R gene variants and brain functional connectivity according to the severity of alcohol use. Sci Rep. 2022;12(1):13027. https://doi.org/10.1038/s41598-022-17190-3.

    Article  CAS  Google Scholar 

  46. Yapici-Eser H, Appadurai V, Eren CY, et al. Association between GLP-1 receptor gene polymorphisms with reward learning, anhedonia and depression diagnosis. Acta Neuropsychiatr. 2020;32(4):218–25. https://doi.org/10.1017/neu.2020.14.

    Article  Google Scholar 

  47. Han E, Park HS, Kwon O, et al. A genetic variant in GLP1R is associated with response to DPP-4 inhibitors in patients with type 2 diabetes. Medicine (Baltimore). 2016;95(44):e5155. https://doi.org/10.1097/MD.0000000000005155.

    Article  CAS  Google Scholar 

  48. Nishiya Y, Daimon M, Mizushiri S, et al. Nutrient consumption-dependent association of a glucagon-like peptide-1 receptor gene polymorphism with insulin secretion. Sci Rep. 2020;10(1):16382. Published 2020 Oct 2. https://doi.org/10.1038/s41598-020-71853-7.

  49. Holst JJ, Gasbjerg LS, Rosenkilde MM. The Role of Incretins on Insulin Function and Glucose Homeostasis. Endocrinology. 2021;162(7):bqab065. https://doi.org/10.1210/endocr/bqab065.

  50. Mashayekhi M, Wilson JR, Jafarian-Kerman S, et al. Association of a glucagon-like peptide-1 receptor gene variant with glucose response to a mixed meal. Diabetes Obes Metab. 2021;23(1):281–6. https://doi.org/10.1111/dom.14216.

    Article  CAS  Google Scholar 

  51. Javorský M, Gotthardová I, Klimčáková L, et al. A missense variant in GLP1R gene is associated with the glycaemic response to treatment with gliptins. Diabetes Obes Metab. 2016;18(9):941–4. https://doi.org/10.1111/dom.12682.

    Article  CAS  Google Scholar 

  52. Chedid V, Vijayvargiya P, Carlson P, et al. Allelic variant in the glucagon-like peptide 1 receptor gene associated with greater effect of liraglutide and exenatide on gastric emptying: a pilot pharmacogenetics study. Neurogastroenterol Motil. 2018;30(7):e13313. https://doi.org/10.1111/nmo.13313.

    Article  CAS  Google Scholar 

  53. de Luis DA, Diaz Soto G, Izaola O, Romero E. Evaluation of weight loss and metabolic changes in diabetic patients treated with liraglutide, effect of RS 6923761 gene variant of glucagon-like peptide 1 receptor. J Diabetes ications. 2015;29(4):595–8. https://doi.org/10.1016/j.jdiacomp.2015.02.010.

    Article  Google Scholar 

  54. Yu M, Wang K, Liu H, Cao R. GLP1R variant is associated with response to exenatide in overweight Chinese Type 2 diabetes patients. Pharmacogenomics. 2019;20(4):273–7. https://doi.org/10.2217/pgs-2018-0159.

    Article  CAS  Google Scholar 

  55. Kwak SH, Kim SH, Cho YM, et al. A genome-wide association study of gestational diabetes mellitus in Korean women. Diabetes. 2012;61(2):531–41. https://doi.org/10.2337/db11-1034.

    Article  CAS  Google Scholar 

  56. Zeng Z, Huang SY, Sun T. Pharmacogenomic studies of current antidiabetic agents and potential new drug targets for precision medicine of diabetes. Diabetes Ther. 2020;11(11):2521–38. https://doi.org/10.1007/s13300-020-00922-x.

    Article  Google Scholar 

Download references

Acknowledgements

None.

Funding

None.

Author information

Authors and Affiliations

Authors

Contributions

Study design: PL, BM. Data collection and analysis: PL, YF, YX, HF. Manuscript writing: PL, ZY, CZ. Manuscript revision: PL, BM. All authors read and approved the final version of the manuscript submitted for publication.

Corresponding author

Correspondence to Bing Mei.

Ethics declarations

Ethics approval and consent to participate

All subjects provided informed consent for participation, and the study was approved by the ethics committees of Jingzhou Hospital Affiliated to Yangtze University and Gongan County Maternal and Child Health Care Hospital. This study followed the principles of the Declaration of Helsinki.

Consent for publication

All authors gave consent for publication.

Competing interests

The authors have no relevant financial or nonfinancial interests to disclose.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Luo, P., Fan, Y., Xiong, Y. et al. Genetic variants of the GLP-1R gene affect the susceptibility and glucose metabolism of gestational diabetes mellitus: a two-center nested case‒control study. Diabetol Metab Syndr 14, 190 (2022). https://doi.org/10.1186/s13098-022-00963-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13098-022-00963-1

Keywords