Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Dietary Factors Impact on the Association between CTSS Variants and Obesity Related Traits

  • Henri Hooton ,

    Contributed equally to this work with: Henri Hooton, Lars Ängquist

    henri.hooton@upmc.fr

    Affiliation Institut national de la santé et de la recherché médicale (INSERM), U872, Nutriomique, Paris, France; Université Pierre et Marie Curie-Paris Paris, France, 6, Centre de Recherche des Cordeliers, U872, Paris, France; Université Paris Descartes, Paris, France

  • Lars Ängquist ,

    Contributed equally to this work with: Henri Hooton, Lars Ängquist

    Affiliation Institute of Preventive Medicine, Copenhagen University Hospital, Copenhagen, Denmark

  • Claus Holst,

    Affiliation Institute of Preventive Medicine, Copenhagen University Hospital, Copenhagen, Denmark

  • Jorg Hager,

    Affiliation Centre national de genotypage (CNG), Paris, France

  • Francis Rousseau,

    Affiliation INTEGRAGEN, Paris, France

  • Rikke D. Hansen,

    Affiliation Danish Cancer Society, Institute of Cancer Epidemiology, Copenhagen, Denmark

  • Anne Tjønneland,

    Affiliation Danish Cancer Society, Institute of Cancer Epidemiology, Copenhagen, Denmark

  • Nina Roswall,

    Affiliation Danish Cancer Society, Institute of Cancer Epidemiology, Copenhagen, Denmark

  • Daphne L. van der A,

    Affiliation National Institute for Public Health and the Environment (RIVM), Bilthoven, The Netherlands

  • Kim Overvad,

    Affiliations Department of Cardiology, Aalborg Hospital, Aarhus University Hospital, Aalborg, Denmark, Department of Clinical Epidemiology, Aarhus University Hospital, Aalborg, Denmark

  • Marianne Uhre Jakobsen,

    Affiliation Department of Cardiology, Aalborg Hospital, Aarhus University Hospital, Aalborg, Denmark

  • Heiner Boeing,

    Affiliation Department of Epidemiology, German Institute of Human Nutrition, Potsdam, Germany

  • Karina Meidtner,

    Affiliation Department of Epidemiology, German Institute of Human Nutrition, Potsdam, Germany

  • Domenico Palli,

    Affiliation Molecular and Nutritional Epidemiology Unit, Cancer Research and Prevention Institute (ISPO), Florence, Italy

  • Giovanna Masala,

    Affiliation Molecular and Nutritional Epidemiology Unit, Cancer Research and Prevention Institute (ISPO), Florence, Italy

  • Nabila Bouatia-Naji,

    Affiliations Medical Research Council (MRC) Epidemiology Unit, Institute of Metabolic Science, Addenbrooke’s Hospital, Cambridge, United Kingdom, Universite Paris-Descartes, Paris, France, Institut national de la santé et de la recherché médicale (INSERM) U970 Paris Cardiovascular Research Centre, Paris, France

  • Wim H. M. Saris,

    Affiliation Department of Human Biology, Nutrition and Toxicology Research Institute of Maastricht (NUTRIM), Maastricht, The Netherlands

  • Edith J. M. Feskens,

    Affiliation Division of Human Nutrition, Wageningen University, Wageningen, The Netherlands

  • Nicolas J. Wareham,

    Affiliation Medical Research Council (MRC) Epidemiology Unit, Institute of Metabolic Science, Addenbrooke’s Hospital, Cambridge, United Kingdom

  • Karani S. Vimaleswaran,

    Affiliations Medical Research Council (MRC) Epidemiology Unit, Institute of Metabolic Science, Addenbrooke’s Hospital, Cambridge, United Kingdom, Centre for Paediatric Epidemiology and Biostatistics and MRC Centre of Epidemiology for Child Health, UCL Institute of Child Health, London, United Kingdom

  • Dominique Langin,

    Affiliations Institut national de la santé et de la recherché médicale (INSERM), U1048, Obesity Research Laboratory, Team 4, I2 MC, Institute of Metabolic and Cardiovascular Diseases, Toulouse, France, University of Toulouse, U1048, Paul Sabatier University, Toulouse, France, Clinical Biochemistry Department, Toulouse University Hospitals, Toulouse, France

  • Ruth J. F. Loos,

    Affiliation Medical Research Council (MRC) Epidemiology Unit, Institute of Metabolic Science, Addenbrooke’s Hospital, Cambridge, United Kingdom

  • Thorkild I. A. Sørensen,

    Affiliations Institute of Preventive Medicine, Copenhagen University Hospital, Copenhagen, Denmark, The Novo Nordisk Foundation Center for Basic Metabolic Research, Universiy of Copenhagen, Copenhagen, Denmark

  •  [ ... ],
  • Karine Clément

    Affiliations Institut national de la santé et de la recherché médicale (INSERM), U872, Nutriomique, Paris, France; Université Pierre et Marie Curie-Paris Paris, France, 6, Centre de Recherche des Cordeliers, U872, Paris, France; Université Paris Descartes, Paris, France, Assistance Publique-Hôpitaux de Paris, Hôpital Pitié-Salpêtrière, Département de Nutrition, Paris, France; Centre de Recherche en Nutrition Humaine-Ile de France, Paris, France

  • [ view all ]
  • [ view less ]

Abstract

Background/Aims

Cathepsin S, a protein coded by the CTSS gene, is implicated in adipose tissue biology–this protein enhances adipose tissue development. Our hypothesis is that common variants in CTSS play a role in body weight regulation and in the development of obesity and that these effects are influenced by dietary factors–increased by high protein, glycemic index and energy diets.

Methods

