Adiponectin gene polymorphisms and risk of gestational diabetes mellitus:A meta-analysis

Lin-Ting Huang,Shi-Lan Wu,Xin Liao,Shu-Juan Ma,Hong-Zhuan Tan

Abstract

Key words: Gestational diabetes mellitus; Single nucleotide; Polymorphism; Adiponectin;Gene; Meta-analysis

INTRUDUCTION

Gestational diabetes mellitus (GDM) is a condition of impaired glucose tolerance during pregnancy in women without a previous diagnosis of diabetes.It is associated with serious complications for both mother and child in the pre- and postnatal periods[1].Many kinds of risk factors contribute to GDM,such as ethnicity,genetics,family history,dietary habits,and physical activity[1].Obesity is a usual risk factor for GDM and can cause insulin resistance.Many biochemical mediators compounded in the adipose tissue and secreted in the circulatory system,such as resistin,adiponectin(ADIPOQ),and leptin,are suggested to correlate with obesity and insulin resistance[2].

ADIPOQ is produced in the adipose tissue and modulates various metabolic processes,including lipid metabolism,glucose and fatty acid oxidation.This hormone decreases insulin resistance,improves lipid metabolism,and exerts anti-inflammatory properties.Decreased plasma ADIPOQ levels were observed in patients with type 2 diabetes (T2D),metabolic syndrome,and obesity[3].During normal pregnancy,ADIPOQ levels progressively decline,with its plasma concentration reaching even lower in GDM women[1].Previous studies suggested thatADIPOQgene single nucleotide polymorphisms (SNPs) could influence the concentration of plasma ADIPOQ and subsequently insulin sensitivity[4-6].

Studies have paid attention to the SNPs+45T/Gin exon 2 and+276G/C inintron 2,-11391G/A,and-11377C/Gin the promoter region.The twoADIPOQlinkage disequilibrium blocks are where these four variants located within.Block 1 comprises the promoter sequence spanning the region -14811 to -4120,and block 2 encompasses the exons in the region -450 to +4545[2].The conclusions of these studies have been disputed regarding whether the metabolic phenotypes of GDM are influenced by the variability at this locus and which polymorphisms contribute to this effect.For example,Lowet al[7]reported that a significant association was found betweenSNP45T/Gand GDM,and normal patients with the TT genotype had significantly higher plasma ADIPOQ levels compared to those with the TG or GG genotype.Beltchevaet al[1]reported that-11377C/Gis associated with GDM.According to Daheret al[8],GDM is not associated with+45T/Gand-11377C/Gpolymorphisms.Reasons for the conflicting results are small sample sizes in a single study and the hereditary difference of ethnicity.

As the results are discrepant and+45T/Gwas the only polymorphism which participated in the meta-analysis,our study was meant to evaluate whether and to what extentADIPOQgene polymorphisms contribute to GDM.

MATERIALS AND METHODS

Literature search strategy

Electronic databases PubMed,EMBASE,Web of Science,WANFANG DATA,and CNKI were used to search possibly association articles on human genetic studies ofADIPOQand GDM that had been published up to 20 October 2018.The search terms used were:“Gestational diabetes mellitus” or “GDM” and “adiponectin” or“ADIPOQ” and “single nucleotide polymorphism” or “polymorphism”.

Selection criteria and data extraction

We searched the database and identified 87 articles.The selection criteria of the publications were as follow:(1) well-designed case control studies on genetic association of ADIPOQ and the risk of GDM; (2) clear diagnostic criteria for GDM; (3)independent and sufficient genotype data must be contained in the original papers and the data can calculate odds ratios (ORs) and 95% confidence intervals (CIs); and(4) there should be at least two articles that have studied each polymorphism that we used in our meta-analysis.During the selection,we removed 25 articles for duplicate publication and excluded 49 articles for review,animal studies,case reports with unrelated outcomes or other diseases,and articles with noADIPOQ +45T/G,+276G/T,or-11377C/Greported.

Finally,13 articles were adopted in this meta-analysis.Among them,nine studies investigated+45T/G,five studies investigated+276G/T,and five studies investigated-11377C/G.A flow diagram of study selection is presented in Figure 1.

Data were extracted by two researchers independently.We have extracted the following information from every included study:first author,year of publication,country,ethnicity,matching criteria,genotyping method,numbers of cases and controls,minor allele frequency in controls,and Hardy-Weinberg equilibrium (HWE)status.We obtained the HWE status of controls by calculating from genotype distributions using STATA12.0.The Newcastle-Ottawa quality assessment scale(NOS) was used for quality assessment of primary studies.The study would be regarded as a high-quality study when it had an NOS scores ≥ 6 (Table 1)[9].

