Skip to main content

A metabolome-wide case-control study of african american breast cancer patients



Breast cancer survivors face long-term sequelae compared to the general population, suggesting altered metabolic profiles after breast cancer. We used metabolomics approaches to investigate the metabolic differences between breast cancer patients and women in the general population, aiming to elaborate metabolic changes among breast cancer patients and identify potential targets for clinical interventions to mitigate long-term sequelae.


Serum samples were retrieved from 125 breast cancer cases recruited from the Chicago Multiethnic Epidemiologic Breast Cancer Cohort (ChiMEC), and 125 healthy controls selected from Chicago Multiethnic Prevention and Surveillance Study (COMPASS). We used liquid chromatography-high resolution mass spectrometry to obtain untargeted metabolic profiles and partial least squares discriminant analysis (PLS-DA) combined with fold change to select metabolic features associated with breast cancer. Pathway analyses were conducted using Mummichog to identify differentially enriched metabolic pathways among cancer patients. As potential confounders we included age, marital status, tobacco smoking, alcohol drinking, type 2 diabetes, and area deprivation index in our model. Random effects of residence for intercept was also included in the model. We further conducted subgroup analysis by treatment timing (chemotherapy/radiotherapy/surgery), lymph node status, and cancer stages.


The entire study participants were African American. The average ages were 57.1 for cases and 58.0 for controls. We extracted 15,829 features in total, among which 507 features were eventually selected by our criteria. Pathway enrichment analysis of these 507 features identified three differentially enriched metabolic pathways related to prostaglandin, leukotriene, and glycerophospholipid. The three pathways demonstrated inconsistent patterns. Metabolic features in the prostaglandin and leukotriene pathways exhibited increased abundances among cancer patients. In contrast, metabolic intensity in the glycerolphospholipid pathway was deregulated among cancer patients. Subgroup analysis yielded consistent results. However, changes in these pathways were strengthened when only using cases with positive lymph nodes, and attenuated when only using cases with stage I disease.


Breast cancer in African American women is associated with increase in serum metabolites involved in prostaglandin and leukotriene pathways, but with decrease in serum metabolites in glycerolphospholipid pathway. Positive lymph nodes and advanced cancer stage may strengthen changes in these pathways.

Peer Review reports


Breast cancer survivorship has been improved since 1990s due to advances in early detection and cancer therapy [1]. Over 90% of all female breast cancer patients survived 5 years after their initial diagnosis [2]. Currently, more than 3.8 million women in the US are estimated to have a history of breast cancer, comprising the largest group of cancer survivors [1]. Meanwhile, breast cancer patients generally report poorer health compared to the general population even after successful treatment [3, 4]. Higher risks for cardiovascular disease (CVD) [5,6,7], pulmonary disease [8], fatigue [9], chronic pain [10], and cognitive decline [11] are documented among breast cancer survivors. Studies have suggested that CVD is the largest single cause of death among breast cancer survivors, exceeding cancer-related causes [12].

The long-term sequelae of breast cancer are usually attributed to the side effects of cancer therapies, because radiotherapy, chemotherapy, and trastuzumab therapy are able to induce cardiotoxicity and have been linked to higher risks for cardiovascular and pulmonary diseases [13,14,15,16,17,18,19,20,21,22,23]. However, not all sequelae can be explained by cancer therapies. Genetic susceptibility [24, 25] and shared risk factors including age, obesity, and physical inactivity [26,27,28,29,30], also contribute to the development of long-term sequelae after breast cancer diagnosis. Among all potential risk factors, biological changes induced by breast cancer among survivors should not be overlooked.

Given that these long-term sequelae appear to manifest nearly 5 years after initial diagnosis of breast cancer [5, 7], there exists a potential window period for interventions mitigating disease burdens. Meanwhile, minority groups such as African American face a disproportionately large burden of both breast cancer and chronic diseases [31]. Therefore, understanding the cause and etiology of these sequelae among breast cancer patients will be crucial for the development of such medical interventions, especially among minority groups.

As disturbances in metabolic activities underlie most diseases, the study of the metabolome can provide important insight into the etiology of breast cancer sequelae as well as offer the potential to identify metabolic pathways for clinical interventions [32]. More recently, metabolomics technologies give us the ability to measure thousands of metabolites in biological samples, assisting researchers in the investigation of metabolic changes. Metabolomics has demonstrated an emerging and promising role in the diagnosis and prognosis of chronic diseases and clinical interventions [33].

Within this context, we aim to investigate the metabolic profile of breast cancer patients soon after diagnosis and identify potential biological pathways using the state-of-the-art metabolomics technologies among a case-control study with 125 breast cancer cases and 125 controls that were frequency-matched by age. All the participants were African American women residing in Chicago, offering the opportunity to mitigate the larger burdens of both breast cancer and chronic disease in this group [31].


Study population

The Chicago Multiethnic Epidemiologic Breast Cancer Cohort (ChiMEC) was initiated as a hospital-based case-control study to facilitate research on the effects of high-penetrance susceptibility genes, common genetic variants, and environmental risk factors for breast cancer [34,35,36,37]. Breast cancer cases were followed for survival, disease recurrence, and other outcomes to form the ChiMEC cohort, with a current sample size of 5097 patients [38]. Patients diagnosed or treated at the University of Chicago Hospitals were ascertained through the cancer risk clinic and breast center. Clinical, pathological, and treatment data were collected via electronic medical records. Epidemiological risk factor data were collected via a questionnaire. A biobank was established, collecting blood and tumor samples. This analysis randomly selected female African American patients who were ≥ 18 years of age at diagnosis, lived in Chicago, enrolled in the ChiMEC study between 2012 and 2018, had histologically diagnosed non-metastatic invasive breast cancer, and had serum samples available. Stage IV distant metastatic patients were excluded as these patients might have tumor cells in circulation and experience dramatic metabolic changes.

