Skip to main content

ORIGINAL RESEARCH article

Front. Genet., 13 June 2022
Sec. Applied Genetic Epidemiology

GDF-15 as a Therapeutic Target of Diabetic Complications Increases the Risk of Gallstone Disease: Mendelian Randomization and Polygenic Risk Score Analysis

  • 1Department of Big Data in Health Science School of Public Health, Center of Clinical Big Data and Analytics of The Second Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, China
  • 2Unit of Cardiovascular and Nutritional Epidemiology, Institute of Environmental Medicine, Karolinska Institutet, Stockholm, Sweden
  • 3Unit of Medical Epidemiology, Department of Surgical Sciences, Uppsala University, Uppsala, Sweden
  • 4Centre for Global Health Research, Usher Institute, University of Edinburgh, Edinburgh, United Kingdom

Growth differentiation factor 15 (GDF-15) levels have been revealed as a robust biomarker for metformin use. We conducted Mendelian randomization (MR) analysis to explore the association between GDF-15 and gallstone disease to inform potential therapeutic effects targeting GDF-15. Four genetic variants associated with GDF-15 levels at p < 5 × 10–8 were selected as instrumental variables from a genome-wide association meta-analysis including 21,758 individuals. Two-sample MR analysis was conducted using summary-level data from UK Biobank (10,520 gallstone cases and 350,674 controls) and FinnGen consortium (19,023 gallstone cases and 195,144 controls). Polygenic risk score analysis using individual-level data in UK biobank was performed to complement the MR findings by examining the non-linearity of the association. Diabetic complications were taken as positive controls to validate the therapeutic effect of targeting GDF-15. Linear and nonlinear associations between genetically predicted GDF-15 levels and gallstones were estimated with stratification by the diabetic status. In the two-sample MR analysis, the odds ratio (OR) of gallstones was 1.09 (95% confidence interval (CI), 1.03–1.15; p = 0.001) for one standard deviation increase in genetically predicted GDF-15 levels in the meta-analysis of two datasets. Polygenic risk score analysis found this association to be U-shaped (p = 0.037). The observed association was predominantly seen in nondiabetic population (OR = 1.11, 95% CI: 1.01–1.21; p = 0.003). An inverse association between genetically predicted GDF-15 levels and diabetic complications (OR = 0.77, 95% CI: 0.62–0.96; p = 0.023) was observed, validating the potential therapeutic effects of targeting GDF-15 levels. This MR study indicates that the increased risk of gallstone disease should be taken into account when considering GDF-15 as a therapeutic target for diabetic complications.

Introduction

Growth differentiation factor 15 (GDF-15) is one of the members of transforming growth factor beta (TGF-β), specifically produced by macrophages and epithelial cells, as well as by adipocytes, with the basic effects on inhibition of macrophage activation, anti-inflammation, and maintenance of homeostasis (Bootcov et al., 1997; Unsicker et al., 2013).

Experimental data on mice found that GDF-15 controlled appetite (Johnen et al., 2007; Macia et al., 2012), reduced body weight and fat mass (Johnen et al., 2007; Macia et al., 2012; Chrysovergis et al., 2014), increased metabolism rate (Chrysovergis et al., 2014), and improved glycemic profiles (Li et al., 2013). A high expression of GDF-15 has been proposed to be a predicting factor (Bao et al., 2019) and a potential treatment target for diabetes mellitus along with complications (Nakayasu et al., 2020; Xiao et al., 2022).

Evidence from some population-based studies has also been provided that increased GDF-15 might be a cause mediator or predicting biomarker in health-adverse outcomes, such as insulin resistance (Vila et al., 2011; Kempf et al., 2012) and cardiovascular disease (Brown et al., 2002; Wollert et al., 2017). Recently, GDF-15 levels have been revealed as a robust biomarker for metformin use (Gerstein et al., 2017; Day et al., 2019), which is well-recognized as the first-line medication for diabetes (Suissa and Azoulay, 2012; Griffin et al., 2017). The medication can lower diabetic complications and has been suggested to have potential benefits in cardiovascular disease and cancer prevention in diabetic individuals, although corresponding evidence remains inconclusive. However, its secondary effects on gastrointestinal complaints, including gallstones, are scarcely evaluated due to lack of high-quality data (Schlender et al., 2017).

Gallstone disease is defined by the occurrence of symptoms or complications caused by gallstones in the gallbladder and/or the bile ducts (Lammert et al., 2016), which has been proposed to be a potential secondary disease of diabetes (Lieber, 1952; Everhart and Ruhl, 2009), obesity (Yuan et al., 2021), and insulin resistance (Nervi et al., 2006), possibly through increased hepatic cholesterol uptake and synthesis, gallbladder hypomotility, etc. Of note, some studies have observed a higher expression of GDF-15 in patients with benign biliary diseases like gallstone disease, cholecystitis, etc. than in healthy controls (Li et al., 2020; Ozkan et al., 2011), providing a new perspective for exploration of the novel role of GDF-15 in the progress of gallstone disease. In spite of limited direct evidence on GDF-15 in relation to gallstones, the associations of GDF-15 or gallstones, respectively, with diabetes indirectly support and enrich the hypothesis.