Statistical analysis

We used Stata12.0 software for statistical analyses.A fixed-effects or random-effects model was used to merge OR and 95%CI based on allelic models,recessive models,and dominant models to evaluate the association between each genetic variant and the risk of GDM.TheZ-test was used for determining the significance of the merged OR.P＜ 0.05 was considered statistically significant.

We used the CochranQtest to assess the heterogeneity among the studies and HigginsI²statistic for quantifying the heterogeneity.We used the random-effects model as the merging method when the variant association presented significant interstudy heterogeneity (Qtest,P-value ＜ 0.05,orI²> 50%),otherwise,we used the fixed-effects model.Subgroup analysis was performed based on the ethnicity of the study population to evaluate ethnic-specific effects.Publication bias was tested by Begg’s funnel plot.

The statistical methods of this study were reviewed by Shi-Min Hu from Department of Epidemiology and Health Statistics,Xiangya School of Public Health,Central South University.

RESULTS

Main characteristics of all the included studies

Figure 1 Flow diagram of study selection.

Thirteen studies were adopted in this meta-analysis; among them,nine studies were about the association of+45T/Gand GDM[2,4,7,8,10-14],five studies were about the association of+276G/Tand GDM[4,10,12-14],and five studies were about the association of-11377C/Gand GDM[1,3,8,15,16].In total,1667 GDM cases and 1682 controls were included; ten studies were from Asian descendants[2,4,7,10-15],two studies were from European descendants[1,3],and one study was from a South American descendant[8].Detailed characteristics of all studies included are shown in Table 2.

Association between +45T/G and GDM

Nine articles evaluated the association of+45T/Gand GDM; eight of them were conducted in Asia and 1 in South America,with a total of 1024 GDM cases and 1059 controls.Heterogeneity test revealedP> 0.05 andI² ＜ 50%; the fixed-effects model was used.

The pooled results suggested a significant association between+45T/Gand GDM(allelic model:OR=1.45,95%CI:1.26-1.67; dominant model:OR=1.50,95%CI:1.25-1.79; recessive model:OR=2.00,95%CI:1.42-2.84).Ethnicity-based subgroup analysis showed that+45T/Gwas associated with GDM in Asians (allelic model:OR=1.47,95%CI:1.27-1.70; dominant model:OR=1.54,95%CI:1.27-1.85; recessive model:OR=2.00,95%CI:1.43-2.85).However,there was no association of+45T/Gwith the risk of GDM in South Americans (allelic model:OR=1.21,95%CI:0.68-2.41; dominant model:OR=1.13,95%CI:0.59-2.15; recessive model:OR=2.18,95%CI:0.43-11.07)(Figure 2 and Table 3).

Association between +276G/T and GDM

The association between+276G/Tand the risk of GDM was investigated by five studies,including 590 GDM cases and 595 controls.Heterogeneity test revealedP>0.05 andI² ＜ 50%,and the fixed-effects model was used.

The results showed that+276G/Twas not associated with the risk of GDM (allelic model:OR=0.88,95%CI:0.74-1.05; dominant model:OR=0.91,95%CI:0.65-1.26;recessive model:OR=0.82,95%CI:0.64-1.05) (Figure 3 and Table 3).

Association between -11377C/G and GDM

The association between-11377C/Gand the risk of GDM was investigated by five studies,of which two were conducted in Asia,two conducted in Europe,and one in South America,with a total of 722 GDM cases and 791 controls.Heterogeneity test revealedP＜ 0.05 andI² > 50%; the random-effects model was used.

The results showed that-11377C/Gwas not associated with the risk of GDM (allelic model:OR=0.96,95%CI:0.72-1.26; dominant model:OR=1.00,95%CI:0.73-1.37;recessive model:OR=0.90,95%CI:0.61-1.32).Ethnicity-based subgroup analysis also showed that-11377C/Gwas not associated with GDM in Asian (allelic model:OR=1.04,95%CI:0.77-1.41; dominant model:OR=1.09,95%CI:0.79-1.50; recessive model:OR=0.97,95%CI:0.51-1.86),European (allelic model:OR=0.94,95%CI:0.45-1.96;dominant model:OR=1.00,95%CI:0.42-2.33; recessive model:OR=0.87,95%CI:0.51-1.49) and South American populations (allelic model:OR=0.80,95%CI:0.50-1.29;dominant model:OR=0.77,95%CI:0.44-1.36; recessive model:OR=0.81,95%CI:0.51-1.86) (Figure 4 and Table 3).

Table 1 Quality assessment of included case control studies using the Newcastle-Ottawa scale

Publication bias

We used the Egger regression asymmetry test and Begg’s funnel plot to assess the public bias of the studies.The evidence of publication bias cannot be found in the meta-analysis of+45T/G(allelic model:continuity correctedP-value=1.000,Egger regression asymmetry testt=-0.62,P=0.554; recessive model:continuity correctedPvalue=0.466,Egger regression asymmetry testt=-0.15,P=0.883),+276G/T(allelic model:continuity correctedP-value=0.26,Egger regression asymmetry testt=-1.24,P=0.282) and-11377C/G(allelic model:continuity correctedP-value=0.221,Egger regression asymmetry testt=-2.48,P=0.089)(Figure 5).

Sensitivity analysis

To assess the stability of the results,we performed the sensitivity analysis by sequentially excluding individual studies for each meta-analysis.For the association between+45T/G,+276G/T,or-11377C/Gand GDM,there was no significant change of corresponding merged ORs when one study was sequentially excluded from every meta-analysis.Therefore,the results of our meta-analysis are stable and reliable(Figure 6).

DISCUSSION

ADIPOQ has been considered an important factor in regulating glucose and lipid metabolism.It is secreted by adipose tissue,which has a negative correlation with insulin resistance,T2D,and metabolic syndrome.ADIPOQ can increase insulin sensitivity,anti-inflammation,and anti-atherosclerosis,promote glucose uptake in muscle tissue,and inhibit intrahepatic synthetic glucose[17].

In the chromosomal region where theADIPOQgene is located,there aresusceptible sites of T2D and metabolic syndrome,and its SNPs can affect the level of ADIPOQ in blood,leading to obesity,insulin resistance,and the occurrence of T2D[15].

Table 2 Detailed characteristics of all eligible studies for the association between ADIPOQ single nucleotide polymorphism and gestational diabetes mellitus

Plasma ADIPOQ levels gradually decreased with gestational week during pregnancy,consistent with the gradual decrease of insulin sensitivity[3],and plasma adiponectin levels decreased more significantly in GDM women.This phenomenon is closely related to the decreased transcriptional activity ofADIPOQduring pregnancy.Previous studies revealed the association ofADIPOQSNPs,such as+45T/G[1,2,4,7,8,10-14],+276G/T[1,4,10,12,13,18],-11377C/G[1,3,8,15,16],-3971A/G[13,19],and-11426A/G[15,19],and the risk of GDM.A total of 66.7% (6 of 9) of the studies adopted in this meta-analysis reported that the+45T/Gpolymorphism increased the risk of GDM,and 40% (2 of 5) reported that the-11377C/Gpolymorphism was associated with GDM.A higher prevalence of the G allele was observed among women with GDM.All studies regarding+276G/Treported that this polymorphism had no association with the risk of GDM.

Thirteen studies were included in our study; nine studies were about+45T/G,with 1024 cases and 1059 controls,five studies were about+276G/T,with 590 cases and 595 controls,and five studies were about-11377C/G,with 722 cases and 791 controls.We not only had a larger sample size than previous studies but also performed subgroup analysis based on the ethnicity of the study population to evaluate ethnic-specific effects.+45T/Gwas proved by our meta-analysis to be a risk factor for GDM (allelic model:OR=1.45,95%CI:1.26-1.67,P=0.000),and 66.7% (6 of 9) of studies reported a positive result[2,4,7,10,12,14].Subgroup analysis showed that+45T/Gwas associated with GDM in Asians (allelic model:OR=1.47,95%CI:1.27-1.70,P=0.000) but not in South Americans.In addition,no association of+276G/Tor-11377C/Gand the risk of GDM was observed.

