Next Article in Journal
Incorporation of Suppression of Tumorigenicity 2 into Random Survival Forests for Enhancing Prediction of Short-Term Prognosis in Community-ACQUIRED Pneumonia
Previous Article in Journal
Osteoarthritis: New Insight on Its Pathophysiology
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
JCM
Volume 11
Issue 20
10.3390/jcm11206014
jcm-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Effect of Body Fat Distribution on Systemic Sclerosis

by
Gonzalo Villanueva-Martin
1,*,
Marialbert Acosta-Herrera
2,
Martin Kerick
1,
Elena López-Isac
1,
Carmen P. Simeón
3,
José L. Callejas
4,
Shervin Assassi
5,
Lorenzo Beretta
6,
International SSc Group
†,
Australian Scleroderma Interest Group (ASIG)
†,
Yannick Allanore
7,
Susanna M. Proudman
8,
Mandana Nikpour
8,
Carmen Fonseca
9,
Christopher P. Denton
9,
Timothy R. D. J. Radstake
10,
Maureen D. Mayes
5,
Xia Jiang
11,12,13,
Javier Martin
1,‡ and
Lara Bossini-Castillo
14,15,*,‡
1
Department of Cell Biology and Immunology, Institute of Parasitology and Biomedicine López-Neyra, CSIC, 18016 Granada, Spain
2
Systemic Autoimmune Disease Unit, Hospital Clínico San Cecilio, Instituto de Investigación Biosanitaria Ibs. GRANADA, 18016 Granada, Spain
3
Department of Internal Medicine, Valle de Hebrón Hospital, 08035 Barcelona, Spain
4
Department of Internal Medicine, Hospital San Cecilio, 18016 Granada, Spain
5
Department of Rheumatology, The University of Texas Health Science Center, Houston, TX 77030, USA
6
Referral Center for Systemic Autoimmune Diseases, Fondazione IRCCS Ca’ Granda Ospedale Maggiore Policlinico di Milano, 20122 Milan, Italy
7
Department of Rheumatology A, Hospital Cochin, 75005 Paris, Île-de-France, France
8
Department of Medicine, The University of Melbourne at St. Vincent’s Hospital, Melbourne, VIC 3010, Australia
9
Center for Rheumatology, Royal Free and University College Medical School, London 41125, UK
10
Department of Rheumatology and Clinical Immunology, University Medical Center Utrecht, 3584 CX Utrecht, The Netherlands
11
Department of Clinical Neuroscience, Center for Molecular Medicine, Karolinska Institutet, 17176 Stockholm, Sweden
12
Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, MA 02115, USA
13
West China School of Public Health, West China Fourth Hospital, Sichuan University, Chengdu 610041, China
14
Departamento de Genética e Instituto de Biotecnología, Centro de Investigación Biomédica (CIBM), Universidad de Granada, 18100 Granada, Spain
15
Advanced Therapies and Biomedical Technologies (TEC-14), Instituto de Investigación Biosanitaria Ibs. GRANADA, 18016 Granada, Spain
 
add Show full affiliation list
 
remove Hide full affiliation list
*
Authors to whom correspondence should be addressed.
†
See Supplementary Materials.
‡
Authors share senior authorship.
J. Clin. Med. 2022, 11(20), 6014; https://doi.org/10.3390/jcm11206014
Submission received: 16 September 2022 / Revised: 6 October 2022 / Accepted: 10 October 2022 / Published: 12 October 2022
(This article belongs to the Section Immunology)
Browse Figures
Versions Notes

Abstract

:
Obesity contributes to a chronic proinflammatory state, which is a known risk factor to develop immune-mediated diseases. However, its role in systemic sclerosis (SSc) remains to be elucidated. Therefore, we conducted a two-sample mendelian randomization (2SMR) study to analyze the effect of three body fat distribution parameters in SSc. As instrumental variables, we used the allele effects described for single nucleotide polymorphisms (SNPs) in different genome-wide association studies (GWAS) for SSc, body mass index (BMI), waist-to-hip ratio (WHR) and WHR adjusted for BMI (WHRadjBMI). We performed local (pHESS) and genome-wide (LDSC) genetic correlation analyses between each of the traits and SSc and we applied several Mendelian randomization (MR) methods (i.e., random effects inverse-variance weight, MR-Egger regression, MR pleiotropy residual sum and outlier method and a multivariable model). Our results show no genetic correlation or causal relationship between any of these traits and SSc. Nevertheless, we observed a negative causal association between WHRadjBMI and SSc, which might be due to the effect of gastrointestinal complications suffered by the majority of SSc patients. In conclusion, reverse causality might be an especially difficult confounding factor to define the effect of obesity in the onset of SSc.

1. Introduction