Four tag SNPs (rs7511673, rs11576175, rs10888390 and rs1136774) were selected to capture all common variation in the CTSS region. Association between these four SNPs and several adiposity measurements (BMI, waist circumference, waist for given BMI and being a weight gainer–experiencing the greatest degree of unexplained annual weight gain during follow-up or not) given, where applicable, both as baseline values and gain during the study period (6–8 years) were tested in 11,091 European individuals (linear or logistic regression models). We also examined the interaction between the CTSS variants and dietary factors–energy density, protein content (in grams or in % of total energy intake) and glycemic index–on these four adiposity phenotypes.

Results

We found several associations between CTSS polymorphisms and anthropometric traits including baseline BMI (rs11576175 (SNP N°2), p = 0.02, β = −0.2446), and waist change over time (rs7511673 (SNP N°1), p = 0.01, β = −0.0433 and rs10888390 (SNP N°3), p = 0.04, β = −0.0342). In interaction with the percentage of proteins contained in the diet, rs11576175 (SNP N°2) was also associated with the risk of being a weight gainer (pinteraction = 0.01, OR = 1.0526)–the risk of being a weight gainer increased with the percentage of proteins contained in the diet.

Conclusion

CTSS variants seem to be nominally associated to obesity related traits and this association may be modified by dietary protein intake.

Introduction

