ORIGINAL RESEARCH article

Front. Physiol., 25 September 2018

Sec. Gastrointestinal Sciences

Volume 9 - 2018 | https://doi.org/10.3389/fphys.2018.01331

A Pilot Study Investigating the Influence of Glucagon-Like Peptide-1 Receptor Single Nucleotide Polymorphisms on Gastric Emptying Rate in Caucasian Men

  • 1. School of Healthcare Science, Manchester Metropolitan University, Manchester, United Kingdom

  • 2. Division of Diabetes, Endocrinology and Gastroenterology, School of Medical Sciences, The University of Manchester, Manchester, United Kingdom

  • 3. Salford Royal Hospitals, Salford, United Kingdom

  • 4. School of Medicine, University of St Andrews, St Andrews, United Kingdom

  • 5. School of Biomedical Sciences, Ulster University, Coleraine, United Kingdom

Article metrics

View details

13

Citations

3,1k

Views

1k

Downloads

Abstract

Gastric emptying rate in humans is subject to large individual variability, but previous research on the influence of genetics is scarce. Variation in the glucagon-like peptide-1 receptor (GLP1R) gene is a plausible candidate gene to partially explain the high variance. This study aimed to investigate the influence of genetic variation in the GLP1R gene on gastric emptying rate of a glucose solution in humans. Forty eight healthy Caucasian males took part in this investigation. Gastric emptying rate of a 6% glucose solution was assessed using the 13C breath test method and a venous blood sample was obtained from each participant. Participants were genotyped for 27 Tag single nucleotide polymorphisms (SNPs) in the GLP1R locus using Sequenom MassARRAY iPLEX GOLD analysis and MALDI-TOF mass spectrometry. The time at which maximal emptying rate occurred (Tlag) was faster in participants with the CC genotype than in TT and TC genotypes for SNP rs742764: [median (quartiles) CC, 35 (30–36) min vs. TT, 43 (39–46) min, and TC, 41 (39–45) min; P < 0.01]. Tlag was also slower in participants with the AA genotype compared to the TT and TA genotypes for SNP rs2254336: [AA, 43 (39–49) min vs. TT, 36 (34–41) min, and TA, 39 (35–42) min; P < 0.05]. Analysis by phenotype also showed differences in half-emptying time (T12) and Tlag for SNPs rs9283907, rs2268657, and rs2254336. Several neighboring Tag SNPs within the GLP1R gene were found to be associated with gastric emptying rate, and should be further investigated.

Introduction

An increasing number of genetic variants within the genes coding for gastrointestinal hormones and their receptors are being found to be associated with obesity phenotypes, appetite, and food regulation. As many of the gastrointestinal hormones involved in appetite regulation also play a role in the regulation of gastric emptying, it is possible that genetic variations in gastrointestinal hormones associated with obesity may affect gastric emptying rate. However, research on the potential influence of genetics on gastric emptying rate is scarce. Acosta et al. (2014) reported a common genetic variant rs17782313 in the Melanocortin 4 receptor (MC4R) gene, which has previously been found to be strongly associated with common obesity (Loos et al., 2008; Vogel et al., 2011), to be associated with reduced gastric emptying rate and satiation. In addition, Cremonini et al. (2005) found an association between the 779T > C polymorphism in the cholecystokinin (CCK) gene and slower gastric emptying rate. On the other hand, Jones et al. (2010) found no effect of common genetic polymorphisms of the CCK or CCK-1 receptor genes on gastric emptying rate.

Genetic variation within a gene related to the action of the gastrointestinal hormone glucagon-like peptide-1 (GLP-1) presents a plausible area of investigation, as GLP-1 has been shown to influence gastric emptying rate (Wettergren et al., 1993; Wishart et al., 1998). Furthermore, large inter-individual variations in gastric emptying rate of a glucose solution as well as the large differences in GLP-1 hormone responses to carbohydrate ingestion have been observed (Yau et al., 2014, 2017). GLP-1 exerts its effects via a G-protein coupled receptor called the GLP-1 receptor (GLP1R). Stimulation of this receptor by GLP-1 triggers cAMP production as the primary signal transduction pathway (Mayo et al., 2003). The GLP1R recognizes GLP-1 specifically despite the hormone having strong sequence homology to other hormones within the glucagon-related family of peptides, and furthermore, does not bind to a number of other related peptides, including secretin (Fehmann et al., 1994).

Genetic polymorphisms of the GLP1R gene have previously been investigated in relation to insulin secretion (Sathananthan et al., 2010) and the pathogenesis of diabetes (Tanizawa et al., 1994; Tokuyama et al., 2004; Beinborn et al., 2005). Only one study has previously investigated the influence of GLP1R genetic variation on gastric emptying rate, but in mice (Kumar et al., 2008). The presence of a non-synonymous cysteine to tyrosine substitution at amino acid 416 of the GLP1R gene found in CAST strain mice was associated with reduced GLP1R expression and significantly faster (by 20%) gastric emptying rate compared to B6 strain mice. This was also seen for the congenic strains where gastric emptying rate was significantly higher in B6.CAST-17 congenic mice compared to homozygous B6 controls. Furthermore, administration of the GLP1R antagonist extendin-(9-36) resulted in no increase in gastric emptying rate compared to a 10% increase in the homozygous B6 control mice. The GLP1R gene is therefore a plausible candidate gene for a genetic association study on gastric emptying rate in humans. The human GLP1R gene consists of 13 exons interrupted by 12 introns (Wilmen et al., 1998) and is situated on chromosome 6, band p21.1 (Stoffel et al., 1993). A major transcription start point and a minor transcription start point 42 base pairs (bp) and 360 bp upstream of the translation initiation site, respectively, have been reported (Lankat-Buttgereit and Goke, 1997). Three putative Sp1 binding sites have also been located in the proximal 5′ flanking sequence at -108, -173, and -389 bp from the translation initiation codon (Lankat-Buttgereit and Goke, 1997). The receptor is 463 amino acids in length and is highly conserved between species with 90% being identical to rat GLP1R (Dillon et al., 1993) and approximately 95% homology with mice GLP1R. The primary aim of this study was to investigate the influence of genetic variation in the GLP1R gene on gastric emptying rate of a glucose solution in humans.