Employing genetic variants as instrumental variables for exposure (e.g., GDF-15 levels) and Mendelian randomization (MR) analysis can strengthen causal inference by minimizing unobserved confounding and diminishing reverse causality (Burgess and Thompson, 2015). Therefore, we conducted a two-sample Mendelian randomization (MR) study and polygenic risk score analysis, here, to determine the association between genetically predicted GDF-15 levels, a biomarker of metformin use, and the risk of gallstones.

Methods

Study Design and Data Sources

Figure 1 shows the overview of the study design. We first performed a two-sample MR study based on the UK Biobank study and the FinnGen study to determine the causal association between genetically predicted GDF-15 and gallstone disease. We then conducted a polygenic risk score analysis based on individual-level data from the UK Biobank study to further explore the diabetes-specific association and the nonlinearity of the association. UK Biobank received ethical approval from the North West Multi-centre Research Ethics Committee, the National Information Governance Board for Health and Social Care in England and Wales, and the Community Health Index Advisory Group in Scotland. All participants provided written informed consent. The ethical permit for MR analyses based on summary-level data was unnecessary.

FIGURE 1
www.frontiersin.org

FIGURE 1. Schematic representation of the study design. MR, Mendelian randomization; IVW, inverse variance weighted; GRS, weighted genetic risk score; LD, linkage disequilibrium; GWAS, genome-wide association study.

Genetic Instrument Selection

Genetic variants associated with GDF-15 levels at the genome-wide significance threshold (p < 5 × 10–8) were identified from a meta-analysis of 13 genome-wide association studies including up to 21,758 individuals of the European ancestry (Folkersen et al., 2020). Linkage disequilibrium was calculated based on the 1,000 genomes European reference panel (Clarke et al., 2012). Four single nucleotide polymorphisms (rs2517481, rs1227734, rs60164552, and rs112253475) without linkage disequilibrium (r2 ≤ 0.01 and clump window > 10,000 kb) were selected as instrumental variables for GDF-15 (Table 1). These SNPs explain approximately 8.4% of phenotypic variance. The association test was adjusted for population structure and study-specific parameters (Folkersen et al., 2020). Standard sample-level quality control (QC) exclusions were based on call rate filters, sex mismatch, population outliers, heterozygosity, and cryptic relatedness. (Folkersen et al., 2020).

TABLE 1
www.frontiersin.org

TABLE 1. Genetic instruments for growth differentiation factor 15 levels and their associations with gallstones.

Two-Sample MR Analysis

Two-sample MR analyses were performed using summary-level data for the association of GDF-15–associated SNPs with gallstone disease from the UK Biobank study (10,520 cases and 350,674 controls) and the FinnGen study (19,023 cases and 195,144 controls). We used the Neale Lab second wave results of the genome-wide analysis in UK Biobank after the removal of individuals of the non-European ancestry, closely related individuals, individuals with sex chromosome aneuploidies, and individuals who had withdrawn consent. The association test was adjusted for age, sex, and up to 20 genetic principal components. The R5 release of FinnGen consortium data was used, where individuals with ambiguous gender, high genotype missingness, excess heterozygosity, and non-Finnish ancestry had been excluded. The association test in FinnGen was adjusted for age, sex, 10 genetic principal components, and genotyping batch.

Polygenic Risk Score Analysis

Polygenic risk score analysis was conducted in the UK Biobank study including 16,463 gallstone patients and 318,308 controls after exclusion of individuals with high missingness or heterozygosity, outlying short runs of homozygosity, sex mismatch, and non-British ancestry (to minimize population bias). Incident and prevalent gallstone cases were defined by corresponding codes in the International Classification of Diseases (ICD) versions 9 and 10 with information from national medical records (including inpatient hospital episode records, primary general practitioner data, cancer registry, and death registry). Diabetic complications, including type 2 diabetes with ketoacidosis, renal manifestations, ophthalmic manifestations, and neurological manifestations (1,935 cases and 314,181 controls) were taken as positive controls to validate the effect of GDF-15 as a therapeutic target for diabetic complications. Individual-level UK Biobank data was accessed up to 6th March 2021 under the approved project application (ID: 66345). Polygenic risk scores for GDF-15 levels were constructed by summing the number of risk alleles (each SNP was recoded as 0, 1, or 2) carried by each participant and weighted by the effect size estimates (i.e., beta coefficient) of the four genetic instrument variables (Table 1) selected for GDF-15.

Statistical Analysis

For two-sample MR analysis, the inverse-variance weighted method with random effects was used as the main method and supplemented by four sensitivity analyses, including the weighted median (Bowden et al., 2016), MR-Egger regression (Bowden et al., 2015), MR-PRESSO (Verbanck et al., 2018), and leave-one-out analysis. The weighted median method can provide consistent causal estimates if more than 50% of the weight comes from valid instrumental variables (Bowden et al., 2016). MR-Egger can generate estimates after correcting for horizontal pleiotropy; however, this method compromises statistical power (Bowden et al., 2015). MR-PRESSO can detect outlying instrumental variables and provide causal estimates after the removal of these outliers (Verbanck et al., 2018). The leave-one-out analysis can detect whether the association is driven by a certain SNP. In addition, we made a scatterplot to show the association for each SNP. Cochrane’s Q value was used to assess the heterogeneity among estimates of genetic instruments, and the p-value for the intercept in MR-Egger was used to detect horizontal pleiotropy (Bowden et al., 2015). Possible pleiotropic associations with GDF-15-associated SNPs were searched in PhenoScanner V2, a database of human genotype–phenotype associations (Kamat et al., 2019).

