Introduction

There are marked social inequalities in the occurrence of coronary heart disease (CHD) [1, 2]. The first Whitehall study of British Civil Servants showed that these health inequalities were not confined to persons in poverty [3]. Similarly, in the Whitehall II study of white-collar employees in relatively secure employment, the lower the position in the occupational hierarchy, the higher the risk of CHD [4]. In other European countries [5, 6] and the USA [2], there is also a graded relation between socioeconomic position and CHD incidence [7]. The biological mechanisms linking social position and CHD have not been clear [8]. In the first Whitehall study the combination of blood pressure, plasma total cholesterol, body mass index, plasma glucose and a health behaviour (smoking) accounted for less than one-third of the social gradient in CHD mortality rates [9]. In the British Regional Heart Study a slightly different combination of risk factors accounted for 39% of the social class differences in CHD [10].

A central rationale of the Whitehall II study was to extend the range of potential biological and behavioural explanations. The metabolic syndrome is a promising candidate. It predicts CHD in addition to the predictive power of the established risk factors [11]; moreover, components of the metabolic syndrome and inflammatory markers show a social gradient, evidenced by progressively higher levels as the social hierarchy is descended [12, 13]. Here we wished to test the hypothesis that variables associated with the metabolic syndrome, i.e. inflammation (C-reactive protein [CRP], IL-6, fibrinogen) and height, may mediate the relation between socioeconomic position and CHD. We also sought to determine whether this effect was in addition to the role of adult health behaviours.

Methods

Participants

The Whitehall II cohort of non-industrial civil servants was first examined between 1985 and 1988 (phase 1) [14, 15]. All men and women aged 35 to 55 years and working in the London offices of 20 departments were invited to participate in this study. The cohort consisted of 10,308 persons with an overall response rate of 73%. Postal questionnaires were administered at six further study phases, with phases 3 (1991–1993), 5 (1997–1999) and 7 (2002–2004) including a screening examination. Written informed consent to follow-up sickness absence and clinical records was obtained from participants.

Phase 3 provides the baseline for the analyses reported here, when the metabolic syndrome variables were first measured. The maximum follow-up from phases 3 to 7 was 13.1 years (median 12.2 years). South Asians and African-Caribbeans together make up less than 10% of the cohort. There is insufficient power to analyse them separately. Ethnic group was included as an adjustment variable in all analyses, which yielded similar results when conducted with whites alone. Analyses were confined to men (n = 5312), as too few incident major coronary events occurred among women during follow-up.

Socioeconomic position

Current employment grade within the civil service at phase 3 was used as a measure of socioeconomic position in six categories. For presentation, a three-level variable for grade was used, consisting of clerical and office support staff (low), the three executive officer categories combined (medium) and the two upper administrator categories combined (high).

Behavioural factors

Cigarette smoking was assessed at phase 3 and categorised as: never, ex-smoker, pipe/cigar only or current cigarette smoker with the number of cigarettes smoked. Self-reported leisure-time physical activity at phases 1 and 3 was categorised as vigorous (≥1 h of vigorous activity/week), moderate (≥1 h moderate but <1 h vigorous activity/week) and none/mild (<1 h of vigorous and moderate activity/week). Food frequency questionnaires at phase 3, validated against biomarkers and 7 day diet diaries [16], were cluster-analysed to identify dietary patterns (unhealthy, healthy, Mediterranean-type, sweet) [17]. Alcohol consumption (units/week, 1 unit = 10 ml ethanol) was categorised in four levels (0, 1–14, 15–21 and≥22).

Biological factors

Waist and hip circumferences and BMI were measured at phase 3 using standard protocols [18]. Systolic and diastolic blood pressure were measured twice after 5 min rest using a Hawksley random zero sphygmomanometer (Hawksley and Sons, Lancing, UK). Oral glucose tolerance tests were administered at phase 3, following an overnight fast or in the afternoon after no more than a light fat-free breakfast taken before 08:00 hours. Known diabetic patients did not participate in this part of the screening. Total cholesterol, HDL-cholesterol and triacylglycerol were measured within 72 h in serum stored at 4°C using enzymatic colorimetric methods [12]. Fibrinogen, factor VIIc and von Willebrand factor were determined in citrated plasma, and CRP and IL-6 were determined in serum; samples were stored at −80°C. Plasma glucose was measured using an electrochemical glucose oxidase method and serum insulin by radioimmunoassay using a polyclonal antiserum [12, 13].

The definition of metabolic syndrome used by the National Cholesterol Education Program Adult Treatment Panel III (ATPIII) is based on waist circumference, HDL-cholesterol, triacylglycerol, blood pressure and fasting plasma glucose. We used these variables but instead of using cut-off points to define caseness, we treated the variables as individual continuous variables.

Follow-up of CHD events

About 99.9% of the men were flagged at the NHS Central Registry, who notified us of the date and cause of death. Coronary death was defined by underlying cause (International Classification of Diseases 9, code 410–414) [19]. Potential cases of non-fatal myocardial infarction were ascertained by questionnaire items on chest pain [20], doctors’ diagnoses and hospitalisations. Twelve-lead resting ECGs (Mingorec, Siemens Healthcare, Erlangen, Germany) were performed at each clinic phase and assigned Minnesota codes [21]. Details of physician diagnoses and investigation results were sought from clinical records for all potential cases. Based on all available data from questionnaires, study and hospital ECGs, and cardiac enzymes, definite non-fatal myocardial infarction was defined following MONICA criteria [22]. Myocardial infarction was defined as positive, if a questionnaire or clinical record of diagnosed myocardial infarction was obtained in the presence of an ischaemic ECG; it was defined as negative when self-reported only. Classification of myocardial infarction was carried out independently by two trained coders, who were blind to other study data, with adjudication by a third coder in the rare event of disagreement.

Statistical analysis

Continuous variables were log transformed if skewness was >1.0. Least squares regressions were used to produce age-adjusted least squares means (geometric means for log–normal variables) by employment grade using the age at entry distribution (phase 3) of the analysis sample as the standard population. Differences in means across grades were tested using a linear trend term. Tests for trends in age-adjusted proportions across grades were based on the Mantel–Haenszel test.

Mortality follow-up was available to 30 September 2004 and morbidity follow-up to the phase of most recent participation between phases 3 and 7. Tracing through clinical records identified non-fatal CHD events that occurred after the most recent contact. To allow for this and to overcome differential follow-up for mortality and morbidity rates, non-fatal and fatal follow-up was censored on 30 September 2004 for the 81% (of the 5312 men) who attended phase 7 or on the mid-point date of the phase after the last phase of contact for the remaining 19%, whose last contact had been prior to phase 7. For each coronary outcome, we only considered the time to the first event and excluded prevalent cases of myocardial infarction (n = 86). Event rates were calculated using person-years at risk and standardised for age at entry by the direct method. The associations of employment grade and other risk factors with incident CHD were analysed using Cox models with the time since phase 3 as the underlying time scale. The proportional hazards assumption was tested by fitting exposure by log follow-up time interaction terms and found to be not violated (all p values > 0.25).

The effect of employment grade on CHD incidence was summarised using the relative index of inequality (RII) [23, 24]. The procedure transforms employment grade, a categorical hierarchical variable, into a summary measure scaled from zero (highest socioeconomic position) to one (lowest position), which is weighted to reflect the proportion of the sample in each grade category. We refer to this as the employment grade score. The score is used in a Cox or logistic regression model to estimate the RII: the ratio of instantaneous event rates or log odds between extremes (bottom versus top). The employment grade index was computed within 5-year age strata using six categories of employment grade, since attained grade was associated with age. Using Cox regression, the ability of risk factors to explain grade differences in CHD incidence rates was calculated by computing the proportional (%) change in the RII from models fitted without and with the risk factors of interest. The attenuating effects of traditional risk factors (smoking, blood pressure, serum cholesterol) with respect to the RII for CHD incidence rates were examined, as were those of behavioural factors, metabolic syndrome components and inflammatory markers. Attenuation is calculated as the per cent change in the log hazard ratio or log odds ratio compared with the reference model, controlling for age and ethnic group.