Materials and Methods

Ethical Approval

Forty eight British and South-Western European Caucasian males aged 18–35 years (mean ± SD, age 23 ± 5 years, height 178 ± 7 cm, body mass 75.8 ± 11.2 kg, body mass index 23.9 ± 3.3 kg.m-2, and estimated body fat percentage 19 ± 6%) volunteered for this study. All participants were non-smokers, were not taking medication with any known effect on gastrointestinal function, and had no known history of chronic gastrointestinal disease or other medical conditions as assessed by a medical screening questionnaire. All participants provided written informed consent prior to participation. The study was granted prior ethical approval from the Manchester Metropolitan University Ethical Advisory committee and was performed in accordance with the Declaration of Helsinki (2008).

Experimental Trial

All participants attended the laboratory for one experimental trial. In the 24 h preceding the trial, participants were requested to refrain from alcohol and caffeine consumption, and the performance of strenuous physical activity. Participants arrived at the laboratory in the morning following an overnight fast from 2100 h, with the exception of drinking 500 mL of water 90 min prior to arrival at the laboratory, in an attempt to ensure a consistent and adequate hydration status.

Upon arrival at the laboratory, participants were asked to completely empty their bladder. Body mass was then recorded. Following this, a whole blood sample was obtained by venipuncture of an antecubital vein. Once this procedure was completed, participants ingested 595 mL of a 6% glucose solution as quickly as they were comfortably able over a maximum period of 2 min. Test solutions were prepared fresh on the morning of the trial containing 39.6 g glucose monohydrate 1 dissolved in water and made up to a volume of 600 mL. In addition, solutions contained 100 mg of 13C sodium acetate (Cambridge Isotopes Laboratories, Inc., Andover, MA, United States) for the assessment of gastric emptying rate. The drink was served at ambient temperature and a 5 mL sample retained for later analysis of osmolality. Participants remained seated throughout a 60 min post-ingestion sampling period where end expiratory breath samples were collected at baseline then every 10 min after ingestion into foil bags (Wagner Analyzen-Technik, Breman, Germany) by exhalation through a mouthpiece. Bags were then sealed with a plastic stopper and stored for later analysis. Breath samples were analyzed for the ratio of 13CO2 to 12CO2 by non-dispersive infrared spectroscopy, with the difference in the ratio from baseline to post-breath samples expressed as delta over baseline (DOB). Half emptying time (T½) and the time at which maximal emptying rate occurred (Tlag) were calculated using the manufacturer’s integrated software evaluation embedded with equations previously described (Ghoos et al., 1993). Each participant’s own physiological CO2 production assumed as 300 mmol CO2 per m2 body surface per hour was set as default and body surface area was calculated by the integrated software according to a previously described formula (Haycock et al., 1978). Following all sample collections at 60 min, participants were free to leave the laboratory.

Genotyping

Genomic DNA was extracted from 3 mL of whole blood using Flexigene DNA Kit (Qiagen, West Sussex, United Kingdom), according to the manufacturer’s instructions, except for the final step of resuspension in 300 μL of nuclease-free ultrapure water. Extracted DNA was subsequently quantified using a NanoDrop 2000 (Thermo Scientific, 116 Loughborough, United Kingdom) and normalized to a concentration of 40 ng/μL then stored at -20°C before working concentrations of 20 ng/μL were prepared in 96-well plate format prior to use.

Twenty-eight tag single nucleotide polymorphisms (SNPs) in the GLP1R locus incorporating 10,000 bp upstream and downstream of the major transcription initiation site and the last exon, respectively (Chr6: 39114595….39173498) were selected from HapMap2 (HapMap Data release 27/phase II + III, Feb09, on NCBI B36 assembly, dbSNP b126). The Tagger algorithm for multi-marker tagging with r2 > 0.8 and minor allele frequency (MAF) > 0.1 was used (de Bakker et al., 2005). Furthermore, three additional non-synonymous SNPs were selected based on previous literature. One of two SNPs found to be associated with insulin secretion in response to exogenous infusion of GLP-1 by Sathananthan et al. (2010), which was not already in the generated list of 28 tag SNPs was genotyped. The other two additional SNPs genotyped were two in very close proximity (two amino acids upstream and three amino acids downstream) to the equivalent locus of the SNP found in mice to be associated with gastric emptying rate (Kumar et al., 2008). All 31 SNPs were genotyped using Sequenom MassARRAY iPLEX GOLD analysis. Forward and reverse primers as well as extension probes were designed using Sequenom Assay Design Suite (v1.0) which produced two appropriate assay plexes; plex 1 containing 24 SNPs and plex 2 containing seven SNPs (Table 1). Primers were purchased from Metabion International AG (Martinsreid, Germany).

Table 1