For polygenic risk score analysis, logistic regression was first used to explore the linear association between the polygenic score of GDF-15 and gallstone risk. We have three levels of models: 1) the univariate regression model without adjustment for any covariates; 2) the multivariable regression model with adjustment for basic covariates, including age, sex, body mass index (BMI), and the first 20 genetic principal components; and 3) the multivariable regression model with adjustment for a comprehensive list of covariates, including age, sex, BMI and the first 20 genetic principal components, deprivation index, smoking, drinking, physical activity, blood low-density lipoprotein cholesterol (LDL-C) levels, blood glucose levels, blood HbA1c levels, and metformin use. Considering confounding effects from type 2 diabetes (Yuan et al., 2021) and most diabetic patients had taken metformin, in whom the genetically predicted GDF-15 cannot proxy the true level of exposure, we stratified the association by diabetic status. We also examined the nonlinearity of the association by using logistic regression with restricted cubic splines. Spline knot from nonlinearity regression was used as a cutoff to stratify the study population into genetically predicted low and high GDF-15 groups, and the linear association between the polygenic score of GDF-15 and gallstone risk was assessed within each group separately. The nonlinearity regression of polygenic score analysis was conducted using the rms package, and MR analyses were performed using the TwoSampleMR and MR-PRESSO packages (Hemani et al., 2018; Verbanck et al., 2018) in R Software 4.0.2.

Results

Two-Sample MR Analysis

Higher genetically predicted GDF-15 levels were associated with an increased risk of gallstone disease (Figure 2). For one standard deviation increase in genetically predicted GDF-15 levels, the odds ratio (OR) of gallstones was 1.14 (95% confidence interval [95% CI]: 1.01, 1.29; p = 0.028) in UK Biobank, 1.08 (95% CI: 1.02, 1.14; p = 0.012) in FinnGen, and 1.09 (95% CI: 1.03, 1.15; p = 0.001) in the meta-analysis of the two datasets. Four GDF-15–associated SNPs showed consistent associations with gallstones in FinnGen (Supplementary Figure S1). The association was driven by rs1227734 which explains the majority of explained variance of GDF-15 (Figure 2; Supplementary Figure S2). The association was consistent in sensitivity analyses (Table 2). We detected moderate heterogeneity between estimates of SNPs, no indication of horizontal pleiotropy in MR-Egger (p for intercept >0.05), and no outlier in MR-PRESSO analysis (Table 2). Rs2517481 but no other SNPs showed possible pleiotropic associations with other phenotypes (Supplementary Table S1).

FIGURE 2
www.frontiersin.org

FIGURE 2. Association of genetically predicted growth differentiation factor 15 levels with gallstone disease risk in Mendelian randomization analysis. CI, confidence interval; IVW, inverse variance weighted; OR, odds ratio; SNPs, single nucleotide polymorphisms; UKBB, UK Biobank. Estimates of analyses based on all SNPs were derived from the random-effects inverse-variance weighted method.

TABLE 2
www.frontiersin.org

TABLE 2. Association of genetically predicted growth differentiation factor 15 levels with gallstone disease risk in Mendelian randomization sensitivity analysis.

Polygenic Risk Score Analysis

Characteristics of 339,197 participants in the UK Biobank are shown in Table 3. The population with gallstones was more likely to be females with higher body mass index, diabetes, and metformin use. A positive association between genetically predicted GDF-15 levels and gallstone disease was observed in the polygenic risk score analysis (Table 4). The OR of gallstones was 1.11 (95% CI: 1.03, 1.19; p = 0.008) for one standard deviation increase in genetically predicted GDF-15 levels in a univariate model, and the effect estimate remained significant (OR = 1.10, 95% CI: 1.03, 1.19; p = 0.027) after adjustment for additional full list of covariates. In the stratification analysis by diabetic status, the association between genetically predicted GDF-15 and gallstones was consistently observed in nondiabetic participants (OR = 1.11, 95% CI: 1.01, 1.21; p = 0.003), indicating the observed association was not driven by diabetes (Table 4). In diabetic individuals, the association was insignificant (OR = 1.06, 95% CI: 0.85, 1.32; p = 0.605), and the power calculation indicated insufficient statistical power (13%) to detect this effect among diabetic individuals. Furthermore, we found a nonlinear U-shaped association by a regression model with cubic splines (the turning point of polygenic risk score of -0.06, p for nonlinearity = 0.037) (Figure 3). When restricting to nondiabetic individuals, genetically predicted GDF-15 levels showed an inverse but non-significant association with gallstone risk (OR = 0.82, 95% CI: 0.59, 1.13; p = 0.228) when it is below the turning point, while genetically predicted high levels of GDF-15 (above the turning point) showed a significant positive association with an increased risk of gallstones (OR = 1.28, 95% CI: 1.90, 1.43; p = 0.002).