Although data completeness for each individual variable was high, data for one or more variables were missing for approximately one-quarter of participants in the main analysis (Table 3). To retain all participants, multiple imputed values were generated for the missing data from the variables used in the analysis, by means of PROC MI (SAS 8.2, SAS Institute, Cary, NC, USA). Five datasets were randomly selected, and analyses conducted on each of these imputed datasets gave similar results. The mean of these estimates was presented using SAS PROC MIANALYSE. Among the categorical (behavioural) variables, missing values were confined to diet. A complete case analysis of the effect of diet on the social gradient in CHD gave a similar result to that using imputed values.

The impact of measurement error on the contribution of metabolic syndrome variables to the social gradient in CHD was investigated. Within-person variation was estimated using short-term repeat measures obtained at phase 3 (n = 215). This sample size was large enough to ensure that sampling error in the intraclass correlation coefficients estimates would have little impact on the measurement error adjustments. The regression calibration method was used to adjust for measurement error [25]. Since methods for correcting for measurement error in Cox models were not available, logistic regression was used for the measurement error analysis (STATA 8.2, StataCorp, College Station, TX, USA), producing similar results to the Cox models reported in the tables.

Results

In 13.1 years follow-up we documented 201 cases of coronary death or non-fatal myocardial infarction. The overall event rate was 3.3 per 1000 person-years. Figure 1 shows the hazard ratios for CHD across the employment grades, which are well described by a strong linear trend (p < 0.001). Using the relative index of inequality method, the hazard ratio (95% CI) for being of low employment grade was 2.43 (1.48–4.00).

Fig. 1
figure 1

Age-adjusted HR (95% CI) for coronary death/non-fatal myocardial infarction by employment grade, Whitehall II men. The number of CHD events is shown above error bars for each grade. p < 0.001 for linear trend; p = 0.63 for departures from linear trend

Adverse health behaviours and a majority of biological risk factors were inversely associated with employment grade at baseline in 1991 to 1993 (Table 1). Hypotensive medication was related to low employment grade (high 4.9, medium 5.4, low 8.8%, p = 0.02 for trend), whereas use of lipid-lowering drugs was weakly related to high employment grade (high 0.8, medium 0.6, low 0.0%, p = 0.10 for trend). Table 2 shows the age-adjusted relation of risk factors to 13 year incidence of major CHD events. For continuous variables, standardised effects (hazard ratio for 1 SD risk factor difference) are presented. Considered singly, all risk factors except factor VIIc showed the expected association with CHD. The contribution of each variable to the social gradient in CHD incidence was calculated using the RII for employment grade, adjusted for age and ethnicity. Although the majority of this population, even in the lowest grade, were non-smokers, current smoking accounted for 19% of the social gradient in incident CHD. Traditional risk factors (smoking, blood pressure, serum cholesterol) explained, statistically, 30% of the CHD gradient—a finding similar to that observed in the first Whitehall study [3]. Based on a logistic model, additional correction for measurement error increased the explanation by a further 8 percentage points.

Table 1 Prevalence and means of demographic factors, health-related behaviours and risk factors in 5,312 men at phase 3 by employment grade
Table 2 Association between phase 3 risk factors and incident CHD, and the effect of adjustment for single risk factors on the relative index of inequality for grade

Diet and physical activity accounted for 15 and 8% of the social gradient in CHD respectively (6 and 6% after adjustment for smoking). A complete case analysis of the effect of diet on the social gradient in CHD gave estimates and standard errors similar to those derived using imputed values. Components of the metabolic syndrome and inflammatory markers contributed substantially to the social gradient in CHD. BMI, which highly correlated with waist circumference, contributed to the gradient (12%). WHR accounted for more of the gradient (24%) than did waist circumference, but we retained waist, following ATPIII. Height contributed 9% to the social gradient in CHD incidence.