SNP IDPos. from major transcription site (bp)Forward primerReverse primerExtension probe
1rs7738586-9,033ACGTTGGATGACACCCAGACTGACGTTATCACGTTGGATGTTGTGCAAAAGCAGCCCAAGTCAAAGTGATTGTCACCATAAG
2rs9380825-5,610ACGTTGGATGTGAGCCAGGAAGTGATGTTCACGTTGGATGTCTCATCCAGGCAGCAAGTAcccCTGCACCAGAGCCCCT
3rs9296274-1,455ACGTTGGATGTTTGGATATCGTGGCTGGAGACGTTGGATGTGCCTGGGTCTTCTAGCTTCgggtGTGGAAATCTTGGACCA
4rs9266742,930ACGTTGGATGACTCACACATACCTGGGAACACGTTGGATGATACAGACAGGTAGTCTGAGtagcAGGGAGTAGGCTATATGA
5rs22686573,969ACGTTGGATGGTTCTGCCGTCCATAAAATGACGTTGGATGTACAGGGCTTGAGAAGTCACggaagTGGGCATATCATTCTTCTCA
6rs132023696,125ACGTTGGATGTGATCCACCAGGACTTGCTCACGTTGGATGTAACAGCTGCAAAGGTGTTGgacaTCCTCAGCTGTGGCTAAT
7rs37997076,937ACGTTGGATGTTTGGTTGCTGTGTCAGAGGACGTTGGATGGCACTCACTTACAGATGCACAGATGCACTCAACACA(Inosine)C
8rs103054327,057ACGTTGGATGTCTTTGTAGCCCTGAACGCCACGTTGGATGAGTCCTTCAGATCAGTGACCCCGCACACCTTGCGA
9rs928390710,130ACGTTGGATGGCTCCTATCATCACACCTTGACGTTGGATGCCAGAGCATAACCTCATGCCaagaA(Inosine)CAACTGGCCCAGAA
10rs74276413,261ACGTTGGATGTCTCAGCTTCTGCATCTGTAACGTTGGATGTGTGGAGAGCTGCTCATGAAggagcGCTGCTCATGAATCCATTA
11rs225433616,262ACGTTGGATGCTGGTCTAAAAGGAGTACACACGTTGGATGAAGGTAGGAGCTGGGTATTCggAGGCTGCATACGACCA
12rs91016316,847ACGTTGGATGTGAGCCTCAGCCCAGAATAGACGTTGGATGCAGGGATAGCCCTCAGAATGtATGGGGAGGAAGGGG
13rs376546717,022ACGTTGGATGAACCCCGCCTCAACTCACTCACGTTGGATGTGCAGAAGGACAACTCCAGCAACTCACTCT(Inosine)TCCCCT
14rs692376117,499ACGTTGGATGTTCTCTGCTCTGGTTATCGCACGTTGGATGGAGTTAGGATGAAGCAGCCCGGGCCACCTTACCTGAAGC
15rs776666319,209ACGTTGGATGAGCAATAGGTCTGCATGTGGACGTTGGATGCTGGTCTCCAATCTCTGGCCCTAGCTAATTGAGAGGC
16rs93244325,761ACGTTGGATGGTGAATGAAGGAGTGGCAAGACGTTGGATGTCCTCCCTAATCTGCCATTGggattGTGTGGAACAGGAAAACTC
17rs226864626,943ACGTTGGATGCCAACTGTGTCAGAGTCCTAACGTTGGATGGTGATGCCAGGAGGCCTTGcccaTTCCACTTGCACATGAA
18rs230061432,400ACGTTGGATGTCCAAACCTAGGGCAGGTTCACGTTGGATGCCCTGCTAAAATTCTTATTTCTCTTTCTGATCTTCAGTGTT
19rs226864133,693ACGTTGGATGCTGGGTCCTCTAAGACCTGTACGTTGGATGCAAAGAGTGGCCCATAAATGggggAAGACCTGTCCCAGGA
20rs226864033,811ACGTTGGATGTGCACTTCCTCGTTTGCATCACGTTGGATGTCCTCCTCCACTGCCATATCCACTGCCATATCCTCAAAATGA
21rs226863934,049ACGTTGGATGGATAGAGAAGTGAGAAACGGACGTTGGATGATGAGGAGCAGAGGCCTGTAggcaCTGCTGCCACCTTGTCATCT
22rs220694234,866ACGTTGGATGAATTGGGAAGCTCATTCACCACGTTGGATGGGCAAGTCATTTTGCCTCCCggcGCTCATTCACCTTCATTTAC
23rs289442036,523ACGTTGGATGAACAGGGATCCTGGCTGACACGTTGGATGAGTGAGGGCTTCTCAACTGCACTGCAGTGTCTCTCT
24rs199796313†37,124ACGTTGGATGTCCTTTTCCCATGGAAGGTCACGTTGGATGTGGATGTGCAAGTGCTCAAGGGTCCAGCTGGAATTT
25rs200691429†37,140ACGTTGGATGAAGGTCCAGCTGGAATTTCGACGTTGGATGTGGATGTGCAAGTGCTCAAGgGGAAGAGCTGGGAGC
26rs471421139,117ACGTTGGATGAAGGCACCCCTTATTTGCTGACGTTGGATGACTTGCACCAGCACTGTTTCccccTTTGCTGTCTCTTCGT
27rs1030552539,577ACGTTGGATGAATGGCACTGCACTCTTTCCACGTTGGATGCATTGCATTCAATAGTTCCCaATTCAATAGTTCCCAGACCT
28rs929629140,356ACGTTGGATGGATGGTGAAAGTGTCATCTCACGTTGGATGAAGACAAGGATGAATGAAGgTGAATGAAGTACCAGTGT
29rs996888643,905ACGTTGGATGACAGTGAGGTTTTCCCCATCACGTTGGATGTCATGTAGTCCAGCTTGTGCGCTTGTGCTGCTAGTT
30rs214373344,923ACGTTGGATGCATCTAATCGATGGGTAGCACGTTGGATGCAGAACCCTTCTGAACCTTCGAAATTGAATTTACAGCTTTAATAAA
31rs929629246,417ACGTTGGATGTCACAATATGTTTGGCACTGACGTTGGATGGCTTTGTTTTGCAGAGCTTGTGGCACTGCCAAACT

List of SNPs analyzed and primers used.

Unshaded SNPs analyzed in plex 1, shaded SNPs in plex 2. Additional missense SNP associated with insulin secretion in response to exogenous infusion of GLP-1 by Sathananthan et al. (2010). †Additional missense SNPs in close proximity to equivalent SNP associated with gastric emptying rate in mice (Kumar et al., 2008).

Briefly, an initial locus-specific amplification was performed using polymerase chain reaction (PCR; GeneAmp PCR System 9700, Applied Biosystems) carried out in a total volume of 5 μL containing 40 ng of DNA and final concentrations of 1.25x buffer, 1.0 mM MgCl2, 500 μM dNTP mix, 100 nM of each forward and reverse primers, and 0.1 U/μL Hotstart Taq. PCR conditions consisted of an initial denaturing step at 94°C for 5 min, followed by 40 cycles of denaturing at 94°C for 20 s, annealing at 56°C for 30 s and extension at 72°C for 1 min, and then a final extension step at 72°C for 3 min. Unincorporated dNTPs within the PCR products were dephosphorylated by incubation for 40 min at 37°C followed by 5 min at 85°C with 2 μL Shrimp Alkaline Phosphatase (SAP) mixture composed of 1.7 μL 1x hME buffer and 0.3 μL SAP 1.7 U/μL. A single base extension step was then performed with the addition of 2 μL reaction mix containing 50 μM of each ddNTP, 3.3:6.6 μM extension probe from low:high mass (2.8:5.6 μM for plex 2), 1.0x buffer, and 1.25 U Thermo Sequenase. Probe extension conditions consisted of a two-step 200 short cycle program involving initial denaturing at 94°C for 30 s, 40 cycles of denaturing at 94°C for 5 s, annealing at 52°C for 5 s, and extension at 80°C for 5 s where the annealing and extension steps also cycled five times, before a final extension at 72°C for 3 min. Resulting products were diluted with 20 μL nuclease-free ultrapure water and desalted with 6 mg resin by gentle inversion for 10 min before centrifugation at 4000 rpm for 5 min. Products were subsequently nano-dispensed (Samsung Sequenom MassARRAY Nanodispenser) onto a 384-element SpectroCHIP II bioarray and analyzed by MALDI-TOF mass spectrometry.

Data and Statistical Analysis

Genotype results were visually checked and appropriate manual genotype assignment was made upon spectrum inspection where no automated genotype calling had successfully taken place. Differences in gastric emptying T½ and Tlag were examined between genotype and phenotype groups. For each SNP analyzed, it was hypothesized that there was a significant difference in gastric emptying rate between genotype groups or phenotype groups. Normality tests indicated that the majority of data were not normally distributed. Furthermore, as group sizes were unequal, non-parametric statistical analysis comprising Kruskal–Wallace and Mann–Whitney U-tests were utilized accordingly. Where appropriate, post hoc tests of Mann–Whitney U with Bonferroni correction were applied. All data were analyzed using SPSS Statistics for Windows (IBM, New York, NY, United States). Statistical significance was accepted at the 5% level and results presented as median and quartiles unless stated otherwise.

Results

SNP Genotyping

Twenty-seven of the 31 SNPs were successfully analyzed for variants. No variants occurred among the participants for three SNPs, Tag SNP 3 within the promoter region (rs9296274) and two of the three selected missense polymorphism, SNPs 24 (rs199796313) and 25 (rs200691429). Genotyping failed in all participants for SNP 13 (rs3765467), one of the three additional missense SNPs. The occasional failure to successfully genotype one participant occurred in five SNPs, reducing the total participant number (n) to 47 instead of 48. Frequencies of each genotype are shown in Table 2.

Table 2

Genotype
SNP IDHomozygous minor allele
Heterozygous
Homozygous major allele
GenotypenGenotypenGenotypenTotal
1rs7738586AA0CA10CC3848
2rs9380825AA5AG28GG1548
4rs926674TT3TC6CC3948
5rs2268657GG9AG27AA1147
6rs13202369GG2AG21AA2548
7rs3799707TT3GT19GG2648
8rs10305432CC3CT19TT2648
9rs9283907AA1AG9GG3848
10rs742764CC8TC25TT1447
11rs2254336TT9TA23AA1648
12rs910163CC3TC18TT2748
14rs6923761AA7GA28GG1348
15rs7766663GG10GT23TT1548
16rs932443GG5AG21AA2248
17rs2268646AA0AG10GG3848
18rs2300614TT5CT21CC2248
19rs2268641AA4AG21GG2247
20rs2268640CC3TC18TT2748
21rs2268639TT3TA22AA2248
22rs2206942AA6AG24GG1848
23rs2894420AA9AG27GG1147
26rs4714211GG7AG27AA1347
27rs10305525AA0CA9CC3948
28rs9296291CC3TC17TT2848
29rs9968886AA0GA12GG3648
30rs2143733GG7GT27TT1448
31rs9296292CC4CT20TT2448

SNP genotype frequencies.

Gastric Emptying Rate

Mean ± SD gastric emptying rate for all participants was 68 ± 16 and 41 ± 8 min for T½ and Tlag, respectively. Median (quartiles) gastric emptying rate for all participants was 63 (55–78) and 40 (36–44) min for T½ and Tlag, respectively.

By Genotype

Results for Tlag and T½ according to genotype are shown in Tables 3, 4, respectively. Differences in median gastric emptying Tlag were seen for SNP 10 rs742764 and SNP 11 rs2254336. For SNP 10 rs742764, gastric emptying Tlag was faster in genotype CC compared to genotype TT (P = 0.006) and TC (P = 0.006) by 15%. T½ was also tending to significance with differences of 18% (P = 0.061). For SNP 11 rs2254336, gastric emptying Tlag was slower in genotype AA compared to genotype TT (P = 0.04) and TA (P = 0.036) by 19 and 10%, respectively. T½ showed a slight trend to be slower than both other groups by 20 and 15% but this did not reach statistical significance (P = 0.138). No significant differences in gastric emptying rate were seen between genotypes for all other SNPs, although T½ tended toward significance for SNP 9 rs9283907 (P = 0.054) and Tlag tended to significance for SNP 5 rs2268657 (P = 0.087) and SNP 15 rs7766663 (P = 0.076).

Table 3

