Open Access

Promoter polymorphisms in the chitinase 3-like 1 gene influence the serum concentration of YKL-40 in Danish patients with rheumatoid arthritis and in healthy subjects

  • Kaspar R Nielsen1Email author,
  • Rudi Steffensen1,
  • Martin Boegsted2,
  • John Baech1,
  • Soeren Lundbye-Christensen3,
  • Merete L Hetland4,
  • Sophine B Krintel4, 5,
  • Hans E Johnsen2,
  • Mette Nyegaard2 and
  • Julia S Johansen5, 6
Arthritis Research & Therapy201113:R109

DOI: 10.1186/ar3391

Received: 30 January 2011

Accepted: 29 June 2011

Published: 29 June 2011

Abstract

Introduction

The present study investigates the association between single nucleotide polymorphisms (SNPs) in the chitinase 3-like 1 (CHI3L1) gene and serum concentrations of YKL-40 in Danish patients with rheumatoid arthritis (RA) and healthy controls as well as the association with RA in the Danish population. The CHI3L1 gene is located on chromosome 1q32.1 and encodes the YKL-40 glycoprotein. YKL-40 concentrations are elevated in the serum of patients with RA compared to healthy subjects, and YKL-40 has been suggested to be an auto-antigen and may play a role in development of RA and in inflammation.

Methods

Eight SNPs in the CHI3L1 gene and promotor were genotyped in 308 patients with RA and 605 controls (healthy blood donors) using TaqMan allele discrimination assays. Serum concentrations of YKL-40 were determined by an enzyme-linked immunosorbent assay (ELISA).

Results

We found significant association between the serum concentrations of YKL-40 and polymorphism in the CHI3L1 gene among both patients with RA and controls. The g.-131(C > G) polymorphism (rs4950928) was most strongly associated with age adjusted serum concentrations of YKL-40 in patients with RA (P < 2.4e-8) and controls (P < 2.2e-16). No significant allelic- or genotypic association with RA was found in this Danish cohort.

Conclusions

We suggest that the g.-131(C > G) promoter polymorphism has a substantial impact on serum concentrations of YKL-40 in patients with RA and healthy subjects. However, the polymorphism does not seem to confer risk to RA itself. The effect of CHI3L1 polymorphism on clinical outcome or the response to treatment in patients with RA remains to be investigated.

Introduction

Rheumatoid arthritis (RA) is a systemic autoimmune inflammatory disorder, affecting approximately 1% in western populations. The disease is primarily characterised by chronic polyarthritis [1, 2]. The aetiology of RA remains unknown, although it is estimated that the contribution of genetic factors is about 50 to 60% [3, 4]. The strongest genetic association is with polymorphic alleles within the human leukocyte antigen HLA-DRB1 locus on chromosome 6p21.3 and a single nucleotide polymorphism (SNP) in the PTPN22 gene on chromosome 1p13.2 [5]. Another proposed potential loci is on chromosome 1q32.1 harbouring the chitinase 3-like 1 (CHI3L1) gene encoding the YKL-40 protein [6]. YKL-40 is a 40 kDa heparin- and chitin-binding glycoprotein, and a member of chitinase like proteins. YKL-40 is expressed by a variety of cells, including macrophages, neutrophils, synovial cells, arthritic chondrocytes and cancer cells [710]. As YKL-40 contains HLA-DR4 binding motifs, it has been suggested to function as an auto antigen in RA [1115].

A high serum concentration of YKL-40 is emerging as a new biomarker of severe disease activity and poor prognosis in patients with diseases characterized by inflammation and ongoing tissue remodelling such as RA, inflammatory bowel disease, asthma and cancer [8, 10, 1626]. The exact biological function of the YKL-40 protein is still largely elusive. YKL-40 is a trans-membrane protein in which cleavaged components bind to an unidentified receptor and the expression of YKL-40 is regulated by various inflammatory cytokines and hormones [2730]. It is suggested that YKL-40 plays a role in cell proliferation, differentiation and protection against apoptotic signals, and has an effect on extracellular tissue remodelling [31, 32]. Two recent studies have explored the effect of YKL-40 as a stimulator of angiogenesis in tumours, suggesting that anti-YKL-40 antibodies could have a place in cancer treatment [33, 34].

The proximal promoter region of the CHI3L1 gene contains a highly polymorphic area, suggesting a possibility for several functional variants of the gene. Rehli et al. [35] demonstrated that binding of the SP1 transcription factor to the most proximal part of the CHI3L1 gene affected gene transcription. This finding was supported by Zhao et al. [36] reporting functional variants based on the binding of the MYC/MAX transcription factors to the proximal promoter region. The relationships between CHI3L1 polymorphisms and YKL-40 production have been studied in a small number of patients with various inflammatory disorders, such as sarcoidosis, asthma, hepatitis, schizophrenia and diabetes [3744]. These studies suggest that serum concentrations of YKL-40 are, at least partly, regulated by polymorphisms in the proximal promotor region. The findings have been somewhat contradictory and the exact position of the regulatory site or sites remains to be demonstrated. Allele frequencies differ significantly between Caucasian, African and Asian populations, and possibly even within these populations, thereby making direct comparison of the reported studies difficult [45].

Only one small study has evaluated CHI3L1 polymorphisms in patients with RA [46]. In 182 Hungarian patients with RA and 194 healthy controls there were no significant differences in genotype frequencies for the g.-131(C > G) or the g.-329(C > T) polymorphisms between the two groups. This study did not evaluate the functional properties of these polymorphisms. Several questions remain unanswered, namely the relationship between CHI3L1 polymorphisms and serum concentrations of YKL-40 in patients with RA, the association of CHI3L1 promoter genotypes to risk of RA and the Linkage Disequilibrium (LD) properties in different populations.

We aimed to investigate these questions in a cohort of well defined Danish patients with RA and a group of healthy Danish controls. Our hypothesis was that polymorphisms in the proximal promoter region of CHI3L1, most likely the g.-131(C > G) polymorphism (rs4950928), are associated with serum concentrations of YKL-40 in both patients with RA and healthy controls. Moreover, we hypothesized that these polymorphisms could be associated with the risk of developing RA and possibly also associated to IgM rheumatoid factor (RF), since YKL-40 seems to play a role in the pathogenesis and immunomodulation in RA.