Reductions of the social gradient in CHD incidence by adjustments for metabolic syndrome risk factors in logistic models without (and with) adjustment for measurement error were similar: waist −14% (−14%), systolic blood pressure −8% (−10%), diastolic blood pressure −8% (−12%), HDL-cholesterol −13% (−15%), triacylglycerol −18% (−22%) and fasting glucose 0% (0%). Adjustment for measurement error increased the explanatory effect of fibrinogen (24 versus 15%), CRP (27 versus 22%), IL-6 (29 versus 17%), smoking (25 versus 20%) and diet (24 versus 16%). Adjustment for antihypertensive medication, for which no measurement error estimate was available, explained approximately a further 2% of the RII for employment grade in all models. Adjustment for lipid-lowering medication did not change the RII for employment grade in any of the models.

Health behaviours together explained 30% of the employment grade gradient (Table 3). Metabolic syndrome variables explained 27% (model B). Metabolic syndrome and behavioural variables together reduced the social gradient in CHD by 51% (model C), indicating a degree of overlap of the two sets of explanatory effects. Comparison of attenuating effects in models A, B and C showed that about one quarter of the contribution of metabolic syndrome to the disease gradient arises through social differences in health behaviours. If inflammatory markers are treated as components of the metabolic syndrome cluster, the attenuation of the disease gradient is higher (42%, model B1). Based on models A, B1 and C1, where inflammatory mechanisms are assumed to play a causal role in CHD development, the estimate for the overlap in the explanatory effects of health behaviours and metabolic syndrome on the disease gradient was higher.

Table 3 Effect of adjustment for phase 3 risk factors on the relationship between incident CHD death/non-fatal myocardial infarction and the relative index of inequality for grade

The two groups of variables, behavioural and biological, together with height as a marker of early life influences, both genetic and environmental, increased the explanation of the social gradient in CHD to about 60%. Using WHR rather than waist circumference, attenuation in the full explanatory model including height was the same.

The contribution to the social gradient in CHD of combinations of risk factors was estimated in analyses based on logistic regression models without (and with) adjustment for measurement error. Compared with the corresponding estimates from Cox models, those derived from logistic models were similar for behavioural factors (−31% [−41%]) and metabolic syndrome variables (−27% [−29%]). The contribution to the gradient increased to −44% (−48%) when inflammatory markers were added to the metabolic variables and to −59% (−61%) when health behaviours were also added. In the final model with health behaviours, biological variables and height, allowance for measurement error added 3% (72 versus 69%) to the explanation.

Discussion

Our study set out to determine the extent to which metabolic syndrome variables, inflammatory markers and health behaviours account for the substantial social gradient in incident major coronary events observed in the Whitehall II cohort. We found that, together and without adjustment for measurement error, they account for about half of the social gradient, and indeed more if measurement error is taken into account. The addition of height as a measure of early life influences increased the proportion explained to about 60%, markedly more than that explained by traditional risk factors in the original Whitehall study and consistent with findings from the Kuopio study [26].

Specificity of biological findings

A statistical explanation does not necessarily equate with a biological one. It is nevertheless helpful in testing our prior hypotheses. Earlier, based on cross-sectional analysis, we postulated that von Willebrand factor, a marker of endothelial damage, would be related to the social gradient in CHD [27]. We found little support for that mechanism. On the other hand, there is support for our previously stated hypothesis that social differentials in metabolic syndrome variables provide an important part of the explanation for the social gradient in CHD incidence [12, 18, 28]. In Whitehall II, serum total and LDL-cholesterol, although strongly related to CHD incidence, are not inversely related to social position and are thus not major explanations of the social gradient in CHD.

The nature of the relation between inflammatory markers and CHD remains unresolved [2932]. Since it is not certain that IL-6, CRP and fibrinogen are markers of causal processes or measures of disease status, final analyses are presented with and without inclusion of inflammatory markers. When metabolic syndrome and behavioural variables were included in the model, inflammatory markers added modestly to the explanation of the CHD gradient.

Behavioural factors

