Clinical risk factors and blood protein biomarkers of 10-year pneumonia risk
============================================================================

* Ming-Ming Lee
* Yi Zuo
* Katrina Steiling
* Joseph P. Mizgerd
* Bindu Kalesan
* Allan J. Walkey

## ABSTRACT

**Background** Chronic inflammation may increase susceptibility to pneumonia.

**Research Question** To explore associations between clinical comorbidities, serum protein immunoassays, and long-term pneumonia risk.

**Methods** Framingham Heart Study Offspring Cohort participants ≥65 years were linked to their Centers for Medicare Services claims data. Clinical data and 88 serum protein immunoassays were evaluated for associations with 10-year incident pneumonia risk using Fine-Gray models for competing risks of death and least absolute shrinkage and selection operators for covariate selection.

**Results** We identified 1,370 participants with immunoassays and linkage to Medicare data. During 10 years of follow up, 428 (31%) participants had a pneumonia diagnosis. Chronic pulmonary disease [subdistribution hazard ratio (SHR) 1.87; 95% confidence interval (CI), 1.33-2.61], current smoking (SHR 1.79, CI 1.31-2.45), heart failure (SHR 1.74, CI 1.10-2.74), atrial fibrillation/flutter (SHR 1.43, CI 1.06-1.93), diabetes (SHR 1.36, CI 1.05-1.75), hospitalization within one year (SHR 1.34, CI 1.09-1.65), and age (SHR 1.06 per year, CI 1.04-1.08) were associated with pneumonia. Three baseline serum protein measurements were associated with pneumonia risk independent of measured clinical factors: growth differentiation factor 15 (SHR 1.32; CI 1.02-1.69), C-reactive protein (SHR 1.16, CI 1.06-1.27) and matrix metallopeptidase 8 (SHR 1.14, CI 1.01-1.30). Addition of C-reactive protein to the clinical model improved prediction (Akaike information criterion 4950 from 4960; C-statistic of 0.64 from 0.62).

**Conclusions** Clinical comorbidities and serum immunoassays were predictive of pneumonia risk. C-reactive protein, a routinely-available measure of inflammation, modestly improved pneumonia risk prediction over clinical factors. Our findings support the hypothesis that prior inflammation may increase the risk of pneumonia.

**Summary or “take home” message** Here, we show that long-term pneumonia risks are associated with clinical characteristics and serum markers of inflammation. Our findings offer preliminary support to the hypothesis that chronic inflammation may be a marker of increased pneumonia risk.

## INTRODUCTION

Pneumonia is the leading infectious cause of death globally [1] and is a substantial public health issue. In the U.S., pneumonia accounts for an estimated 1.5 million adult hospitalizations annually with 100,000 in-hospital deaths [2]. Furthermore, pneumonia survivors experience numerous deleterious consequences, such as reduced pulmonary function, new cardiovascular complications, decline in cognitive function, depression, decreased functional status and quality of life, and shortened lifespans [3].

Prior studies have investigated clinical risk factors for pneumonia and found associations with age (both the youngest and oldest), chronic conditions, socioeconomic status and lifestyle factors such as smoking and alcohol ingestion [4]. Given associations between pneumonia, pre-existing comorbid conditions, and future health decline, some have hypothesized that pneumonia represents an acute manifestation of a chronic susceptibility state [5]. Although the mechanisms of increased pneumonia susceptibility are multi-factorial and incompletely understood, prior studies have posited links between chronic inflammation, immunosenescence and pneumonia risk [6-8]. However, few studies have evaluated associations between comorbidities, inflammation, and pneumonia in humans.

We explored the associations between clinical comorbidities, serum protein markers of organ function and inflammation, and long-term pneumonia risk in a well-characterized, community-dwelling, prospective cohort. Furthermore, we evaluated whether routinely-measured serum markers of inflammation such as C-reactive protein may be associated with future risk of pneumonia. Thus, we hypothesize that clinical risk factors and serum markers of inflammation would be associated with long-term risk for the development of pneumonia.

## METHODS

### Data Source and Cohort

