Thorax 2005;60:320-325
© 2005 BMJ Publishing
Group Ltd & British Thoracic Society
CYSTIC FIBROSIS |
1 Department for Human Genetics, KULeuven,
Herestraat 49, O&N6, 3000 Leuven, Belgium
2 Institute of
Molecular Genetics and Genetic Engineering, Vojvode Stepe 444a, Belgrade, Serbia
and Montenegro
3 Department of Pediatrics, UZ Gasthuisberg,
Herestraat 49, 3000 Leuven, Belgium
4 Charles University Prague
and University Hospital Motol, V Uvalu 84, CZ 15006 Prague, Czech Republic
Correspondence to:
Dr H Cuppens
Department for Human Genetics, KULeuven, Herestraat 49, O&N6, 3000
Leuven, Belgium; harry.cuppens{at}med.kuleuven.ac.be
Received 19 March 2004
Accepted for
publication 25 October 2004
|
|
ABSTRACT |
|---|
Methods: Three polymorphic loci in the promoter (851c/t,
308g/a, 238g/a) and one polymorphic locus in intron
1 (+691g ins/del) of the TNF
gene were typed by a single nucleotide primer extension
assay in CF patients and healthy controls. Spirometric data and first
age of infection with Pseudomonas aeruginosa were collected
retrospectively from patients medical records.
Results: An association was found between the TNF
+691g ins/del polymorphic locus and severity of CF
lung disease. Patients heterozygous for +691g ins and +691g
del were more likely to have better pulmonary function (mean (SD)
forced expiratory volume in 1 second (FEV1) 79.7 (12.8)%
predicted) than patients homozygous for +691g ins (mean (SD)
FEV1 67.5 (23.0)% predicted; p = 0.008, mean difference
12.2%, 95% CI 3.5 to 21.0). Also, patients heterozygous for +691g
ins and +691g del were more likely to have an older first
age of infection with P aeruginosa (mean (SD) 11.4 (6.0)
years) than patients homozygous for +691g ins (mean (SD) 8.3
(4.6) years; p = 0.018, mean difference 3.1 years, 95% CI 0.5 to
5.6). An association was also found with the 851c/t
polymorphic locus. In the group of patients with more severe
FEV1% predicted, a higher proportion of patients were
homozygous for the 851c allele than in the other group of
patients (p = 0.04, likelihood ratio
2, odds ratio = 2.4).
Conlusion: TNF
polymorphisms are associated with the severity of CF lung
disease in Czech and Belgian patients with CF.
Abbreviations: CF, cystic fibrosis; CFTR, cystic fibrosis
transmembrane conductance regulator; FEV1, forced expiratory volume
in 1 second; MBL, mannose binding lectin; TNF
, tumour necrosis factor
Keywords: tumour necrosis factor
; cystic fibrosis; genetics
Cystic fibrosis (CF) is a common lethal autosomal recessive disease affecting the white population with an incidence of about 1 in 2500.1,2 Certain phenotypes of CF are clearly associated with different CF transmembrane conductance regulator (CFTR) genotypes such as pancreatic insufficiency, but pulmonary disease is strongly influenced by other factors.3 There are many genes that may affect the CF phenotype. These genes may affect secondary symptoms, such as genes encoding proteins involved in lung defence (immune system, inflammation processes, protection against oxidative stress) or genes encoding proteins which directly or indirectly interact with CFTR. Phenotypic differences among patients with the same CFTR genotype are most likely caused by differences in the function of such proteins and/or by environmental factors. Previous studies have shown an association between variant alleles and the severity of the CF phenotypefor example, it was found that polymorphisms in the mannose binding lectin (MBL) protein, which result in low or no functional MBL, contribute to more severe CF lung disease.4,5
Tumour necrosis factor
(TNF
), an endogenous factor of the lung involved in host
defence,6 is therefore a possible modulator of CF
pulmonary disease severity. TNF
is a proinflammatory cytokine largely produced by
macrophages,7 whose gene is located on chromosome
6p21.3.8 It is prothrombotic, promotes leucocyte adhesion
and migration, modulates haematopoiesis and lymphocyte
development, induces other cytokines, and has an important role in
macrophage activation and the immune response in tissues.9
Even though there are many polymorphisms in the promoter region of
TNF
, much attention has been given to the TNF
308g/a polymorphic locus. In a previous study of CF
patients from a UK CF clinical centre, an association was found
between the 308a allele and a more severe CF
phenotype.10 A modulating role of a genetic factor with
disease, based on association studies, should however only be treated
as tentative until the finding of an association has been replicated
in other studies. Moreover, we broadened the search with additional
polymorphisms for an association between TNF
and severity of lung function in CF patients. We therefore
investigated whether the TNF
308g/a polymorphic locus, together with the
851c/t, 238g/a promoter polymorphic loci and the
polymorphic locus +691g ins/del in intron 111, are
associated with pulmonary function, body mass index (BMI), and
susceptibility to Pseudomonas aeruginosa infection in Belgian
and Czech patients with CF.
|
|
METHODS |
|---|
For Hardy-Weinberg equilibrium testing, all unrelated CF patients
were included. For the actual association study of TNF
polymorphisms with the CF pulmonary phenotype, only 113 CF
patients (41 Belgian and 72 Czech) between the ages of 12 and 15
years (mean age 13.4) were included; 49% were male. Indeed, a
modulating gene may only become penetrant by age so it will only be
detected in association studies using a well defined group of CF
patients of a particular age or age range. All clinical data were
retrieved from the patients case records. The BMI was measured
and percentage predicted values of forced expiratory volume in
1 second (FEV1) were calculated according to Knudson et
al.12 Statistical analysis of lung function was
performed on the last FEV1 % predicted value between the
ages of 12 and 15 years.
The age of first infection with P aeruginosa was available from the case records of 119 of the total group of 180 CF patients. Patients were followed at the outpatient clinic at a minimum of 3 monthly intervals. At every visit, sputum samples (or throat swabs in patients who did not expectorate) were taken. The age at which P aeruginosa was isolated for the first time is referred to as first infection with P aeruginosa.
The association studies were approved by the ethical committees of both universities.
PCR and single nucleotide
extension assay
We developed multiplex single nucleotide
primer extension assays in which the TNF
single nucleotide polymorphisms (SNPs) were included. This
study is part of a much larger study involving semi-high throughput
analysis of tens of SNPs, so gap filling experiments were needed for
the samples in which a genotype could not be obtained for the first
time. However, from a practical point of view, complete gap filling
is not possible in such a large study because the DNA of some samples
becomes exhausted. Thus, the number of genotypes/alleles obtained for
each SNP may differ between different SNPs. For the single
nucleotide primer extension assay, the SNP detection kit (SNaPshot;
Applied Biosystems) was used. The primers used for amplification
of the genomic regions covering the SNPs in TNF
and the single nucleotide extension oligonucleotides are
shown in tables S1 and S2 respectively, available on the
Thorax website at http://www.thoraxjnl.com/supplemental.
Statistical
analysis
Database management and statistical analysis were
performed with the SAS statistical package release 8.01 (SAS
Institute Inc, Cary, NC, USA) and the SYSTAT package release 7.0
(SPSS Inc, Chicago, IL, USA). Statistical tests were considered
significant when their type I error was less than 0.05. Differences
between the mean FEV1 % predicted values of Czech and
Belgian populations were assessed using a two-sample t test
after checking whether the variances of the populations were equal or
not. Comparison between one independent and one dependent qualitative
variables were assessed using the
2 likelihood ratio test or by the Fisher exact
probability test when more than 20% of the cells had an expected
count of less than 5. When the independent variables were qualitative
and the dependent one was quantitative, a parametric analysis of
variance (ANOVA) was done. Pearsons correlation coefficient was
calculated to test the goodness of fit of one independent
quantitative value and one dependent quantitative value in a linear
regression model.
|
|
RESULTS |
|---|
FEV1 % predicted values of CF patients decline with age13 and females have a steeper decline in pulmonary function than males.14 FEV1 % predicted values calculated according to Knudson are corrected for sex and age.12 Nevertheless, we determined whether the FEV1 % predicted values were truly independent of sex and age in our age selected cohort of patients and found no influence of age (Pearson correlation coefficient r2 = 0.006) or sex (p = 0.21, mean difference 5.4%, 95% CI 13.7 to 3.0, two-sample t test) on FEV1 % predicted.
We then determined whether the values of FEV1 % predicted in
CF patients were dependent on the distribution of TNF
genotypes. Given that a modifier effect may only become
penetrant by age, as has been shown for
MBL2,5,15 only CF patients aged 1215 years
were included in the actual association study. The mean
FEV1 % of the age selected cohort of CF patients was
approximately 70% predicted, so patients with an FEV1 less
than 70% predicted were compared with those with an FEV1
of more than 70% predicted. The distribution of +691g ins/del
was found to be significantly different between the two groups, with
a higher proportion of those in whom the FEV1 was more
than 70% predicted being heterozygous for the two alleles at the
+691g ins/del locus (table 1
). More importantly, the mean FEV1 %
predicted values differed significantly between the two groups,
patients with the +691g ins/+691g del genotype having a higher
mean FEV1 % predicted value than those with the +691g
ins/+691g ins genotype (table 2
). Both these tests were significant even when a
conservative Bonferroni correction (correction for multiple testing
for four SNPs) was performed.
The distribution of 851c/t genotypes did not differ
significantly between patients with FEV1 values above and
below 70% predicted (table 1
). However, in a dominant/recessive model, a significantly
higher proportion of patients with an FEV1 lower than 70%
predicted were homozygous for the 851c allele than in the
other group of patients (table 1
). The mean FEV1 % predicted values were
not significantly different between the three genotypes of
851c/t separately or in a dominant/recessive model (table
2
).
The 308g/a and 238g/a polymorphic loci were not
significantly associated with FEV1 values above or below
70% predicted (table 1
). Also, the mean FEV1 % predicted values
were not significantly different between the three genotypes,
separately or in a dominant/recessive model, for 308g/a
or 238g/a (table 2
).
Records of the age at first onset of P aeruginosa infection
were available for 119 patients with CF. We examined whether
the age of first infection with P aeruginosa was associated
with SNP alleles or genotypes in TNF
and found that +691g ins/+691g ins homozygotes were
colonised with P aeruginosa at a significantly earlier age
than +691g ins/+691g del heterozygotes (table 3
). The mean first age of infection with P
aeruginosa for the +691g ins/del genotypes did not differ
between the Belgian and Czech CF cohorts (p = 0.21, ANOVA),
indicating that there was no skewing to either of the two centres. No
other variants were significantly associated with the first age of
infection with P aeruginosa when the genotypes were tested
separately or in a dominant/recessive model.
We also tested whether the BMI of the age selected cohort of
patients was associated with SNP alleles or genotypes in TNF
. None of the TNF
variants was significantly associated with the BMI of these
patients (see table S4 available on the Thorax website at
http://www.thoraxjnl.com/supplemental).
|
|
DISCUSSION |
|---|
Apart from the association of +691g del with better CF pulmonary
function, this variant was also significantly associated with
an older first recorded/observed age of infection with P
aeruginosa. The most common pulmonary infection in the lungs of
patients with CF is P aeruginosa,19 and progression
of lung disease unequivocally accelerates after colonisation with
this organism.20 Therefore, the younger the first age of
infection with P aeruginosa, the more rapid the decline in
pulmonary function. Mucoid P aeruginosa is a key factor in
accelerating the decline in pulmonary function in patients with
CF.21 TNF
plays an important role in the innate resistance to P
aeruginosa and its clearance from the respiratory
tract.22,23 Polymorphisms in the TNF
gene may lead to changes in levels of TNF
16 which may, in turn, increase the ability of
the lung to clear P aeruginosa. On the other hand, increased
concentrations of TNF
have been found in lung secretions of CF
patients7,24 which may contribute directly to
neutrophil influx and elastase activity that eventually destroy the
CF lung.7 Furthermore, if a certain TNF
allele has a quantitative or qualitative effect on the
TNF
protein, it would most likely have a knock on effect on the
secretion and synthesis of other members of the cytokine
cascade.9 There is probably a balance between its activity
in protecting the immune system against pathogens and its destructive
side effects. It is feasible that this balance may be tipped in one
or other direction when a patient has a certain polymorphic variant
in an important mediator of the immune system such as the TNF
gene. Interestingly, the variability in FEV1 %
predicted values in patients with the +691g del variant is
small and in the mild range, but not the most mild (fig 1
). This may be due to the balance between the
beneficial and detrimental properties of the TNF
protein.
|
Both variants of the +691g ins/del and 851c/t polymorphic
loci were found to be associated with the degree of severity of
CF disease. Variants of the +691g ins/del locus were associated
with the severity of CF lung disease in three separate tests.
In contrast to +691g ins/del, 851c/t was associated
with pulmonary function in only one test and was not associated
with pulmonary defence. The association of +691g ins/del
variants with the severity of CF pulmonary disease therefore seems
to be more important than the 851c/t variants. Of
course, it is possible that these observed associations with CF
disease severity are caused by another mutation in linkage
disequilibrium, either in this gene or in neighbouring genes. In this
regard, TNF
alleles have been shown to be in linkage disequilibrium
with HLA allelesfor example, +691gdel is linked to
DRB*11*13 and DQB1*301.25
In a previous study by Hull and Thomson10 in which an
association was found between the 308a allele and a more
severe CF phenotype, no other TNF
polymorphisms were tested. There have also been a number of
investigations into the possible association of the TNF
308g/a polymorphic locus with susceptibility to
chronic obstructive pulmonary disease, but their conclusions were
conflicting.26 In the present study no significant association
was found between 308g/a polymorphic variants and severity
of CF disease. Furthermore, we found that the 308g allele
tended to be associated with more severe CF disease, although
not significantly. Possible explanations for these dissimilar
results might be explained by the fact that we selected only
F508del homozygotes in order to exclude disease variability
because of the CFTR genotype, while Hull and Thomson10
included patients who were either F508del homozygotes or
heterozygotes. Also, we tested more than twice as many patients, and
our patients were all aged between 12 and 15 years with a mean age of
13.4 years compared with the previous study where the patients
were all aged 8 years.
The allele frequencies of the TNF
polymorphisms in this study were very similar to those
found among French controls by Herrmann et al.11 In
that study, six different haplotypes from these polymorphisms could
be constructed. One haplotype contained the +691g del allele,
the 851c allele, and the 308a allele. The +691g
del allele is only observed on this haplotype while the
308a allele appears on this haplotype and an additional
haplotype.26 This may explain the conflicting data of the
association of TNF
308a in CF lung disease in the study by Hull and
Thomson,10 and possibly the association of TNF
308a in COPD.2730 Indeed, it is
feasible that only one of the two haplotypes on which 308a
resides is the causal one, but which of them cannot be
discriminated when typing the 308g/a SNP alone. This might
also explain why several studies have found that the rare
308a allele is associated with increased levels of TNF
while others have found that 308a does not affect
levels of TNF
.16
Our findings might be of interest from a pharmacogenetic point of
view. Our results show that particular polymorphic TNF
loci are associated with better lung function.
Anti-inflammatory agents are currently being used for the treatment
of CF31for example, ibuprofen was found to reduce the
rate of decline in FEV1 values in CF patients under 13
years of age.32 Since ibuprofen seems to lessen the
effects of TNF
,33 it may be important to determine the genotype
of TNF
polymorphic variants of a patient before deciding on
treatment with anti-inflammatory agents. Indeed, administration of
such drugs to CF patients may be of therapeutic value in some
patients but detrimental in others, based on the TNF
genotype. It should be noted that ibuprofen was not
administered to any of the Belgian CF patients in this study but it
was given to six Czech CF patients (+691g ins homozygotes) for
a few weeks only.
In conclusion, we did not find any evidence that the 308g/a
TNF
polymorphic locus is associated with CF lung disease in
Belgian and Czech patients with CF. However, our results indicate
an involvement of TNF
in the modulation of CF disease since the +691g
ins/del and 851c/t loci are associated with lung disease
severity. Functional studies are needed to confirm this association
and to unravel the mechanism of modulation.
|
|
FOOTNOTES |
|---|
|
|
REFERENCES |
|---|