Obvious heterogeneity was detected among the-11377C/Gstudies (allelic model:I²=64.0%,P=0.03,dominant model:I²=55.0%,P=0.06,recessive model:I²=32.0%,P=0.21).We used subgroup analysis based on ethnicity,and the heterogeneity could not be reduced,indicating that a small sample size and other reasons may have influenced the heterogeneity.The association of-11377C/Gwith the risk of GDM remains to be verified by further studies.No heterogeneity was found in the studies of+45T/G(allelic model:I²=14.3%,P=0.32,dominant model:I²=18.1%,P=0.28,recessive model:I²=0.0%,P=0.90) and+276G/T(allelic model:I²=0.0%,P=0.74,dominant model:I²=0.0%,P=0.78,recessive model:I²=0.0%,P=0.83),so the conclusion that+45T/Ghas,but+276G/Thas no,association with the risk of GDM is relatively reliable.Begg’s funnel plot was used to test publication bias.The test showed that there was no publication bias among the studies.Sensitivity analysis indicated that the results are stable and reliable.

Figure 2 Forest plot for the association of ADIPOQ +45T/G polymorphism and gestational diabetes mellitus under the allelic model.

The SNP+45T/Gis a synonymous mutation (GGTGGG,Gly15Gly) at exon 2.The results of Yanget al[20]indicated that+45T/Gpolymorphism may influence the expression ofADIPOQby influencing RNA splicing and stability.Some studies reported that the G allele of+45T/Gpolymorphism in theADIPOQgene is associated with obesity,insulin resistance,and T2D in several populations.Very few studies have investigated the association ofADIPOQ +45T/Gpolymorphism with GDM and the results of these studies were controversial.

As far as we know,only+45T/Gwas involved in a previous meta-analysis that reported no association ofADIPOQ +45T/Gpolymorphism with the risk of GDM(allelic model:OR=1.17,95%CI:0.79-1.76; dominant model:OR=0.86,95%CI:0.50-1.48; recessive model:OR=1.21,95%CI:0.62-2.33)[21].The reason for this controversy is most likely the following:(1) the small sample size (case number=875,control number=884); (2) obvious heterogeneity (allP-values for heterogeneity less than 0.01); and (3) false HWE status of 25% studies (2 of 8) involved in that meta-analysis will cause insufficient power which may lead to the false-negative results.

The role of+276G/Tin the pathogenesis of metabolic syndrome and diabetes mellitus has also been reported to be contradictory.Commonly,the T allele has an association with a higher adiponectin level and protection against T2D[22,23],but some studies showed that T carriers have a higher risk of obesity and diabetes[24,25],or+276G/T polymorphism is not associated with T2D or GDM[26].

All the studies on+276G/Tpolymorphism included in our study came from the Chinese population.This result suggests that+276G/Tin the Chinese population may associate with the risk of GDM.The studies about the association between+276G/Tand GDM from other countries could not be found.Geographical,environmental,and genetic factors of different ethnic groups lead to different susceptibility to diabetes;therefore,we need more studies about the association between+276G/Tand GDM of other ethnic groups to reach reliable conclusions.

Zhanget al[27]found thatADIPOQgene promoter region has four transcription stimulatory protein (SP1) binding sites,while the G allele of- 11377 C/Gin the promoter region can change the DNA sequence of one of the SP1 binding sites,leading to the loss of binding force to SP1.This may reduce theADIPOQgene transcription activity,inhibit the expression of genes,and lead to lower plasma ADIPOQ,which could associate with glucolipid metabolic abnormalities and insulin resistance.Consistent with the results of Vasseuret al[28],Petroneet al[29]reported that-11377Ghaplotype is associated with low plasma ADIPOQ levels and T2D.However,due to ethnic and geographical differences,the results of studies on the association betweenADIPOQgene-11377C/Gpolymorphism and diabetes mellitus are not completely consistent[30].

Table 3 Main results of the pooled odds ratios in meta-analysis for the association between ADIPOQ polymorphisms and gestational diabetes mellitus

According to the literatures on-11377C/Gpolymorphism included in our studies,Asians accounted for 40%,Europeans for 40%,and South Americans for 20%.Due to the small sample size and large heterogeneity of each ethnic subgroup,our results,which are inconsistent with the previous studies,are unreliable.

The limitations of this study should be considered.First,the number of cases and controls involved in the meta-analysis for exploring the association ofADIPOQand GDM in different ethnicities may have little power,and studies with larger sample sizes and multiple ethnicities are needed.Second,GDM has complicated cases,with genetic susceptibility,environmental triggers,and acquired dispositions,such as age,gestational weeks,condition of nutrition,and physique.In this meta-analysis,we failed to conduct a multivariate analysis of confounders.Therefore,further comprehensive studies with strict matching criteria for cases and controls are needed.Third,few studies have reported the association between polymorphisms and serum ADIPOQ levels,so genotype-phenotype analysis was prevented[21].