Materials and methods

Patients with rheumatoid arthritis

Three-hundred and eight patients with RA treated at the Department of Rheumatology, Hvidovre Hospital, Hvidovre, Denmark were included in the study. The patients had RA according to the ACR 1987 criteria [47]. The patients with available blood samples were identified in the DANBIO Registry (The Copenhagen Cohort). DANBIO is a Danish nationwide registry that prospectively collects clinical data on patients with rheumatic diseases receiving medical treatment [48]. The blood samples (serum and whole blood) were collected at the time of diagnosis or at the time of starting treatment with TNFα inhibitors. All patients provided informed consents for inclusion in the study population. The study was approved by the local ethics committee. Table 1 summarizes the demographic data for the patients with RA and the controls.
Table 1

Characteristics of the study population

 

Group

 

All RA

(n = 308)

IgM RFpos RA

(n = 178)

IgM RFneg RA

(n = 130)

Controls

(n = 605)

Age in years (mean ± SD and range)

54.5 ± 14.7

(22 to 93)

56.2 ± 14.0

(22 to 86)

52.4 ± 15,4

(23 to 93)

42.6 ± 12.8

(19 to 65)

Male/Female

74/234

47/131

27/103

367/238

Serum YKL-40 ng/ml (median and 95% CI)

86 (79 to 94)

91 (81 to 102)

80 (70 to 91)

46 (44 to 48)

CI, confidence interval; RA, rheumatoid arthritis; RF, rheumatoid factor.

Healthy controls

Six-hundred and five healthy blood donors from the Aalborg Hospital Blood Bank, Aalborg, Denmark were included in the study. The donors were known not to take any medication and were clinically healthy at the time of blood drawing. The over-2representation of female controls was a random phenomenon. The samples were handled anonymously and all donors gave consent to the blood being used for this purpose and the sampling was approved by the local ethics committee.

Handling of blood samples

From the patients with RA and blood donors Ethylenediaminetetraacetic acid (EDTA)-stabilised whole blood and blood samples without anticoagulants were drawn. Serum was isolated from coagulated whole blood within three hours and stored at -80°C until analysis of YKL-40 and IgM-RF was performed. Genomic DNA was prepared from EDTA-stabilised blood samples using a Maxwell 16 blood DNA purification kit (Promega, Madison, WI, USA).

Biochemical analysis

Serum concentration of YKL-40 was measured by a commercial two-site sandwich type ELISA (Quidel, Mountain View, CA, USA) [49]. The detection limit was 10 ng/ml. The intra-assay coefficient of variations (CV) was 5% and the inter-assay CV was < 6%. IgM-RF was measured using an ELIA fluorescence immunoassay on a Unicap250 system (Phadia AB, Uppsala, Sweden). A validated diagnostic cut off (< 17 kI U/l) was used to classify patients as IgM-RF negative or IgM-RF positive.

Genotyping

A total of eight SNPs located within the promoter or coding regions of the CHI3L1 gene was analysed. Genotyping was performed using real-time polymerase chain reaction (rt-PCR) with TaqMan® SNP Genotyping Assays (Applied Biosystems, Foster City, CA, USA). Applied Biosystems: assay identification numbers are reported as SNP identification. DNA amplification was carried out in a 5 μl volume containing 20 ng DNA, 0.9 μM primers and 0.2 μM probes (final concentrations). The product was amplified using TaqMan Universal PCR Master Mix (Applied Biosystems). Reactions were performed in 384-well plates with the following protocol on a GeneAmp PCR 9700 or a 7900 HT Sequence Detection System: 95°C for 10 minutes, followed by 40 cycles at 95°C for 15 seconds and 60°C for 1 minute. To determine genotypes, end-point fluorescence was read on the 7900 HT Sequence Detection Systems using SDS version 2.3 software (QIAGEN Inc. 27220 Turnberry Lane, CA 91355, USA).

Statistical analysis

The genotype distribution among patients with RA and controls was tested for deviation from Hardy-Weinberg equilibrium and haplotypes were estimated using the Helix Tree SNP analysis software package (Golden Helix Software, Bozeman, MT, USA). The degree of LD between the SNPs was determined using the SHEsis software (Bio-X Center, Shanghai Jiao Tong University, 1954 Huashan Road, Shanghai 200030, China) [50]. Serum concentrations of YKL-40 were log-normally distributed and, therefore, log-transformed before analysis. Statistical analysis was performed using the statistical software system R, version 2.12.1 [51]. The initial non-linear association between serum concentrations of YKL-40 and age was modelled by a restricted cubic spline function, using the user-contributed package design [52] integrated in R. Analysis of variance based on multiple linear regression models was used to investigate the association between age, gender, case-control status, genotypes and serum YKL-40. Prior to SNP-wise association analysis with serum YKL-40, all serum concentrations of YKL-40 were age adjusted to 44.4 years (mean age for the total sample of controls and cases age 65 years and below) using a linear model. Genotypic associations with age-adjusted serum concentrations of YKL-40 were carried out for cases (age 65 years and below) and controls separately using a multiple linear regression model. For association analysis with RA, allelic and genotypic association was performed using Fisher:s exact test including all patients (n = 308) and controls (n = 605) and using a significance level of 0.05.

Results

No deviations from Hardy-Weinberg equilibrium were found for any of the eight SNPs in the patient or control group. Age stratification into one-year age groups did not reveal deviations from Hardy-Weinberg equilibrium in any of the age groups.

Prior to the SNP association analysis, the effect of age and case-control status on serum YKL-40 was tested using a multiple linear regression model, with serum YKL-40 as dependent variable and case-control status and a non-linear function of age included as covariate. Strong significant association of the serum concentration of YKL-40 with age (P < 2.0e-16) and case-control status (P < 2.0e-16) was observed (Figure 1). Moreover an apparent increase in serum YKL-40 with age was found for the older patients in the case group. To avoid a potential bias due to the high influence of individuals older than 65 years in the RA group, we excluded in all further analysis patients with RA older than 65 years.
https://static-content.springer.com/image/art%3A10.1186%2Far3391/MediaObjects/13075_2011_Article_3136_Fig1_HTML.jpg
Figure 1

Non-linear association between age and serum concentrations of YKL-40. Restricted cubic spline model with six knots applied for patients with rheumatoid arthritis (RA) and controls.

