Skip to main content

Proteomic risk markers for coronary heart disease and stroke: validation and mediation of randomized trial hormone therapy effects on these diseases

Abstract

Background

We previously reported mass spectrometry-based proteomic discovery research to identify novel plasma proteins related to the risk of coronary heart disease (CHD) and stroke, and to identify proteins with concentrations affected by the use of postmenopausal hormone therapy. Here we report CHD and stroke risk validation studies for highly ranked proteins, and consider the extent to which protein concentration changes relate to disease risk or provide an explanation for hormone therapy effects on these outcomes.

Methods

Five proteins potentially associated with CHD (beta-2 microglobulin (B2M), alpha-1-acid glycoprotein 1 (ORM1), thrombospondin-1(THBS1), complement factor D pre-protein (CFD), and insulin-like growth factor binding protein 1 (IGFBP1)) and five potentially associated with stroke (B2M, IGFBP2, IGFBP4, IGFBP6, and hemopexin (HPX)) had high discovery phase significance level ranking and an available ELISA assay, and were included in case-control validation studies within the Women’s Health Initiative (WHI) hormone therapy trials. Protein concentrations, at baseline and 1 year following randomization, were assessed for 358 CHD cases and 362 stroke cases, along with corresponding disease-free controls. Disease association, and mediation of estrogen-alone and estrogen plus progestin effects on CHD and stroke risk, were assessed using logistic regression.

Results

B2M, THBS1, and CFD were confirmed (P <0.05) as novel CHD risk markers, and B2M, IGFBP2, and IGFBP4 were confirmed as novel stroke disease risk markers, while the assay for HPX proved to be unreliable. The change from baseline to 1 year in B2M was associated (P <0.05) with subsequent stroke risk, and trended similarly with subsequent CHD risk. Change from baseline to 1 year in IGFBP1 was also associated with CHD risk, and this change provided evidence of hormone therapy effect mediation.

Conclusions

Plasma B2M is confirmed to be an informative risk marker for both CHD and stroke. The B2M increase experienced by women during the first year of hormone therapy trial participation conveys cardiovascular disease risk. The increase in IGFBP1 similarly conveys CHD risk, and the magnitude of the IGFBP1 increase following hormone therapy may be a mediator of hormone therapy effects. Plasma THBS1 and CFD are confirmed as CHD risk markers, and plasma IGFBP4 and IGFBP2 are confirmed as stroke risk markers.

Clinical trials registration

ClinicalTrials.gov identifier: NCT00000611

Background

Cardiovascular disease (CVD), particularly coronary heart disease and stroke, remains the leading cause of death in the United States among both women and men in all racial and ethnic groups. CHD is the designated cause for about 25% of all deaths and stroke for an additional 12% [1]. Risk factor epidemiology has played a crucial role in attempts to understand CVD mechanisms and pathways, and has led to the identification of effective approaches to disease prevention, for example through the treatment of hypertension [2], hypercholesterolemia [3], and arguably chronic inflammation [4]. Risk factor data, such as those arising from the Framingham Study cohort, have been effectively used to develop risk prediction models for CHD [5, 6] and for stroke [7, 8].

Both the age-incidence pattern and the strength of association for some risk factors differ between women and men, and it has been standard procedure to study CVD risk factors and risk prediction models in a sex-specific manner. Some studies have examined the ability of non-traditional risk factors to add to discrimination between CVD cases and controls. For example, an ischemic stroke study [9] showed that the estimated area under the receiver-operator-characteristic curve (AUC) among women in the Atherosclerosis Risk in Communities cohort increased from 0.83 to only 0.84 when certain non-traditional risk markers were included. Corresponding numbers were 0.76 and 0.80 among men. A study of CHD risk prediction models [10] in the Women’s Health Initiative (WHI) postmenopausal hormone therapy (HT) trial cohort found that the AUC increased from 0.73 to 0.75 when certain non-traditional risk factors were added. While these analyses imply an ability to assign CHD and stroke risk estimates that vary by several-fold among individuals, there is still a limited ability to identify individuals who are highly likely to develop disease, say, in the next 5 years. Additional blood-based biomarkers may lead to improvements in risk discrimination and risk prediction.

Blood biomarkers also have potential to provide biological insights into the effects of interventions on CVD outcomes. In particular, the pathways influenced and the key mediators of the effects of postmenopausal estrogen (E-alone) and estrogen plus progestin (E + P) on CVD remain substantially unknown. In the WHI randomized controlled trials these effects include an early elevation of CHD risk with E + P [11] that was less apparent for E-alone [12], and sustained elevations in stroke risk, of a similar magnitude with E + P [13] and E-alone [14].

We have carried out proteomic discovery work using an Intact Protein Analysis System [15] to compare, for about 370 proteins, pre-diagnostic plasma concentrations between CHD cases and matched controls, and between stroke cases and matched controls, drawn from the WHI Observational Study cohort [16]. We confirmed [17] the associations of beta-2-microglobulin (B2M) with short-term CHD risk, and the association of insulin-like growth factor-binding protein 4 (IGFBP4) with short-term stroke risk by comparing baseline blood concentrations for these proteins among women developing these diseases during the first year of participation in the WHI HT trials to 1-1 matched controls, using enzyme-linked immunosorbent assays (ELISAs). The association of these markers with longer-term CHD and stroke risk has yet to be studied. Importantly, there were several other proteins having empirical support for CHD or stroke risk association in our discovery research, with commercially available ELISAs. Here we report on these proteins collectively in relation to CHD and stroke incidence in the WHI HT trials, both as novel disease risk markers and as potential mediators of corresponding postmenopausal HT effects on CHD and stroke risk.

Methods

Study subjects and outcome ascertainment