Systemic sclerosis (SSc) is an immune-mediated disease (IMD), characterized by abnormal immunological activation, vascular damage and fibrosis of the skin [1]. SSc represents a major challenge for clinicians, as it has a deep impact on the life quality and life expectancy of the affected patients [1]. Recent efforts in the study of the genetic factors that contribute to the onset and progression of SSc, such as several large-scale genetic association studies and genome-wide association studies (GWAS) [2], have contributed to identifying genetic susceptibility markers both in the human leukocyte antigen (HLA) locus and outside this highly polymorphic region [3]. The largest GWAS to date comprised more than 9000 patients with SSc and allowed the identification of 19 non-HLA loci associated with the disease [2]. Moreover, recent studies have identified specific HLA-DQA1 alleles exclusively associated with different clinical subtypes of SSc [3]. Therefore, the number of relevant loci that have been firmly associated with this condition has remarkably increased over the last decade. Although the use of genetic risk factors to predict the risk of developing SSc was explored in a recent genomic risk score (GRS) [4], the involvement of these genetic risk factors in the disease pathogenesis and the affected biological pathways have not been fully established yet [5].
Despite the advances in the identification of the genetic factors contributing to the heritability of SSc, the complex nature of this disorder is an intrinsic obstacle to studying the pathological mechanisms that lead to the disruption of the immune homeostasis and to the onset of fibrotic processes in affected individuals. Well-established environmental triggers for SSc are silica and solvents, extreme or long-term exposure to which is related to the disease’s development [6,7]. Moreover, demographic and clinical characteristics, such as sex, age, ethnical origin, hormone levels, etc., have been pointed out as risk factors for SSc [6,8]. However, the roles of lifestyle and environmental triggers in the manifestation and prognosis of SSc are still elusive.
Mendelian randomization (MR) uses SNPs as instrumental variants (IVs) in order to determine if they are acting on a disease or outcome through a risk factor or exposure [9,10]. The principle of the method is that alleles are randomly distributed during gametogenesis and their presence pre-exists the disease. These genetic facts mimic the random distribution of clinical trials and take away the causality of the disease on the variable, reducing confounding factors [11]. For a genetic variant to be considered an IV, it is assumed that it is associated with exposure. However, an IV cannot be associated with any confounding factor related to the risk factor or the outcome, either directly or indirectly. Additionally, the effects of the IV on the outcome should only be mediated by the exposure [9]. Therefore, only when genetic polymorphisms are relevant, independent and have a restricted effect on the outcome, can they be considered IVs. In a classical MR study, the allele effects on outcome and exposure are obtained from the same individuals [9,10]. However, detailed information for multiple traits is difficult to obtain in a large population. Two-sample MR (2SMR) methods allow us to combine the estimations of the IV allele effects relying only on GWAS summary statistics for the outcome and for the exposure from independent studies. The implementation of these methods has improved the statistical power to detect causal associations between risk factors and disease, which has shown promising results in several conditions [12].
Obesity-related diseases are becoming a public health issue in Western countries [13], since obesity rates are increasing due to unhealthy lifestyles. Obesity is defined by an excess of fat in the body and body fat distribution can be measured by a variety of methods, for instance body mass index (BMI) and waist to hip ratio (WHR). BMI is the most common body fat proxy and it is the gold standard for obesity. BMI is measured as the body weight normalized by height squared (kg/m2) [14], and it is known that BMI > 25 kg/m2 is associated with an increased risk of suffering from chronic diseases such as cardiovascular disease, type II diabetes or specific cancers [15]. Nevertheless, BMI has certain limitations, and anthropometric measures of abdominal obesity, such as WHR, seem to be better indicators of excessive fat mass [16]. Since WHR measures both visceral and gluteal fat, it stands out among other anthropometric traits [17]. If WHR is adjusted for BMI (WHRadjBMI), it is possible to obtain an anthropometric measure which is independent from the overall adiposity, and to combine the most standardized measure of obesity and the anthropometric measure that best captures the distribution of body fat [17,18]. Taking advantage of the publicly available GWAS results, MR approaches have been successful in identifying risk factors for IMDs, such as obesity-related traits [19,20]. Excess of fat has been associated with a low but persistent proinflammatory state that is believed to promote IMDs [13,21]. However, in the case of SSc, the relationship between body fat distribution and SSc remains to be explored.
Consequently, in order to analyze the effect of nutritional status on SSc risk, we applied the novel 2SMR methods on the largest GWAS of SSc patients [2] with European ancestry and the biggest GWAS meta-analysis for anthropometric fat distribution measures to date [22].

2. Materials and Methods

2.1. Instrumental Variables

The study design of the 2SMR study of SSc and three obesity-related traits is summarized in Figure 1. The outcome instrumental variables (IV-outcome), i.e., the selected genetic variants and their effect sizes in SSc, were obtained from the largest SSc GWAS meta-analysis, which included 9846 SSc patients and 18,333 healthy controls from 14 different cohorts with European ancestry [2]. Patient demographic data are shown in Table S1. Additionally, SNP effect sizes after stratification by sex, serological and clinical subtype, as reported elsewhere [3], were also analyzed. Finally, we performed sex-specific analyses, including only either the female or the male individuals from the different cohorts and following the previously described analysis framework [2].
In the case of the exposures, we obtained the IVs (IV-exposure) from a recent GWAS meta-analysis between the cohorts included in the Genetic Investigation of Anthropometric Traits consortium (GIANT) project and those recruited for the UK Biobank (UKBB) repository for different anthropometric measures [23]. We only used the summary statistics comprising individuals with European ancestry, which included 806,810 individuals and 27,381,302 SNPs for BMI, a classical obesity parameter, and for two parameters that assess body fat distribution: WHR, comprising 697,734 individuals and 27,376,273 SNPs, and WHRadjBMI, covering 694,649 individuals and 27,375,636 SNPs [23]. None of the participants recruited in the SSc studies overlapped with the exposure GWASs to the best of our knowledge.

2.2. Genomic Association Analysis

Genetic correlation. To determine causality between obesity risk factors and SSc, we calculated the total genomic correlation between them. First, we performed an approximation implemented in the linkage disequilibrium regression score (LDSC) software [24]. Then, to study the contribution of specific regions (pairwise local genetic correlation), we used the methods supported by ρ-HESS software (version 0.5.3) [25]. Briefly, ρ-HESS software splits the genome into 1703 small regions through the chromosomes and uses LD matrices to create eigenvectors and to project the GWAS effect sizes. Then, the local SNP heritability per trait is calculated and, finally, genetic covariance between traits is estimated. We adjusted our significance thresholds for multiple testing, i.e., 1.1 × 10−3 (0.05/45) for LDSC and 2.9 × 10−5 (0.05/1703) for ρ-HESS.
Mendelian randomization analysis. In order to assess if there was a causal relationship between body fat distribution and SSc or any of the stratified sets of patients, we performed a 2SMR study implementing the R package “TwoSampleMR” [10]. Considering the complex linkage disequilibrium (LD) patterns and the strong genetic associations described in the HLA locus SSc [2,3,26], the extended HLA region (chromosome 6: 20,000,000–40,000,000 bp) was excluded from the MR analyses in order to prevent biases.
The selected IVs were based on the original independent signal analysis reported by Pulit et al. [23]. Briefly, the independent signals from results from the inverse variance meta-analysis (p < 5 × 10−9) were identified by LD-based clumping (r2 > 0.05 and ±5 Mb). Secondary signals were also defined by conditional analyses (p < 5 × 10−9) and locus LD clumping. We extracted the association estimates for these SNPs or the best available proxy (according to the LD patterns observed in the UKBB cohort) that was present in the SSc dataset. The number of shared SNPs between SSc and the exposures reached 533, 247 and 262 for BMI, WHR and WHRadjBMI, respectively (Table S2).
Three gold-standard 2SMR methods were selected. A random effects inverse-variance weight (IVW) approach was taken, which pools the effects of each IV and balances the global pleiotropy to zero by assuming the validity or invalidity of all the SNPs [10]. An MR-Egger regression method [27] was applied, which is able to estimate causality even when all IVs are weak or invalid and to calculate horizontal pleiotropy. Although these methods are very robust for MR analysis, both of them have limitations in dealing with outlier IVs. For that reason, we also applied the MR pleiotropy residual sum and outlier (MR-PRESSO) method [28]. The MR-PRESSO algorithm detects outlier IVs that exert horizontal pleiotropy in a multi-instrument Mendelian randomization analysis. Moreover, MR-PRESSO provides outlier-free causality estimates.
Additionally, to estimate the effect of the IVs controlling for their effect on other exposures, we performed a multivariable Mendelian randomization analysis (MVMR), implemented in the TwoSampleMR package [29]. This analysis included a set of unique LD-clumped IV exposures for both BMI and WHR, which were regressed against SSc together, weighting for the inverse variance of SSc for these IVs.
The Benjamini–Hochberg false discovery rate (FDR) correction was applied, and we considered p < 0.05 as significant [30].

2.3. Sensitivity Analysis

The statistical power of our analyses was calculated using the algorithm described by Brion et al. for MR studies [31]. Aiming to control for the effect of potential confounding factors, we removed any of the SNPs with reported associations with known obesity-related confounding factors (Table S3) from the MR analysis, as reported by the GWAS catalog [32], SNPnexus [33] and ClinVar [34]. We studied the contribution of each SNP to the observed effects by carrying out a leave-one-out sensitivity analysis, implemented in the “TwoSampleMR” package [10]. By these means, we observed that the exclusion of one SNP at a time did not affect the observed results.

3. Results

Leveraging Mendelian randomization as a novel methodological strategy, we studied for the first time the causal contribution of body fat distribution to the risk of suffering from SSc (Figure 1). Here we used the GWAS summary statistics of the largest SSc meta-analysis [2] as an outcome and three obesity-related traits GWAS comprising thousands of European ancestry individuals as exposures.

3.1. Genomic Correlation. Only the HLA Locus Harbours Local Genetic Correlation between SSc and Body Fat Distribution

At a genomic scale, we observed a strong genome-wide correlation between BMI and WHR (rg = 0.59, [95% CI −0.016–0.051]) and between WHR and WHRadjBMI (rg = 0.78, [95% CI −0.01–0.03]), but not between WHRadjBMI and BMI (rg = −4.02 × 10−2, [95% CI −0.016–0.049]), as previously described [19] (Figure 2). However, our results show no evidence of correlation between SSc and the three tested obesity-related traits (BMI rg = −0.039 [95% CI −0.033–0.102]; WHR rg = −0.054, [95% CI −0.035–0.106]; WHRadjBMI rg = −0.041, [95% CI −0.04–0.122], all observed p > 0.05) (Figure 2).
Even when there is no correlation between traits at a genome-wide level, it is possible that the traits show local correlation at specific loci. To address this potential correlation, we performed a local genetic correlation analysis between BMI, WHR, WHRadjBMI and SSc (Figure A1). The local correlation observed in these regions reached local-rg = 8.5 × 10−4 and local-rg = 2.6 × 10−4 (Figure A1).

3.2. The Analysis of the Causal Relationship between Obesity-Related Traits and Systemic Sclerosis Is Limited by Confounding Factors

Despite the limited genetic correlation found, we explored the possible causal relationship between body fat distribution and SSc. Considering the complex LD patterns in the HLA regions and the local genetic correlation found only in this locus, it was excluded from the following MR analyses. The available SSc dataset was powerful enough to detect associations of 25% increased risk of SSc with BMI (99%), WHR (83%) and WHRadjBMI (92%) (Table S4), considering an explained phenotypic variance of 2.5–5% and the complete set comprising 28,179 individuals (34.9% cases). We were confident about the statistical power estimated for the largest subsets of patients, for instance, females (BMI power = 79%, WHR power = 82% and WHRadjBMI = 87%), lcSSc (BMI power = 94%, WHR power = 70% and WHRadjBMI = 81%) and ACA+ (BMI power = 83%, WHR power = 53% and WHRadjBMI = 65%). However, the analyses for the less frequent patient groups, i.e., males (BMI power = 30%, WHR power = 8% and WHRadjBMI = 10%) and ATA+ (BMI power = 14%, WHR power = 10% and WHRadjBMI = 12%) were clearly insufficient to identify true causal relationships (Table S4).
As reported in Table 1 and Table S5, classical MR methods showed no significant evidence of causality for BMI or WHR on SSc, whether including only the index SNPs or considering both the index SNPs and the secondary signals. The results for BMI under the random effects IVW model show a suggestive positive association with BMI, but this association did not reach statistical significance (OR under random effects IVW = 1.15 [95% CI 0.67–1.98]). A trend of negative association considering index and secondary signals was only observed in the case of the random effects IVW model for WHR (Table 1). All the remaining models showed p > 0.05 and the ORs were in the range of 0.93–1.15 for BMI and 0.27–0.82 for WHR. In the case of WHRadjBMI (WHR after regressing out the effect of BMI), a negative association with SSc reached statistical significance in the three tested models (OR under random effects IVW = 0.73 [95% CI 0.56–0.94], MR-Egger = 0.43 [95% CI 0.20–0.90], MR-PRESSO = 0.77 [95% CI 0.60–0.99]). These associations with WHRadjBMI remained negative in the analyses that included only index signals, but only the MR-Egger model was significant after multiple testing correction (OR under MR-Egger = 0.69 [95% CI 0.51–0.93], see Table S5).
We carried out a sensitivity analysis, which implied the removal of SNPs associated with known obesity-related confounders (Table S3), to address the effect of these confounders on the lack of significance for the BMI models and the negative relationships with WHR and WHRadjBMI. As shown in Table 2 and Table S6, the confounder-free models did not change the observed negative relationship and none of them reached a significant result after FDR correction. Although we observed effect size heterogeneity for the different genetic variants (Table S7), the analyses of the intercept parameter in the MR-Egger models did not reveal any signs of horizontal pleiotropy, and the effects were not affected by the removal of the outlier SNPs identified by the MR-PRESSO algorithm (Table 1 and Table 2, Tables S5 and S6). Furthermore, leave-one-out analyses did not highlight that these effects were influenced only by one variant (Figure A2).
We decided to implement an MVMR model, considering the significant associations observed for WHRadjBMI and the limitations of the univariate models, to test for the combined influence of several exposures and to control for the effect of confounding factors. This analysis allowed us to directly test the association of BMI and WHR with SSc, controlling for the effects of both parameters at the same time. As expected, the results of these analyses show an effect for WHR (MVMR OR 0.80 [95% CI 0.57–1.13]) that is similar to the previously identified effect for WHRadjBMI (Table 3). Nevertheless, no significant association of BMI with SSc was revealed (MVMR OR 1.03 [95% CI 0.79–1.33]) (Table 3). These findings might point towards a negative or inexistent effect of WHR in SSc and, if any, a very modest risk effect for BMI.
Considering the well-known clinical and genetic differences between the SSc subsets of patients [35], we explored subset-specific effects for the selected exposures. Several associations remained significant in the stratified analyses, especially in the largest and more powerful subsets, such as lcSSc (Table S8). However, the direction and magnitude of the exposure effects were consistent in all the subsets (Table S8), which suggested a uniform effect, if any, in all the patients. There were no significant differences between the models with and without the secondary signals (Table S8). Moreover, taking into account the higher frequency of SSc in females (9 female: 1 male ratio) [8], we performed sex-specific analyses too. In these analyses, we relied on female-only and male-only GWAS summary statistics for both SSc and the obesity-related risk factors. Once more, although the risk effect of BMI, WHR and WHRadjBMI seemed more evident in men, these effects did not reach statistical significance (Table S8).