SNP IDHomozygous minor allele
Heterozygous
Homozygous major allele
P-value
GenotypeMedianQuartilesGenotypeMedianQuartilesGenotypeMedianQuartiles
1rs7738586CA3734–41CC4037–450.274
2rs9380825AA4261–80AG4137–45GG3736–400.349
4rs926674TT3737–41TC3836–40CC4136–460.448
5rs2268657GG4237–44AG3837–46AA3734–390.087
6rs13202369GG3431–36AG4137–49AA3936–420.295
7rs3799707TT4236–43GT4035–43GG4037–450.696
8rs10305432CC3635–37CT4139–43TT4035–460.267
9rs9283907AA2929–29AG3735–41GG4137–440.158
10rs742764CC3530–36TC4139–45TT4137–460.008
11rs2254336TT3634–41TA3935–42AA43†39–490.031
12rs910163CC3737–39TC4036–45TT4136–440.788
14rs6923761AA4341–48GA4037–43GG3734–410.222
15rs7766663GG3631–40GT4038–44TT4237–480.076
16rs932443GG3736–41AG4038–45AA4135–440.634
17rs2268646AG4037–48GG4035–430.493
18rs2300614TT3635–37CT4037–44CC4236–460.210
19rs2268641AA3737–38AG4035–44GG4135–440.724
20rs2268640CC3737–39TC4034–45TT4137–440.741
21rs2268639TT3737–39TA4035–44AA4136–440.816
22rs2206942AA3937–51AG4036–43GG4132–460.878
23rs2894420AA4339–49AG3936–42GG3736–480.470
26rs4714211GG3633–39AG4037–45AA4137–440.244
27rs10305525CA3736–41CC4036–450.369
28rs9296291CC3737–39TC4035–45TT4137–430.812
29rs9968886GA3736–45GG4136–430.543
30rs2143733GG3737–48GT3935–43TT4239–480.549
31rs9296292CC3734–38CT3934–44TT4138–450.166

Gastric emptying Tlag results according to genotype.

Values are minutes. Significantly faster than other two genotypes (P < 0.01). †Significantly slower than other two genotypes (P < 0.05). Bold values indicate significant difference.

Table 4

SNP IDHomozygous minor allele
Heterozygous
Homozygous major allele
P-value
GenotypeMedianQuartilesGenotypeMedianQuartilesGenotypeMedianQuartiles
1rs7738586CA6256–72CC6455–790.919
2rs9380825AA7161–80AG6456–80GG6055–700.650
4rs926674TT5956–62TC6055–71CC6456–820.569
5rs2268657GG7259–83AG6156–73AA5953-670.203
6rs13202369GG7058–81AG6254–77AA6456–730.935
7rs3799707TT6157–66GT5954–83GG6659–780.452
8rs10305432CC5956–71CT6557–75TT6255–800.852
9rs9283907AA4141–41AG5654–62GG6759–820.054
10rs742764CC5453–59TC6656–80TT6659–810.061
11rs2254336TT5953–68TA6254–75AA7162–850.138
12rs910163CC5956–60TC6654–82TT6257–730.431
14rs6923761AA7266–81GA6256–74GG5954–830.297
15rs7766663GG5753–66GT6254–79TT7161–810.168
16rs932443GG6059–68AG6554–83AA6257–720.820
17rs2268646AG6757–87GG6354–760.485
18rs2300614TT5453–59CT6554–83CC6658–730.120
19rs2268641AA5754–59AG6554–80GG6358–730.283
20rs2268640CC5956–60TC6654–82TT6257–730.453
21rs2268639TT5956–60TA6555–79AA6356–730.475
22rs2206942AA6057–86AG6555–78GG6355–730.959
23rs2894420AA7264–92AG6255–72GG5954–810.187
26rs4714211GG5454–60AG6456–75AA6859–890.233
27rs10305525CA5955–68CC6455–790.412
28rs9296291CC5956–60TC6654–83TT6256–720.386
29rs9968886GA6156–92GG6554–740.934
30rs2143733GG5955–77GT6455–79TT6559–730.791
31rs9296292CC5652–59CT6554–81TT6660–770.125

Gastric emptying T½ results according to genotype.

Values are minutes.

By Phenotype

Results for Tlag and T½ according to phenotype are shown in Tables 5, 6, respectively. A significant effect of the minor allele on median gastric emptying T½ was seen for SNP 9 rs9283907. Participants with one or two A alleles had faster T½ than participants homozygous for the allele G by 18% (P = 0.033). A significant effect of the minor allele on median gastric emptying Tlag was also seen for SNP 5 rs2268657 and SNP 11 rs2254336. Participants with one or two G alleles had slower Tlag than participants homozygous for the allele A by 11% (P = 0.028), and participants with one or two T alleles had faster Tlag than participants homozygous for the allele A by 13% (P = 0.012), respectively. No differences in gastric emptying rate were seen between phenotypes for all other SNPs, though T½ tended toward significance for SNP 11 rs2254336 (P = 0.055) and SNP 15 rs7766663 (P = 0.097).

Table 5

SNP IDHomozygous major allele
Heterozygous/homozygous minor allele
P-value
PhenotypenMedianQuartilesPhenotypenMedianQuartiles
1rs7738586CC384037–45A103734–410.274
2rs9380825GG153736–40A334137–440.157
4rs926674CC394136–46T93736–400.239
5rs2268657AA113734–39G364137–450.028
6rs13202369AA253936–42G234136–480.627
7rs3799707GG264037–45T224034–430.396
8rs10305432TT264035–46C224137–420.844
9rs9283907GG384137–44A103633–410.137
10rs742764TT144137–46C334035–430.492
11rs2254336AA164339–49T323934–410.012
12rs910163TT274136–44C213936–440.950
14rs6923761GG133734–41A354037–440.197
15rs7766663TT154237–48G333935–420.136
16rs932443AA224135–44G264036–440.983
17rs2268646GG384035–43A104037–480.493
18rs2300614CC224236–46T263935–420.330
19rs2268641GG224135–44A253936–420.856
20rs2268640TT274137–44C213935–440.546
21rs2268639AA224136–44T253935–420.701
22rs2206942GG184132–46A304036–440.806
23rs2894420GG113736–48A364036–430.930
26rs4714211AA134137–44G344035–440.497
27rs10305525CC394036–45A93736–410.369
28rs9296291TT284137–43C204035–440.908
29rs9968886GG364136–43A123736–450.543
30rs2143733TT144239–48G343935–440.276
31rs9296292TT244138–45C243834–430.107

