Introduction

Over the last decades, advances in diagnostic imaging techniques have been impressive. At the same time, this success has resulted in an increase in unrequested information during the diagnostic work-up of patients. Unrequested information can be defined as imaging abnormalities of potential clinical relevance that are unexpectedly discovered and unrelated to the clinical indication. Although unrequested information occurs across different imaging techniques and anatomical regions, computed tomography (CT) of the chest provides numerous exemplary cases of unrequested information: vascular calcifications, pulmonary emphysema, pulmonary nodules, or skeletal osteoporosis. Due to lack of follow-up, the clinical importance of unexpected abnormalities is largely unknown. Performing follow-up studies is the only way to determine which imaging abnormalities have or do not have clinical relevance. Alternatively, unrequested information could also contain valuable prognostic information on a variety of prevalent diseases in an ageing western population providing a new means of identifying patients that could profit from primary prevention efforts.

Non-invasive assessment of coronary artery calcium (CAC) by various CT techniques has been widely used as a prognostic marker of cardiovascular disease (CVD) in large prospective screening studies [15]. Vascular calcifications have a high prevalence in the normal population [6] and will thus be detected frequently when reading diagnostic CT in a clinical care population. To date, no follow-up studies have been performed to investigate whether vascular calcifications—as a prevalent example of unrequested information—can be used as prognostic marker in a population receiving routine clinical care.

In this study we examined whether subclinical coronary and/or aortic calcification detected unexpectedly by routine diagnostic chest CT is associated with future CVD events.

Methods

Population

All subjects are participants of the PROVIDI study, a multicentre retrospective cohort study among persons 40 years of age and older who underwent chest CT in one of 8 participating hospitals between January 2002 and December 2005. In patients with more than one chest CT examination during follow-up, only the first CT scan was used for analysis. The original cohort consisted of 23,443 patients. Based on CT reports, a research physician excluded all patients (n = 9077) with a poor prognosis defined as (1) diagnosis of primary lung cancer (including mesothelioma), or (2) diagnosis of distant metastatic disease from other types of cancer (not including haematological malignancies). This exclusion criterion was applied as it is highly unlikely that in patients with a poor prognosis detection of subclinical disease markers will alter clinical decision making. Similarly, subjects with suspected cardiovascular disease or patients known to have symptomatic cardiovascular disease were excluded (n = 2303). Starting May 2008, CT data of all eligible patients were transferred to the primary study site to be evaluated for the presence of unrequested information.

The institutional review boards of all participating centres approved this study, and written informed consent was waived because of the retrospective design of this study; a privacy protocol was implemented.

After exclusion, a total of 12,063 patients were eligible for this study. All 1653 patients from one randomly chosen study centre were additionally excluded to be used as an external validation sample for future research. Chest CT examinations were thus available for 10,410 patients. A flow-chart detailing the inclusion and exclusion steps for this research project is provided as Appendix to this study.

Study design

We used a case-cohort design in which a representative subcohort (n = 1285) is randomly sampled from the baseline cohort (n = 10,410) at the start of the study. This random sample should be >10% of the size of the baseline cohort [7]. All participants in the baseline cohort are followed-up for the outcome-of-interest. Incident cases plus the random subcohort define the actual population under study. In the data analysis, all members of the subcohort are weighted by the inverse sampling fraction to reach valid estimates for the full cohort. This design provides an efficient way to avoid scoring CT characteristics for all 10,410 subjects of the baseline cohort, but instead only in incident cases and a representative baseline sample.

Chest CT protocols

To mimic day-to-day clinical practice and to increase the generalisability of our results, we chose to include CT from a wide spectrum of different CT machines to imaging protocols. All CT were obtained with dual, 4-, 8-, 16-, 40-, or 64-slice systems from different vendors. The field-of-view (FOV) of all eligible chest CT protocols had to include the heart and the full length of the thoracic aorta. Slice thickness varied according to the CT indication and corresponding protocol: e.g. 3 mm for CT pulmonary embolism, 5 mm for standard chest CT and 10 mm for high-resolution CT of the lung parenchyma.—CT examinations with and without intravenous contrast enhancement were used. One research physician abstracted and classified the CT indication on the basis of information from CT reports; examples include suspected pulmonary disease, cardiovascular disease, traumatic injury, or mediastinal disease.

Scoring of CT characteristics

Computed tomography data sets were scored with regard to coronary and aortic calcification by three readers: one board certified radiologist with 10 years of experience and two research physicians with 2 and 3 years of experience in reading chest CT. All readers were blinded for patient characteristics and outcome status. CT images were stored in DICOM format and read using a commercially available software program (Dicomworks, version 1.3.5). All CT were analysed in the axial plane; no other reformatted planes or thin slice reconstructions were included for analysis. Readings were performed using a standardised score form including information on the type of CT protocol (section thickness, tube load (mAs) and voltage (kVp), and the use of contrast agent) and image quality (defined as good, adequate or poor).

The extent of vascular calcification was scored using simple visual grading (0–3) as previously described [8]. Table 1 provides exact definitions. Briefly, the four main coronary artery branches were scored separately (0–3) before being summed into a single score for the coronary artery tree (0–12). Calcification of the thoracic aorta was scored in three locations: calcification in the ascending (0–3) and descending (0–3) aorta and calcification in the supra-aortic arteries branching of the aortic arch (0–2), which were summed into a single score for the thoracic aorta (0–8). The lower cut-off level for the descending aorta was defined as the mid-level of the 11th thoracic vertebra. Figure 1 provides a visual example of different grades of CAC.

Table 1 Definitions used for visual grading of calcified plaques on routine diagnostic chest CT
Fig. 1
figure 1

Axial 16-slice CT images of different grades of coronary artery calcification in the left anterior descending (LAD) branch (arrows demonstrate the example area): a grade 1 (2 focal calcifications); b grade 2 (1 large calcification); c grade 3 (extensive calcification covering a large segment of the LAD branch)

Use of the score form was trained under the supervision of an experienced board certified chest radiologist using a training set of 50 randomly selected patients. Reproducibility of visual grading was evaluated. Briefly, weighted kappa for inter -and intra-observer variability of aortic wall abnormalities (including calcifications) were 0.72 and 0.88, for calcifications of the supra-aortic arteries 0.89 and 0.96, and for calcifications in the coronary arteries 0.77 and 0.91, respectively.

Follow-up

We recorded incident fatal and non-fatal CVD events for a mean of 17.8 months. End-point status was obtained through linkage of patients with the National Death Registry and the National Registry of Hospital Discharge Diagnoses from January 2002 to December 2005. Database linkage was performed with a validated probabilistic method [912]. In these databases, cause of death and the occasion of hospitalisation are coded according to the International Classification of Diseases, 9th and 10th revision (ICD-9, ICD-10) [13]. Correct designation of causes of death has been established with in a comparison study with patient medical records [14]. The incidence of CVD events in the full cohort of 10,410 participants was n = 515.

A trained research physician abstracted all valid end points from these databases on the basis of ICD codes. Fatal and non-fatal CVD events were defined by ICD-9 codes as coronary heart disease (CHD) (codes 410–414), heart failure (code 428), peripheral arterial disease (PAOD) (codes 440, 443–444), aortic aneurysm (code 441), cerebrovascular disease (codes 430–438), or non-rheumatic valvular disease (code 424). All other codes included in the ICD-9 paragraph on diseases of the circulatory system were not considered as valid end points for this study. In the case of multiple valid end points in the same patient, cause of death prevailed over hospital admissions or else the first hospital discharge diagnosis was used. Secondary analyses were performed for a selection of cases from the primary end point: all fatal and non-fatal coronary events (codes 410–414) and all fatal and non-fatal non-cardiac events (PAOD, codes 440, 443–444; aortic aneurysm, code 441; and cerebrovascular events, codes 430–438).

Data analysis

Differences in risk factors and CT protocols between subjects from the subcohort to various case-groups were tested with Chi-squared tests (discrete variables) and one-way ANOVA (continuous variables). Next, we calculated annualised event rates for all three end points stratified according to coronary and aortic calcium scores.

We used Cox proportional-hazards regression to estimate hazard ratios for any CVD event, coronary events, and any non-cardiac event according to calcium scores. Sum scores for visually graded coronary and aortic calcifications were analysed as categorical and as continuous variables (per SD). Categories of the aortic calcification sum score (0–8) were chosen as: no calcification (reference category), mild calcification (sum score 1–2), moderate calcification (sum score 3–4) and severe calcification (sum score 5 and higher). In the case of coronary calcifications (0–12), category cut-offs were chosen as: no calcification (reference category), mild calcification (sum score 1–2), moderate calcification (sum score 3–5) and severe calcification (sum score 5 and higher).

After estimating the hazard ratios for all non-cardiac events according to the aortic sum score, we performed additional regression analyses for all cases of stroke, looking at calcifications in each anatomically distinct part of the thoracic aorta (ascending, descending, arch). Finally, we investigated the multivariate associations when coronary and aortic sum scores were modelled simultaneously.

Apart from crude associations, all models were adjusted for age, sex, indication for CT, image quality to type of medical centre in which the CT was performed (tertiary/secondary). We also tested the interaction of aortic and coronary sum scores with section thickness (continuous) and the use of contrast agent (yes/no) by entering interaction terms into the model. We had missing values for CT characteristics in <3%. We used regression methods implemented in SPSS software (SPSS 14.0, Chicago, IL, USA) to impute missing values. All other analyses were performed with R software, version 6.2.

Results

A total of 116 subjects (39 of 515 cases (7.6%), and 77 of 1285 subjects from the subcohort (6.0%)) were excluded from analyses because CT images could not be retrieved from CT databases. A total of 515 subjects had a CVD event during follow-up: fatal CVD events occurred in 240 patients, and 275 patients had a non-fatal event. Table 2 summarises the baseline characteristics of the subcohort and the various case-groups.

Table 2 Baseline characteristics

Table 3 shows a gradual increase in annualised event rate for any CVD event with increasing CAC and TAC score categories.

Table 3 Annualised event rates for any CVD event (n = 515) according to coronary artery calcium score and thoracic aorta calcium score categories

Subjects with coronary and non-cardiac events were older than people in the subcohort, but only subjects with coronary events were more often male.

The prevalence of CAC (sum score >0) was 67% for people in the subcohort compared with 88% in subjects with CVD events, 90% for coronary events, and 85% for non-cardiac events (P < 0.001 for all three groups). The corresponding figures for TAC >0 were 61% (subcohort), 82%, 80% and 86%, respectively (P < 0.001 for all three groups).

Table 4 shows that the risk of any CVD event increases with an increase in category of TAC and CAC scores after adjustment for age, sex, CT indication, CT image quality and type of medical centre. A statistically significant increase in hazard ratio occurs already in the case of mild calcification (TAC and CAC scores = 1–2). The risk of any CVD event was increased by a factor of 2.7 (95% CI, 2.0–3.7) and 3.7 (95% CI, 2.7–5.2), respectively, among patients with severe TAC or CAC (TAC score ≥5; CAC score ≥6) compared with those without any calcium. In order to compare the strength of the associations between the two calcium measures, Table 4 also shows the increase in risk of any CVD event associated with CAC and TAC as continuous measures (per 1 SD increase in the TAC and CAC sum scores). After adjustment, an increase of 1 SD in TAC and CAC sum scores resulted in a 46% and 41% increase in the risk of any CVD event.

Table 4 Risk of any CVD event, coronary events, and non-cardiac events associated with increasing coronary artery calcium score and thoracic aorta calcium score in clinical care patients undergoing routine diagnostic chest CT

Looking at coronary and non-cardiac events separately (Table 4, lower half), we could demonstrate a more differentiated pattern in the associations with TAC and CAC scores. Risk of coronary events increased 5-fold among subjects with severe CAC (CAC score ≥6) compared with those without any coronary calcium, whereas having severe TAC (TAC score ≥5) resulted in an increase in risk by a factor of 2.4. Conversely, risk of non-cardiac events (stroke, aortic aneurysm/dissection, PAOD) increased by a factor of 3.3 for subjects with severe TAC compared with patients without any aortic calcium, whereas people with severe CAC showed an increase in risk by a factor of 2.3. No statistically significant interaction was found between section thickness and TAC (P = 0.46) or CAC (P = 0.89) and the use of intravenous contrast and TAC (P = 0.40) or CAC (P = 0.81) in relation to the risk of CVD events.

Table 5 shows the risk of stroke associated with increasing amounts of calcium in each of the three anatomically distinct parts of the thoracic aorta. Remarkably, strongest associations (and the only to reach statistical significance) were demonstrated for calcifications in the ascending aorta. Moderate calcification of the ascending aorta increased the risk of stroke by a factor 3.0 compared with no calcification, whereas moderate calcification of the descending aorta or supra-aortic branches increased the risk by a factor of 1.8 and 1.4, respectively.

Table 5 Risk of stroke (n = 58) associated with calcium in three anatomically distinct parts of the thoracic aorta