To test the effect of genotypes on serum concentrations of YKL-40 in the RA group (age 65 and below) and control group, a multiple linear regression model including serum YKL-40 as dependent variable and a non-linear function of age, case-control status, genotypes and gender as well as the interaction between case-control status and genotype with age as independent variables was applied (Table 2). From this analysis a strong association was observed with case-control status (P < 2.0e-16) (as before), age (P < 2.0e-16) (as before) and genotype (P < 2.0e-16).
Table 2

Sequential analysis of variance table for the regression model for patients ≤ 65 years of age.

Effects

Df

SSQ

F value

Pr(> F)

Main Effects

    

   Age

1

41.844

180.7462

< 2e-16

   Status

1

23.888

103.1865

< 2e-16

   Genotype

14

51.14

15.7787

< 2e-16

   Sex

1

0.516

2.2304

0.1357

Interactions effects

    

   Genotype * Status

13

4.202

1.3961

0.1551

   (Status + Genotype) * Age

13

1.671

0.5552

0.8896

   (Status + Genotype) * RCS(Age)

55

14.932

1.1727

0.1892

Error

    

   Residual

744

172.24

  

RCS, (restricted cubic spline denotes a non-linear relation with age)

Regarding the age-dependent increase in the serum concentrations of YKL-40, no significant difference was found between a non-linear and a linear model for the age-dependence in both the case group (age 65 and below) and control group (P = 0.19) suggesting that the linear model can be used for age adjustment of the serum concentrations of YKL-40 in both groups. The linear model was fitted and depicted in Figure 2.
https://static-content.springer.com/image/art%3A10.1186%2Far3391/MediaObjects/13075_2011_Article_3136_Fig2_HTML.jpg
Figure 2

Linear association between age and serum concentrations of YKL-40. Linear model applied for patients with rheumatoid arthritis (RA) ≤ 65 years of age (n = 238) and controls (n = 605) is sufficient to explain the age dependent variation (P = 0.96). The y-axis represents serum concentrations of YKL-40. Dotted lines represent 95% confidence intervals.

Serum concentrations of YKL-40 were not associated with gender (P = 0.16). There were no interaction effects between case-control status or genotype and age (P = 0.89) and no association between serum YKL-40 and the interaction effect between genotype and case-control status (P = 0.16) (Table 2). This suggests that age, case-control status and genotypes are all strong independent factors affecting serum concentrations of YKL-40.

To test the association of each SNP on age-adjusted serum YKL-40 in the RA group (age 65 and below) and control group, a linear age-adjustment was applied and genotypes were included one-by-one as dependent variables in a multiple linear regression analysis. The g.-131(C > G) genotype was found to be most strongly associated with age-adjusted serum concentrations of YKL-40 in both the patient (P = 2.4e-08) and control group (P < 2.2e-16) (Table 3). Consistently within both groups, the rare GG genotype was associated with low serum YKL-40, the CG genotype with intermediate serum concentrations of YKL-40, and the common CC genotype with high serum YKL-40 (Figure 3). With respect to genotypes, the RA patients had significantly higher serum YKL-40 than controls for both the CC and CG group. For the rare GG group, the difference was not significant, most likely because of low statistical power due to the limited number of individuals in the GG groups.
Table 3

HI3L1 genotypes and the effect on serum YKL-40 levels in patients with rheumatoid arthritis and controls.

SNP

CHI3L1position

Controls (n = 605)

F df = 603

P-value

RA age ≤ 65 years (n = 238)

F

Df = 236

P-value

rs6691378

g.-1371G/A

Serum YKL-40 ng/ml

  

serum YKL-40 ng/ml

  
 

G/G

46 (44 to 48)

6.63

0.0014

67 (61 to 74)

2.85

0.06

 

G/A

54 (50 to 59)

  

77 (62 to 95)

  
 

A/A

47 (34 to 68)

  

151 (71 to 323)

  

rs10399931

g.-329C/T

      
 

C/C

56 (53 to 58)

83.13

< 2.2e-16

80 (71 to 89)

12.34

8.0e-06

 

C/T

40 (38 to 42)

  

61 (53 to 70)

  
 

T/T

25 (22 to 29)

  

35 (24 to 49)

  

rs10399805

g.-247G/A

      
 

G/G

45 (44 to 47)

7.47

6.3e-3

67 (61 to 74)

2.69

0.07

 

G/A

54 (50 to 59)

  

76 (61 to 94)

  
 

A/A

48 (34 to 68)

  

151 (71 to 323)

  

rs4950928

g.-131C/G

      
 

C/C

56 (53 to 58)

102.32

< 2.2e-16

81 (73 to 90)

18.91

2.4e-08

 

C/G

38 (38 to 40)

  

59 (51 to 68)

  
 

G/G

25 (22 to 28)

  

31 (22 to 43)

  

rs7515776

g.+48A/T

      
 

A/A

46 (44 to 48)

7.02

9.7e-3

67 (61 to 74)

1.76

0.17

 

A/T

54 (50 to 59)

  

77 (62 to 95)

  
 

T/T

48 (33 to 70)

  

119 (62 to 231)

  

rs1538372

g.+1219G/A

      
 

G/G

56 (53 to 59)

46.93

< 2.2e-16

78 (68 to 89)

7.58

6.5e-3

 

G/A

43 (41 to 45)

  

69 (60 to 78)

  
 

A/A

34 (30 to 37)

  

45 (34 to 57)

  

rs2071579

g.+2117G/C

      
 

C/C

56 (52 to 59)

21.91

6.5e-10

72 (62 to 84)

6.31

0.0022

 

C/G

47 (44 to 49)

  

76 (67 to 86)

  
 

G/G

39 (36 to 42)

  

52 (43 to 62)

  

rs880633

g.+2950C/T

      
 

C/C

56 (52 to 59)

21.91

6.7e-10

72 (62 to 84)

6.31

0.0022

 

C/T

47 (44 to 49)

  

76 (67 to 86)

  
 

T/T

39 (36 to 42)

  

52 (43 to 62)

  

Serum concentrations of YKL-40 are given as median ± 95% CI.

CHI3L1, chitinase 3-like 1 gene; CI, confidence interval; RA, rheumatoid arthritis; SNP, single nucleotide polymorphism.

https://static-content.springer.com/image/art%3A10.1186%2Far3391/MediaObjects/13075_2011_Article_3136_Fig3_HTML.jpg
Figure 3

Association between the CHI3L1 g.-131(C > G) polymorphism and age adjusted serum concentrations of YKL-40. Box plot illustrating the association in 238 patients with rheumatoid arthritis (RA) ≤ 65 years of age (P < 2.0e-16 and 605 healthy controls (P < 1.1e-8). The x-axis represents CHI3L1 g.-131(C > G) genotypes. The y-axis represents serum YKL-40, horizontal bars represents median serum YKL-40 and quartiles for patients with RA and controls.

When the g.-131 C/G was used as a covariate to determine the influence of the remaining seven SNPs on serum concentrations of YKL-40 none of the other SNPs contributed significantly to the association supporting the isolated highly significant effect of the g.-131 C/G polymorphism on serum concentrations of YKL-40 (Table 4).
Table 4

g.-131(C/G) used as a covariate to determine the influence of the remaining 7 SNP:s on s-YKL-40

SNP

CHI3L1position

P- value

rs4950928

g.-131C/G

< 2.2e-16

rs6691378

g.-1371G/A

0.21

rs10399931

g.-329C/T

0.88

rs10399805

g.-247G/A

0.19

rs7515776

g.+48A/T

0.25

rs1538372

g.+1219G/A

0.57

rs2071579

g.+2117G/C

0.72

rs880633

g.+2950C/T

0.72

CHI3L1, chitinase 3-like 1 gene; SNP, single nucleotide polymorphism, s-YKL-40 serum concentrations of YKL-40.

Haplotype analysis did not add further information as all the haplotypes associated with low serum concentrations of YKL-40 carried the g.-131G allele and no further increase in association was seen with any of the haplotypes (data not shown). LD analysis of the eight genotyped SNPs revealed that both the proximal promoter and the distal part of the gene contained blocks of high or moderate LD (Figure 4) explaining the effect of all the included polymorphisms on serum YKL-40 when analysed individually. In particular the -131 C/G polymorphism displayed moderate LD with g.-329C/T (R2 0.78) indicating that the effect on serum concentrations of YKL-40 with g.-329C/T is caused by LD. These findings are in line with CEU HapMap data (Figure 5).
https://static-content.springer.com/image/art%3A10.1186%2Far3391/MediaObjects/13075_2011_Article_3136_Fig4_HTML.jpg
Figure 4

Linkage disequilibrium in the Danish control individuals (R 2 values). SNPs are defined by RefSNP number. SNP: single nucleotide polymorphism.

https://static-content.springer.com/image/art%3A10.1186%2Far3391/MediaObjects/13075_2011_Article_3136_Fig5_HTML.jpg
Figure 5

Linkage disequlibrium between SNPs in the CHI3L1 gene in the CEU HapMap population. All SNPs are defined by RefSNP number. CHI3L1, chitinase 3-like 1 gene; SNP, single nucleotide polymorphism.

To investigate the association of the eight SNPs with case-control status, allelic and genotypes were tested for association with RA using Fishers exact test. No association was found with alleles or genotypes for any of the eight SNPs (Table 5) indicating that these SNPs do not confer risk to the development of RA itself.
Table 5

Association of CHI3L1 SNPs with rheumatoid arthritis

SNP

SNP location

CHI3L1position

RA (n = 308)

Controls (n = 605)

P-value

rs6691378

Promoter

g.-1371G/A

n (%)

n (%)

 
  

G/G

249 (80.8)

468 (77.4)

0.48

  

G/A

56 (18.2)

130 (21.5)

 
  

A/A

3 (1.0)

7 (1.2)

 
  

Allele

   
  

G

554 (89.9)

1,066 (88.1)

0.24

  

A

62 (10.1)

144 (11.9)

 

rs10399931

Promoter

g.-329C/T

   
  

C/C

175 (56.8)

360 (59.5)

0.61

  

C/T

116 (37.7)

208 (34.4)

 
  

T/T

17 (5.5)

37 (6.1)

 
  

Allele

   
  

C

466 (75.6)

928 (76.7)

0.62

  

T

150 (24.4)

282 (23.3)

 

rs10399805

Promoter

g.-247G/A

   
  

G/G

248 (80.5)

464 (76.7)

0.38

  

G/A

56 (18.2)

134 (22.1)

 
  

A/A

4 (1.3)

7 (1.2)

 
  

Allele

   
  

G

552 (89.6)

1,062 (87.8)

0.25

  

A

64 (10.4)

148 (12.2)

 

rs4950928

Promoter

g.-131C/G

   
  

C/C

190 (61.7)

390 (64.5)

0.61

  

C/G

101 (32.8)

179 (29.5)

 
  

G/G

17 (5.5)

36 (6.0)

 
  

Allele

   
  

C

481 (78.1)

959 (79.3)

0.56

  

G

135 (21.9)

251 (20.7)

 

rs7515776

Intron 1/ exon 1

g.+48A/T

   
  

A/A

247 (80.2)

467 (77.2)

0.20

  

A/T

55 (17.9)

132 (21.8)

 
  

T/T

6 (1.9)

6 (1.0)

 
  

Allele

   
  

A

549 (89.1)

1,066 (88.1)

0.52

  

T

67 (10.9)

144 (11.9)

 

rs1538372

Intron 2/ exon 3

g.+1219G/A

   
  

G/G

137 (44.5)

277 (45.8)

0.93

  

G/A

137 (44.5)

262 (43.3)

 
  

A/A

34 (11.0)

66 (10.9)

 
  

Allele

   
  

G

411 (66.7)

816 (67.4)

0.76

  

A

205 (33.3)

394 (32.6)

 

rs2071579

intron 4/ exon 4

g.+2117G/C

   
  

C/C

63 (20.5)

123 (20.3)

0.89

  

C/G

149 (48.4)

302 (49.9)

 
  

G/G

96 (31.2)

180 (29.8)

 
  

Allele

   
  

C

275 (44.6)

548 (45.3)

0.79

  

G

341 (55.4)

662 (54.7)

 

rs880633

exon 5

g.+2950C/T

   
  

C/C

96 (31.2)

180 (29.8)

0.88

  

C/T

149 (48.4)

303 (50.0)

 
  

T/T

63 (20.4)

122 (20.2)

 
  

Allele

   
  

C

341 (55.4)

663 (54.8)

0.82

  

T

275 (44.6)

547 (45.2)

 

CHI3L1, chitinase 3-like 1 gene; RA, rheumatoid arthritis; SNP, single nucleotide polymorphism.

The high producer genotypes were not more frequent in the IgM-RF positive subgroup and no difference was found in geno- or phenotype distribution between seropositive and seronegative patients with RA (data not shown).

Discussion

This study aimed to investigate eight polymorphic sites in the CHI3L1 gene with possible functional properties in both patients with RA and healthy individuals. We focused on the g.-131(C > G) allele and closely related polymorphisms described in Caucasian populations [26, 3639, 43, 44, 46]. The g.1219(G > A) polymorphism was also included as one study reported an individual functional property of this polymorphism [43]. Serum concentrations of YKL-40 were strongly associated with age and case-control status. After adjustment of the serum concentrations of YKL-40 for these two variables, serum YKL-40 was found to be significantly associated with SNPs in the CHI3L1 gene. The strongest association was with the g.-131(C > G) promoter polymorphism. The association of serum YKL-40 with the remaining SNPs could be explained by LD. Our findings indicate that serum concentrations of YKL-40 are under the influence of genetic variability in the CHI3L1 gene in both patients with RA and healthy controls, and the effect of genotypes seems to be the same in both groups. Though our results indicate that CHI3L1 polymorphisms are not involved in the pathogenesis of RA, we do not know if high producer genotypes results in a more severe clinical phenotype.

Several other studies have suggested the g.-131(C > G) is a strong candidate for a functional promoter polymorphism influencing the serum concentrations of YKL-40 [36, 4245]. The promoter SNP g.-131(C > G) in the CHI3L1 gene was associated with elevated serum YKL-40, asthma, bronchial hyper responsiveness and pulmonary function [44, 45], and with elevated serum YKL-40 and the severity of hepatitis C virus-induced liver fibrosis [43]. This indicates a functional role of YKL-40 in these diseases. An association is also found between schizophrenia and haplotypes within the promoter region of the CHI3L1 gene suggesting that polymorphisms in an area starting from base pair position -180 could have functional properties [36, 42]. Our findings support these earlier studies.

Zhao et al. [36] investigated Chinese patients with schizophrenia and found lower activity of the transcription factor MYC/MAX and decreased CHI3L1 gene expression related to the low frequency G allele for the g.-131(C > G) SNP. Ober et al. [44] studied 443 patients with asthma and 491 healthy controls from a genetically preserved group of Americans of European descent. They found that serum concentrations of YKL-40 were associated with many alleles in the promoter region including the g.-131(C > G) and g.-329(C < T) polymorphisms. This supports the g.-131(C > G) polymorphism as a site of genetic regulation in both healthy controls and patients with asthma. Conclusions were complicated by strong LD in the promoter region in the population studied. They also showed a strong association to the g.-1219(G > A) polymorphisms, which was not in LD with the promoter polymorphisms. This individual effect on serum YKL-40 with g.-1219(G > A) was not supported in our study as we found this phenomena related to LD in the Danish population. In contrast, Sohn et al. [40] demonstrated a functional effect of the g.-247 (G > A) polymorphisms in a study of 295 atopic children and 180 healthy controls from a Korean population. They concluded that this polymorphism was responsible for most of the genetic effects on YKL-40 production, and that the g.-131C allele was associated to low promoter activity. These results are complicated by the fact that the g.-131(C > G) and g.-247(G > A) polymorphisms showed no LD in the Asian populations, contradictory to our finding which suggest a high degree of LD in this part of the proximal promoter. In the Danish population the region on chromosome 1 bearing the g.-131(C > G) polymorphism was in strong LD, illustrated by the occurrence of just 8 frequent haplotypes (f > 1%). The g.-131(G > C) allele was found to be in LD with several other loci in the CHI3L1 gene, and three haplotypes could be defined as low producer haplotypes, all including the g.-131G allele. It must be emphasized that ethnicity seems to play an important role in the genetic regulation of YKL-40 production, and results from our and similar studies can only be considered valid in ethnically similar populations. Further studies in different populations are awaited.

Serum concentrations of YKL-40 increased with increasing age in both healthy controls and patients with RA making age a possible confounding variable. The cause of this remains unknown, but the phenomena has been explained by a higher level of general inflammation and apoptosis in the elderly, which is well known for other inflammatory mediators [53]. Similar increases in plasma YKL-40 with age have recently been described in a large group of 8,899 subjects from the general Danish population [20]. We initially decided to fit a non-linear model to explain the effect of age on serum concentrations of YKL-40. Our control group did not include any persons above the age of 65, but below this age we were able to fit a linear model explaining the relationship between age and serum concentrations of YKL-40. This supports the findings by Kruit et al. [38], who suggested a linear relationship between age and serum YKL-40. In the patients with RA it seems as if serum concentrations of YKL-40 rise more rapidly with age above age 65, indicating that elevated serum YKL-40 in this age group needs careful interpretation. It is possible that high serum YKL-40 is associated to co-morbidity or a latent malignant disease [10, 2026, 53, 54]

It remains unknown whether high serum YKL-40 affects a person:s risk of autoimmune disease in the long term. YKL-40 expression is stimulated by the inflammatory cytokines TNF-α, IL-6 [30] and IL-1β, whereas YKL-40 inhibits cellular responses induced by IL-1 and TNF-α, suggesting an autocrine feed-back mechanism [9, 28]. YKL-40 is strongly expressed by macrophages in the synovial membrane of RA patients possibly activated by a pro-inflammatory IFNγ-mediated immune response, and elevated YKL-40 can stimulate local production of anti-inflammatory IL-10 [32]. In inflammatory diseases such as RA, the excessive YKL-40 production may also have the opposite effect stimulating a continuous pro-inflammatory state and stimulation of VEGF and angiogenesis [3234].

Conclusions

In conclusion, this study reports a strong association between the g.-131(C > G) allele and serum concentrations of YKL-40 in both patients with RA and healthy controls. Our findings indicate that the g.-131(C > G) polymorphism is the main contributor to the inter-individual variation of serum YKL-40 in Caucasian patients with RA, and that the effect of other polymorphic sites in this region is related to a high degree of LD in this area of the genome.

Abbreviations

CHI3L1

chitinase 3-like 1 gene

CV: 

coefficient of variations

LD: 

linkage disequilibrium

RA: 

rheumatoid arthritis

RF: 

IgM rheumatoid factor

RS: 

RefSNP number

RT-PCR: 

real-time polymerase chain reaction

SNP: 

single nucleotide polymorphism

Declarations

Acknowledgements

The study was funded by grants from The Obel Family Foundation, Badehusvej 59000 Aalborg, Denmark

Authors’ Affiliations

(1)
Department of Clinical Immunology, Aalborg Hospital, Aarhus University Hospital
(2)
Department of Haematology, Aalborg Hospital Science and Innovation Center, Aarhus University Hospital
(3)
Department of Cardiology, Center for Cardiovascular Research, Aalborg Hospital, Aarhus University Hospital
(4)
Department of Rheumatology, Copenhagen University Hospital, Hvidovre and Glostrup
(5)
Department of Medicine, Copenhagen University Hospital
(6)
Department of Oncology, Copenhagen University Hospital

References

  1. Klareskog L, Catrina AI, Paget S: Rheumatoid arthritis. Lancet. 2009, 373: 659-672. 10.1016/S0140-6736(09)60008-8.View ArticlePubMedGoogle Scholar
  2. Firestein GS: Evolving concepts of rheumatoid arthritis. Nature. 2003, 423: 356-361. 10.1038/nature01661.View ArticlePubMedGoogle Scholar
  3. Oliver JE, Worthington J, Silman A: Genetic epidemiology of rheumatoid arthritis. Curr Opin Rheumatol. 2006, 18: 141-146. 10.1097/01.bor.0000209425.84775.23.View ArticlePubMedGoogle Scholar
  4. Robert M, Plengea B: Recent progress in rheumatoid arthritis genetics:one step towards improved patient care. Curr Opin Rheumatol. 2009, 21: 262-271. 10.1097/BOR.0b013e32832a2e2d.View ArticleGoogle Scholar
  5. Karlson EW: Associations between HLA, PTPN22, CTLA4 genotypes and RA phenotypes of autoantibody status, age at diagnosis, and erosions in a large cohort study. Ann Rheum Dis. 2008, 67: 358-363.View ArticlePubMedGoogle Scholar
  6. Rehli M, Krause SW, Andreesen R: Molecular characterization of the gene for human cartilage gp-39 (CHI3L1), a member of the chitinase protein family and marker for late stages of macrophage-differentiation. Genomics. 1997, 43: 221-225. 10.1006/geno.1997.4778.View ArticlePubMedGoogle Scholar
  7. Johansen JS, Williamson MK, Rice JS, Price PA: Identification of proteins secreted by human osteoblastic cells in culture. J Bone Miner Res. 1992, 7: 501-512.View ArticlePubMedGoogle Scholar
  8. Johansen JS: Studies on serum YKL-40 as a biomarker in diseases with inflammation, tissue remodelling, fibroses and cancer. Dan Med Bull. 2006, 53: 172-209.PubMedGoogle Scholar
  9. Eurich K, Segawa M, Toei-Shimizu S, Mizoguchi E: Potential role of chitinase 3-like-1 in inflammationassociated carcinogenic changes of epithelial cells. World J Gastroenterol. 2009, 15: 5249-5259. 10.3748/wjg.15.5249.View ArticlePubMedGoogle Scholar
  10. Johansen JS, Schultz NA, Jensen BV: Plasma YKL-40: a potential new cancer biomarker?. Future Oncol. 2009, 5: 1065-1082. 10.2217/fon.09.66.View ArticlePubMedGoogle Scholar
  11. Verheijden GF, Rijnders AW, Bos E, Coenen-de Roo CJ, van Staveren CJ, Miltenburg AM, Meijerink JH, Elewaut D, de Keyser F, Veys E, Boots AM: Human cartilage glycoprotein-39 as a candidate autoantigen in rheumatoid arthritis. Arthritis Rheum. 1997, 40: 1115-1125. 10.1002/art.1780400616.View ArticlePubMedGoogle Scholar
  12. Cope AP, Patel SD, Hall F, Congia M, Hubers HA, Verheijden GF, Boots AM, Menon R, Trucco M, Rijnders AW, Sønderstrup G: T cell responses to a human cartilage autoantigen in the context of rheumatoid arthritis-associated and nonassociated HLA-DR4 alleles. Arthritis Rheum. 1999, 42: 1497-1507. 10.1002/1529-0131(199907)42:7<1497::AID-ANR25>3.0.CO;2-#.View ArticlePubMedGoogle Scholar
  13. Vos K, Miltenburg AM, van Meijgaarden KE, van den Heuvel M, Elferink DG, van Galen PJ, van Hogezand RA, van Vliet-Daskalopoulou E, Ottenhoff TH, Breedveld FC, Boots AM, de Vries RR: Cellular immune response to human cartilage glycoprotein-39 (HC gp-39)-derived peptides in rheumatoid arthritis and other inflammatory conditions. Rheumatology. 2000, 39: 1326-1331. 10.1093/rheumatology/39.12.1326.View ArticlePubMedGoogle Scholar
  14. Patil NS, Hall FC, Drover S, Spurrell DR, Bos E, Cope AP, Sonderstrup G, Mellins ED: Autoantigenic HCgp39 epitopes are presented by the HLADM-dependent presentation pathway in human B cells. J Immunol. 2001, 166: 33-41.View ArticlePubMedGoogle Scholar
  15. Kavanaugh A, Genovese M, Baughman J, Kivitz A, Bulpitt K, Olsen N, Weisman M, Matteson E, Furst D, van Vollenhoven R, Anderson J, Cohen S, Wei N, Meijerink J, Jacobs C, Mocci S: Allele and antigen-specific treatment of rheumatoid arthritis: a double blind, placebo controlled phase 1 trial. J Rheumatol. 2003, 30: 449-454.PubMedGoogle Scholar
  16. Johansen JS, Christoffersen P, Møller S, Price PA, Henriksen JH, Garbarsch C, Bendtsen F: Serum YKL-40 is increased in patients with hepatic fibrosis. J Hepatol. 2000, 32: 911-920. 10.1016/S0168-8278(00)80095-1.View ArticlePubMedGoogle Scholar
  17. Koutroubakis IE, Petinaki E, Dimoulios P, Vardas E, Roussomoustakaki M, Maniatis AN, Kouroumalis EA: Increased serum levels of YKL-40 in patients with inflammatory bowel disease. Int J Colorectal Dis. 2003, 18: 254-259.PubMedGoogle Scholar
  18. Harvey S, Whaley J, Eberhardt K: The relationship between serum levels of YKL-40 and disease progression in patients with early rheumatoid arthritis. Scand J Rheumatol. 2000, 29: 391-393. 10.1080/030097400447606.View ArticlePubMedGoogle Scholar
  19. Matsumoto T, Tsurumoto T: Serum YKL-40 levels in rheumatoid arthritis: correlations between clinical and laborarory parameters. Clin Exp Rheumatol. 2001, 19: 655-660.PubMedGoogle Scholar
  20. Johansen JS, Bojesen SE, Tybjaerg-Hansen A, Mylin AK, Price PA, Nordestgaard BG: Plasma YKL-40 and total and disease-specific mortality in the general population. Clin Chem. 2010, 56: 1580-1591. 10.1373/clinchem.2010.146530.View ArticlePubMedGoogle Scholar
  21. Kjaergaard AD, Bojesen SE, Johansen JS, Nordestgaard BG: Elevated plasma YKL-40 levels and ischemic stroke in the general population. Ann Neurol. 2010, 68: 672-680. 10.1002/ana.22220.View ArticlePubMedGoogle Scholar
  22. Kastrup J, Johansen JS, Winkel P, Hansen JF, Hildebrandt P, Jensen GB, Jespersen CM, Kjøller E, Kolmos HJ, Lind I, Nielsen H, Gluud C, CLARICOR Trial Group: High serum YKL-40 concentration is associated with cardiovascular and all-cause mortality in patients with stable coronary artery disease. Eur Heart J. 2009, 30: 1066-1072. 10.1093/eurheartj/ehp049.View ArticlePubMedGoogle Scholar
  23. Nielsen AR, Erikstrup C, Johansen JS, Fischer CP, Plomgaard P, Krogh-Madsen R, Taudorf S, Lindegaard B, Pedersen BK: Plasma YKL-40: a BMI-independent marker of type 2 diabetes. Diabetes. 2008, 57: 3078-3082. 10.2337/db08-0182.View ArticlePubMedGoogle Scholar
  24. Létuvé S, Kozhich A, Arouche N, Grandsaigne M, Reed J, Dombret MC, Kiener PA, Aubier M, Coyle AJ, Pretolani M: YKL-40 is elevated in patients with chronic obstructive pulmonary disease and activates alveolar macrophages. J Immunol. 2008, 181: 5167-5173.View ArticlePubMedGoogle Scholar
  25. Fontana RJ, Dienstag JL, Bonkovsky HL, Sterling RK, Naishadham D, Goodman ZD, Lok AS, Wright EC, Su GL, HALT-C Trial Group: Serum fibrosis markers are associated with liver disease progression in non-responder patients with chronic hepatitis C. Gut. 2010, 59: 1401-1409. 10.1136/gut.2010.207423.View ArticlePubMedGoogle Scholar
  26. Chupp GL, Lee CG, Jarjour N, Shim YM, Holm CT, He S, Dziura JD, Reed J, Coyle AJ, Kiener P, Cullen M, Grandsaigne M, Dombret MC, Aubier M, Pretolani M, Elias JA: A chitinase-like protein in the lung and circulation of patients with severe asthma. N Engl J Med. 2007, 357: 2016-2027. 10.1056/NEJMoa073600.View ArticlePubMedGoogle Scholar
  27. Recklies AD, Ling H, White C: Inflammatory cytokines induce production of chitinase-3-like protein 1 by articular chondrocytes. J Biol Chem. 2005, 280: 41213-41221. 10.1074/jbc.M510146200.View ArticlePubMedGoogle Scholar
  28. Ling H, Recklies AD: The chitinase 3-like protein human cartilage glycoprotein 39 inhibits cellular responses to the inflammatory cytokines interleukin-1 and tumour necrosis factor-alpha. Biochem J. 2004, 380: 651-659. 10.1042/BJ20040099.View ArticlePubMedGoogle Scholar
  29. Mizoguchi E, Mizoguchi A: Is the sugar always sweet in intestinal inflammation?. Immunol Res. 2007, 37: 47-60. 10.1007/BF02686089.View ArticlePubMedGoogle Scholar
  30. Nielsen AR, Plomgaard P, Krabbe KS, Johansen JS, Pedersen BK: IL-6, but not TNF-α, increases plasma YKL-40 in human subjects. Cytokine. 2011, 55: 152-155. 10.1016/j.cyto.2011.03.014.View ArticlePubMedGoogle Scholar
  31. Recklies AD, White C, Ling H: The chitinase 3-like protein human cartilage glycoprotein 39 (HC-gp39) stimulates proliferation of human connective-tissue cells and activates both extracellular signal-regulated kinase- and protein kinase B mediated signalling pathways. Biochem J. 2002, 365: 119-126. 10.1042/BJ20020075.View ArticlePubMedGoogle Scholar
  32. Lee CG, Da Silva CA, Dela Cruz CS, Ahangari F, Ma B, Kang MJ, He CH, Takyar S, Elias JA: Role of chitin and chitinase/chitinase-like proteins in inflammation, tissue remodeling and injury. Ann Rev Physiol. 2011, 73: 479-501. 10.1146/annurev-physiol-012110-142250.View ArticleGoogle Scholar
  33. Faibish M, Francescone R, Bentley B, Yan W, Shao R: A YKL-40-neutralizing antibody blocks tumor angiogenesis and progression: a potential therapeutic agent in cancers. Mol Cancer Ther. 2011, 10: 742-751. 10.1158/1535-7163.MCT-10-0868.View ArticlePubMedGoogle Scholar
  34. Francescone RA, Scully S, Faibish M, Taylor SL, Oh D, Moral L, Yan W, Bentley B, Shao R: Role of YKL-40 in the angiogenesis, radioresistance, and progression of glioblastoma. J Biol Chem. 2011, 286: 15332-15343. 10.1074/jbc.M110.212514.View ArticlePubMedGoogle Scholar
  35. Rehli M, Niller HH, Ammon C: Transcriptional regulation of CHI3L1, a marker gene for late stages of macrophage differentiation. J Biol Chem. 2003, 278: 44058-44067. 10.1074/jbc.M306792200.View ArticlePubMedGoogle Scholar
  36. Zhao X, Tang R, Gao B: Functional variants in the promoter region of chitinase 3-like 1 (CHI3L1) and susceptibility to schizophrenia. Am J Hum Genet. 2007, 80: 12-18. 10.1086/510438.View ArticlePubMedGoogle Scholar
  37. Rathcke CN, Holmkvist J, Jørgensen T, Borch-Johnsen K, Hansen T, Pedersen OB, Vestergaard H: Variation in CHI3L1 in relation to type 2 diabetes and related quantitative traits. Plos One. 2009, 4: e5469-10.1371/journal.pone.0005469.View ArticlePubMedGoogle Scholar
  38. Kruit A, Grutters JC, Ruven HJ, van Moorsel CC, van den Bosch JM: CHI3L1 gene polymorphism is associated with serum levels of YKL-40, a novel sarcoidosis marker. Respir Med. 2007, 101: 1563-1571. 10.1016/j.rmed.2006.12.006.View ArticlePubMedGoogle Scholar
  39. Rathcke CN, Holmkvist J, Husmoen LL, Hansen T, Pedersen O, Vestergaard H, Linneberg A: Association of polymorphisms of the CHI3L1 gene with asthma and atopy: a populations-based study of 6514 danish adults. PLoS One. 2009, 4: e6106-10.1371/journal.pone.0006106.View ArticlePubMedGoogle Scholar
  40. Sohn MH, Lee JH: Genetic variation in the promoter region of chitinase 3-like 1 is associated with atopy. Am J Respir Crit Care Med. 2009, 179: 449-456. 10.1164/rccm.200809-1422OC.View ArticlePubMedGoogle Scholar
  41. Ohi K, Hashimoto R, Yasuda Y: The chitinase 3-like 1 gene and schizophrenia: evidence from a multi-center case-control study and meta-analysis. Schizophr Res. 2010, 116: 126-132. 10.1016/j.schres.2009.12.002.View ArticlePubMedGoogle Scholar
  42. Yang MS, Morris DW, Donohoe G, Kenny EO, Dushalaine CT, Schwaiger S, Nangle JM, Clarke S, Scully P, Quinn J, Meagher D, Baldwin P, Crumlish NO, Callaghan E, Waddington JL, Gill M, Corvin A: Chitinase-3-like 1 (CHI3L1) gene and schizophrenia: genetic association and a potential functional mechanism. Biol Psychiatry. 2008, 64: 98-103. 10.1016/j.biopsych.2007.12.012.View ArticlePubMedGoogle Scholar
  43. Berres ML, Papen S, Pauels K, Schmitz P, Zaldivar MM, Hellerbrand C, Mueller T, Berg T, Weiskirchen R, Trautwein C, Wasmuth HE: A functional variation in CHI3L1 is associated with severity of liver fibrosis and YKL-40 serum levels in chronic hepatitis C infection. J Hepatol. 2009, 50: 370-376. 10.1016/j.jhep.2008.09.016.View ArticlePubMedGoogle Scholar
  44. Ober C, Tan Z, Sun Y, Possick JD, Pan L, Nicolae R, Radford S, Parry RR, Heinzmann A, Deichmann KA, Lester LA, Gern JE, Lemanske RF, Nicolae DL, Elias JA, Chupp GL: Effect of variation in CHI3L1 on serum YKL-40 level, risk of asthma, and lung function. N Engl J Med. 2008, 358: 1682-1691. 10.1056/NEJMoa0708801.View ArticlePubMedGoogle Scholar
  45. Ober C, Chupp GL: The chitinase and chitinase-like proteins: a review of genetic and functional studies in asthma and immune-mediated diseases. Curr Opin Allergy Clin Immun. 2009, 9: 401-408. 10.1097/ACI.0b013e3283306533.View ArticleGoogle Scholar
  46. Srivastava SK, Antal P, Gál J: Lack of evidence for association of two functional SNPs of CHI3L1 gene (HC-gp39) with rheumatoid arthritis. Rheumatol Int. 2011, 31: 1003-1007. 10.1007/s00296-010-1396-3.View ArticlePubMedGoogle Scholar
  47. Arnett FC, Edworthy SM, Bloch DA, McShane DJ, Fries JF, Cooper NS, Healey LA, Kaplan SR, Liang MH, Luthra HS: The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum. 1988, 31: 315-324. 10.1002/art.1780310302.View ArticlePubMedGoogle Scholar
  48. Hetland ML: DANBIO - powerful research database and electronic patient record. Rheumatology. 2011, 50: 69-77. 10.1093/rheumatology/keq309.View ArticlePubMedGoogle Scholar
  49. Harvey S, Weisman MO, Dell J: Chondrex: new marker of joint disease. Clin Chem. 1998, 44: 509-516.PubMedGoogle Scholar
  50. Yong Y, He L: SHEsis, a powerful software platform for analyses of linkage disequilibrium, haplotype construction, and genetic association at polymorphism loci. Cell Research. 2005, 15: 97-98. 10.1038/sj.cr.7290272.View ArticleGoogle Scholar
  51. R Development Core Team: R A language and environment for statistical computing. R Foundation for Statistical Computing Vienna, Austria, ISBN 3-900. 2004, [http://www.r-project.org]Google Scholar
  52. Harrell F: Regression Modeling Strategies. 2001, New York: Springer VerlagView ArticleGoogle Scholar
  53. Johansen JS, Pedersen AN, Schroll M: High serum YKL-40 level in a cohort of octogenarians is associated with increased risk of all-cause mortality. Clin Exp Immunol. 2008, 151: 260-266.View ArticlePubMedGoogle Scholar
  54. Johansen JS, Bojesen SE, Mylin AK: Elevated plasma YKL-40 predicts increased risk of gastrointestinal cancer and decreased survival after any cancer diagnosis in the general population. J Clin Oncol. 2009, 27: 572-578.View ArticlePubMedGoogle Scholar

Copyright

© Nielsen et al.; licensee BioMed Central Ltd. 2011

This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.