We conducted a retrospective cohort study including participants in the Framingham Heart Study (FHS) Offspring Cohort with FHS data linked to Centers for Medicare and Medicaid Services (CMS) claims data. The design of the FHS Offspring Cohort has been previously described [9], as well as linkage of FHS to CMS data [10]. In brief, FHS Offspring participants underwent standardized examinations and questionnaires with a study physician every 4 to 8 years. We used data from examinations 5 (1991-1995), 6 (1995-1998) and 7 (1998-2001) to characterize patient characteristics potentially associated with pneumonia based upon prior literature reviews [4]. Data missing from more recent examinations were carried forward if present at prior examinations. Data were missing if a participant did not answer a particular questionnaire item or the question was not asked at an examination. In addition to clinical characteristics collected during examinations, FHS data contained measurements of 88 serum protein immunoassays drawn at examination 7 as part of the Systems Approach to Biomarker Research in Cardiovascular Disease Initiative, a sub-study of the FHS (see **Table E1** in the online data supplement) [11]. Samples were obtained between 7 to 9 a.m. after an overnight fast and stored at -80°C after centrifugation and until assay. Luminex™ xMAP multiplex immunoassay technology was used to measure the biomarkers [12]. Fourteen proteins were analyzed as binary variables defined as above or below a detection limit due to a high rate of missing data above the detection limit [11]. CMS inpatient, outpatient, skilled nursing facility, home health agency, hospice, carrier and durable medical equipment files were used to identify pneumonia diagnoses. Participants were included if they were 65 years old or older, had FHS data linked to CMS data, had no history of pneumonia prior to exam 7, and had a protein immunoassay panel measurement at examination 7 (**Figure 1**). Subjects were followed for up to 10 years from the time of examination 7 until diagnosis of pneumonia or death. Data for this study was de-identified, with a waiver of informed consent and protocol approval by the Boston University Medical Center Institutional Review Board (H-35870).

![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/12/09/2023.12.07.23299678/F1.medium.gif)

[Figure 1:](http://medrxiv.org/content/early/2023/12/09/2023.12.07.23299678/F1)

Figure 1: 

### Outcomes

The primary outcome was pneumonia defined by International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM) codes identified from CMS claims within 10 years of follow up. Based upon prior studies [13], we explored several ICD-9-CM coding algorithms for identifying acute pneumonia episodes from claims data. **Table E2** shows that clinical variables associated with pneumonia were similar regardless of the pneumonia definition based on various ICD-9CM codes for pneumonia and in various clinical settings, thus we chose to use the most sensitive definition capturing any pneumonia diagnosis in any clinical setting for all analyses. Death was defined from Medicare denominator files and loss to follow-up from Medicare files was likely rare.

### Covariates

We abstracted baseline demographics, clinical comorbidities, vital signs, anthropometric measurements (i.e., body mass index [BMI]), smoking and alcohol use and measures of functional and cognitive status from FHS clinical exam files. Clinical comorbidities included chronic pulmonary disease, hypertension, atrial fibrillation/flutter, valvular disease, myocardial infarction, heart failure, stroke, diabetes, chronic kidney disease, peripheral vascular disease, pulmonary embolus/deep vein thrombosis, cancer, hip fracture and depression as defined by FHS clinician assessments. Measures of functional and cognitive status included the Kannel physical activity index [14], mini-mental state examination [15], Katz index of independence in activities of daily living [16], physical disability [17] and overnight hospitalization in the year prior to the FHS examination. Detailed definitions of the covariates validated by the FHS protocols are provided in **Table E3**.

### Statistical methods

We used a least absolute shrinkage and selection operator (LASSO) machine approach to select model covariates [18]. LASSO imposes a constraint to the model that minimizes prediction error and causes some regression coefficients to shrink toward 0, excluding those variables, and thus selects variables most strongly associated with the outcome of interest. Because of the presence of missing data deemed missing at random, we used a multiple imputation approach to missing data. We imputed missing data across 5 datasets and used LASSO methods to select covariates averaged across the pooled imputed data sets, where the tuning parameter lamba was selected using 5-fold cross-validation. Because of the potential for competing risk of death for pneumonia outcomes, we used Fine-Gray competing risk models to estimate the subdistribution hazard ratio of each selected predictor of pneumonia [19]. In order to explore contributions from comorbidities and serum protein markers to pneumonia risk, we estimated 4 risk models: 1) a model including only clinical variables from the FHS clinical exam data; 2) a model including only serum protein immunoassay measurements; 3) a model including both clinical variables and serum protein immunoassay measurements and 4) an exploratory model including clinical variables and routinely available protein biomarker C-reactive protein (CRP) (to assess how CRP may add to risk prediction above comorbidities). Models were assessed for fit by measured change in C-statistic (increase) and Akaike Information Criterion (AIC; decrease). Model calibration was assessed by plotting predicted risk of pneumonia against observed frequencies of pneumonia. Because optimal choice of model-building strategy is an area of debate [20] and strongly determines covariates included in the model, we performed sensitivity analyses building clinical models using multiple different strategies (e.g., p<0.2 based selection, use of Cox proportional hazards models that did not account for competing risk of death) to assess the robustness of the multiple imputation LASSO selection approach (**Tables E4 and E5**). We performed additional sensitivity analyses using the four aforementioned risk models restricted to only pneumonia in the inpatient setting to assess for changes in model performance among patients with more severe pneumonia. All statistical analyses were performed with R software version 3.4.3.