TABLE 3
www.frontiersin.org

TABLE 3. Characteristics of participants in the UK Biobank.

TABLE 4
www.frontiersin.org

TABLE 4. Association of genetically predicted growth differentiation factor 15 levels with gallstone disease risk in UK Biobank.

FIGURE 3
www.frontiersin.org

FIGURE 3. Nonlinear regression of association between the polygenic risk scores of growth differentiation factor 15 and gallstones in overall participants (A), nondiabetic participants (B) and diabetic participants (C). Associations were adjusted for a comprehensive list of covariates, including age, sex, BMI and the first 20 genetic principal components, deprivation index, smoking, drinking, physical activity, blood LDL-c levels, blood glucose levels, blood HbA1c levels, and metformin intake.

Secondary Analysis of Positive Controls

Higher genetically predicted GDF-15 levels were associated with a decreased risk of diabetic complications, validating the therapeutic effect on diabetic complications targeting GDF-15 levels; the corresponding OR was 0.77 (95% CI: 0.62–0.96; p = 0.023) for one standard deviation increase in genetically predicted GDF-15 levels.

Discussion

In this first study on GDF-15 in relation to gallstones, we found an association between genetically predicted GDF-15 levels and the risk of gallstone disease. This association was confined to nondiabetic individuals and appeared to be U-shaped even though a significant association was observed for high genetically predicted GDF-15 levels but not for low genetically predicted GDF-15 levels. We observed an inverse association between GDF-15 levels and diabetic complications, which strengthened the potential usage of GDF-15 levels as a therapeutic target for diabetic complications. Taken together, our findings implied a possible adverse outcome, gallstones, when considering GDF-15 as a potential therapeutic target for diabetic complications.

Our study found that genetically predicted GDF-15 levels were significantly associated with gallstones, which support our hypothesis. Also, evidence from previous studies has explored the association between GDF-15 and gallstone disease. A case-control study observed a higher level of GDF-15 in patients with benign biliary diseases (acute cholecystitis, acute cholangitis, choledocholithiasis, and cholelithiasis), suggesting GDF-15 to be a potential biomarker in biliary diseases (Ozkan et al., 2011). Likewise, another case-control study found that the levels of GDF-15 were significantly higher in non-cirrhotic primary biliary cholangitis (PBC) patients than in healthy controls, along with a higher expression of total bilirubin (TBIL) and direct bilirubin (DBIL) (Li et al., 2020), the typical markers of cholestasis, showing the notable role of GDF-15 as a risk predictor for cholestasis and cholelithiasis.

There were also some studies indicating opposite findings to ours. An experiment in mice found that metformin that elevates circulating GDF-15 levels (Coll et al., 2020) reduced the risk of gallstones induced by a high-fat diet (Dorvash et al., 2018). In a population-based study involving 36 patients with polycystic ovary syndrome and 20 healthy controls, metformin treatment was associated with an improved metabolic and hormonal imbalance associated with polycystic ovary syndrome and increased gallbladder motility (Isik et al., 2012) that is associated with a lower risk of cholesterol gallstones (Pauletzki et al., 1996). A cohort study in diabetes also observed a reduced risk of gallstones in individuals with >180 cumulative defined daily dose of metformin compared to those with <29 cumulative defined daily dose of metformin (Liao et al., 2017). Even though these findings nominally are conflicting with our findings, several issues in previous studies deserve consideration, including lack of direct GDF-15 measurement, unadjusted important confounders (e.g., smoking status), no examination of the nonlinearity of the association, and certain special features of the diabetic population and patients with polycystic ovary syndrome. In addition, time-related biases, such as immortal time bias, time-window bias, and time-lag bias, have been acknowledged as predominant methodological issues in observational studies on metformin use (Suissa and Azoulay, 2012), which may mislead the association between metformin use and gallstones.

In stratification analysis, the significant association between higher genetically predicted GDF-15 and gallstones was only observed in nondiabetic individuals. The insignificant association in diabetic patients could be due to the lack of study power (13%) and the possibility that genetic polymorphisms may explain a smaller variance of GDF-levels, given that metformin treatment is much more frequent in diabetic patients. In addition, we observed a U-shaped association between GDF-15 and gallstone disease. The nonlinearity of the observed association has important clinical implications on the dose of metformin use which should maintain GDF-15 at a reasonable level to balance the therapeutic effects on diabetic complications and the risk of gallstone disease.

Even though the underlying mechanisms that explain the positive association between GDF-15 and gallstone disease have not been expounded, several possible pathways have been proposed. Increased GDF-15 levels are mostly linked to inflammation (Wischhusen et al., 2020), such as higher levels of interleukin 6 (Breit et al., 2011; Wischhusen et al., 2020; Georgakis et al., 2021), and such inflammation can facilitate cholelithogenesis (Liu et al., 2018). GDF-15 can influence hepatic lipid metabolism and secretion and therefore alters biliary cholesterol secretion and bile concentration per se and via inflammation (Maurer et al., 2009; Luan et al., 2019). In addition, GDF-15 may be a proxy for immobilization (Conte et al., 2020), impaired fasting glucose (Hong et al., 2014), and insulin resistance (Vila et al., 2011; Kempf et al., 2012), which are risk factors for the development of gallstones.

