2.1 Study Design and population
A prospective cohort study was conducted using the Atherosclerosis Risk in Communities (ARIC) study population from four communities in the United States17. A total of 15,792 participants, aged between 45 and 64 years, were recruited in 1987–1989 (Visit 1) and followed up every three years. Visit 2 was conducted in 1990–1992, Visit 3 in 1993–1995 and Visit 4 in 1996–1998. The ARIC study was approved by the institutional review committee of all participating institutions with the informed consent from all participants. Functional limitations were assessed for the first time in Visit 4. In our study, Visit 1 was regarded as the baseline, part of the data from Visit 2 was collected, and Visit 4 was regarded as the research outcome. Among participants in Visit 4 (n = 15,028), exclusion criteria were as follows: absence of an inflammatory biomarkers assessment (n = 336), absence of one or more covariates (n = 1111), and absence of an assessment on functional limitations (n = 3537). Finally, a total of 10,044 participants were enrolled in this study.
2.2 Inflammatory biomarkers
The physiological response of biomarkers to stress varies across individuals. One biomarker may be compensated by the effective function of other regions, and may not accurately describe the overall physiological function and the body's response to stress18. It is unclear whether a single measurement can adequately capture inflammation chronicity. Thus, the use of an inflammation composite score alleviated this possibility. Therefore, four biomarkers in plasma collected from Visit 1 were used to describe the state of systemic inflammation, including white blood cell count (after log-transformed for corrected skewness), von Willebrand factor, fibrinogen and Factor VIII. The inflammation composite score refers to the average value of four biomarkers after standardization to z-score. The selection of these markers is largely dependent on their availability in the ARIC study and that the inflammation composite score was often used in other studies19.
In this study, hypersensitive-C-reactive proteins (hs-CRP) in blood sample from Visit 1 and 2 was used as a sub-indicator of systemic inflammation. There was no difference between the relevant conditions of the two measurements. The hs-CRP was classified as “low” or “elevated”, with a cutoff value of 3.0 mg/L20. According to the hs-CRP values from Visit 1 to Visit 2, the following four longitudinal patterns of change were established: (1) Consistent low hs-CRP: both were at low level; (2) Ascending hs-CRP: first at low level and then at elevated level; (3) Descending hs-CRP: first at elevated level and then at low level; (4) Consistent elevated hs-CRP: both were at elevated level.
2.3 Functional limitations
At Visit 4, the functional status assessment was completed in the form of a self-reported questionnaire containing the designated activities, as was done in the third National Health and Nutrition Examination Survey4. These activities included activities of daily living (eating, dressing, getting up, and walking)21, instrumental activities of daily living (cooking, housework, and financial management)22 and lower limb functions (standing up from a chair without arm support, bending or kneeling, walking 1/4 miles, going upstairs, and carrying 10 pounds)23. Participants were required to indicate the level of difficulty in conducting these activities. When participants could not answer for themselves, proxies of participants could answer those questions with high degree of reliability. The final results were reported as "no limits" or "function impaired". Functional limitations were defined as three categories: impaired activities of daily living (ADLs), impaired instrumental activities of daily living (IADLs) and impaired lower limb function (LEF).
2.4 Assessment of Multiple Organ Function
Organ function was assessed at Visit 4, including cardiac function, brain function, lung function, liver function, and kidney function. Cardiac function was assessed using N-Terminal Pro-Brain Natriuretic Peptide (NT-proBNP), which was measured in plasma using electrochemiluminescence immunoassay on an automatic Cobase411 analyzer (Roche diagnosis). Brain function was assessed using the global composite cognition score, including delayed word recall test (DWRT), digit symbol substitution test (DSST), and word fluency test (WFT). The composite cognition score was the average value of the three tests after standardization to z-score24–27. Lung function was assessed using the ratio of forced expiratory volume in one second / forced vital capacity (FEV1/FVC). The values were measured by the standardized Collins Survey II spirometer28. Liver function was assessed using plasma alaninetransaminase (ALT). Kidney function was assessed using estimated glomerular filtration rate (eGFR) which was calculated using creatinine in the Eq. 29.
2.5 Baseline Covariates
At baseline, the following covariates were assessed based on self-report or medical record evidence from participants: age, gender, race, income, education, hours of metabolic equivalent of task (MET-hour) per week, prevalence of medical comorbidity (hypertension, diabetes, coronary heart disease, heart failure, cancer, and chronic obstructive pulmonary disease), and medication use. Total cholesterol, high density lipoprotein, low density lipoprotein and triglycerides were measured by enzymatic analysis.
Hypertension was defined as blood pressure ≥ 140/90 mmHg, or a history of medication. Diabetes was defined as fasting glucose ≥ 126.0 mg/dl or non-fasting glucose ≥ 200.0 mg/dl, or a history of medication or insulin therapy. Coronary heart disease was defined as acute coronary syndrome, or chronic coronary artery disease, or a history of medication or coronary revascularization. Chronic obstructive pulmonary disease was identified as FEV1/FVC < 70% after bronchodilator, or a history of medication. Heart failure and cancer were diagnosed with the evidence from medical records.
2.6 Statistical analysis
In this study, continuous variables were represented with median (25th − 75th) or mean ± standard deviation, and categorical variables were represented with numbers (percentages). Continuous variables were compared using Mann–Whitney U test or ANOVA, and categorical variables were compared using chi-square test.
Logistic regression models were conducted to assess the association between the inflammation composite score at Visit 1, hs-CRP at Visit 1, and longitudinal pattern of hs-CRP (from Visit 1 to Visit 2) with functional limitations (impaired ADLs, IADLs and LEF), respectively. Two regression models were established to examine the independent role of these relationships. Model 1 was adjusted according to age, gender, race, education, income, and MET-hour/week. Model 2 was additionally adjusted according to the prevalent of hypertension, diabetes, coronary heart disease, heart failure, cancer, chronic obstructive pulmonary disease, total cholesterol, high density lipoprotein, low density lipoprotein, triglycerides, cholesterol-lowering and anti-inflammatory medications. Furthermore, a logistic regression analysis in Model 2 was conducted to assess the effect modification between the inflammation composite score and functional limitations of age, gender and race subgroups. The value of multiplicative interaction term < 0.05 was considered statistically significant.
According to the procedure recommended by Hayes30, Pathway analysis was carried out through the structural equation Model 4 to calculate the indirect effects of the inflammation composite score on functional limitations, including cardiac function (NT-proBNP was log-transformed to correct for skewness), brain function (composite cognition score), lung function (FEV1/FVC), liver function (ALT), and kidney function (eGFR). In structural equation modeling, cross-products from estimated mediation effects were considered statistically significant when confidence intervals (CIs) did not include zero. Our study used the standardized regression coefficients (β) to report point estimates of direct and indirect effects.
A two-tailed P value of < 0.05 was considered statistically significant. SPSS Statistics (version 36.0, IBM Corp, Armonk, USA), PROCESS (version 3.5 for SPSS) and R software (version 3.5.0, Vienna, Austria) were used for data analysis.