The rs10229583 polymorphism near paired box gene 4 is associated with gestational diabetes mellitus in Chinese women

Abstract

Objective

The rs10229583 polymorphism near paired box gene 4 (PAX4) is associated with insulin resistance and type 2 diabetes. Mutations in the PAX4 gene may be associated with impaired differentiation/development of pancreatic islet beta cells during fetal development and, consequently, a compromised insulin response to high blood glucose. To ascertain whether this polymorphism plays a role in gestational diabetes mellitus (GDM), we investigated the genotypic and allele frequency differences between GDM and normal pregnancies.

Methods

A total of 310 GDM and 440 normal pregnancies were evaluated. Allele and genotype frequencies of rs10229583 were determined for all participants with Sanger sequencing and SNaPshot. Association of the allele and genotypes of the single nucleotide polymorphism with the disease was analyzed using Pearson’s χ² test and OR (odds ratio).

Results

The G allele was more frequent in patients with GDM compared with controls (OR = 1.47, 95% confidence interval (CI): 1.12–1.939). The GG genotype frequency of rs10229583 was significantly different between subjects with GDM and normal controls (OR = 1.411, 95% CI: 1.032–1.928). The OR of the GA + GG genotype was 3.182 (95% CI: 1.294–7.826) for patients with GDM compared with controls.

Conclusion

The present study suggests that rs10229583 is associated with GDM.

Keywords

Gestational diabetes mellitus polymorphism insulin resistance prenatal health mutation

Introduction

Gestational diabetes mellitus (GDM) is defined as glucose intolerance first diagnosed during pregnancy and is associated with increased risk of preeclampsia, perinatal mortality of the fetus and type 2 diabetes (T2D) development.¹ During pregnancy, women experience increased peripheral insulin resistance, which is explained, in part, by increased production of estrogen, human placental lactogen and prolactin.^2–4 To maintain glucose homeostasis, insulin secretion is increased. GDM occurs in cases where insulin secretion fails to compensate for insulin resistance. Women with a history of GDM are at a higher risk of developing T2D later in life.⁵ GDM and T2D are multifactorial disorders with a common pathogenesis. The diseases share several features, including pancreatic beta-cell dysfunction, insulin resistance and association with similar risk factors. A variety of risk factors, including obesity, family history, lifestyle, ethnicity, genetics and physical inactivity, contribute to the development of both GDM and T2D.⁶ Earlier studies strongly supported the possibility that GDM and T2D share the same genetic background. Over the past few decades, numerous studies have focused on the identification of accumulating genetic variants that are strongly associated with T2D.⁷

The rs10229583 polymorphism is located downstream of the paired box gene 4 (PAX4) gene, which is located on chromosome 7q32 and belongs to a family of transcription factors. The gene plays an important role in the differentiation and development of pancreatic islet beta cells during fetal development.^8,9 A number of investigations to date have identified rs10229583 near PAX4 as a locus for T2D.¹⁰ However, little is known about its potential association with GDM. Based on these findings, we hypothesized that the rs10229583 single nucleotide polymorphism (SNP) is associated with predisposition to GDM. Therefore, PAX4 may be a novel target candidate gene for GDM therapy.

In the current study, we examined the genetic, clinical and biochemical differences between GDM and normal pregnancies according to the Strengthening the Reporting of Genetic Association Studies initiative.¹¹ We aimed to clarify the relationship between rs10229583 and GDM as well as its association with clinical and biochemical parameters. Our results may provide better insights into the mechanisms underlying PAX4 gene variations in the development of GDM.

Materials and methods

Subjects