Our study has several strengths, including MR design that strengthened causal inference, consistent associations in two independent populations, the confined study population of the European ancestry that minimized the population structure bias, stratification analysis by the diabetic status that diminished unobserved confounding from diabetes, and examination of nonlinearity. In addition, our study validated the use of GDF-15 as a potential therapeutic target to reduce the risk of diabetic complications which could be regarded as positive controls. Nevertheless, limitations need attention when interpreting our results. We observed pleiotropic associations for rs2517481 but not for other SNPs at the genome significance level. However, this pleiotropy was supposed to generate a limited impact on our findings since this SNP explains a small variance in GDF-15 levels, and the observed association between GDF-15 and gallstones was mainly driven by rs1227734 located in the GDF15 gene region. The population confinement to European individuals may limit the generalizability of our findings to other populations, such as Asians and Africans. More research on GDF-15 in relation to gallstone disease in other populations is needed. In addition, the validity of genetic instruments for GDF-15 could not be examined in outcome data, given no phenotype information was available in UK Biobank and FinnGen. However, we obtained these SNPs from a hither to the largest genome-wide association meta-analysis study on GDF-15 (Folkersen et al., 2020), and the observed association appeared to be steered by the SNP located in the protein-encoding GDF15 gene. Also, the confirmation of the shared causal variants between GDF-15 circulating levels and gallstone disease could not be examined due to insufficient data. Given that our findings are novel, a more functional study and randomized controlled trial (RCT) with causal inference potentials are needed to verify this association, and the corresponding underlying mechanisms also deserve further investigation.

In conclusion, this MR study suggests a detrimental effect of high GDF-15 levels on gallstone formation for the first time. This finding may indicate gallstone disease as a potential adverse outcome, when taking GDF-15 as a therapeutic target for diabetic complications as indicated by the metformin use.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.

Ethics Statement

The studies involving human participants were reviewed and approved by the North West Multi-centre Research Ethics Committee, the National Information Governance Board for Health and Social Care in England and Wales, and the Community Health Index Advisory Group in Scotland. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

LY and YZ had full access to all the data in the study and took responsibility for the integrity of the data and the accuracy of the data analysis. SY and XL conceived and designed the study. LY and YZ undertook the statistical analyses. LY wrote the first draft of the manuscript. All authors made critical revisions of the manuscript for important intellectual content.

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.

Publisher’s Note

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.

Acknowledgments

This research was conducted using the UK Biobank study under the application number 66354. The authors thank the Neale Lab and FinnGen consortium for sharing the summary-level data on gallstones.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.814457/full#supplementary-material

References

Bao, X., Borné, Y., Muhammad, I. F., Nilsson, J., Lind, L., Melander, O., et al. (2019). Growth Differentiation Factor 15 Is Positively Associated with Incidence of Diabetes Mellitus: the Malmö Diet and Cancer-Cardiovascular Cohort. Diabetologia 62 (1), 78–86. doi:10.1007/s00125-018-4751-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Bootcov, M. R., Bauskin, A. R., Valenzuela, S. M., Moore, A. G., Bansal, M., He, X. Y., et al. (1997). MIC-1, a Novel Macrophage Inhibitory Cytokine, Is a Divergent Member of the TGF-β Superfamily. Proc. Natl. Acad. Sci. U.S.A. 94 (21), 11514–11519. doi:10.1073/pnas.94.21.11514

PubMed Abstract | CrossRef Full Text | Google Scholar

Bowden, J., Davey Smith, G., and Burgess, S. (2015). Mendelian Randomization with Invalid Instruments: Effect Estimation and Bias Detection through Egger Regression. Int. J. Epidemiol. 44 (2), 512–525. doi:10.1093/ije/dyv080

PubMed Abstract | CrossRef Full Text | Google Scholar

Bowden, J., Davey Smith, G., Haycock, P. C., and Burgess, S. (2016). Consistent Estimation in Mendelian Randomization with Some Invalid Instruments Using a Weighted Median Estimator. Genet. Epidemiol. 40 (4), 304–314. doi:10.1002/gepi.21965

PubMed Abstract | CrossRef Full Text | Google Scholar

Breit, S. N., Johnen, H., Cook, A. D., Tsai, V. W. W., Mohammad, M. G., Kuffner, T., et al. (2011). The TGF-β Superfamily Cytokine, MIC-1/GDF15: A Pleotrophic Cytokine with Roles in Inflammation, Cancer and Metabolism. Growth factors. 29 (5), 187–195. doi:10.3109/08977194.2011.607137

PubMed Abstract | CrossRef Full Text | Google Scholar