We also calculated multivariate-adjusted models with TAC and CAC sum scores entered simultaneously. Both scores were independently associated with the risk of (1) any CVD event increasing the risk by 25% (95% CI, 13–38%) per 1 SD increase in CAC score and 30% (95% CI, 17–44%) per 1 SD increase in TAC score, respectively; and (2) coronary events increasing the risk by 42% (95% CI, 25–61%) per 1 SD increase in CAC score and 19% (95% CI, 4–36%) per 1 SD increase in TAC score, respectively. Only TAC was independently associated with the risk of non-cardiac events (45% increased risk per 1 SD increase in TAC score).

Discussion

We examined whether unrequested information from routine diagnostic chest CT in a clinical care population could be used as predictor of future disease. Severe coronary (CAC) or aortic (TAC) calcification—based on visual grading of calcified plaques—increased the risk of any CVD event almost 4-fold and 3-fold, respectively, after a mean follow-up of close to one and a half years. CAC was found to be a stronger predictor of coronary events, whereas TAC was more strongly associated with non-cardiac events.

More than 60 million diagnostic CT are performed annually in the USA and this number will likely increase as the population continues to age [15]. These diagnostic images contain unrequested, subclinical findings that could contain important prognostic information. Obtaining this information is ‘free’, as it is contained in routine care and comes at no additional exposure to ionising radiation. Currently, no data are available on the natural course and prognosis of these findings. The prognostic value (or lack thereof) of unrequested subclinical findings is important to indicate new ways of extracting clinically relevant information. Results from this study suggest that CT images obtained as part of routine clinical care do indeed convey important ancillary information.

We investigated the prognostic value of coronary and aortic calcifications as these are highly prevalent [6] and long recognised as prognostic markers of CVD in asymptomatic screening populations [15]. The prevalence of any CAC in our study (67%) is comparable with these previous studies specifically aimed at screening for CAC and showing ranges between 53 and 69% [2, 16, 17]. Like these screening studies we show that CAC is a strong predictor of future CVD events. However, to our knowledge, our results are the first to demonstrate this association in a routine clinical care population. Furthermore, we used a simple visual grading to score calcifications on a heterogeneous set of CT protocols. Existing calcium scoring typically involves non-contrast-enhanced, ECG-gated CT protocols with semi-automatic calcium scoring software. The scoring method we used consists of simple definitions, is reproducible, can be performed with a minimum of additional reviewing time and can be applied to a whole range of different CT protocols (including non-gated CT). Prognostic information on CVD risk can be appended to the CT report. Together with all relevant clinical data available from a patient’s medical record, the clinician can decide whether additional work-up is required. This will result in a more consistent and clinically useful way of reporting unrequested subclinical findings.

Cardiovascular disease is one of the leading causes of death and hospitalisations in western society [18]. Major advances in knowledge about CVD risk factors have been made in recent decades and international guidelines for the primary prevention of CVD have been actively developed [1921]. However, for a variety of reasons substantial groups of people at intermediate or high risk of CVD events according to existing guidelines are currently not recognised and treated as such [22, 23]. This has prompted expert panels to state that the search for new strategies to detect patients who would benefit most from intensive primary prevention efforts is a clinically important objective [24, 25]. The present study is an attempt at defining such a new strategy.

Unfortunately, no information from patients’ medical records was available for this study. This might have influenced our results in two ways. First, exclusion of patients with symptomatic CVD could now only be performed based on information from CT reports, possibly resulting in a number of symptomatic CVD patients being missed. This would have caused an overestimation of our risk ratios. However, we excluded approximately 16% of patients from the baseline clinical care cohort based on information from CT reports (2,303/14,366) which is slightly more than the prevalence of symptomatic CVD (11.2%) in the general white population of the USA [26]. Second, lack of information from medical records prevented us from adjusting our risk estimates for traditional CVD risk factors (hypertension, hyperlipidaemia, diabetes, smoking). Based on previous reports we expect that only limited attenuation of our age- and sex-adjusted risk estimates will occur when adjusting additionally for other CVD risk factors [3].

Conclusion

Unrequested subclinical findings—such as coronary and aortic calcifications—detected on routine care, diagnostic chest CT in a clinical care population can be used as a prognostic marker of future CVD events. This may be achieved with a simple visual scoring method that is easy to use in daily clinical practice. This offers a novel approach in trying to employ potentially valuable prognostic information contained in routine diagnostic imaging and may be used to identify patients who could benefit from primary preventive efforts.