## RESULTS

### Cohort Characteristics

We identified 1,370 FHS participants aged 65 years or older with linked CMS claims data and protein immunoassay panel measurements from examination 7. Subject characteristics are summarized in **Table 1**. During the 10 years of follow up, 428 (31%) subjects had a diagnosis of pneumonia. Characteristics of subjects with and without a diagnosis of pneumonia are shown in **Table 2**. Patients with pneumonia more frequently were older, had higher rates of comorbid conditions, physical disability, and smoking.

View this table:
[Table 1.](http://medrxiv.org/content/early/2023/12/09/2023.12.07.23299678/T1)

Table 1. 
Baseline characteristics of participants (N=1,370)

View this table:
[Table 2.](http://medrxiv.org/content/early/2023/12/09/2023.12.07.23299678/T2)

Table 2. 
Distribution of baseline characteristics by pneumonia outcome (N=1,370)

### Missing data

Of the 1,370 participants, the following data regarding baseline characteristics were missing: diabetes (n=120 [9%]), peripheral vascular disease (n=134 [10%]), Kannel physical activity index (n=260 [19%]), valve disease (n=134 [10%]), and BMI (n=138 [10%]). WNK lysine deficient protein kinase 1 was excluded from analyses due to substantial missing data (n=1,254 [92%]). Of the remaining 87 biomarkers, an average of 156 (11%) of samples were missing; fourteen biomarkers (**Table E1**), including interleukin-6 (IL-6), were analyzed as binary variables defined as above or below a detection limit due to missing data regarding levels above the detection limit.

### Factors associated with pneumonia

In the primary analysis using multiple imputation, Fine-Gray competing risk regression models, and LASSO variable selection, clinical variables associated with pneumonia were chronic pulmonary disease (subdistribution hazard ratio [SHR], 1.87; 95% CI, 1.33-2.61), current smoking (SHR, 1.79; 95% CI, 1.31-2.45), heart failure (SHR, 1.74; 95% CI, 1.10-2.74), atrial fibrillation/flutter (SHR, 1.43; 95% CI, 1.06-1.93), diabetes (SHR, 1.36; 95% CI, 1.05-1.75), hospitalization within one year (SHR, 1.34; 95% CI, 1.09-1.65), and age (SHR, 1.06; 95% CI, 1.04-1.08) (**Figure 2a**). In a model including only immunoassay measurements, five of 88 baseline serum protein measurements were predictive of pneumonia, including growth differentiation factor 15 (GDF-15; SHR, 1.65; 95% CI, 1.29-2.12), insulin-like growth factor binding protein 2 (IGFBP-2; SHR, 1.27; 95% CI, 1.00-1.60), matrix metallopeptidase 8 (MMP-8; SHR, 1.20; 95% CI, 1.06-1.37), C-reactive protein (CRP; SHR, 1.15; 95% CI, 1.04-1.27) and N-terminal prohormone of brain natriuretic peptide (NT-proBNP; SHR, 1.15; 95% CI, 1.04-1.28) (**Figure 2b**). In a model adjusted for clinical variables, baseline GDF-15 (SHR, 1.32; 95% CI, 1.02-1.69), MMP-8 (SHR, 1.14; 95% CI, 1.01-1.30) and CRP (SHR, 1.16; 95% CI, 1.06-1.27) remained significant predictors of pneumonia (**Figure 2c**). The addition of CRP (SHR, 1.17; 95% CI, 1.06-1.28) to the clinical model improved the clinical model fit (AIC value 4960 to 4950) and discrimination of pneumonia events (C-statistic 0.62 to 0.64) (**Figure 2d**). Calibration plots showed good calibration across the models, except at extreme levels of pneumonia risk, with the exploratory clinical and CRP model showing optimal calibration (**Figure E1**). Sensitivity analyses of clinical variables in complete case analysis showed similar results to the multiple imputed LASSO approach (**Table E6**). Sensitivity analyses using the risk models restricted to pneumonia in the inpatient setting did not substantively change model performance (C-statistic of the model including only clinical variables was 0.67; C-statistic of the model including only serum protein immunoassays was 0.64; C-statistic of the model including clinical variables and serum protein immunoassays was 0.69; C-statistic of the exploratory model including clinical covariates and C-reactive protein was 0.68 (**Table E7**).

![Figure 2A:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/12/09/2023.12.07.23299678/F2.medium.gif)

[Figure 2A:](http://medrxiv.org/content/early/2023/12/09/2023.12.07.23299678/F2)

Figure 2A: 

![Figure 2B:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/12/09/2023.12.07.23299678/F3.medium.gif)

[Figure 2B:](http://medrxiv.org/content/early/2023/12/09/2023.12.07.23299678/F3)

Figure 2B: 

![Figure 2C:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/12/09/2023.12.07.23299678/F4.medium.gif)

[Figure 2C:](http://medrxiv.org/content/early/2023/12/09/2023.12.07.23299678/F4)

Figure 2C: 

![Figure 2D:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/12/09/2023.12.07.23299678/F5.medium.gif)

[Figure 2D:](http://medrxiv.org/content/early/2023/12/09/2023.12.07.23299678/F5)

Figure 2D: 

## DISCUSSION

We explored clinical and immunoassay predictors of pneumonia risk in a cohort of community dwelling participants in the FHS. In addition to observing multiple clinical variables associated with pneumonia, we identified 3 baseline serum protein immunoassay measures (GDF-15, MMP-8 and CRP), after adjusting for clinical variables, that were associated with pneumonia events up to 10 years prior to the pneumonia diagnosis. Furthermore, the addition of routinely-available CRP to a model including clinical variables improved model performance for discrimination of pneumonia events. Models performed similarly when restricted to more severe pneumonia requiring hospitalization. Our findings suggest that long-term pneumonia risks are associated not only with previously described clinical characteristics, but also with serum markers of inflammation. Taken together, our findings provide preliminary support to the hypothesis that prior inflammation may be a marker of increased pneumonia risk and inform future efforts at pneumonia risk prediction.

We caution against interpreting the association between the clinical variables and protein immunoassays and future pneumonia as casual and acknowledge that these predictors may correlate with unmeasured causal variables. In addition, single measurements of protein immunoassays may not represent a chronic inflammatory state. However, our results expand upon prior studies of inflammation, clinical comorbidities, and pneumonia risk. Yende et al. [8] identified associations between increased baseline levels of inflammatory cytokines (IL-6 and tumor necrosis factor-alpha [TNF-α]) and future pneumonia risk in a cohort of well-functioning, community-dwelling elderly subjects. Unfortunately, a high rate of missing data in our study restricted analysis of IL-6, and tumor necrosis factor-alpha was not included in the immunoassay panel. However, we found other markers of inflammation independent of clinical comorbidities associated with long-term pneumonia risk, including routinely-available clinical measures such as CRP, an acute-phase protein synthesized by the liver. Although the utility of CRP in the diagnosis of acute pneumonia and as a predictor of more severe disease with poor outcomes has been reported [21], the association of baseline CRP with future pneumonia risk has not been described.

Inflammation may increase pneumonia susceptibility through various biological mechanisms. The highest incidence of pneumonia in adults is amongst the elderly, and aging is associated with inflammation and immunosenescence [22,23]. Aged lungs experience greater pathogen adhesion and host cell invasion, a delayed and diminished immune response during early infection, attenuated phagocytic capabilities of macrophages and neutrophils via Toll-like receptor dysfunction, and a prolonged immune response to infection [6,7,24,25]. These pre-clinical findings suggest age-related chronic inflammation results in immune dysregulation that may contribute to increased pneumonia susceptibility. But, the exact causal pathways between the clinical factors and biomarkers and future pneumonia risk remain unclear.

In addition to CRP, baseline GDF-15 and MMP-8 were also associated with future pneumonia risk independent of clinical comorbidities. The association of these inflammatory markers with future pneumonia risk has not been previously reported. Although our methodological approach was structured to maximize pneumonia risk prediction, and not establish causal pathways, we review below previously described biological roles of the biomarkers associated with pneumonia risk to guide potential future research. GDF-15 is a transforming growth factor-ß superfamily cytokine that is expressed in various cell types, and its expression increases in response to tissue injury and inflammation [26]. In subjects with chronic obstructive pulmonary disease (COPD) and in response to cigarette smoke exposure, GDF-15 is upregulated in airway epithelial cells and activates pathways that promote mucin production and cellular senescence, altering mucosal immunity and leading to chronic airway inflammation, respectively. Furthermore, GDF-15 overexpression has been shown to increase susceptibility to and severity of respiratory viral infections via increased viral replication and exaggerated virus-induced inflammation [27]. Baseline serum GDF-15 levels have been found to be an independent predictor of adverse outcomes in COPD patients, including accelerated pulmonary function decline, higher yearly exacerbation rates and increased mortality [28]. Finally, increased GDF-15 expression in the muscle of COPD subjects suggests it contributes to muscle wasting in COPD [29]. Although these studies report the association of GDF-15 with increased susceptibility to viral infection, chronic inflammation and senescence in the context of COPD and cigarette smoke exposure, they introduce possible mechanisms by which baseline GDF-15 increases future pneumonia risk. MMP-8 is a member of the matrix metalloproteinase family of proteases and is involved in the breakdown of extracellular matrix and cleavage of chemokines [30]. MMPs contribute to the pathogenesis of multiple inflammatory pulmonary conditions such as emphysema, COPD, bronchiectasis and acute respiratory distress syndrome through primarily increased neutrophilic activity [31-34]. Hartog et al. [35] found elevated levels of MMP-8 in the bronchoalveolar lavage fluid of patients with hospital-acquired pneumonia, with levels correlated to clinical severity. Our findings add to Hartog et al.’s by demonstrating that baseline MMP-8 levels may predict future pneumonia. Potential mechanisms linking baseline levels of CRP, GDF-15 and MMP-8 to pneumonia risk require further study.

Our study also examined the association of clinical comorbidities and lifestyle factors with pneumonia susceptibility, and our finding that age, atrial fibrillation, chronic pulmonary disease, diabetes, heart failure, prior overnight hospitalization and smoking were associated with pneumonia risk is similar to that of previous studies [4,23,36] providing face validity to our results. In contrast to prior studies [4,37], we did not observe chronic renal disease, cancer, BMI, alcohol use and nursing home residence to be independently predictive of future pneumonia risk. The adjustment for additional clinical covariates and potential confounders in our study as compared with prior studies, which may be partially collinear with some comorbidities, as well as the use of competing risk statistical models, may explain some of the differences between our results and prior studies of pneumonia risk factors. Despite studies showing an association between cognitive impairment and pneumonia [4,38], we interestingly did not find mini-mental status scores to be predictive of future pneumonia. Regardless, clinical comorbidities increase pneumonia susceptibility through a variety of mechanisms, in addition to chronic inflammation, that increase pathogen burden and impair immune responses. A more detailed description of links between clinical comorbidities and pneumonia risk can be found in Quinton et al. [3].

Our study has several limitations. Our current prediction model is not intended to be used for clinical decision-making, however it highlights a novel association between inflammatory markers and future pneumonia. We chose to study FHS participants because they are a cohort of well characterized, community-dwelling residents with available serum protein immunoassay measurements and linked CMS data; however, we acknowledge the lack of generalizability of our findings to other dissimilar populations. Thus, our findings of associations between individual immunoassays and pneumonia risk are hypothesis-generating and should be replicated in additional datasets prior to further exploration of potential mechanisms between these biomarkers and pneumonia susceptibility. Furthermore, because we used the most sensitive definition of pneumonia for analysis, we found biomarkers that predicted future pneumonia of various types (i.e., community-acquired, hospital-acquired, bacterial, viral, etc.). We also lacked data on the pneumonia-causing pathogens. Thus, future studies examining more specific definitions of pneumonia and causative pathogens might improve our understanding of the pathogenesis of different types of pneumonia. We also did not include the effect of vaccination, such as vaccination against pneumococcus or influenza, immunosuppressive medications, other infections or time-varying factors on pneumonia risk. In addition, studies have shown that the microbiome plays an important role in host immunity and that dysbiosis leads to impaired immune defenses [39] and potentially chronic inflammation. Future studies should explore potential relationships between dysbiosis, inflammation and pneumonia risk. Missing data on specific biomarkers (e.g., IL-6), use of multiple imputation and multiple hypothesis tests may increase risk of both systematic and random bias. However, we used methods designed to minimize bias such as multiple imputation, and LASSO covariate selection, and multiple sensitivity analyses demonstrated robustness of our models. Finally, misclassification of pneumonia diagnoses by ICD-9-CM codes may increase risk of bias.

In conclusion, we found clinical and serum protein predictors of pneumonia risk. Our findings showing associations between biomarkers of inflammation and pneumonia support the hypothesis that pneumonia is, in part, a disease of persistent or increased susceptibility. The exact mechanisms by which aging, inflammation and clinical comorbidities increase pneumonia susceptibility remain fertile areas for further investigation and potential targets for novel interventions to reduce pneumonia risk.

## Data Availability

Data cannot be shared publicly because of pre-existing data use rules surrounding Center of Medicare and Medicaid Services and Framingham Heart Study data.

## Author Contributions

A.J.W., K.S., and J.P.M. contributed to the concept and design of the work. A.J.W., K.S., Y.Z., and B.K. contributed to the acquisition and analysis of data. A.J.W., K.S., J.P.M. and M.L. contributed to the interpretation of data. M.L. and A.J.W. wrote the initial draft. All authors contributed to revising the manuscript and approved the final manuscript.

## Financial/nonfinancial Disclosures

None of the authors have any conflicts of interest to disclose.

## Funding Source(s)

Funding for this project was provided through Boston University Evans Department of Medicine Affiliate Research Collaborative; J.P.M. was supported by NIH grant R35 HL135756, A.J.W. was supported by NHLBI R01 HL136660 and R01 HL139751.

## Disclaimer

The FHS is conducted and supported by the National Heart, Lung, and Blood Institute (NHLBI) in collaboration with Boston University (Contract No. N01-HC-25195 and HHSN268201500001I). This manuscript was not prepared in collaboration with investigators of the FHS and does not necessarily reflect the opinions or views of the FHS, Boston University, or NHLBI. Additional funding for SABRe was provided by Division of Intramural Research, NHLBI, and Center for Population Studies, NHLBI. Funding support for the Framingham Multiplexed Immunoassay Panel 1-17 (SABRe project) datasets was provided by FHS contract supplement, NHLBI Intramural Funding.

## ACKNOWLEDGMENTS

M.L., A.J.W., Y.Z. and B.K. are responsible for the content of the manuscript, including the integrity of the data and the accuracy of analysis.

*   Received December 7, 2023.
*   Revision received December 7, 2023.
*   Accepted December 9, 2023.


*   © 2023, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/)

## REFERENCES

1.  1.WHO. Global Health Estimate: Life expectancy and leading causes of death and disability. [who.int/data/gho/data/themes/mortality-and-global-health-estimates](https://who.int/data/gho/data/themes/mortality-and-global-health-estimates). Date last accessed: November 5, 2021.
    
    

2.  2.Ramirez JA, Wiemken TL, Peyrani P, et al. Adults hospitalized with pneumonia in the United States: incidence, epidemiology, and mortality. Clin Infect Dis. 2017;65(11):1806– 1812.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/cix647&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F09%2F2023.12.07.23299678.atom) 

3.  3.Quinton LJ, Walkey AJ, Mizgerd JP. Integrative physiology of pneumonia. Physiol Rev. 2018;98(3):1417–1464.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1152/physrev.00032.2017&link_type=DOI) 

4.  4.Torres A, Peetermans WE, Viegi G, Blasi F. Risk factors for community-acquired pneumonia in adults in Europe: a literature review. Thorax. 2013;68(11):1057–1065.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToidGhvcmF4am5sIjtzOjU6InJlc2lkIjtzOjEwOiI2OC8xMS8xMDU3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMTIvMDkvMjAyMy4xMi4wNy4yMzI5OTY3OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

5.  5. Dela Cruz CS, Wunderink RG, Christiani DC, et al. Future research directions in pneumonia. NHLBI working group report. Am J Respir Crit Care Med. 2018;198(2):256– 263.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1164/rccm.201801-0139WS&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F09%2F2023.12.07.23299678.atom) 

6.  6.Boyd AR, Orihuela CJ. Dysregulated inflammation as a risk factor for pneumonia in the elderly. Aging Dis. 2011;2(6):487–500.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22288022&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F09%2F2023.12.07.23299678.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000208950700006&link_type=ISI) 

7.  7.Hinojosa E, Boyd AR, Orihuela CJ. Age-associated inflammation and toll-like receptor dysfunction prime the lungs for pneumococcal pneumonia. J Infect Dis. 2009;200(4):546– 554.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/600870&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19586419&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F09%2F2023.12.07.23299678.atom) 

8.  8.Yende S, Tuomanen EI, Wunderink R, et al. Preinfection systemic inflammatory markers and risk of hospitalization due to pneumonia. Am J Respir Crit Care Med. 2005;172(11):1440–1446.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1164/rccm.200506-888OC&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16166617&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F09%2F2023.12.07.23299678.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000233697300015&link_type=ISI) 

9.  9.Feinleib M, Kannel WB, Garrison RJ, McNamara PM, Castelli WP. The Framingham offspring study. Design and preliminary data. Prev Med. 1975;4(4):518–525.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0091-7435(75)90037-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1208363&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F09%2F2023.12.07.23299678.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1975BD40000012&link_type=ISI) 

