Research article | Open | Open Peer Review | Published:
DNA methylation of the glucagon-like peptide 1 receptor (GLP1R) in human pancreatic islets
BMC Medical Geneticsvolume 14, Article number: 76 (2013)
Insulin secretion is enhanced upon the binding of Glucagon-like peptide-1 (GLP-1) to its receptor (GLP1R) in pancreatic β cells. Although a reduced expression of GLP1R in pancreatic islets from type 2 diabetic patients and hyperglycaemic rats has been established, it is still unknown if this is caused by differential DNA methylation of GLP1R in pancreatic islets of type 2 diabetic patients.
In this study, DNA methylation levels of 12 CpG sites close to the transcription start site of GLP1R were analysed in pancreatic islets from 55 non-diabetic and 10 type 2 diabetic human donors as well as in β and α cells isolated from human pancreatic islets. DNA methylation of GLP1R was related to GLP1R expression, HbA1c levels and BMI. Moreover, mRNA expression of MECP2, DNMT1, DNMT3A and DNMT3B was analysed in pancreatic islets of the non-diabetic and type 2 diabetic donors.
One CpG unit, at position +199 and +205 bp from the transcription start site, showed a small increase in DNA methylation in islets from donors with type 2 diabetes compared to non-diabetic donors (0.53%, p=0.02). Furthermore, DNA methylation levels of one CpG site located 376 bp upstream of the transcription start site of GLP1R correlated negatively with GLP1R expression (rho=−0.34, p=0.008) but positively with BMI and HbA1c (rho=0.30, p=0.02 and rho=0.30, p=0.03, respectively). This specific CpG site is located in an area with known SP1 and SP3 transcription factor binding sites. Moreover, when we compared the DNA methylation of the GLP1R promoter in isolated human β and α cells, we found that it was higher in α- compared with β-cells (p=0.009). Finally, there was a trend towards decreased DNMT3A expression (p=0.056) in type 2 diabetic compared with non-diabetic islets.
Together, our study shows that while BMI and HbA1c are positively associated with DNA methylation levels of GLP1R, its expression is negatively associated with DNA methylation of GLP1R in human pancreatic islets.
Glucagon-like peptide-1 (GLP-1) is an incretin hormone that is secreted by gastrointestinal L-cells in response to meal intake. The peptide is produced by the posttranslational modification of proglucagon . The glucagon-like peptide-1 receptor (GLP1R) is a G protein-coupled receptor, which is expressed in the pancreas, lungs, heart, kidney, stomach, and brain [1–4]. When GLP-1 binds to its receptor in pancreatic β cells, an intracellular signalling cascade is initiated, resulting in the activation of adenylate cyclase and the formation of cAMP. The increase in cAMP enhances the secretion of insulin, providing a mechanism for GLP-1 to regulate insulin secretion in humans. Furthermore, while several studies have shown that GLP1R is expressed in pancreatic β cells, its expression is low or absent in pancreatic α cells [5, 6].
Type 2 diabetes is a multifactorial polygenic disease characterised by chronic hyperglycaemia due to impaired insulin secretion and action. Dissecting the mechanisms that contribute to insufficient insulin secretion in type 2 diabetes patients is an important goal of understanding the disease. Disease susceptibility is affected by genetic and non-genetic factors and a combination thereof. However, epigenetic factors, including DNA methylation and histone modifications, also participate in type 2 diabetes . Indeed, our group has previously demonstrated that DNA methylation correlates negatively with type 2 diabetes candidate gene expression in human pancreatic islets and skeletal muscle [7–12]. Increased DNA methylation and decreased expression of PPARGC1A, INS and PDX1 in pancreatic islets of type 2 diabetic patients is further associated with decreased insulin secretion [8, 11, 13]. However, knowledge about the role of epigenetic mechanisms in the growing incidence of type 2 diabetes is still limited and additional studies analysing epigenetics in humans are hence needed.
Interestingly, GLP1R expression is decreased in pancreatic islets from patients with type 2 diabetes and hyperglycaemic rats [14–16]. Although the GLP1R promoter is GC rich and cytosine-residues in CG-dinucleotides are targets for DNA methylation, no previous study has analysed DNA methylation of the GLP1R promoter in human pancreatic islets. The aim of this study was therefore to analyse the levels of DNA methylation of the GLP1R promoter in human pancreatic islets from 55 non-diabetic organ donors and 10 donors with type 2 diabetes. DNA methylation of GLP1R was further correlated to gene expression, HbA1c levels and BMI. We also tested if genes coding for proteins involved in epigenetic processes are differentially expressed in pancreatic islets from patients with type 2 diabetes compared with non-diabetic donors.
Pancreatic islets from 55 non-diabetic and 10 type 2 diabetic deceased organ donors were obtained from the Human Tissue Laboratory at Lund University Diabetes Centre (Table 1). Islets were prepared by collagenase digestion and density gradient purification. After isolation, islets were cultured free floating in CMRL 1066 culture medium (ICN Biomedicals, Costa Mesa, CA, USA) supplemented with 10 mmol/l HEPES, 2 mmol/l l-glutamine, 50 μg/ml gentamicin, 0.25 μg/ml Fungizone (GIBCO, BRL, Gaithersburg, MD, USA), 20 μg/ml ciprofloxacin (Bayer Healthcare, Leverkusen, Germany), and 10 mmol/l nicotinamide at 37°C (5% CO2) prior to RNA and DNA preparation. Glucose-stimulated insulin release from the human islets was measured in duplicate during dynamic glucose perifusion (Brandel, London, UK) in order to calculate the stimulation index (SI), which was defined as the ratio between the areas under the curves that were calculated for the low (1.67 mM) and high (16.7 mM) glucose concentrations as previously described . The donor before death or her/his relatives upon admission to Intensive Care Unit (ICU) had given their consent to donate organs and the local ethics committees approved the protocols.
Total RNA was isolated with the AllPrep DNA/RNA Mini Kit (Qiagen GmbH, Hilden, Germany). RNA quality and concentration was measured using an Agilent 2100 bioanalyzer (Agilent Technologies, Inc., Santa Clara, CA, USA) and Nanodrop ND-1000 equipment (NanoDrop Technologies, Wilmington, DE), respectively. Gene expression was analysed using the Human Gene 1.0 ST Array (Affymetrix, Santa Clara, CA, USA) analysis following the Affymetrix standard protocol. The array data was summarised and normalised with Robust Multi-array Analysis (RMA) method using the software “Expression Console” (Affymetrix).
β and α cell purification
β and α cells were purified from pancreatic islets of three human donors (aged 54, 55 and 74 years old, with BMI 21.5-23.1 kg/m2), different from the donors described in Table 1, using a method previously described [18, 19]. In short, dissociation of islet cells was achieved by incubation with constant agitation for 3 minutes at 37°C in 0.05 % trypsin-EDTA (Life Technologies Ltd, Paisley, UK) supplemented with 3 mg/ml DNAse I (Roche, Basel, Switzerland) followed by vigorous pipetting. Labelling and FACS sorting of the β- and α-cell fractions was performed as previously described . Sorted α and β-cells were applied to microscope slides and co-immunostained for insulin and glucagon in order to detect the amount of α-cells in the β-cell fraction, and vice versa. Using this method, a β-cell purity of 89 ± 9 (mean ± SD) was achieved .
500 ng of genomic DNA was bisulfite treated using the EZ DNA Methylation kit (Zymo Research, Orange, CA, USA). DNA methylation analysis was performed with EpiTYPER using Sequenom MassARRAY system (Sequenom, Inc., San Diego, CA, USA) as previously described . Two EpiTYPER assays were designed using the online EpiDesigner tool , covering a total of 18 CpG sites in the region upstream or downstream of the transcription start site of the GLP1R gene. Primer information is given in Additional file 1: Table S1. Due to either high or low mass of the cleavage product, no data was generated for 6 CpG sites. Because of the base specific cleavage of the EpiTYPER method, two CpG sites positioned downstream of the transcription start site of the GLP1R gene were analysed as a unit and therefore called CpG site +199/+205.
Statistical analyses were performed using PASW Statistics 18 for Windows (SPSS, Chicago, IL, USA). Non-parametric two samples test, Mann–Whitney U test, was performed to analyse differences between type 2 diabetes and non-diabetic donors. Correlations were analysed using the non-parametric Spearman correlation using all individuals in the study. Paired samples t-test was used to analyse the difference in methylation between α and β cells. The p-values presented in this study have not been corrected for multiple testing. All data is presented as mean ± sd and data in figures as mean ± SEM.
Pancreatic islets from 55 non-diabetic and 10 type 2 diabetic human donors were analysed in this study. Donor characteristics are described in Table 1. Pancreatic islets from type 2 diabetic organ donors showed a decrease in glucose-stimulated insulin secretion compared with islets from non-diabetic donors. In agreement with previous studies [14, 16], we found that GLP1R expression was reduced in pancreatic islets from donors diagnosed with type 2 diabetes compared with non-diabetic donors (type 2 diabetic 213±76.6 vs. non-diabetic 390.4±170.2, p=0.0006). Furthermore, GLP1R mRNA expression in the human islets correlated positively with the stimulation index (SI) of glucose-stimulated insulin secretion (Table 1) (rho=0.33, p=0.015). We next quantified DNA methylation levels of 12 CpG sites of the GLP1R gene, including 5 CpG sites upstream of and 7 CpG sites downstream of the transcription start site (Figure 1A). Two of the studied CpG sites, +199 and +205 bp from the transcription start site, were analysed as a CpG unit, due to the sequence characteristic of GLP1R. The majority of the analysed CpG sites displayed a degree of DNA methylation below 10% (Additional file 2: Table S2). The CpG unit including CpG sites +199 and +205 showed increased DNA methylation in islets from donors with type 2 diabetes compared with non-diabetic donors (Figure 1B). However, the observed difference was small (0.53%, p=0.02) and did not persist after correction for multiple testing. Because previous studies have shown that increased DNA methylation may be associated with decreased gene expression , we next examined if DNA methylation of GLP1R correlates with its gene expression in the human islets. The CpG site located 376 bp upstream of the transcription start site showed negative correlation with GLP1R expression (rho=−0.34, p=0.008) (Figure 2A). In addition, the same CpG site at position −376 correlated positively with both BMI (rho=0.30, p=0.02) and HbA1c levels (rho=0.30, p=0.03) (Figure 2B and C). There were no significant correlations between the degree of DNA methylation of the CpG unit including CpG sites +199 and +205 and GLP1R mRNA expression, BMI or HbA1c levels (p>0.05). Neither did DNA methylation of the other analysed CpG sites correlate negatively with GLP1R mRNA expression or positively with BMI and HbA1c levels (p>0.05). Transcription factor binding sites around position −376 were queried in silico with a web-based tool, TFsearch . CpG site −376 is located in a putative SP1 transcription factor binding site as marked in Figure 1A.
Moreover, since it has been shown that GLP1R mainly is expressed in β cells of pancreatic islets , we tested if DNA methylation of the GLP1R promoter differed in α and β cells isolated from three human pancreatic islet donors. We found a significant increase in DNA methylation of CpG site −376 in α cells compared with β cells (α cells 14±5 vs. β cells 8±4, p=0.009) (Figure 3). However, there were no differences in DNA methylation between α and β cells for the other analysed CpG sites (Additional file 3: Table S3).
Finally, we tested if the methyl binding protein MECP2, which is known to control gene expression through the interaction with transcriptional repressors , as well as three DNA methyl transferases; DNMT1, DNMT3A and DNMT3B, show differential expression in pancreatic islets from patients with type 2 diabetes. There were no differences in MECP2, DNMT1 and DNMT3B expression between diabetic and non-diabetic donors (Figure 4A, B, D). However, there was a trend towards decreased expression of DNMT3A (p=0.056) in diabetic versus non-diabetic islets (Figure 4C).
A number of studies have shown differential DNA methylation of a number of type 2 diabetic candidate genes in human pancreatic islets pointing to a potential role of DNA methylation in the pathogenesis of the disease [8, 11–13, 23, 24]. The present study suggests that DNA methylation of a CpG site located 376 bp upstream from the transcription start site in the promoter of GLP1R may affect the expression of this gene. We found an inverse correlation between DNA methylation of this CpG site and GLP1R gene expression, suggesting that methylation of this specific CpG site could have a negative effect on GLP1R gene expression. Wildhage et al. have reported that the GLP1R promoter contains at least three SP1 binding sites which SP1 can bind to. The related transcription factor SP3 can bind to only one of the sites . CpG −376 is located in the binding site to which both SP1 and SP3 can bind. While SP1 in most cases serves as a transcriptional activator, SP3 has been reported to act as both an activator and as a repressor, depending on cell and tissue type. The role of SP3 in regulating gene expression depends on several different factors, including relative levels of SP1 and SP3, the number of SP1 binding sites in a promoter, cell and/or tissue type, interaction with other proteins, and chemical modification of the transcription factor.  Increased DNA methylation of promoter regions is associated with transcriptional silencing , either through preventing the binding of specific transcription factors or through recruitment of methyl CpG binding proteins, e.g. MeCP2, that promote recruitment of histone deacetyltransferases and/or co-repressors. DNA methylation of CpG −376 could potentially repress the binding of SP1 to the GLP1R promoter and hence result in transcriptional silencing. However, additional functional studies are required to examine the specific role of GLP1R DNA methylation on gene transcription.
It has previously been established that both protein and mRNA expression of GLP1R is reduced in islets from patients with type 2 diabetes compared with non-diabetic controls [14, 16] as well as in rodent islets exposed to hyperglycaemia . Our study identified one CpG site showing a small increase in GLP1R DNA methylation in islets from type 2 diabetes patients. However, the GLP1R promoter is highly GC rich making assay design for DNA methylation analysis in large parts of the GLP1R promoter inaccessible for analysis due to the requirement for CpG-free primers when using standard techniques e.g. EpiTYPER or Pyrosequencing. Hence, we cannot exclude that other areas of the GLP1R may also show differential DNA methylation in type 2 diabetes islets.
Hyperglycaemia, a high fat diet and obesity have previously been associated with differential DNA methylation in target tissues for type 2 diabetes such as pancreatic islets, skeletal muscle and liver [11, 13, 27–30]. For instance, HbA1c has been shown to correlate positively with DNA methylation of INS and PDX-1 in human pancreatic islets and hyperglycaemia was associated with increased DNA methylation of the same genes in clonal β cells cultured in vitro. Based on the positive correlations between both HbA1c and BMI and DNA methylation of CpG site −376, the present study suggests that hyperglycaemia and/or obesity may affect DNA methylation of GLP1R in human islets. Interestingly, hyperglycaemia has previously been shown to increase the expression of a methyl transferase, Dnmt1, in clonal β cells . Elevated levels of this methyl transferase may be a mechanism behind glucose induced DNA methylation. However, we did not find any significant mRNA expression differences of DNMT1, DNMT3A and DNMT3B in pancreatic islets from patients with type 2 diabetes. There was a trend towards a decreased expression of DNMT3A in type 2 diabetics compared with non-diabetics, but the difference was small and not significant (200.8±13.8 and 212.7±19.1 respectively, p=0.056). It is further possible that the levels of methyl donors affect the degree of DNA methylation in diabetic patients. Indeed, the level of S-adenosylmethionine, a methyl group donor, has been reported to be decreased in the erythrocytes of patients with type 2 diabetes .
Overall, our data suggest that DNA methylation influences gene expression of GLP1R in human pancreatic islets. Decreased binding of the transcription factors SP1 and SP3 due to increased methylation may be involved in reducing GLP1R expression. Furthermore, reduced islet GLP1R levels associate with lower insulin secretion, which is seen in patients with type 2 diabetes. The results of this study again indicate that epigenetic mechanisms may contribute to type 2 diabetes, however additional studies are needed to fully understand the mechanisms involved in this regulation.
Holst JJ: The physiology of glucagon-like peptide 1. Physiol Rev. 2007, 87 (4): 1409-1439.
Wei Y, Mojsov S: Tissue-specific expression of the human receptor for glucagon-like peptide-I: brain, heart and pancreatic forms have the same deduced amino acid sequences. FEBS Lett. 1995, 358 (3): 219-224.
Mayo KE, Miller LJ, Bataille D, Dalle S, Goke B, Thorens B, Drucker DJ: International Union of Pharmacology. XXXV. The glucagon receptor family. Pharmacol Rev. 2003, 55 (1): 167-194.
Doyle ME, Egan JM: Mechanisms of action of glucagon-like peptide 1 in the pancreas. Pharmacol Ther. 2007, 113 (3): 546-593.
De Marinis YZ, Salehi A, Ward CE, Zhang Q, Abdulkader F, Bengtsson M, Braha O, Braun M, Ramracheya R, Amisten S, et al: GLP-1 inhibits and adrenaline stimulates glucagon release by differential modulation of N- and L-type Ca2+ channel-dependent exocytosis. Cell Metab. 2010, 11 (6): 543-553.
Tornehave D, Kristensen P, Romer J, Knudsen LB, Heller RS: Expression of the GLP-1 receptor in mouse, rat, and human pancreas. J Histochem Cytochem. 2008, 56 (9): 841-851.
Ling C, Groop L: Epigenetics: a molecular link between environmental factors and type 2 diabetes. Diabetes. 2009, 58 (12): 2718-2725.
Ling C, Del Guerra S, Lupi R, Ronn T, Granhall C, Luthman H, Masiello P, Marchetti P, Groop L, Del Prato S: Epigenetic regulation of PPARGC1A in human type 2 diabetic islets and effect on insulin secretion. Diabetologia. 2008, 51 (4): 615-622.
Ling C, Poulsen P, Simonsson S, Ronn T, Holmkvist J, Almgren P, Hagert P, Nilsson E, Mabey AG, Nilsson P, et al: Genetic and epigenetic factors are associated with expression of respiratory chain component NDUFB6 in human skeletal muscle. J Clin Invest. 2007, 117 (11): 3427-3435.
Ronn T, Poulsen P, Hansson O, Holmkvist J, Almgren P, Nilsson P, Tuomi T, Isomaa B, Groop L, Vaag A, et al: Age influences DNA methylation and gene expression of COX7A1 in human skeletal muscle. Diabetologia. 2008, 51 (7): 1159-1168.
Yang BT, Dayeh TA, Kirkpatrick CL, Taneera J, Kumar R, Groop L, Wollheim CB, Nitert MD, Ling C: Insulin promoter DNA methylation correlates negatively with insulin gene expression and positively with HbA(1c) levels in human pancreatic islets. Diabetologia. 2011, 54 (2): 360-367.
Dayeh TA, Olsson AH, Volkov P, Almgren P, Ronn T, Ling C: Identification of CpG-SNPs associated with type 2 diabetes and differential DNA methylation in human pancreatic islets. Diabetologia. 2013, 56 (5): 1036-1046.
Yang BT, Dayeh TA, Volkov PA, Kirkpatrick CL, Malmgren S, Jing X, Renstrom E, Wollheim CB, Nitert MD, Ling C: Increased DNA methylation and decreased expression of PDX-1 in pancreatic islets from patients with type 2 diabetes. Mol Endocrinol. 2012, 26 (7): 1203-1212.
Shu L, Matveyenko AV, Kerr-Conte J, Cho JH, McIntosh CH, Maedler K: Decreased TCF7L2 protein levels in type 2 diabetes mellitus correlate with downregulation of GIP- and GLP-1 receptors and impaired beta-cell function. Hum Mol Genet. 2009, 18 (13): 2388-2399.
Xu G, Kaneto H, Laybutt DR, Duvivier-Kali VF, Trivedi N, Suzuma K, King GL, Weir GC, Bonner-Weir S: Downregulation of GLP-1 and GIP receptor expression by hyperglycemia: possible contribution to impaired incretin effects in diabetes. Diabetes. 2007, 56 (6): 1551-1558.
Taneera J, Lang S, Sharma A, Fadista J, Zhou Y, Ahlqvist E, Jonsson A, Lyssenko V, Vikman P, Hansson O, et al: A systems genetics approach identifies genes and pathways for type 2 diabetes in human islets. Cell Metab. 2012, 16 (1): 122-134.
Stahle MU, Brandhorst D, Korsgren O, Knutson F: Pathogen inactivation of human serum facilitates its clinical use for islet cell culture and subsequent transplantation. Cell Transplant. 2011, 20 (5): 775-781.
Parnaud G, Bosco D, Berney T, Pattou F, Kerr-Conte J, Donath MY, Bruun C, Mandrup-Poulsen T, Billestrup N, Halban PA: Proliferation of sorted human and rat beta cells. Diabetologia. 2008, 51 (1): 91-100.
Kirkpatrick CL, Marchetti P, Purrello F, Piro S, Bugliani M, Bosco D, de Koning EJ, Engelse MA, Kerr-Conte J, Pattou F, et al: Type 2 diabetes susceptibility gene expression in normal or diabetic sorted human alpha and beta cells: correlations with age or BMI of islet donors. PLoS One. 2010, 5 (6): e11053-
Bird A: Perceptions of epigenetics. Nature. 2007, 447 (7143): 396-398.
Sandovici I, Smith NH, Nitert MD, Ackers-Johnson M, Uribe-Lewis S, Ito Y, Jones RH, Marquez VE, Cairns W, Tadayyon M, et al: Maternal diet and aging alter the epigenetic control of a promoter-enhancer interaction at the Hnf4a gene in rat pancreatic islets. Proc Natl Acad Sci USA. 2011, 108 (13): 5449-5454.
Volkmar M, Dedeurwaerder S, Cunha DA, Ndlovu MN, Defrance M, Deplus R, Calonne E, Volkmar U, Igoillo-Esteve M, Naamane N, et al: DNA methylation profiling identifies epigenetic dysregulation in pancreatic islets from type 2 diabetic patients. EMBO J. 2012, 31 (6): 1405-1426.
Wildhage I, Trusheim H, Goke B, Lankat-Buttgereit B: Gene expression of the human glucagon-like peptide-1 receptor is regulated by Sp1 and Sp3. Endocrinology. 1999, 140 (2): 624-631.
Li L, He S, Sun JM, Davie JR: Gene regulation by Sp1 and Sp3. Biochem Cell Biol. 2004, 82 (4): 460-471.
Brons C, Jacobsen S, Nilsson E, Ronn T, Jensen CB, Storgaard H, Poulsen P, Groop L, Ling C, Astrup A, et al: Deoxyribonucleic acid methylation and gene expression of PPARGC1A in human muscle is influenced by high-fat overfeeding in a birth-weight-dependent manner. J Clin Endocrinol Metab. 2010, 95 (6): 3048-3056.
Franks PW, Ling C: Epigenetics and obesity: the devil is in the details. BMC Med. 2010, 8: 88-
Jiang M, Zhang Y, Liu M, Lan MS, Fei J, Fan W, Gao X, Lu D: Hypermethylation of hepatic glucokinase and L-type pyruvate kinase promoters in high-fat diet-induced obese rats. Endocrinology. 2011, 152 (4): 1284-1289.
Jiang MH, Fei J, Lan MS, Lu ZP, Liu M, Fan WW, Gao X, Lu DR: Hypermethylation of hepatic Gck promoter in ageing rats contributes to diabetogenic potential. Diabetologia. 2008, 51 (8): 1525-1533.
Poirier LA, Brown AT, Fink LM, Wise CK, Randolph CJ, Delongchamp RR, Fonseca VA: Blood S-adenosylmethionine concentrations and lymphocyte methylenetetrahydrofolate reductase activity in diabetes mellitus and diabetic nephropathy. Metabolism. 2001, 50 (9): 1014-1018.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/14/76/prepub
We thank SCIBLU at Lund University for analysing gene expression. We thank the Nordic Network for Clinical Islet Transplantation (JDRF award 31-2008-413), the tissue isolation teams and Human Tissue Laboratory within EXODIAB/Lund University Diabetes Center. We are grateful to Prof. P. Marchetti of the University of Pisa and Dr Bosco of Geneva University Hospital for the provision of human islets for the β- and α-cell sorting and to Dr Jalal Taneera for help with extracting RNA and DNA from human islets. This work was supported by grants from the Swedish Research Council, the Wallenberg Foundation, ALF, Linné grant (B31 5631/2006), the Novo Nordisk Foundation, UMAS Fonder, the Lund University Diabetes Centre (LUDC), Tore Nilsson, Söderberg, Syskonen Svenssons Foundation, Diabetes foundation, Kungliga Fysiografiska Sällskapet, SSMF and Påhlsson foundation.
The authors declare that they have no competing interests.
EH designed and conducted the study, collected, analysed and interpreted data and wrote the manuscript. TD collected and interpreted data and reviewed and edited manuscript. CLK collected data and reviewed and edited manuscript. CBW interpreted data and reviewed and edited the manuscript. MDN interpreted data and reviewed and edited the manuscript. CL designed and conducted the study, interpreted data and reviewed and edited the manuscript. EH and CL are guarantors of this work and, as such, had full access to all of the data in the study and take responsibility for the integrity of the data. All authors read and approved the final manuscript.