Hua Xu1,2†
Chuandi Jin
Qingbo Guan- 1Department of Endocrinology, Shandong Provincial Hospital Affiliated to Shandong University, Jinan, China
- 2Shandong Clinical Medical Center of Endocrinology and Metabolism, Institute of Endocrinology and Metabolism, Shandong Academy of Clinical Medicine, Jinan, China
- 3Institute for Medical Dataology, Shandong University, Jinan, China
Overall and abdominal obesity were significantly associated with insulin resistance and type 2 diabetes mellitus (T2DM) risk in observational studies, though these associations cannot avoid the bias induced by confounding effects and reverse causation. This study aimed to test whether these associations are causal, and it compared the causal effects of overall and abdominal obesity on T2DM risk and glycemic traits by using a two-sample Mendelian randomization (MR) design. Based on summary-level statistics from genome-wide association studies, the instrumental variables for body mass index (BMI), waist-to-hip ratio (WHR), and WHR adjusted for BMI (WHRadjBMI) were extracted, and the horizontal pleiotropy was analyzed using MR–Egger regression and the MR–pleiotropy residual sum and outlier (PRESSO) method. Thereafter, by using the conventional MR method, the inverse-variance weighted method was applied to assess the causal effect of BMI, WHR, and WHRadjBMI on T2DM risk, Homeostatic model assessment of insulin resistance (HOMA-IR), fasting insulin, fasting glucose, and Hemoglobin A1c (HbA1c). A series of sensitivity analyses, including the multivariable MR (diastolic blood pressure, systolic blood pressure, high-density lipoprotein cholesterol, and low-density lipoprotein cholesterol as covariates), MR–Egger regression, weighted median, MR–PRESSO, and leave-one-out method, were conducted to test the robustness of the results from the conventional MR. Despite the existence of horizontal pleiotropy, consistent results were found in the conventional MR results and sensitivity analyses, except for the association between BMI and fasting glucose, and WHRadjBMI and fasting glucose. Each one standard deviation higher BMI was associated with an increased T2DM risk [odds ratio (OR): 2.741; 95% confidence interval (CI): 2.421–3.104], higher HbA1c [1.054; 1.04–1.068], fasting insulin [1.202; 1.173–1.231], and HOMA-IR [1.221; 1.187–1.255], similar to findings for causal effect of WHRadjBMI on T2DM risk [1.993; 1.704–2.33], HbA1c [1.061; 1.042–1.08], fasting insulin [1.102; 1.068–1.136], and HOMA-IR [1.127; 1.088–1.167]. Both BMI (P = 0.546) and WHRadjBMI (P = 0.443) were unassociated with fasting glucose in the multivariable MR analysis. In conclusion, overall and abdominal obesity have causal effects on T2DM risk and insulin resistance but no causal effect on fasting glucose. Individuals can substantially reduce their insulin resistance and T2DM risk through reduction of body fat mass and modification of body fat distribution.
Introduction
Type 2 diabetes mellitus (T2DM) is a chronic metabolic disease characterized by hyperglycemia secondary to insulin resistance and pancreatic β-cell failure (Alejandro et al., 2015). The findings of human epidemiologic studies indicate that the global prevalence of T2DM is increasing rapidly, and this increase parallels the increase in the prevalence of obesity (Sampath Kumar et al., 2019). The body mass index (BMI) is routinely used to quantify the overall obesity although body fat distribution of individuals can vary substantially. The waist-to-hip ratio (WHR) and WHR adjusted for BMI (WHRadjBMI) are frequently used surrogate measures of abdominal obesity. Many observational epidemiologic studies, including case–control and cohort studies, have demonstrated that higher WHR and BMI are two important risk factors for developing T2DM (Vazquez et al., 2007; Lv et al., 2017). Moreover, cohort and cross-sectional studies (Wang et al., 2018; Benites-Zapata et al., 2019) demonstrated that BMI and WHR were associated with glycemic traits, including fasting insulin, Hemoglobin A1c (HbA1c), and insulin resistance [measured by Homeostatic model assessment of insulin resistance (HOMA-IR)]. Longitudinal and cross-sectional studies have found an association between increased risk of T2DM and higher genetic predisposition to both BMI and WHRadjBMI in European and East Asian populations (Robiou-du-Pont et al., 2013; Zhu et al., 2014; Huang et al., 2015).
However, these observational studies cannot avoid the bias induced by the confounding effect and reverse causation and, therefore, are incapable of confirming whether these associations are causal (Smith and Ebrahim, 2003). Mendelian randomization (MR) is an approach that is used to unbiasedly test or estimate the causal relationship between an exposure and an associated outcome by using data on inherited genetic variants that influence exposure status in the presence of unmeasured confounding (Didelez and Sheehan, 2007; Lawlor et al., 2008). In the past few years, MR has been extensively used in epidemiology and other related areas of population science (Smit et al., 2019; Wainberg et al., 2019; Yang et al., 2019).
Previous MR analyses have demonstrated that per 1 standard deviation (SD) higher WHRadjBMI and BMI were causally associated with T2DM risk in European populations (Dale et al., 2017; Emdin et al., 2017). Wang et al. (2018) conducted an MR to further investigate the causal effect of both BMI and WHR on glycemic traits, and they found that BMI had a causal relevance for insulin secretion, whereas neither WHR and BMI was causally associated with HOMA-IR in a conventional MR in a Chinese Han population. However, there is a dearth of MR studies for testing and comparing the causal effect of both overall and abdominal obesity on glycemic in the European population. Epidemiologic studies have found differences in T2DM epidemiologic characteristics between the Asian and European population wherein, in comparison with South Asians, Europeans have a lower T2DM risk, typically develop T2DM 5–10 years later, and have a slower disease progression (Gujral et al., 2013; Admiraal et al., 2014; Gupta and Misra, 2016; Meeks et al., 2016; Banerjee and Shah, 2018). Moreover, Europeans developed many metabolic abnormalities, including hyperglycemia and elevated triacylglycerol and low high-density lipoprotein cholesterol (HDL-C), at a higher BMI and age (Raji et al., 2001; Razak et al., 2007). Therefore, the estimation and comparison of the causal effects of overall and abdominal obesity on glycemic traits could provide insights into the obesity-related mechanism of T2DM.
In this study, a two-sample Mendelian randomization (TSMR) with a large sample size was conducted to determine whether a genetic predisposition to increased BMI, WHR, and WHRadjBMI was causally associated with T2DM and glycemic traits, including HOMA-IR, fasting insulin, fasting glucose, and HbA1c. The causal effects were further compared to identify differences in the effect of overall and abdominal obesity on T2DM development and glycemic traits.
Materials and Methods
Data Source
This study aimed to explore the causal effect of WHR, BMI, and WHRadjBMI on the risk of T2DM and glycemic traits (HOMA-IR, fasting insulin, fasting glucose, and HbA1c) in an European population, and used diastolic blood pressure (DBP), systolic blood pressure (SBP), HDL-C, and low-density lipoprotein cholesterol (LDL-C) as the covariates. The genome-wide association study summary statistics datasets used in this study were obtained from Zenodo1 for WHR (Censin et al., 2019), BMI (Censin et al., 2019), and WHRadjBMI (Censin et al., 2019); the Program in Complex Trait Genomics2 for T2DM (Xue et al., 2018); MAGIC Consortium3 for HOMA-IR (Dupuis et al., 2010), fasting glucose (Manning et al., 2012), fasting insulin (Manning et al., 2012), and HbA1c (Wheeler et al., 2017); the MRBase platform4 for HDL-C (Kettunen et al., 2016) and LDL-C (Kettunen et al., 2016); and the MRC-IEU Consortium5 for SBP and DBP. Detailed information of the summary statistics datasets are displayed in Table 1. We obtained the β-coefficients and standard errors for the per allele association of each single-nucleotide polymorphism (SNP) as well as all exposures and outcomes from these data sources.
Table 1. Summary statistics of data source.
Selection of Genetic Instrumental Variables
In the TSMR analysis conducted in this study, the genetic variants for exposures (BMI, WHR, and WHRadjBMI) were used as instrumental variables (IVs) and were obtained by two steps. Firstly, SNPs that are strongly associated with exposures (P < 5.0 × 10–8) were extracted. Secondly, we pruned these extracted SNPs by linkage disequilibrium (LD; r2 = 0.001, clumping distance = 10,000 kb) to ensure that each IV was independent of the others. To test the strength of the IVs, the F-statistics were calculated as previously described (Xu and Hao, 2017). F-statistics >10 are considered adequately strong to mitigate against any bias of the causal IV estimate.
Heterogeneity and Horizontal Pleiotropic Analysis
In MR, heterogeneity in the causal estimate may indicate that a variant has an effect on the outcome outside of its effect on the exposure (known as horizontal pleiotropy), and this can cause severe bias (Davey Smith and Hemani, 2014). Mendelian randomization–Egger (MR–Egger) regression was undertaken to assess the horizontal pleiotropy of the IVs, where a regression intercept that significantly differed from zero (P < 0.05) indicated the presence of horizontal pleiotropy exists or that the InSIDE (INstrument Strength Independent of Direct Effect) assumption was violated (Bowden et al., 2015). Heterogeneity between IVs in the conventional MR, with the inverse-variance weighted (IVW) method, was estimated by Cochran’s Q statistic. The MR pleiotropy residual sum and outlier (MR–PRESSO) method can be used to test horizontal pleiotropic outliers and can obtain the corrected causal effect after removal of these outliers in MR (Verbanck et al., 2018). In the present study, both MR–Egger regression and MR–PRESSO tests were conducted using the TwoSampleMR and MRPRESSO R package in R (version 3.6.1), respectively.
Mendelian Randomization
Mendelian randomization can test and estimate the causal effect of an exposure on an outcome by using genetic variants as the IVs (Zheng et al., 2017). Firstly, Wald ratios were calculated for each IV by dividing the per-allele log-odds ratio (or beta) of that variant in the outcome data by the log-odds ratio (or beta) of the same variant in the exposure data. Then, the random-effects IVW method was applied to estimate the association between exposures and outcomes. In IVW, the Wald ratio for each SNP was weighted by its inverse variance, and the effect estimates were meta-analyzed using random effects.
Sensitivity Analysis
Sensitivity analysis was used to test the disproportionate effects of variants and the pleiotropy in the MR analysis (Mokry et al., 2016). A series of sensitivity analyses were conducted to test the robustness of the conventional MR results.
Multivariable IVW, which included the DBP, SBP, HDL-C, and LDL-C as covariates, was carried out in accordance with the method proposed by Rees et al. (2017) that was used to account for possible horizontal pleiotropy arising from the association of the instrument with these variables.
The MR–Egger regression and weighted median method are two pleiotropy-robust MR methods that are used to estimate consistent causal effects against unknown directional pleiotropy under the InSIDE assumptions (Bowden et al., 2015). In the MR–Egger regression method, the regression line fitted to the data is not constrained to pass through the origin, and the intercept represents the horizontal pleiotropic effect that may bias the IVW estimate, whereas the slope represents pleiotropy-corrected causal estimates. The weighted median method has considerable robustness to individual genetics with strongly outlying causal estimates and could provide a consistent causal estimate when the valid IVs exceed 50%.
The MR–PRESSO method was used to identify potential outliers in the conventional MR testing, and provided a robust estimate with outlier correction. Moreover, testing of significant distortion in the IVW causal estimate before and after MR–PRESSO correction, was undertaken and served as a sensitivity analysis.
The leave-one-out sensitivity analysis was conducted to ascertain whether the association was being disproportionately influenced by a single SNP. In this analysis, the random-effects IVW was repeated by leaving out each SNP in turn, and the overall analysis including all SNPs was used for the comparison. The variation of the results from before and after the removal of each SNP reflects the sensitivity of this SNP.
Results
Genetic IVs
A total of 546, 356, and 330 IVs were identified for BMI, WHR, and WHRadjBMI, respectively. Some IVs were absent in the outcome data; however, the F statistics for BMI-IVs (86.250–89.078), WHR-IVs (67.502–67.991), and WHRadjBMI-IVs (90.758–96.860) that were used for MR were more than 10, which indicated that the weak instrument bias was negligible. Detailed information of IVs used in this study are shown in Supplementary Table S1.
Horizontal Pleiotropy and Heterogeneity Analysis
The MR–Egger regression intercepts obtained in this study (Table 2) showed that horizontal pleiotropy (P = 0.029) was only found in the MR with WHRadjBMI as exposure and fasting insulin as the outcome. Heterogeneity (Table 2) was observed between IVs; therefore, the random-effect IVW method was used in the subsequent stages of the research analysis. The MR–PRESSO test showed that horizontal pleiotropy was found in all IVW analyses in this study, and the horizontal pleiotropic outliers were identified and removed (Supplementary Table S2). After the removal of these outliers, the F-statistics of BMI-IVs (85.704–89.033), WHR-IVs (62.240–67.991), and WHRadjBMI-IVs (84.668–96.860) continued to remain well powered to estimate the causal effect of the exposure on the outcome.
Table 2. Heterogeneity and horizontal pleiotropy analysis.
Causal Effect of WHR, BMI, and WHRadjBMI on T2DM and Glycemic Traits
Table 3 and Figure 1 show the causal effect estimates of WHR, BMI, and WHRadjBMI on T2DM and glycemic traits. The TSMR analysis by the IVW method showed a significant causal effect, wherein each SD of genetically higher BMI was associated with an increased T2DM risk [OR: 2.741; 95% confidence interval (CI): 2.421–3.104], higher fasting glucose [1.073; 1.048–1.099], higher fasting insulin [1.202; 1.173–1.231], higher HOMA-IR [1.221; 1.187–1.255], and higher HbA1c [1.054; 1.04–1.068]. Each SD of genetically higher WHR was associated with increased T2DM risk [3.12; 2.653–3.668], higher fasting glucose [1.087; 1.054–1.12], higher fasting insulin [1.193; 1.153–1.234], higher HOMA-IR [1.203; 1.155–1.252], and higher HbA1c [1.075; 1.056–1.095]. Each SD of genetically higher WHRadjBMI was associated with increased T2DM risk [1.993; 1.704–2.33], higher fasting glucose [1.039; 1.012–1.067], higher fasting insulin [1.102; 1.068–1.136], higher HOMA-IR [1.127; 1.088–1.167], and higher HbA1c [1.061; 1.042–1.08].
Table 3. Mendelian randomization results.
Figure 1. Scatterplot and causal effect of body mass index (BMI), waist-to-hip ratio (WHR), and WHR adjusted for BMI (WHRadjBMI) on type 2 diabetes mellitus (T2DM) and glycemic traits including HOMA-IR, fasting insulin, fasting glucose and glycate hemoglobin (HbA1c). (A) Causal effect of WHRadjBMI on T2DM, (B) Causal effect of WHRadjBMI on HOMA-IR, (C) Causal effect of WHRadjBMI on fasting glucose, (D) Causal effect of WHRadjBMI on fasting insulin, (E) Causal effect of WHRadjBMI on HbA1c, (F) Causal effect of WHR on T2DM, (G) Causal effect of WHR on HOMA-IR, (H) Causal effect of WHR on fasting glucose, (I) Causal effect of WHR on fasting insulin, (J) Causal effect of WHR on HbA1c, (K) Causal effect of BMI on T2DM, (L) Causal effect of BMI on HOMA-IR, (M) Causal effect of BMI on fasting glucose, (N) Causal effect of BMI on fasting insulin, and (O) Causal effect of BMI on HbA1c. The x-axis presents the single nucleotide polymorphism (SNP) effect on exposure, and the y-axis presents the SNP effect on outcome. The light blue, dark blue and green regression line represents the inverse-variance weighted (IVW), Mendelian randomization (MR)–Egger, and weighted median estimate, respectively.
Sensitivity Analysis
In the leave-one-out sensitivity analysis, no single SNP strongly or reversely drove the overall effect of exposure on outcome in the IVW (Supplementary Figure S1). Consistent results were observed in the IVW after the MR–PRESSO correction, MR–Egger regression, and the weighted median method, with the exception of the causal estimates of WHR on HbA1c (P = 0.058) and fasting glucose (P = 0.098) in the MR–Egger regression. The MR–Egger regression could obtain pleiotropy-corrected causal estimates, although this method had less statistical power than an equivalent IVW method, and the CIs were wider and included the null value (Bowden, 2017; Weng et al., 2018). Because the intercept of the MR–Egger regression indicates that there was no horizontal pleiotropy in the MR–Egger regression between WHR and both HbA1c (P = 0.695) and fasting glucose (P = 0.935), the causal estimate was more convincing in the IVW. In the multivariable IVW (DBP, SBP, HDL-C, and LDL-C as covariates), BMI (P = 0.546) and WHRadjBMI (P = 0.443) were not causally associated with fasting glucose, whereas other multivariable IVW results persisted with that in the univariable IVW.
Taken together, the causal effect estimates of BMI and WHRadjBMI on fasting glucose in conventional MR might be biased by the horizontal pleiotropy of SBP, DBP, HDL-C, and LDL-C, while no significant bias was found in other causal effect estimates despite the existence of horizontal pleiotropy and heterogeneity.
Discussion
Numerous observational studies indicated that obesity was strongly associated with T2DM risk and glycemic traits (Lv et al., 2017), however, a causal effect cannot be ascertained from these studies due to residual confounding or reverse causality. This present study utilized a TSMR design that was applied to the summary-level data from a large-scale genome-wide association study to address the potential causal role of overall obesity (measured by BMI) and abdominal obesity (measure by WHRadjBMI) on the risk of T2DM and glycemic traits. The well-powered conventional MR (random-effect IVW method) confirmed that genetic predisposition to higher BMI, WHR, and WHRadjBMI are causally associated with higher fasting glucose, fasting insulin, HOMA-IR, HbA1c, and increased risk of T2DM in the European population.
However, heterogeneity and horizontal pleiotropy was found in the conventional MR analysis, a series of sensitivity analyses that included the multivariable MR (DBP, SBP, HDL-C, and LDL-C as covariates), MR–Egger regression, weighted median method, MR–PRESSO method, and leave-one-out analysis to test the robustness of the conventional MR results. The causal effect of BMI and WHRadjBMI on T2DM risk, HbA1c, fasting insulin, and HOMA-IR in the conventional MR were consistent with that in all the sensitivity analyses, which suggested that the causal estimate was robust and unbiased.
Each SD of genetically higher BMI [2.741; 2.421–3.104] and WHRadjBMI [1.993; 1.704–2.33] was associated with increased T2DM risk. Human epidemiologic studies have considered obesity to be a major risk factor of T2DM, and the substantial increase in the incidence of obesity contributes to the current T2DM epidemic (Sampath Kumar et al., 2019). Using the MR method in the European descendants, Emdin et al. confirmed that a 1 SD genetic increase in WHRadjBMI was associated with a higher risk of T2DM [1.77; 1.57–2.00] (Emdin et al., 2017), and Dale et al. revealed that each SD higher BMI was associated with increased T2DM risk [1.98; 1.41–2.78] (Dale et al., 2017). The results of this study are in agreement with those of previous observational studies (Lv et al., 2017; Sampath Kumar et al., 2019) and MR studies (Dale et al., 2017; Emdin et al., 2017) which suggested that both overall and abdominal obesity play a causal role on T2DM risk in the European population. In addition, our MR studies suggested that the causal effect of overall obesity on T2DM risk was greater than that of abdominal obesity. Moreover, both BMI [1.054; 1.04–1.068] and WHRadjBMI [1.061; 1.042–1.08] were found to have a causal effect on HbA1c, which suggested that overall and abdominal obesity have a similar but small causal effect on the HbA1c.
Insulin resistance refers to a decreased physiological response of peripheral tissues to insulin action, which implies an impaired effect of insulin in lowering the blood glucose (Gelaye et al., 2010). This serves as the key mechanism and a major global driver of the T2DM condition (Roglic, 2016; Czech, 2017). The accumulation of body fat and abdominal body fat are risk factors for increased insulin resistance (Kohrt et al., 1993; Gobato et al., 2014), and high BMI and WHR were found to be positively correlated with insulin resistance in observational epidemiological studies (Gobato et al., 2014; Benites-Zapata et al., 2019; Lin et al., 2019). Wang et al. reported that higher BMI was causally correlated with increased Stumvoll first- and second-phase insulin secretion and HOMA-IR, whereas no causal relationship between WHR and HOMA-IR was found in a conventional MR study in the Chinese Han population (Wang et al., 2018). In Europeans, a previous MR study found that higher WHRadjBMI was causally associated with higher fasting insulin levels (Emdin et al., 2017). The present MR study provides a similar conclusion with regard to the European population, each SD of genetically higher WHRadjBMI [1.102; 1.068–1.136] and BMI [1.202; 1.173–1.231] was found to play a positive causal effect on higher fasting insulin. Furthermore, each SD of genetically higher BMI [1.221; 1.187–1.255] and WHRadjBMI [1.127; 1.088–1.167] was causally associated with the HOMA-IR. These results suggested that higher overall and abdominal obesity serve as causal risk factors of fasting insulin and insulin resistance in the European population. The findings of the present study are supported by experimental studies as well. Obesity could stimulate the formation of lipid metabolites, hormones, and cytokines, which involves changes in the insulin signaling pathway and the accelerated progression of insulin resistance (Patel and Abate, 2013; Balsan et al., 2015). Moreover, the causal effect of overall obesity on fasting insulin and insulin resistance is slightly greater than that of abdominal obesity. Thus, we highlighted that both mass and distribution of body fat play a causal role on insulin resistance and T2DM risk. This indicates that the development of therapies to modify the mass and distribution of body fat to reduce overall and abdominal obesity might contribute to the prevention and alleviation of T2DM and insulin resistance-related diseases.
Furthermore, although higher BMI and WHRadjBMI was found to be causally associated with higher fasting glucose in our conventional MR in the European population, no statistical significance was found between BMI and fasting glucose (P = 0.546) or with WHRadjBMI and fasting glucose (P = 0.443) in the multivariable MR (DBP, SBP, HDL-C, and LDL-C as covariates). The casual estimates of BMI and WHRadjBMI on fasting glucose in conventional MR might be biased by the horizontal pleiotropy of DBP, SBP, HDL-C, and LDL-C. These negative results warrant further investigation.
Through a comparison of the causal estimates of BMI and WHRadjBMI on glycemic traits (fasting glucose, fasting insulin, HOMA-IR, and HbA1c), this study further emphasizes that overall and abdominal obesity might increase the T2DM risk mainly via elevation of insulin resistance.
In conclusion, overall and abdominal obesity have a causal effect on the T2DM risk and insulin resistance, and overall obesity may have stronger effects, whereas they may have no causal effect on the fasting glucose. These results suggest that individuals can substantially reduce their insulin resistance and T2DM risk through reduction of body fat mass and modification of body fat distribution.
Data Availability Statement
Publicly available datasets were analyzed in this study. This data can be found here: The Zenodo (https://zenodo.org), the Program in Complex Trait Genomics (https://cnsgenomics.com), the MAGIC Consortium (http://www.magicinvestigators.org/), the MRBase platform (http://www.mrbase.org), and the MRC-IEU Consortium (http://www.bristol.ac.uk/integrative-epidemiology/).
Author Contributions
HX and CJ collected, analyzed, and interpreted the data. All authors conceived and designed the project and wrote and approved the final version fo the manuscript.
Funding
This research was supported by Key Technology Research and Development Program of Shandong Province (2017CXGC1214).
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.
Acknowledgments
We would like to thank Editage (www.editage.com) for English language editing.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.00603/full#supplementary-material
FIGURE S1 | Leave-one-out sensitive analysis.
Footnotes
- ^ https://zenodo.org
- ^ https://cnsgenomics.com
- ^ http://www.magicinvestigators.org/
- ^ http://www.mrbase.org
- ^ http://www.bristol.ac.uk/integrative-epidemiology/
References
Admiraal, W. M., Holleman, F., Snijder, M. B., Peters, R. J., Brewster, L. M., Hoekstra, J. B., et al. (2014). Ethnic disparities in the association of impaired fasting glucose with the 10-year cumulative incidence of type 2 diabetes. Diabetes Res. Clin. Pract. 103, 127–132. doi: 10.1016/j.diabres.2013.10.014
Alejandro, E. U., Gregg, B., Blandino-Rosano, M., Cras-Méneur, C., and Bernal-Mizrachi, E. (2015). Natural history of β-cell adaptation and failure in type 2 diabetes. Mol. Aspects Med. 42, 19–41. doi: 10.1016/j.mam.2014.12.002
Balsan, G. A., Vieira, J. L., Oliveira, A. M., and Portal, V. L. (2015). Relationship between adiponectin, obesity and insulin resistance. Rev. Assoc. Med. Bras. 61, 72–80. doi: 10.1590/1806-9282.61.01.072
Banerjee, A. T., and Shah, B. R. (2018). Differences in prevalence of diabetes among immigrants to Canada from South Asian countries. Diabet. Med. 35, 937–943. doi: 10.1111/dme.13647
Benites-Zapata, V. A., Toro-Huamanchumo, C. J., Urrunaga-Pastor, D., Guarnizo-Poma, M., Lazaro-Alcantara, H., Paico-Palacios, S., et al. (2019). High waist-to-hip ratio levels are associated with insulin resistance markers in normal-weight women. Diabetes Metab. Syndr. 13, 636–642. doi: 10.1016/j.dsx.2018.11.043
Bowden, J. (2017). Misconceptions on the use of MR-Egger regression and the evaluation of the InSIDE assumption. Int. J. Epidemiol. 46, 2097–2099. doi: 10.1093/ije/dyx192
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, 512–525. doi: 10.1093/ije/dyv080
Censin, J. C., Peters, S. A. E., Bovijn, J., Ferreira, T., Pulit, S. L., Magi, R., et al. (2019). Causal relationships between obesity and the leading causes of death in women and men. PLoS Genet. 15:e1008405. doi: 10.1371/journal.pgen.1008405
Czech, M. P. (2017). Insulin action and resistance in obesity and type 2 diabetes. Nat. Med. 23, 804–814. doi: 10.1038/nm.4350
Dale, C. E., Fatemifar, G., Palmer, T. M., White, J., Prieto-Merino, D., Zabaneh, D., et al. (2017). Causal associations of adiposity and body fat distribution with coronary heart disease, stroke subtypes, and Type 2 diabetes mellitus: a mendelian randomization analysis. Circulation 135, 2373–2388. doi: 10.1161/circulationaha.116.026560
Davey Smith, G., and Hemani, G. (2014). Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum. Mol. Genet. 23, R89–R98. doi: 10.1093/hmg/ddu328
Didelez, V., and Sheehan, N. (2007). Mendelian randomization as an instrumental variable approach to causal inference. Stat. Methods Med. Res. 16, 309–330. doi: 10.1177/0962280206077743
Dupuis, J., Langenberg, C., Prokopenko, I., Saxena, R., Soranzo, N., Jackson, A. U., et al. (2010). New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk. Nat. Genet. 42, 105–116. doi: 10.1038/ng.520
Emdin, C. A., Khera, A. V., Natarajan, P., Klarin, D., Zekavat, S. M., Hsiao, A. J., et al. (2017). Genetic association of waist-to-hip ratio with cardiometabolic traits. type 2 diabetes, and coronary heart disease. JAMA 317, 626–634. doi: 10.1001/jama.2016.21042
Gelaye, B., Revilla, L., Lopez, T., Suarez, L., Sanchez, S. E., Hevner, K., et al. (2010). Association between insulin resistance and c-reactive protein among Peruvian adults. Diabetol. Metab. Syndrome 2:30.
Gobato, A. O., Vasques, A. C., Zambon, M. P., Barros Filho Ade, A., and Hessel, G. (2014). Metabolic syndrome and insulin resistance in obese adolescents. Rev. Paul. Pediatr. 32, 55–62. doi: 10.1590/s0103-05822014000100010
Gujral, U. P., Pradeepa, R., Weber, M. B., Narayan, K. M., and Mohan, V. (2013). Type 2 diabetes in South Asians: similarities and differences with white Caucasian and other populations. Ann. N. Y. Acad. Sci. 1281, 51–63. doi: 10.1111/j.1749-6632.2012.06838.x
Gupta, R., and Misra, A. (2016). Epidemiology of microvascular complications of diabetes in South Asians and comparison with other ethnicities. J. Diabetes 8, 470–482. doi: 10.1111/1753-0407.12378
Huang, T., Qi, Q., Zheng, Y., Ley, S. H., Manson, J. E., Hu, F. B., et al. (2015). Genetic predisposition to central obesity and risk of type 2 diabetes: two independent cohort studies. Diabetes Care 38, 1306–1311. doi: 10.2337/dc14-3084
Kettunen, J., Demirkan, A., Wurtz, P., Draisma, H. H., Haller, T., Rawal, R., et al. (2016). Genome-wide study for circulating metabolites identifies 62 loci and reveals novel systemic effects of LPA. Nat. Commun. 7:11122. doi: 10.1038/ncomms11122
Kohrt, W. M., Kirwan, J. P., Staten, M. A., Bourey, R. E., King, D. S., and Holloszy, J. O. (1993). Insulin resistance in aging is related to abdominal obesity. Diabetes Metab. Res. Rev. 42, 273–281. doi: 10.2337/diab.42.2.273
Lawlor, D., Harbord, R., Jonathan, A. C. S., Timpson, N., and Davey-Smith, G. (2008). Mendelian randomization: using genes as instruments for making causal inferences in epidemiology. Stat. Med. 27, 1133–1163. doi: 10.1002/sim.3034
Lin, C., Chen, K., Zhang, R., Fu, W., Yu, J., Gao, L., et al. (2019). The prevalence, risk factors, and clinical characteristics of insulin resistance in Chinese patients with schizophrenia. Compr. Psychiatry 96:152145. doi: 10.1016/j.comppsych.2019.152145
Lv, J., Yu, C., Guo, Y., Bian, Z., Yang, L., Chen, Y., et al. (2017). Adherence to a healthy lifestyle and the risk of type 2 diabetes in Chinese adults. Int. J. Epidemiol. 46, 1410–1420. doi: 10.1093/ije/dyx074
Manning, A. K., Hivert, M. F., Scott, R. A., Grimsby, J. L., Bouatia-Naji, N., Chen, H., et al. (2012). A genome-wide approach accounting for body mass index identifies genetic variants influencing fasting glycemic traits and insulin resistance. Nat. Genet. 44, 659–669. doi: 10.1038/ng.2274
Meeks, K. A., Freitas-Da-Silva, D., Adeyemo, A., Beune, E. J., Modesti, P. A., Stronks, K., et al. (2016). Disparities in type 2 diabetes prevalence among ethnic minority groups resident in Europe: a systematic review and meta-analysis. Intern. Emerg. Med. 11, 327–340. doi: 10.1007/s11739-015-1302-9
Mokry, L. E., Ross, S., Timpson, N. J., Sawcer, S., Davey Smith, G., and Richards, J. B. (2016). Obesity and multiple sclerosis: a mendelian randomization study. PLoS Med. 13:e1002053. doi: 10.1371/journal.pmed.1002053
Patel, P., and Abate, N. (2013). Body fat distribution and insulin resistance. Nutrients 5, 2019–2027. doi: 10.3390/nu5062019
Raji, A., Seely, E. W., Arky, R. A., and Simonson, D. C. (2001). Body fat distribution and insulin resistance in healthy Asian Indians and Caucasians. J. Clin. Endocrinol. Metab. 86, 5366–5371. doi: 10.1210/jcem.86.11.7992
Razak, F., Anand, S. S., Shannon, H., Vuksan, V., Davis, B., Jacobs, R., et al. (2007). Defining obesity cut points in a multiethnic population. Circulation 115, 2111–2118. doi: 10.1161/circulationaha.106.635011
Rees, J. M. B., Wood, A. M., and Burgess, S. (2017). Extending the MR-Egger method for multivariable Mendelian randomization to correct for both measured and unmeasured pleiotropy. Stat. Med. 36, 4705–4718. doi: 10.1002/sim.7492
Robiou-du-Pont, S., Bonnefond, A., Yengo, L., Vaillant, E., Lobbens, S., Durand, E., et al. (2013). Contribution of 24 obesity-associated genetic variants to insulin resistance, pancreatic beta-cell function and type 2 diabetes risk in the French population. Int. J. Obes. 37, 980–985. doi: 10.1038/ijo.2012.175
Roglic, G. (2016). WHO Global report on diabetes: a summary. Int. J. Noncommun. Dis. 1:3. doi: 10.4103/2468-8827.184853
Sampath Kumar, A., Maiya, A. G., Shastry, B. A., Vaishali, K., Ravishankar, N., Hazari, A., et al. (2019). Exercise and insulin resistance in type 2 diabetes mellitus: A systematic review and meta-analysis. Ann. Phys. Rehabil. Med. 62, 98–103. doi: 10.1016/j.rehab.2018.11.001
Smit, R. A. J., Trompet, S., Leong, A., Goodarzi, M. O., Postmus, I., Warren, H., et al. (2019). Statin-induced LDL cholesterol response and type 2 diabetes: a bidirectional two-sample Mendelian randomization study. Pharmacogenom. J. 20, 462–470. doi: 10.1038/s41397-019-0125-x
Smith, G. D., and Ebrahim, S. (2003). ‘Mendelian randomization’: can genetic epidemiology contribute to understanding environmental determinants of disease? Int. J. Epidemiol. 32, 1–22. doi: 10.1093/ije/dyg070
Vazquez, G., Duval, S., Jacobs, D. R. Jr., and Silventoinen, K. (2007). Comparison of body mass index, waist circumference, and waist/hip ratio in predicting incident diabetes: a meta-analysis. Epidemiol. Rev. 29, 115–128. doi: 10.1093/epirev/mxm008
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, 693–698. doi: 10.1038/s41588-018-0099-7
Wainberg, M., Mahajan, A., Kundaje, A., McCarthy, M. I., Ingelsson, E., Sinnott-Armstrong, N., et al. (2019). Homogeneity in the association of body mass index with type 2 diabetes across the UK Biobank: A Mendelian randomization study. PLoS Med. 16:e1002982. doi: 10.1371/journal.pmed.1002982
Wang, T., Zhang, R., Ma, X., Wang, S., He, Z., Huang, Y., et al. (2018). Causal Association of Overall Obesity and Abdominal Obesity with Type 2 Diabetes: A Mendelian Randomization Analysis. Obesity 26, 934–942. doi: 10.1002/oby.22167
Weng, L. C., Roetker, N. S., Lutsey, P. L., Alonso, A., Guan, W., Pankow, J. S., et al. (2018). Evaluation of the relationship between plasma lipids and abdominal aortic aneurysm: A Mendelian randomization study. PLoS One 13:e0195719. doi: 10.1371/journal.pone.0195719
Wheeler, E., Leong, A., Liu, C. T., Hivert, M. F., Strawbridge, R. J., Podmore, C., et al. (2017). Impact of common genetic determinants of Hemoglobin A1c on type 2 diabetes risk and diagnosis in ancestrally diverse populations: A transethnic genome-wide meta-analysis. PLoS Med. 14:e1002383. doi: 10.1371/journal.pmed.1002383
Xu, L., and Hao, Y. T. (2017). Effect of handgrip on coronary artery disease and myocardial infarction: a Mendelian randomization study. Sci. Rep. 7:954. doi: 10.1038/s41598-017-01073-z
Xue, A., Wu, Y., Zhu, Z., Zhang, F., Kemper, K. E., Zheng, Z., et al. (2018). Genome-wide association analyses identify 143 risk variants and putative regulatory mechanisms for type 2 diabetes. Nat. Commun. 9:2941. doi: 10.1038/s41467-018-04951-w
Yang, Q., Millard, L. A. C., and Davey Smith, G. (2019). Proxy gene-by-environment Mendelian randomization study confirms a causal effect of maternal smoking on offspring birthweight, but little evidence of long-term influences on offspring health. Int. J. Epidemiol. doi: 10.1093/ije/dyz250 [Epub ahead of print]
Zheng, J., Baird, D., Borges, M. C., Bowden, J., Hemani, G., Haycock, P., et al. (2017). Recent developments in mendelian randomization studies. Curr. Epidemiol. Rep. 4, 330–345. doi: 10.1007/s40471-017-0128-6
Keywords: type 2 diabetes mellitus, insulin resistance, abdominal obesity, body fat mass, body fat mass distribution, Mendelian randomization
Citation: Xu H, Jin C and Guan Q (2020) Causal Effects of Overall and Abdominal Obesity on Insulin Resistance and the Risk of Type 2 Diabetes Mellitus: A Two-Sample Mendelian Randomization Study. Front. Genet. 11:603. doi: 10.3389/fgene.2020.00603
Received: 26 January 2020; Accepted: 18 May 2020;
Published: 02 July 2020.
Edited by:
Lide Han, Vanderbilt University Medical Center, United StatesReviewed by:
Dan Zhou, Vanderbilt University Medical Center, United StatesPing Zeng, Xuzhou Medical University, China
Copyright © 2020 Xu, Jin and Guan. 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: Qingbo Guan, ZG9jdG9yZ3VhbnFpbmdib0AxNjMuY29t
†These authors have contributed equally to this work