Brown, D. A., Breit, S. N., Buring, J., Fairlie, W. D., Bauskin, A. R., Liu, T., et al. (2002). Concentration in Plasma of Macrophage Inhibitory Cytokine-1 and Risk of Cardiovascular Events in Women: a Nested Case-Control Study. Lancet 359 (9324), 2159–2163. doi:10.1016/s0140-6736(02)09093-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Burgess, S., and Thompson, S. G. (2015). Mendelian Randomization: Methods for Using Genetic Variants in Causal Estimation. Boca Raton: CRC Press.

Google Scholar

Chrysovergis, K., Wang, X., Kosak, J., Lee, S.-H., Kim, J. S., Foley, J. F., et al. (2014). NAG-1/GDF-15 Prevents Obesity by Increasing Thermogenesis, Lipolysis and Oxidative Metabolism. Int. J. Obes. 38 (12), 1555–1564. doi:10.1038/ijo.2014.27

PubMed Abstract | CrossRef Full Text | Google Scholar

Clarke, L., Zheng-Bradley, X., Zheng-Bradley, X., Smith, R., Kulesha, E., Xiao, C., et al. (2012). The 1000 Genomes Project: Data Management and Community Access. Nat. Methods 9 (5), 459–462. doi:10.1038/nmeth.1974

PubMed Abstract | CrossRef Full Text | Google Scholar

Coll, A. P., Chen, M., Taskar, P., Rimmington, D., Patel, S., Tadross, J. A., et al. (2020). GDF15 Mediates the Effects of Metformin on Body Weight and Energy Balance. Nature 578 (7795), 444–448. doi:10.1038/s41586-019-1911-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Conte, M., Martucci, M., Mosconi, G., Chiariello, A., Cappuccilli, M., Totti, V., et al. (2020). GDF15 Plasma Level Is Inversely Associated with Level of Physical Activity and Correlates with Markers of Inflammation and Muscle Weakness. Front. Immunol. 11, 915. doi:10.3389/fimmu.2020.00915

PubMed Abstract | CrossRef Full Text | Google Scholar

Day, E. A., Ford, R. J., Smith, B. K., Mohammadi-Shemirani, P., Morrow, M. R., Gutgesell, R. M., et al. (2019). Metformin-induced Increases in GDF15 Are Important for Suppressing Appetite and Promoting Weight Loss. Nat. Metab. 1 (12), 1202–1208. doi:10.1038/s42255-019-0146-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Dorvash, M. R., Khoshnood, M. J., Saber, H., Dehghanian, A., Mosaddeghi, P., and Firouzabadi, N. (2018). Metformin Treatment Prevents Gallstone Formation but Mimics Porcelain Gallbladder in C57Bl/6 Mice. Eur. J. Pharmacol. 833, 165–172. doi:10.1016/j.ejphar.2018.06.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Everhart, J. E., and Ruhl, C. E. (2009). Burden of Digestive Diseases in the United States Part III: Liver, Biliary Tract, and Pancreas. Gastroenterology 136 (4), 1134–1144. doi:10.1053/j.gastro.2009.02.038

PubMed Abstract | CrossRef Full Text | Google Scholar

Folkersen, L., Gustafsson, S., Wang, Q., Hansen, D. H., Hedman, Å. K., Schork, A., et al. (2020). Genomic and Drug Target Evaluation of 90 Cardiovascular Proteins in 30,931 Individuals. Nat. Metab. 2 (10), 1135–1148. doi:10.1038/s42255-020-00287-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Georgakis, M. K., Malik, R., Li, X., Gill, D., Levin, M. G., Vy, H. M. T., et al. (2021). Genetically Downregulated Interleukin-6 Signaling Is Associated with a Favorable Cardiometabolic Profile. Circulation 143 (11), 1177–1180. doi:10.1161/circulationaha.120.052604

PubMed Abstract | CrossRef Full Text | Google Scholar

Gerstein, H. C., Pare, G., Hess, S., Ford, R. J., Sjaarda, J., Raman, K., et al. (2017). Growth Differentiation Factor 15 as a Novel Biomarker for Metformin. Diabetes Care 40 (2), 280–283. doi:10.2337/dc16-1682

PubMed Abstract | CrossRef Full Text | Google Scholar

Griffin, S. J., Leaver, J. K., and Irving, G. J. (2017). Impact of Metformin on Cardiovascular Disease: a Meta-Analysis of Randomised Trials Among People with Type 2 Diabetes. Diabetologia 60 (9), 1620–1629. doi:10.1007/s00125-017-4337-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Hemani, G., Zheng, J., Elsworth, B., Wade, K. H., Haberland, V., Baird, D., et al. (2018). The MR-Base Platform Supports Systematic Causal Inference across the Human Phenome. Elife 7, e34408. doi:10.7554/eLife.34408

PubMed Abstract | CrossRef Full Text | Google Scholar

Hong, J. H., Chung, H. K., Park, H. Y., Joung, K.-H., Lee, J. H., Jung, J. G., et al. (2014). GDF15 Is a Novel Biomarker for Impaired Fasting Glucose. Diabetes Metab. J. 38 (6), 472–479. doi:10.4093/dmj.2014.38.6.472

PubMed Abstract | CrossRef Full Text | Google Scholar