Women who were postmenopausal and in the age range of 50-79 years enrolled in the WHI HT trials during 1993 to 1998, including 10,739 women who were post-hysterectomy in the E-alone trial and 16,608 women with uterus in the E + P trial. Of these, women who experienced CHD or stroke through February 2001 were included in a CVD biomarker case-control study [18, 19]. Controls who were free of CVD through this date were matched 1-1 to cases on age, randomization date, hysterectomy status, and prevalent study disease at baseline in each trial. These cases and controls were well-characterized in terms of traditional risk factors. One hundred incident stroke cases arising during subsequent HT trial follow-up, and 1-1 matched controls using the same matching criteria, were subsequently added to enhance the ability to study the two diseases separately. A small number of the selected controls developed the disease of their matched case by the end of the planned trial intervention phase (8 April 2005) and are included in the case group here, giving a total of 358 CHD cases and a corresponding 352 controls, and a total of 362 stroke cases and a corresponding 346 matched controls for whom baseline plasma protein concentrations were assessed. Of these, 106 CHD cases and 68 stroke cases had their disease events during the first year following randomization. For all other cases and controls plasma protein concentrations were also assessed from blood drawn at 1 year following randomization.

CHD in the HT trials is defined as non-fatal myocardial infarction (MI) or death due to coronary heart disease. Disease event ascertainment involved physician adjudication based on the review of pertinent documents at each clinical center, and further adjudication by a central committee [20] with agreement rates of 90% for MI and 97% for death due to coronary heart disease, between local and central adjudication. Cases of hospitalized stroke were based on rapid neurologic deficit attributable to arterial obstruction or rupture, or a demonstrable lesion compatible with acute stroke [13]. Central neurologists reviewed all stroke cases, as well as transient ischemic attacks and self-reports of stroke. Strokes were classified as ischemic or hemorrhagic, as well as according to various outcome scales.

All participants provided written informed consent for their HT trial and their overall WHI participation. The related protocols were approved by the Institutional Review Board of the Fred Hutchinson Cancer Research Center and each of the 40 participating clinical centers. The research was conducted in accordance with the Helsinki Declaration and with pertinent local legislation.

Specimen preparation and analysis

Fasting blood specimens were obtained at baseline in WHI as a part of eligibility screening and at 1 year following randomization, for clinical trial women. Serum and plasma were sent to a central laboratory and stored at -70°C. Plasma specimens for this project were plated so that case and matched control specimens were analyzed together. For cases occurring after the first year from randomization, the 1-year plasma specimens were plated with baseline specimens for concurrent ELISA analyses. Assays for each plasma sample followed ELISA kit manufacturer (R and D Systems, Minneapolis, MN, USA, for IGFBP1, IGFBP2, IGFBP4, IGFBP6, ORM1, THBS1, and CHD; CalBiotech, Spring Valley, CA, USA, for B2M; Abnova, Taipei, Taiwan, for HPX) recommendations. All samples were assayed with sample characteristics blinded, and in duplicate. Quality control activities included 5% blind duplicate analyses, using plasma from postmenopausal women outside of the HT trial cohorts. The reliability of the ELISA measurements was assessed by examining intra-class correlations between blind duplicates. Each protein concentration was reliably measured, with intra-class correlations ranging from 0.79 to 0.97, with the exception of HPX where the intra-class correlation was 0.38. Linear dilution curves were examined and a dilution was selected for each analyte that was above the detection threshold and below saturation. The dilutions applied were 1:50, 1:10,000, 1:1,500, 1:3,000, and 1:25 for B2M, ORM1, THSB1, CFD, and IGFBP1, respectively; and 1:250, 1:50, 1:500, and 1:400 for IGFBP2, IGFBP4, IGFBP6, and HPX, respectively.

Proteomic biomarker selection

The in-depth proteomic discovery methodology [15, 21] led to nine proteins in a false discovery rate (FDR) bin [22] of less than 20% for CHD and 11 such proteins for stroke [17]. Of these B2M, ORM1, THBS1, CFD, and IGFBP1 were selected on the basis of not being established as CHD risk markers, and having commercially available ELISA assays. The same criteria applied to the novel stroke candidates led to the selection of IGFBP2, IGFBP4, and IGFBP6, and HPX. B2M was also assessed for stroke cases and controls, based on its CHD association and a nominal P value of 0.03, even though FDR bin was higher (0.31) for this protein.

The CVD proteomic discovery work was complemented by additional IPAS analyses comparing blood protein concentrations at baseline and at 1 year following randomization for 50 women who adhered to active intervention during the first year of the E-alone trial [23], and 50 women who adhered to active intervention in the E + P trial [24]. These analyses suggested many proteomic changes following 1 year of use of these preparations, with 169 (44.7%) of the 378 proteins quantified having some evidence (nominal P <0.05) of change for one or both of E-alone or E + P [24]. Proteins with changed concentrations contributed to multiple biologic pathways relevant to the observed clinical effects of HT, including coagulation, inflammation, immune response, metabolism, cell adhesion, growth factors, and osteogenesis. The estimated 1-year versus baseline concentration ratios were very similar for E-alone and E + P for most highly ranked proteins, supporting the notion of combining proteomic analyses across the two trials.

Of the risk marker candidates selected here, B2M, CFD, and IGFBP1 for CHD, and each of the five proteins selected for stroke were among the proteins whose concentrations were observed to change (nominal P <0.05) as a result of HT [24].

Statistical methods

Principal association analysis estimated CHD or stroke odds ratios (ORs) as a function of log-transformed baseline biomarker values using binary logistic regression of case (1) versus control (0) status [25]. For either disease, the logistic regression model also included systolic and diastolic blood pressure, cigarette smoking, diabetes, prior HT use, and body mass index, as well as the case-control matching factors (age, hysterectomy status, randomization year, prior history of study disease).

Analyses to examine the extent to which treatment-related changes in protein concentrations between baseline and 1 year following randomization can provide an explanation for E-alone and E + P effects on CHD and stroke, also relied on binary logistic regression. Analyses of the type just described, but based on cases occurring after the first year from HT trial enrollment and all controls for the specific trial, were carried out to estimate HT effects on disease OR and to examine the possibility of an OR dependence on the (log-transformed) baseline level of a biomarker under study. A biomarker mediation analysis then proceeded by adding the (log-transformed) year 1 biomarker value or equivalently, the logarithm of the year 1-to-baseline biomarker ratio, to the regression model, and examining the change in the HT OR following the year-1 biomarker addition.

Results

Table 1 shows some characteristics of contributing cases and controls, separately for CHD and stroke, and separately for the E-alone and E + P trials. Compared to E + P trial women, E-alone trial women tended to have higher BMI, and to be more likely to have used postmenopausal hormones prior to trial enrollment.

Table 1 Characteristics of CHD cases and controls and stroke cases and controls from the WHI postmenopausal HT trials of estrogen plus progestin (E + P) and estrogen-alone (E-alone): Mean ± SE for continuous baseline covariates, N (%) for categorical baseline covariates

Table 2 shows geometric means and 95% confidence intervals (CIs) for each selected analyte at baseline, for cases and controls separately, along with P values for their comparison. P values for testing equality of case-control differences between the two trials are also shown, and no between-trial differences were suggested. Combined trial comparisons show CHD case-control differences (P <0.05) for B2M and CFD, and stroke case-control differences for B2M and IGFBP4.

Table 2 Geometric mean and 95% CI for baseline and year-1 biomarker values and for the year 1 to baseline ratio (change) for cases and controls from the WHI postmenopausal HT trials of estrogen plus progestin (E + P) and estrogen-alone (E-alone)

Case-control comparisons based on blood drawn at 1 year following randomization (year 1) are also shown in Table 2, excluding cases occurring in the first year of trial participation. For B2M, combined trial case-control differences at 1 year were evident for both CHD and stroke; and for IGFBP4, were evident for stroke.

More refined baseline analyte comparisons, using logistic regression, are shown in Table 3. These use all cases and controls for each disease, include indicator variables for treatment (active vs. placebo) and trial (E + P vs. E-alone), all control matching variables, and several other CVD risk factors (listed above). The ORs shown are for a 30% increment in the proteomic marker, a value well within the observed range of values for the biomarkers studied. These logistic regression analyses demonstrate positive associations of B2M and CFD, a modest inverse association of THBS1, and a possible positive association of IGFBP1 with CHD risk, each of which concurs in direction with the preceding proteomic discovery results. These analyses also imply positive associations of stroke risk with IGFBP2, IGFBP4, and B2M. Disease risk did not differ between the two trial cohorts, after controlling for the listed factors.

Table 3 ORs (95% CIs) for a 30% baseline biomarker increment from analyses that include data from both WHI postmenopausal HT trials

We previously reported [17] positive associations of baseline B2M with CHD risk and baseline IGFBP4 with stroke risk during the first year from randomization in the WHI HT trials. With the longer-term data analyzed here the estimated CHD OR (95% CI) for a 30% increment in baseline B2M is 1.28 (1.08, 2.54) following year 1, and the stroke OR (95% CI) for a 30% increment in baseline IGFBP4 is 1.16 (1.03, 1.30) following year 1. These are similar to the ORs for the first year [17], and for the overall time period as shown in Table 3.

Additional analyses (not shown) examined ORs for a 30% increase in baseline biomarker, separately in the placebo and treatment groups for each clinical trial.

The B2M associations with risk were evident in both placebo and treatment groups for both diseases. The CFD associations with CHD were primarily evident in the active treatment groups in both trials; whereas the IGFBP4 association with stroke was most evident in the placebo groups for both trials.

There were some noteworthy correlations among the protein concentrations considered for each disease. For CHD, baseline and year 1 B2M concentrations correlated positively with corresponding CFD concentrations in both the placebo and active hormone groups, in both the E-alone and E + P trials. For stroke, positive correlations of B2M with each of IGFBP2, IGFBP4, and IGFBP6 were also evident at baseline and year 1 in both treatment groups, and both trials. Table 4 shows ORs for the five analytes jointly, for each disease, in analyses that otherwise include the same regression variables as Table 3. The strongest associations in these analyses are for B2M and CHD and for IGFBP4 and stroke, while an inverse association of CHD risk with THBS1 and a positive association of stroke risk with IGFBP2 are also observed.

Table 4 ORs (95% CIs) for a 30% baseline biomarker increment, from analyses that combine trials and treatment arms in the WHI postmenopausal HT trials, with biomarkers simultaneously modeled

In Table 5, the year-1 protein concentrations are included in analyses like those in Table 3, but based on cases occurring after year 1 and their disease-specific controls. Estimated ORs are shown for a 30% increment in baseline analyte, and for a 30% increment in the ratio of year 1 to baseline concentration (that is, 30% ‘change’). Toward assessing mediation of HT effects, treatment ORs are shown with only baseline biomarkers included in the analysis, along with corresponding ORs when the analyte change is added to the OR model. For CHD, the treatment OR was essentially unchanged when B2M change was added to the analysis, though there was a suggestion of higher risk (OR of 1.13) among women having a positive concentration change. Interestingly, even though evidence for an association of baseline IGFBP1 with CHD risk was weak, there was a nearly significant association (P = 0.06) of risk with change from baseline to year 1 in IGFBP1 concentration, and the OR for treatment was null (OR = 0.96) after controlling for IGFBP1 change.

Table 5 ORs (95% CIs) for a 30% increment in baseline biomarker, and in year 1 to baseline biomarker ratio (change), based on cases occurring after year 1 and disease-specific controls, from the WHI estrogen plus progestin and estrogen-alone HT trials

For stroke, baseline to 1 year B2M change was positively associated with risk, though the treatment OR was not affected by including B2M change in the analysis. Similarly, OR for treatment was little altered when changes in any of the IGF binding proteins or HPX was added to the analytic model.

Additional analyses (not shown) examined treatment effect mediation by each of these proteins when allowing for technical measurement error in protein assessment, with the blind duplicate data used to estimate the measurement error variance for the (log-transformed) protein concentrations. The treatment ORs were little changed from Table 5 after making this measurement error correction. For example, the treatment effect OR (95% CI), in analyses that include measurement error corrected baseline and 1-year IGFBP1 values, was 0.89 (0.54, 1.45). Other work (also not shown) repeated the Table 5 analyses with cases and controls restricted to women who adhered to their assigned medications during the first year of HT trial participation (at least 80% of pills taken), with little change in HT mediation findings. Further analyses repeated Table 3 excluding women who were being treated for diabetes at baseline (Table 1) with very little change in the ORs for the proteomic markers.

Discussion

Table 6 presents correlations of each of the proteomic markers with plasma measures of lipids, inflammatory factors, and thrombotic factors, as well as insulin, glucose, white blood cell count, and blood pressure, to facilitate the integration of the associations just described with knowledge about CHD and stroke pathogenesis. These measures were available [18, 19] at baseline and 1 year, for all cases and controls, except the 100 stroke case-control pairs that were added late to the case-control study. All variables were log-transformed in calculating these correlations (95% CIs). The correlations shown in Table 6, and ORs for a 30% increment in the proteomic marker from analyses like those shown in Table 3, but with each of the variables on the left side of Table 6 included in the logistic regression model, will be described below for each of the proteomic markers in turn.

Table 6 Pearson correlation (bolded) and 95% CIs for baseline log-transformed proteomic markers and baseline log-transformed variables involved in CHD and stroke pathophysiology

This study confirms plasma B2M to be a risk marker for both CHD and stroke, over an average follow-up period of about 4 years for CHD, and longer for stroke. B2M is an amyloidogenic protein that is elevated in hemodialysis patients [26, 27], and has been reported to be positively associated with CVD risk factors [28] and with CVD events among patients having chronic kidney disease [29], asymptomatic carotid atherosclerosis [30], or peripheral arterial disease in a healthy elderly population [31]. B2M was found to be increased by about 15% by both types of HT both in the control groups studied here and in our preceding proteomics discovery research [23, 24].

From Table 6 it can be seen that B2M has a moderate inverse correlation with HDL-cholesterol, and moderate positive correlations with coagulation factors FVIII and vWF, and with insulin and diastolic blood pressure in this population of postmenopausal women. When the Table 6 variables were added (left side) to the regression model, the OR (95% CI) for a 30% increment in B2M was 1.21 (1.06, 1.37) for CHD, in close proximity to that shown in Table 3 without the addition of these variables, and was 1.46 (1.21, 1.78) for stroke, noticeably stronger than that given in Table 3.

While baseline B2M relates rather clearly to the risk of both CHD and stroke (Table 3), when baseline to 1 year B2M change was added to the regression model (Table 5), it was B2M change that conveyed the greater disease risk, especially for stroke. Recent B2M change deserves consideration as a stroke and possibly a CHD risk marker. However, the HT treatment OR estimates changed little when the B2M change was added to the regression model. Evidently, B2M change is not correlated strongly enough with randomization assignment in the WHI trials for evidence of important mediation of HT effects on these diseases to emerge.

The Table 3 analyses also suggest CFD (adipsin) to be a CHD risk marker, although its association with CHD risk is not significant in analyses that include B2M and the other CHD risk marker candidates (Table 4). CFD correlates inversely with HDL-cholesterol and positively with vWF, insulin, glucose, and diastolic blood pressure. When the Table 6 variables are added to the regression model, the CHD OR (95% CI) for a 30% increment in CFD becomes a non-significant 1.11 (0.96, 1.29). CFD is secreted by adipocytes into the bloodstream. Such adipocytes have been reported to impact multiple functions (blood pressure, lipid metabolism, and hemostasis) to be linked to CVD [32] and to be elevated among obese persons. There is a suggestion (Table 5) that CFD change may relate positively to CHD risk, but the association is not significant.

IGF1 and IGFBP1 have been found to associate positively with all cause and ischemic heart disease mortality in the elderly Rancho Bernardo cohort [33]. Here, baseline IGFBP1 was positively correlated with HDL-cholesterol and inversely associated with glucose. After including these and the other Table 6 measures in the regression analysis, the CHD OR (95% CI) for a 30% IGFBP1 increment was 1.08 (1.02, 1.14). Also, the IGFBP1 change from baseline to year 1 was nearly significant (P = 0.06; Table 5) in its association with disease risk. Moreover, the HT treatment OR was reduced to a null value (0.96) after allowing for the IGFBP1 change. Hence, IGFBP1 deserves consideration among biomarkers that may help to explain HT effects on CHD. This mediation possibility may be dampened, however, by the suggestively larger IGFBP1 changes with E-alone versus E + P [24], whereas CHD ORs were somewhat larger for E + P [11] than for E-alone [12]. The HT regimes studied here are taken orally, and the first-pass hepatic metabolism is known to stimulate a wide variety of proteins. IGFBP1 is recognized as a liver-selective protein [34], so that this protein may be unaffected by transdermal estrogens, that are being increasingly used in clinical practice to treat menopausal symptoms, since these bypass the liver.

THBS1 has a modest inverse association with CHD risk, consistent with our discovery work [17]. This protein has little correlation with the Table 6 factors, though there is some positive association with systolic blood pressure. The OR (95% CI) for a 30% increment in THBS1 is 0.95 (0.91, 1.00) after including all measures on the left side of Table 6 in the analysis. Thrombin, an important factor in relation to inflammation, coagulation, and wound healing, has been shown to regulate THBS1 expression in endothelial cells [35].

The analyses presented here provide little support for ORM1 as a risk marker for CHD among postmenopausal women. It relates inversely to HDL-cholesterol, and positively to IL-6 and MMP-9, D-dimer, vWF, and glucose. When the Table 6 variables are included in the regression analysis, the OR (95% CI) for a 30% increment in ORM1 is 0.84 (0.77, 1.00), suggesting a possible weak inverse incremental association.

There is an extensive literature, reviewed in [36], on the key role of the insulin-like growth factor system in central nervous system development, function, and repair in animal models. The six IGFBPs coordinate and regulate the biologic activity of IGF1 and IGF2. In spite of sequence homology, the IGFBPs may have quite different biological activity, due to differing abundances, with IGFBP-2, -4, and -5 predominating in the brain, and due to post-translational modifications [36, 37]. Also, IGF1 and IGF2 levels have been found to associate inversely with ischemic stroke risk in a Danish case-control study [38], and also with stroke outcome in human studies [39, 40].

Here, our validation exercises confirm a positive association of plasma IGFBP2 and IGFBP4 with stroke risk among postmenopausal women. An association was not confirmed for IGFBP6, which has low abundance in the brain in animal models [36, 37]. From Table 6 one sees that baseline IGFBP4 and IGFBP6 have quite similar correlation patterns with cardiovascular risk biomarkers, with a negative correlation with HDL-cholesterol and positive correlation with IL-6, MMP-9, Factor VIII, insulin, glucose, and diastolic blood pressure. These patterns are also very similar to those for B2M. In contrast, IGFBP2 has a fairly strong positive correlation with HDL-cholesterol, and negative correlations with MMP-9, WBC, and systolic blood pressure. When the Table 6 variables are included in the analysis, the OR (95% CI) for a 30% increment in IGFBP2 is 1.21 (1.08, 1.35), larger than that given in Table 3, while corresponding ORs (95% CIs) are 1.14 (0.99, 1.32) for IGFBP4 and 1.06 (0.95, 1.17) for IGFBP6. In conjunction with Table 4, one can infer that IGFBP4 and IGFBP2 are associated with stroke risk beyond that attributable to the established biomarkers considered here, while there is little evidence for further association with IGFBP6. However, there was limited evidence of important mediation of HT effects on stroke by either IGFBP4 or IGFBP2, and only a weak suggestion of a positive association between IGFBP4 change from baseline to year 1 and stroke risk.

Hemopexin had been shown in mice to be neuroprotective through high-affinity binding of the pre-oxidant free heme [41]. The baseline HPX measures obtained here correlated negatively with HDL-cholesterol and positively with IL-6, Factor 1.2, and glucose, and the OR (95% CI) for a 30% increment in HPX was 1.05 (0.87, 1.28) after including each of the Table 6 variables in the analysis. Overall, neither baseline nor change in HPX was related to stroke risk in these analyses, but this could be due to the poor reliability of the HPX ELISA assay used here.

The proteins studied here provide some interesting leads concerning the pathophysiology of both CHD and stroke. Their ability to enhance discrimination between cases and controls, however, appears to be limited. For example, when baseline values of these proteins were added to a model that includes the variables used here to control confounding (Table 3 analyses), the AUC for CHD did not increase from its value of 0.670 with 95% CI of (0.615, 0.726) without such addition. For stroke there was some modest AUC increase from 0.645 (0.590, 0.700) without any such proteins added, to 0.663 (0.604, 0.718) when B2M was added, and to 0.665 (0.605, 0.722) when IGFBP4 was added, based on analyses that randomly divided the data into training and validation subsets with the model fitted in the training set used to estimate AUC in the validation set.

Using the same split sample approach, one can estimate positive (PPV) and negative (NPV) predictive values for the proteins for which there is evidence of disease risk association at, say, a specificity of 80%. The PPV and NPV estimates for CHD are 0.70 and 0.58, respectively, without inclusion of the proteins studied here, and are essentially unchanged when any of B2M, CFD, THBS1, or IGFBP1 is added to the model. The corresponding estimated PPV and NPV values of 0.62 and 0.54 without the proteins evaluated increased slightly to 0.68 and 0.58 when B2M was added, to 0.66 and 0.57 when IGFBP4 was added, and to 0.64 and 0.56 when IGFBP2 was added to the regression model. The highly overlapping distribution of the proteomic markers between cases and controls (Table 2) prevents these measures from adding much to case versus control discrimination or, presumably, to personalized risk assessment. Even though the discovery and validation phases of this work took place in distinct cohorts, the WHI observational study and clinical trial, respectively, the two cohorts were drawn from essentially the same catchment population, and evaluation of these findings in other populations will be useful.

The set of proteins studied here was limited by our requirement of a commercially available ELISA and not recognized as CHD or stroke risk markers. There were several other proteins within FDR <0.20 bins that could be evaluated, perhaps using multiple reaction monitoring mass spectrometry. These proteins are listed in Tables 1 and 2 of [17].

Conclusions

Proteomic discovery work has led to the identification of plasma B2M, CFD, and THBS1 as novel risk markers for CHD. Additionally, an increase in IGFBP1 over a 1-year period also appears to convey CHD risk, and may be relevant to the early elevation in CHD risk among women initiating oral HT. This work has also led to the identification of plasma B2M, IGFBP2, and especially IGFBP4 as novel risk markers for stroke risk among postmenopausal women.

Abbreviations

AUC:

Area under receiver-operator curve

B2M:

Beta-2-microglobulin

CFD:

Complement factor D pre-protein

CHD:

Coronary heart disease

CI:

Confidence interval

CRP:

C-reactive protein

CVD:

Cardiovascular disease

E + P:

Estrogen plus progestin

E-alone:

Estrogen-alone

ELISA:

Enzyme-linked immunosorbent assay

Factor 1.2:

Prothrombin fragment 1.2

Factor VIII:

Clotting factor eight

FDR:

False discovery rate

HDL-C:

High-density lipoprotein cholesterol

HPX:

Hemopexin

HT:

(postmenopausal) hormone therapy

IGFBP:

Insulin-like growth factor binding protein

IL-6:

Interleukin 6

LDL-C:

Low-density lipoprotein cholesterol

MMP-9:

Matrix metallaprotinase-9

OR:

Odds ratio

ORM1:

Alpha-1-acid glycoprotein 1

THBS1:

Thrombospondin-1

vWF:

Von villebrand factor

WBC:

White blood cell count

WHI:

Women’s Health Initiative.

References

  1. King HC, Hoyert DL, Xu JQ, Murphy SL: Deaths: Final Data for 2005. National Vital Statistics Reports. 2008, Hyattsville, MD: National Center for Health Statistics, 56:

    Google Scholar 

  2. Law MR, Morris JK, Wald NJ: Use of blood pressure lowering drugs in the prevention of cardiovascular disease: meta-analysis of 147 randomised trials in the context of expectations from prospective epidemiological studies. BMJ. 2009, 338: b1665-10.1136/bmj.b1665.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Shepherd J, Cobbe SM, Ford I, Isles CG, Lorimer AR, MacFarlane PW, McKillop JH, Packard CJ, West of Scotland Coronary Prevention Study Group: Prevention of coronary heart disease with pravastatin in men with hypercholesterolemia. N Engl J Med. 1995, 333: 1301-1307. 10.1056/NEJM199511163332001.

    Article  CAS  PubMed  Google Scholar 

  4. Ridker PM, Danielson E, Fonseca FA, Genest J, Gotto AM, Kastelein JJ, Koenig W, Libby P, Lorenzatti AJ, MacFadyen JG, Nordestgaard BG, Shepherd J, Willerson JT, Glynn RJ, JUPITER Study Group: Rosuvastatin to prevent vascular events in men and women with elevated C-reactive protein. N Engl J Med. 2008, 359: 2195-2207. 10.1056/NEJMoa0807646.

    Article  CAS  PubMed  Google Scholar 

  5. Gordon T, Kannel WB: Multiple risk function for predicting coronary heart disease: the concept, accuracy, and application. Am Heart J. 1982, 103: 1031-1039. 10.1016/0002-8703(82)90567-1.

    Article  CAS  PubMed  Google Scholar 

  6. D’Agostino RB, Grundy S, Sullivan LM, Wilson P, for the CHD Risk Prediction Group: Validation of the Framingham coronary heart disease prediction scores: results of a multiple ethnic groups investigation. JAMA. 2001, 286: 180-187. 10.1001/jama.286.2.180.

    Article  PubMed  Google Scholar 

  7. Wolf PA, D’Agostino RB, Belanger AJ, Kannel WB: Probability of stroke: a risk profile from the Framingham Study. Stroke. 1991, 22: 312-318. 10.1161/01.STR.22.3.312.

    Article  CAS  PubMed  Google Scholar 

  8. D’Agostino RB, Wolf PA, Belanger AJ, Kannel WB: Stroke risk profile: adjustment for antihypertensive medication. Stroke. 1994, 25: 40-43. 10.1161/01.STR.25.1.40.

    Article  PubMed  Google Scholar 

  9. Chambless LE, Heiss G, Shahar E, Earp MJ, Toole J: Prediction of ischemic stroke risk in the atherosclerosis risk in communities study. Am J Epidemiol. 2004, 160: 259-269. 10.1093/aje/kwh189.

    Article  PubMed  Google Scholar 

  10. Kim HC, Greenland P, Rossouw JE, Manson JE, Cochrane BB, Lasser NL, Limacher MC, Lloyd-Jones DM, Margolis KL, Robinson JG: Multimarker prediction of coronary heart disease risk. J Am Coll Cardiol. 2010, 55: 2080-2091. 10.1016/j.jacc.2009.12.047.

    Article  CAS  PubMed  Google Scholar 

  11. Manson JE, Hsia J, Johnson KC, Rossouw JE, Assaf AR, Lasser NL, Trevisan M, Black HR, Heckbert SR, Detrano R, Strickland OL, Wong ND, Crouse JR, Stein E, Cushman M, Women’s Health Initiative Investigators: Estrogen plus progestin and the risk of coronary heart disease. N Engl J Med. 2003, 349: 523-534. 10.1056/NEJMoa030808.

    Article  CAS  PubMed  Google Scholar 

  12. Hsia J, Langer RD, Manson JE, Kuller L, Johnson KC, Hendrix SL, Pettinger M, Heckbert SR, Greep N, Crawford S, Eaton CB, Kostis JB, Caralis P, Prentice R, Women’s Health Initiative Investigators: Conjugated equine estrogens and coronary heart disease: The Women’s Health Initiative. Arch Intern Med. 2006, 166: 357-365. 10.1001/archinte.166.3.357.

    Article  CAS  PubMed  Google Scholar 

  13. Wassertheil-Smoller S, Hendrix SL, Limacher M, Heiss G, Kooperberg C, Baird A, Kotchen T, Curb JD, Black H, Rossouw JE, Aragaki A, Safford M, Stein E, Laowattana S, Mysiw WJ, WHI Investigators: Effect of estrogen plus progestin on stroke in postmenopausal women: the Women’s Health Initiative: A randomized trial. JAMA. 2003, 289: 2673-2684. 10.1001/jama.289.20.2673.

    Article  CAS  PubMed  Google Scholar 

  14. Hendrix SL, Wassertheil-Smoller S, Johnson KC, Howard BV, Kooperberg C, Rossouw JE, Trevisan M, Aragaki A, Baird AE, Bray PF, Buring JE, Criqui MH, Herrington D, Lynch JK, Rapp SR, Torner J, WHI Investigators: Effects of conjugated equine estrogen on stroke in the Women’s Health Initiative. Circulation. 2006, 113: 2425-2434. 10.1161/CIRCULATIONAHA.105.594077.

    Article  CAS  PubMed  Google Scholar 

  15. Faca V, Coram M, Phanstiel D, Glukhova V, Zhang Q, Fitzgibbon M, McIntosh M, Hanash S: Quantitative analysis of acrylamide labeled serum proteins by LC-MS/MS. J Proteome Res. 2006, 5: 2009-2018. 10.1021/pr060102+.

    Article  CAS  PubMed  Google Scholar 

  16. The Women’s Health Initiative Study Group: Design of the Women’s Health Initiative clinical trial and observational study. Control Clin Trials. 1998, 19: 61-109.

    Article  Google Scholar 

  17. Prentice RL, Paczesny SJ, Aragaki A, Amon LM, Chen L, Pitteri SJ, McIntosh M, Wang P, Busald Buson T, Hsia J, Jackson RD, Rossouw JE, Manson JE, Johnson K, Eaton C, Hanash SM: Novel proteins associated with risk for coronary heart disease or stroke among postmenopausal women identified by in-depth plasma proteome profiling. Genome Med. 2010, 2: 48-60. 10.1186/gm169.

    Article  PubMed Central  PubMed  Google Scholar 

  18. Rossouw JE, Cushman M, Greenland P, Lloyd-Jones DM, Bray P, Kooperberg C, Pettinger M, Robinson J, Hendrix S, Hsia J: Inflammatory, lipid, thrombotic, and genetic markers of coronary heart disease in the Women’s Health Initiative trials of hormone therapy. Arch Intern Med. 2008, 168: 2245-2253. 10.1001/archinte.168.20.2245.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Kooperberg C, Cushman M, Hsia J, Robinson JG, Aragaki AK, Lynch JK, Band AE, Johnson KC, Kuller LH, Beresford SAA, Rodriguez B: Can biomarkers identify women at increased stroke risk? The Women’s Health Initiative hormone trials. PLoS Clin Trials. 2007, 2: e28-10.1371/journal.pctr.0020028.

    Article  PubMed Central  PubMed  Google Scholar 

  20. Curb JD, McTiernan A, Heckbert SR, Kooperberg C, Stanford J, Nevitt M, Johnson KC, Proulx-Burns L, Pastore L, Criqui M, Daugherty S, WHI Morbidity and Mortality Committee: Outcomes ascertainment and adjudication methods in the Women’s Health Initiative. Ann Epidemiol. 2003, 13: 122-128. 10.1016/S1047-2797(03)00048-6.

    Article  Google Scholar 

  21. Rauch A, Bellew M, Eng J, Fitzgibbon M, Holzman T, Hussey P, Igra M, Maclean B, Lin CW, Detter A, Fang R, Faca V, Gafken P, Zhang H, Whiteaker J, States D, Hanash S, Paulovich A, McIntosh MW: Computational Proteomics Analysis System (CPAS): an extensible, open-source analytic system for evaluating and publishing proteomic data and high throughput biological experiments. J Proteome Res. 2006, 5: 112-121. 10.1021/pr0503533.

    Article  CAS  PubMed  Google Scholar 

  22. Benjamini Y, Hochberg Y: Controlling the False Discovery Rate: a practical and powerful approach to multiple testing. J Roy Stat Soc B (Methodological). 1995, 57: 289-300.

    Google Scholar 

  23. Katayama H, Pacznesny S, Prentice RL, Aragaki A, Faca VM, Pitteri SJ, Zhang Q, Wang H, Silva M, Kennedy J, Rossouw J, Jackson R, Hsia J, Chlebowski R, Manson J, Hanash S: Application of serum proteomics to the Women’s Health Initiative conjugated equine estrogens trial reveals a multitude of effects relevant to clinical findings. Genome Med. 2009, 1: 47.1-47.16.

    Article  Google Scholar 

  24. Pitteri SJ, Hanash SM, Aragaki A, Amon L, Chen L, Busald Buson T, Paczesny S, Katayama H, Wang H, Johnson MM, Zhang Q, McIntosh M, Wang P, Kooperberg C, Rossouw JE, Jackson R, Manson JE, Hsia J, Liu S, Martin L, Prentice RL: Postmenopausal estrogen and progestin effects on the serum proteome. Genome Med. 2009, 1: 121.1-121.14.

    Article  Google Scholar 

  25. Prentice RL, Pyke R: Logistic disease incidence models and case-control studies. Biometrika. 1979, 66: 403-412. 10.1093/biomet/66.3.403.

    Article  Google Scholar 

  26. Polat H, Yeksan M, Dalmaz M, Kaptanoglu B, KoÅŸar A, AkkuÅŸ I: Serum amyloid A protein levels in haemodialysis patients. Nephrol Dial Transplant. 1996, 11: 1492-1493.

    Article  CAS  PubMed  Google Scholar 

  27. Saijo Y, Utsugi M, Yoshioka E, Horikawa N, Sato T, Gong Y, Kishi R: Relationship of B2-Microglobulin to arterial stiffness in Japanese subjects. Hypertens Res. 2005, 28: 505-511. 10.1291/hypres.28.505.

    Article  CAS  PubMed  Google Scholar 

  28. Wilson AM, Kimura E, Harada RK, Nair N, Narasimhan B, Meng X-Y, Zhang F, Beck KR, Olin JW, Fung ET, Cooke JP: B2-Microglobulin as a biomarker in peripheral arterial disease: proteomic profiling and clinical studies. Circulation. 2007, 116: 1396-1403. 10.1161/CIRCULATIONAHA.106.683722.

    Article  CAS  PubMed  Google Scholar 

  29. Liabeuf S, Lenglet A, Desjardins L, Neirynck N, Glorieux G, Lemke H-D, Vanholder R, Diouf M, Choukroun G, Massy ZA, on behalf of the European Uremic Toxin Work Group (EUTox): Plasma beta-2 microglobulin is associated with cardiovascular disease in uremic patients. Kidney Int. 2012, 82: 1297-1303. 10.1038/ki.2012.301.

    Article  CAS  PubMed  Google Scholar 

  30. Amighi J, Hoke M, Mlekusch W, Schlager O, Exner M, Haumer M, Pernicka E, Koppensteiner R, Minar E, Rumpold H, Schillinger M, Wagner O: Beta 2 microglobulin and the risk for cardiovascular events in patients with asymptomatic carotid atherosclerosis. Stroke. 2011, 42: 1826-1833. 10.1161/STROKEAHA.110.600312.

    Article  CAS  PubMed  Google Scholar 

  31. Shinkai S, Chaves PHM, Fujiwara Y, Watanabe S, Shibata H, Yoshida H, Suzuki T: B2-Microglobulin for risk stratification of total mortality in the elderly population. Comparison with Cystatin C and C-Reactive Protein. Arch Intern Med. 2008, 168: 200-206. 10.1001/archinternmed.2007.64.

    Article  CAS  PubMed  Google Scholar 

  32. Ronti T, Lupattelli G, Mannarino E: The endocrine function of adipose tissue: an update. Clin Endocrinol. 2006, 64: 355-365.

    CAS  Google Scholar 

  33. Laughlin GA, Barrett-Connor E, Criqui MH, Kritz-Silverstein D: The prospective association of serum Insulin-like Growth Factor I (IGF-1) and IGF-Binding Protein-1 levels with all cause and cardiovascular disease mortality in older adults: the Rancho Bernardo Study. J Clin Endocrinol Metab. 2004, 89: 114-120. 10.1210/jc.2003-030967.

    Article  CAS  PubMed  Google Scholar 

  34. Hsiao LL, Dangong F, Yoshida T, Hong R, Jensen RV, Misra J, Dillon W, Lee KF, Clark KE, Haverty P, Weng Z, Mutter GL, Frosch MP, Macdonald ME, Milford EL, Crum CP, Bueno R, Pratt RE, Mahadevappa M, Warrington JA, Stephanoupoulos G, Gullans SR: A compendium of dene expression in normal human tissues. Physiol Genomics. 2001, 7: 97-104.

    Article  CAS  PubMed  Google Scholar 

  35. McLaughlin JN, Mazzoni MR, Cleator JH, Earls L, Perdigoto AL, Brooks JD, Muldowney JAS, Vaughan DE, Hamm HE: Thrombin modulates the expression of a set of genes including thrombospondin-1 in human microvascular endothelial cells. J Biol Chem. 2005, 280: 22172-22180. 10.1074/jbc.M500721200.

    Article  CAS  PubMed  Google Scholar 

  36. Russo VC, Gluckman PD, Feldman EL, Werther GA: The insulin-like growth factor system and its pleiotropic functions in brain. Endocr Rev. 2005, 26: 916-943. 10.1210/er.2004-0024.

    Article  CAS  PubMed  Google Scholar 

  37. Kooijman R, Sarre S, Michotte Y, De Keyser J: Insulin-like growth factor I: a potential neuroprotective compound for the treatment of acute ischemic stroke?. Stroke. 2009, 40: e83-e88. 10.1161/STROKEAHA.108.528356.

    Article  PubMed  Google Scholar 

  38. Johnsen SP, Hundborg HH, Sorensen HT, Orskov H, Tjonneland A, Overvad K, Jorgensen JOL: Insulin-like Growth Factor (IGF) I, -II, and IGF Binding Protein-3 and risk of ischemic stroke. J Clin Endocrinol Metab. 2005, 90: 5937-5941. 10.1210/jc.2004-2088.

    Article  CAS  PubMed  Google Scholar 

  39. Denti L, Annoni V, Cattadori E, Salvagnini MA, Visioli S, Merli MF, Corradi F, Ceresini G, Valenti G, Hoffman AR, Ceda GP: Insulin-like growth factor 1 as a predictor of ischemic stroke outcome in the elderly. Am J Med. 2004, 117: 312-317. 10.1016/j.amjmed.2004.02.049.

    Article  CAS  PubMed  Google Scholar 

  40. Bondanelli M, Ambrosio MR, Onofri A, Bergonzoni A, Lavezzi S, Zatelli MC, Valle D, Basaglia N, degli Uberti EC: Predictive value of circulating insulin-like growth factor I levels in ischemic stroke outcome. J Clin Endocrinol Metab. 2006, 91: 3928-3934. 10.1210/jc.2006-1040.

    Article  CAS  PubMed  Google Scholar 

  41. Li RC, Saleem S, Zhen G, Cao W, Zhuang H, Lee J, Smith A, Altruda F, Tolosano E, Dore S: Heme-hemopexin complex attenuates neuronal cell death and stroke damage. J Cereb Blood Flow Metab. 2009, 29: 953-964. 10.1038/jcbfm.2009.19.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

A list of key investigators involved in WHI follows:

Program Office: (National Heart, Lung, and Blood Institute, Bethesda, Maryland) Jacques Rossouw, Shari Ludlam, Dale Burwen, Joan McGowan, Leslie Ford, and Nancy Geller.

Clinical Coordinating Center: Clinical Coordinating Center: (Fred Hutchinson Cancer Research Center, Seattle, WA) Garnet Anderson, Ross Prentice, Andrea LaCroix, and Charles Kooperberg.

Investigators and Academic Centers: (Brigham and Women’s Hospital, Harvard Medical School, Boston, MA) JoAnn E. Manson; (MedStar Health Research Institute/Howard University, Washington, DC) Barbara V. Howard; (Stanford Prevention Research Center, Stanford, CA) Marcia L. Stefanick; (The Ohio State University, Columbus, OH) Rebecca Jackson; (University of Arizona, Tucson/Phoenix, AZ) Cynthia A. Thomson; (University at Buffalo, Buffalo, NY) Jean Wactawski-Wende; (University of Florida, Gainesville/Jacksonville, FL) Marian Limacher; (University of Iowa, Iowa City/Davenport, IA) Robert Wallace; (University of Pittsburgh, Pittsburgh, PA) Lewis Kuller; (Wake Forest University School of Medicine, Winston-Salem, NC) Sally Shumaker.