The healthy controls were selected from the Chicago Multiethnic Prevention and Surveillance Study (COMPASS), a large scale, longitudinal cohort study with a current sample size of 7728 participants from 72 of the 77 Chicago community areas [39]. Residents of the greater Chicago area were eligible for COMPASS if they were: (1) 18 or older at the time of enrollment; (2) able to give consent and provide survey data in English or Spanish; (3) willing to provide blood, urine, saliva samples, and access to medical records. Recruitment strategies to increase minority enrollment have included a predominantly minority interviewer team and focus on recruitment in census tracts with minority and diverse populations as the primary sampling unit. Enrollment entails the completion of a 1-hour long survey, consenting for past and future medical records from all sources, the collection of clinical and physical measurement data and the on-site collection of biological samples including blood, urine and saliva. On collection, all biological samples are processed and liquated within 24 h before long-term storage and subsequent analysis. Participants completed an extensive survey providing information on medical history, socioeconomic status, psychosocial variables, lifestyle behaviors, social environment, immune status, use of medical services, medication use, among other covariates. Blood collection occurred at the same time as consent and the in-person interview.

Informed consent

has been obtained from all ChiMEC and COMPASS participants. For the current analysis, 125 African American participants were randomly selected from ChiMEC and COMPASS, respectively, frequency-matched by age. In addition to data collection through questionnaire interview and electronic medical records, we also conducted geocoding analysis based on residence addresses of participants in both studies to collect neighborhood-level characteristics, and we calculated area deprivation index (ADI) [40].

High-resolution metabolomics

The precipitation of proteins were performed using the Ostro 96-well plate by following the manufacturer’s protocol for each serum sample. In brief, 100uL of serum were placed in the well with 100uL of surrogate standard Loratadine. Then 200uL of cold solvent (acetonitrile:formic acid 99:1) were added, followed by gentle mixing before the filtration by manifold processor. The procedure was repeated with 400ul cold solvent (acetonitrile/water/formic acid 3:1:1%). The resulting solution were dried by nitrogen flow and reconstituted in 200uL solvent (acetonitrile/water 1:1) with spiked-in internal standards (Celecobix and 4-Aminobiphenyl).

For testing the serum extraction procedure, a quality control sample was prepared by mixing same volumes of all serum samples. The same extraction protocol was performed on these quality control samples randomly placed in the well plate. The quality control extracts together with the extraction blanks were injected during the sample analysis.

Liquid chromatography mass spectrometry (LC-MS) analysis of the metabolite extracts was performed using an Agilent 6545 Q-TOF and 1290UPLC system controlled by the Agilent Mass Hunter acquisition software. The mass spectrometer was operated in 2 GHz extended dynamic range mode employing precursor ion analysis for relative quantification experiments in positive/negative ion modes. Internal references were used for calibration. Mobile phase A was 0.1% formic acid in H2O and B was 0.1% formic acid in acetonitrile. The samples were loaded onto a 2.1 × 100 mm Agilent Poroshell C18, 1.9 μm column (Agilent Technologies Inc., Santa Clara, CA, USA) and separation performed using an Agilent 1290 Ultra Performance Liquid Chromatography system at a flow rate of 300 µL/ min. The gradient started with 2%B and was increased to 65% in 11 min, then to 95%B in 2 min. The gradient was held at 95%B for 4 min. A post column equilibration time of 4 min was used for all runs. The post time was set to 3 min to re-equilibrate after each run. Source parameters were as follows: gas temp (300 °C), drying gas (11 L/min), nebulizer (35 psi), sheath gas temp (350 °C), sheath gas flow (12 L/min ), VCap (3000 V), and fragmentor (145 V). The column compartment was held at 35 °C. Data was collected for relative quantification using a scan speed of 3 MS spectra per second. Pooled quality control samples were injected across every 10 sample injections. Each extract was separated and quantitated only once.

The surrogate standard was used to track the recovery after extraction and reconstitution. Each sample was used within ± 15% recovery. This protocols were widely used in metabolomics studies [41,42,43]. As suggested by the literature, the protocols used in this study offered the greatest extraction performance, though the survivability of certain metabolites cannot be evaluated [41].

Statistical analysis

Metabolic features that were present in at least 80% of one group and > 50% of all samples were filtered and maintained in following analyses. After filtering, missing values in these features were imputed by one-half of the lowest signal detected for that feature across all samples as suggested by prior studies [44]. Log2-transformed feature intensities were used for analyses. To control for confounding from different aspects, we used residuals of intensities derived from linear mixed effects regression models against potential confounders including age, marital status, tobacco smoking, alcohol drinking, type 2 diabetes, and area deprivation index of residential address [40]. To further eliminate potential confounding arising from spatial variations, we included residential zipcode as random intercepts in the mixed effects models.

Partial least squares discriminative analysis (PLS-DA) was employed to identify differential features between cases and controls. PLS-DA is a supervised multivariate analytical approach for dimensionality reduction that maximizes covariance between cancer status and feature intensities [45]. A Variable Importance in Projection (VIP) score is generated for each feature in the PLS-DA model. The VIP score estimates the importance of each feature in the model and is commonly used for feature selection. Features with a VIP score > 2 in PLS-DA and fold change (FC) > 2 were considered clinically importance. Moreover, to evaluate the performance of selected features, we conducted 10-fold cross-validation tests utilizing the support vector machine and calculated the classification accuracy of the selected features. All feature selection approaches were implemented with the R package mixOmics v6.3.1.

Pathway analysis

The enriched metabolic pathways differential between controls and cases were identified using Mummichog (v. 1.0.10) [46]. All features were included in the pathway enrichment analysis, while the selected features were denoted as significant. Mummichog is a novel and reliable algorithm for pathway enrichment analysis designed specifically for high-resolution LC-MS. This algorithm has been proven to be valid and reflect real biological activity [47,48,49,50]. To reduce false positive match rate, the Mummichog algorithm requires all annotated metabolites to present in at least their primary adduct (M + H or M-H for positive and negative mode, respectively). We used a P value of 0.05 as the threshold. Only enriched pathways with at least 3 overlapping metabolites were kept for further interpretation.

Subgroup analysis

Clinical characteristics related to breast cancer had the potential to modify metabolic activities among patients. To evaluate the influence of these variables, we conducted subgroup analysis based on selected clinic characteristics. Specifically, we repeated aforementioned metabolomics analysis using cases stratified by treatment, lymph node, and cancer stage. To eliminate the influence of cancer treatment on metabolomics in breast cancer cases, we included only cases without any treatment before serum collection in the analysis. To examine the impact of cancer extent on metabolomics, we also conducted subgroup analysis based on cancer stage and lymph node. In summary, we have a total of seven scenarios in the analysis of this study: (1) all breast cancer cases as described above; (2) cases without any cancer treatment (chemotherapy, radiotherapy, and surgery) before serum collection; (3) cases with one of treatment modalities (chemotherapy, radiotherapy, and surgery) before serum collection; (4) cases with negative lymph nodes; (5) cases with positive lymph nodes; (6) cases with stage I cancer; and (7) cases with stages II or III cancer.

Study protocols of both studies (ChiMEC and COMPASS) were approved by the Institutional Review Board at the University of Chicago. All participant recruitment, human sample experiment, and data analysis in this study were carried out in accordance with Helsinki Declaration guidelines.


The average age was 57.1 years for controls and 58.0 years for cases (Table 1). There were more current tobacco smokers (40.8% vs. 13.6%) and alcohol drinkers (61.6% vs. 44.8%) among controls compared to cases. The prevalence of type 2 diabetes were comparable between these two groups (13.6% vs. 11.2%). Both cases and controls lived in similar communities, as indicated by zip codes, in south side of Chicago. The areas that cases and controls lived had similar ranks of area deprivation index (median value: 71 vs. 66). Compared to all ChiMEC population [38], the breast cancer cases in this study demonstrated consistent distributions of demographic characteristics because of the random selection process, as among the original ChiMEC African American patients the average age was 59.6 years and 15.8% reported current smoking. Similarly, the demographic characteristics of healthy controls in this study also exhibited no substantial difference from the original COMPASS cohort, as the average age was 53.7 years, 47.8% reported current smoking, and 10.1% had type 2 diabetes in the original population.

Table 1 Distribution of selected characteristics among participants

Among breast cancer patients in this study, the median time from diagnosis to serum collection was 72 days. There were 49, 58, and 18 patients having breast stage I, II, and III breast cancer, respectively. More than half cancer patients had lymph node tumor negative disease (57.6%). Almost all patients received breast surgery (mastectomy or lumpectomy), and 25 of them received surgery before serum collection (median 98 days). Most patients (62.4%) received chemotherapy, and 46 of them received chemotherapy before serum collection (median 158 days after initiation of chemotherapy). There were 65.6% receiving radiotherapy, and only 9 of them before serum collection. There were 57 (45.6%) patients receiving any treatment (chemotherapy, radiotherapy, and surgery) before serum collection, while the rest 68 patients had sera collected after date of diagnosis and before any cancer treatment.

In total, we detected 15,829 features, among which 12,023 remained after filtering for missing values. When we compared controls with all cases, 518 unique metabolic features were selected using PLS-DA after adjusting for aforementioned confounders, while setting the VIP scores > 2 (Table 2; Fig. 1, Scenario 1). The classification error rate obtained from the 10-fold cross-validation was 4.8%, suggesting the set of discriminatory features effectively separate the cases and controls. Among the 518 features, 507 (97.9%) features had FC > 2 and were eventually used for enriched pathway. In our subgroup analysis, the classification error rates ranged from 2.9 to 8.4%, and the final numbers of selected metabolic features from 478 to 612. The scenario where breast cancer cases had positive lymph node had the lowest classification error rate and the largest number of selected metabolic features. In contrast, the scenario where cases with stage I cancer had the smallest number of selected metabolic features. See supplemental Table S1 for VIP scores and FC of each feature in each scenario in this study.

Table 2 Model performance for feature selection
Fig. 1
figure 1

Identification of metabolic features associated with breast cancer. The positive fold change (log2) indicates higher feature intensity among cases. (1) Scenario 1, all breast cancer cases (n = 125); (2) Scenario 2, cases without any treatments (chemotherapy, radiotherapy, and surgery) before serum collection (n = 68); (3) Scenario 3, cases with any treatments (chemotherapy, radiotherapy, and surgery) before serum collection (n = 57); (4) Scenario 4, cases with negative lymph nodes (n = 72); (5) Scenario 5, cases with positive lymph nodes (n = 52); (6) Scenario 6, cases with stage I cancer (n = 49); (7) Scenario 7, cases with stage II/III cancer (n = 76)

We evaluated whether metabolic features selected in this study were enriched in specific metabolic pathways using Mummichog (Table 3). When using all breast cancer cases (Scenario 1), enriched pathway analysis indicated that three metabolic pathways were differentially enriched between cases and controls with a P-value < 0.05 and had at least 3 overlapped metabolites when we compared breast cancer cases with healthy controls. The three metabolic pathways were prostaglandin formation from arachidonate (8 overlap metabolites), leukotriene metabolism (8 overlap metabolites), and glycerophospholipid metabolism (6 overlap metabolites). The first two pathways, prostaglandin formation from arachidonate and leukotriene metabolism, indicate an inflammatory response. The last pathways is lipid-related metabolic pathway, which is closely related to essential biological structures and functions. Tentative annotation results of metabolites in each differentially enriched pathways and FCs of these features are present in Table 4. Notably, the three pathways demonstrated different patterns. In prostaglandin and leukotriene pathways, FCs of all significant metabolic features were positive, meaning that these features in the two pathways exhibited increased abundance among breast cancer patients. In contrast, in the glycerophospholipid pathway, all significant metabolic features had negative FC and thus showed a decreased intensity among breast cancer patients.

Table 3 Enriched metabolic pathways associated with breast cancer status
Table 4 Tentative match of the serum metabolic features associated with breast cancer to the metabolites within mummichog enriched pathways

The association directions of features in the three pathways were all consistent across all scenarios (Table 4). In the subgroup analysis using data from cases with sera collected after breast diagnosis and before any cancer treatment, we found that the three pathways were still significant or marginal significant (Table 3, Scenario 2). The directions and FC magnitudes were also consistent in overall and the subgroup analysis (Table 4, Scenarios 1 and 2). Interestingly, we found slightly larger number of overlap metabolic features (Table 3) and generally larger effect sizes (Table 4) in prostaglandin and leukotriene pathways when comparing the subgroup analysis of lymph node positive cases (Scenarios 5) to the subgroup analysis of lymph node negative cases (Scenarios 4), while there were similar effect sizes in the glycerophospholipid pathway between the two scenarios. Similarly results were found when comparing the subgroup analysis of stage II/III cancer cases (Tables 3 and 4, Scenario 7) to the subgroup analysis of stage I cancer cases (Scenarios 6). Except for the aforementioned three metabolic pathways, no other significant pathways were observed in the subgroup analyses.


In this study, we used untargeted high-resolution metabolomics approach to investigate the metabolic changes soon after diagnosis among breast cancer patients compared with healthy controls. With more than 15,000 metabolic features detected from serum of the study samples, we found that metabolic pathways related to inflammatory reactions and lipid metabolism were differentially enriched among breast cancer patients, providing clues for the long-term sequelae manifested after breast cancer diagnosis. Specifically, two inflammatory reactions pathways, prostaglandin formation from arachidonate and leukotriene metabolism, were differentially enriched between controls and cases. Metabolic features in these two pathways showed an increased abundance among breast cancer patients. The changes in these pathways were strengthened when only considering breast cancer cases with positive lymph nodes or stage II/III cancer. By contrast, in the enriched lipid-related pathway, glycerophospholipid metabolism, the intensity of metabolic feature were mainly downregulated among cancer patients.

Glycerophospholipids are the most abundant lipids in virtually all mammalian membranes. Dysfunction in glycerophospholipid metabolism is among the most prominent metabolic alterations in cancer [51], as cancer cells need to continually regulate glycerophospholipids for membrane production, energy acquisition, and molecular signaling [52]. Therefore, it is likely that the observed changes in glycerophospholipid metabolism in this study are the direct consequence of breast cancer and thus suggest disruptions among breast cancer patients. The altered lipid metabolism further disrupts cholesterol homeostasis [53, 54], which has been linked to the progression of various chronic disease [55,56,57,58]. Given the crucial role of glycerophospholipid in cell membranes and surfactant, it is natural and intuitive that dysfunction of glycerophospholipid metabolism results in a series of chronic diseases [59]. However, the long-term impact of the disruptions in glycerophospholipid metabolism is under-investigated, warranting more studies.

In addition to serving as key structural components of membranes, glycerophospholipids can also provide arachidonic acid, the precursor of prostaglandins and leukotrienes [60]. As a key inflammatory intermediate, arachidonic acid is released from cell membranes controlled by glycerophospholipids and then converted to eicosanoids through pathways including: the lipoxygenase (LOX) pathway, where arachidonic acids is dioxygenated to produce hydroperoxyeicosatetraenoic acid (HPETE) and then further converted to leukotrienes; and the cyclooxygenase (COX) pathway that produces prostaglandin [61]. Prostaglandins and leukotrienes are two major pro-inflammatory mediators, and their metabolic pathways were observed to be differentially enriched among breast cancer patients in this study.

The most common type of prostaglandins is prostaglandin E2 (PGE2) whose link to CVD has been extensively documented in the literature [62]. PGE2 is involved in a variety of biological functions in the cardiovascular system, including vascular tone regulation, cardiac remodeling, and cardiac inflammation [62]. While clinical studies are still lacking, evidence from animal experiments reveals that PGE2 is associated with higher risks of heart failure, reperfusion injury, arrhythmias, hypertension, hypertensive heart disease, and atherosclerosis [62]. In this study, we observed that PGE2 was increased among breast cancer patients.

Leukotrienes have been linked to increased risks for CVD and other diseases by prior studies [63]. The overproduction of leukotrienes is a major cause of inflammation [64]. Leukotriene D4 (LTD4) is one of leukotrienes that retains biological activity. In our study, we observed that LTD4 was increased intensity among breast cancer patients, revealing another potential mechanistic pathway of the ling-term sequelae after breast cancer diagnosis. Recently, given the close relationship between CVD and leukotrienes, there is an emerging interest develop cardiovascular drugs based on leukotriene modifiers, which are expected to mitigate stroke, myocardial infarction, and atherosclerosis [65].

In the subgroup analysis, though consistent with our main findings, the changes in these differentially enriched pathways were somewhat strengthened in the scenario where cases had positive lymph nodes or stage II/III cancer, while attenuated in scenarios where cases had negative lymph nodes or stage I cancer. The observation of more extent of cancer with stronger associations with the aforementioned three pathways is reasonable and anticipated, as ancillary lymph nodes are considered the most important route for breast cancer to spread to blood and then other organs.

The study has several strengths. First, to our knowledge, this study is the first one that employ metabolomics approaches to analyze the difference in metabolism between breast cancer cases and healthy controls. Findings from this study identified potential metabolic target for clinical intervention that mitigates long-term sequelae among breast cancer patients. Second, substantial efforts have been made to ensure the homogeneity of the study population with an aim to preclude the influence from external factors. Specifically, the cases and controls were matched by age and all of them were African Americans, so that variations arising from genetics and age could be reduced. Moreover, all participants were recruited from similar communities in south side of Chicago and random effects of spatial variations were further included in the regression model. Therefore, impacts of environmental factors, such as air pollution and neighborhood disadvantages, could be minimized. All these approaches were expected to increases the validity and biological reliability of our results, though we cannot completely rule out external influences.

Our study also has some limitations that warrant caution in interpretation. First, the breast cancer samples are heterogeneous regarding time of serum collection and cancer stage. Chemotherapy and radiotherapy can have short and long-term influence on metabolic profile in circulation, while breast surgery can have short-term impact on inflammatory markers. In the subgroup analysis using samples collected after diagnosis and before any of these treatment, the main study findings remained, suggesting that the differentially regulated pathways identified in this study may be due to intrinsic cancer metabolism. The heterogeneity by cancer stage, as discussed above, suggests that these differentially regulated pathways depend on extent of breast cancer. Second, without further metabolite identification using tandem MS, we could only tentatively annotate the extracted features and enriched pathways using computational approaches. This is an inherent challenge in untargeted metabolomics studies. We could not rule out false matches that influenced our interpretation. Observations from the current study only provide helpful, but not conclusive, suggestions for metabolic changes related to breast cancer. Future studies are warranted to improve the identification of metabolites using either tandem MS or internal standards. Third, the study population were all African Americans, therefore, generalizability to other populations might be weakened. Particularly, genetic variations across racial/ethnic groups could prevent us generalizing conclusions from African Americans to other groups. However, given African American females face higher disease burdens of both breast cancer and chronic diseases, our results can also provide valuable insights regarding public health and clinical interventions.

In summary, we applied high-resolution metabolomics to identify perturbations in the serum metabolome associated with breast cancer. We observed metabolic pathways consistent with inflammatory reactions and lipid metabolism that may contribute to long-term sequelae among breast cancer patients. Particularly, we observed differentially enriched metabolic pathways of prostaglandins and leukotrienes with increased abundances of metabolic features among cancer patients. Changes in these two pathways appear to be strengthened by positive lymph nodes, and weakened at negative lymph nodes. We also observed decreased metabolic intensities in the pathway of glycerolphospholipid among breast cancer patients. The findings from this study may provide insights to identify clinical intervention targets that mitigate long-term sequelae among breast cancer patients.

Data Availability

The data that support the findings of this study are available on request from the corresponding author, D.H. The data are not publicly available due to their containing information that could compromise the privacy of research participants.


  1. Miller KD, Nogueira L, Mariotto AB, Rowland JH, Yabroff KR, Alfano CM, et al. Cancer treatment and survivorship statistics, 2019. Cancer J Clin. 2019;69(5):363–85.

    Article  Google Scholar 

  2. Siegel RL, Miller KD, Fuchs HE, Jemal A, Cancer statistics. 2021. CA: a cancer journal for clinicians. 2021;71(1):7–33.

  3. Yabroff KR, Lawrence WF, Clauser S, Davis WW, Brown ML. Burden of illness in cancer survivors: findings from a population-based national sample. J Natl Cancer Inst. 2004;96(17):1322–30.

    Article  PubMed  Google Scholar 

  4. Weaver KE, Forsythe LP, Reeve BB, Alfano CM, Rodriguez JL, Sabatino SA, et al. Mental and Physical Health–Related Quality of Life among US Cancer Survivors: Population estimates from the 2010 National Health interview SurveyHealth-Related Quality of Life among US Cancer Survivors. Cancer Epidemiol Biomarkers Prev. 2012;21(11):2108–17.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Park N-J, Chang Y, Bender C, Conley Y, Chlebowski RT, Van Londen G, et al. Cardiovascular disease and mortality after breast cancer in postmenopausal women: results from the women’s Health Initiative. PLoS ONE. 2017;12(9):e0184174.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Patnaik JL, Byers T, DiGuiseppi C, Dabelea D, Denberg TD. Cardiovascular disease competes with breast cancer as the leading cause of death for older females diagnosed with breast cancer: a retrospective cohort study. Breast Cancer Res. 2011;13(3):1–9.

    Article  Google Scholar 

  7. Bradshaw PT, Stevens J, Khankari N, Teitelbaum SL, Neugut AI, Gammon MD. Cardiovascular disease mortality among breast cancer survivors. Epidemiol (Cambridge Mass). 2016;27(1):6.

    Article  Google Scholar 

  8. Carver JR, Shapiro CL, Ng A, Jacobs L, Schwartz C, Virgo KS, et al. American Society of Clinical Oncology clinical evidence review on the ongoing care of adult cancer survivors: cardiac and pulmonary late effects. J Clin Oncol. 2007;25(25):3991–4008.

    Article  CAS  PubMed  Google Scholar 

  9. Berger AM, Gerber LH, Mayer DK. Cancer-related fatigue: implications for breast cancer survivors. Cancer. 2012;118(S8):2261–9.

    Article  PubMed  Google Scholar 

  10. Peuckmann V, Ekholm O, Rasmussen NK, Groenvold M, Christiansen P, Møller S, et al. Chronic pain and other sequelae in long-term breast cancer survivors: nationwide survey in Denmark. Eur J Pain. 2009;13(5):478–85.

    Article  CAS  PubMed  Google Scholar 

  11. Vannorsdall TD. Cognitive changes related to cancer therapy. Med Clin. 2017;101(6):1115–34.

    Google Scholar 

  12. Gulati M, Mulvagh SL. The connection between the breast and heart in a woman: breast cancer and cardiovascular disease. Clin Cardiol. 2018;41(2):253–7.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Hozumi Y, Suemasu K, Takei H, Aihara T, Takehara M, Saito T, et al. The effect of exemestane, anastrozole, and tamoxifen on lipid profiles in japanese postmenopausal early breast cancer patients: final results of National Surgical Adjuvant Study BC 04, the TEAM Japan sub-study. Ann Oncol. 2011;22(8):1777–82.

    Article  CAS  PubMed  Google Scholar 

  14. Amir E, Seruga B, Niraula S, Carlsson L, Ocaña A. Toxicity of adjuvant endocrine therapy in postmenopausal breast cancer patients: a systematic review and meta-analysis. J Natl Cancer Inst. 2011;103(17):1299–309.

    Article  CAS  PubMed  Google Scholar 

  15. Cuppone F, Bria E, Verma S, Pritchard KI, Gandhi S, Carlini P, et al. Do adjuvant aromatase inhibitors increase the cardiovascular risk in postmenopausal women with early breast cancer: Meta-analysis of randomized trials. Cancer: Interdisciplinary International Journal of the American Cancer Society. 2008;112(2):260–7.

    Article  CAS  Google Scholar 

  16. Albini A, Pennesi G, Donatelli F, Cammarota R, De Flora S, Noonan DM. Cardiotoxicity of anticancer drugs: the need for cardio-oncology and cardio-oncological prevention. J Natl Cancer Inst. 2010;102(1):14–25.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Bowles EJA, Wellman R, Feigelson HS, Onitilo AA, Freedman AN, Delate T, et al. Risk of heart failure in breast cancer patients after anthracycline and trastuzumab treatment: a retrospective cohort study. J Natl Cancer Inst. 2012;104(17):1293–305.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Darby SC, Ewertz M, McGale P, Bennet AM, Blom-Goldman U, Brønnum D, et al. Risk of ischemic heart disease in women after radiotherapy for breast cancer. N Engl J Med. 2013;368(11):987–98.

    Article  CAS  PubMed  Google Scholar 

  19. Du XL, Xia R, Liu CC, Cormier JN, Xing Y, Hardy D, et al. Cardiac toxicity associated with anthracycline-containing chemotherapy in older women with breast cancer. Cancer. 2009;115(22):5296–308.

    Article  PubMed  Google Scholar 

  20. Hooning MJ, Botma A, Aleman BM, Baaijens MH, Bartelink H, Klijn JG, et al. Long-term risk of cardiovascular disease in 10-year survivors of breast cancer. J Natl Cancer Inst. 2007;99(5):365–75.

    Article  PubMed  Google Scholar 

  21. Khouri MG, Douglas PS, Mackey JR, Martin M, Scott JM, Scherrer-Crosbie M, et al. Cancer therapy–induced cardiac toxicity in early breast cancer: addressing the unresolved issues. Circulation. 2012;126(23):2749–63.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Pinder MC, Duan Z, Goodwin JS, Hortobagyi GN, Giordano SH. Congestive heart failure in older women treated with adjuvant anthracycline chemotherapy for breast cancer. J Clin Oncol. 2007;25(25):3808–15.

    Article  CAS  PubMed  Google Scholar 

  23. Giordano SH, Hortobagyi GN, Kau S-WC, Theriault RL, Bondy ML. Breast cancer treatment guidelines in older women. J Clin Oncol. 2005;23(4):783–91.

    Article  PubMed  Google Scholar 

  24. van Westerop L, Arts-de Jong M, Hoogerbrugge N, De Hullu J, Maas A. Cardiovascular risk of BRCA1/2 mutation carriers: a review. Maturitas. 2016;91:135–9.

    Article  PubMed  Google Scholar 

  25. Sajjad M, Fradley M, Sun W, Kim J, Zhao X, Pal T, et al. An exploratory study to determine whether BRCA1 and BRCA2 mutation carriers have higher risk of cardiac toxicity. Genes. 2017;8(2):59.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Eliassen AH, Colditz GA, Rosner B, Willett WC, Hankinson SE. Adult weight change and risk of postmenopausal breast cancer. JAMA. 2006;296(2):193–201.

    Article  CAS  PubMed  Google Scholar 

  27. McTigue KM, Chang YF, Eaton C, Garcia L, Johnson KC, Lewis CE, et al. Severe obesity, heart disease, and death among white, african american, and hispanic postmenopausal women. Obesity. 2014;22(3):801–10.

    Article  PubMed  Google Scholar 

  28. Nichols HB, Trentham-Dietz A, Egan KM, Titus-Ernstoff L, Holmes MD, Bersch AJ, et al. Body mass index before and after breast cancer diagnosis: associations with all-cause, breast cancer, and cardiovascular disease mortality. Cancer Epidemiol Prev Biomarkers. 2009;18(5):1403–9.

    Article  Google Scholar 

  29. Hildebrand JS, Gapstur SM, Campbell PT, Gaudet MM, Patel AV. Recreational physical activity and leisure-time sitting in relation to postmenopausal breast cancer risk. Cancer Epidemiol Prev Biomarkers. 2013;22(10):1906–12.

    Article  Google Scholar 

  30. Reddigan JI, Ardern CI, Riddell MC, Kuk JL. Relation of physical activity to cardiovascular disease mortality and the influence of cardiometabolic risk factors. Am J Cardiol. 2011;108(10):1426–31.

    Article  PubMed  Google Scholar 

  31. Mehta LS, Watson KE, Barac A, Beckie TM, Bittner V, Cruz-Flores S, et al. Cardiovascular disease and breast cancer: where these entities intersect: a scientific statement from the American Heart Association. Circulation. 2018;137(8):e30–e66.

    Article  PubMed  PubMed Central  Google Scholar 

  32. McGarrah RW, Crown SB, Zhang G-F, Shah SH, Newgard CB. Cardiovascular metabolomics. Circul Res. 2018;122(9):1238–58.

    Article  CAS  Google Scholar 

  33. Ussher JR, Elmariah S, Gerszten RE, Dyck JR. The emerging role of metabolomics in the diagnosis and prognosis of cardiovascular disease. J Am Coll Cardiol. 2016;68(25):2850–70.

    Article  CAS  PubMed  Google Scholar 

  34. Nanda R, Schumm LP, Cummings S, Fackenthal JD, Sveen L, Ademuyiwa F, et al. Genetic testing in an ethnically diverse cohort of high-risk women: a comparative analysis of BRCA1 and BRCA2 mutations in american families of european and african ancestry. JAMA. 2005;294(15):1925–33.

    Article  CAS  PubMed  Google Scholar 

  35. Huo D, Senie RT, Daly M, Buys SS, Cummings S, Ogutha J, et al. Prediction of BRCA mutations using the BRCAPRO model in clinic-based african american, hispanic, and other minority families in the United States. J Clin Oncol. 2009;27(8):1184.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Huo D, Feng Y, Haddad S, Zheng Y, Yao S, Han Y-J, et al. Genome-wide association studies in women of african ancestry identified 3q26. 21 as a novel susceptibility locus for oestrogen receptor negative breast cancer. Hum Mol Genet. 2016;25(21):4835–46.

    CAS  PubMed  PubMed Central  Google Scholar 

  37. Huo D, Zheng Y, Ogundiran TO, Adebamowo C, Nathanson KL, Domchek SM, et al. Evaluation of 19 susceptibility loci of breast cancer in women of african ancestry. Carcinogenesis. 2012;33(4):835–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Shang L, Hattori M, Fleming G, Jaskowiak N, Hedeker D, Olopade OI, et al. Impact of post-diagnosis weight change on survival outcomes in Black and White breast cancer patients. Breast Cancer Res. 2021;23(1):1–15.

    Article  Google Scholar 

  39. Aschebrook-Kilfoy B, Kibriya MG, Jasmine F, Stepniak L, Gopalakrishnan R, Craver A, et al. Cohort profile: the ChicagO Multiethnic Prevention and Surveillance Study (COMPASS). BMJ open. 2020;10(9):e038481.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Kind AJH, Buckingham WR. Making Neighborhood-Disadvantage Metrics Accessible - The Neighborhood Atlas. N Engl J Med. 2018;378(26):2456–8.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Tulipani S, Mora-Cubillos X, Jáuregui O, Llorach R, García-Fuentes E, Tinahones FJ, et al. New and vintage solutions to enhance the plasma metabolome coverage by LC-ESI-MS untargeted metabolomics: the not-so-simple process of method performance evaluation. Anal Chem. 2015;87(5):2639–47.

    Article  CAS  PubMed  Google Scholar 

  42. Graham SF, Chevallier OP, Elliott CT, Hölscher C, Johnston J, McGuinness B, et al. Untargeted metabolomic analysis of human plasma indicates differentially affected polyamine and L-arginine metabolism in mild cognitive impairment subjects converting to Alzheimer’s disease. PLoS ONE. 2015;10(3):e0119452.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Anesi A, Berding K, Clarke G, Stanton C, Cryan JF, Caplice N, et al. Metabolomic workflow for the Accurate and High-Throughput Exploration of the Pathways of Tryptophan, Tyrosine, Phenylalanine, and branched-chain amino acids in human biofluids. J Proteome Res. 2022;21(5):1262–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Yan Q, Liew Z, Uppal K, Cui X, Ling C, Heck JE, et al. Maternal serum metabolome and traffic-related air pollution exposure in pregnancy. Environ Int. 2019;130:104872.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Wold S, Sjöström M, Eriksson L. PLS-regression: a basic tool of chemometrics. Chemometr Intell Lab Syst. 2001;58(2):109–30.

    Article  CAS  Google Scholar 

  46. Li S, Park Y, Duraisingham S, Strobel FH, Khan N, Soltow QA, et al. Predicting network activity from high throughput metabolomics. PLoS Comput Biol. 2013;9(7):e1003123.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Li S, Dunlop AL, Jones DP, Corwin EJ. High-resolution metabolomics: review of the field and implications for nursing science and the study of preterm birth. Biol Res Nurs. 2016;18(1):12–22.

    Article  PubMed  Google Scholar 

  48. Li S, Todor A, Luo R. Blood transcriptomics and metabolomics for personalized medicine. Comput Struct Biotechnol J. 2016;14:1–7.

    Article  PubMed  Google Scholar 

  49. Li S, Sullivan NL, Rouphael N, Yu T, Banton S, Maddur MS, et al. Metabolic phenotypes of response to vaccination in humans. Cell. 2017;169(5):862–77. e17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Uppal K, Walker DI, Liu K, Li S, Go Y-M, Jones DP. Computational metabolomics: a framework for the million metabolome. Chem Res Toxicol. 2016;29(12):1956–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Bian X, Liu R, Meng Y, Xing D, Xu D, Lu Z. Lipid metabolism and cancer. Journal of Experimental Medicine. 2021;218(1).

  52. Dolce V, Rita Cappello A, Lappano R, Maggiolini M. Glycerophospholipid synthesis as a novel drug target against cancer. Curr Mol Pharmacol. 2011;4(3):167–75.

    Article  CAS  PubMed  Google Scholar 

  53. Sankaram MB, Thompson TE. Interaction of cholesterol with various glycerophospholipids and sphingomyelin. Biochemistry. 1990;29(47):10670–5.

    Article  CAS  PubMed  Google Scholar 

  54. Ohvo-Rekilä H, Ramstedt B, Leppimäki P, Slotte JP. Cholesterol interactions with phospholipids in membranes. Prog Lipid Res. 2002;41(1):66–97.

    Article  PubMed  Google Scholar 

  55. Kelly RS, Virkud Y, Giorgio R, Celedón JC, Weiss ST, Lasky-Su J. Metabolomic profiling of lung function in Costa-Rican children with asthma. Biochimica et Biophysica Acta (BBA)-Molecular Basis of Disease. 2017;1863(6):1590–5.

  56. Zhu T, Li S, Wang J, Liu C, Gao L, Zeng Y, et al. Induced sputum metabolomic profiles and oxidative stress are associated with chronic obstructive pulmonary disease (COPD) severity: potential use for predictive, preventive, and personalized medicine. EPMA J. 2020;11(4):645–59.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Deng P, Barney J, Petriello MC, Morris AJ, Wahlang B, Hennig B. Hepatic metabolomics reveals that liver injury increases PCB 126-induced oxidative stress and metabolic dysfunction. Chemosphere. 2019;217:140–9.

    Article  CAS  PubMed  Google Scholar 

  58. Dang VT, Zhong LH, Huang A, Deng A, Werstuck GH. Glycosphingolipids promote pro-atherogenic pathways in the pathogenesis of hyperglycemia-induced accelerated atherosclerosis. Metabolomics. 2018;14(7):1–10.

    Article  CAS  Google Scholar 

  59. Cruickshank-Quinn CI, Jacobson S, Hughes G, Powell RL, Petrache I, Kechris K, et al. Metabolomics and transcriptomics pathway approach reveals outcome-specific perturbations in COPD. Sci Rep. 2018;8(1):1–18.

    Article  CAS  Google Scholar 

  60. Hermansson M, Hokynar K, Somerharju P. Mechanisms of glycerophospholipid homeostasis in mammalian cells. Prog Lipid Res. 2011;50(3):240–57.

    Article  CAS  PubMed  Google Scholar 

  61. Tam VC, Quehenberger O, Oshansky CM, Suen R, Armando AM, Treuting PM, et al. Lipidomic profiling of influenza infection identifies mediators that induce and resolve inflammation. Cell. 2013;154(1):213–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Bryson TD, Harding P. Prostaglandin E2 EP receptors in cardiovascular disease: an update. Biochem Pharmacol. 2022;195:114858.

    Article  CAS  PubMed  Google Scholar 

  63. Poeckel D, Funk CD. The 5-lipoxygenase/leukotriene pathway in preclinical models of cardiovascular disease. Cardiovascular Res. 2010;86(2):243–53.

    Article  CAS  Google Scholar 

  64. Nelson DL, Lehninger AL, Cox MM. Leukotrienes. Lehninger principles of biochemistry. 5th ed. Macmillan; 2008. p. 359.

  65. Funk CD. Leukotriene modifiers as potential therapeutics for cardiovascular disease. Nat Rev Drug Discovery. 2005;4(8):664–72.

    Article  CAS  PubMed  Google Scholar 

Download references


We thank the dedicated ChiMEC and COMPASS field staff and community partners for their support of this work. We also thank the UChicago Medicine for its infrastructure support.


This research was funded by a Chicago Center for Health and the Environment CACHET pilot grant (2P30ES027792) and in part by National Cancer Institute (P20-CA233307).

Author information

Authors and Affiliations



B.A. and D.H. conceived the idea and acquired the funding. J.L. analyzed the data, interpreted the results, wrote the main manuscript text, and prepared the tables and figures. M.G.K. and H.C. conducted the experiments. H.A. helped supervise the project. O.I.O. and C.S.O provided necessary clinical support. All authors reviewed the manuscript and contributed to the discussion to shape the article.

Corresponding authors

Correspondence to Briseis Aschebrook-Kilfoy or Dezheng Huo.

Ethics declarations

Conflict of Interest

The authors declare no conflict of interests.

Ethical approval

Study protocols were approved by the Institutional Review Board at the University of Chicago. All participant recruitment, human sample experiment, and data analysis in this study were carried out in accordance with Helsinki Declaration guidelines. Informed consent has been obtained from all ChiMEC and COMPASS participants.

Consent for publication

Not applicable.

Competing interests

The authors have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Luo, J., Kibriya, M.G., Chen, H. et al. A metabolome-wide case-control study of african american breast cancer patients. BMC Cancer 23, 183 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • breast cancer
  • metabolome
  • metabolomics
  • prostaglandin
  • leukotriene
  • glycerophospholipid