Isik, S., Ozcan, H. N., Ozuguz, U., Berker, D., Tutuncu, Y., Akbaba, G., et al. (2012). Impaired Gallbladder Motility and the Effect of Metformin Therapy in Patients with Polycystic Ovary Syndrome. Clin. Endocrinol. (Oxf) 76 (3), 373–378. doi:10.1111/j.1365-2265.2011.04210.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnen, H., Lin, S., Kuffner, T., Brown, D. A., Tsai, V. W.-W., Bauskin, A. R., et al. (2007). Tumor-induced Anorexia and Weight Loss Are Mediated by the TGF-β Superfamily Cytokine MIC-1. Nat. Med. 13 (11), 1333–1340. doi:10.1038/nm1677

PubMed Abstract | CrossRef Full Text | Google Scholar

Kamat, M. A., Blackshaw, J. A., Young, R., Surendran, P., Burgess, S., Danesh, J., et al. (2019). PhenoScanner V2: an Expanded Tool for Searching Human Genotype-Phenotype Associations. Bioinformatics 35, 4851–4853. doi:10.1093/bioinformatics/btz469

PubMed Abstract | CrossRef Full Text | Google Scholar

Kempf, T., Guba-Quint, A., Torgerson, J., Magnone, M. C., Haefliger, C., Bobadilla, M., et al. (2012). Growth Differentiation Factor 15 Predicts Future Insulin Resistance and Impaired Glucose Control in Obese Nondiabetic Individuals: Results from the XENDOS Trial. Eur. J. Endocrinol. 167 (5), 671–678. doi:10.1530/eje-12-0466

PubMed Abstract | CrossRef Full Text | Google Scholar

Lammert, F., Gurusamy, K., Ko, C. W., Miquel, J.-F., Méndez-Sánchez, N., Portincasa, P., et al. (2016). Gallstones. Nat. Rev. Dis. Prim. 2, 16024. doi:10.1038/nrdp.2016.24

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Yang, L., Qin, W., Zhang, G., Yuan, J., and Wang, F. (2013). Adaptive Induction of Growth Differentiation Factor 15 Attenuates Endothelial Cell Apoptosis in Response to High Glucose Stimulus. PLoS One 8 (6), e65549. doi:10.1371/journal.pone.0065549

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Z., Liu, Y., Li, X., Wu, Y., Yang, F., Mo, Q., et al. (2020). Association between Circulating Growth Differentiation Factor 15 and Cirrhotic Primary Biliary Cholangitis. BioMed Res. Int. 2020, 1–11. doi:10.1155/2020/5162541

PubMed Abstract | CrossRef Full Text | Google Scholar

Liao, K.-F., Chuang, H.-Y., and Lai, S.-W. (2017). Metformin Use Correlates with Reduced Risk of Gallstones in Diabetic Patients: A 12-Year Follow-Up Study. Front. Pharmacol. 8, 765. doi:10.3389/fphar.2017.00765

PubMed Abstract | CrossRef Full Text | Google Scholar

Lieber, M. M. (1952). The Incidence of Gallstones and Their Correlation with Other Diseases. Ann. Surg. 135 (3), 394–405. doi:10.1097/00000658-195203000-00011

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Z., Kemp, T. J., Gao, Y.-T., Corbel, A., McGee, E. E., Wang, B., et al. (2018). Association of Circulating Inflammation Proteins and Gallstone Disease. J. Gastroenterology Hepatology 33 (11), 1920–1924. doi:10.1111/jgh.14265

CrossRef Full Text | Google Scholar

Luan, H. H., Wang, A., Hilliard, B. K., Carvalho, F., Rosen, C. E., Ahasic, A. M., et al. (2019). GDF15 Is an Inflammation-Induced Central Mediator of Tissue Tolerance. Cell 178 (5), 1231–1244. doi:10.1016/j.cell.2019.07.033

PubMed Abstract | CrossRef Full Text | Google Scholar

Macia, L., Tsai, V. W.-W., Nguyen, A. D., Johnen, H., Kuffner, T., Shi, Y.-C., et al. (2012). Macrophage Inhibitory Cytokine 1 (MIC-1/GDF15) Decreases Food Intake, Body Weight and Improves Glucose Tolerance in Mice on Normal & Obesogenic Diets. PLoS One 7 (4), e34868. doi:10.1371/journal.pone.0034868

PubMed Abstract | CrossRef Full Text | Google Scholar

Maurer, K. J., Carey, M. C., and Fox, J. G. (2009). Roles of Infection, Inflammation, and the Immune System in Cholesterol Gallstone Formation. Gastroenterology 136 (2), 425–440. doi:10.1053/j.gastro.2008.12.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakayasu, E. S., Syed, F., Tersey, S. A., Gritsenko, M. A., Mitchell, H. D., Chan, C. Y., et al. (2020). Comprehensive Proteomics Analysis of Stressed Human Islets Identifies GDF15 as a Target for Type 1 Diabetes Intervention. Cell Metab. 31 (2), 363–374. doi:10.1016/j.cmet.2019.12.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Nervi, F., Miquel, J. F., Alvarez, M., Ferreccio, C., García-Zattera, M. J., González, R., et al. (2006). Gallbladder Disease Is Associated with Insulin Resistance in a High Risk Hispanic Population. J. Hepatology 45 (2), 299–305. doi:10.1016/j.jhep.2006.01.026

CrossRef Full Text | Google Scholar

Özkan, H., Demirbaş, S., Ibiş, M., Akbal, E., and Köklü, S. (2011). Diagnostic Validity of Serum Macrophage Inhibitor Cytokine and Tissue Polypeptide-specific Antigen in Pancreatobiliary Diseases. Pancreatology 11 (3), 295–300. doi:10.1159/000328963

PubMed Abstract | CrossRef Full Text | Google Scholar

Pauletzki, J., Althaus, R., Holl, J., Sackmann, M., and Paumgartner, G. (1996). Gallbladder Emptying and Gallstone Formation: a Prospective Study on Gallstone Recurrence. Gastroenterology 111 (3), 765–771. doi:10.1053/gast.1996.v111.pm8780583

PubMed Abstract | CrossRef Full Text | Google Scholar

Schlender, L., Martinez, Y. V., Adeniji, C., Reeves, D., Faller, B., Sommerauer, C., et al. (2017). Efficacy and Safety of Metformin in the Management of Type 2 Diabetes Mellitus in Older Adults: a Systematic Review for the Development of Recommendations to Reduce Potentially Inappropriate Prescribing. BMC Geriatr. 17 (Suppl. 1), 227. doi:10.1186/s12877-017-0574-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Suissa, S., and Azoulay, L. (2012). Metformin and the Risk of Cancer: Time-Related Biases in Observational Studies. Diabetes Care 35 (12), 2665–2673. doi:10.2337/dc12-0788

PubMed Abstract | CrossRef Full Text | Google Scholar

Unsicker, K., Spittau, B., and Krieglstein, K. (2013). The Multiple Facets of the TGF-β Family Cytokine Growth/differentiation Factor-15/macrophage Inhibitory Cytokine-1. Cytokine & Growth Factor Rev. 24 (4), 373–384. doi:10.1016/j.cytogfr.2013.05.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Verbanck, M., Chen, C.-Y., Neale, B., and Do, R. (2018). Detection of Widespread Horizontal Pleiotropy in Causal Relationships Inferred from Mendelian Randomization between Complex Traits and Diseases. Nat. Genet. 50 (5), 693–698. doi:10.1038/s41588-018-0099-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Vila, G., Riedl, M., Anderwald, C., Resl, M., Handisurya, A., Clodi, M., et al. (2011). The Relationship between Insulin Resistance and the Cardiovascular Biomarker Growth Differentiation Factor-15 in Obese Patients. Clin. Chem. 57 (2), 309–316. doi:10.1373/clinchem.2010.153726

PubMed Abstract | CrossRef Full Text | Google Scholar

Wischhusen, J., Melero, I., and Fridman, W. H. (2020). Growth/Differentiation Factor-15 (GDF-15): From Biomarker to Novel Targetable Immune Checkpoint. Front. Immunol. 11, 951. doi:10.3389/fimmu.2020.00951

PubMed Abstract | CrossRef Full Text | Google Scholar

Wollert, K. C., Kempf, T., and Wallentin, L. (2017). Growth Differentiation Factor 15 as a Biomarker in Cardiovascular Disease. Clin. Chem. 63 (1), 140–151. doi:10.1373/clinchem.2016.255174

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiao, Q.-A., He, Q., Zeng, J., and Xia, X. (2022). GDF-15, a Future Therapeutic Target of Glucolipid Metabolic Disorders and Cardiovascular Disease. Biomed. Pharmacother. 146, 112582. doi:10.1016/j.biopha.2021.112582

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, S., Gill, D., Giovannucci, E. L., and Larsson, S. C. (2022). Obesity, Type 2 Diabetes, Lifestyle Factors, and Risk of Gallstone Disease: A Mendelian Randomization Investigation. Clin. Gastroenterology Hepatology 20, e529–e537. doi:10.1016/j.cgh.2020.12.034

CrossRef Full Text | Google Scholar

Keywords: gallstones, growth differentiation factor 15, metformin, Mendelian randomization, diabetic complication

Citation: Yu L, Zhou Y, Wang L, Zhou X, Sun J, Xiao J, Xu X, Larsson SC, Yuan S and Li X (2022) GDF-15 as a Therapeutic Target of Diabetic Complications Increases the Risk of Gallstone Disease: Mendelian Randomization and Polygenic Risk Score Analysis. Front. Genet. 13:814457. doi: 10.3389/fgene.2022.814457

Received: 13 November 2021; Accepted: 02 May 2022;
Published: 13 June 2022.

Edited by:

Tao-Hsin Tung, Taizhou Hospital of Zhejiang Province Affiliated to Wenzhou Medical University, China

Reviewed by:

Nicholas B. Larson, Mayo Clinic, United States
Gregory Livshits, Ariel University, Israel

Copyright © 2022 Yu, Zhou, Wang, Zhou, Sun, Xiao, Xu, Larsson, Yuan and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Xue Li, xueli157@zju.edu.cn

These authors share first authorship

These authors share last authorship

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.