Gastric emptying Tlag results according to phenotype.

Values are minutes. Significantly faster emptying rate compared to other phenotype (P < 0.05). Bold values indicate significant difference.

Table 6

SNP IDPhenotypenMedianQuartilesPhenotypenMedianQuartilesP-value
1rs7738586CC386455–79A106256–720.919
2rs9380825GG156055–70A336456–800.367
4rs926674CC396456–82T95954–650.355
5rs2268657AA115953–67G366456–780.119
6rs13202369AA256456–73G236254–800.757
7rs3799707GG266659–78T226054–760.218
8rs10305432TT266255–80C226555–760.836
9rs9283907GG386759–82A105553–610.033
10rs742764TT146659–81C336254–770.478
11rs2254336AA167162–85T326054–740.055
12rs910163TT276257–73C216454–800.950
14rs6923761GG135954–83A356457–750.516
15rs7766663TT157161–81G336154–770.097
16rs932443AA226257–72G266554–820.702
17rs2268646GG386354–76A106757–870.485
18rs2300614CC226658–73T266254–790.521
19rs2268641GG226358–73A256254–770.733
20rs2268640TT276257–73C216454–800.755
21rs2268639AA226555–79T256254–770.806
22rs2206942GG186355–73A306355–790.774
23rs2894420GG115954–81A366456–740.715
26rs4714211AA136859–89G346254–730.385
27rs10305525CC396455–79A95955–680.412
28rs9296291TT286256–72C206554–810.875
29rs9968886GG366554–74A126156–920.934
30rs2143733TT146559–73G346354–790.691
31rs9296292TT246660–77C246054–780.239

Gastric emptying T½ results according to phenotype.

Values are minutes. Significantly faster emptying rate compared to other phenotype (P < 0.05). Bold values indicate significant difference.

Body Mass Index and Body Fat Percentage

Median (quartiles) BMI for all participants was 23.1 (21.6–25.3) kg.m-2. No effect of genotype on BMI was found for all SNPs although SNP 15 rs7766663 tended to significance [GG, 24.5 (21.6–26.2) kg.m-2 vs. GT, 24.0 (22.2–26.0) kg.m-2 vs. TT, 22.1 (21.1–23.3) kg.m-2; P = 0.091]. Analysis by phenotype, however, revealed median BMI was significantly higher in participants with one or two minor alleles compared to homozygotes of the major allele for SNP 12 rs910163 [TT, 22.2 (21.3–24.2) kg.m-2 vs. CC/CT, 24.6 (22.6–26.7) kg.m-2; P = 0.039] and SNP 15 rs7766663 [TT, 22.1 (21.1–23.3) kg.m-2 vs. GG/GT, 24.0 (22.1–26.3) kg.m-2; P = 0.028].

Median (quartiles) body fat for all participants was 18.3 (14.1–22.2)%. No effect of genotype nor phenotype was seen for any SNPs, although SNP 28 rs9296291 tended to significance for phenotype [TT, 17.2 (13.5–19.7)% vs. CC/CT, 20.8 (15.5–24.9)%; P = 0.061].

Discussion

The results of this study showed several Tag SNPs of the GLP1R gene are significantly associated with the gastric emptying rate of a glucose solution in healthy young Caucasian men. Two neighboring polymorphisms, SNPs 10 rs742764 and 11 rs2254336, were found to be significantly associated with gastric emptying rate by genotype. In addition, SNPs 5 rs2268657, 9 rs9283907, and 15 rs7766663 tended to significance. Furthermore, three SNPs, 5 rs2268657, 9 rs9283907, and 11 rs2254336, were found to be significantly associated with gastric emptying rate by phenotype. SNP 15 rs7766663 also tended to significance. Thus, significant associations between genetic variation and gastric emptying rate were found for four Tag SNPs by one or more measures of genetic association. Tag SNP 5 rs2268657 is situated in intron one and the other three 9, 10, 11 are situated in intron 3.

The aforementioned Tag SNPs identified to be associated with gastric emptying rate signify a region where a causative variant is most likely to reside (Xia and Grant, 2013). As three of the associated SNPs are neighboring Tag SNPs, this presents a genomic region that warrants further investigation with particularly high interest. Approximately 3000 bp exist between the locations of SNP 9 rs9283907 and 10 rs742764 and similarly between SNP 10 rs742764 and 11 rs2254336. Spanning this whole 6132 bp “hot spot” region, a total of 129 SNPs have been sequenced to date. No known functional SNPs within the GLP1R gene have currently been identified within this particular region or within introns 1, 3, or 5. However, it is widely known that SNPs in regulatory elements residing within intronic regions can alter silencing, enhancer, or splicing events (Chorev and Carmel, 2012; Lee et al., 2012). Further work by in silico analysis or multi-array analysis of all 129 known SNPs within this region to narrow down on the precise SNP or SNPs responsible for the observed differences in gastric emptying rate should therefore be conducted. The genetic variants surrounding SNP 5 rs2268657 in intron one should also be further investigated as SNPs within intron one of several genes have been shown to influence gene transcription events (Murani et al., 2009; Berulava and Horsthemke, 2010).

Variants in the proximity of SNP 15 rs7766663 also provide a directive area of further research in gastric emptying regulation as it tended toward significance but was also significantly associated with BMI by phenotype. This association may signal a link between gastric emptying rate and BMI but further participants are required to confirm these concurrent associations. Indeed, further potential “links” or associations may also be identified with a much larger sample size which will increase the power of future studies.

The two additional missense SNPs in close proximity to the locus of the variant seen in the mice model by Kumar et al. (2008) did not show any variants in the sample population of this present study. It may be that the variant is rare or non-existent in the Caucasian population. The MAFs for these two missense SNPs are unknown. The other additional missense SNP selected for its previous association with insulin secretion (Sathananthan et al., 2010) also showed no variants in this participant group. This was mostly likely due to its small MAF of 0.0646 indicating the SNP is somewhat rare. Tag SNPs 6 rs13202369 and 14 rs6923761, located in intron 1 and exon 5, respectively, have previously been reported to alter insulin secretion responses to intravenous GLP-1 (Vella et al., 2009; Sathananthan et al., 2010). These were not found to be associated with gastric emptying rate in this study, however.

Conclusion

The results of this targeted gene study to investigate the potential influence of GLP1R genetic variation on gastric emptying rate in humans revealed several Tag SNPs to be associated with gastric emptying rate of a glucose solution in healthy Caucasian men. This suggests that genetic variation within the GLP1R gene may influence gastric emptying rate in humans. Further work should be undertaken to identify the precise SNP or SNPs responsible and functional analysis conducted. In addition, this experiment should be replicated with the use of foods more representative of a typical meal. Furthermore, this association study should be repeated with a larger population sample to independently confirm the detected associations between GLP1R genetic variation and gastric emptying rate.

Statements

Author contributions

Experimental trials and DNA extractions were performed at Manchester Metropolitan University. All other genotyping work was performed at CIGMR in Manchester University. AY, GE, and JA conceived and designed the experiments. AY and JA performed the experiments. AY analyzed the data and wrote the paper with contributions from GE, JA, JM, RM, and WG. All authors read and approved the final manuscript and agreed to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.

Funding

AY was supported by a Manchester Metropolitan University Ph.D. studentship. This research received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors.

Acknowledgments