Smoking is important as a cause of CHD. However, more than 85% of participants were non-smokers and 81% (75% after adjustment for measurement imprecision) of social differences in disease remained after taking account of smoking. Smoking measured at one point in time may be an inaccurate measure of exposure. We therefore repeated the analyses using accumulated pack-years in current and ex-smokers, based on the first three study phases. The degree of attenuation of the social gradient was unchanged. Both physical activity and dietary pattern [17, 33], but not alcohol consumption, made further modest independent contributions to the explanation, after the effect of smoking had been allowed for. Health behaviours overall accounted for just under a third of the socioeconomic differential in CHD. This exceeds the estimate derived from a representative sample of the US population, which found that health behaviours accounted for 12% of the predictive effect of income on total mortality rates [34]. In the latter study, attenuation of the income effect by behavioural factors was similar for cardiovascular and tumour deaths analysed separately.

Relation between behavioural and biological factors

If indeed these are the relevant biological pathways, it raises the question of how social gradients in the metabolic syndrome variables and in a marker of inflammation are generated. There may be three types of influence: health behaviours, early life and psychosocial factors. Diet, physical activity, smoking and alcohol may all affect metabolic syndrome variables and fibrinogen [12, 28, 35, 36]. Our multivariate analysis indicates that behaviours, excluding alcohol, account for a part of the social gradient in metabolic syndrome. Moreover, the metabolic syndrome variables appear to mediate the effect of social position in part independently of their association with these health behaviours.

Early life

Early life is related to risk of CHD in adulthood [37]. Adult height is inversely related to metabolic syndrome variables [38], plasma fibrinogen [28] and CHD incidence in Whitehall II [39], consistent with an effect of early life on these biological pathways. A further effect of early life is possible, as addition of height to the full model increased the explanation of social differences in CHD.

Psychosocial factors

Psychosocial influences are also possible. Activity of the hypothalamic–pituitary–adrenal axis is related to the metabolic syndrome [40]. We have also shown links between metabolic syndrome variables and both sympathoadrenal and cardiac autonomic activity [41, 42]. Psychosocial stress is related to plasma fibrinogen [28] and prospectively to development of central obesity, metabolic syndrome and CHD [4345]. Among Swedish women [46], a combination of psychosocial and biochemical factors explained much of the social gradient in CHD.

Civil servants and generalisability

The scientific challenge is to understand why CHD follows a social gradient in people who are not poor: the lower the position in the hierarchy the higher the risk [1, 2]. Civil servants are an illuminating population to study. They exclude the richest and poorest members of society, as well as manual occupations and the unemployed. They are men and women with higher than average job security. Despite this, there is a marked social gradient in coronary risk [3, 4]. Similar social gradients have been observed in USA [7, 47], Europe [5, 6] and Australasia [48].

Potential limitations

Non-responders at the phase 3 examination had higher rates of all-cause mortality than those who attended, but the social gradient in mortality rates was similar (p = 0.22). Loss to follow-up may lead to under-identification of non-fatal CHD, but the impact on the grade gradient is likely to be small [49]. Measurement error may lead to underestimation of the contribution of risk factors, particularly of diet, to the social gradient. It had a smaller effect on components of the metabolic syndrome. Furthermore, adjustment for measurement error did not alter our conclusion that health behaviours and biological risk factors make independent contributions to explaining the social gradient in CHD. This observation does not preclude the possibility that explanatory power might be enhanced by multiple measurements of exposure over the life-course. Finally, the analyses here are confined to men; as events accumulate in women, we shall test the degree to which explanations are common across the sexes.

Implications for reducing social inequalities in CHD

Our study provides further evidence that differences in traditional risk factors provide a poor explanation for the socioeconomic gradient in CHD. Two implications flow from this observation. First, it adds to our understanding that cardiovascular risk scores such as QRISK, which include an area measure of deprivation, are better discriminators of risk across the socioeconomic hierarchy than those which do not include such a measure [50]. Second, it highlights the possibility that a prevention strategy focusing on smoking cessation and both blood pressure and cholesterol lowering may not be optimal for reducing social inequalities in CHD.

Conclusion

Our findings highlight the importance of specific biological pathways involving the metabolic syndrome in mediating the social gradient in CHD.