Conceived and designed the experiments: MDW AN. Performed the experiments: SK ALC BN AN. Analyzed the data: SK GCE ALC AN. Contributed reagents/materials/analysis tools: MDW DJP AN. Assisted in the drafting and editing of the manuscript and figures: SK MDW GCE ALC BN DJP WNR AN.
After 9/11/2001, some Fire Department of New York (FDNY) workers had excessive lung function decline. We hypothesized that early serum matrix metalloproteinases (MMP) expression predicts World Trade Center-Lung Injury (WTC-LI) years later.
This is a nested case-control analysis of never-smoking male firefighters with normal pre-exposure Forced Expiratory Volume in one second (FEV1) who had serum drawn up to 155 days post 9/11/2001. Serum MMP-1, 2,3,7,8, 9, 12 and 13 were measured. Cases of WTC-LI (N = 70) were defined as having an FEV1 one standard deviation below the mean (FEV1≤77%) at subspecialty pulmonary evaluation (SPE) which was performed 32 months (IQR 21–53) post-9/11. Controls (N = 123) were randomly selected. We modeled MMP's ability as a predictor of cases status with logistic regression adjusted for time to blood draw, exposure intensity, weight gain and pre-9/11 FEV1.
Each log-increase in MMP-3 and MMP-12 showed reduced odds of developing WTC-LI by 73% and 54% respectively. MMP-3 and MMP-12 consistently clustered together in cases, controls, and the cohort. Increasing time to blood draw significantly and independently increased the risk of WTC-LI.
Elevated serum levels of MMP-3 and MMP-12 reduce the risk of developing WTC-LI. At any level of MMP-3 or 12, increased time to blood draw is associated with a diminished protective effect.
During the events of September 11th, 2001, the collapse of the World Trade Center (WTC) pulverized much of the building's glass and metal structure and released an estimated 10 million tons of particulate matter (WTC-PM).
Approximately 92% of the exposed rescue workers were subsequently enrolled in the Fire Department of New York-World Trade Center-Medical Monitoring and Treatment Program.
We have previously investigated serum inflammatory, metabolic syndrome and cardiovascular biomarkers and found them to be predictors of FEV1 decline and WTC-LI.
The role of proteases has been studied in the setting of many diseases including cancer, coronary disease, chronic obstructive pulmonary disease (COPD), and cigarette-induced chronic lung diseases.
This study investigates the levels of MMPs present in serum after exposure to WTC dust in rescue workers as potential systemic biomarkers predicting susceptibility to WTC-LI. It tests the hypothesis that proteases present within three months of 9/11/2001 can predict airway injury years later.
This study was approved by both Montefiore Medical Center and New York University IRB. All subjects signed informed Institutional Review Board-approved consent at the time of enrollment allowing analysis of their information and samples for research (Montefiore Medical Center; #07-09-320 and New York University; #11-00439). De-identified data devoid of sensitive information acquired through WTC National Institute of Occupational Safety and Health (NIOSH) funding is available to health researchers and others in accordance with the Zadroga Act and in full compliance with the Centers for Disease Control and Prevention (CDC) and the Agency for Toxic Substances and Disease Registry (ATSDR) Policy on Releasing and Sharing Data by request to the WTC Health Program directors.
Prior to 9/11, all subjects had spirometry as part of an annual physical. At medical monitoring entry (MME), serum samples and the first post-9/11 spirometry were obtained (Portascreen Spirometry; S&M Instruments). The population was evaluated and monitored, and N = 13,234 symptomatic firefighters were referred to subspecialty pulmonary evaluation (SPE) if they became symptomatic.
Similar to previous studies a nested case-control design was used to measure the association of serum analytes with case status.
Subjects were included in the study if they were never-smoking male firefighters who had reliable National Health and Nutrition Examination Survey (NHANES) normative data for predicted FEV1, post-9/11 FDNY PFTs within 200 days of 9/11, and pre-9/11 FEV1 >75% pred (n = 801 (47%) out of 1720).
Blood drawn at MME was allowed to stand for 1 hour at room temperature, and then centrifuged at 1,800g for 10 minutes. Serum was stored at −80°C (Bio-Reference Laboratories, Inc. Elmwood Park, NJ), thawed once at 4°C, and assayed using MMP panel (Procarta/Affymetrix), which had a detection range of 2.44–40,000 ng/mL. Panels were read on Luminex 200-IS (Luminex Corporation, Austin, TX) and analyzed with MasterPlexQT software (Ver. 1.2; MiraiBio, Inc., San Diego, CA). MMP values were internally validated using manufacturer controls and standards of known concentrations. Each plate contained a 1∶2 ratio of cases to controls to control for batch artifact.
We tested normality using the Shapiro-Wilk test and Q-Q plots. We used Mann–Whitney U test for between group comparisons, as appropriate. Pearson's Chi-squared test was used for inferences on proportions.
Odds ratios were calculated by multivariable logistic regression using case definition as the outcome variable. Since biomarker levels were skewed, values were log-transformed and included as continuous covariates in the model. The Hosmer-Lemeshow goodness-of-fit test was used to assess the calibration of the model. The model discrimination was quantified using the receiver operating characteristic (ROC) area under the curve (AUC).
Hierarchical clustering was performed using Cluster 3.0 (Ver. 1.47, Michael Eisen: Stanford University; Michiel de Hoon: University of Tokyo) and Java Treeview (Ver.1.1, Saldanha).
Data are expressed as median (interquartile range, IQR) or Odds Ratio (95% confidence interval), unless otherwise stated. A two sided p value less than 0.05 was considered significant. Database management and analyses were performed using SPSS 20 (IBM, Armonk, NY) and STATA/SE 12 (StataCorp, College Station, TX).
Derivation of the baseline cohort, cases and controls are described in
Derivation of Study Cohort from N = 1720 symptomatic FDNY firefighters who presented for subspecialty pulmonary testing. Serum available and inclusion criteria met for N = 70/100 cases and N = 123/153 controls.
| Date/Event | Baseline Cohort | Susceptible Cases | Sub-Cohort Controls | p | |
| N | 801 | 70 | 123 | ||
| High | 197(25%) | 18(26%) | 21(17%) | 0.151 | |
| Intermediate | 604(75.4%) | 52(74%) | 102(83%) | ||
| MME | 2.7(2–4) | 2.7(2–4) | 2.5(2–3) | 0.145 | |
| SPE | 33.8(25–57) | 32.6(21–53) | 35.5(26–55) | 0.327 | |
| MME | 28.0(26–30) | 29.0(27–31) | 28.0(26–31) | 0.106 | |
| SPE | 28.9(27–31) | 29.6(27–34) | 29.0(27–31) | 0.018 | |
| Change | 0.8(–0.1–1.8) | 1.1(0–2.3) | 0.6(–0.4–1.7) | 0.003 | |
| 9/11/01 | 13(7–19) | 15(8–18) | 14(7–18) | 0.907 | |
| 9/11/01 | 40(36–45) | 40(36–45) | 42(37–46) | 0.764 |
Median (IQR).
WTC, World Trade Center; MME, Medical Monitoring Exam; SPE, Subspecialty Pulmonary Exam; PFT, Pulmonary Function Test; BMI, Body Mass Index.
Three sequential pulmonary functions tests (PFT) were performed on this cohort; prior to exposure (Pre-9/11), at MME, and at SPE [
Median and IQR are represented by box plots, with median drawn in the middle of the box. The extremes of the error bars (whiskers) represent 10–90% percentile.
| Time | Variable | Cases | N | Controls | N | p |
| 88(81–96) | 70 | 104(92–113) | 123 | <0.001 | ||
| 81.7(78–86) | 70 | 84.9(81–88) | 123 | 0.001 | ||
| 78(71–89) | 70 | 93(84–99) | 123 | <0.001 | ||
| 81.4(76–86) | 70 | 83.8(80–87) | 123 | 0.016 | ||
| 72(66–75) | 70 | 96(88–104) | 123 | <0.001 | ||
| 71.1(65–77) | 70 | 77.1(73–81) | 123 | <0.001 | ||
| 15(7–29) | 46 | 5(2–8) | 46 | <0.001 | ||
| 30(65) | 46 | 10(22) | 46 | <0.001 | ||
| 0.24(0.06–1.78) | 34 | 0.05(0.03–0.11) | 102 | 0.001 | ||
| 17(50) | 34 | 19(19) | 102 | <0.001 | ||
| 96(83–106) | 47 | 103(98–109) | 53 | 0.002 | ||
| 130(109–157) | 47 | 123(111–140) | 53 | 0.525 | ||
| 96(85–107) | 46 | 107(101–116) | 52 | <0.001 | ||
| 83(75–89) | 35 | 94(87–101) | 43 | <0.001 | ||
| 122(113–134) | 36 | 117(105–122) | 40 | 0.038 | ||
| 18(45) | 40 | 25(39) | 64 | 0.550 | ||
| 22(55) | 40 | 27(42) | 64 | 0.203 | ||
| 13(33) | 40 | 23(36) | 64 | 0.833 |
Median (IQR).
FEV1, Forced Expiratory Volume in one second; FVC, Forced Vital capacity; BDR, Bronchodilator Response; MCT, Methacholine Challenge Testing; PC20, Provocative concentration of methacholine that results in a 20% drop in FEV1; TLC, Total Lung Capacity; RV, Residual Volume; DLco, Diffusing Capacity of the Lung for Carbon Monoxide; VA, Alveolar Ventilation; BWT, Bronchial Wall Thickening.
Although cases had lower pre-9/11 FEV1 (88%) than the controls (104%), all FDNY workers had normal FEV1 by design. The controls and cases suffered similar decline of FEV1 in the interval from pre-9/11 to MME (11% for the controls and cases). Importantly, controls regained some FEV1 between MME and SPE, whereas cases continued to lose lung function over this interval.
Cases had more decline in FEV1/FVC ratio between MME and SPE (0.81 to 0.71) than controls (0.84 to 0.77) [
MMP-1,3,8 and 12 were significantly lower in cases that would progress to WTC-LI than controls who would subsequently recover some of the lung function lost after WTC exposure. There was no difference in MMP-2,7,9 and 13 expression in cases and control [
| Analyte | Cases | Controls | p |
| pg/mL | N = 70 | N = 123 | |
| 387(116–864) | 775(296–1368) | 0.001 | |
| 2840(1281–5130) | 3020(1815–4640) | 0.268 | |
| 3194(1962–7542) | 7653(3320–13765) | <0.001 | |
| 293(67–396) | 222(96–320) | 0.382 | |
| 2(2–20) | 2(2–174) | 0.030 | |
| 25610(12222–74000) | 23490(11196–47084) | 0.548 | |
| 35(7–218) | 66(22–313) | 0.008 | |
| 58(9–109) | 75(3–141) | 0.355 |
p<0.05; All values shown as median (IQR).
Hierarchical clustering of serum biomarkers tested coordinate expression of the MMPs [
Clustering of MMPs in the cohort, N = 193 showed that MMP3 and MMP-12 clustered together.
We then tested if the MMPs that were different between cases and controls predicted future lung function with logistic regression. All MMP levels remained as continuous variables and were log-transformed to approximate a normal distribution [
| Model | OR (95% CI) | AUC | HL |
| 0.476 (0.261–0.866) | 0.819 (0.755–0.882) | 0.034 | |
| 1.013 (1.000–1.027) | |||
| 0.267 (0.121–0.589) | 0.832 (0.771–0.893) | 0.182 | |
| 1.014 (1.000–1.027) | |||
| 0.675 (0.410–1.111) | 0.810 (0.745–0.874) | 0.075 | |
| 1.010 (0.997–1.024) | |||
| 0.462 (0.260–0.821) | 0.818 (0.755–0.881) | 0.506 | |
| 1.013(1.000–1.027) |
Each Model Includes: ΔBMI between MME and SPE, Exposure Group, Pre-9/11 FEV1% Predicted, Time to MME (days) and MMP (Log10 pg/mL).
OR, Odds Ratio; CI, Confidence Interval; AUC, Area Under the Curve; HL, Hosmer Lemeshow.
Four models were developed using case definition as the outcome variable, and assessed for quality using ROC analysis and Hosmer-Lemeshow goodness of fit statistic. MMP-8 was not a significant predictor of lung function after adjustments. Although MMP-1, MMP-3, and MMP-12 demonstrated significant predictive ability, MMP-1 was discounted from further analysis because it failed Hosmer-Lemeshow goodness of fit test. Time to serum draw was also a significant covariate for MMP-3 and MMP-12; every day post-exposure increased the odds of developing WTC-LI by 1.4% and 1.3% respectively. The odds ratios for MMP-3 and MMP-12 were 0.267 and 0.462 respectively, indicating a protective function of elevated levels of the proteases. MMP-3 and MMP-12 demonstrated excellent predictive ability with AUC of 0.832 and 0.818 respectively. The MMP-12 and MMP-3 models have a sensitivity of 58.6% and specificity of 86.2%.
To better understand the effect of time on the association between serum MMP level and WTC-LI, we utilized a contour plot [
Contour Plots express probability isopleths for the development of WTC-LI with all other covariates held constant. When either MMP-3 (A) or MMP-12 (B) increases, the probability of lung injury decreases. As time to blood draw increases, the probability of lung Injury increases.
We report that early elevated expression of MMP-3 and MMP-12 in serum within 200 days after WTC exposure predicts protected lung function over the subsequent seven years (between 2001–2008). We found that MMP-3 and MMP-12 levels of the cohort were comparable to other published patient populations, including healthy controls and patient cohorts with emphysema or rheumatoid arthritis.
This cohort provides a unique opportunity identify biomarkers of susceptibility to dust induced lung injury. First, FEV1 was measured prior to exposure. Second, after the WTC collapse produced a massive acute dust exposure, serum was drawn early during the evolution process to disease. Longitudinal FEV1 measurements allowed us to define case status years after the serum was drawn. The cohort is also large enough to study only never-smokers, eliminating a major confounder of reduced FEV1. The development of FEV1 loss in obstructive lung disease is a multifactorial process. MMPs are central in the disease pathogenesis of obstructive lung disease, and are thus of interest in this population. MMPs are involved in parenchymal destruction and lung remodeling, and their quantification allows us to start to understand their role in the development of FEV1 loss in this WTC exposed population. Furthermore, they are involved in separate pathways of disease attenuation than those of inflammation, cardiovascular disease, and metabolic syndrome that were the focus of our prior work.
Cases and controls differed in several ways that impacted future analyses. Incident cases gained significantly more weight than controls by the time of SPE. BMI is a known cause of decreased lung function, and is thus a possible confounder. Because our case definition of WTC-LI relies on FEV1 at SPE, it is heavily impacted by pre-9/11 FEV1 and BMI at SPE, which are both significantly different between cases and controls prior to model building. For this reason, we adjust for pre-9/11 FEV1 and change in BMI (as a measure of weight gain) as confounders in the logistic regression model. The change in BMI in the logistic regression is simply used as a marker of weight gain. The model remains robust whether the weight gain is measured as change in BMI or BMI at SPE.
Cases and controls suffered similar acute reductions in FEV1 from pre-9/11 to MME. This is consistent with the similar WTC exposure intensities in both case definitions. However, controls had higher pre-9/11 FEV1 than susceptible cases (102 vs 88%). The differences in pre-9/11 FEV1 could reflect predisposition to lung injury after non-tobacco smoke exposure during the population's 18 years of service in fighting fires. This is suggested by continued deterioration of cases in FEV1 compared to the relative recovery of controls years after the acute WTC exposure. Because cases were defined by a threshold of FEV1 of <77% at SPE, those with pre-9/11 FEV1 close to threshold values were more likely to meet case definition. Cases also had markedly diminished airway function across multiple measures after exposure, including FEV1, lower FEV1/FVC ratio, more bronchodilator response, and steeper methacholine response slope.
We observe that increased MMP-3 and 12 expression reduces the odds of subsequent lung injury. Also, the time to medical monitoring enrollment and serum draw post-exposure was an independent predictor of lung injury. Those who had serum obtained two months post-9/11 (early presenters) had lower risk of lung injury than those who had serum obtained four months post-9/11 (late presenters). This result could be due to selection bias with individuals entering medical monitoring later at a more advanced stage of disease. Alternately, delay in early interventional treatment for those who enrolled in medical monitoring later may have inhibited their ability to recover from the exposure. We investigated this possibility, but did not find differences between early and late presenters in known confounders such as exposure intensity or pre-9/11 FEV1. There was a significant increase in BMI in the late presenters, but the effect of time did not change with BMI as a covariate in the logistic models. Time to serum draw was not a significant covariate with all other serum biomarkers of risk or protection that we have reported.
The protective role of MMP-12, also known as macrophage elastase, in this cohort was surprising since most studies on patients with established lung disease demonstrate elevated MMP-12 as a biomarker of disease severity in COPD patients.
This is the first paper that describes a protective function of elevated MMP-3, stromelysin-1, in humans. This novel correlation needs to be shown in other lung injury cohorts prior to investigating the relationship of progressive injury pathways and MMP-3 expression. MMP-3 has been shown to be protective in murine models of atherosclerosis.
This nested case control study has several limitations. The cohort was narrowly defined, reducing generalizability of these results. Specifically, the vast majority of FDNY workers were exposed at the time of the WTC destruction. The few unexposed workers were markedly less healthy than the exposed workers pre-9/11 and most were not on active duty at the time of the event. Therefore, all the study participants were highly exposed to WTC dust and products of combustion. Since there is no unexposed control group with comparable lung function, years of service, age, racial background and serum available it is therefore not feasible to construct an unexposed control. Also, using FEV1 to define case status may produce misclassification of disease. However, the authors maintain that using FEV1 is the simplest, most robust, and best validated single measure of airway injury in an obstructed cohort. Since the design reduced the potential for selection bias, misclassification of disease is non-differential. The reported associations therefore likely underestimate the true strength of association between the serum biomarkers and subsequent FEV1. Although the serum biomarkers are expressed years before the disease is defined, the results are correlations and do not imply causation. This restricts our ability to assess the impact of WTC exposure to the observed biomarker disease relationship. Biomarker associations with lung injury must also be investigated in independent cohorts.
WTC exposure impacted multiple distinct injury and repair pathways. One interpretation of the findings is that biomarkers that reduce the risk of WTC-LI reflect effective repair after an acute dust induced lung injury. We report that elevated MMP-3 and MMP-12 are associated with a reduction in the odds of developing WTC-LI years later. Finally, we show that time to serum sample acquisition is independently associated with the probability of lung injury and highlights the utility of serum sampling efforts continuing for some time. We will investigate in future studies about the relationship between genetic polymorphisms that may improve or change the predictive ability of MMP-3 and MMP-12. Also, we will analyze longitudinal samples of serum on the same cohort to understand how chronicity of the biomarkers impact disease development.
The authors would like to thank the firefighters and rescue workers for their participation in this study and for their selfless contributions.