Abstract
BACKGROUND/OBJECTIVES
Limited evidence suggests that the inflammatory state of aging is a risk factor for accelerated renal function (RF) decline. We examined this hypothesis using inflammatory biomarkers and RF collected over a 9-year follow-up in relatively healthy individuals enrolled in the InCHIANTI study.
DESIGN
Longitudinal
SETTING
Community
PARTICIPANTS
Participants included 687 adults age 60 and older with baseline estimated glomerular filtration rate (eGFR)≥60 ml/min/1.73m2 and no diabetes mellitus (DM).
MEASURES
eGFR, proxy for RF, was determined using the CKD-EPI equation at baseline and then at 3-, 6-, and 9-year follow-ups. Incident chronic kidney disease (CKD) was defined as new-onset eGFR<60 ml/min/1.73m2 at each follow-up. Predictors included baseline and time-dependent inflammatory biomarkers: soluble tumor necrosis factor alpha receptors (sTNFα-R1 and -R2), interleukin(IL)-6, IL-18, IL-1β, IL-1 receptor antagonist, and high sensitivity C-reactive protein.
Results
Higher baseline sTNFα-R1 was significantly associated with lower eGFR over 9 years, independent of DM or blood pressure:(Baseline:β̂=−0.39, p = 0.001; 3-year: β̂=−0.26, p=0.001; 6-year: β̂=−0.36, p=0.001; 9-year: β̂=−0.47, p=0.001). The rate of TNFα-R1 change was significantly associated with rate of eGFR change (β̂=−0.18, p=0.001). Baseline sTNFα-R1 predicted incident CKD [per 1-SD increment]: 3-year: relative risk (RR)=1.3 (95% CI: 1.1–1.5); 6-year: RR=1.5 (1.1–2.2); 9-year RR=1.6 (1.1–2.2,). Similar results were found for sTNFα-R2.
Conclusion
Baseline and the rate of change in TNFα receptors were significantly associated with faster RF decline and incident CKD among older adults independent of DM or blood pressure.
Keywords: Cytokines, inflammation, renal function, CKD, aging
Introduction
Chronic kidney disease (CKD) is a major health burden worldwide.1 Diabetes mellitus (DM) and essential hypertension (HTN) are established risk factors for CKD,2 yet they do not fully explain the increased risk of renal function (RF) decline and CKD with aging.3, 4 The role of high concentrations of pro-inflammatory cytokines has drawn considerable attention because studies have shown that patients with known CKD and high inflammation experience accelerated declines in RF.5, 6 RF also declines with aging in the absence of major comorbidities.7–10 Moreover, chronic inflammation is known as a feature of aging.11 There is body of evidence showing associations between inflammatory biomarkers, which activate two main inflammatory pathways, and age-related diseases such as functional decline, mobility and cognititive impairments, and mortality.11–14 Examples of such biomarkers include interleukin-6 (IL-6) and tumor necrosis factor α (TNF α) that activate nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB), and IL-1β and IL-18 that are related to the nucleotide-binding oligomerization domain (NOD)-like receptor containing pyrin domain 3 (NLRP3) pathways.
To date, the results of studies that show inflammation predicts incident CKD are still controversial because of methodological limitations, including limited number of inflammatory biomarkers, not considering potentially informative censoring and loss to follow-up, or limited follow-up time.14–18 Moreover, most of these studies were focused on patients with DM. In this study we focused on IL-6 and two soluble TNFα receptors (sTNFα-R1 and sTNFα-R2) from NF-κB pathway, IL-18 and IL-1β from NLRP3 pathway, IL-1 receptor antagonist, and HS-CRP.14
Using InCHIANTI data,19 we tested if baseline levels of these inflammatory biomarkers and their increased levels over 9 years were prospectively associated with accelerated RF decline and CKD development in nondiabetic participants who initially had normal RF.
Methods
Study Design and Setting
Analyses were performed using 9 years of prospective data from the InCHIANTI (Invecchiare in Chianti) study. Participants were enrolled from 1998–2000 and had follow-up visits every 3 years. The 3-year follow-up occurred in 2001–2003, the 6-year in 2004–2006, and the 9-year in 2007–2009. All participants explicitly consented that their collected data could be used for research purposes.
Participants
Of 1,453 individuals, 1,326 had measured baseline serum creatinine. Exclusion criteria for these analyses included age<60 years (n=237); baseline estimated glomerular filtration rate (eGFR)<60 ml/min/1.73m2 (n=232); DM defined as fasting plasma glucose≥126 mg/dL, history of DM or antidiabetic medication (n=117); history of steroid utilization (n=15); high sensitive C-reactive protein (HS-CRP) ≥30 g/L (n=9) to exclude persons with severe infections or autoimmune diseases that resulted in exclusion of the outliers of other biomarkers as well; and missing baseline TNFα-R1 (n=29). Because some individuals met multiple exclusion criteria, the baseline sample size was 687.
Measurements
Inflammatory biomarkers
Measured serum levels of inflammatory biomarkers included sTNFα-R1, sTNFα-R2, IL-18, and IL-6 at baseline and all three follow-ups, while IL-1 receptor antagonist (IL-1RA) and IL-1β were measured only at baseline, using previously described assays.20, 21 IL-1β and IL-1RA were measured in duplicate using high-sensitivity enzyme-linked immunoabsorbent assays (ELISA) (BIOSOURCE International, Camarillo, CA). IL-18 was measured using highly sensitive quantitative sandwich assays (Quantikine HS, R&D Systems, Minneapolis, MN). sTNFα-R1 and sTNFα-R2 were measured utilizing ELISA (BIOSOURCE International, Camarillo, CA). HS-CRP was measured using colorimetric competitive immunoassay with purified protein and polyclonal anti-CRP antibodies (Calbiochem, San Diego, CA). We calculated the inter-assay and intra-assay coefficient of variation (CV) at baseline and all follow-ups (Table S1). Additionally, to understand “within-subject” variations of the inflammatory biomarkers in the total study population and in those who never developed CKD, we performed a mixed-effect model for each longitudinal biomarker over 9-year follow-ups and reported random estimated within-subject variations with 95% confidence interval (95% CI) (Table S2). From InCHIANTI genotyped data we also extracted a genetic variant for R92Q mutation (rs4149584) within the tumor necrosis factor receptor superfamily member 1A (TNFRSF1A), a gene that encodes TNFα-R1,22 as this variant increases the affinity of TNFα to TNFα-R1.23
Covariates
Demographic covariates included age; sex; smoking, reported as “never,” “former,” or “current”; marital status, reported as “not married,” “married,” or “widowed,”; and education (number of school years). Systolic blood pressure (SBP), measured three times in both arms, was then determined as the average of the three maximum measures. Similarly, diastolic blood pressure (DBP) was determined. Use of angiotensin converting enzyme inhibitor (ACEI), angiotensin II receptor blockers (ARBs), or nonsteroidal anti-inflammatory drugs and presence of cancer were self-reported. We adjusted for new-onset DM over the time course. A diagnosis of congestive heart failure (CHF) was adjudicated based on self-reported history and/or current treatment for CHF. Baseline urinary protein (normal range for proteinuria was 0–150mg/dL)24 was measured using a Modular P800 Hitachi Analyzer and normalized to urine creatinine concentration, and the details of urine creatinine measure is described elsewhere.25, 26 Creatinine levels depend on muscle mass, which often declines with aging; therefore, we included calf muscle cross-sectional area (proxy for muscle mass) derived from computerized tomography and changes in cross-sectional area over time.27
Outcomes
Serum creatinine was quantified using an enzymatic assay and calibrated to standardized isotope dilution mass spectrometry (IDMS).28 Glomerular filtration rate (eGFR), a proxy for RF, was estimated from IDSM-traceable serum creatinine, age, and sex using the Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) equation (henceforth called eGFR).29, 30 We also used the rate of eGFR change over 9 years as an outcome (slope, see Statistical Analysis). Moreover, we defined incident CKD as the development of eGFR<60 ml/min/1.73m2 in individuals who had normal eGFR (≥60 ml/min/1.73m2) at the previous visit.
Statistical analyses
Participants’ characteristics are reported as means±standard deviation (SD) or median and intesrquartile range (IQR), percentages of total study population, and by incident CKD status (Table 1).
Table 1.
Baseline characteristics of total study participants, those with incident CKD (eGFR <60 anytime during the follow-up) vs. no-CKD (eGFR ≥60 ml/min/1.73 m2).
Baseline | Total n = 687 |
No CKD ≥ 60, n = 434a |
Incident CKD < 60, n = 188 |
p-Value |
---|---|---|---|---|
Demographic variables | ||||
Age, mean yr ± SD | 72.9 ± 7.1 | 71.4 ± 6.4 | 73.5 ± 6.2 | 0.001 |
Sex, % | ||||
Male | 46.0 | 45.6 | 43.6 | 0.64 |
Education, mean yr ± SD | 5.6 ± 3.4 | 5.8 ± 3.2 | 5.6 ± 3.7 | 0.42 |
Marital status, % | ||||
Not married | 8.3 | 8.5 | 8.5 | |
Married | 66.1 | 71.4 | 59.6 | 0.002 |
Widowed | 25.6 | 20.1 | 31.9 | |
Smoking, % | ||||
Never | 58.1 | 58.1 | 59.0 | 0.77 |
Former | 25.8 | 25.3 | 26.6 | |
Current | 16.1 | 16.6 | 14.4 | |
NSAIDs, % | ||||
Yes | 41.9 | 41.4 | 43.3 | 0.65 |
| ||||
Physical exam variables | ||||
Systolic BP, mean ± SD | 158.8 ± 19.2 | 157.5 ± 19.4 | 162.4 ± 18.5 | 0.003 |
Diastolic BP, mean ± SD | 84.0 ± 8.4 | 83.6 ± 8.5 | 85.1 ± 8.2 | 0.04 |
ACEI, column % | 19.21 | 18.3 | 19.6 | 0.60 |
ARBs, column % | 1.86 | 1.8 | 2.39 | 0.50 |
CHF, column % | 1.03 | 0.81 | 1.44 | 0.35 |
Base muscle mass area calf mm2 ± SD | 71.1 ± 3.5 | 6286 ± 1252 | 6461 ± 1254 | 0.12 |
| ||||
Laboratory | ||||
Serum glucose, mg/dL | 88.1 ± 10.7 | 88.1 ± 10.5 | 88.7 ± 10.9 | 0.50 |
Serum creatinine, mg/dL | 0.9 ± 0.1 | 0.8 ± 0.1 | 0.9 ± 0.1 | 0.001 |
Baseline eGFR ml/min/1.73m2 | 76.4 ± 10.0 | 78.8 ± 9.8 | 71.8 ± 8.9 | 0.001 |
Annual eGFR slope ml/min/1.73 m2 | −0.32 ± 0.1 | 0.9 ± 0.1 | −7.9 ± 0.3 | 0.001 |
Inflammatory biomarkers Median (IQR) | −0.50 ± 0.6 | −0.2 ± 0.06 | −2.1 ± 0.2 | <0.001 |
sTNFα-R1 pg/mL | 1253 (1037–1526) | 1174 (977–1432) | 1402 (1135–1571) | 0.001 |
sTNFα-R2 pg/mL | 2499 (2205–2914) | 2377 (2091–2757) | 2587 (2297–2961) | 0.001 |
IL-1RA pg/mL | 130.8 (95.4–176.3) | 127.0 (94.0–169.0) | 134.3 (94.2–180.4) | 0.27 |
IL-1β pg/mL | 0.12 (0.08–0.20) | 0.12 (0.08–0.18) | 0.13 (0.09–0.2) | 0.10 |
IL-6 pg/mL | 2.7 (1.9–3.7) | 2.5 (1.9–3.5) | 2.8 (1.9–3.7) | 0.03 |
IL-18 pg/mL | 367.7 (283.6–457.2) | 360.4 (280.3–450.7) | 65.8 (284.3–440.0) | 0.73 |
HS-CRP mg/dL | 2.5 (1.3–5.1) | 2.2 (1.1–4.6) | 2.8 (1.4–5.3) | 0.03 |
Estimated glomerular filtration rate based on the Chronic Kidney Disease Epidemiology Collaboration(CKD-EPI)equation
Of 687, 42 individuals died at 3-year follow-up and 23 missed 3-year follow-up and died at 6-year follow-up, therefore, not possible to quantify incident CKD for.
p-value based on non-parametric test: Kolmogorov-Smirnov or chi-square tests
NSAIDs: Nonsteroidal anti-inflammatory drugs
BP: blood pressure
ACEI: Angiotensin converting enzyme inhibitor
ARBs: Angiotensin II receptor blockers
CHF: Congestive heart failure
sTNFα-R1: soluble tumor necrosis factor alpha-1 receptor
sTNFα-R2: soluble tumor necrosis factor alpha-2 receptor
IL: Interleukin
HS-CRP: High sensitivity C-reactive protein
IL-1RA: Interleukin-1 receptor antagonist,
HS-CRP: high sensitive C-reactive protein
Column % is reported
IQR: Inter Quartile Range
The skewed inflammatory biomarkers were natural log-transformed. Then, the biomarkers, SBP, muscle cross-sectional area, and eGFR were standardized to mean =0 and SD=1 to permit comparison of their incremental predictive associations. We used inverse probability weighted generalized estimating equations (W-GEE) to address data missing at random, after accounting for measured covariates and truncation by death.31 We assessed the relationship between standardized baseline inflammatory biomarkers and standardized eGFR by separate W-GEE models adjusted for muscle cross-sectional area, SBP, new-onset DM, and other covariates. Results are regression coefficients relating baseline biomarkers with eGFR at each follow-up. To assess which inflammatory biomarker had a stronger association with time-dependent eGFR or rate of change of eGFR over time, we—after evaluating collinearity (Pearson correlation, data not shown)—included all significant biomarkers in a single model adjusting for demographics, CHF, SBP, new-onset DM, and other covariates (see Methods). Then, a parsimonious model was generated by applying backward stepwise elimination of variables with p-value>0.1. We also examined possible threshold effects by comparing quartiles 2, 3, and 4 against quartile 1(reference group).
We computed the rates of change inflammatory biomarkers over time using mixed-effect models and operationalized time as a continuous variable obtaining individual time-slopes as annual rate of change. Similarly, we modeled eGFR change using a mixed-effect model and operationalized time as a continuous variable to obtain annual rate of eGFR (slope), a common RF outcome used in longitudinal studies. We reported the mean±SD of slope in incident CKD versus non-CKD. We determined if baseline inflammatory biomarkers or their standardized slopes predicted the rate of RF change after accounting for time varying conventional risk factors (SBP and DM) and baseline RF. We adjusted for baseline eGFR to account for inter-individual variability.
The slopes of biomarkers and eGFR were normally distributed. We used standardized slopes (using Z-transformed slopes) for all continuous measures in the models.
The frequently-used single measure of eGFR<60 was considered as an operational definition of CKD. For individuals who had eGFR<60 at one visit and eGFR≥60 at the next visit, we allowed them to re-enter the analyses in the next visit to calculate incidence rate for the subsequent round of analysis. Therefore, we treated CKD as a recurrent event. We then examined if the inflammatory biomarkers significantly predicted incident CKD within the subsequent follow-up visits using weighted generalized linear models (Poisson distribution model, lag-time) accounting for missing values, with results expressed as relative risk (RR) and 95%CI. All models accounted for covariates described above. We performed all statistical analysis in STATA 14.1(College Station, Texas 77845 USA).
Results
Participant characteristics
Baseline characteristics of the 687 participants are summarized in Table 1. The mean age was 72.9 years (range 60–102 years) and 54% were women. Particularly, participants who experienced incident CKD were older, more likely to be widowed, and had higher baseline levels of SBP, DBP, sTNFα-R1, sTNFα-R2, IL-6, and HS-CRP than participants who did not experience CKD (p<0.05)
There were less within-subject variations in estimated inflammatory biomarkers in those with no CKD experience versus those with incident CKD (Table S2). The CKD participants started with lower baseline eGFR than participants who did not experience incident CKD (p=0.001). Mean urine protein to creatinine ratio at baseline was 7.3±1.9 mg/g. Mean eGFR at baseline was 71.8±8.9 ml/min/1.73m2. The mean eGFR over 9 years for CKD was 51.9±6.6, but 77.2±9.7 ml/min/1.73m2 in non-CKD groups. The mean annual rate of change was −0.09±0.13 ml/min/1.73m2 in all participants, while it was −7.6±0.3 in CKD and 1.2±0.1 ml/min/1.73m2 in non-CKD groups (p=0.001). Further characteristics are shown in Table 1. Only two individuals developed severe loss of kidney function (eGFR 15–29ml/min/1.73m2) and 23 developed new onset DM over the 9-year follow-up. No participant carried a TNFRSF1A R92Q variant.
Baseline inflammatory biomarker and RF over time
In separate analyses, after adjustment for covariates, we found that higher baseline serum levels of sTNFα-R1 and TNFα-R2 were associated with significantly lower baseline eGFR (β̂=−0.34, p<0.001; β̂=−0.23, p=0.001, respectively; Figure 1). The association between baseline levels of sTNFα-R1 and lower eGFR remained strongly significant over the 9 -year of follow-up (p<0.001; Figure 1). A similar association pattern of baseline sTNFα-R2 with lower RF was also observed (p=0.001 for all time points; Figure 1).
Figure 1.
Effect size (standardized beta) from separate regression models examining associations of each log-transformed baseline inflammatory biomarker and longitudinal log-transformed renal function (eGFR), adjusted for covariates. sTNFα-R1 & 2: soluble tumor necrosis factor alpha-receptor 1 and 2, IL-6 (Interleukin-6), IL-1β: Interleukin-1 beta, IL-1RA: Interleukin-1 receptor antagonist, IL-18: Interleukin-18, CRP: C-reactive protein.
Other inflammatory biomarkers had less consistent associations with RF. We found that baseline IL-6, IL1-β, IL-1RA, and IL-18 were not associated with eGFR at baseline (p>0.05), but were associated with eGFR at the 6- or 9-year follow-up visits (p<0.05; Figure 1). We examined the possibility of a threshold effect by comparing quartiles 2, 3, and 4 against quartile 1 as the reference group. Compared with the first quartile (Q1), participants in the Q2, Q3, and Q4 of baseline sTNFα-R1 and sTNFα-R2 were all significantly associated with steeper declines in eGFR: (sTNFα-R1 Q2: β̂= −0.12, p=0.001; Q3: β̂ =−0.20, p<0.001; Q4: β̂==−0.21, p=0.001), and (sTNFα-R2 Q2: β̂=−0.24, p=0.001; Q3: β̂=−0.30, p<0.001; Q4: β̂=−0.28, p=0.004). Similarly, higher baseline TNFα-R1 and TNFα-R2, when considered as a continuous variable (dose-response), were associated with a steeper slope of eGFR decline (β̂=−0.27, p=0.001; β̂=−0.28, p=0.001, respectively). In contrast, the association of baseline IL-6 with the slope of eGFR was not significant when considered as either quartiles (Q2: β̂=−0.03, p=0.5; Q3: β̂=−0.003, p=0.8; Q4: β̂=0.02, p=0.6) or as a continuous variable. Similar consideration holds for the other nonsignificant biomarkers (results not shown).
In a model that included sTNFα-R1, IL-6, IL-1RA, and HS-CRP, higher levels of baseline sTNFα-R1 remained significantly associated with lower RF over the course of 9-year follow up (3-year follow-up β̂=−0.25, p=0.001; 6-year follow-up: β̂=−0.37, p=0.001; 9-year follow-up: β̂=−0.45, p=0.001). The associations of baseline IL-6, IL-1RA, and HS-CRP with RF were less consistent than that of sTNFα-R1 with RF (Table S3). In a similar fashion, we replaced sTNFα-R1 with sTNFα-R2, and baseline sTNFα-R2 remained significantly associated with lower eGFR over the course of follow-ups (3-year follow-up β̂=−0.11, p=0.03; 6-year follow-up: β̂=−0.26, p=0.001; 9-year follow-up: β̂=−0.37, p=0.001). We also found that IL-6, IL-1RA, and HS-CRP remained significant. Other inflammatory biomarkers were not included due to the backward stepwise elimination.
We estimated participant-specific time-slopes of RF (henceforth called “rate” of change; see Methods). In separate W-GEE models, including each baseline biomarker, the eGFR slope was regressed on the baseline biomarkers and adjusted for baseline RF and other covariates. Of note, baseline sTNFα-R1 and IL1-RA were significantly associated with faster declines of RF (β̂=−0.10, p<0.001; β̂=−0.05, p=0.001, respectively; Table S4). Additionally, a sensitivity analysis that removed baseline eGFR from the models resulted in even stronger estimates. In the model evaluating association of baseline TNFα-R1 and slope of eGFR unadjusted for baseline eGFR, the standardized beta coefficient was −0.25 (p<0.001).
Rates of inflammatory biomarkers and RF change (slope)
We estimated participant-specific time-slopes of inflammatory biomarkers and RF change and performed separate W-GEE models regressing standardized eGFR slope on the standardized slopes of biomarkers adjusted for baseline RF and other covariates. The accelerated increase of sTNFα-R1 was strongly associated with the accelerated decline of eGFR (β̂=0.25, p<0.001). A similar pattern was observed for the association between the rate of increase of sTNFα-R2, IL-6, and HS-CRP; and the rate of eGFR change (β̂=−0.18, p<0.001; β̂=−0.06, p=0.01; β̂=−0.07, p=0.002, respectively) in separate models (Table S5).
After adjustment for baseline eGFR and other covariates, we simultaneously entered the significant biomarkers’ slopes in the model. sTNFα-R1 and sTNFα-R2 slopes were highly correlated (r=0.98); therefore, we excluded sTNFα-R2. The rate of sTNFα-R1 change remained significantly associated with the rate of eGFR change, so that 1 SD increase of TNFα-R1 is associated with 0.25 SD faster decline in eGFR (β̂=−0.25, p<0.001; Table 2). Having CHF and using angiotensin-converting enzyme inhibitor or angiotensin II receptor blockers were not significantly associated with the rate of eGFR change (p>0.7). When we performed similar analyses for non-hypertensive individuals, both baseline association and the association of rate of sTNFα-R1 change with rate of change of eGFR remained significant (baseline: β̂=−0.23, p<0.001; rate of change: β̂=−0.15, p<0.001). Similarly, in these sensitivity analyses, the baseline sTNFα-R2 levels and the rate of sTNFα-R2 change were significantly associated with the rate of eGFR change (β̂=−0.06, p=0.001; β̂=−0.13 p<0.001, respectively).
Table 2.
Association between the rates of all significant inflammatory biomarkers simultaneously entered in the model (slope of change of biomarkers over time) and the rate of renal function (slope of eGFR over time).
Slope of eGFR over 9 yearsa |
||
---|---|---|
β̂ | p | |
TNFα-R1 | −0.25(0.03) | <0.001 |
IL-6 | −0.03(0.03) | 0.50 |
HS_CRP | −0.01(0.02) | 0.49 |
DM | −0.08(0.05) | 0.57 |
SBP | −0.05(0.054) | 0.02 |
Slope of eGFR: slope of estimated rate of glomerular filtration, the model is adjusted for baseline eGFR, age, sex, marital status, education, site of study, cancer, NSAIDs use, calf muscle area, SBP, new-onset DM, baseline systolic blood pressure (SBP); slopes of eGFR and all inflammatory biomarkers, calf muscle area and SBP are standardized; β̂: standardized beta regression coefficient, presented per SD unit increase; SD: standard deviation
Inflammatory biomarkers and incident CKD
We assessed whether sTNFα-R1 at previous visits predicted incident CKD at 3-, 6-, and 9-year follow-ups. Participants with high levels of sTNFα-R1 at a given visit (baseline, 3-year, or 6-year) had a significantly higher risk of incident CKD at a subsequent visit. Additionally, relative risks (RR) for incident CKD per SD increment of baseline TNFα-R1 were: 3-year RR: 1.3 (95%CI: 1.1 to 1.5), 6-year RR: 1.5 (95%CI: 1.1 to 2.2), and 9-year follow-ups RR: 1.6, (95%CI: 1.1 to 2.3). Similarly, SD increments of TNFα-R1 measured at the 3-year and 6-year visits were significantly associated with incident CKD at later visits (Table 3). We found similar results for TNFα-R2 and incident CKD. No other inflammatory biomarker significantly predicted incident CKD (Table S6).
Table 3.
Association of TNFα-R1 at each visit with incident CKD at subsequent Follow-up visits.
Time of TNFα-R1 Measurement |
Year 3 Incident CKD RR (95% C |
Year 6 Incident CKD RR (95% CI) |
Year 9 Incident CKD RR (95% CI) |
---|---|---|---|
Baseline | 1.3(1.1–1.5) | 1.5(1.1–2.2) | 1.6(1.1–2.3) |
Year 3 | - | 2.0(1.5–2.6) | 1.5(1.1–2.2) |
Year 6 | - | - | 1.3(1.07–1.7) |
Standardized TNFα-R1, RR: relative risk, CKD is defined as new development of eGFR <60 ml/min/1.73 m2 at subsequent follow-up visit, adjusted for sex, marital status, education, site of study, and time relevant age, cancer, NSAID use, calf muscle area, SBP, new-onset DM
Discussion
Our study shows that in the group of nondiabetic community-dwelling older adults with initially normal RF, higher baseline sTNFα receptors are associated with accelerated RF decline, and a greater likelihood of incident CKD. This association remained strong over the 9-year follow-up and was independent of SBP and new-onset DM. We also showed that greater increases of these inflammatory biomarkers were associated with faster decline of RF, independent of baseline RF, SBP, new-onset DM, and CHF.
Our results are similar to other studies showing the predictive roles of TNFα receptors for incident CKD. In contrast, a study conducted in a Norwegian population showed CRP as the inflammatory biomarker most strongly associated with incident CKD.18 That study, however, was limited to a younger population (middle age) with shorter follow-up time. Additionally, the data were collected from only two follow-up visits, which might have resulted in regression toward mean error.18 These authors also suggested that the previously reported association of TNFα-R2 with incident CKD could be secondary to decline in renal excretion. Our findings, however, are not consistent with this interpretation, as we evaluated the association of the baseline TNFα receptors and the rate of RF decline, independent of baseline RF. Moreover, our results were confirmed in lag-time models where the incident CKD was estimated independent of baseline RF and based only on high inflammation states in the previous visits. These analytical methods should minimize the interference of renal excretion of inflammatory biomarkers.
Of note, most longitudinal studies that reported associations of TNFα-R1 or TNFα-R2 with the rate of RF decline were performed in participants who were significantly older than the Norwegian participants.15–18 While CRP is a marker of general inflammatory status, sTNFα-R2 is a tissue-specific biomarker usually expressed in endothelial cells and leukocytes32 while its expression in other tissues is increased in disease states.33
Some studies have demonstrated that TNF-α and/or its circulating receptors are associated with reduced RF, but these studies have focused mainly on diabetic patients. One such investigation demonstrated that sTNFα receptors strongly predict stage 3 CKD in individuals with type-1 DM without proteinuria.34 There is also evidence that sTNFα receptors predict end-stage renal disease in individuals with type 2 DM.35 Furthermore, previous studies have found an association between high levels of sTNFα-R1 and diabetic nephropathy.34, 36, 37 Our results suggest that, independent of other risk factors, micro-inflammation emerges at a very early stage of RF decline and should be considered as a potential independent risk factor in older adults.15–17, 38 Among the inflammatory biomarkers tested, we found that sTNFα-R1 and sTNFα-R2 were the strongest and most consistent predictors of CKD. Reports indicate that TNFα induces cleavage of the TNFα receptors from the membrane receptor.39 TNFα-R1 activates pathways downstream to NF-κB, such as kinases of mitogen-activated protein kinase, or “death domain molecules,” resulting in IL-6 synthesis or apoptosis.40–42 sTNFα-R1 is considered an indicator of long-term activation of TNFα.34 It has been shown that sTNFα-R1also activates the apoptosis pathway through transforming growth factor beta-1.43 Allugor et al. have reported that the R92Q mutation (rs4149584) increases the affinity of TNF to its receptor. In our study, all individulas were normal homozygotes for this variant.23
Consistent with previous research, we found that increased levels of HS-CRP were associated with the prospective, accelerated decline of RF and CKD development, and increased CRP at the 6-year follow-up predicted incident CKD at the 9-year.18 These associations, however, were attenuated after adjustment for sTNFα-R1.
In our study, we found that sTNFα-R1 and sTNFα-R2 predict incident CKD. We used a categorical threshold of eGFR<60 ml/min/1.73m2 as representing CKD, consistent with consensus definitions commonly utilized in large epidemiologic studies. We further adjusted for baseline eGFR to minimize bias from variation in initial RF.
Although IL-6 and TNFα have been frequently reported in association with age-related outcomes.44–45 The mechanisms that link circulating inflammatory biomarkers and CKD risk, however, remain hypothetical. sTNFα receptors have been shown to play a pivotal role in the development of tubulointerstitial fibrosis and, therefore, nephropathy.46 Moreover, inflammatory-induced atherosclerosis or oxidative stress detected by these biomarkers might mediate RF decline, with the latter linked to increased sTNFα-R2.45 Additionally, age-related chronic inflammation might have an important role in the development of glomerulosclerosis. In a study of kidney donors, Kremers et al. reported that older age was associated with global glomerulosclerosis and higher odds of sclerosis in hypertensive donors.47 Although these investigators did not measure inflammatory biomarkers, it is known that aging is associated with the development of inflammatory states. Mayr et al. demonstrated that angiotensin-II infusion in mice introduces a fibroinflammatory response that is followed by kidney tubulointerstitial collagen deposition mediated through the TNF-R1 signaling pathway and speculating that TNFα-R1 triggers collagen deposition.48 Additionally, the significant association of lipocalin 2 in mesangial cells with sTNFα-R2 levels and impaired kidney function has been reported.38 Therefore, a similar persistent age-related chronic inflammation might result in structural modification of glomerular and interstitial components of the kidneys, with resulting decline in RF.
Strengths and Limitations
Strengths of our study include the longitudinal design, with multiple repeated measures of RF and inflammatory biomarkers over a 9-year follow-up, and use of inverse probability weighting in the analysis to account for random missing values due to missing visits and truncation by death. Furthermore, due to the uncertainty of eGFR thresholds in aging, which identify true disease from age-related decline, we performed multiple analyses using both continuous measures of RF and an epidemiological-based threshold definition for CKD. These analyses demonstrate that the link between sTNFα receptors and RF is robust and independent of baseline RF. Lack of proteinuria measures in follow-up visits is a limitation. Consensus guidelines, however, define eGFR<60ml/min/1.73m2 for more than 3 months as defining stage 3 or greater CKD, irrespective of proteinuria. In addition, albuminuria with preserved RF in older adults is primarily a manifestation of diabetic or hypertensive renal disease. Excluding the individuals with these morbidities in our study, the associations between TNFα receptors and the time-dependent slope and other RF outcomes remained significant.26 One limitation of our study is that the follow-up intervals are 3 years apart. We operationalized time as a continuous variable to overcome this limitation for the annual slope of eGFR; incident CKD, though, is probably underestimated.
Conclusion
sTNFα receptors were the inflammatory biomarkers with the strongest association with the age-related decline of RF both at baseline and over time. Temporal changes in sTNFα receptors were associated with the rate of RF change. Moreover, these biomarkers predicted incident CKD independent of the main conventional risk factors and baseline RF. In contrast, in multivariate analyses, other inflammatory biomarkers, including IL-6, IL-18, and IL-1β, were not consistently significantly associated with RF outcomes. These results suggest that higher levels of sTNFα receptors, in the absence of other chronic conditions such as DM, might be an early sign of impending, accelerated RF decline. Furthermore, these biomarkers could be used to identify individuals at elevated risk of CKD who could be candidates for interventions aimed at preventing progression of RF decline.
Supplementary Material
Table S1. Table S1. Inter- and intra-assay coefficient variation percentage for inflammatory biomarkers at baseline and three follow-ups (FU)ma.
Table S2. Within-individual’s variations of inflammatory biomarkers the 9-year follow-ups.
Table S3. Association of baseline inflammatory biomarker levels measured at baseline with renal function (eGFR) at baseline and each of 3 follow-ups.
Table S4. Association between each baseline inflammatory biomarkers levels and the rate of renal function change (slope of eGFR over time).
Table S5. Association between the rate of each inflammatory biomarkers separately (slope of change of biomarkers over time) and the rate of renal function change (slope of eGFR over time).
Table S6. Association of the inflammatory biomarkers at each visit with incident CKD at subsequent follow-up visits.
Acknowledgments
Dr. Salimi is supported by NIH training grant T32 AG00262. The InCHIANTI study baseline (1998e2000) was supported as a targeted project (ICS110.1/RF97.71) by the Italian Ministry of Health and in part by the National Institute on Aging (NIA) (contracts 263 MD 9164 and 263 MD 821336). Follow-up 1 (2001–2003) was funded by NIA contracts N.1-AG-1– 1 and N.1-AG-1–2111 and follow-ups 2 and 3 (2004–2010) by NIA contract N01-AG-5–0002. Additional support came in part by the Intramural Research Program of the NIA, National Institutes of Health, Baltimore, Maryland.
Footnotes
Conflict of interest
None.
The results presented in this paper have not been published previously in whole or part, except in abstract format.
Author Contributions
All authors contributed to this paper.
Sponsor’s Role
None.
References
- 1.Mills KT, Xu Y, Zhang W, et al. A systematic analysis of worldwide population-based data on the global burden of chronic kidney disease in 2010. Kidney Int. 2015;88:950–957. doi: 10.1038/ki.2015.230. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 2.Fox CS, Larson MG, Leip EP, et al. Predictors of new-onset kidney disease in a community-based population. JAMA. 2004;291:844–850. doi: 10.1001/jama.291.7.844. [DOI] [PubMed] [Google Scholar]
- 3.Fox CS, Gona P, Larson MG, et al. A multi-marker approach to predict incident CKD and microalbuminuria. J Am Soc Nephrol. 2010;21:2143–2149. doi: 10.1681/ASN.2010010085. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 4.Kshirsagar AV, Bang H, Bomback AS, et al. A simple algorithm to predict incident kidney disease. Arch Intern Med. 2008;168:2466–2473. doi: 10.1001/archinte.168.22.2466. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 5.Amdur RL, Feldman HI, Gupta J, et al. Inflammation and Progression of CKD: The CRIC Study. Clin J Am Soc Nephrol. 2016;11:1546–1556. doi: 10.2215/CJN.13121215. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 6.Gupta RK, Kuppusamy T, Patrie JT, et al. Association of plasma des-acyl ghrelin levels with CKD. Clin J Am Soc Nephrol. 2013;8:1098–1105. doi: 10.2215/CJN.09170912. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 7.Epstein M. Aging and the kidney. J Am Soc Nephrol. 1996;7:1106–1122. doi: 10.1681/ASN.V781106. [DOI] [PubMed] [Google Scholar]
- 8.Glassock RJ, Rule AD. The implications of anatomical and functional changes of the aging kidney: with an emphasis on the glomeruli. Kidney Int. 2012;82:270–277. doi: 10.1038/ki.2012.65. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 9.Lindeman RD, Tobin J, Shock NW. Longitudinal studies on the rate of decline in renal function with age. J Am Geriatr Soc. 1985;33:278–285. doi: 10.1111/j.1532-5415.1985.tb07117.x. [DOI] [PubMed] [Google Scholar]
- 10.Yang H, Fogo AB. Cell senescence in the aging kidney. J Am Soc Nephrol. 2010;21:1436–1439. doi: 10.1681/ASN.2010020205. [DOI] [PubMed] [Google Scholar]
- 11.Franceschi C. Inflammaging as a major characteristic of old people: can it be prevented or cured? Nutr Rev. 2007;65:S173–176. doi: 10.1111/j.1753-4887.2007.tb00358.x. [DOI] [PubMed] [Google Scholar]
- 12.Youm YH, Grant RW, McCabe LR, et al. Canonical Nlrp3 inflammasome links systemic low-grade inflammation to functional decline in aging. Cell Metab. 2013;18:519–532. doi: 10.1016/j.cmet.2013.09.010. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 13.Franceschi C, Campisi J. Chronic inflammation (inflammaging) and its potential contribution to age-associated diseases. J Gerontol A Biol Sci Med Sci. 2014;69(Suppl 1):S4–9. doi: 10.1093/gerona/glu057. [DOI] [PubMed] [Google Scholar]
- 14.Michaud M, Balardy L, Moulis G, et al. Proinflammatory cytokines, aging, and age-related diseases. J Am Med Dir Assoc. 2013;14:877–882. doi: 10.1016/j.jamda.2013.05.009. [DOI] [PubMed] [Google Scholar]
- 15.Shankar A, Sun L, Klein BE, et al. Markers of inflammation predict the long-term risk of developing chronic kidney disease: a population-based cohort study. Kidney Int. 2011;80:1231–1238. doi: 10.1038/ki.2011.283. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 16.Carlsson AC, Nordquist L, Larsson TE, et al. Soluble Tumor Necrosis Factor Receptor 1 Is Associated with Glomerular Filtration Rate Progression and Incidence of Chronic Kidney Disease in Two Community-Based Cohorts of Elderly Individuals. Cardiorenal Med. 2015;5:278–288. doi: 10.1159/000435863. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 17.Medenwald D, Girndt M, Loppnow H, et al. Inflammation and renal function after a four-year follow-up in subjects with unimpaired glomerular filtration rate: results from the observational, population-based CARLA cohort. PLoS One. 2014;9:e108427. doi: 10.1371/journal.pone.0108427. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 18.Schei J, Stefansson VT, Eriksen BO, et al. Association of TNF Receptor 2 and CRP with GFR Decline in the General Nondiabetic Population. Clin J Am Soc Nephrol. 2017 doi: 10.2215/CJN.09280916. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 19.Ferrucci L, Bandinelli S, Benvenuti E, et al. Subsystems contributing to the decline in ability to walk: bridging the gap between epidemiology and geriatric practice in the InCHIANTI study. J Am Geriatr Soc. 2000;48:1618–1625. doi: 10.1111/j.1532-5415.2000.tb03873.x. [DOI] [PubMed] [Google Scholar]
- 20.Bandeen-Roche K, Walston JD, Huang Y, et al. Measuring systemic inflammatory regulation in older adults: evidence and utility. Rejuvenation Res. 2009;12:403–410. doi: 10.1089/rej.2009.0883. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 21.Ferrucci L, Corsi A, Lauretani F, et al. The origins of age-related proinflammatory state. Blood. 2005;105:2294–2299. doi: 10.1182/blood-2004-07-2599. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 22.Melzer D, Perry JR, Hernandez D, et al. A genome-wide association study identifies protein quantitative trait loci (pQTLs) PLoS Genet. 2008;4:e1000072. doi: 10.1371/journal.pgen.1000072. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 23.Agullo L, Malhotra S, Fissolo N, et al. Molecular dynamics and intracellular signaling of the TNF-R1 with the R92Q mutation. J Neuroimmunol. 2015;289:12–20. doi: 10.1016/j.jneuroim.2015.10.003. [DOI] [PubMed] [Google Scholar]
- 24.Lauretani F, Semba RD, Bandinelli S, et al. Plasma polyunsaturated fatty acids and the decline of renal function. Clin Chem. 2008;54:475–481. doi: 10.1373/clinchem.2007.095521. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 25.Zamora-Ros R, Rabassa M, Cherubini A, et al. Comparison of 24-h volume and creatinine-corrected total urinary polyphenol as a biomarker of total dietary polyphenols in the Invecchiare InCHIANTI study. Anal Chim Acta. 2011;704:110–115. doi: 10.1016/j.aca.2011.07.035. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 26.Kidney Disease: Improving Global Outcomes (KDIGO) CKD Work Group. KDIGO 2012 Clinical Practice Guideline for the Evaluation and Management of Chronic Kidney Disease. Kidney inter, Suppl. 2013;3:1–150. [Google Scholar]
- 27.Cesari M, Pahor M, Lauretani F, et al. Skeletal muscle and mortality results from the InCHIANTI Study. J Gerontol A Biol Sci Med Sci. 2009;64:377–384. doi: 10.1093/gerona/gln031. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 28.Pizzarelli F, Lauretani F, Bandinelli S, et al. Predictivity of survival according to different equations for estimating renal function in community-dwelling elderly subjects. Nephrol Dial Transplant. 2009;24:1197–1205. doi: 10.1093/ndt/gfn594. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 29.Levey AS, Inker LA, Coresh J. GFR estimation: from physiology to public health. Am J Kidney Dis. 2014;63:820–834. doi: 10.1053/j.ajkd.2013.12.006. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 30.Levey AS, Stevens LA, Schmid CH, et al. A new equation to estimate glomerular filtration rate. Ann Intern Med. 2009;150:604–612. doi: 10.7326/0003-4819-150-9-200905050-00006. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 31.Shardell M, Hicks GE, Ferrucci L. Doubly robust estimation and causal inference in longitudinal studies with dropout and truncation by death. Biostatistics. 2015;16:155–168. doi: 10.1093/biostatistics/kxu032. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 32.Bradley JR. TNF-mediated inflammatory disease. J Pathol. 2008;214:149–160. doi: 10.1002/path.2287. [DOI] [PubMed] [Google Scholar]
- 33.Al-Lamki RS, Wang J, Vandenabeele P, et al. TNFR1- and TNFR2-mediated signaling pathways in human kidney are cell type-specific and differentially contribute to renal injury. FASEB J. 2005;19:1637–1645. doi: 10.1096/fj.05-3841com. [DOI] [PubMed] [Google Scholar]
- 34.Gohda T, Niewczas MA, Ficociello LH, et al. Circulating TNF receptors 1 and 2 predict stage 3 CKD in type 1 diabetes. J Am Soc Nephrol. 2012;23:516–524. doi: 10.1681/ASN.2011060628. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 35.Niewczas MA, Gohda T, Skupien J, et al. Circulating TNF receptors 1 and 2 predict ESRD in type 2 diabetes. J Am Soc Nephrol. 2012;23:507–515. doi: 10.1681/ASN.2011060627. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 36.Navarro-Gonzalez JF, Mora-Fernandez C, Muros de Fuentes M, et al. Inflammatory molecules and pathways in the pathogenesis of diabetic nephropathy. Nat Rev Nephrol. 2011;7:327–340. doi: 10.1038/nrneph.2011.51. [DOI] [PubMed] [Google Scholar]
- 37.Fernandez-Real JM, Vendrell J, Garcia I, et al. Structural damage in diabetic nephropathy is associated with TNF-alpha system activity. Acta Diabetol. 2012;49:301–305. doi: 10.1007/s00592-011-0349-y. [DOI] [PubMed] [Google Scholar]
- 38.Hashikata A, Yamashita A, Suzuki S, et al. The inflammation-lipocalin 2 axis may contribute to the development of chronic kidney disease. Nephrol Dial Transplant. 2014;29:611–618. doi: 10.1093/ndt/gft449. [DOI] [PubMed] [Google Scholar]
- 39.Lantz M, Malik S, Slevin ML, et al. Infusion of tumor necrosis factor (TNF) causes an increase in circulating TNF-binding protein in humans. Cytokine. 1990;2:402–406. doi: 10.1016/1043-4666(90)90048-x. [DOI] [PubMed] [Google Scholar]
- 40.Cabal-Hierro L, Lazo PS. Signal transduction by tumor necrosis factor receptors. Cell Signal. 2012;24:1297–1305. doi: 10.1016/j.cellsig.2012.02.006. [DOI] [PubMed] [Google Scholar]
- 41.Fujita K, Srinivasula SM. Ubiquitination and TNFR1 signaling. Results Probl Cell Differ. 2009;49:87–114. doi: 10.1007/400_2009_18. [DOI] [PubMed] [Google Scholar]
- 42.Rangan G, Wang Y, Harris D. NF-kappaB signalling in chronic kidney disease. Front Biosci (Landmark Ed) 2009;14:3496–3522. doi: 10.2741/3467. [DOI] [PubMed] [Google Scholar]
- 43.Waetzig GH, Rosenstiel P, Arlt A, et al. Soluble tumor necrosis factor (TNF) receptor-1 induces apoptosis via reverse TNF signaling and autocrine transforming growth factor-beta1. FASEB J. 2005;19:91–93. doi: 10.1096/fj.04-2073fje. [DOI] [PubMed] [Google Scholar]
- 44.De Martinis M, Franceschi C, Monti D, et al. Inflammation markers predicting frailty and mortality in the elderly. Exp Mol Pathol. 2006;80:219–227. doi: 10.1016/j.yexmp.2005.11.004. [DOI] [PubMed] [Google Scholar]
- 45.Svenungsson E, Fei GZ, Jensen-Urstad K, et al. TNF-alpha: a link between hypertriglyceridaemia and inflammation in SLE patients with cardiovascular disease. Lupus. 2003;12:454–461. doi: 10.1191/0961203303lu412oa. [DOI] [PubMed] [Google Scholar]
- 46.Guo G, Morrissey J, McCracken R, et al. Role of TNFR1 and TNFR2 receptors in tubulointerstitial fibrosis of obstructive nephropathy. Am J Physiol. 1999;277:F766–772. doi: 10.1152/ajprenal.1999.277.5.F766. [DOI] [PubMed] [Google Scholar]
- 47.Kremers WK, Denic A, Lieske JC, et al. Distinguishing age-related from disease-related glomerulosclerosis on kidney biopsy: the Aging Kidney Anatomy study. Nephrol Dial Transplant. 2015;30:2034–2039. doi: 10.1093/ndt/gfv072. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 48.Mayr M, Duerrschmid C, Medrano G, et al. TNF/Ang-II synergy is obligate for fibroinflammatory pathology, but not for changes in cardiorenal function. Physiol Rep. 2016;4 doi: 10.14814/phy2.12765. [DOI] [PMC free article] [PubMed] [Google Scholar]
Associated Data
This section collects any data citations, data availability statements, or supplementary materials included in this article.
Supplementary Materials
Table S1. Table S1. Inter- and intra-assay coefficient variation percentage for inflammatory biomarkers at baseline and three follow-ups (FU)ma.
Table S2. Within-individual’s variations of inflammatory biomarkers the 9-year follow-ups.
Table S3. Association of baseline inflammatory biomarker levels measured at baseline with renal function (eGFR) at baseline and each of 3 follow-ups.
Table S4. Association between each baseline inflammatory biomarkers levels and the rate of renal function change (slope of eGFR over time).
Table S5. Association between the rate of each inflammatory biomarkers separately (slope of change of biomarkers over time) and the rate of renal function change (slope of eGFR over time).
Table S6. Association of the inflammatory biomarkers at each visit with incident CKD at subsequent follow-up visits.