4. Discussion

This report addressed the risk effect of body fat distribution in SSc for the first time. We exhaustively exploited public GWAS summary statistics for both SSc and for anthropometric traits and the development of novel MR methods. We did not observe global genomic correlation between the outcome and any of the exposures. Moreover, local genetic correlation was only found in the HLA locus, a highly complex region. Different MR methods were then applied to identify possible causal relationships between the obesity traits and SSc. However, no significant causal risk effect of the exposures was found in this case.
Although our results do not support the causal relation between exposures and outcome, it should be noted that the statistical power of the SSc dataset is modest compared to similar studies performed to date in other IMDs, such as RA or IBD [36] (Table S2). SSc is a rare IMD and, despite recent advances [1,2,26], the recruitment of large patient cohorts remains challenging. Therefore, future efforts to enlarge the size or to complement the available SSc GWAS information might help to identify causal risk factors.
We found that the effect of confounders might be more severe in the case of SSc than in other IMDs. Gastrointestinal involvement (GI), which affects more than 70% of SSc patients [36], hinders food ingestion, and patients are mostly thin [37]. In fact, weight loss has been used as one of the SSc diagnostic markers [35]. This direct effect of the onset symptoms in the exposures is known as reverse causality, and it is a remarkably difficult confounding factor to control for [38]. Reverse causality might be the cause behind both the lack of significant risk effects of BMI in SSc and the reported negative relationship between WHR and SSc, which becomes more evident when the effect of BMI is subtracted in the analysis of WHRadjBMI (Table 1 and Table 2, Table S5 and Table S6).
Bad diet habits and obesity are associated with an increased risk to suffer from IMDs such as RA and IBD [19,20,39]. Higher BMI has been associated with increased risk to Crohn’s disease (CD) and Rheumatoid Arthritis (RA), but negative associations with BMI have been reported for ulcerative colitis (UC), and a recent study found reverse causality between WHR and RA [19,20,39]. IMDs are often present as comorbidities and share altered molecular pathways, environmental triggers and genetic risk factors [40]. Furthermore, the role of adipocytes in the activation of the immune system is prominent, especially due to the release of adipokines [41]. Adipokines are molecules known to be involved in the “obesity–autoimmunity” relationship [13,42], such as lectins or cytokines, especially adiponectin, but also interleukins and tumor necrosis factor alpha (TNFɑ) [13]. Interestingly, patients with SSc and a high BMI have been shown to have higher lectin levels than healthy controls [43], and it has been established that subcutaneous adipocytes can act as progenitor cells for fibroblasts [44,45]. These fibroblasts may eventually transdifferentiate into myofibroblasts [46], activated profibrotic fibroblasts that are characteristic of the fibrotic lesions observed in SSc patients, and recent evidence has shown that the activation of adipocyte-derived mesenchymal cells from SSc skin biopsies to myofibroblasts is possible using soluble molecules present in the skin microenvironment in SSc [47].
In order to rule out the role of obesity as a risk factor for SSc, body fat distribution measures from the patients before the onset of GI or BMI-matched case-control sets would be very valuable resources.
The negative association that is observed for WHR might be due to additional confounding factors that are inherent to SSc and that affect body fat distribution, for example, sex or lipid profiles [15]. Remarkably, WHR is different in women than in men and there is a clear sex bias in SSc [35]. Therefore, we hypothesized that there could be a sex-specific association and performed stratified analyses with the female and male cohorts separately. Our results show significant causal associations with SSc only in females, but considering the statistical power differences and the similarity between the effect sizes, the lack of significance for the male group may be likely due to the reduced sample size (Table S4). The key role of sample size as a limitation of our study to identify weak risk effects was also clear in other stratified analyses, as we found consistent ORs for all the tested clinical subtypes of SSc patients, but the models reached statistical significance only in the largest subsets (Table S8).
In conclusion, this study found no significant evidence that supported the role of body fat distribution as a causal risk factor for SSc using 2SMR methods. Nevertheless, the current GWAS have a limited statistical power to identify modest contributions to SSc risk and the intrinsic nature of the SSc clinical complications might be acting as potential constraints in this study. Consequently, further analyses are needed to rule out the role of obesity in the onset of Ssc.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/jcm11206014/s1. Table S1. Demographic data from SSc meta GWAS. Table S2. GWAS selection summary. Table S3. List of discarded SNPs due to association with known obesity-related confounders. Table S4. Statistical power calculation for each stratification and outcome-exposure combination. Table S5. Association between genetically-predicted obesity-related traits and risk of SSc. Analysis including index and secondary signals for the obesity-related traits and excluding the HLA region. Table S6. Association between genetically-predicted obesity-related traits and risk of SSc. Analysis including index signals for the obesity-related traits and excluding the HLA region and known obesity-related confounder SNPs. Table S7. Effect sizes in the exposures and in the outcome of the BMI, WHR and WHRadjBMI index and secondary SNPs. Table S8. Association between genetically-predicted obesity-related traits and risk of SSc. Analysis including index and secondary signals for the obesity-related traits and excluding the HLA region and known obesity-related confounder SNPs. SSc stratifications comprise sex (females and males), clinical subtypes (lcSSc and dcSSc) and serological subsets (ACA+ and ATA+). File S1: Collaborators.

Author Contributions

G.V.-M.: data analysis, manuscript drafting, revision and approval; M.A.-H.: data interpretation, manuscript revision and approval; M.K.: data interpretation, manuscript revision and approval; E.L.-I.: data interpretation, manuscript revision and approval; C.P.S.: data acquisition, manuscript revision and approval; J.L.C.: data acquisition; S.A.: data acquisition; L.B.: data acquisition, manuscript revision and approval; International SSc Group: data acquisition; Australian Scleroderma Interest Group (ASIG): data acquisition; Y.A.: data acquisition, manuscript revision and approval; S.M.P.: data acquisition, manuscript revision and approval; M.N.: data acquisition; C.F.: data acquisition, manuscript revision and approval, C.P.D.: data acquisition, manuscript revision and approval, T.R.D.J.R.: data acquisition, manuscript revision and approval, M.D.M.: data acquisition, manuscript revision and approval; X.J.: study design, data interpretation, manuscript revision and approval; J.M.: study design, manuscript drafting, revision and approval; L.B.-C.: study design, manuscript drafting, revision and approval. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by grant RTI2018101332-B-100 funded by MCIN/AEI/10.13039/501100011033 and by “ERDF A way of making Europe” funded by the European Union. Red de Investigación en Inflamación y Enfermedades Reumáticas (RIER) from Instituto de Salud Carlos III (RD16/0012/0013). 115565. LB-C was funded by Grant IJC2018-038026-I funded by MCIN/AEI/10.13039/501100011033. MA-H is a recipient of a Miguel Servet fellowship (CP21/00132) from Instituto de Salud Carlos III (Spanish Ministry of Science and Innovation). EL-I was funded by Grant IJC2019-040080-I funded by MCIN/AEI/10.13039/501100011033. GV-M was funded by Grant PRE2019-087586 funded by MCIN/AEI/10.13039/501100011033 and by “ESF Investing in your future”.

Institutional Review Board Statement

An ethical protocol was prepared, consensus was reached across all partners, academic and industrial, translated into all participants’ languages and approved by each of the local ethical committees of the clinical recruitment centers. The studies adhered to the standards set by the International Conference on Harmonization and Good Clinical Practice (ICH-GCP), and to the ethical principles that have their origin in the Declaration of Helsinki (2013). The protection of the confidentiality of records that could identify the included subjects is ensured as defined by the EU Directive 2001/20/EC and the applicable national and international requirements relating to data protection in each participating country. For the systemic sclerosis data, as previously described (Acosta-Herrera M, et al. Comprehensive analysis of the major histocompatibility complex in systemic sclerosis identifies differential HLA associations by clinical and serological subtypes. Ann Rheum Dis. 2021. https://doi.org/10.1136/annrheumdis-2021-219884) (López-Isac E, et al. GWAS for systemic sclerosis identifies multiple risk loci and highlights fibrotic and vasculopathy pathways. Nat Commun. 2019. https://doi.org/10.1038/s41467-019-12760-y). CSIC’s ethics committee approved the study protocol, and written informed consent was obtained in accordance with the tenets of the Declaration of Helsinki. In the case of the obesity-related traits, GWAS summary statistics data were public at the time of the study and accessible as described in the data accessibility section.

Informed Consent Statement

Not applicable.

Acknowledgments

This research is part of the doctoral degree awarded to G.V.-M., within the Biomedicine program from the University of Granada, entitled “Deciphering the genetic basis of systemic sclerosis”. We would like to thank Matthew Brown, from the NIHR Biomedical Research Centre, Guy’s and Saint Thomas’ NHS Foundation Trust and King’s College, London, UK, for his contributions to sample and data acquisition. We would like to thank Sofia Vargas and Gema Robledo for their excellent technical assistance. We also appreciate the controls and the affected individuals who generously provided the samples for these studies.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Jcm 11 06014 g0a1 550
Figure A1. Local genetic correlation, local genetic variance and local SNP heritability between SSc and (a) BMI, (b) WHR and (c) WHRadjBMI.
Figure A1. Local genetic correlation, local genetic variance and local SNP heritability between SSc and (a) BMI, (b) WHR and (c) WHRadjBMI.
Jcm 11 06014 g0a1
Jcm 11 06014 g0a2 550
Figure A2. MR leave-one-out sensitivity analyses for: (A) BMI, (B) WHR and (C) WHRadjBMI.
Figure A2. MR leave-one-out sensitivity analyses for: (A) BMI, (B) WHR and (C) WHRadjBMI.
Jcm 11 06014 g0a2

References

  1. Radstake, T.R.D.J.; Gorlova, O.; Rueda, B.; Martin, J.-E.; Alizadeh, B.Z.; Palomino-Morales, R.; Coenen, M.J.; Vonk, M.C.; Voskuyl, A.E.; Schuerwegh, A.J.; et al. Genome-wide association study of systemic sclerosis identifies CD247 as a new susceptibility locus. Nat. Genet. 2010, 42, 426–429. [Google Scholar] [CrossRef] [Green Version]
  2. López-Isac, E.; Acosta-Herrera, M.; Kerick, M.; Assassi, S.; Satpathy, A.T.; Granja, J.; Mumbach, M.R.; Beretta, L.; Simeón, C.P.; Carreira, P.; et al. GWAS for systemic sclerosis identifies multiple risk loci and highlights fibrotic and vasculopathy pathways. Nat. Commun. 2019, 10, 4955. [Google Scholar] [CrossRef] [Green Version]
  3. Acosta-Herrera, M.; Kerick, M.; Lopéz-Isac, E.; Assassi, S.; Beretta, L.; Simeón-Aznar, C.P.; Ortego-Centeno, N.; International SSc Group; Proudman, S.M.; Australian Scleroderma Interest Group (ASIG); et al. Comprehensive analysis of the major histocompatibility complex in systemic sclerosis identifies differential HLA associations by clinical and serological subtypes. Ann. Rheum. Dis. 2021, 80, 1040–1047. [Google Scholar] [CrossRef]
  4. Bossini-Castillo, L.; Villanueva-Martin, G.; Kerick, M.; Acosta-Herrera, M.; López-Isac, E.; Simeón, C.P.; Ortego-Centeno, N.; Assassi, S.; International SSc Group; Australian Scleroderma Interest Group (ASIG); et al. Genomic Risk Score Impact on Susceptibility to Systemic Sclerosis. Ann. Rheum. Dis. 2021, 80, 118–127. [Google Scholar] [CrossRef]
  5. Bossini-Castillo, L.; López-Isac, E.; Mayes, M.D.; Martín, J. Genetics of systemic sclerosis. Semin. Immunopathol. 2015, 37, 443–451. [Google Scholar] [CrossRef]
  6. Barnes, J.; Mayes, M.D. Epidemiology of Systemic Sclerosis: Incidence, Prevalence, Survival, Risk Factors, Malignancy, and Environmental Triggers. Curr. Opin. Rheumatol. 2012, 24, 165–170. [Google Scholar] [CrossRef] [PubMed]
  7. Calderon, L.M.; Pope, J.E. Scleroderma epidemiology update. Curr. Opin. Rheumatol. 2021, 33, 122–127. [Google Scholar] [CrossRef] [PubMed]
  8. Allanore, Y.; Simms, R.; Distler, O.; Trojanowska, M.; Pope, J.; Denton, C.P.; Varga, J. Systemic Sclerosis. Nat. Rev. Dis. Primers 2015, 1, 15002. [Google Scholar] [CrossRef]
  9. Burgess, S.; Butterworth, A.; Thompson, S.G. Mendelian Randomization Analysis with Multiple Genetic Variants Using Summarized Data. Genet. Epidemiol. 2013, 37, 658–665. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Bowden, J.; Davey Smith, G.; Burgess, S. Mendelian randomization with invalid instruments: Effect estimation and bias detection through Egger regression. Int. J. Epidemiol. 2015, 44, 512–525. [Google Scholar] [CrossRef]
  11. Smith, G.D.; Ebrahim, S. “Mendelian Randomization”: Can Genetic Epidemiology Contribute to Understanding Environmental Determinants of Disease? Int. J. Epidemiol. 2003, 32, 1–22. [Google Scholar] [CrossRef] [Green Version]
  12. Hartwig, F.P.; Davies, N.M.; Hemani, G.; Smith, G.D. Two-sample Mendelian randomization: Avoiding the downsides of a powerful, widely applicable but potentially fallible technique. Int. J. Epidemiol. 2016, 45, 1717–1726. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Versini, M.; Jeandel, P.-Y.; Rosenthal, E.; Shoenfeld, Y. Obesity in autoimmune diseases: Not a passive bystander. Autoimmun. Rev. 2014, 13, 981–1000. [Google Scholar] [CrossRef] [PubMed]
  14. Borga, M.; West, J.; Bell, J.D.; Harvey, N.C.; Romu, T.; Heymsfield, S.B.; Leinhard, O.D. Advanced body composition assessment: From body mass index to body composition profiling. J. Investig. Med. 2018, 66, 1–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Huxley, R.; Mendis, S.; Zheleznyakov, E.; Reddy, S.; Chan, J. Body Mass Index, Waist Circumference and Waist:hip Ratio as Predictors of Cardiovascular Risk--a Review of the Literature. Eur. J. Clin. Nutr. 2010, 64, 16–22. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Cao, Q.; Yu, S.; Xiong, W.; Li, Y.; Li, H.; Li, J.; Li, F. Waist-Hip Ratio as a Predictor of Myocardial Infarction Risk: A Systematic Review and Meta-Analysis. Medicine 2018, 97, e11639. [Google Scholar] [CrossRef] [PubMed]
  17. Pulit, S.L.; Karaderi, T.; Lindgren, C.M. Sexual dimorphisms in genetic loci linked to body fat distribution. Biosci. Rep. 2017, 37, BSR20160184. [Google Scholar] [CrossRef] [PubMed]
  18. Randall, J.C.; Winkler, T.W.; Kutalik, Z.; Berndt, S.I.; Jackson, A.U.; Monda, K.L.; Kilpeläinen, T.O.; Esko, T.; Mägi, R.; Li, S.; et al. Sex-stratified Genome-wide Association Studies Including 270,000 Individuals Show Sexual Dimorphism in Genetic Loci for Anthropometric Traits. PLoS Genet. 2013, 9, e1003500. [Google Scholar] [CrossRef] [PubMed]
  19. Tang, B.; Shi, H.; Alfredsson, L.; Klareskog, L.; Padyukov, L.; Jiang, X. Obesity-Related Traits and the Development of Rheumatoid Arthritis: Evidence from Genetic Data. Arthritis Rheumatol. 2021, 73, 203–211. [Google Scholar] [CrossRef] [PubMed]
  20. Carreras-Torres, R.; Ibáñez-Sanz, G.; Obón-Santacana, M.; Duell, E.J.; Moreno, V. Identifying environmental risk factors for inflammatory bowel diseases: A Mendelian randomization study. Sci. Rep. 2020, 10, 19273. [Google Scholar] [CrossRef] [PubMed]
  21. Alwarawrah, Y.; Kiernan, K.; Maciver, N.J. Changes in Nutritional Status Impact Immune Cell Metabolism and Function. Front. Immunol. 2018, 9, 1055. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Pulit, S.L. Summary-Level Data from Meta-Analysis of Fat Distribution Phenotypes in UK Biobank and Giant. Zenodo 2018. [Google Scholar] [CrossRef]
  23. Pulit, S.L.; Stoneman, C.; Morris, A.P.; Wood, A.R.; Glastonbury, C.A.; Tyrrell, J.; Yengo, L.; Ferreira, T.; Marouli, E.; Ji, Y.; et al. Meta-analysis of genome-wide association studies for body fat distribution in 694,649 individuals of European ancestry. Hum. Mol. Genet. 2019, 28, 166–174. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Bulik-Sullivan, B.K.; Loh, P.R.; Finucane, H.K.; Ripke, S.; Yang, J.; Schizophrenia Working Group of the Psychiatric Genomics Consortium; Patterson, N.; Daly, M.J.; Price, A.L.; Neale, B.M. LD Score regression distinguishes confounding from polygenicity in genome-wide association studies. Nat. Genet. 2015, 47, 291–295. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Shi, H.; Mancuso, N.; Spendlove, S.; Pasaniuc, B. Local Genetic Correlation Gives Insights into the Shared Genetic Architecture of Complex Traits. Am. J. Hum. Genet. 2017, 101, 737–751. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Mayes, M.D.; Bossini-Castillo, L.; Gorlova, O.; Martin, J.E.; Zhou, X.; Chen, W.V.; Assassi, S.; Ying, J.; Tan, F.K.; Arnett, F.C.; et al. Immunochip Analysis Identifies Multiple Susceptibility Loci for Systemic Sclerosis. Am. J. Hum. Genet. 2014, 94, 47–61. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Burgess, S.; Thompson, S.G. Interpreting Findings from Mendelian Randomization Using the MR-Egger Method. Eur. J. Epidemiol. 2017, 32, 377–389. [Google Scholar] [CrossRef] [Green Version]
  28. Verbanck, M.; Chen, C.-Y.; Neale, B.; Do, R. Publisher Correction: Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat. Genet. 2018, 50, 1196. [Google Scholar] [CrossRef]
  29. Hemani, G.; Zheng, J.; Elsworth, B.; Wade, K.H.; Haberland, V.; Baird, D.; Laurin, C.; Burgess, S.; Bowden, J.; Langdon, R.; et al. The MR-Base platform supports systematic causal inference across the human phenome. eLife 2018, 7, e34408. [Google Scholar] [CrossRef]
  30. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B 1995, 57, 289–300. [Google Scholar] [CrossRef]
  31. Brion, M.-J.A.; Shakhbazov, K.; Visscher, P. Calculating statistical power in Mendelian randomization studies. Int. J. Epidemiology 2013, 42, 1497–1501. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Buniello, A.; MacArthur, J.A.L.; Cerezo, M.; Harris, L.W.; Hayhurst, J.; Malangone, C.; McMahon, A.; Morales, J.; Mountjoy, E.; Sollis, E.; et al. The NHGRI-EBI GWAS Catalog of published genome-wide association studies, targeted arrays and summary statistics 2019. Nucleic Acids Res. 2019, 47, D1005–D1012. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Oscanoa, J.; Sivapalan, L.; Gadaleta, E.; Ullah, A.Z.D.; Lemoine, N.R.; Chelala, C. SNPnexus: A web server for functional annotation of human genome sequence variation (2020 update). Nucleic Acids Res. 2020, 48, W185–W192. [Google Scholar] [CrossRef] [PubMed]
  34. Landrum, M.J.; Lee, J.M.; Benson, M.; Brown, G.R.; Chao, C.; Chitipiralla, S.; Gu, B.; Hart, J.; Hoffman, D.; Jang, W.; et al. ClinVar: Improving access to variant interpretations and supporting evidence. Nucleic Acids Res. 2018, 46, D1062–D1067. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Denton, C.P.; Khanna, D. Systemic Sclerosis. Lancet 2017, 390, 1685–1699. [Google Scholar] [CrossRef]
  36. Miller, J.B.; Gandhi, N.; Clarke, J.; McMahan, Z. Gastrointestinal Involvement in Systemic Sclerosis: An Update. J. Clin. Rheumatol. 2018, 24, 328–337. [Google Scholar] [CrossRef]
  37. Hughes, M.; Heal, C.; Siegert, E.; Hachulla, E.; Airó, P.; Riccardi, A.; Distler, O.; Matucci-Cerinic, M. EUSTAR Collaborators Significant Weight Loss in Systemic Sclerosis: A Study from the EULAR Scleroderma Trials and Research (EUSTAR) Database. Ann. Rheum. Dis. 2020, 79, 1123–1125. [Google Scholar] [CrossRef]
  38. Davies, N.M.; Holmes, M.V.; Smith, G.D. Reading Mendelian randomisation studies: A guide, glossary, and checklist for clinicians. BMJ 2018, 362, k601. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Julià, A.; Martínez-Mateu, S.H.; Domènech, E.; Cañete, J.D.; Ferrándiz, C.; Tornero, J.; Gisbert, J.P.; Fernández-Nebro, A.; Daudén, E.; Barreiro-de Acosta, M.; et al. Food groups associated with immune-mediated inflammatory diseases: A Mendelian randomization and disease severity study. Eur. J. Clin. Nutr. 2021, 75, 1368–1382. [Google Scholar] [CrossRef] [PubMed]
  40. González-Serna, D.; Villanueva-Martin, G.; Acosta-Herrera, M.; Márquez, A.; Martín, J. Approaching Shared Pathophysiology in Immune-Mediated Diseases through Functional Genomics. Genes 2020, 11, 1482. [Google Scholar] [CrossRef]
  41. Francisco, V.; Pino, J.; Campos-Cabaleiro, V.; Ruiz-Fernández, C.; Mera, A.; Gonzalez-Gay, M.A.; Gómez, R.; Gualillo, O. Obesity, Fat Mass and Immune System: Role for Leptin. Front. Physiol. 2018, 9, 640. [Google Scholar] [CrossRef] [Green Version]
  42. Żółkiewicz, J.; Stochmal, A.; Rudnicka, L. The Role of Adipokines in Systemic Sclerosis: A Missing Link? Arch. Dermatol. Res. 2019, 311, 251–263. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Iannone, F.; Praino, E.; Rotondo, C.; Natuzzi, D.; Bizzoca, R.; Lacarpia, N.; Fornaro, M.; Cacciapaglia, F. Body mass index and adipokines/cytokines dysregulation in systemic sclerosis. Clin. Exp. Immunol. 2021, 206, 153–160. [Google Scholar] [CrossRef] [PubMed]
  44. Driskell, R.R.; Lichtenberger, B.M.; Hoste, E.; Kretzschmar, K.; Simons, B.D.; Charalambous, M.; Ferron, S.R.; Herault, Y.; Pavlovic, G.; Ferguson-Smith, A.C.; et al. Distinct fibroblast lineages determine dermal architecture in skin development and repair. Nature 2013, 504, 277–281. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Salzer, M.C.; Lafzi, A.; Berenguer-Llergo, A.; Youssif, C.; Castellanos, A.; Solanas, G.; Peixoto, F.O.; Attolini, C.S.-O.; Prats, N.; Aguilera, M.; et al. Identity Noise and Adipogenic Traits Characterize Dermal Fibroblast Aging. Cell 2018, 175, 1575–1590.e22. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Zhang, Z.; Shao, M.; Hepler, C.; Zi, Z.; Zhao, S.; An, Y.A.; Zhu, Y.; Ghaben, A.L.; Wang, M.-Y.; Li, N.; et al. Dermal adipose tissue has high plasticity and undergoes reversible dedifferentiation in mice. J. Clin. Investig. 2019, 129, 5327–5342. [Google Scholar] [CrossRef] [PubMed]
  47. Taki, Z.; Gostjeva, E.; Thilly, W.; Yaseen, B.; Lopez, H.; Mirza, M.; Hassuji, Z.; Vigneswaran, S.; Abdi, B.A.; Hart, A.; et al. Pathogenic Activation of Mesenchymal Stem Cells Is Induced by the Disease Microenvironment in Systemic Sclerosis. Arthritis Rheumatol. 2020, 72, 1361–1374. [Google Scholar] [CrossRef] [PubMed]
Jcm 11 06014 g001 550
Figure 1. Schematic representation of the study design. The study is divided into several phases, i.e., selection of the instrumental variables for the outcome and the exposures, data harmonization and generation of different Mendelian randomization models.
Figure 1. Schematic representation of the study design. The study is divided into several phases, i.e., selection of the instrumental variables for the outcome and the exposures, data harmonization and generation of different Mendelian randomization models.
Jcm 11 06014 g001
Jcm 11 06014 g002 550
Figure 2. Pairwise global genetic correlation observed between the three obesity-related exposures and SSc. * = p > 0.05 (suggestive for statistical significance); ** = p > 0.00625 (Bonferroni-corrected).
Figure 2. Pairwise global genetic correlation observed between the three obesity-related exposures and SSc. * = p > 0.05 (suggestive for statistical significance); ** = p > 0.00625 (Bonferroni-corrected).
Jcm 11 06014 g002
Table
Table 1. Association between genetically predicted obesity-related traits and risk of SSc. Analysis including index and secondary signals for the obesity-related traits and excluding the HLA region. BMI: body mass index, WHR: waist to hip ratio, WHRadjBMI: WHR adjusted for BMI, MR: Mendelian randomization, nSNPs: number of single nucleotide polymorphisms, OR: Odds Ratio, CI: confidence interval, p: p value, p adj: p value after FDR correction for multiple testing, IVW: inverse-variance weight, PRESSO: pleiotropy residual sum and outlier, NA: not applicable.
Table 1. Association between genetically predicted obesity-related traits and risk of SSc. Analysis including index and secondary signals for the obesity-related traits and excluding the HLA region. BMI: body mass index, WHR: waist to hip ratio, WHRadjBMI: WHR adjusted for BMI, MR: Mendelian randomization, nSNPs: number of single nucleotide polymorphisms, OR: Odds Ratio, CI: confidence interval, p: p value, p adj: p value after FDR correction for multiple testing, IVW: inverse-variance weight, PRESSO: pleiotropy residual sum and outlier, NA: not applicable.
Index and secondary SNPs (p < 5 × 10−9)
MR ApproachnSNPsOR (95% CI)ppadjpfor Heterogeneity or Pleiotropy
BMIMR-Egger5331.0575(0.6403–1.7466)0.8270.82730.6005
Random-effects IVW0.9326(0.7787–1.117)0.4490.4485<0.001
MR-PRESSO (1) *0.943(0.7892–1.1269)0.5189NANA
WHRMR-Egger2470.2698(0.0914–0.7965)0.01850.03840.0519
Random-effects IVW0.7564(0.5567–1.0278)0.07430.11145<0.001
MR-PRESSO (3) *0.7809(0.5907–1.0324)0.0838NANA
WHRadjBMIMR-Egger2620.4251(0.2014–0.8971)0.02560.03840.1344
Random-effects IVW0.7269(0.5603–0.9431)0.01630.0489<0.001
MR-PRESSO (1) *0.77(0.6015–0.9857)0.039NANA
* Number of outlier SNPs detected by MR-PRESSO.
Table
Table 2. Association between genetically predicted obesity-related traits and risk of SSc. Analysis including index and secondary signals for the obesity-related traits and excluding the HLA region and known obesity-related confounder SNPs. BMI: body mass index, WHR: waist to hip ratio, WHRadjBMI: WHR adjusted for BMI, MR: Mendelian randomization, nSNPs: number of single nucleotide polymorphisms, OR: odds ratio, CI: confidence interval, p: p value, p adj: p value after FDR correction for multiple testing, IVW: inverse-variance weight, PRESSO: pleiotropy residual sum and outlier, NA: not applicable.
Table 2. Association between genetically predicted obesity-related traits and risk of SSc. Analysis including index and secondary signals for the obesity-related traits and excluding the HLA region and known obesity-related confounder SNPs. BMI: body mass index, WHR: waist to hip ratio, WHRadjBMI: WHR adjusted for BMI, MR: Mendelian randomization, nSNPs: number of single nucleotide polymorphisms, OR: odds ratio, CI: confidence interval, p: p value, p adj: p value after FDR correction for multiple testing, IVW: inverse-variance weight, PRESSO: pleiotropy residual sum and outlier, NA: not applicable.
Index and Secondary SNPs (p < 5 × 10−9)
MR ApproachnSNPsOR (95% CI)ppadjpfor Heterogeneity or Pleiotropy
BMIMR-Egger4831.422(0.721–2.803)0.31030.31030.1769
Random-effects IVW0.909(0.741–1.115)0.35980.35980.0011
MR-PRESSO (1) *0.922(0.754–1.128)0.4288NANA
WHRMR-Egger2210.301(0.086–1.060)0.06290.094350.1391
Random-effects IVW0.752(0.535–1.057)0.10070.15105< 0.001
MR-PRESSO (2) *0.764(0.559–1.044)0.0927NANA
WHRadjBMIMR-Egger2370.335(0.137–0.819)0.01720.05160.0772
Random-effects IVW0.716(0.534–0.961)0.02610.0783< 0.001
MR-PRESSO (1) *0.769(0.582–1.015)0.0651NANA
* Number of outlier SNPs detected by MR-PRESSO.
Table
Table 3. Multivariable MR (MVMR) model including BMI. WHR and risk of SSc. Analysis including index and secondary signals for the obesity-related traits and excluding the HLA region, with and without known obesity-related confounding SNPs. BMI: body mass index, WHR: waist to hip ratio, MR: Mendelian randomization, nSNPs: number of single nucleotide polymorphisms, OR: odds ratio, CI: confidence interval, p: p value.
Table 3. Multivariable MR (MVMR) model including BMI. WHR and risk of SSc. Analysis including index and secondary signals for the obesity-related traits and excluding the HLA region, with and without known obesity-related confounding SNPs. BMI: body mass index, WHR: waist to hip ratio, MR: Mendelian randomization, nSNPs: number of single nucleotide polymorphisms, OR: odds ratio, CI: confidence interval, p: p value.
Before Confounder SNP Removal.After Confounder SNP Removal
Index and secondary SNPs (p < 5 × 10−9)Index and secondary SNPs (p < 5 × 10−9)
OutcomeExposurenSNPOR (95% CI)pOutcomeExposurenSNPOR (95% CI)p
SScBMI6661.026(0.79–1.331)0.849SScBMI6101.027(0.760–1.387)0.863
WHR6660.804(0.573–1.128)0.207WHR6100.812(0.552–1.195)0.291
Index SNPs (p < 5 × 10−9)Index SNPs (p < 5 × 10−9)
OutcomeExposurenSNPOR (95% CI)pOutcomeExposurenSNPOR (95% CI)p
SScBMI5810.99(0.749–1.309)0.946SScBMI5241.013(0.726–1.412)0.941
WHR5810.876(0.607–1.263)0.477WHR5240.881(0.574–1.352)0.561
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Villanueva-Martin, G.; Acosta-Herrera, M.; Kerick, M.; López-Isac, E.; Simeón, C.P.; Callejas, J.L.; Assassi, S.; Beretta, L.; SSc Group, I.; , A.S.I.G.; et al. The Effect of Body Fat Distribution on Systemic Sclerosis. J. Clin. Med. 2022, 11, 6014. https://doi.org/10.3390/jcm11206014

AMA Style

Villanueva-Martin G, Acosta-Herrera M, Kerick M, López-Isac E, Simeón CP, Callejas JL, Assassi S, Beretta L, SSc Group I, ASIG, et al. The Effect of Body Fat Distribution on Systemic Sclerosis. Journal of Clinical Medicine. 2022; 11(20):6014. https://doi.org/10.3390/jcm11206014

Chicago/Turabian Style

Villanueva-Martin, Gonzalo, Marialbert Acosta-Herrera, Martin Kerick, Elena López-Isac, Carmen P. Simeón, José L. Callejas, Shervin Assassi, Lorenzo Beretta, International SSc Group, Australian Scleroderma Interest Group (ASIG), and et al. 2022. "The Effect of Body Fat Distribution on Systemic Sclerosis" Journal of Clinical Medicine 11, no. 20: 6014. https://doi.org/10.3390/jcm11206014

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop