Postmenopausal estrogen and progestin effects on the serum proteome

Background Women's Health Initiative randomized trials of postmenopausal hormone therapy reported intervention effects on several clinical outcomes, with some important differences between estrogen alone and estrogen plus progestin. The biologic mechanisms underlying these effects, and these differences, have yet to be fully elucidated. Methods Baseline serum samples were compared with samples drawn 1 year later for 50 women assigned to active hormone therapy in both the estrogen-plus-progestin and estrogen-alone randomized trials, by applying an in-depth proteomic discovery platform to serum pools from 10 women per pool. Results In total, 378 proteins were quantified in two or more of the 10 pooled serum comparisons, by using strict identification criteria. Of these, 169 (44.7%) showed evidence (nominal P < 0.05) of change in concentration between baseline and 1 year for one or both of estrogen-plus-progestin and estrogen-alone groups. Quantitative changes were highly correlated between the two hormone-therapy preparations. A total of 98 proteins had false discovery rates < 0.05 for change with estrogen plus progestin, compared with 94 for estrogen alone. Of these, 84 had false discovery rates <0.05 for both preparations. The observed changes included multiple proteins relevant to coagulation, inflammation, immune response, metabolism, cell adhesion, growth factors, and osteogenesis. Evidence of differential changes also was noted between the hormone preparations, with the strongest evidence in growth factor and inflammation pathways. Conclusions Serum proteomic analyses yielded a large number of proteins similarly affected by estrogen plus progestin and by estrogen alone and identified some proteins and pathways that appear to be differentially affected between the two hormone preparations; this may explain their distinct clinical effects.

B Ba ac ck kg gr ro ou un nd d Postmenopausal hormone therapy was shown to have multiple effects of public-health importance in the Women's Health Initiative (WHI) randomized, placebo-controlled hormone-therapy trials of 0.625 mg/day conjugated equine estrogen (E-alone) [1] or of this same estrogenic preparation plus 2.5 mg/day medroxyprogesterone acetate (E+P) [2], over respective average intervention periods of 7.1 and 5.6 years. The observed effects were similar for the two preparations for some outcomes, including stroke [3,4] and hip fracture [5,6]; whereas E+P effects were unfavorable (P < 0.05) compared with those for E-alone for other outcomes, including coronary heart disease (CHD) [7,8], breast cancer [9,10], and venous thromboembolism (VT) [11,12], and a global index [1,2] that was designed to summarize major health benefits versus risks [13].
Several of the articles just cited formally examined whether interactions occurred between the hormone-therapy hazard ratios and baseline study-subject characteristics. Although some moderate variations were detected (for example, for E-alone and breast cancer [10]), these tended to provide limited insight into the biologic mechanisms and pathways involved in the observed clinical effects. A cardiovascular disease nested case-control study also was conducted to relate baseline values of candidate biomarkers and postrandomization biomarker changes to observed hormonetherapy effects. This study confirmed baseline biomarker disease associations and identified some pertinent biomarker changes after hormone-therapy initiation, but identified few interactive or explanatory biomarkers for either CHD [14] or stroke [15], although the E+P hazard ratio elevation for CHD appeared to be smaller among women having relatively low baseline low-density lipoprotein cholesterol [14].
It follows that much remains to be explained about the pattern of biologic changes induced by these hormone-therapy preparations in relation to the outcome effects mentioned earlier. Proteomic discovery work has the potential to identify biomarkers that may help to explain E+P or E-alone clinical effects or differences in effects between the two preparations. Hence, we applied a comprehensive quantitative proteomic approach designated Intact Protein Analysis System (IPAS) [16][17][18][19] to compare the serum proteome at 1 year after randomization to baseline for 50 women assigned to E+P and for 50 women assigned to E-alone, in the WHI hormone-therapy trials. These women were selected to be free of major disease outcomes through the WHI clinical trial intervention phase and were selected to be adherent to their assigned hormone regimen over the first year of treatment, but were otherwise randomly selected from women assigned to active treatment in the trial cohorts. The IPAS approach involves extensive fractionation followed by tandem mass spectrometry and is capable of identifying proteins over seven orders of abundance. For reasons of throughput, serum pools were formed from 10 E+P women (five baseline and five 1-year pools), or from 10 E-alone women, before proteomic analysis.
We recently reported [20] proteomic changes from the E-alone component of this project. An impressive 10.5% of proteins had false discovery rates, for a change, of <0.05. The affected proteins had relevance to multiple pathways, including coagulation, metabolism, osteogenesis, and inflammation, among others. Ten of 14 protein changes tested were confirmed with enzyme-linked immunosorbent assays (ELISAs) in the original samples, and in serum samples from 50 nonoverlapping randomly chosen women, selected by using the same criteria, from the E-alone trial treatment group.
Here, we sought to uncover proteins and pathways that are differentially affected by E+P therapy relative to E-alone that would provide leads for the comparatively unfavorable effects with E+P observed in these trials.

M Me et th ho od ds s S St tu ud dy y c co oh ho or rt ts s
The use of human samples was approved by the Fred Hutchinson Cancer Research Center Institutional Review Board. Fifty study subjects were randomly selected from the 8,506 women assigned to active E+P in the WHI clinical trial, which also included 8,102 women assigned to placebo. All women were postmenopausal, with a uterus, and in the age range from 50 to 79 years, at recruitment during 1993 through 1998. The selected women were required to have been adherent to study medication (80% or more of pills taken) over the first year after randomization, and without a major clinical event (CHD, stroke, VT, breast or colorectal cancer, or hip fracture) over the intervention and follow-up period (through March 2005). A second nonoverlapping subset of E+P women was selected, by using the same criteria, for replication studies with ELISA. As previously reported [20], the same selection criteria were used for the E-alone discovery and replication phases of the study. Women enrolled in the E-alone trial (10,739) satisfied the same eligibility criteria as E+P enrollees, but were posthysterectomy at randomization. Women who used hormone therapy before trial enrollment had mostly stopped such treatment, months or years before enrollment, and were otherwise required to undergo a 3-month washout before randomization. Serum samples, collected at baseline and 1 year, were stored at -80° C until proteomic analyses.
S Sa am mp pl le e p pr re ep pa ar ra at ti io on n, , p pr ro ot te ei in n f fr ra ac ct ti io on na at ti io on n, , a an nd d m ma as ss s s sp pe ec ct tr ro om me et tr ry y a an na al ly ys si is s These methods were previously described [20] in detail and are only briefly summarized here. As in the E-alone project component, pools formed from 30 µl of serum for 10 randomly selected women from the 50 E+P group women were formed from baseline and 1-year specimens.
After immunodepletion of the six most abundant proteins (albumin, IgG, IgA, transferrin, haptoglobin, and antitrypsin), pools were concentrated, and intact proteins having cysteine residues were isotopically labeled with acrylamide (baseline pools received the 'light' C12 acrylamide; 1-year pools the 'heavy' C13 acrylamide). The baseline and 1-year pools were than mixed together for further analysis.
The combined sample was diluted, and each sample was separated into 12 subsamples by using anion exchange chromatography, and each subsample was further separated into 60 fractions by using reversed-phase chromatography, giving a total of 720 fractions for each original mixed sample. Aliquots of 200 µl from each fraction, corresponding to about 200 µg of protein, were separated for mass spectrometry 'shotgun' analysis.
Lyophilized aliquots from the 720 individual fractions were subjected to in-solution trypsin digestion, and individual digested fractions, four to 60 from each reversed-phase run, were combined into 11 pools, giving a total of 132 (12 × 11) fractions for analysis from each original mixed baseline and 1-year pool. Tryptic peptides were analyzed with an LTQ-FT mass spectrometer. Spectra were acquired in a data-dependent mode in a mass/charge range of 400 to 1,800, and the five most abundant +2 or +3 ions were selected from each spectrum for tandem mass spectrometry (MS/MS) analysis. P Pr ro ot te ei in n i id de en nt ti if fi ic ca at ti io on n a an nd d b ba as se el li in ne e v ve er rs su us s 1 1--y ye ea ar r c co on nc ce en nt tr ra at ti io on n a as ss se es ss sm me en nt t The acquired LC-MS/MS data were automatically processed by the Computational Proteomics Analysis System [21]. Database searches were performed by using X!Tandem against the human International Protein Index (IPI) by using tryptic search. Database search results were analyzed by using PeptideProphet [22] and ProteinProphet [23].
The relative quantitation of 1-year to baseline concentration for cysteine-containing peptides identified by MS/MS was extracted by using a script designated Q3 ProteinRatioParser [16], which calculates the relative peak areas of heavy to light acrylamide-labeled peptides. Peaks with zero area were reset to a background value to avoid singularities. Peptides having PeptideProphet ≥0.75, Tandem expect score <0.10, and mass deviation <20 ppm were considered for quantification. Proteins were identified as those having ProteinProphet scores ≥0.90, and their ratios were calculated by taking the geometric mean of all the associated peptide ratios. Proteins from all 10 IPAS experiments were aligned by their protein group number, assigned by ProteinProphet, to identify master groups of indistinguishable proteins across experiments. Ratios for these protein groups were logarithmically transformed and median-centered at zero. The following protein groups were removed in this analysis: groups that had fewer than five peptide ratios across all 10 experiments; groups that contained proteins that were targeted for depletion; and groups in which all proteins had been annotated as 'defunct' by IPI.
S St ta at ti is st ti ic ca al l a an na al ly ys si is s o of f 1 1--y ye ea ar r v ve er rs su us s b ba as se el li in ne e p pr ro ot te ei in n c co on nc ce en nt tr ra at ti io on ns s Protein log-(concentration) ratios were analyzed by first normalizing further, so that the median of the log-ratios is zero for all the proteins identified from a mixed baseline and 1-year sample. Concentration changes after E+P use were identified by testing the hypothesis that the mean of the logratios across the (up to 5) mixed samples is zero, by using a weighted moderated t statistic [24] implemented in the R package LIMMA [25]: the log-ratios were weighted by the number of quantified peptides for each protein, and a matrix of weights was included in the linear model. The variance was estimated by using the sum of the sample variances from the E+P and E-alone data, with the requirement of at least one degree of freedom for variance estimation. Benjamini and Hochberg's method [26] was used to accommodate multiple testing for the large number of proteins quantified, through the calculation of estimated false discovery rates (FDRs).
The same method was used to identify proteins for which the 1-year to baseline change in concentration differed between E+P and E-alone. Specifically, a moderated t statistic was used to test for a difference in means between the log-ratios for E+P and those for E-alone, with common log-ratio variance for the two preparations.
B Bi io ol lo og gi ic c p pa at th hw wa ay y a an na al ly ys si is s We developed a regularized Hotelling T 2 procedure (Chen LS, Prentice RL, and Wang P, submitted for publication, 2009) to identify sets of proteins, defined by biologic pathways, that change concentration with E+P, or that change differentially in the E+P and E-alone project components. This testing procedure takes advantage of the correlation structure among the log-ratios for proteins in a given set. Protein sets were defined by using the KEGG database [27,28].
To accommodate multiple hypotheses testing issues, the significance for individual proteins or for biologic pathways, is based on a 5% FDR criterion.
E EL LI IS SA A--b ba as se ed d v va al li id da at ti io on n ELISAs are commercially available for some of the proteins for which evidence emerged of change after E+P use, or of differential change between E+P and E-alone. ELISA tests were applied according to manufacturer's protocols for individual baseline and 1-year serum samples from an additional randomly selected nonoverlapping 50 E+P and 50 E-alone women, for independent validation of leads from the proteomic discovery work. P values were obtained by applying t tests to log-transformed 1-year-to-baseline concentration ratios. Log-ratios from ELISA and IPAS were compared to assess discovery platform signals.

R Re es su ul lt ts s
The average age at enrollment for the selected 50 E+P women is 63.2, similar to that for the trial cohort as a whole. Other study-subject characteristics were generally similar also to those for the entire trial cohort [2], as was also the case for the 50 selected E-alone women [20]. Subject characteristics for both studies are shown in Table 1. Some characteristics varied among the pools of size 10, as expected with the random assignment of women to pools. For example, the average baseline age (standard deviation) for the five E+P pools was 60.6 (8.4), 65.8 (5.3), 63.5 (8.5), 63.2 (7.1), and 62.8 (7.0), respectively. The project generated 2,576,869 spectra from the E+P pools, as compared with 2,458,506 from the E-alone pools. These led to the identification of 3,669 IPI-based proteins from the E+P pools compared with 4,679 from the E-alone analyses; and 942 IPI-based relative protein concentrations for E+P, versus 1,054 for E-alone, including 698 that were quantified in both E+P and E-alone analyses.
S Se er ru um m p pr ro ot te ei in n c co on nc ce en nt tr ra at ti io on n r ra at ti io os s Protein concentration ratios were further filtered and curated by using stringent standards (see Methods) for protein identification, including a requirement that a protein is quantified in at least two of the 10 IPAS experiments leading to a focus on 378 proteins (IPIs), all but 10 of which were quantified for both E+P and E-alone. A remarkable 169 (44.7%) of these showed evidence (nominal P < 0.05) of change from baseline to 1-year with E+P or E-alone, or with both. For E+P, 371 proteins were quantified under these quality standards, of which 132 (35.6%) had P < 0.05 as compared with 18.6 expected by chance, and 98 (26.4%) had FDRs <0.05 compared with 94 for E-alone. Of these, 84 had FDR <0.05 for both preparations. Table S1 in Additional file 1 shows estimated 1-year-to-baseline concentration log ratios for all 378 proteins ranked according to the minimum of P values for change with E+P or change with E-alone. Significance levels (P values) are also given for a test of equality of the E+P and E-alone ratios. Table 2 lists proteins for which strong evidence (FDR < 0.01) exists of changed concentration with E+P or with E-alone, according to biologic pathways that were found to be associated with E-alone use in [20]. Nominal P values and FDRs for change also are provided. Of note, five proteins involved in the insulin growth factor pathway are represented in Table 2. Five of these proteins have 1.25-fold or greater changes in their concentrations with E-alone or E+P treatment or both. Protein NOV homologue (NOV) and insulin-like growth factor 1 (IGF1) were both decreased with E-alone and E+P. Insulin-like growth factor-binding protein 1 (IGFBP1) level was increased with both E-alone and E+P. Strong evidence of the E+P effect exists on each of blood coagulation and inflammation, metabolism, osteogenesis, complement and immune response, and cell adhesion. Moreover, the changes (base 2 logarithm of 1-year-to-baseline concentration ratio in Table 2) are mostly quantitatively very similar between E+P and E-alone, attesting to major effects of conjugated equine estrogens on the serum proteome. Table 3 presents the differences in quantitative ratios for E+P minus E-alone that were nominally significantly different (P < 0.05) from each other. For this analysis, we tested the 368 proteins meeting our identification and quantification criteria that were common to both E+P and E-alone. Twenty-six proteins were identified with nominal P values <0.05 for differential change between E+P and E-alone. The list includes proteins involved in insulin growth factor binding and inflammation. Insulin growth factor-binding proteins (IGFBPs) may be affected differently by E+P compared with E-alone. Specifically, three proteins show nominally statistically significant differences with E+P compared with E-alone. IGFBP1 is increased (log 2 ratio, 1.27) with E-alone, but the increase is mitigated with E+P (log 2 ratio, 0.528). IGFBP4 is increased with E-alone (log 2 ratio, 0.511), but not with E+P (log ratio, 0.179). NOV, also known as IGFBP9, is decreased with E-alone (log 2 ratio, -0.344) and decreased to a greater extent with E+P (log 2 ratio, -0.759).

P Pr ro ot te ei in n--s se et t a an na al ly ys se es s
In addition to the protein classifications presented earlier, the 368 proteins quantified for both E+P and E-alone were subjected to protein set analysis. In total, 41 KEGG human disease pathways were represented by at least two proteins in this group. Each protein has been quantified in at least T Ta ab bl le e 1 1 B Ba as se el li in ne e c ch ha ar ra ac ct te er ri is st ti ic cs s a am mo on ng g w wo om me en n i in nc cl lu ud de ed d i in n h ho or rm mo on ne e--t th he er ra ap py y p pr ro ot te eo om mi ic cs s p pr ro oj je ec ct t ( (n n = = 5 50 0 f fo or r E E+ +P P a an nd d f fo or r E E--a al lo on ne e t tr ri ia al ls s) ) E+P E-alone two IPAS experiments in both E+P and E-alone. Table 4 indicates pathways that show a baseline versus 1-year difference for E+P and E-Alone at FDR < 0.05 by using a regularized Hotelling T 2 test. We also tested the equality of log-concentration ratios between the two regimens for proteins in these pathways by using the same test statistic. Two pathways had FDR < 0.05 (Table 5). The gonadotropinreleasing hormone (GnRH) signaling pathway, known to be regulated by estrogen, was represented by two proteins (matrix metalloproteinase 2 (MMP2) and phospholipase A 2 (PLA2G1B)), and a pathway associated with bladder cancer was represented by three proteins (MMP2, thrombospondin 1 (THBS1), and vascular endothelial growth factor C (VEGFC)). Both pathways had nominal P values of 0.002, with corresponding FDRs of 0.041. MMP2, a collagenase with the ability to break down extracellular matrix proteins, was common to both pathways, and substantially explains the difference between the two regimens for these pathways.

E EL LI IS SA A--b ba as se ed d p pr ro ot te ei in n a as ss sa ay ys s i in n a an n i in nd de ep pe en nd de en nt t s se et t o of f s su ub bj je ec ct ts s
The 26 proteins with nominal P < 0.05 for differential change between E+P and E-alone based on IPAS mass spectrometry findings each had FDR > 0.3, so that many of these may be attributable to chance. We sought to determine whether concordant changes in these proteins can be demonstrated in an independent set of subjects and with independent methods. Figure 1 shows 1-year-to-baseline concentration log-ratios (95% confidence intervals (CIs)) from IPAS mass spectrometry data along with corresponding values from ELISA evaluation of 1-year-to-baseline ratios for an independent set of 50 women selected from the active-treatment group in the E+P trial. Corresponding IPAS and ELISA information also is provided for E-alone. We observed concordance of IPAS and ELISA data between the two sets of subjects for six of eight proteins assayed. The lack of replication for ceruloplasmin (CP) and ICAM1 may be due to multiple comparison effects in the discovery component or other factors, notably distinct epitope targets by ELISA assays compared with quantified peptides by mass spectrometry.  independent set of samples. ELISA assays of six of the eight proteins showed changes concordant with the mass spectrometry data. The correlation of initial concentration ratios by mass spectrometry with ELISA ratios from an independent set of samples supports the reliability of the protein changes observed. Our previous report on E-alone [20] provided a detailed discussion of the proteins that changed after treatment with conjugated equine estrogens, in which 19% of proteins were changed after 1 year of treatment. Findings for 10 proteins were confirmed and validated by ELISA assays. Proteins altered with E-alone therapy had relevance to processes such as coagulation, inflammation, growth factors, osteogenesis, metabolism, and cell adhesion, among others. The most striking feature of the present E+P analysis is the similarity in these quantitative proteomic changes when medroxyprogesterone acetate is added to the daily conjugated equine estrogen. For changes with E+P, 98 proteins had FDR < 0.05 compared with 94 proteins for E-alone. Of these, 84 proteins had FDR < 0.05 for both preparations, and corresponding intensity ratios tended to be quite similar between the two regimens for most of these proteins. Hence, our prior discussion [20] of proteins and pathways that were changed after E-alone is largely applicable to the E+P hormone preparation as well. The 1 year of aging between the baseline and 1-year blood-sample collection could have some influence on the serum proteome, but any such influence should be absent for the comparison of E+P versus E-alone changes, because age-related changes would apply equally for the two regimens.

D Di is sc cu us ss si io on n
When we specifically sought proteins for which the change with E+P differed from that for E-alone, a number of potential proteins (Table 3)    be ruled out as an explanation for any particular protein.
To check whether these suggested differences could be attributable to differences in the E+P and E-alone study cohorts (Table 1) we repeated the Table 3 analyses with the mean age and mean BMI at baseline in each pool as adjustment factors. The log intensity ratio differences were not appreciably affected by this adjustment, although P values tended to become less significant because of reduction in 'degrees of freedom' for the moderated t tests. Interestingly, after this adjustment, the FDR for NOTCH2 decreased to 0.02. None of the other FDRs in this sensitivity analysis were <0.05. Aberrant NOTCH signaling has been implicated in tumorigenesis and has been reported to play an oncogenic role in breast cancer [29,30]. A decrease of NOTCH2 serum levels with E+P, but not E-alone, could be related to alteration of signaling and increased risk of breast cancer with E+P but not with E-alone. A differential change between E+P and E-alone in IGFBP1 was supported by ELISA data in an independent set of 50 subjects for each regimen ( Figure 2) and may provide an important lead to understanding clinical effects that differ between the two preparations, including breast cancer. As elaborated later, Table 3 also contains proteins that are associated with atherogenesis.
First, consider proteins involved in the insulin growth factor-signaling pathway. The overall pattern (Table 3) is a greater increase in insulin growth factor-binding proteins (IGFBP1, IGFBP4) with E alone compared with E+P, whereas the decrease in NOV was relatively greater with E+P. ELISA testing produced trends in these same directions for all three proteins, but only that for IGFBP1 approached statistical significance in the independent set. Collectively, these analyses suggest that progestin may attenuate some of the estrogen-induced increases in IGFbinding proteins.
It has been previously suggested that medroxyprogesterone acetate has only a weak degree of opposition to the estrogeninduced decrease of total IGF-1 (which is primarily of hepatic origin), in agreement with our study findings for IGF-1 levels [31]. However, given reduced levels of IGFbinding proteins, it would be expected that less IGF is bound, possibly increasing the availability of free IGF. IGF-signaling pathway plays a role in cell proliferation, tissue development, and tumorigenesis. The IGF pathway has been linked to colorectal malignancy [32,33], and serum levels of IGF1 have been associated with colon cancer risk [34]. Changes in the IGF pathway with E+P compared with E-alone could potentially explain some of the differences in the clinical outcomes. In particular, IGF-1 is a strong mitogen, and varying levels of free IGF-1 between E+P and E-alone treatment could help explain the increased risk of breast cancer with E+P.
In addition to proteins related to the IGF pathway, several other proteins of biologic interest in the context of E+P versus E-alone effects on the serum proteome are presented in Table 3. Expression of the α 2 -glycoprotein 1 zinc (AZGP1) gene is regulated predominantly androgens and progestins [35,36]. Our data suggest an increase in AZGP1 protein levels with E+P and a decrease with E-alone. AZGP1 has been identified as a potential prognostic marker for earlystage breast cancer and a useful immunohistochemical marker of apocrine cell differentiation in human breast tissue [37]. Increased levels of circulating AZGP1 in E+P compared with E-alone may be associated with increased risk of breast cancer in the former group. Circulating levels of extracellular matrix proteins (collagen α-1 chain (COL1A1), lumican (LUM), and extracellular matrix protein 1 (ECM1) may also be differentially affected by E+P compared with Ealone, whereas MMP2, a metalloproteinase that breaks down COL1A1, may be increased with E+P compared with E. The extracellular matrix plays a variety of physiological roles, many of which are related to cancer, including tumor invasion. Changes in the extracellular matrix with E+P compared with E-alone could also help explain the differences in cancer risks associated with these treatments.
Several proteins listed in Table 3 have been linked to atherogenesis, and thus may suggest avenues for exploring mechanisms underlying a more-substantial early increase in CHD risk with E+P than with E-alone in the randomized trials. For example, matrix metalloproteinases (for example, MMP2) are thought to participate in atherogenic inflammation [38]. The role of innate immunity in atherogenesis is less well established; PGLYRP1 participates in recognition of bacteria by neutrophils, but is independently  [39]. The difference in PGLYRP1 concentration among women taking E+P versus E-alone suggests a possible mechanistic link.
Other extracellular matrix proteins (Table 3) may also shed light on the relation between E+P and CHD events. One such protein, lumican (LUM), contributes to variation in proteoglycan composition of arterial intima by location within the human vasculature. Enhanced deposition of lumican has been observed in the intima of the atherosclerosis-prone internal carotid artery compared with the internal thoracic artery, a relatively atherosclerosis-resistant vessel [40]. The relation between serum and tissue proteoglycan levels, the impact of proteoglycan composition on plaque stability, and the clinical significance of lumican all remain to be determined.
Our proteomic comparisons (Table 3) may also provide insight into the greater elevation in venous thromboembolic event risk when progestin was added to conjugated estrogens. For example, coagulation factor V (F5) binds to multimerin 1 (MMFN1), with high affinity for storage in human platelet granules, and may modulate thrombosis [41].
Some important considerations exist in assessing the effects of estrogens and progestins broadly on the serum proteome. The preparations considered here are conjugated equine estrogen and medroxyprogesterone acetate. Further studies would be needed to determine whether the changes reported here also arise for the other estrogen (for example, 17βestradiol) and progesterone (norethisterone acetate or levonorgestrel) treatments. Related to this, these substances are taken orally, and the first-pass hepatic metabolism of oral estrogens is known to stimulate a wide variety of proteins, synthesized in the liver. Of the 378 proteins reported here, 73 are included in the liver-specific gene set listed in Hsiao and colleagues [42] (Table S1 in Additional file 1). For example, of 66 significantly upregulated proteins (FDR < 0.05), 35 are in the liver-specific gene path, as were 28 of 71 for E-alone.
Of the proteins emphasized in the preceding discussion, IGFBP1, AZGP1, and F5, but not others, are part of the liverspecific list. Given that transdermal estrogen, which is being increasingly used in clinical practice to treat menopausal symptoms, bypasses the liver, these proteins may not be affected when estrogen is administered transdermally.

C Co on nc cl lu us si io on ns s
In summary, E+P, like E-alone, has a profound effect on the serum proteome and affects multiple pathways that are relevant to observed clinical effects on cancer, cardiovascular disease, and fractures, among others. The addition of 2.5 mg/d medroxyprogesterone acetate to 0.625 mg/d conjugated equine estrogen may have an impact on the IGF pathway proteins and may affect circulating levels of extracellular matrix proteins (for example, MMP2) of potential relevance to the less-favorable E+P effects, compared with those for E-alone, on breast cancer, and CHD. Similarly the addition of medroxyprogesterone acetate may also augment the effects of conjugated estrogens on coagulation factors (for example, factor V), of potential relevance to a relatively greater elevation in venous thromboembolism with E+P. These and other leads from our proteomic study will benefit from further testing in women who experienced major clinical outcomes and in matched controls from the WHI hormone therapy trials, to evaluate more directly the potential of these protein-concentration changes to contribute to a biologic explanation for observed trial-outcome patterns.
C Co om mp pe et ti in ng g i in nt te er re es st ts s The authors declare that they have no competing interests.
A Au ut th ho or rs s' ' c co on nt tr ri ib bu ut ti io on ns s SJP, SMH, CK, JR, RDJ, JEM, JH, SL, LM, and RLP participated in drafting the manuscript. Data acquisition was performed by HW and HK. Data were analyzed and interpreted by SJP, SMH, LA, LC, SP, HK, QZ, MM, PW, and RLP. Immunoassays were performed by TBB and MMJ. SMH and RLP were responsible for the study design. Statistical analysis was performed by AA, LC, MM, PW, and RLP.
A Ad dd di it ti io on na al l f fi il le es s The following additional files for this article are available online: Additional file 1 contains Table S1, which shows year 1 to baseline log-transformed concentration ratios after estrogen plus progestin (E+P) or estrogen (E-Alone) exposure for all 378 quantified proteins.
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 comprising WHI investigators that included NHLBI representatives.