Women’s Health Initiative Memory Study: (Wake Forest University School of Medicine, Winston-Salem, NC) Sally Shumaker.

For a list of all the investigators who have contributed to WHI science, please visit: https://cleo.whi.org/researchers/Documents%20%20Write%20a%20Paper/WHI%20Investigator%20Long%20List.pdf.

Funding/support

This work was supported by National Heart, Lung and Blood Institute, National Institutes of Health, U.S. Department of Health and Human Services grant R21HL105927 (Prentice). The Women’s Health Initiative Program (all investigators) is also supported by the National Heart, Lung and Blood Institute under contracts HHSN268201100046C, HHSN268201100001C, HHSN268201100002C, HHSN268201100003C, HHSN268201100004C, and HHSN271201100004C. Dr. Prentice’s work was also partially supported by grant P01CA53996 from the National Cancer Institute.

Clinical Trials Registration: ClinicalTrials.gov identifier: NCT00000611

Role of the sponsor

Decisions concerning study design, data collection and analysis, interpretation of the results, the preparation of the manuscript, or the decision to submit the manuscript for publication resided with committees comprised of WHI investigators that included National Heart, Lung and Blood Institute representatives.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Ross L Prentice.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

RLP and SMH secured funding, oversaw the data collection, and drafted the manuscript. MJ conducted the ELISA assays and related quality assurance activities. AA assembled project data and provided ongoing data analyses support. SZ was the principal data analyst for manuscript development, and shared in the initial drafting of the manuscript. JH, RDJ, JER, and JEM served as consultants throughout the project, providing cardiovascular disease expertise and input to analyte selection. All authors provided review and critical input to manuscript development and finalization. All authors read and approved the final manuscript.

Rights and permissions

Reprints and permissions

About this article

Cite this article

Prentice, R.L., Zhao, S., Johnson, M. et al. Proteomic risk markers for coronary heart disease and stroke: validation and mediation of randomized trial hormone therapy effects on these diseases. Genome Med 5, 112 (2013). https://doi.org/10.1186/gm517

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/gm517

Keywords