In conclusion,our meta-analysis reveals the association of theADIPOQ +45T/Gpolymorphism and the risk of GDM; this polymorphism increases GDM risk in Asian populations.Another two polymorphisms,+276G/Tand-11377C/G,seem to have no association with the risk of GDM.Prospective studies of high quality with larger sample sizes are required to reveal the association ofADIPOQpolymorphisms with GDM,the existence of ethnicity-specific factors,and the role that ADIPOQ polymorphisms play in pathology.

Figure 3 Forest plot for the association of ADIPOQ +276G/T polymorphism and gestational diabetes mellitus under the allelic model.

Figure 4 Forest plot for the association of ADIPOQ -11377C/G polymorphism and gestational diabetes mellitus under the allelic model.

Figure 5 Begg’s funnel plots for testing publication bias.

Figure 6 Sensitivity analysis between ADPIOQ polymorphisms and gestational diabetes mellitus in all studies.

ARTICLE HIGHLIGHTS

Research background

Many biochemical mediators that are synthesized in the adipose tissue and secreted in the circulation,such as leptin,adiponectin (ADIPOQ),and resistin,are thought to be involved in obesity and insulin resistance.ADIPOQ is produced in the adipose tissue and regulates a variety of metabolic processes such as lipid metabolism,glucose and fatty acid oxidation.This hormone can reduce insulin resistance,improve lipid metabolism,and exert anti-inflammatory effects.Plasma ADIPOQ levels are decreased in patients with type 2 diabetes,metabolic syndrome,and obesity.Previous studies have shown thatADIPOQsingle nucleotide polymorphisms can affect plasma ADIPOQ concentrations,which in turn affect insulin sensitivity.

Research motivation

Previous studies have evaluated the relationship betweenADIPOQpolymorphisms and gestational diabetes mellitus (GDM),but the results of the association betweenADIPOQpolymorphisms and GDM is uncertain.

Research objectives

We evaluated the association betweenADIPOQ +45T/G,+276G/T,and-11377C/Gpolymorphisms and GDM with a bigger sample size and less heterogeneity.Moreover,subgroup analysis was performed by ethnicity.

Research methods

Potentially related articles on human fat metabolism and GDM gene research published before October 20,2018 were retrieved through the electronic databases EMBASE,Web of Science,PubMed,WANFANG DATA,and China National Knowledge Infrastructure.A fixed-effects or random-effects model was used to calculate pooled odds ratios (ORs) with 95% confidence intervals (CIs),based on the between-study heterogeneity to evaluate the association betweenAIDPOQ +45T/G,+276G/T,and-11377C/Gpolymorphisms and the risk of GDM.

Research results

Nine+45T/Gstudies included 1024 GDM cases and 1059 controls,five+276G/Tstudies included 590 GDM cases and 595 controls,and five-11377C/Gstudies included 722 GDM cases and 791 controls.Pooled ORs showed that+45T/Gincreased Asian GDM risk (allele model OR=1.47,95%CI:1.27-1.70,P=0.000; dominant model OR=1.54,95%CI:1.27-1.85,P=0.000;recessive mode:OR=2.00,95%CI:1.43-2.85,P=0.000),but not in South Americans (equal pattern:OR=1.21,95%CI:0.68-2.41,P=0.510; dominant model OR=1.13,95%CI:0.59-2.15,P=0.710; recessive mode OR=2.18,95%CI:0.43-11.07,P=0.350).There was no significant correlation between+276G/T(allele model OR=0.88,95%CI:0.74-1.05,P=0.158; dominant model OR=0.91,95%CI:0.65-1.26,P=0.561; recessive mode:OR=0.82,95%CI:0.64-1.05,P=0.118) or-11377C/G(equal pattern:OR=0.96,95%CI:0.72-1.26,P=0.750; dominant model OR=1.00,95%CI:0.73-1.37,P=0.980; recessive model:OR=0.90,95%CI:0.61-1.32,P=0.570) and GDM risk.

Research conclusions

Our meta-analysis reveals the association of theADIPOQ +45T/Gpolymorphism and the risk of GDM; this polymorphism increases GDM risk in Asian populations.

Research perspectives

In order to reveal the association ofADIPOQpolymorphisms with GDM,the existence of ethnicity-specific factors,and the role thatADIPOQpolymorphisms play in pathology,studies focused on delineating ethnicity-specific factors with larger sample sizes are needed.