The 310 patients with GDM and 440 normal pregnancy subjects selected for study were of Han Chinese ethnic origin. Between 2010 and 2014, subjects were recruited from the International Peace Maternity and Child Health Hospital affiliated with the Medical School of Shanghai Jiao Tong University and Bao’an Maternal and Child Health Hospital, Shenzhen, China. No previous diagnosis of chronic hypertension, type 1 or type 2 diabetes, infection, renal dysfunction, rheumatoid arthritis or systemic lupus erythematosus was reported in any case. GDM diagnosis was performed according to the International Association of the Diabetes and Pregnancy Study Groups criteria,¹² which later was also implemented by the World Health Organization (http://apps.who.int/iris/bitstream/10665/85975/1/WHO_NMH_MND_13.2_eng.pdf). Our study protocol was in accordance with institutional ethic guidelines and national laws. Informed written consent was obtained from all subjects. This study was approved by the Ethics Committees of the International Peace Maternity and Child Health Hospital affiliated with the Medical School of Shanghai Jiao Tong University and Bao’an Maternal and Child Health Hospital and conducted in accordance with the principles of the Declaration of Helsinki.

Study protocol

All subjects underwent a 75-g oral glucose tolerance test between 24 and 28 weeks of gestation to screen for GDM. Subjects were examined in the morning after a 12-h overnight fast between 08:30 and 09:30. Plasma glucose levels were measured at 0, 1 and 2 h during the glucose challenge. GDM was diagnosed in cases where one of the glucose values met or exceeded the threshold. Body mass index (BMI) for all pregnancies before delivery was calculated as the ratio of weight (kg) to the square of height (m²).¹³ Large, appropriate and small for gestational age rates were recorded. Biochemical data, including plasma glucose, glycosylated hemoglobin A1c (HbA1c) and triglyceride levels, were additionally measured in all cases.

Genotyping

Genomic DNA samples were extracted from peripheral whole blood using the TIANamp Blood DNA kit (TIANGEN Biotech, Beijing, China). Genotyping was performed using SNaPshot (ABI Prism SNaPShot multiplex kit; Applied Biosystems) on an ABI PRISM 3730 genetic analyzer according to the manufacturer’s instructions and evaluated using ABI GeneScan software. Approximately 15% of the samples in our study were sequenced using an ABI PRISM 3730 genetic analyzer to reconfirm genotype. Reproducibility of data was 100%. The proportion of technically successful genotyping out of all samples assessed for genotype was 98%.

Primers for amplification were designed using Primer3 (http://frodo.wi.mit.edu/). All primers were synthesized by Sangon Biotech Co., Ltd. (Shanghai, China). Fragments of the variant were amplified from 100 ng DNA in a 20-µl reaction mixture containing 50 mmol MgCl₂, 50 mmol dNTP, 1 µM forward primer, 1 µM reverse primer and 0.5 U Taq polymerase (TaKaRa). A touchdown PCR program was used according to the manufacturer’s instructions. The primer sequences for PCR amplification were as follows: forward 5′-CAACAGTGCTGCCTATTCCC-3′ and reverse 5′-TAGGAAAGGAGATGCCCTCA-3′.

Statistical analysis

Statistical analyses were performed using the Statistical Package for the Social Sciences Windows version 17.0 (SPSS Inc., Chicago, IL, USA). Continuous variables were expressed as mean ± standard deviation and compared using the Student’s t-test. Association of the allele and genotype of the SNP with the disease was analyzed using Pearson’s χ² test and OR (odds ratio). Deviation from Hardy–Weinberg equilibrium was assessed using the Chi-square test. All tests were two-sided, and a value of p < 0.05 was accepted as statistically significant.

Results

Clinical and biochemical characteristics of subjects

The basic clinical and biochemical characteristics of all patients with GDM and controls are presented in Table 1. No significant differences in gestational age or neonatal birth weight were observed between subjects with GDM and normal controls. BMI, maternal age and concentrations of plasma glucose, HbA1c and triglycerides were significantly higher in subjects with GDM compared with those of normal controls (p < 0.05).

Table 1.

Demographic characteristics of the subjects

	Cases	Controls	P value
	(n = 310)	(n = 440)	P value
Age (mean ± sd, years)	30.63 ± 4.58	29.57 ± 4.43	0.002
0-h plasma glucose (mmol/l)	4.84 ± 0.77	4.33 ± 0.64	<0.001
1-h plasma glucose (mmol/l)	10.12 ± 1.81	7.23 ± 1.61	<0.001
2-h plasma glucose (mmol/l)	8.84 ± 1.64	6.53 ± 1.25	<0.001
Gestational weeks	38.13 ± 2.17	38.18 ± 2.68	0.761
HbA1c (%)	5.6 ± 0.42	5.25 ± 0.29	<0.001
Triglycerides (mmol/L)	2.87 ± 1.76	2.26 ± 1.17	<0.001
BMI (kg/m²)	27.01 ± 3.3	26.4 ± 3.25	0.025
Large gestational age	45 (14.52%)	65 (14.77%)
Appropriate gestational age	256 (82.58%)	355 (80.68%)	0.106
Small gestational age	9 (2.90%)	20 (4.55%)

HbA1c: hemoglobin A1c; BMI: body mass index

Genotyping and allelic frequencies

Genotypes were in Hardy–Weinberg equilibrium for all subjects. The distributions of genotypes and allelic frequencies are presented in Table 2. The genotypes of rs10229583 were AA, GA and GG, with frequencies of 1.93% (6/310), 27.10% (84/310) and 70.96% (220/310) in subjects with GDM and 5.91% (26/440), 30.68% (135/440) and 63.41% (279/440) in normal controls, respectively. The A and G allele frequencies were 15.48% (96/620) and 84.52% (524/620) in patients with GDM and 21.25% (187/880) and 78.75% (693/880) in control subjects, respectively.

Table 2.

The rs10229583 polymorphism genotypes and the risk of developing GDM

Genotypes	Cases	Controls	χ²	P value	OR (95% CI)
Genotypes	n = 310 (%)	n = 440 (%)	χ²	P value	OR (95% CI)
GG (vs. GA)	70.96%	63.41%	2.048	0.152	1.267 (0.916–1.753)
GG (vs. AA)	70.96%	63.41%	7.897	0.005	3.417 (1.382–8.447)
GG + GA (vs. AA)	98.06%	94.09%	7.03	0.008	3.182 (1.294–7.826)
GG (vs. GA + AA)	70.96%	63.41%	4.666	0.031	1.411 (1.032–1.928)

CI: confidence interval

The G allele was associated with GDM risk (OR = 1.47, 95% CI: 1.12–1.939, p < 0.005) (Table 2). The OR of the GA + GG genotype was 3.182 (95% CI: 1.294–7.826, p = 0.008), and that of the GG genotype was 1.411 (95% CI: 1.032–1.928, p = 0.031) for patients with GDM compared with controls. These data indicate that the G allele is a risk factor for GDM development.

Association of rs10229583 with clinical and biochemical characteristics in subjects with GDM

Association analysis of the rs10229583 SNP with clinical and biochemical characteristics in subjects with GDM revealed no significant differences between the AA, GA and GG genotypes in terms of BMI, large/appropriate/small for gestational age status or plasma glucose, HbA1c or triglyceride levels (Table 3).

Table 3.

Association of SNP rs10229583 with clinical and biochemical characteristics in subjects with GDM

Genotype	Cases (n)	0-h plasma glucose (mmol/l)	1-h plasma glucose (mmol/l)	2-h plasma glucose (mmol/l)	Gestational weeks	HbA1c (%)	Triglycerides (mmol/L)	BMI (kg/m²)	Large gestational age (n = 45)	Appropriate gestational age (n = 256)	Small gestational age (n = 9)
AA + GA	90	4.85 ± 0.7	10.12 ± 1.87	8.95 ± 1.7	38.05 ± 2.44	5.65 ± 0.39	3.04 ± 1.85	27.03 ± 3.63	11	76	3
GG	220	4.84 ± 0.8	10.12 ± 1.79	8.79 ± 1.62	38.16 ± 2.06	5.57 ± 0.43	2.81 ± 1.73	26.89 ± 3.56	34	180	6
t or χ²		0.139	−0.031	0.705	−0.362	1.209	0.871	0.293	0.594
P		0.889	0.975	0.482	0.718	0.228	0.385	0.770	0.743

HbA1c: hemoglobin A1c; BMI: body mass index

Discussion

GDM is a multifactorial disease caused by the complex interplay of multiple genetic, environmental and immunologic factors.^14–17 As GDM appears to share the same genetic background as T2D, candidate genes for GDM have been identified based on those associated with T2D. A number of genetic variants related to T2D are also reported to be associated with GDM.¹⁸ Several genome-wide association studies have suggested that mutations or polymorphisms in PAX4 are associated with both type 1 and 2 diabetes.^19–20 To clarify whether PAX4 is linked to GDM, we evaluated the possible association between the rs10229583 locus and GDM for the first time in this study.

Consistent with earlier findings, BMI and plasma glucose, HbA1c and triglyceride concentrations were significantly higher in subjects with GDM compared with those of normal controls. However, neonatal birth weight did not differ between subjects with GDM and normal controls. All pregnant women diagnosed with GDM in the clinic underwent diet and exercise therapy to ensure that weight before delivery and newborn birth weight were effectively controlled. Pharmacological therapy was also indicated in some cases to provide appropriate blood glucose control.

Our data revealed a significant association of the rs10229583 SNP with GDM. In particular, the G allele may be a risk factor for development of the disease. The relative risk for GA + GG was 3.182 (95% CI: 1.294–7.826, p = 0.008), while that for GG was 1.411 (95% CI: 1.032–1.928, p = 0.031) for patients with GDM. The rs10229583 polymorphism is located downstream of PAX4. While several studies to date have focused on the effects of PAX4 in the pathogenesis of T2D, to the best of our knowledge, none have addressed its potential association with GDM. PAX4 is a transcription factor belonging to the PAX family of proteins that contains a paired box domain, a paired-type homeodomain and an octapeptide involved in fetal development and tumor growth. PAX4 is expressed in early pancreatic endocrine cells, which is necessary for differentiation, development, proliferation and insulin secretion of pancreatic islet beta cells. The protein stimulates the expression of PDX1, a key transcription factor for islet development.^21,22 Targeted disruption of PAX4 in mice has been shown to reduce beta-cell mass at birth. Shimajiri et al. reported that the R121W mutation in PAX4 is significantly associated with T2D.¹¹ Ma and colleagues identified an association of rs10229583 near PAX4 with elevated fasting plasma glucose and impaired beta-cell function in controls and earlier age at diagnosis in Chinese populations.²³ Moreover, the rs10229583 locus was associated with T2DM in a Mongolian population in China.²⁴

To determine the functional implications of the rs10229583 SNP in GDM, we performed an association analysis between these variants and biochemical characteristics. Notably, the GG genotype of rs10229583, the risk variant for GDM, was not associated with BMI, neonatal birth weight or plasma glucose, HbA1c or triglyceride levels. Further research is required to elucidate the precise mechanisms by which these loci contribute to GDM susceptibility.

In conclusion, rs10229583 near PAX4 has been identified as a novel locus associated with GDM in Chinese populations, providing further insights into the pathogenesis of this disease.

Footnotes

Acknowledgements

We are grateful to all the subjects for their ongoing participation.

Declaration of conflicting interests

The authors declare that there is no conflict of interest.

Funding

This work was supported by a grant from the Shenzhen Municipal Science and Technology Innovation Council (grant no. JCYJ20150403105513698).

References

Metzger

Buchanan

Coustan

et al.

Summary and recommendations of the Fifth International Workshop-Conference on Gestational Diabetes Mellitus. Diabetes Care 2007; 30: S251–S260.

Kim

Toyofuku

Lynn

et al.

Serotonin regulates pancreatic beta cell mass during pregnancy. Nat Med 2010; 16: 804–808.

Di Cianni

Miccoli

Volpe

et al.

Intermediate metabolism in normal pregnancy and in gestational diabetes. Diabetes Metab Res Rev 2003; 19: 259–270.

Karnik

Chen

McLean

et al.

Menin controls growth of pancreatic beta-cells in pregnant mice and promotes gestational diabetes mellitus. Science 2007; 318: 806–809.

Linné

Barkeling

Rössner

. Natural course of gestational diabetes mellitus: long term follow up of women in the SPAWN study. BJOG 2002; 109: 1227–1231.

Konig

Shuldiner

. The genetic interface between gestational diabetes and type 2 diabetes. J Matern Fetal Neonatal Med 2012; 25: 36–40.

Hindorff

Sethupathy

Junkins

et al.

Potential etiologic and functional implications of genome-wide association loci for human diseases and traits. Proc Natl Acad Sci USA 2009; 106: 9362–9367.

Brun

Franklin

St-Onge

et al.

The diabetes-linked transcription factor PAX4 promotes {beta}-cell proliferation and survival in rat and human islets. J Cell Biol 2004; 167: 1123–1135.

Nagai

Ohno

et al.

Aberrant DNA demethylation in promoter region and aberrant expression of mRNA of PAX4 gene in hematologic malignancies. Leuk Res 2006; 30: 1547–1553.

10.

Mauvais-Jarvis

Smith

Le May

et al.

PAX4 gene variations predispose to ketosis-prone diabetes. Hum Mol Genet 2004; 13: 3151–3159.

11.

Little

Higgins

Ioannidis

et al.

Strengthening the reporting of genetic association studies (STREGA): an extension of the strengthening the reporting of observational studies in epidemiology (STROBE) statement. J Clin Epidemiol 2009; 62: 597–608.

12.

Metzger

Gabbe

et al.

International Association of Diabetes and Pregnancy Study Groups Consensus Panel. International association of diabetes and pregnancy study groups recommendations on the diagnosis and classification of hyperglycemia in pregnancy. Diabetes Care 2010; 33: 676–682.

13.

Tarquini

Picchiassi

Centra

et al.

Body mass index associated to rs2021966 ENPP1 polymorphism increases the risk for gestational diabetes mellitus. Gynecol Endocrinol 2015; 31: 83–86.

14.

Catalano

Huston

Amini

et al.

Longitudinal changes in glucose metabolism during pregnancy in obese women with normal glucose tolerance and gestational diabetes mellitus. Am J Obstet Gynecol 1999; 180: 903–916.

15.

Kusunoki

Katsuno

Nakae

et al.

Insulin resistance and β-cell function influence postprandial blood glucose levels in Japanese patients with gestational diabetes mellitus. Gynecol Endocrinol 2015; 31: 929–933.

16.

Kwak

Kim

Cho

et al.

A genome-wide association study of gestational diabetes mellitus in Korean women. Diabetes 2012; 61: 531–541.

17.

Rosta

Al-Aissa

Hadarits

et al.

Association Study with 77 SNPs Confirms the Robust Role for the rs10830963/G of MTNR1B variant and identifies two novel associations in gestational diabetes mellitus development. PLoS One 2017; 12: e0169781–e0169781.

18.

Ayaz

Karakaş Çelik

Cayan

. The G1057D polymorphism of insulin receptor substrate-2 associated with gestational diabetes mellitus. Gynecol Endocrinol 2014; 30: 165–168.

19.

Tokuyama

Matsui

Ishizuka

et al.

The Arg121Trp variant in PAX4 gene is associated with beta-cell dysfunction in Japanese subjects with type 2 diabetes mellitus. Metabolism 2006; 55: 213–216.

20.

Biason-Lauber

Boehm

Lang-Muritano

et al.

Association of childhood type 1 diabetes mellitus with a variant of PAX4: possible link to beta cell regenerative capacity. Diabetologia 2005; 48: 900–905.

21.

Habener

Kemp

Thomas

et al.

Minireview: transcriptional regulation in pancreatic development. Endocrinology 2005; 146: 1025–1034.

22.

Wang

Elghazi

Parker

et al.

The concerted activities of Pax4 and Nkx2.2 are essential to initiate pancreatic beta-cell differentiation. Dev Biol 2004; 266: 178–189.

23.

Tam

Zhang

et al.

Genome-wide association study in a Chinese population identifies a susceptibility locus for type 2 diabetes at 7q32 near PAX4. Diabetologia 2013; 56: 1291–1305.

24.

Bai

Liu

Suyalatu

et al.

Association Analysis of Genetic Variants with Type 2 Diabetes in a Mongolian Population in China. J Diabetes Res 2015; 2015: 613236–613236.