Obesity is caused by a large number of factors that can be summarized as an interaction between an unhealthy environment and a predisposing genetic background. There is a wide spectrum of obesity-susceptibility ranging from strictly genetically determined obesity to fully environmentally determined obesity with most individuals containing a complex mix of these factors–i.e. many individual effects of genes, environmental influences and the interaction between these two. While epidemiological approaches, including twin studies, have shown that genetic factors may account for as much as 57 to 86% of body mass index (BMI) variations [1], there have been more than 450 genes referenced in the national canter for biotechnology information “NCBI gene” database (http://www.ncbi.nlm.nih.gov/sites/entrez) as being associated with obesity, but each of these genes individually has a small effect on BMI variance. This is very well illustrated by the results of genome wide association studies on large populations which have investigated the implication of several hundreds of thousands of single nucleotide polymorphisms (SNPs) in BMI variance [2][4]. New variants influencing BMI have been discovered, yet these variants, even in combination, only explain a very small part of the observed BMI variations and therefore it seems that most of the causal variants remain to be discovered. It is most likely that the largest part of the variance of BMI or other adiposity related traits attributable to genetic factors is due to a large number of variants, each of which has a very small effect [5][7]. Nevertheless, it seems that genes influencing BMI and those influencing waist circumference and adiposity may belong to a different pool of genes [8][10]. Furthermore, the genes involved in BMI and changes in BMI over time might also belong to a different set of genes [1]. However, the possibility that some genes may have a pleiotropic effect on several adiposity phenotypes should not be excluded.

Understanding the mechanisms that underlie development of adipose tissue will contribute to the identification of novel candidate genes involved in BMI and fat mass variations during life. Using a large scale transcriptomic approach in human adipose tissue, we previously identified Cathepsin S as a putative novel biomarker of adiposity [11] produced by adipose tissue. Expression of the CTSS gene, encoding for Cathepsin S, in adipose tissue correlates with BMI in obese and lean subjects. Clinical studies also revealed that Cathepsin S circulating levels were correlated with BMI and triglycerides [11], [12]. Furthermore, both CTSS adipose tissue expression and Cathepsin S systemic circulating levels were significantly modulated by weight variations either induced by dietary change or bariatric surgery in independent studies [12], [13]. In vitro studies showed that this protease also has a local role on adipose tissue. In particular, Cathepsin S contributes to the stimulation of adipocyte differentiation by degrading fibronectin, one of the main components of extra cellular matrix [14]. In vitro studies also showed that CTSS expression and Cathepsin S secretion in adipose tissue were induced by LPS, TNF-α, and IL-1β, proinflammatory factors that are secreted by cells such as macrophages or smooth muscle cells [11]. In addition Cathepsin S belongs to a family of cystein protease that includes other proteases involved in the development of obesity. In particular, CTSK−/− [15] and CTSL−/− [16] mice are protected against diet induced obesity. These animals also have improved glucose metabolism related parameters [17].

While the metabolic phenotype of CTSS−/− mice is currently under investigation, it is not known whether CTSS variants could influence obesity-related phenotypes.

We recently found an association between obesity related phenotypes and rs2424577 [18], a variant located in CST3, the gene that encodes Cystatin C, which is the main endogenous enzymatic inhibitor of Cathepsins [19][21].

A genetic study carried out by our team showed an association between several SNPs in CTSS and metabolic features in women. Rs11576175 was found to be associated with Apo A1 and HDL levels in a group of lean women from the SUVIMAX [22], [23] study; rs10888390, rs10888394 and rs1136774 were found to be associated with Apo A1 circulating levels in a group of obese women [24]. However no consistent association was found between CTSS variants and BMI in this study, irrespective of the relatively large sample size (N = 2368 unrelated lean and obese individuals).

Based on these findings, we hypothesize that genetic variation at the CTSS locus might influence obesity related phenotypes and their variation over time. We investigated four distinct phenotypes–BMI (body mass index measured in kilograms per squared meters), body fat distribution (measured by waist circumference and waist circumference for given BMI– based on sex-study stratified initial regressions of waist vs. BMI), change in weight during follow-up (either as a quantitative outcome or as a binary weight gainers indicator–experiencing the greatest degree of unexplained annual weight gain during follow-up or not), and change in body fat distribution during follow-up (measured by change in waist circumference and waist circumference for given BMI during follow-up). These four types of phenotypes were chosen since they might be influenced by different sets of genes, although they are all in some way related to adiposity. We addressed this question in a subset of the EPIC [25] cohorts, within the DiOGenes [26], [27] (Diet Obesity and genes) project where both possible associations corresponding to CTSS-SNP main effects and some SNP-dietary interactions (GI, protein intake and energy density) were investigated. These dietary factors were chosen since several studies have suggested that diets high in protein and low in GI were beneficial for obesity prevention and weight control by enhancing satiety leading to a decreased energy intake [28], [29].

Methods

Ethics Statement

EPIC study has been approved by local review board of all participating institutions, namely the Florence Local Health Authority Ethical Committee (Italy), the Ethics Committee of the Norwich District Health Authority (UK), the Medical Ethics Committee of TNO (Netherlands Organisation for Applied Scientific Research) (the Netherlands), the Ethics Committee of the Medical Association of the State of Brandenburg (Germany), and the Danish National Committee on Biomedical Research Ethics (Denmark). Written informed consent has been obtained from all participants before joining EPIC study.

Participants

Participants came from cohorts established in eight regions within five European countries (Italy, UK, the Netherlands, Germany, Denmark) participating in the European Prospective Investigation into Cancer and Nutrition (EPIC) study [25].The cohorts were those in the EPIC that had a follow-up program including reassessment of anthropometry completed. Individuals were eligible if the following inclusion criteria were met: younger than 60 years of age at baseline and younger than 65 years at follow-up, blood sample available, had baseline information on diet, weight and height and follow-up information on weight, stable smoking habits, no cancer, cardiovascular diseases (CVD), and diabetes, and an annual weight change not more than 5 kg/year. A reported or recorded weight gain exceeding 5 kg/year is very unlikely to correspond to increased adiposity but rather much more likely to either being caused by an error in the measurement, in the data or emergence of a disease that induces water retention as oedema or ascites. As such, a total of 50,293 men and women out of 146,543 initially recruited participants were eligible to participate in our study.

Cases were defined as those individuals who had experienced the greatest degree of unexplained annual weight gain during follow-up (with an average duration of 6–8 years). They were identified by using the residuals from a regression model of annual weight change on baseline values of age, weight and height, smoking status (current/former/never smokers), and follow-up time. Regression models were run separately for each sex-country strata. For each of the five countries, except Italy, we selected 600 male and 600 female cases. As the Italian cohort consisted of a general population-based sample and of a women-only sample (population-based breast cancer screening program), men were underrepresented (27%). Approximately consistent with the sex-ratio in the Italian cohort, we selected 300 male and 900 female cases. In addition to this a random subcohort (RSC) sample was selected comprising a random sample of the total eligible cohort, drawn in such a way that the total number of noncases should generally equal the number of cases (with respect to number of individuals and sex-strata distribution). Since the original case-group sized stratified random samples resulted in some overlap of cases, in practice this was performed by random oversampling of noncases, except in Denmark where overlap between cases and subcohort was negligible (n  =  79). In total, 11,921 participants were included in the present genetic association study: 6,000 cases and a subcohort of 7,061 individuals, of which 5,921 were noncases. The demographic, anthropometric and dietary characteristics of cases, noncases and random subcohort are presented in Table 1. We used both a case-noncase group and a random subcohort group to be able to test for associations with different obesity related parameters–a categorical (dichotomous) variable in the case-control group and quantitative variables in the random sub cohort.

thumbnail
Table 1. Characteristics of participants of cases, noncases and subcohort.

https://doi.org/10.1371/journal.pone.0040394.t001

Measurements of Diet, Anthropometrics and Smoking Status

Validated country-specific food frequency questionnaires (FFQs) were used to collect dietary information at baseline [25] on GI and protein intake, assessed using the methods described earlier [30], [31], and energy density.

Details of the anthropometric measurements have been described previously [30], [31]. In brief, at baseline all participants were measured for weight and height using standard study protocols [32]. At follow-up, participants in the UK and one center in the Netherlands (Doetinchem) were measured again by trained technicians, while all other participants measured their weight at home according to the guidance provided. Therefore participants from Doetinchem were analyzed separately from other Dutch participants. As such, we analyzed the data from six study centers in five countries. The cohort from Norfolk and one of the Dutch cohorts used objective measurement, but in the analyses of the various associations between baseline factors and weight gain, there were no indication that the association results were systematically different from the cohorts in which self-measured weight were used [27], [30], [31], [33].

Information on smoking status (never, former, or current smoker) was collected via self-administered questionnaires at baseline and at follow-up. Only those who had not changed their smoking habits during follow-up were included in the analyses.

Selection of Candidate Genes and tagSNPs

We originally used the International HapMap data for European ancestry (CEU) (release 20, NCBI Build 35) to select SNPs such that full coverage of the common genetic variation in the CTSS gene (+/−5 kb) was ensured. We later checked that we still covered 100% of genetic variability with the latest HapMap version (HapMap Data Rel 27 Phase II + III, Feb 09 on NCBI B36 assembly, dbSNP B126).

The Haploview software V3.3 was used to assess the linkage disequilibrium (LD) structure between SNPs [34]. Tagger software was used to select tagSNPs with the ‘pairwise tagging only’ option and an LD measure r2 threshold of 0.8. In total, 4 SNPs were selected using the pre-requisite criteria based on the minor allele frequency (MAF) and Hardy-Weinberg Equilibrium (HWE): MAF≥5%, pHWE > 0.01.

DNA Extraction and Genotyping

Genomic DNA was extracted from the buffy coats with a salting out method [35], except for participants from the UK, for whom whole-genome amplified DNA was used. Genomic and amplified DNA samples were quality-checked, quantified and normalized to approximately 100 ng/ml and 2.0 µg before genotyping. High throughput SNP genotyping was carried out using the Illumina® GoldenGate Genotyping System at IntegraGen, France.

We subjected all SNPs to country-specific HWE genotype distribution-tests. Significant deviations from equilibrium were defined as pHWE ≤ 0.001. This threshold was chosen in order to be concordant with other genetic studies carried out in the DiOGenes project. All four SNPs passed the tests for each country and were successfully genotyped for 11,091 participants. The case group included 5,584 participants and the random subcohort included 6,566 participants of whom 5,507 were noncases (Table 1).

Genetic Variability at the CTSS Loci

Four tag SNPs were selected in order to obtain a full coverage of the common variability at the CTSS locus +/−5 kb (chromosome 1, 1q21, position 148964178 to 149009929) in the HapMap CEU population. According to the latest HapMap Data Rel 27 Phase II + III, Feb 09 on NCBI B36 assembly dbSNP B126, rs7511673 (SNP N°1) captured 7 other SNPs–rs1415148, rs12089989, rs7418501, rs7521898, rs7540874, rs12086472 and rs11587444; rs11576175 (SNP N°2) captured no other SNP, rs10888390 (SNP N°3) captured 6 other SNPs–rs2275235, rs11204722, rs16827671, rs3768018, rs4537557 and rs10888391; and rs1136774 (SNP N°4) captured 2 other SNPs–rs12568757 and rs11204725. Figure 1 shows the LD pattern for the 4 selected tag SNPs in cases and RSC respectively. There seems to be no difference in the LD pattern at the CTSS locus between the cases and RSC. Two tag SNPs–rs7511373 and rs10888390 (SNP N°3)–are in strong LD in these two groups (r2 = 0.83 in both groups). Table 2 provides Hardy-Weinberg P-values, frequencies and counts for genotypes and alleles for the 4 SNPs investigated in this study both for the cases and the RSC. None of these SNPs significantly deviated from Hardy-Weinberg equilibrium in both the cases and the RSC (all pHWE>0.05).

thumbnail
Figure 1. Linkage disequilibrium (LD) plot of the CTSS locus in cases and random subcohort.

This Figure shows LD (linkage disequilibrium) values (r2) between each tag SNP in (A) cases and (B) subcohort. Each diamond contains the LD value (r2) between the two SNPs that face each of the upper sides of the diamond, ex: the LD between rs10888390 (SNP N°3) and rs1136774 (SNP N°4) is r2 = 0.62; the darker the diamond, the higher the LD value. There seems to be no difference in the LD pattern at the CTSS locus between the cases and the subcohort.

https://doi.org/10.1371/journal.pone.0040394.g001

thumbnail
Table 2. Description of CTSS variability in subcohort and cases.

https://doi.org/10.1371/journal.pone.0040394.t002

Statistical Methods

Each SNP was coded 0, 1 and 2 according to the number of minor alleles an individual carries (0 for those homozygous for the common allele, 1 for heterozygote and 2 for those homozygous for the minor allele).

First, the association between each SNP and each quantitative phenotype was tested using linear regression, assuming an additive effect of the minor allele. Second, we tested for SNP-dietary interaction associations with quantitative phenotype in the same manner. Third, case-noncase (CNC) logistic regression analyses were run, investigating possible SNP main effects on case-status (i.e. based on the risk of being a weight-gainer in the sense outlined above). These logistic regression analyses were then repeated as described above but including SNP-dietary effects.

CNC analyses of main effects were not adjusted, whereas RSC analyses were adjusted for variables that had been included in the case-status defining model (i.e. baseline values of age, height, sex, smoking status, and follow-up time) to reduce the residual variation and potential confounding. SNP-dietary variable interaction analyses were performed by including the corresponding interaction term as well as the complementary dietary main effect term in the model. Finally, change-based analyses were additionally adjusted for corresponding baseline values (additionally including baseline BMI when considering waist circumference for given BMI), and follow-up time was not used for adjustment when considering the cross-sectional (baseline) analyses.

All association analyses were first conducted for each study center separately and then effect-estimates were meta-analyzed. We used random effects to account for the possible heterogeneity across study centers, which presence was tested for using the Cochran Q-test [36].

Nominally significant associations (p<0.05) were retested assuming a dominant and a recessive model in the same way as described above.

All association analyses were conducted using Stata 9.2/11.1 for Windows (StataCorp LP, Texas, USA). The descriptive analyses were performed with SAS 9.1 for Windows (SAS Institute, Cary, NC).

Power calculations were performed using QUANTO software, Version 1.2.4 (May 2009) [37]. In the CNC analysis, the minimum detectable main effects, at 80% power, were ORs (odds ratios) 1.08 for rs7511673 (SNP N°1), 1.13 for rs11576175 (SNP N°2), and 1.08 for both rs10888390 (SNP N°3) and rs1136774 (SNP N°4). In the RSC analysis, the minimum detectable main effects, at 80% power, for weight change during the study, were regression coefficients (β) 40 g/y for rs7511673 (SNP N°1), 66 g/y for rs11576175 (SNP N°2), 41 g/y for rs10888390 (SNP N°3) and 40 g/y for rs1136774 (SNP N°4).

Results

Association between CTSS SNPs and BMI at Baseline

We found that the minor allele of rs11576175 (SNP N°2) was associated with lower BMI at baseline (p = 0.02, β = −0.24, Figure S1, Table 3). When tested assuming a dominant model, the association was also significant (p = 0.01, β = −0.29, Table S1).

thumbnail
Table 3. Associations between CTSS SNPs, BMI and body fat distribution at baseline.

https://doi.org/10.1371/journal.pone.0040394.t003

Association between CTSS SNPs and Body Fat Distribution at Baseline

No significant association between studied SNPs and body fat distribution were found (Table 3).

Association between CTSS SNPs and Annual Weight Change

The interaction between rs11576175 (SNP N°2) and the percentage of proteins contained in the diet was significantly associated to case-status (interaction p = 0.01, OR = 1.05, Table 4). For each additional minor allele, the estimated risk of being a weight gainer increases by 1.05 odds per extra one percent of proteins in the diet. This association was also significant in this population when assuming a dominant model (p = 0.004, OR = 1.06, Table S1).

thumbnail
Table 4. Association between CTSS SNPs and weight change during the study.

https://doi.org/10.1371/journal.pone.0040394.t004

Association between CTSS SNPs and Annual Body Fat Distribution Change

Both rs7511673 (SNP N°1) and rs10888390 (SNP N°3) were associated with annual waist change (p = 0.01, β = −0.04, Figure S2, Table 5 and p = 0.04, β = −0.03, Figure S3, Table 5 respectively). Rs7511673 (SNP N°1) was associated with a change in waist circumference of 0.04 cm per year and per minor allele and rs10888390 (SNP N°3) was associated with a change in waist circumference of 0.03 cm per year and per minor allele. Nevertheless these two SNPs are in strong LD in our populations (r2 = 0.83, Figure 1). The association between rs7511673 (SNP N°1) and waist gain was significant when assuming a dominant model (p =  0.02, β = −0.06, Table S1). Rs7511673 (SNP N°1) was also associated with change in waist circumference for given BMI (p = 0.03, β = −0.03, Figure S4, Table 5)–rs7511673 (SNP N°1) was associated with a change in waist circumference of 0.03 cm per year and per minor allele. This association was significant when assuming a dominant model (p = 0.02, β = −0.04, Table S1).

thumbnail
Table 5. Association between CTSS SNPs and body fat distribution change during the study.

https://doi.org/10.1371/journal.pone.0040394.t005

Discussion

In this study we found several associations between CTSS polymorphisms and anthropometric parameters including baseline BMI (rs11576175 (SNP N°2)), waist change over time (rs7511673 (SNP N°1) and rs10888390 (SNP N°3)). Although this waist change (0.03–0.04 cm/yr) is unlikely to have clinical relevance if considered on its own, this association should rather be considered in combination with other risk factors. Importantly rs11576175 (SNP N°2) was also associated with the risk of being a weight gainer, and this association was under the influence of the percentage of proteins contained in the diet. Rs7511673 (SNP N°1) captured 7 other SNPs and rs10888390 (SNP N°3) captured 6 other SNPs, besides this, two tag SNPs–rs7511673 (SNP N°1) and rs10888390 (SNP N°3)–are in LD in both of our study groups (r2 = 0.83 in each group), which means that any association with one of these variants could be caused by one of at least 14 other SNPs. There is a controversy regarding the role of fat intake on obesity related phenotypes–some studies found that fat intake had an important role [38] whereas others found that it had no importance at all [39][49]. Furthermore a study carried out in the EPIC cohorts, which investigated the role of fat intake on body weight change yielded no significant association between the type or amount of dietary fat and weight change [50]. For this reason, we decided not to investigate the interaction between CTSS SNPs and the type or amount of dietary fat in our study.

Many statistical tests have been performed therefore the question of multiple testing should be raised. The p-values presented in our study are uncorrected in order to avoid conservative corrections and loss of power (after correcting by an FDR adjustment [51] (data not shown) none of the p-values were significant). A further – although largely overlapping – motivation for not restricting the presentation and discussion to p-values adjusted for multiple comparisons is that our study is exploratory; therefore our results will need to be replicated in large independent cohorts (for related discussion, see e.g. [52], [53]).

Our group has previously published an association between CTSS variants and lipid metabolism related parameters [24]. In addition, we identified an association between a genetic variant located in CST3, a gene coding for an endogenous inhibitor of Cathepsin S, and BMI measured repeatedly during lifetime in independent European populations [18]. These observations suggest that potential alterations of Cathepsin pathway, eventually genetically induced, might contribute to changes in corpulence over time and are therefore consistent with the observations reported in this present paper. The obesity related phenotypes of CTSK−/− [15] and CTSL−/− [16] mice are also in agreement with this hypothesis [17]. Fontanesi et al [54] found an association between a CTSS polymorphism and feed:gain ratio and average daily gain in a group of Italian large white pigs. These findings seem to be in agreement with ours.

Noteworthy, CTSS has not been identified as associated to obesity related parameters by the large GWAS [2], [9]. However this may be due to the fact that these studies focus on one time point and do not investigate longitudinal data, therefore the genes that influence changes in corpulence may not be detectable by these approaches. Moreover, these studies do not account for dietary habits. Finally, it might be that these associations were not identified by GWAS simply because of the small effect size of the associations–although GWAS include many more individuals than in our study, the significance level that is generally applied in GWAS is much lower than the one applied in our study (0.05). We cannot exclude that these associations are caused by one or several variants acting on a gene nearby CTSS. CTSK, the gene that codes for Cathepsin K, an enzyme that is also involved in obesity [17], is located in the same genomic region as CTSS (1q21) [55][57]. In the HapMap CEU population, CTSSrs11576175 (SNP N°2) is in perfect LD with CTSKrs4379678 (r2 = 1), which means that the associations we found with rs11576175 (SNP N°2) might actually reflect an association with rs4379678. Furthermore we have identified a complex association between rs11576175 (SNP N°2) and the risk of being a weight gainer–the interaction between rs11576175 (SNP N°2) and the percentage of proteins in the diet was associated with the risk of being a weight gainer. A potential link between high protein diet and improved weight and fat loss has been reported [58]. These observations may be explained by the fact that proteins might be more satiating than fat or carbohydrate [59]. Very little is known concerning the molecular mechanisms underlying this process and especially regarding the potential link between Cathepsins, and in particular Cathepsin S, and dietary protein intake. The possibility that dietary changes could influence the expression of Cathepsins has been highlighted by the outcomes of both animal models and in vitro studies. In mice, after infection by Paracoccidioides brasiliensis (a fungus that causes Paracoccidioidomycosis, a systemic mycosis), a very high protein diet was associated with a greater increase in spleen and liver Cathepsin G mRNA than a low protein diet [60]. Furthermore, in vitro, pyridoxal phosphate, a coenzyme form of vitamin B6, strongly inhibits Cathepsin B activity and weakly inhibits Cathepsin S and K activities [61].

In conclusion, we have identified nominally significant associations between several CTSS variants and obesity related parameters. One of these associations seems to be influenced by dietary protein intake. However this link needs to be further investigated in order to gain knowledge on the mechanisms governing weight homeostasis.

Supporting Information

Figure S1.

BMI at baseline according to rs11576175 (SNP N°2). Mean +/− SEM of BMI at baseline according to rs11576175 genotypes (G/G n = 5341, G/A n = 1155, and A/A n = 70) in the subcohort, n = 6566. Rs11576175 was associated with a decrease of 0.24 kg/m2 per A allele (p = 0.02, β = −0.24).

https://doi.org/10.1371/journal.pone.0040394.s001

(TIF)

Figure S2.

Annual waist gain according to rs7511673 (SNP N°1). Mean +/- SEM of annual waist gain according to rs7511673 genotypes (A/A n = 2382, A/T, n = 3142, and T/T, n = 1041) in the subcohort, n = 6566. In the regression analysis rs7511673 was associated with a decrease in waist circumference of 0.04 cm per year and per T allele (p = 0.01, β = −0.04). This association was also significant when assuming a dominant model (p = 0.02, β = −0.06), A/T and T/T carriers gained 0.06 cm per year less than A/A carriers.

https://doi.org/10.1371/journal.pone.0040394.s002

(TIF)

Figure S3.

Annual waist gain according to rs10888390 (SNP N°3). Mean +/− SEM of annual waist gain according to rs10888390 genotypes (G/G n = 2721, G/A n = 2999, and A/A n = 844) in the subcohort, n = 6566. rs10888390 was associated with a decrease in waist circumference of 0.03 cm per year and per A allele (p = 0.04, β = −0.03).

https://doi.org/10.1371/journal.pone.0040394.s003

(TIF)

Figure S4.

Annual waist for given BMI gain per year according to rs7511673 (SNP N°1). Mean +/− SEM of annual waist gain for given BMI according to rs7511673 genotypes (A/A n = 2382, A/T, n = 3142, and T/T, n = 1041) in the subcohort, n = 6566. Rs7511673 was associated with a decrease in waist circumference for given BMI of 0.03 cm per year and per T allele (p = 0.03, β = −0.03). This association was also significant when assuming a dominant model (p = 0.02, β = −0.04).

https://doi.org/10.1371/journal.pone.0040394.s004

(TIF)

Table S1.

Dominant and recessive models for associations which were significant when assuming an additive model.

https://doi.org/10.1371/journal.pone.0040394.s005

(DOC)

Author Contributions

Conceived and designed the experiments: KC TIAS NJW RJFL DL. Performed the experiments: HH. Analyzed the data: LÄ CH. Contributed reagents/materials/analysis tools: KC HH DLvdA NJW AT KO HB DP TIAS WHMS NB-N NR EJMF RJFL DL RDH MUJ KM GM KSV JH FR. Wrote the paper: HH. Read and corrected the manuscript: KC TIAS LÄ KSV JH WHMS CH RJFL EJMF.

References

  1. 1. Silventoinen K, Kaprio J (2009) Genetics of tracking of body mass index from birth to late middle age: evidence from twin and family studies. Obes Facts 2: 196–202.
  2. 2. Speliotes EK, Willer CJ, Berndt SI, Monda KL, Thorleifsson G, et al. (2010) Association analyses of 249,796 individuals reveal 18 new loci associated with body mass index. Nat Genet. Accessed 15 October 2010.
  3. 3. Willer CJ, Speliotes EK, Loos RJF, Li S, Lindgren CM, et al. (2009) Six new loci associated with body mass index highlight a neuronal influence on body weight regulation. Nat Genet 41: 25–34.
  4. 4. Thorleifsson G, Walters GB, Gudbjartsson DF, Steinthorsdottir V, Sulem P, et al. (2009) Genome-wide association yields new sequence variants at seven loci that associate with measures of obesity. Nat Genet 41: 18–24.
  5. 5. Liu Y-J, Guo Y-F, Zhang L-S, Pei Y-F, Yu N, et al. (2010) Biological Pathway-Based Genome-Wide Association Analysis Identified the Vasoactive Intestinal Peptide (VIP) Pathway Important for Obesity. Obesity (Silver Spring). Accessed 26 July 2010.
  6. 6. Moonesinghe R, Liu T, Khoury MJ (2010) Evaluation of the discriminative accuracy of genomic profiling in the prediction of common complex diseases. Eur J Hum Genet 18: 485–489.
  7. 7. Naukkarinen J, Surakka I, Pietiläinen KH, Rissanen A, Salomaa V, et al. (2010) Use of Genome-Wide Expression Data to Mine the “Gray Zone” of GWA Studies Leads to Novel Candidate Obesity Genes. PLoS Genet 6: e1000976.
  8. 8. Hasselbalch AL, Benyamin B, Visscher PM, Heitmann BL, Kyvik KO, et al. (2008) Common genetic components of obesity traits and serum leptin. Obesity (Silver Spring) 16: 2723–2729.
  9. 9. Heid IM, Jackson AU, Randall JC, Winkler TW, Qi L, et al. (2010) Meta-analysis identifies 13 new loci associated with waist-hip ratio and reveals sexual dimorphism in the genetic basis of fat distribution. Nat Genet. Accessed 20 October 2010.
  10. 10. Lindgren CM, Heid IM, Randall JC, Lamina C, Steinthorsdottir V, et al. (2009) Genome-wide association scan meta-analysis identifies three Loci influencing adiposity and fat distribution. PLoS Genet 5: e1000508.
  11. 11. Taleb S, Lacasa D, Bastard J-P, Poitou C, Cancello R, et al. (2005) Cathepsin S, a novel biomarker of adiposity: relevance to atherogenesis. FASEB J 19: 1540–1542.
  12. 12. Naour N, Rouault C, Fellahi S, Lavoie M-E, Poitou C, et al. (2010) Cathepsins in human obesity: changes in energy balance predominantly affect cathepsin s in adipose tissue and in circulation. J Clin Endocrinol Metab 95: 1861–1868.
  13. 13. Taleb S, Cancello R, Poitou C, Rouault C, Sellam P, et al. (2006) Weight loss reduces adipose tissue cathepsin S and its circulating levels in morbidly obese women. J Clin Endocrinol Metab 91: 1042–1047.
  14. 14. Taleb S, Cancello R, Clément K, Lacasa D (2006) Cathepsin s promotes human preadipocyte differentiation: possible involvement of fibronectin degradation. Endocrinology 147: 4950–4959.
  15. 15. Yang M, Sun J, Zhang T, Liu J, Zhang J, et al. (2008) Deficiency and inhibition of cathepsin K reduce body weight gain and increase glucose metabolism in mice. Arterioscler Thromb Vasc Biol 28: 2202–2208.
  16. 16. Yang X, Peters C, Shi G-P, Kitamoto S, Sukhova GK, et al. (2007) Cathepsin L deficiency reduces diet-induced atherosclerosis in low-density lipoprotein receptor-knockout mice. Circulation 115: 2065–2075.
  17. 17. Lafarge J-C, Naour N, Clément K, Guerre-Millo M (2010) Cathepsins and cystatin C in atherosclerosis and obesity. Biochimie. Accessed 8 September 2010.
  18. 18. Hooton H, Dubern B, Henegar C, Paternoster L, Nohr EA, et al. (2011) Association between CST3 rs2424577 polymorphism and corpulence related phenotypes during lifetime in populations of European ancestry. Obes Facts 4: 131–144.
  19. 19. Clarke GM, Cardon LR (2010) Aspects of observing and claiming allele flips in association studies. Genet Epidemiol 34: 266–274.
  20. 20. Lin P-I, Vance JM, Pericak-Vance MA, Martin ER (2007) No gene is an island: the flip-flop phenomenon. Am J Hum Genet 80: 531–538.
  21. 21. Zaykin DV, Shibata K (2008) Genetic flip-flop without an accompanying change in linkage disequilibrium. Am J Hum Genet 82: 794–796; author reply 796–797. https://doi.org/10.1016/j.ajhg.2008.02.001
  22. 22. Hercberg S, Preziosi P, Briançon S, Galan P, Triol I, et al. (1998) A primary prevention trial using nutritional doses of antioxidant vitamins and minerals in cardiovascular diseases and cancers in a general population: the SU.VI.MAX study–design, methods, and participant characteristics. SUpplementation en VItamines et Minéraux AntioXydants. Control Clin Trials 19: 336–351.
  23. 23. Dolley G, Bertrais S, Frochot V, Bebel J-F, Guerre-Millo M, et al. (2008) Promoter adiponectin polymorphisms and waist/hip ratio variation in a prospective French adults study. Int J Obes (Lond) 32: 669–675.
  24. 24. Spielmann N, Mutch DM, Rousseau F, Tores F, Hager J, et al. (2008) Cathepsin S genotypes are associated with Apo-A1 and HDL-cholesterol in lean and obese French populations. Clin Genet 74: 155–163.
  25. 25. Riboli E, Hunt KJ, Slimani N, Ferrari P, Norat T, et al. (2002) European Prospective Investigation into Cancer and Nutrition (EPIC): study populations and data collection. Public Health Nutr 5: 1113–1124.
  26. 26. Saris WHM, Harper A (2005) DiOGenes: a multidisciplinary offensive focused on the obesity epidemic. Obes Rev 6: 175–176.
  27. 27. Du H, Vimaleswaran KS, Angquist L, Hansen RD, van der A DL, et al. (2011) Genetic Polymorphisms in the Hypothalamic Pathway in Relation to Subsequent Weight Change - The DiOGenes Study. PLoS ONE 6: e17436.
  28. 28. Du H, Van der A DL, Feskens EJM (2006) Dietary glycaemic index: a review of the physiological mechanisms and observed health impacts. Acta Cardiol 61: 383–397.
  29. 29. St Jeor ST, Howard BV, Prewitt TE, Bovee V, Bazzarre T, et al. (2001) Dietary protein and weight reduction: a statement for healthcare professionals from the Nutrition Committee of the Council on Nutrition, Physical Activity, and Metabolism of the American Heart Association. Circulation 104: 1869–1874.
  30. 30. Du H, van der A DL, van Bakel MME, Slimani N, Forouhi NG, et al. (2009) Dietary glycaemic index, glycaemic load and subsequent changes of weight and waist circumference in European men and women. Int J Obes (Lond) 33: 1280–1288.
  31. 31. Halkjær J, Olsen A, Overvad K, Jakobsen MU, Boeing H, et al. (2010) Intake of total, animal and plant protein and subsequent changes in weight or waist circumference in European men and women: the Diogenes project. Int J Obes (Lond). Accessed 10 February 2011.
  32. 32. Rinaldi S, Key TJ, Peeters PHM, Lahmann PH, Lukanova A, et al. (2006) Anthropometric measures, endogenous sex steroids and breast cancer risk in postmenopausal women: a study within the EPIC cohort. Int J Cancer 118: 2832–2839.
  33. 33. Vimaleswaran KS, Angquist L, Hansen RD, van der A DL, Bouatia-Naji N, et al. (2012) Association Between FTO Variant and Change in Body Weight and Its Interaction With Dietary Factors: The DiOGenes Study. Obesity (Silver Spring, Md). Accessed 3 May 2012.
  34. 34. Barrett JC, Fry B, Maller J, Daly MJ (2005) Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21: 263 –265. https://doi.org/10.1093/bioinformatics/bth457
  35. 35. Miller SA, Dykes DD, Polesky HF (1988) A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res 16: 1215.
  36. 36. Higgins JPT, Thompson SG, Deeks JJ, Altman DG (2003) Measuring inconsistency in meta-analyses. BMJ 327: 557–560.
  37. 37. Gauderman WJ, Morrison JM (2006) QUANTO 1.1: A computer program for power and sample size calculations for gene-epidemiology studies. Available. http://hydra.usc.edu/gxe.
  38. 38. Bray GA, Popkin BM (1998) Dietary fat intake does affect obesity! Am J Clin Nutr 68: 1157–1173.
  39. 39. Willett WC (2002) Dietary fat plays a major role in obesity: no. Obes Rev 3: 59–68.
  40. 40. Willett WC, Leibel RL (2002) Dietary fat is not a major determinant of body fat. Am J Med 113: 47S–59S.
  41. 41. Pirozzo S, Summerbell C, Cameron C, Glasziou P (2003) Should we recommend low-fat diets for obesity? Obes Rev 4: 83–90.
  42. 42. Pirozzo S, Summerbell C, Cameron C, Glasziou P (2002) Advice on low-fat diets for obesity. Cochrane Database Syst Rev. CD003640 p. https://doi.org/10.1002/14651858.CD003640
  43. 43. Carmichael HE, Swinburn BA, Wilson MR (1998) Lower fat intake as a predictor of initial and sustained weight loss in obese subjects consuming an otherwise ad libitum diet. J Am Diet Assoc 98: 35–39.
  44. 44. Sheppard L, Kristal AR, Kushi LH (1991) Weight loss in women participating in a randomized trial of low-fat diets. Am J Clin Nutr 54: 821–828.
  45. 45. Westerterp-Plantenga MS, Wijckmans-Duijsens NE, Verboeket-van de Venne WP, de Graaf K, van het Hof KH, et al. (1998) Energy intake and body weight effects of six months reduced or full fat diets, as a function of dietary restraint. Int J Obes Relat Metab Disord 22: 14–22.
  46. 46. Astrup A, Grunwald GK, Melanson EL, Saris WH, Hill JO (2000) The role of low-fat diets in body weight control: a meta-analysis of ad libitum dietary intervention studies. Int J Obes Relat Metab Disord 24: 1545–1552.
  47. 47. Sherwood NE, Jeffery RW, French SA, Hannan PJ, Murray DM (2000) Predictors of weight gain in the Pound of Prevention study. Int J Obes Relat Metab Disord 24: 395–403.
  48. 48. Rock CL, Thomson C, Caan BJ, Flatt SW, Newman V, et al. (2001) Reduction in fat intake is not associated with weight loss in most women after breast cancer diagnosis: evidence from a randomized controlled trial. Cancer 91: 25–34.
  49. 49. Brehm BJ, Spang SE, Lattin BL, Seeley RJ, Daniels SR, et al. (2005) The role of energy expenditure in the differential weight loss in obese women on low-fat and low-carbohydrate diets. J Clin Endocrinol Metab 90: 1475–1482.
  50. 50. Forouhi NG, Sharp SJ, Du H, van der A DL, Halkjaer J, et al. (2009) Dietary fat intake and subsequent weight change in adults: results from the European Prospective Investigation into Cancer and Nutrition cohorts. Am J Clin Nutr 90: 1632–1641.
  51. 51. Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society Series B, 57, 289–300.
  52. 52. Rothman KJ (1990) No adjustments are needed for multiple comparisons. Epidemiology 1: 43–46.
  53. 53. Bender R, Lange S (2001) Adjusting for multiple testing–when and how? J Clin Epidemiol 54: 343–349.
  54. 54. Fontanesi L, Speroni C, Buttazzoni L, Scotti E, Costa LN, et al. (2010) Association between cathepsin L (CTSL) and cathepsin S (CTSS) polymorphisms and meat production and carcass traits in Italian Large White pigs. Meat Sci 85: 331–338.
  55. 55. Gelb BD, Shi GP, Heller M, Weremowicz S, Morton C, et al. (1997) Structure and chromosomal assignment of the human cathepsin K gene. Genomics 41: 258–262.
  56. 56. Rood JA, Van Horn S, Drake FH, Gowen M, Debouck C (1997) Genomic organization and chromosome localization of the human cathepsin K gene (CTSK). Genomics 41: 169–176.
  57. 57. Shi GP, Webb AC, Foster KE, Knoll JH, Lemere CA, et al. (1994) Human cathepsin S: chromosomal localization, gene structure, and tissue distribution. J Biol Chem 269: 11530–11536.
  58. 58. Halton TL, Hu FB (2004) The effects of high protein diets on thermogenesis, satiety and weight loss: a critical review. J Am Coll Nutr 23: 373–385.
  59. 59. Westerterp-Plantenga MS (2003) The significance of protein in food intake and body weight regulation. Curr Opin Clin Nutr Metab Care 6: 635–638.
  60. 60. Oarada M, Igarashi M, Tsuzuki T, Kamei K, Hirasaka K, et al. (2010) Effects of a high-protein diet on host resistance to Paracoccidioides brasiliensis in mice. Biosci Biotechnol Biochem 74: 620–626.
  61. 61. Katunuma N, Matsui A, Endo K, Hanba J, Sato A, et al. (2000) Inhibition of intracellular cathepsin activities and suppression of immune responses mediated by helper T lymphocyte type-2 by peroral or intraperitoneal administration of vitamin B6. Biochem Biophys Res Commun 272: 151–155.