The authors would like to acknowledge Mr. Dave Maskew (Manchester Metropolitan University) for his technical support in the laboratory, the staff at Salford Royal Hospital’s Gastrointestinal Physiology department for their cooperation with breath sample analysis, and Miss Hazel Platt (CIGMR, Manchester University) for her technical support with the genetics analysis.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  • 1

    AcostaA.CamilleriM.ShinA.CarlsonP.BurtonD.O’NeillJ.et al (2014). Association of melanocortin 4 receptor gene variation with satiation and gastric emptying in overweight and obese adults.Genes Nutr.9:384. 10.1007/s12263-014-0384-8

  • 2

    BeinbornM.WorrallC. I.McBrideE. W.KopinA. S. (2005). A human glucagon-like peptide-1 receptor polymorphism results in reduced agonist responsiveness.Regul. Pept.13016. 10.1016/j.regpep.2005.05.001

  • 3

    BerulavaT.HorsthemkeB. (2010). The obesity-associated SNPs in intron 1 of the FTO gene affect primary transcript levels.Eur. J. Hum. Genet.1810541056. 10.1038/ejhg.2010.71

  • 4

    ChorevM.CarmelL. (2012). The function of introns.Front. Genet.3:55. 10.3389/fgene.2012.00055

  • 5

    CremoniniF.CamilleriM.McKinzieS.CarlsonP.CamilleriC. E.BurtonD.et al (2005). Effect of CCK-1 antagonist, dexloxiglumide, in female patients with irritable bowel syndrome: a pharmacodynamic and pharmacogenomic study.Am. J. Gastroenterol.100652663. 10.1111/j.1572-0241.2005.41081.x

  • 6

    de BakkerP. I.YelenskyR.Pe’erI.GabrielS. B.DalyM. J.AltshulerD. (2005). Efficiency and power in genetic association studies.Nat. Genet.3712171223. 10.1038/ng1669

  • 7

    DillonJ. S.TanizawaY.WheelerM. B.LengX. H.LigonB. B.RabinD. U.et al (1993). Cloning and functional expression of the human glucagon-like peptide-1 (glp-1) receptor.Endocrinology13319071910. 10.1210/endo.133.4.8404634

  • 8

    FehmannH. C.JiangJ. W.SchweinfurthJ.WheelerM. B.BoydA. E.GokeB. (1994). Stable expression of the rat glp-1 receptor in cho cells - activation and binding characteristics utilizing glp-1(7-36)-amide, oxyntomodulin, exendin-4, and exendin(9-39).Peptides15453456. 10.1016/0196-9781(94)90204-6

  • 9

    GhoosY. F.MaesB. D.GeypensB. J.MysG.HieleM. I.RutgeertsP. J.et al (1993). Measurement of gastric-emptying rate of solids by means of a carbon-labeled octanoic-acid breath test.Gastroenterology10416401647. 10.1016/0016-5085(93)90640-X

  • 10

    HaycockG.SchwartzG.WisotskyD. (1978). Geometric method for measuring body surface area: a height-weight formula validated in infants children and adults.J. Paediatr.936266. 10.1016/S0022-3476(78)80601-5

  • 11

    JonesR. B.PaytonA.OilierW. E.DockrayG. J.ThompsonD. G. (2010). Genetic polymorphism of the cck and cck1 receptor genes does not affect gastric function or reported satiety.Gut59:A138. 10.1136/gut.2009.209023y

  • 12

    KumarK. G.ByerleyL. O.VolaufovaJ.DruckerD. J.ChurchillG. A.LiR.et al (2008). Genetic variation in Glp1r expression influences the rate of gastric emptying in mice.Am. J. Physiol. Regul. Integr. Comp. Physiol.294R362R371. 10.1152/ajpregu.00640.2007

  • 13

    Lankat-ButtgereitB.GokeB. (1997). Cloning and characterization of the 5′ flanking sequences (promoter region) of the human GLP-1 receptor gene.Peptides18617624. 10.1016/S0196-9781(97)00001-6

  • 14

    LeeY.GamazonE. R.RebmanE.LeeY.LeeS.DolanM. E.et al (2012). Variants affecting exon skipping contribute to complex traits.PLoS Genet.8:e1002998. 10.1371/journal.pgen.1002998

  • 15

    LoosR. J. F.LindgrenC. M.LiS.WheelerE.ZhaoJ. H.PropenkoI. (2008). Common variants near MC4R are associated with fat mass, weight and risk of obesity.Nat. Genet.40768775. 10.1038/ng.140

  • 16

    MayoK. E.MillerL. J.BatailleD.DalleS.GokeB.ThorensB.et al (2003). International union of pharmacology. XXXV. The glucagon receptor family.Pharmacol. Rev.55167194. 10.1124/pr.55.1.6

  • 17

    MuraniE.PonsuksiliS.SeyfertH.-M.ShiX.WimmersK. (2009). Dual effect of a single nucleotide polymorphism in the first intron of the porcine secreted phosphoprotein 1 gene: allele-specific binding of C/EBP beta and activation of aberrant splicing.BMC Mol. Biol.10:96. 10.1186/1471-2199-10-96

  • 18

    SathananthanA.ManC. D.MichelettoF.ZinsmeisterA. R.CamilleriA.GieslerP. D.et al (2010). Common genetic variation in GLP1R and insulin secretion in response to exogenous GLP-1 in nondiabetic subjects a pilot study.Diabetes Care3320742076. 10.2337/dc10-0200

  • 19

    StoffelM.EspinosaR.LebeauM. M.BellG. I. (1993). Human glucagon-like peptide-1 receptor gene - localization to chromosome band-6p21 by fluorescence in-situ hybridization and linkage of a highly polymorphic simple tandem repeat dna polymorphism to other markers on chromosome-6.Diabetes Metab. Res. Rev.4212151218.

  • 20

    TanizawaY.RiggsA. C.ElbeinS. C.WhelanA.Donis-KellerH.PermuttM. A. (1994). Human glucagon-like peptide-1 receptor gene in niddm - identification and use of simple sequence repeat polymorphisms in genetic-analysis.Diabetes Metab. Res. Rev.43752757.

  • 21

    TokuyamaY.MatsuiK.EgashiraT.NozakiO.IshizukaT.KanatsukaA. (2004). Five missense mutations in glucagon-like peptide 1 receptor gene in Japanese population.Diabetes Res. Clin. Prac.666369. 10.1016/j.diabres.2004.02.004

  • 22

    VellaA.SathananthanA.MichelettoF.Dalla ManC.ToffoloG.CobelliC.et al (2009). A pilot study examining the effect of GLP-1R polymorphisms in insulin secretion in response to glucose and GLP-1 infusion in non-diabetic subjects.Diabetologia52:S58.

  • 23

    VogelC. I. G.BoesT.ReinehrT.RothC. L.ScheragS.ScheragA.et al (2011). Common variants near MC4R: exploring gender effects in overweight and obese children and adolescents participating in a lifestyle intervention.Obes. Facts46775. 10.1159/000324557

  • 24

    WettergrenA.SchjoldagerB.MortensenP. E.MyhreJ.ChristiansenJ.HolstJ. J. (1993). Truncated glp-1 (proglucagon 78-107-amide) inhibits gastric and pancreatic functions in man.Dig. Dis. Sci.38665673. 10.1007/BF01316798

  • 25

    WilmenA.WalkenbachA.FullerP.Lankat-ButtgereitB.GokeR.GokeB. (1998). The genomic organization of the human GLP-1 receptor gene.Exp. Clin. Endocrinol. Diabetes106299302. 10.1055/s-0029-1211989

  • 26

    WishartJ. M.HorowitzM.MorrisH. A.JonesK. L.NauckM. A. (1998). Relation between gastric emptying of glucose and plasma concentrations of glucagon-like peptide-1.Peptides1910491053. 10.1016/S0196-9781(98)00052-7

  • 27

    XiaQ. H.GrantS. F. A. (2013). The genetics of human obesity.Ann. N. Y. Acad. Sci.1281178190. 10.1111/nyas.12020

  • 28

    YauA. M. W.McLaughlinJ.MaughanR. J.GilmoreW.EvansG. H. (2014). Short-term dietary supplementation with fructose accelerates gastric emptying of a fructose but not a glucose solution.Nutrition3013441348. 10.1016/j.nut.2014.03.023

  • 29

    YauA. M. W.McLaughlinJ.MaughanR. J.GilmoreW.EvansG. H. (2017). The effect of short-term dietary fructose supplementation on gastric emptying rate and gastrointestinal hormone responses in healthy men.Nutrients9:258. 10.3390/nu9030258

Summary

Keywords

GLP-1 receptor, single nucleotide polymorphisms, gastric emptying rate, GLP-1, gastrointestinal motility

Citation

Yau AMW, McLaughlin J, Maughan RJ, Gilmore W, Ashworth JJ and Evans GH (2018) A Pilot Study Investigating the Influence of Glucagon-Like Peptide-1 Receptor Single Nucleotide Polymorphisms on Gastric Emptying Rate in Caucasian Men. Front. Physiol. 9:1331. doi: 10.3389/fphys.2018.01331

Received

13 March 2018

Accepted

04 September 2018

Published

25 September 2018

Volume

9 - 2018

Edited by

Ravinder Abrol, California State University, Northridge, United States

Reviewed by

Manlio Vinciguerra, International Clinical Research Center (FNUSA-ICRC), Czechia; Matthias J. Bahr, Sana Kliniken Lübeck, Germany

Updates

Copyright

*Correspondence: Gethin H. Evans,

This article was submitted to Gastrointestinal Sciences, a section of the journal Frontiers in Physiology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics