This article requires a subscription to view the full text. If you have a subscription you may use the login form below to view the article. Access to this article can also be purchased.
Abstract
Background and objectives Carotid intima-media thickness has been reported to predict kidney function decline. However, whether carotid intima-media thickness is associated with a hard kidney end point, ESRD, has not been investigated.
Design, setting, participants, & measurements We studied 13,197 Atherosclerosis Risk in Communities participants at visit 1 (1987–1989) without history of cardiovascular disease, including coronary heart disease, stroke, and heart failure, at baseline and assessed whether carotid intima-media thickness measured by B-mode ultrasound is associated with ESRD risk using Cox proportional hazards models. Regarding carotid intima-media thickness parameters, we investigated the mean and maximum values of overall and segment-specific (common, bifurcation, and internal carotid arteries) measurements.
Results Mean age was 54.0 (SD=5.7) years old, and there were 3373 (25.6%) blacks and 7370 (55.8%) women. During a median follow-up of 22.7 years, 433 participants developed ESRD (1.4/1000 person-years). After adjusting for shared risk factors for atherosclerosis and kidney disease, including baseline kidney function, carotid intima-media thickness was significantly associated with ESRD risk (hazard ratio [HR] between quartiles 4 and 1, 1.46; 95% confidence interval [95% CI], 1.02 to 2.08 for overall mean intima-media thickness and HR between quartiles 4 and 1, 1.75; 95% CI, 1.24 to 2.48 for overall maximum intima-media thickness). The associations were largely consistent in demographic and clinical subgroups. When we explored segment–specific intima-media thicknesses, the associations with ESRD were most robust for bifurcation carotid (e.g., adjusted HR between quartiles 4 and 1 of mean intima-media thickness, 1.49; 95% CI, 1.04 to 2.13 for bifurcation; adjusted HR between quartiles 4 and 1 of mean intima-media thickness, 1.36; 95% CI, 0.94 to 1.97 for common; and adjusted HR between quartiles 4 and 1 of mean intima-media thickness, 0.93; 95% CI, 0.67 to 1.29 for internal).
Conclusions Carotid intima-media thickness was independently associated with incident ESRD in the general population, suggesting the shared etiology of atherosclerosis and ESRD.
- arteriosclerosis
- end-stage renal disease
- Epidemiology and outcomes
- Atherosclerosis
- Carotid Intima-Media Thickness
- coronary artery disease
- heart failure
- Humans
- Kidney Failure, Chronic
- Stroke
Introduction
Atherosclerosis is characterized by arterial wall thickening caused by the deposition of plaque consisting of lipoprotein particles, macrophages, and smooth muscle cells (1,2). Representative risk factors for atherosclerosis include age, dyslipidemia, hypertension, cigarette smoking, and hyperglycemia (3–7). Atherosclerosis is a leading cause of coronary heart disease (CHD) and stroke. Although the incidence of atherosclerotic cardiovascular disease (CVD) has been decreasing over several decades, of note, it is still a leading cause of death in the United States (8,9).
Atherosclerosis can also affect kidney function through a few mechanisms. Atherosclerotic renal artery stenosis is a representative example and may lead to kidney dysfunction, mainly through hypertension caused by activated renin-angiotensin system (10). Indeed, as many as 14% of patients with ESRD have ESRD attributable to chronic ischemic nephropathy from renal artery stenosis (11). Moreover, some studies suggest the link between atherosclerosis of intrarenal vasculature and glomerulosclerosis (12). However, it is also possible that atherosclerosis and kidney disease merely share common pathogenic mechanisms or risk factors, such as diabetes and hypertension (13,14).
Elsayed et al. (15) have reported that kidney function declines faster among persons with atherosclerotic disease (CHD, stroke, and peripheral artery disease) than in those without it. However, treatment or clinical examinations in those with atherosclerotic disease (e.g., iodinated contrast agent) may have confounded this association. In this connection, measures of subclinical atherosclerosis, such as carotid intima-media thickness (IMT), would be of value, and indeed, several studies have observed that increased carotid IMT predated kidney disease progression (16–19). However, most of them dealt with selected populations of older adults (16), whites (16,17), or Asians (18), leaving uncertainty regarding their generalizability. Most importantly, none of the previous studies have evaluated ESRD, the most severe form of kidney disease with disproportionately high medical expenditure, as a kidney outcome.
Therefore, the primary objective of this study was to evaluate carotid IMT and its relationship to incident ESRD over 20 years of follow-up in the Atherosclerosis Risk in Communities (ARIC) Study, a large community–based, predominantly biethnic cohort of middle-aged people.
Materials and Methods
Study Population
Details of the ARIC Study have been described elsewhere (20). Briefly, the ARIC Study is a prospective cohort of 15,792 individuals ages 45–64 years old from four United States communities (Forsyth County, NC; Jackson, MS; suburban Minneapolis, MN; and Washington County, MD) at visit 1 (1987–1989). There were four follow-up examinations in 1990–1992 (visit 2), 1993–1995 (visit 3), 1996–1998 (visit 4), and 2011–2013 (visit 5). Of 15,792 participants who attended visit 1, we excluded individuals who were neither white nor black (n=48), had missing data for carotid IMT (n=937), or had prevalent kidney failure (eGFR<15 ml/min per 1.73 m2; n=21) or CVD, including CHD, stroke, or heart failure (n=1589), at baseline, leaving 13,197 individuals for this study. Of those 937 individuals without any carotid IMT data at baseline, 340 did not undergo the ultrasound examination, and 597 had unreadable scans. Possible reasons for having unreadable ultrasound scans include limitations in the ultrasound focal characteristics, ultrasound probe size, arterial depth and geometry, scanning and reading protocols, and sonographer and reader proficiency (21). Reflecting some of these possible reasons, participants without carotid IMT data had higher adiposity compared with those with carotid IMT measures, and thus, they were more likely to be diabetic and have higher BP (Supplemental Table 1). They were also more likely to be black. Prevalent CHD and stroke were defined as self-reported history before visit 1. Prevalent CHD also included evidence of a prior myocardial infarction by electrocardiogram (22). Prevalent heart failure was defined as self-reported treatment or Gothenburg stage 3 at visit 1 (23,24). All individuals provided written informed consent, and procedures were approved by the institutional review board at each study center.
Carotid IMT
Carotid IMT was measured from the blood-intimal to the medial-adventitial interface (25). The ultrasound measurements of the ARIC Study are on the basis of the technique validated by Pignoli et al. (25) using a scanning protocol common to the four field centers (26) and standardized central reading of scans (27). Carotid arteries were examined bilaterally at three sites: the level of the common carotid (1 cm proximal to the dilation of the carotid bulb), the carotid bifurcation (1 cm proximal to the flow divider), and the internal carotid artery (1 cm distal to the flow divider). In each 1-cm section, there were up to 11 measurements at 1-mm intervals (28). Carotid measures were adjusted for site–specific reader differences and downward measurement drifts over the baseline visit (28). We analyzed mean carotid IMT and maximum carotid IMT both overall and by carotid segments as exposures. Overall mean carotid IMT was defined as the mean of all carotid measurements along the far walls (farthest from the skin surface) at each of the six sites, whereas segment–specific mean carotid IMT was defined for each segment as the mean of far walls of both the left and right sides. For overall mean IMT, any missing sites were imputed from sex– and race–specific multivariate linear models of mean IMT as a function of age, body mass index (BMI), and arterial depth fitted by maximal likelihood methods as previously conducted in the ARIC Study (28,29). Overall maximum carotid IMT was defined as the maximum value from all carotid measured parameters of the far walls at the six sites, whereas segment–specific maximum carotid IMT was defined for each segment as the maximum far wall measurement from the left and right sides. For segment–specific maximum IMT, we rely on measured values without imputation, and thus, data are available in 13,106 individuals for common carotid IMT, 10,643 individuals for bifurcation IMT, and 8102 individuals for internal IMT.
Covariates
Information on demographics (age, sex, and race), lifestyle, and medical history was collected at visit 1 by trained interviewers using standardized questionnaires. Smoking status and alcohol intake were determined by self-report. BMI was calculated as weight in kilograms divided by the square of height in meters. Certified technicians used a random zero sphygmomanometer to measure BP three times with participants in the sitting position after 5 minutes of rest, and the average of the last two readings was recorded. Hypertension was defined as systolic BP ≥140 mmHg, diastolic BP ≥90 mmHg, or antihypertensive medication use. Diabetes mellitus was defined as a fasting glucose level ≥126 mg/dl, a nonfasting glucose level ≥200 mg/dl, self–reported physician diagnosis of diabetes, or use of glucose-lowering medications. Medication use was verified by inspection of medication bottles. Total cholesterol and HDL cholesterol were determined by enzymatic methods. eGFR was calculated using the Chronic Kidney Disease Epidemiology Collaboration creatinine equation (30).
Outcome Assessment
Incident ESRD was defined as initiation of dialysis therapy, transplantation, or death caused by kidney disease (31). Patients with dialysis therapy and transplantation were identified by linkage to the US Renal Data System (USRDS), which captures information about all Americans with ESRD who receive RRT or are awaiting kidney transplantation (32). Death caused by kidney disease was defined as death with any of the following the 9th International Classification of Diseases (ICD-9) codes (as well as corresponding ICD-10 codes) in the primary code of the death certificate: kidney transplant (V42.0, 55.6, or 996.81), dialysis (V45.1, V56, 39.95, or 54.98), CKD (403, 404, 582, 583, 585, 585.3–585.5, 586, 587, 588, or 593.9), or diabetes with renal complications (250.4). Participants who were free of ESRD by December 31, 2011 were administratively censored.
Statistical Analyses
We first evaluated the Pearson correlation coefficients between carotid IMT and baseline eGFR. Carotid IMT parameters were primarily modeled using their quartiles (28,33). Therefore, baseline characteristics were summarized according to quartiles of overall mean carotid IMT, and differences were assessed by chi-squared test and ANOVA as appropriate. To visualize the potentially nonlinear association between overall mean and maximum carotid IMTs and ESRD, incidence rate of ESRD adjusted for age, sex, race, and center was evaluated using a Poisson regression model with linear spline terms of carotid IMT (three knots corresponding to cutoffs of its quartiles). We then quantified the association of each carotid IMT measure with incident ESRD using Cox proportional hazards models. We used two models with progressive adjustment. Model 1 adjusted for age, sex, race, center, and baseline eGFR. Model 2 also included systolic BP, antihypertensive medication, smoking, alcohol intake, BMI, total and HDL cholesterol, diabetes, and statin use. P value for linear trend was obtained by assigning the median value of each quartile of carotid IMT parameters to individuals in that quartile.
We conducted several sensitivity analyses. First, we examined the association in subgroups defined by age (<55 versus ≥55 years old), sex, race, diabetes, hypertension, obesity (BMI<30 versus ≥30 kg/m2), and smoking. Because of low numbers of patients with ESRD in some of the subgroups, we estimated hazard ratios (HRs) comparing the highest quartile with the lower three quartiles for these subgroup analyses. The likelihood ratio test was used to test interactions. Second, we conducted competing risk analysis with death as a competing end point of ESRD. Third, we included incident CVD (CHD, stroke, and heart failure) events during follow-up as time-varying covariates. Incident CHD and stroke were defined as adjudicated patients after visit 1. Incident heart failure was defined as hospitalization for heart failure after visit 1 (23,24). All statistical analyses and graphical displays were performed using Stata, version 13.0 (StataCorp., College Station, TX).
Results
Mean age was 54.0 years old (SD=5.7), and there were 3373 (25.6%) blacks and 7370 (55.8%) women. Compared with those in the lowest quartile, participants with higher overall mean carotid IMT level were more likely to be older, black, smokers, and diabetic and have higher BMI and BP, abnormal lipid profile, and lower eGFR (Table 1). Overall mean and maximum IMT levels were correlated mildly with baseline eGFR (Pearson correlation coefficients =−0.14 and −0.09, respectively).
Baseline characteristics by quartiles of overall mean intima-media thickness
During a median follow-up of 22.7 years, there were 433 patients with ESRD (crude incidence rate =1.4/1000 person-years). After adjusting for age, sex, race, center, and baseline eGFR, overall, there was a graded association of overall mean IMT with incident ESRD, with risk gradient of two- to fourfold between the range in top and bottom quartiles (Figure 1A). Similar patterns were observed for overall maximum carotid IMT (Figure 1B). Additional adjustment for shared risk factors for atherosclerosis and kidney disease, such as BP, lipids, and diabetes, attenuated the association, but quartile 4 showed significantly higher risk of ESRD compared with quartile 1 (HR, 1.46; 95% confidence interval, 1.02 to 2.08 for overall mean IMT and HR, 1.75; 95% confidence interval, 1.24 to 2.48 for overall maximum IMT [model 2 in Table 2]). When restricting to patients with ESRD determined by USRDS linkage, we found similar results (Supplemental Table 2). The results were largely consistent in demographic and clinical subgroups that we tested, although we observed a significant interaction between races for overall maximum IMT (P=0.02), with a higher relative risk of ESRD according to high IMT in whites than in blacks (Figure 2).
Adjusted incidence rates of ESRD by overall mean carotid intima-media thickness (IMT) and overall maximum carotid IMT. Incidence rate per 1000 person-years was adjusted for age, sex, race, center, and eGFR and trimmed at 0.5% and 99.5%. Gray regions indicate the 95% confidence intervals of adjusted incidence rates. (A) indicates overall mean carotid IMT; (B) indicates overall maximum carotid IMT.
Hazard ratios (95% confidence intervals) of incident ESRD by quartiles of overall intima-media thickness
Adjusted hazard ratio (HR) of ESRD comparing the highest quartile with the lower three quartiles in subgroups. Adjustment was on the basis of model 2 (age, sex, race, center, systolic BP, antihypertensive medication, smoking status, alcohol intake, body mass index (BMI), total and HDL cholesterol levels, diabetes, statin use, and eGFR). 95% CI, 95% confidence interval.
In competing risk analysis, an attenuated association was observed for overall mean IMT, but there was still a significant association for overall maximum IMT (Table 3). When incident CVD events were accounted for as time-varying covariates, the associations were attenuated and became nonsignificant for overall mean IMT, whereas they remained mostly significant for overall maximum IMT (Table 3). Among all CVD subtypes assessed, the adjustment for heart failure attenuated the associations the most.
Hazard ratios (95% confidence intervals) of incident ESRD by top versus bottom quartiles of intima-media thickness on the basis of competing risk model and cardiovascular disease as time-varying covariates
When we explored segment-specific IMTs, the associations with ESRD were most robust for bifurcation IMT, regardless of whether its mean or maximum value was used (Table 4). The mean and maximum common carotid IMTs showed a more evident dose-response relationship than bifurcation IMT in model 1, but the associations were considerably attenuated after accounting for potential confounders, particularly when its mean value was analyzed. The associations of internal IMT and ESRD risk were overall weak and no longer significant in model 2.
Hazard ratios (95% confidence intervals) of incident ESRD by quartiles of intima-media thickness for different carotid segments
Discussion
We found that greater overall carotid IMT levels were associated with incident ESRD in the general population independent of shared risk factors for atherosclerosis and kidney disease. The elevated risk of ESRD was particularly evident in the top quartile of carotid IMT. Although mean and maximum carotid IMTs showed similar results, overall, the associations tended to be stronger for maximum IMT than for mean IMT in our study. These associations were attenuated after controlling for incident CVD events during follow-up, particularly when we accounted for incident heart failure. Although the association between maximum carotid IMT and ESRD risk was greater in whites than in blacks, overall, qualitatively consistent results were seen in several demographic and clinical subgroups. Regarding specific carotid segments, we found that IMT at bifurcation seemed to be most robustly associated with ESRD risk.
Our results are largely consistent with previous studies reporting the associations of increased carotid IMT levels with incident reduced kidney function and kidney function decline (16–18). We extended the results from previous studies by confirming the association with ESRD, a hard kidney outcome with great clinical and social importance. Also, it was important to confirm this association in blacks. Moreover, we comprehensively assessed both mean and maximum values of carotid IMT and systematically investigated segment-specific IMT. We also tried to apply sophisticated statistical approaches, such as spline modeling, competing risk, and time–varying covariate analysis for incident CVD events.
There are several potential mechanisms linking atherosclerosis to elevated ESRD risk. Atherosclerosis and ESRD are known to share risk factors, such as hypertension and diabetes (13,14). However, the associations of carotid IMT with ESRD remained significant after the adjustment for those factors, suggesting the independent contribution of subclinical atherosclerosis to the development of ESRD. This independent association may be because of the property of IMT reflecting the past profile of those risk factors over time. In this context, renal artery stenosis is a representative atherosclerotic disease linked to kidney dysfunction and may damage the kidney through hypertension because of the activated renin-angiotensin system (10,11). Moreover, Kasiske (12) found a significant and independent relationship between intrarenal arterial wall thickness and a number of sclerotic glomeruli on the basis of autopsy data.
Despite these potential mechanisms linking atherosclerosis to kidney disease progression, randomized controlled trials have yielded inconclusive results regarding whether statin can prevent kidney disease progression (34–36). However, we need to keep in mind that most previous statin trials were designed to assess the prevention of cardiovascular outcomes and underpowered to evaluate renoprotection (34–36). Given that the significant association with ESRD risk was observed primarily in the top quartile of carotid IMT, it may be that statin therapy is particularly effective in the prevention of kidney disease progression among people with higher burdens of subclinical atherosclerosis.
Although similar associations were observed in the demographically adjusted model, overall maximum IMT showed a stronger association with ESRD than overall mean when further adjusting for shared risk factors. Overall maximum IMT also displayed a more robust association with ESRD when accounting for incident CVD events as time-varying covariates. None of the previous studies investigating the association of carotid IMT with kidney dysfunction have compared mean and maximum IMTs in this regard (16–18). However, our findings seem consistent with the observations that maximum carotid IMT shows stronger associations with CVD events compared with mean IMT in some studies (37–39). Although this finding warrants confirmation in other studies, it is possible that atherosclerosis does not occur equally in every vascular bed, and maximum IMT may represent the pathologic atherosclerotic process better than mean IMT. It is also possible that maximum IMT reflects the presence of plaque more than mean IMT and that plaque may be important in predicting ESRD risk.
Among all segments assessed, the associations of IMT levels with ESRD were most robust for bifurcation carotid and least robust for internal carotid. This may be partially because of the precision of measurement in each site. Indeed, in the ARIC Study, intraclass correlation coefficients between repeated assessments by different sonographers and different readers were 0.77 for bifurcation, 0.73 for common, and 0.70 for internal far wall IMT (40). Also, it is possible that bifurcation reflects atherosclerotic burden better than other segments. Indeed, the physiologic predisposition of atherosclerosis in carotid bifurcation is well known, possibly because of geometry, velocity profile, and shear stress (41,42). Low shear stress, which is associated with endothelial damage, preferentially affects the outer edges of vessel bifurcations (41). Findings from flow analysis and carotid endarterectomy showed that the greatest plaque thickness occurs in the outer wall of carotid bifurcation (42). Autopsy studies have also shown that the most frequent site for atherosclerosis is carotid bifurcation (43,44).
Among CVD subtypes that we accounted for as a time-varying covariate, incident heart failure attenuated the associations of carotid IMT with ESRD the most. This may suggest that heart failure is an important mediator between subclinical atherosclerosis and ESRD risk. Indeed, ischemic heart disease is a major cause of heart failure, and patients with heart failure are known to be at high ESRD risk because of decreased kidney perfusion from low cardiac output and diuretic therapy and excessive activation of the renin-angiotensin-aldosterone system (45–47). Nevertheless, this finding may highlight the importance of following up kidney function among those with atherosclerosis, particularly when they develop heart failure in the clinical course.
This study has several limitations. First, because of data unavailability at visit 1 in the ARIC Study, we were not able to adjust for albumin-to-creatinine ratio, a potent risk factor for ESRD other than eGFR (48). Second, there were a limited number of ESRD events in participants of certain quartiles of carotid IMT, and thus, we may not have enough power in certain strata in subgroup analysis. Third, as an observational study, this study may not be exempt from residual confounding (for example, related to inflammation and oxidative stress) (49,50). Fourth, as mentioned above, approximately 5% of the ARIC Study participants (937 of 15,792) did not have carotid IMT data usable for this study and had unique characteristics compared with the final study population (e.g., higher likelihood of obesity and blacks). This may somewhat question generalizability of our findings, but it is important that we confirmed consistent results, regardless of these clinical and demographic factors.
In summary, greater carotid IMT levels were independently associated with incident ESRD in a biethnic middle–aged general population. The associations were largely consistent in several demographic and clinical subgroups. These findings suggest the shared etiology of atherosclerosis and ESRD.
Disclosures
None.
Acknowledgments
The authors thank the staff and participants of the Atherosclerosis Risk in Communities (ARIC) Study for their important contributions.
The ARIC Study is carried out as a collaborative study supported by National Heart, Lung, and Blood Institute contracts HHSN268201100005C, HHSN268201100006C, HHSN268201100007C, HHSN268201100008C, HHSN268201100009C, HHSN268201100010C, HHSN268201100011C, and HHSN268201100012C.
Some of the data reported here have been supplied by the US Renal Data System. The interpretation and reporting of these data are the responsibility of the authors and in no way should be seen as an official policy or interpretation of the US Government.
Footnotes
Published online ahead of print. Publication date available at www.cjasn.org.
This article contains supplemental material online at http://cjasn.asnjournals.org/lookup/suppl/doi:10.2215/CJN.11951115/-/DCSupplemental.
- Received November 12, 2015.
- Accepted March 17, 2016.
- Copyright © 2016 by the American Society of Nephrology
References
In this issue