10. 10.Walkey AJ, Pencina KM, Knox D, et al. Five-year risk of mechanical ventilation in community-dwelling adults: the Framingham-intermountain anticipating life support study. J Am Geriatr Soc. 2015;63(10):2082–2088.
    
    

11. 11.Ho JE, Lyass A, Courchesne P, et al. Protein biomarkers of cardiovascular disease and mortality in the community. J Am Heart Assoc. 2018;7(14):e008108.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NToiYWhhb2EiO3M6NToicmVzaWQiO3M6MTI6IjcvMTQvZTAwODEwOCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzEyLzA5LzIwMjMuMTIuMDcuMjMyOTk2NzguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

12. 12.Johannisson A, Jonasson R, Dernfalk J, Jensen-Waern M. Simultaneous detection of porcine proinflammatory cytokines using multiplex flow cytometry by the xMAP technology. Cytometry A. 2006;69A:391–395.
    
    

13. 13.Drahos J, Vanwormer JJ, Greenlee RT, Landgren O, Koshiol J. Accuracy of ICD-9-CM codes in identifying infections of pneumonia and herpes simplex virus in administrative data. Ann Epidemiol. 2013;23(5):291–293.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.annepidem.2013.02.005&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23522903&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F09%2F2023.12.07.23299678.atom) 

14. 14.Kannel WB, Sorlie P. Some health benefits of physical activity. The Framingham study. Arch Intern Med. 1979;139(8):857–861.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/archinte.1979.03630450011006&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=464698&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F09%2F2023.12.07.23299678.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1979HG12100003&link_type=ISI) 

15. 15.Folstein MF, Folstein SE, McHugh PR. “Mini-mental state.” A practical method for grading the cognitive state of patients for the clinician. J Psychiatr Res. 1975;12(3):189– 198.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0022-3956(75)90026-6&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1202204&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F09%2F2023.12.07.23299678.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1975AY27900005&link_type=ISI) 

16. 16.Katz S, Akpom CA. A measure of primary sociobiological functions. Int J Health Serv. 1976;6(3):493–508.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2190/UURL-2RYU-WRYD-EY3K&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=133997&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F09%2F2023.12.07.23299678.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1976BX47800008&link_type=ISI) 

17. 17.Jette AM, Branch LG. The Framingham disability study: II. Physical disability among the aging. Am J Public Health. 1981;71(11):1211–1216.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2105/AJPH.71.11.1211&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7294262&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F09%2F2023.12.07.23299678.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1981MN83800004&link_type=ISI) 

18. 18.Tibshirani R. Regression shrinkage and selection via the lasso. J R Stat Soc B. 1996;58(1):267–288.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18637/jss.v033.i01&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20808728&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F09%2F2023.12.07.23299678.atom) 

19. 19.Fine JP, Gray RJ. A proportional hazards model for the subdistribution of a competing risk. J Am Stat Assoc. 1999;94(446):496–509.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2307/2670170&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000081058500019&link_type=ISI) 

20. 20.Harrell FE, Lee KL, Mark DB. Multivariable prognostic models: issues in developing models, evaluating assumptions and adequacy, and measuring and reducing errors. Stat Med. 1996;15(4):361–387.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/(SICI)1097-0258(19960229)15:4<361::AID-SIM168>3.0.CO;2-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8668867&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F09%2F2023.12.07.23299678.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1996TY77400003&link_type=ISI) 

21. 21.Almirall J, Bolíbar I, Toran P, et al. Contribution of C-reactive protein to the diagnosis and assessment of severity of community-acquired pneumonia. Chest. 2004;125(4):1335– 1342.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1378/chest.125.4.1335&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15078743&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F09%2F2023.12.07.23299678.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000221793700027&link_type=ISI) 

22. 22.Meyer KC, Ershler W, Rosenthal NS, Lu XG, Peterson K. Immune dysregulation in the aging human lung. Am J Respir Crit Care Med. 1996;153(3):1072–1079.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1164/ajrccm.153.3.8630547&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8630547&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F09%2F2023.12.07.23299678.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1996TY52700029&link_type=ISI) 

