|
|
||||||||
1 Division of Pulmonary and Critical Care Medicine, Dept of Medicine, Mount Sinai School of Medicine, New York, 2 Environmental Health Dept (Occupational Health Program), and 3 Dept of Biostatistics, Harvard School of Public Health, and 4 Pulmonary and Critical Care Unit, Dept of Medicine, Massachusetts General Hospital, Harvard Medical School, Boston, MA, USA.
CORRESPONDENCE: D. C. Christiani, Harvard School of Public Health, 665 Huntington Avenue, Boston MA, 02115, USA. Fax: 1 6174323441. E-mail: dchristi@hsph.harvard.edu
Keywords: Acute respiratory distress syndrome, genetic susceptibility, tumour necrosis factor
Received: January 3, 2005
Accepted May 17, 2005
| ABSTRACT |
|---|
|
|
|---|
The above hypothesis was investigated in a nested case-control study of 441 Caucasian controls and 212 cases admitted to an intensive care unit with sepsis, trauma, aspiration or hyper-transfusions.
The 308A and TNFB1 alleles were in linkage disequilibrium. These polymorphisms were not associated with ARDS susceptibility on crude analysis. On subgroup analyses, they were associated with either increased or decreased odds of developing ARDS depending on whether the clinical risk for ARDS results in direct or indirect pulmonary injury. The 308A allele was associated with increased 60-day mortality in ARDS, with the strongest association found among younger patients. There was no association between the TNFB polymorphism and ARDS mortality.
The 308GA, but not the TNFB12, polymorphism was associated with increased mortality in acute respiratory distress syndrome, but their association with acute respiratory distress syndrome susceptibility depended on the site of injury predisposing to acute respiratory distress syndrome.
Understanding why some individuals develop and subsequently die from acute respiratory distress syndrome (ARDS), while others do not, is incomplete. Although clinical predictors such as sepsis and trauma are well recognised, only a minority of patients with these risks develop ARDS 1. It is likely that after the same type and degree of insult, individual differences in susceptibility to developing and dying from ARDS exist. Genetic susceptibility to acute lung injury may explain the observed inter-individual differences in risk and outcomes 24.
Studies on the role of tumour necrosis factor (TNF)-
in ARDS have been conflicting. Some found plasma TNF-
to correlate with the development of, and mortality in, ARDS 57, while others did not 8, 9. These discrepancies may be due to temporal and regional variation in TNF-
release 10, 11 and to differences in assay techniques 12. Clinical heterogeneity and genetic variability in the production of TNF-
may also contribute to such discrepancies 13.
The 308GA promoter polymorphism in the TNF-
gene and the TNFB1/2 NcoI restriction fragment length polymorphism (RFLP) in the TNF-ß gene appear to influence TNF-
level. Carriers of the 308A allele (308A) and homozygotes for the TNFB2 allele (TNFB22) have increased TNF-
production 14, 15 and increased susceptibility to, or increased mortality in, septic shock in some studies 1519, but not in others 20, 21.
In recognition of the clinical heterogeneity in ARDS and the differences between pulmonary and extra-pulmonary ARDS 22, the American European Consensus Conference (AECC) on ARDS suggested dividing aetiologies of lung injury into direct or indirect pulmonary injuries 23. Recently, the 1580C/T missense mutation in the surfactant protein-B gene was linked to ARDS in patients with direct, but not indirect pulmonary injury 2. Thus, it is possible that any associations of the 308GA and TNFB1/2 polymorphisms on the development of ARDS may vary with direct versus indirect pulmonary injury.
The present authors describe a nested case-control study of patients at risk for ARDS and hypothesised that the 308A of the 308GA polymorphism and the TNFB22 of the TNF-ß gene would be associated with increased susceptibility to, and increased mortality in, ARDS. In addition, association was explored to determine if it might differ by direct versus indirect pulmonary injury.
| METHODS |
|---|
|
|
|---|
|
|
Study design
Baseline clinical information was collected on admission to an ICU. Vital signs and laboratory parameters in the first 24 h after ICU admission were collected for calculation of Acute, Physiology, Age and Chronic Health Evaluation (APACHE) III 25. Subjects were screened daily for ARDS, defined by respiratory failure requiring intubation and AECC criteria as previously described 3, 23. Daily chest radiographs were interpreted by two pulmonary and critical care physicians after consensus training on interpretation of infiltrates with disagreements arbitrated by a third physician. All were blinded to the clinical status of the patients. The
-score for agreement between initial interpretation for bilateral infiltrates was 0.75 (95% confidence interval (CI): 0.620.89), comparable to other reports 26.
From the prospective cohort, a case-control study was designed. Cases were those who fulfilled criteria for ARDS during their hospitalisation. All patients who did not develop ARDS during their hospitalisation and had no prior history of ARDS or prior enrolment into the study were selected as controls. All patients were followed for ICU mortality and ARDS cases were followed for all causes of mortality within a 60-day period.
Methods
Blood (10cc) was collected for DNA extraction and PCR amplification. Genotyping was performed with pyrosequencing for 308GA polymorphism (Pyrosequencing AB, Uppsala, Sweden), with the forward primer of 5'-CCAAACACAGGCCTCAGGACTC-3', biotinylated reverse primer of 5'-TCCTCCCTGCTCCGATTCCG-3' and sequencing primer of 5'-AGGCAATAGGTTTTGAGGGGCA-3'. The TNF-ß Nco I polymorphism was detected using a modified PCR-RFLP method that involved published primer sequences and Nco I enzyme digestion (New England BioLabs, Beverly, MA, USA) 19. Results were interpreted by two independent investigators, and genotyping was repeated in a random 5% of samples. Laboratory personnel and research assistants were blinded to the case-control status or genotype of the subjects.
Analysis
Conformity to Hardy Weinberg Equilibrium was determined with a Chi-squared goodness of fit test. Univariate analysis was performed using Fisher's Exact Test, Chi-squared trend test, ANOVA or Wilcoxon Rank Sum tests as appropriate. Variables with p-value
0.2 on univariate analyses were studied in a backwards selection algorithm and eliminated if they did not meet a p-value of
0.2. The final multivariate logistic regression model included the gene effect, results from backwards elimination, significant interactions and clinically relevant parameters such as APACHE III scores for the development of ARDS and septic shock for mortality in ARDS. Possible confounders were adjusted for in the final model, rather than matched between cases and controls on analysis. Given the small number of 308A homozygotes, homozygotes for the 308A allele were grouped with the heterozygotes (308A), and compared with 308G homozygotes in the final model, as in prior reports 15, 16. As the distribution of the TNFB genotype among cases and controls did not clearly suggest a dominant or recessive model, the TNFB genotype was modelled as in previous reports 17, 18 by comparing TNFB2 homozygotes (TNFB22) with homozygotes and heterozygotes for TNFB1. C-statistics (area under the receiver operating characteristic curve) were used to evaluate model fit 27. An a priori decision was made to stratify the analysis by direct versus indirect pulmonary injury. Effect modification was tested with an interaction term. A p-value of 0.05 was considered statistically significant.
Detection of linkage disequlibrium between the two polymorphisms was based on Lewontin's
in controls 28. Haplotypes of the two TNF polymorphisms were generated using the Partition Ligation-Expectation Maximization (PL-EM) version 1.0 29, as in other association studies 30. This software uses an efficient variant of the EM algorithm, to reconstruct individual probabilities for individual phasing accuracy based on unphased genotype data, while providing estimates on the overall haplotype frequencies, as well as their standard errors. Odd ratios (OR) were estimated by comparing the individuals with one or more copies of the haplotype of interest with the individuals without the haplotype.
In the initial study design, assuming an
-error of 0.01 (to allow for multiple comparisons), 80% power and allele frequencies of 16% for 308A 31 and 65% for TNFB2 17, a study with 560 cases and 1,120 controls would have a minimum detectable OR for the development of ARDS of 1.56 for 308A and 1.47 for TNFB22. Assuming 50% of patients had direct pulmonary injury as a risk for ARDS, the minimum detectable OR was 1.9 for 308A and 1.8 for TNFB22. Given the actual study size of 653 subjects and
-error of 0.05, the minimum detectable OR would be 1.78 for 308A, and 1.68 for TNFB22.
| RESULTS |
|---|
|
|
|---|
Clinical risk factors for ARDS and baseline characteristics on admission to the ICU are shown in tables 2
and 3
. ARDS developed a median of 1 day after ICU admission (2575%; 03 days). Results from backwards selection for development of ARDS were; direct pulmonary injury (p<0.001), sepsis (p = 0.05), trauma (p<0.001), age (p = 0.002), female (p = 0.02), diabetes (p = 0.004), platelets
80,000·mm1 (p = 0.005), abnormal bilirubin >2.0 mg·dL1 (p = 0.1) and blood transfusion (p<0.001). A statistically significant interaction was only found between the type of pulmonary injury (direct versus indirect) and the 308A allele or the TNFB22 genotype.
|
|
= 0.90). Applying PL-EM algorithm, the posterior probabilities of individual haplotypes ranged from 0.961.0. Therefore, two haplotypes with the highest posterior probability were assigned to each individual.
The genotype frequencies among the cases and controls are detailed in table 4
. Genotype and haplotype frequencies did not differ between cases and controls. The C-statistic for the multivariate model for ARDS was 0.73 and 0.74 for 308A and TNFB22, respectively. On multivariate analyses, an association with ARDS was found for TNFB22 (adjusted OR: 0.47; 95% CI: 0.250.86; p = 0.01), but not for 308A (adjusted OR: 1.7; 95% CI: 0.893.1).
|
|
2.0 mg·dL1 (p = 0.05), and blood transfusion (p = 0.03). The only significant interaction was between transfusion and 308A (p = 0.05) and TNFB22 (p = 0.02). The C-statistics for the final model for ARDS mortality was 0.85 for both 308A and TNFB22.
ARDS mortality differed significantly depending on the 308GA, but not the TNFB1/2 polymorphism (table 4
) with the numbers of 308A alleles associated with increasing 60-day ARDS mortality (p = 0.01; fig. 3
). The 308A allele was associated with a significantly increased 60-day mortality in ARDS (crude OR: 2.1; 95% CI: 1.13.9; adjusted OR: 3.5; 95% CI: 1.48.6; p = 0.007). This association was strongest among the 117 patients <67 yrs (adjusted OR: 14.9; 95% CI: 3.074; p<0.001). No association with mortality was found among the 95 ARDS patients >67 yrs (p = 0.3). The 308A allele was associated with mortality in both direct (adjusted OR: 4.8; 95% CI: 1.318) and indirect pulmonary injury (adjusted OR: 5.5; 95% CI: 0.9931) with no evidence of effect modification by type of injury (p = 0.2). TNFB22 was not associated with ARDS mortality (adjusted OR: 0.50; 95% CI: 0.211.2). As expected from the linkage disequilibrium, the 308A:TNFB1 haplotype was associated with increased ARDS mortality (adjusted OR: 3.7; 95% CI: 1.59.1; table 3
).
|
| DISCUSSION |
|---|
|
|
|---|
The current study has a number of strengths. First, the prospective determination of ARDS using the AECC definition helps minimise phenotype misclassification given there is no diagnostic gold standard for ARDS. Secondly, clearly defined at-risk controls were used in this study. Using critically ill controls that have the opportunity to develop the outcome is more clinically relevant than using healthy individuals. In addition, this reduces the confounding from any possible association between the gene and the risk condition such as sepsis or pneumonia.
In support of the current authors hypothesis, the reportedly high TNF producing 308A allele of the 308GA polymorphism was associated with increased mortality in ARDS, consistent with other studies in malaria and sepsis 15, 16, 30 This study suggests that genetic heterogeneity among patients may partially explain the differences in prior studies in circulating TNF-
in ARDS and sepsis 510.
The association between 308A and the 308A:TNFB1 haplotype and mortality in ARDS was particularly strong among younger patients even though younger ARDS patients tended to have lower mortality (table 3
). Genetic contribution to complex diseases is greater in diseases with early age of onset 32.
While an association was found between the TNFB22 and ARDS susceptibility on multivariate analysis, this was found to be driven by its differing association with ARDS among subgroups of patients with direct and indirect pulmonary injury. A similar heterogeneous effect was seen with the 308A. Such effect modification by the type of injury on the association between the TNF polymorphisms and ARDS suggests a possible gene-environment interaction in ARDS. The significance of this finding is not clear. The inflammatory response and the radiological, histological and mechanical properties of the lung differ depending on whether the site of infection or the aetiology of ARDS is pulmonary or extra-pulmonary 22, 33, 34. Alternatively, it is possible that these findings are due to the lower sensitivity and specificity of the AECC criteria for ARDS in patients with direct versus indirect pulmonary injury. A recent study found that bilateral pneumonia can often be misdiagnosed as ARDS by AECC criteria 35. While bilateral pneumonia, bilateral aspiration and pulmonary contusion are indistinguishable from ARDS clinically, they have different histology and pathogenesis. Therefore, it is difficult to know whether the high TNF-
producing 308A allele is associated with decreased progression to bilateral pneumonia or a decreased odd of developing the diffuse alveolar damage of ARDS. It is possible that an adequate or high TNF-
response may limit the direct pulmonary injury of pneumonia or aspiration of oral pathogens and decrease the odds of developing bilateral pneumonia. However, even with such misclassification, most patients with direct pulmonary injury who do fulfil AECC criteria do have ARDS histologically, indicating that the ARDS cases may be more homogeneous 24. However, given the reduced sample size in the subgroups, these findings will need to be confirmed in a larger population of patients with clinically homogenous risks for ARDS.
In the present study, the haplotype analysis closely parallels the results of the 308A allele. This suggests that the 308GA polymorphism is more important than the TNFB polymorphism in ARDS mortality, with little contribution from the TNFB polymorphism to the overall findings. However, since the 308GA and TNFB1/2 polymorphisms lie in the highly polymorphic, major histocompatibility region, the possibility that the 308A polymorphism is linked to another polymorphism that is the disease locus cannot be excluded. The current study focused on two particular polymorphisms, specifically because of their previous association with variable function and sepsis, but there are certainly many more polymorphisms in the TNF gene. A formal haplotype analysis after genotyping multiple polymorphisms would be preferable for defining the linkage disequilibrium in the genes and their relationship to ARDS.
The authors acknowledge some other limitations to the present study. The functional significance of these TNF polymorphisms is not confirmed in the patient population. Due to the study design, the results may not be generalised to the community setting, to immunocompromised hosts, to patients without risk factors for ARDS or to patients with different clinical risks for ARDS. In addition, the analyses were restricted to Caucasians, which reduces the possibility of confounding from different genetic make-up, but it does not permit extrapolation of the results to other ethnic groups.
In conclusion, the present study demonstrated an association between the 308A allele and mortality in acute respiratory distress syndrome, especially among younger patients. The association between the 308GA polymorphism and the development of acute respiratory distress syndrome is heterogeneous and may depend upon the type of injury predisposing to acute respiratory distress syndrome. Additional studies are needed to confirm these findings in other populations with other risk factors.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
| REFERENCES |
|---|
|
|
|---|
is associated with 308 TNF-
promoter polymorphism and with clinical severity in chronic beryllium disease. Am J Respir Crit Care Med 2001;164:11921199.
promoter polymorphism, with septic shock susceptibility and mortality. JAMA 1999;282:561568.
concentrations and outcomes of patients with severe sepsis. Crit Care Med 1996;381384.
This article has been cited by other articles:
![]() |
Z. Wang, D. Beach, L. Su, R. Zhai, and D. C. Christiani A Genome-Wide Expression Analysis in Blood Identifies Pre-Elafin as a Biomarker in ARDS Am. J. Respir. Cell Mol. Biol., June 1, 2008; 38(6): 724 - 732. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Agarwal, R. Srinivas, A. Nath, and S. K. Jindal Is the Mortality Higher in the Pulmonary vs the Extrapulmonary ARDS?: A Metaanalysis Chest, June 1, 2008; 133(6): 1463 - 1473. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Kamp, X. Sun, and J. G. N. Garcia Making Genomics Functional: Deciphering the Genetics of Acute Lung Injury Proceedings of the ATS, April 15, 2008; 5(3): 348 - 353. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. L. Lagan, D. D. Melley, T. W. Evans, and G. J. Quinlan Pathogenesis of the systemic inflammatory syndrome and acute lung injury: role of iron mobilization and decompartmentalization Am J Physiol Lung Cell Mol Physiol, February 1, 2008; 294(2): L161 - L174. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. V. Giannoudis, M. van Griensven, E. Tsiridis, and H. C. Pape The genetic predisposition to adverse outcome after trauma J Bone Joint Surg Br, October 1, 2007; 89-B(10): 1273 - 1279. [Abstract] [Full Text] [PDF] |
||||
![]() |
R Zhai, M N Gong, W Zhou, T B Thompson, P Kraft, L Su, and D C Christiani Genotypes and haplotypes of the VEGF gene are associated with higher mortality and lower VEGF plasma levels in patients with ARDS Thorax, August 1, 2007; 62(8): 718 - 722. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |