this study was approved by the institutional ethics committee review boards of chang gung memorial hospital linkou medical center and chang gung university .
we recruited 123 patients , and all patients gave written informed consent to participate in the study .
the screening glucose challenge test for gestational diabetes was performed as previously described ( 24 ) .
genomic dna samples of participants were extracted and purified from anticoagulated blood with the dneasy blood & tissue kit ( qiagen , venlo , netherlands ) .
genotyping of snps was performed using the taqman validated snp genotyping assays ( applied biosystems , foster city , ca ) .
blood samples were collected from patients for hormone measurement 1 h after administration of the oral glucose tolerance test .
total gip , insulin , and c - peptide levels in human serum were measured using sandwich elisa kits ( millipore , billerica , ma ; mercodia , uppsala , sweden ; and calbiotech , spring valley , ca ) .
a 2.15-kb fragment of human gip promoter with the a allele at rs3895874 , rs3809770 , rs3848460 , and rs937301 ( 2073 bp to + 77 bp ) was chemically synthesized ( genescript inc . ,
piscataway , nj ) and subcloned into the pgl4.2 luciferase reporter vector ( promega corp .
promoter fragments with ancestral haplotypes were obtained using the site - directed mutagenesis . in the promoter reporter study , each experiment was conducted at least three times with three or four replicates for each treatment .
patients glucose challenge responses and serum hormone profiles were compared with the  test or the student t test with welch correction .
ratios of patients with glucose levels exceeding the threshold were analyzed with the  test .
all data were presented as mean  sem , and the statistical significance cutoff value was 0.05 .
this study was approved by the institutional ethics committee review boards of chang gung memorial hospital linkou medical center and chang gung university .
we recruited 123 patients , and all patients gave written informed consent to participate in the study .
the screening glucose challenge test for gestational diabetes was performed as previously described ( 24 ) .
genomic dna samples of participants were extracted and purified from anticoagulated blood with the dneasy blood & tissue kit ( qiagen , venlo , netherlands ) .
genotyping of snps was performed using the taqman validated snp genotyping assays ( applied biosystems , foster city , ca ) .
blood samples were collected from patients for hormone measurement 1 h after administration of the oral glucose tolerance test .
total gip , insulin , and c - peptide levels in human serum were measured using sandwich elisa kits ( millipore , billerica , ma ; mercodia , uppsala , sweden ; and calbiotech , spring valley , ca ) .
a 2.15-kb fragment of human gip promoter with the a allele at rs3895874 , rs3809770 , rs3848460 , and rs937301 ( 2073 bp to + 77 bp ) was chemically synthesized ( genescript inc . ,
piscataway , nj ) and subcloned into the pgl4.2 luciferase reporter vector ( promega corp . , madison , wi ) .
promoter fragments with ancestral haplotypes were obtained using the site - directed mutagenesis . in the promoter reporter study , each experiment was conducted at least three times with three or four replicates for each treatment .
patients glucose challenge responses and serum hormone profiles were compared with the  test or the student t test with welch correction . ratios of patients with glucose levels exceeding the threshold were analyzed with the  test .
all data were presented as mean  sem , and the statistical significance cutoff value was 0.05 .
to systematically identify putative metabolic modifiers , we curated and screened snps in 207 gene loci that have previously been implicated in the regulation of diabetes - related and/or obesity - related traits using fst ( supplementary table 1 ) .
the empiric distribution of the fst statistic has been used to detect genomic regions that have rapidly increased in frequency as a result of local selective pressures ( 27 ) . of these 207 genes ,
59 carried genic variants with fst values in the top 5% bracket in comparisons between the corresponding hapmap ii populations yri ( yoruba from ibadan ) , ceu ( u.s .
residents with northern and western european ancestry ) , and asn ( pooled samples of chinese from beijing [ chb ] and japanese from tokyo [ jpt ] ; supplementary table 1 ) .
the genic region in 29 of these 59 genes also contained another indication of local selection : long haplotypes in one of the hapmap populations .
on the basis of the presence of highly divergent allele frequencies between eurasians and africans , the presence of ld and extended haplotypes in eurasians , and a > 30% minor allele frequency in the overall eurasian population , we identified seven snps in the 5 gene region of cdkal1 , cyb5r4 , gad2 , gip , and pparg as potential common metabolic modifiers ( table 1 ) . in earlier gwa studies ,
select variants in cdkal1 , cyb5r4 , gad2 , and pparg were associated with type 2 diabetes - related or obesity - related traits ( 2831 ) . by contrast , there has been no report of linkage of gip variants in gwa studies .
because the three gip variants ( rs3895874 , rs3848460 , and rs937301 ) are highly linked compared with those that appeared alone in other candidate genes ( 32 ) , and because these gip variants are partially linked with a nonsynonymous gip snp ( rs2291725 ) ( 22 ) , these gip variants have a low likelihood of being false positive ( table 1 ) .
snps from the 5 gene region that exhibited signals of partial selection asn , pooled samples of chinese from beijing and japanese from tokyo ; ceu , utah residents with northern and western european ancestry ; gd , gestational diabetes ; ld , linkage - disequilibrium ; ob , obesity ; snp , single nucleotide polymorphism ; t2d , type 2 diabetes ; yri , yoruba in ibadan , nigeria . *
the allele frequency is significantly different among hapmap ii populations ( p < 0.01 ) .
fine mapping of the gip locus showed that these variants are linked with more than three dozen neighboring snps ( from rs9904761 to rs3895874 on chr17 : 44,31144,402 kb ) , and these linked snps exhibited fst values in the top 210% bracket in comparisons between yri and asn ( supplementary fig .
a plot of a 250-kb region of genotypes around gip in the three hapmap ii populations showed that whereas genotypes surrounding gip in yri exhibited a high degree of homozygosity for ancestral alleles ( fig .
1a , upper panel ) , in the same region , genotypes of asn and ceu exhibited a high degree of homozygosity for the derived alleles ( fig .
analysis of genic snps at the gip locus between 11 populations from the hapmap iii project ( 34 ) showed that fst values were the highest for comparisons between east asian and african populations ( fig .
this result is consistent with our recent finding that a nonsynonymous variant in the exon iv of gip ( rs2291725 ) and a 71-kb haplotype block surrounding this variant were positively selected in eurasians in the last 2,000 to 11,800 years ( 22 ) .
differential distribution of polymorphic alleles at the gip locus in human populations . a : visual depiction of genotypes within a 250-kb region around gip in yri , ceu , and asn populations .
polymorphic sites are color coded according to their allelic state . individuals with a homozygous genotype with ancestral allele ( aa ) and a homozygous genotype with derived allele ( aa )
the location of snps in the gip gene region is indicated by a green rectangle .
the position of snps within the genomic region is indicated by blue vertical lines in the bottom panel .
b : average fst for genic snps at the gip locus between 11 populations of the hapmap phase iii dataset .
population descriptors : asw , african ancestry in southwest u.s . ; ceu , utah residents with northern and western european ancestry from the centre de'etude du polymorphism humain ( ceph ) collection ; chb , han - chinese in beijing , china ; chd , chinese in metropolitan denver , co ; gih , gujarati indians in houston , tx ; jpt , japanese in tokyo , japan ; lwk , luhya in webuye , kenya ; mex , mexican ancestry in los angeles , ca ; mkk , maasai in kinyawa , kenya ; tsi , toscans in italy ; yri , yoruba in ibadan , nigeria .
the highest population differentiation was found between asian ( han - chinese and japanese ) and african ( luhya , maasai , or yoruba ) populations .
( a high - quality color representation of this figure is available in the online issue . )
interestingly , analyses of ld and haplotype block diversity of the genomic region that encompassed the three gip variants at the 5 gene region showed that a 91-kb ld block was represented by five inferred haplotypes in ceu and four in and asn chromosomes , respectively ( fig . 2 ; supplementary fig .
in addition , these analyses showed that the high fst values and the major shift in allele frequencies of gip variants at the 5 gene region between eurasian and african populations could be attributed to the increase of a derived haplotype ( haplotype 50 in fig .
because the three 5 gene region variants are completely linked with the positively selected rs2291725 in asn ( 22 ) , we inferred that rs3895874 , rs3848460 , and rs937301 at the human gip promoter region were positively selected in east asian populations as well .
frequencies of the 58 inferred haplotypes found in the 91-kb ld region ( from rs9904761 to rs3895874 ) surrounding gip in yri , ceu , and asn .
most of the ceu and asn chromosomes contained derived haplotypes ( with aaaa residues at positions 7679 ) , whereas yri chromosomes exhibited significant divergence in the haplotype composition .
population descriptors : asn , pooled samples of chinese from beijing and japanese from tokyo ; ceu , utah residents with northern and western european ancestry from the centre de'etude du polymorphism humain ( ceph ) collection ; yri , yoruba in ibadan , nigeria .
( a high - quality color representation of this figure is available in the online issue . )
importantly , we found that the 5 gene region variants ( i.e. , rs3895874 , rs3809770 , rs3848460 , and rs937301 at positions 1920 , 1650 , 1158 , and 320 of gip , respectively ) are located in a haplotype block that is separated from the one containing the nonsynonymous variant rs2291725 by a hotspot for recombination in ceu and yri ( supplementary fig .
a tetranucleotide polymorphism was represented by three inferred haplotypes in all three populations ( referred to as derived gip , ancestral gip , and ancestral gip haplotypes in the following text ; table 2 ) .
the derived gip haplotype , which was found in 18.3% of yri chromosomes , has become the dominant haplotype and has a frequency exceeding 50 and 75% in ceu and asn , respectively ( table 2 ) .
by contrast , the dominant ancestral haplotype in yri chromosomes ( gip , 50.8% ) was only found in 1.1% of asn chromosomes . therefore
, these 5 gene region variants and the nonsynonymous variant rs2291725 could be selected differentially in the overall hapmap populations and represent causal mutations independent of the nonsynonymous variant rs2291725 . frequency of inferred gip haplotypes and rs3895874 genotypes at the 5 gene region in hapmap ii populations asn , pooled samples of chinese from beijing and japanese from tokyo ; ceu , utah residents with northern and western european ancestry from the centre de'etude du polymorphism humain ( ceph ) collection ; yri , yoruba in ibadan , nigeria .
* frequencies of gip haplotypes were significantly different among hapmap ii populations ( p < 0.0001 ) .
because only functional investigations can convincingly demonstrate causal mutations , we sought to obtain direct evidence that these 5 gene region snps represent causal mutations for the population genetics observation .
we constructed and tested gip promoter reporters with each of the three major haplotypes in transfected human embryonic kidney ( hek ) 293 t cells ( fig .
measurement of promoter reporter activities showed that constructs with an ancestral haplotype ( gip or gip ) exhibited luciferase reporter activity 2545% higher than that of a derived haplotype ( gip , p < 0.01 ; fig .
because the gip promoter region contains elements that are important for regulation by transcriptional factors , including pax6 and gata4 ( 23 ) , we also determined whether the gip promoter reporter activity was haplotype - dependent in the presence of these transcription factors . as expected ,
coexpression of pax6 or gata4 increased the basal reporter activities by 1.52.5-fold and 0.70.8-fold , respectively ( fig .
importantly , we found that reporters with an ancestral haplotype consistently exhibited significantly higher activities than those containing the derived haplotype in the presence of pax6 or gata4 ( p < 0.01 ) .
together , these results suggested that derived alleles at rs3895874 , rs3848460 , and rs937301but not rs3809770represent functional mutations .
the gip promoter reporter activity is haplotype - dependent . a : schematic representation of luciferase reporters containing one of the three inferred haplotypes ( derived : gip ; ancestral : gip and gip ) found in a 2.15-kb fragment of the gip promoter .
the genomic fragment from position 2073 to + 77 bp of gip contains four snps ( from rs3895874 to rs937301 ) , and three of them are linked .
genomic fragments that contained each of the three haplotypes were subcloned in the pgl4.2 luciferase reporter .
b : the luciferase reporter activity is shown in hek293 t cells transfected with combinations of reporter and expression vectors for pax6 and gata4 .
the reporter activity is haplotype - dependent in the absence or presence of select transcription factors . to compare gip promoter reporter activities , hek293 t cells in 6-well
culture dishes were transfected with different combinations of transcription factor expression vector and a select gip promoter reporter using lipofectamine 2000 ( life technology inc . ,
, equal amounts of pcmv and pgl4.2 expression vectors and a one - tenth aliquot of a -galactosidase expression vector were transfected . at 48 h after transfection
, cells were free - thawed once in lysis buffer , and luciferase activity in supernatant was assayed using the steady - glo luciferase assay system ( promega corp . ) and a lumimark microplate reader ( bio - rad corp .
the reporter activity has been expressed as the ratio of relative luciferase unit/-galactosidase activity in transfected cells .
( a high - quality color representation of this figure is available in the online issue . ) given that earlier genome studies have not reported an association between gip variants and any trait , we speculated that the partial selection of gip variants in eurasian populations could be associated with adaptation at a life stage that is vulnerable to environmental changes and that has not been specifically studied . because pregnancy represents a critical life stage that subjects individuals to excessive metabolic load , and
because its success has a major impact on reproductive fitness , we hypothesized that studies of phenotypic variation during pregnancy may provide a sensitive model to investigate the role of gip haplotypes . to test for an allelic effect of gip variants , we studied east asian patients for proof - of - concept testing because selection pressures are most likely ongoing in populations that exhibit the most significant evidence .
we genotyped a panel of 123 unrelated han - chinese women who underwent a screening glucose challenge test for gestational diabetes during the 23rd to the 29th weeks of pregnancy , and these patients were assigned to three genotype groups ( gip , gip , or gip ) based on alleles at rs3895874 , rs3848460 , and rs937301 .
the frequency of rs3895874 in these patients was similar to that of the asn population and was in the hardy  weinberg equilibrium ; frequencies of gip , gip , and gip genotypes were 0.089 , 0.480 , and 0.431 , respectively .
in addition , alleles at rs3895874 in these patients were in absolute ld with those at rs3848460 and rs937301 . measurements of serum gip and glucose levels showed that circulating gip and glucose at 1 h after the challenge test were 20.6219.9 pg / ml , and 72230 mg / dl , respectively .
consistent with in vitro promoter reporter assays , gip levels in patients carrying the ancestral gip haplotype ( gip heterozygote and gip homozygote ) were significantly higher than those with a homozygous gip genotype ( fig .
by contrast , serum levels of glucose , insulin , and c - peptide , as well as age and bmi , were not significantly different among patients ( table 3 ) .
gip haplotypes impart a difference in gip response and glucose metabolism after glucose challenge tests .
a : measurements of serum gip and glucose levels in 123 patients during the 23rd to the 29th weeks of pregnancy at 1 h after the 50-g glucose challenge test .
the gip level was significantly higher in patients carrying a gip haplotype compared with those homozygous for the gip haplotype .
b : in 72 patients with the dominant gipr genotype , both serum gip and glucose levels were both significantly different between gip homozygotes and those carrying a gip haplotype .
association between gip haplotypes and gip response after glucose challenge tests in 123 pregnant han - chinese women gip , glucose - dependent insulinotropic polypeptide ; or , odds ratio .
significantly different from patients with a gip haplotype (  test , p < 0.05 ) . because two linked gip receptor ( gipr ) variants rs10423928 and rs1800437 ( referred to as the gipr mutation in table 3 ) have recently been shown to be associated with glucose and insulin levels after oral glucose challenge tests as well as the incretin effect in nondiabetic individuals in gwa studies ( 35,36 ) , we also genotyped these variants and sought to isolate the potential confounding effect of gipr variants .
we found no association between gipr snps and serum glucose or hormone levels , but levels of serum glucose , in addition to gip , were significantly different between gip homozygotes , and heterozygotes and the gip homozygotes combined within the pool of patients with the dominant gipr genotype ( fig .
moreover , we noticed that after controlling for the variation at gipr , the derived gip genotype was associated with increased risk of having a glucose level that exceeded the 140 mg / dl threshold ( odds ratio 3.53 [ 95% ci 1.259.92 ] ; p = 0.015 ; table 3 ) . among patients with a gipr genotype , 48.3% of patients with gip homozygotes exhibited glucose levels that exceeded the threshold , whereas only 20.9% of the remaining patients did
therefore , the homozygous gip genotype could be associated with a reduced gip response and reduced capability of maintaining glucose homeostasis .
to systematically identify putative metabolic modifiers , we curated and screened snps in 207 gene loci that have previously been implicated in the regulation of diabetes - related and/or obesity - related traits using fst ( supplementary table 1 ) .
the empiric distribution of the fst statistic has been used to detect genomic regions that have rapidly increased in frequency as a result of local selective pressures ( 27 ) . of these 207 genes ,
59 carried genic variants with fst values in the top 5% bracket in comparisons between the corresponding hapmap ii populations yri ( yoruba from ibadan ) , ceu ( u.s .
residents with northern and western european ancestry ) , and asn ( pooled samples of chinese from beijing [ chb ] and japanese from tokyo [ jpt ] ; supplementary table 1 ) .
the genic region in 29 of these 59 genes also contained another indication of local selection : long haplotypes in one of the hapmap populations .
on the basis of the presence of highly divergent allele frequencies between eurasians and africans , the presence of ld and extended haplotypes in eurasians , and a > 30% minor allele frequency in the overall eurasian population , we identified seven snps in the 5 gene region of cdkal1 , cyb5r4 , gad2 , gip , and pparg as potential common metabolic modifiers ( table 1 ) . in earlier gwa studies ,
select variants in cdkal1 , cyb5r4 , gad2 , and pparg were associated with type 2 diabetes - related or obesity - related traits ( 2831 ) . by contrast , there has been no report of linkage of gip variants in gwa studies .
because the three gip variants ( rs3895874 , rs3848460 , and rs937301 ) are highly linked compared with those that appeared alone in other candidate genes ( 32 ) , and because these gip variants are partially linked with a nonsynonymous gip snp ( rs2291725 ) ( 22 ) , these gip variants have a low likelihood of being false positive ( table 1 ) .
snps from the 5 gene region that exhibited signals of partial selection asn , pooled samples of chinese from beijing and japanese from tokyo ; ceu , utah residents with northern and western european ancestry ; gd , gestational diabetes ; ld , linkage - disequilibrium ; ob , obesity ; snp , single nucleotide polymorphism ; t2d , type 2 diabetes ; yri , yoruba in ibadan , nigeria . *
the allele frequency is significantly different among hapmap ii populations ( p < 0.01 ) .
fine mapping of the gip locus showed that these variants are linked with more than three dozen neighboring snps ( from rs9904761 to rs3895874 on chr17 : 44,31144,402 kb ) , and these linked snps exhibited fst values in the top 210% bracket in comparisons between yri and asn ( supplementary fig .
a plot of a 250-kb region of genotypes around gip in the three hapmap ii populations showed that whereas genotypes surrounding gip in yri exhibited a high degree of homozygosity for ancestral alleles ( fig .
1a , upper panel ) , in the same region , genotypes of asn and ceu exhibited a high degree of homozygosity for the derived alleles ( fig .
analysis of genic snps at the gip locus between 11 populations from the hapmap iii project ( 34 ) showed that fst values were the highest for comparisons between east asian and african populations ( fig .
this result is consistent with our recent finding that a nonsynonymous variant in the exon iv of gip ( rs2291725 ) and a 71-kb haplotype block surrounding this variant were positively selected in eurasians in the last 2,000 to 11,800 years ( 22 ) .
differential distribution of polymorphic alleles at the gip locus in human populations . a : visual depiction of genotypes within a 250-kb region around gip in yri , ceu , and asn populations .
polymorphic sites are color coded according to their allelic state . individuals with a homozygous genotype with ancestral allele ( aa ) and a homozygous genotype with derived allele ( aa )
the location of snps in the gip gene region is indicated by a green rectangle .
the position of snps within the genomic region is indicated by blue vertical lines in the bottom panel .
b : average fst for genic snps at the gip locus between 11 populations of the hapmap phase iii dataset .
population descriptors : asw , african ancestry in southwest u.s . ; ceu , utah residents with northern and western european ancestry from the centre de'etude du polymorphism humain ( ceph ) collection ; chb , han - chinese in beijing , china ; chd , chinese in metropolitan denver , co ; gih , gujarati indians in houston , tx ; jpt , japanese in tokyo , japan ; lwk , luhya in webuye , kenya ; mex , mexican ancestry in los angeles , ca ; mkk , maasai in kinyawa , kenya ; tsi , toscans in italy ; yri , yoruba in ibadan , nigeria .
the highest population differentiation was found between asian ( han - chinese and japanese ) and african ( luhya , maasai , or yoruba ) populations .
( a high - quality color representation of this figure is available in the online issue . )
interestingly , analyses of ld and haplotype block diversity of the genomic region that encompassed the three gip variants at the 5 gene region showed that a 91-kb ld block was represented by five inferred haplotypes in ceu and four in and asn chromosomes , respectively ( fig . 2 ; supplementary fig .
in addition , these analyses showed that the high fst values and the major shift in allele frequencies of gip variants at the 5 gene region between eurasian and african populations could be attributed to the increase of a derived haplotype ( haplotype 50 in fig .
because the three 5 gene region variants are completely linked with the positively selected rs2291725 in asn ( 22 ) , we inferred that rs3895874 , rs3848460 , and rs937301 at the human gip promoter region were positively selected in east asian populations as well .
frequencies of the 58 inferred haplotypes found in the 91-kb ld region ( from rs9904761 to rs3895874 ) surrounding gip in yri , ceu , and asn .
most of the ceu and asn chromosomes contained derived haplotypes ( with aaaa residues at positions 7679 ) , whereas yri chromosomes exhibited significant divergence in the haplotype composition .
population descriptors : asn , pooled samples of chinese from beijing and japanese from tokyo ; ceu , utah residents with northern and western european ancestry from the centre de'etude du polymorphism humain ( ceph ) collection ; yri , yoruba in ibadan , nigeria .
( a high - quality color representation of this figure is available in the online issue . )
importantly , we found that the 5 gene region variants ( i.e. , rs3895874 , rs3809770 , rs3848460 , and rs937301 at positions 1920 , 1650 , 1158 , and 320 of gip , respectively ) are located in a haplotype block that is separated from the one containing the nonsynonymous variant rs2291725 by a hotspot for recombination in ceu and yri ( supplementary fig .
a tetranucleotide polymorphism was represented by three inferred haplotypes in all three populations ( referred to as derived gip , ancestral gip , and ancestral gip haplotypes in the following text ; table 2 ) .
the derived gip haplotype , which was found in 18.3% of yri chromosomes , has become the dominant haplotype and has a frequency exceeding 50 and 75% in ceu and asn , respectively ( table 2 ) .
by contrast , the dominant ancestral haplotype in yri chromosomes ( gip , 50.8% ) was only found in 1.1% of asn chromosomes . therefore
, these 5 gene region variants and the nonsynonymous variant rs2291725 could be selected differentially in the overall hapmap populations and represent causal mutations independent of the nonsynonymous variant rs2291725 . frequency of inferred gip haplotypes and rs3895874 genotypes at the 5 gene region in hapmap ii populations asn , pooled samples of chinese from beijing and japanese from tokyo ; ceu , utah residents with northern and western european ancestry from the centre de'etude du polymorphism humain ( ceph ) collection ; yri , yoruba in ibadan , nigeria .
* frequencies of gip haplotypes were significantly different among hapmap ii populations ( p < 0.0001 ) .
because only functional investigations can convincingly demonstrate causal mutations , we sought to obtain direct evidence that these 5 gene region snps represent causal mutations for the population genetics observation .
we constructed and tested gip promoter reporters with each of the three major haplotypes in transfected human embryonic kidney ( hek ) 293 t cells ( fig .
measurement of promoter reporter activities showed that constructs with an ancestral haplotype ( gip or gip ) exhibited luciferase reporter activity 2545% higher than that of a derived haplotype ( gip , p < 0.01 ; fig .
because the gip promoter region contains elements that are important for regulation by transcriptional factors , including pax6 and gata4 ( 23 ) , we also determined whether the gip promoter reporter activity was haplotype - dependent in the presence of these transcription factors . as expected ,
coexpression of pax6 or gata4 increased the basal reporter activities by 1.52.5-fold and 0.70.8-fold , respectively ( fig .
importantly , we found that reporters with an ancestral haplotype consistently exhibited significantly higher activities than those containing the derived haplotype in the presence of pax6 or gata4 ( p < 0.01 ) .
together , these results suggested that derived alleles at rs3895874 , rs3848460 , and rs937301but not rs3809770represent functional mutations .
. a : schematic representation of luciferase reporters containing one of the three inferred haplotypes ( derived : gip ; ancestral : gip and gip ) found in a 2.15-kb fragment of the gip promoter .
the genomic fragment from position 2073 to + 77 bp of gip contains four snps ( from rs3895874 to rs937301 ) , and three of them are linked .
genomic fragments that contained each of the three haplotypes were subcloned in the pgl4.2 luciferase reporter .
b : the luciferase reporter activity is shown in hek293 t cells transfected with combinations of reporter and expression vectors for pax6 and gata4 .
the reporter activity is haplotype - dependent in the absence or presence of select transcription factors . to compare gip promoter reporter activities , hek293 t cells in 6-well
culture dishes were transfected with different combinations of transcription factor expression vector and a select gip promoter reporter using lipofectamine 2000 ( life technology inc . , carlsbad , ca ) .
in each well , equal amounts of pcmv and pgl4.2 expression vectors and a one - tenth aliquot of a -galactosidase expression vector were transfected . at 48 h after transfection
, cells were free - thawed once in lysis buffer , and luciferase activity in supernatant was assayed using the steady - glo luciferase assay system ( promega corp . ) and a lumimark microplate reader ( bio - rad corp . , hercules , ca ) .
the reporter activity has been expressed as the ratio of relative luciferase unit/-galactosidase activity in transfected cells .
( a high - quality color representation of this figure is available in the online issue . )
given that earlier genome studies have not reported an association between gip variants and any trait , we speculated that the partial selection of gip variants in eurasian populations could be associated with adaptation at a life stage that is vulnerable to environmental changes and that has not been specifically studied . because pregnancy represents a critical life stage that subjects individuals to excessive metabolic load , and because its success has a major impact on reproductive fitness , we hypothesized that studies of phenotypic variation during pregnancy may provide a sensitive model to investigate the role of gip haplotypes . to test for an allelic effect of gip variants , we studied east asian patients for proof - of - concept testing because selection pressures are most likely ongoing in populations that exhibit the most significant evidence .
we genotyped a panel of 123 unrelated han - chinese women who underwent a screening glucose challenge test for gestational diabetes during the 23rd to the 29th weeks of pregnancy , and these patients were assigned to three genotype groups ( gip , gip , or gip ) based on alleles at rs3895874 , rs3848460 , and rs937301 .
the frequency of rs3895874 in these patients was similar to that of the asn population and was in the hardy  weinberg equilibrium ; frequencies of gip , gip , and gip genotypes were 0.089 , 0.480 , and 0.431 , respectively .
in addition , alleles at rs3895874 in these patients were in absolute ld with those at rs3848460 and rs937301 .
measurements of serum gip and glucose levels showed that circulating gip and glucose at 1 h after the challenge test were 20.6219.9 pg / ml , and 72230 mg / dl , respectively . consistent with in vitro promoter reporter assays , gip levels in patients carrying the ancestral gip haplotype ( gip heterozygote and gip homozygote ) were significantly higher than those with a homozygous gip genotype ( fig .
by contrast , serum levels of glucose , insulin , and c - peptide , as well as age and bmi , were not significantly different among patients ( table 3 ) .
gip haplotypes impart a difference in gip response and glucose metabolism after glucose challenge tests . a : measurements of serum gip and glucose levels in 123 patients during the 23rd to the 29th weeks of pregnancy at 1 h after the 50-g glucose challenge test .
the gip level was significantly higher in patients carrying a gip haplotype compared with those homozygous for the gip haplotype .
b : in 72 patients with the dominant gipr genotype , both serum gip and glucose levels were both significantly different between gip homozygotes and those carrying a gip haplotype .
association between gip haplotypes and gip response after glucose challenge tests in 123 pregnant han - chinese women gip , glucose - dependent insulinotropic polypeptide ; or , odds ratio . *
significantly different from patients with a gip haplotype (  test , p < 0.05 ) . because two linked gip receptor ( gipr ) variants rs10423928 and rs1800437 ( referred to as the gipr mutation in table 3 ) have recently been shown to be associated with glucose and insulin levels after oral glucose challenge tests as well as the incretin effect in nondiabetic individuals in gwa studies ( 35,36 ) , we also genotyped these variants and sought to isolate the potential confounding effect of gipr variants .
we found no association between gipr snps and serum glucose or hormone levels , but levels of serum glucose , in addition to gip , were significantly different between gip homozygotes , and heterozygotes and the gip homozygotes combined within the pool of patients with the dominant gipr genotype ( fig .
moreover , we noticed that after controlling for the variation at gipr , the derived gip genotype was associated with increased risk of having a glucose level that exceeded the 140 mg / dl threshold ( odds ratio 3.53 [ 95% ci 1.259.92 ] ; p = 0.015 ; table 3 ) . among patients with a gipr genotype , 48.3% of patients with gip homozygotes exhibited glucose levels that exceeded the threshold , whereas only 20.9% of the remaining patients did .
therefore , the homozygous gip genotype could be associated with a reduced gip response and reduced capability of maintaining glucose homeostasis .
on the basis of studies of signatures of selection , in vitro promoter assays , and glucose challenge tests in humans , we show that it is possible to identify causal variants related to energy - balance regulation by focusing on genic snps that were subject to environmental selection in a subset of candidate genes .
specifically , we demonstrated that gip variants at the 5 gene region represent metabolic modifiers that contribute to phenotypic variation in gip response and glucose metabolism .
further characterization of these causal variants would open a new venue for understanding the molecular mechanisms that underlie phenotypic variations in energy - balance regulation and improve our ability to stratify and interpret clinical outcomes associated with the gip signaling pathway . for decades
, adaptive selection was assumed to be rare ; however , recent studies have demonstrated that adaptive substitution is pervasive in human genomes ( 1,7,37 ) . despite this progress , it is obvious that population differentiation characteristics of human genomes have yet to be fully explored because physiologic consequences of almost all of these past gene 
environmental interactions remain to be verified experimentally ( 1,2 ) . on the other hand , because many selection pressures could be heterogeneous or reversible in a short time , the signature of selection may have eroded in genes that were responsive to cultural selection pressures ( 12,38 ) compared with those shielded from heterogeneous selection ( e.g. , the adaptation to environmental oxygen levels and altitude ) ( 5,39 ) .
we therefore reasoned that important metabolic modifiers could be hidden in the trove of snps that showed limited evidence of positive selection and that this limitation could be particularly pertinent to modifiers associated with adaptations in response to shifts of subsistence cultures .
consistent with this hypothesis , a survey of earlier studies of genome - wide or chromosome - wide positive selection using the so - called outlier approaches  in which candidate loci are identified in the extreme tails of empiric distributions ( 40)showed that gip variants have not been reported as positively selected ( 1,7 ) .
the selection of gip variants was inferred after we focused the analysis on local genomic regions and assessed the significance of integrated haplotype score using coalescent simulations ( 22 ) .
therefore , our investigation provided a proof - of - concept study for identifying causal mutations that underlie phenotypic variation of complex disease - related traits .
this approach could open new venues for improving the translation of common variant association signals into biologic mechanisms that underlie physiologic variability or disease risk .
neel hypothesized that mismatches between prior adaptations and new environments , or a  conflict of adaptations ,  could lead to changes in fitness or health risks ( 11 ) . because ancient variants could have been selected for the organism s reproductive success but not for its health or longevity , the ancient alleles could confer disease risks as selection pressures change .
therefore , studies of positively selected variants that are associated with adaptations in energy - balance regulation could point not only to novel genotype 
phenotype relationships but also to novel molecular mechanisms that mediate the potential phenotypic variation , thereby providing much - needed insight into how and which phenotypic variations in energy - balance regulation can be attributed to the selected variants . in support of the thrifty genotype hypothesis ,
human capn10 and house - mouse insulin genes have been shown to exhibit characteristics of adaptive evolution after the emergence of agricultural societies ( 41,42 ) .
conceptually , the high gip response associated with the ancestral gip haplotype could have been a beneficial energy - conserving mechanism when the food supply was irregular
. the ancestral haplotype could become deleterious in the last 10 millenniums as agricultural practice became widespread .
one possible deleterious effect of the ancestral gip haplotype in an environment that supplies abundant high - starch food resources is the hypersecretion of insulin and insulin resistance ( 43,44 ) . on the other hand , a reduced gip response associated with the derived gip haplotype
could be protective by decreasing the extent of insulin secretion in the face of oversupply of energy inputs ( 45,46 ) . in support of this speculation , it has been well documented that in the absence of modern medicine , diabetes - associated complications and , possibly , obesity posed detrimental effects on survival when human culture shifted ( 47 ) , even though type 2 diabetes is generally considered a chronic disease in modern society .
alternatively , an elevated glucose level associated with the derived haplotype may have improved the survival of fetuses if the population faced serious famine 
an event frequently experienced by agricultural societies ( 48)despite the reality that an impaired glucose - tolerance response represents a risk to both the mother and fetus under normal circumstances .
moreover , the derived haplotype could have been beneficial by reducing the obesity - promoting effect of gip ( 16 ) . in vitro and in vivo studies
have shown that gip promotes fatty liver and other obesity - associated metabolic disorders , whereas gip antagonists suppress lipid accumulation induced by a high - fat diet ( 23 ) .
therefore , the derived gip haplotype could be selected for its effects on the enteroadipocyte or the enteroinsular axis , or both .
although we speculated that the derived gip haplotype may have provided protective effects in famine - plagued agricultural societies , the observation that the derived haplotype has not been fixed in any population suggests that the selection of the derived gip haplotype(s ) ( e.g. , cycles of famine ) could be opposed ( balance selection ) as populations experienced temporal changes in selection pressure ( e.g. , resumption of population growth with stable food supply ) . alternatively
, the derived gip haplotype could simply be too young to become fixed , or the spread could be limited by the transgenerational effects associated with abnormal gestational glucose metabolism , which raise the risk of macrosomia and diabetes in the offspring ( 11 ) .
although an association between gip variants and glucose - metabolism regulation has not been reported , gipr variants were associated with glucose and insulin levels after challenge tests as well as with bmi in gwa studies that evaluated > 29,000 individuals ( 35,36 ) .
the finding that patients with a homozygous gip genotype have significantly higher glucose levels compared with those carrying an ancestral gip haplotype within the pool of gipr homozygotes suggested that there is a confounding effect stemming from interactions of gip and gipr variants , and that gip and gipr variants represent novel markers for the stratification of the capability to maintain glucose homeostasis during pregnancy .
we also speculate that the significant results observed in pregnant women could be related to the fact that the success of pregnancy has a significant impact on reproductive fitness and that a major fraction of gene  environmental selections probably occurred before birth ( 49 ) .
recent studies have corroborated this idea by showing that associations between many risk alleles and type 2 diabetes can be replicated with smaller sample sizes in patients with gestational diabetes mellitus ( 50,51 ) .
furthermore , given the evolutionary signatures at the gip locus , the plausible molecular mechanism , and the significant results in east asian women , we speculate that the gip variant  mediated phenotypic divergence could also exist in most human populations .
it is also important to note that the selection of gip variants represents a unique example in which the selection process involves regulatory variants that alter the glucose - induced gip response as well as a nonsynonymous variant that affects peptide bioactivity ( 22 ) .
thus , the gip signaling pathway could represent a hotspot for selection in recent human history and play an important role in the manifestation of phenotypic variation in energy - balance regulation among individuals .
in addition to gip variants , our study identified several cdkal1 , cyb5r4 , gad2 , and pparg variants as potential metabolic modifiers .
recent studies have shown that variants in cdkal1 , pparg , and more than two dozen genes are associated with glycemic traits in diabetic patients ( 10 ) . surprisingly
, none of the cdkal1 , cyb5r4 , gad2 , and pparg variants identified here have been implicated in earlier gwa studies , which suggests that these variants could be related to novel energy - balance regulatory mechanisms that operate at certain life stages or under specific physiologic conditions that have not been specifically investigated .
future investigations of these variants could reveal additional metabolic modifiers that have arisen recently and their contributions to phenotypic variation in normal human physiology and metabolic syndrome  related traits . in conclusion
, our data demonstrated a strong association between regulatory gip variants , and gip response and glucose metabolism , reinforcing the indication of an important role of gip signaling in diabetes - related traits from earlier gwa studies of gipr .
importantly , our study also provided a novel approach to reveal metabolic modifiers by studying consequences of previous mismatches of physiologic capabilities and environments .
 objectivealthough recent studies have shown that human genomes contain hundreds of loci that exhibit signatures of positive selection , variants that are associated with adaptation in energy - balance regulation remain elusive . 
 we reasoned that the difficulty in identifying such variants could be due to heterogeneity in selection pressure and that an integrative approach that incorporated experiment - based evidence and population genetics - based statistical judgments would be needed to reveal important metabolic modifiers in humans.research design and methodsto identify common metabolic modifiers that underlie phenotypic variation in diabetes - associated or obesity - associated traits in humans , or both , we screened 207 candidate loci for regulatory single nucleotide polymorphisms ( snps ) that exhibited evidence of gene  environmental interactions.resultsthree snps ( 
 rs3895874 , rs3848460 , and rs937301 ) at the 5 gene region of human gip were identified as prime metabolic - modifier candidates at the enteroinsular axis . 
 functional studies have shown that gip promoter reporters carrying derived alleles of these three snps ( haplotype gip1920a ) have significantly lower transcriptional activities than those with ancestral alleles at corresponding positions ( haplotype gip1920 g ) . 
 consistently , studies of pregnant women who have undergone a screening test for gestational diabetes have shown that patients with a homozygous gip1920a / a genotype have significantly lower serum concentrations of glucose - dependent insulinotropic polypeptide ( gip ) than those carrying an ancestral gip1920 g haplotype . 
 after controlling for a gipr variation , we showed that serum glucose concentrations of patients carrying gip1920a / a homozygotes are significantly higher than that of those carrying an ancestral gip1920 g haplotype ( odds ratio 3.53).conclusionsour proof - of - concept study indicates that common regulatory gip variants impart a difference in gip and glucose metabolism . 
 the study also provides a rare example that identified the common variant - common phenotypic variation pattern based on evidence of moderate gene  environmental interactions . 
RESEARCH DESIGN AND METHODS
Patients.
Genotyping.
GIP, insulin, and C-peptide enzyme-linked immunosorbent assay.
In vitro 
F
Statistical analysis.
RESULTS
SNPs at the 
The 
The 
DISCUSSION
Supplementary Material