23. 23.Jain S, Self WH, Wunderink RG, et al. Community-acquired pneumonia requiring hospitalization among U.S. adults. N Engl J Med. 2015;373:415–427.
    
    

24. 24.Shivshankar P, Boyd AR, Le Saux CJ, Yeh I-T, Orihuela CJ. Cellular senescence increases expression of bacterial ligands in the lungs and is positively correlated with increased susceptibility to pneumococcal pneumonia. Aging Cell. 2011;10(5):798–806.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1474-9726.2011.00720.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21615674&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F09%2F2023.12.07.23299678.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000294864400007&link_type=ISI) 

25. 25.Boyd AR, Shivshankar P, Jiang S, Berton MT, Orihuela CJ. Age-related defects in TLR2 signaling diminish the cytokine response by alveolar macrophages during murine pneumococcal pneumonia. Exp Gerontol. 2012;47(7):507–518.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.exger.2012.04.004&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22548913&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F09%2F2023.12.07.23299678.atom) 

26. 26.Wang TJ, Wollert KC, Larson MG, et al. Prognostic utility of novel biomarkers of cardiovascular stress: the Framingham heart study. Circulation. 2012;126(13):1596–1604.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTQ6ImNpcmN1bGF0aW9uYWhhIjtzOjU6InJlc2lkIjtzOjExOiIxMjYvMTMvMTU5NiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzEyLzA5LzIwMjMuMTIuMDcuMjMyOTk2NzguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

27. 27.Verhamme FM, Freeman CM, Brusselle GG, Bracke KR, Curtis JL. GDF-15 in pulmonary and critical care medicine. Am J Respir Cell Mol Biol. 2019;60(6):621–628.
    
    

28. 28.Husebø GR, Grønseth R, Lerner L, et al. Growth differentiation factor-15 is a predictor of important disease outcomes in patients with COPD. Eur Respir J. 2017;49(3):1601298.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiZXJqIjtzOjU6InJlc2lkIjtzOjEyOiI0OS8zLzE2MDEyOTgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8xMi8wOS8yMDIzLjEyLjA3LjIzMjk5Njc4LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

29. 29.Patel MS, Lee J, Baz M, et al. Growth differentiation factor-15 is associated with muscle mass in chronic obstructive pulmonary disease and promotes muscle wasting in vivo. J Cachexia Sarcopenia Muscle. 2016;7(4):436–448.
    
    

30. 30.Thirkettle S, Decock J, Arnold H, Pennington CJ, Jaworski DM, Edwards DR. Matrix metalloproteinase 8 (collagenase 2) induces the expression of interleukins 6 and 8 in breast cancer cells. J Biol Chem. 2013;288(23):16282–16294.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjEyOiIyODgvMjMvMTYyODIiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8xMi8wOS8yMDIzLjEyLjA3LjIzMjk5Njc4LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

31. 31.Segura-Valdez L, Pardo A, Gaxiola M, Uhal BD, Becerril C, Selman M. Upregulation of gelatinases A and B, collagenases 1 and 2, and increased parenchymal cell death in COPD. Chest 2000;117(3):684–694.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1378/chest.117.3.684&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10712992&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F09%2F2023.12.07.23299678.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000085884000015&link_type=ISI) 

32. 32.Sepper R, Konttinen YT, Ding Y, Takagi M, Sorsa T. Human neutrophil collagenase (MMP-8), identified in bronchiectasis BAL fluid, correlates with severity of disease. Chest. 1995;107(6):1641–1647.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1378/chest.107.6.1641&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7781360&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F09%2F2023.12.07.23299678.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995RD45100031&link_type=ISI) 

33. 33.Torii K, Iida K, Miyazaki Y, et al. Higher concentrations of matrix metalloproteinases in bronchoalveolar lavage fluid of patients with adult respiratory distress syndrome. Am J Respir Crit Care Med. 1997;155(1):43–46.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1164/ajrccm.155.1.9001287&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9001287&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F09%2F2023.12.07.23299678.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1997WC67000009&link_type=ISI) 

34. 34.Betsuyaku T, Nishimura M, Takeyabu K, et al. Neutrophil granule proteins in bronchoalveolar lavage fluid from subjects with subclinical emphysema. Am J Respir Crit Care Med. 1999;159(6):1985–1991.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1164/ajrccm.159.6.9809043&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10351949&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F09%2F2023.12.07.23299678.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000080780000045&link_type=ISI) 

35. 35.Hartog CM, Wermelt JA, Sommerfeld CO, Eichler W, Dalhoff K, Braun J. Pulmonary matrix metalloproteinase excess in hospital-acquired pneumonia. Am J Respir Crit Care Med. 2003;167(4):593–598.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1164/rccm.200203-258OC&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12588713&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F09%2F2023.12.07.23299678.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000181097900017&link_type=ISI) 

36. 36.Krumholz HM. Post-hospital syndrome--an acquired, transient condition of generalized risk. N Engl J Med. 2013;368:100–102.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMp1212324&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23301730&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F09%2F2023.12.07.23299678.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000313238100013&link_type=ISI) 

37. 37.Baik I, Curhan GC, Rimm EB, Bendich A, Willett WC, Fawzi WW. A prospective study of age and lifestyle factors in relation to community-acquired pneumonia in US men and women. Arch Intern Med. 2000;160(20):3082–3088.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/archinte.160.20.3082&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11074737&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F09%2F2023.12.07.23299678.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000165181000007&link_type=ISI) 

38. 38.Shah FA, Pike F, Alvarez K, et al. Bidirectional relationship between cognitive function and pneumonia. Am J Respir Crit Care Med. 2013;188(5):586–592.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1164/rccm.201212-2154OC&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23848267&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F09%2F2023.12.07.23299678.atom) 

39. 39.Schuijt TJ, Lankelma JM, Scicluna BP, et al. The gut microbiota plays a protective role in the host defence against pneumococcal pneumonia. Gut. 2016;65(4):575–583.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ3V0am5sIjtzOjU6InJlc2lkIjtzOjg6IjY1LzQvNTc1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMTIvMDkvMjAyMy4xMi4wNy4yMzI5OTY3OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=)