Next Article in Journal
Use of Clopidogrel and Proton Pump Inhibitors Alone or in Combinations in Persons with Diabetes in Denmark; Potential for CYP2C19 Genotype-Guided Drug Therapy
Next Article in Special Issue
Pre-Diagnostic Circulating Metabolites and Colorectal Cancer Risk in the Cancer Prevention Study-II Nutrition Cohort
Previous Article in Journal
Comparative Effects of Di-(2-ethylhexyl)phthalate and Di-(2-ethylhexyl)terephthalate Metabolites on Thyroid Receptors: In Vitro and In Silico Studies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Metabolomics Analysis of Postmenopausal Breast Cancer Risk in the Cancer Prevention Study II

1
Division of Cancer Epidemiology and Genetics, National Cancer Institute, Rockville, MD 20895, USA
2
Department of Population Science, American Cancer Society, Atlanta, GA 30303, USA
3
Cancer Control and Population Sciences, Huntsman Cancer Institute, Salt Lake City, UT 84112, USA
4
Department of Nutrition and Integrative Physiology, University of Utah, Salt Lake City, UT 84112, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally.
Metabolites 2021, 11(2), 95; https://doi.org/10.3390/metabo11020095
Submission received: 21 January 2021 / Revised: 2 February 2021 / Accepted: 2 February 2021 / Published: 10 February 2021
(This article belongs to the Special Issue Metabolomics Meets Epidemiology)

Abstract

:
Breast cancer is the most common cancer in women, but its incidence can only be partially explained through established risk factors. Our aim was to use metabolomics to identify novel risk factors for breast cancer and to validate recently reported metabolite-breast cancer findings. We measured levels of 1275 metabolites in prediagnostic serum in a nested case-control study of 782 postmenopausal breast cancer cases and 782 matched controls. Metabolomics analysis was performed by Metabolon Inc using ultra-performance liquid chromatography and a Q-Exactive high resolution/accurate mass spectrometer. Controls were matched by birth date, date of blood draw, and race/ethnicity. Odds ratios (ORs) and 95% confidence intervals (CIs) of breast cancer at the 90th versus 10th percentile (modeled on a continuous basis) of metabolite levels were estimated using conditional logistic regression, with adjustment for age. Twenty-four metabolites were significantly associated with breast cancer risk at a false discovery rate <0.20. For the nine metabolites positively associated with risk, the ORs ranged from 1.75 (95% CI: 1.29–2.36) to 1.45 (95% CI: 1.13–1.85), and for the 15 metabolites inversely associated with risk, ORs ranged from 0.59 (95% CI: 0.43–0.79) to 0.69 (95% CI: 0.55–0.87). These metabolites largely comprised carnitines, glycerolipids, and sex steroid metabolites. Associations for three sex steroid metabolites validated findings from recent studies and the remainder were novel. These findings contribute to growing data on metabolite-breast cancer associations by confirming prior findings and identifying novel leads for future validation efforts.

1. Introduction

Breast cancer is the most frequently diagnosed cancer among women, with 2 million cases globally in 2018, and rates rising rapidly [1]. Since incidence rates have historically been two to three times higher in developed than developing nations, it has been thought that “westernization” and accompanying lifestyle changes underlie much of breast cancer risk [2,3]. By the early 2000s, epidemiologists had identified risk factors—such as shorter durations of breastfeeding—that explain up to half the excess risk within developed nations [4,5]. Explanations for the remaining half, however, have proven elusive.
Several lines of evidence indicate that metabolism may be important to breast cancer etiology, including the established associations of physical inactivity and obesity with breast cancer risk [6,7], and the general role that perturbed metabolism plays as a hallmark of cancer [8,9]. Only recently, however, have studies begun to evaluate metabolism more systematically in relation to breast cancer risk, primarily using metabolomics [10,11,12,13,14,15]. These studies reported that levels of many different metabolites were associated with the risk of breast cancer, including metabolites related to mitochondrial and branched-chain amino acid metabolism [14,15], sex steroid hormones [13,14], and various dietary factors [13]. However, some studies were small (e.g., 621 or fewer cases) [10,11,13,14] and therefore in need of replication, while others examined only a small set of metabolites (e.g., <200) [11,12,15] relative to what is currently possible on some metabolomics platforms.
In the present study, we examined levels of 1256 metabolites in prediagnostic serum in a nested case-control study of 782 postmenopausal breast cancer cases and 782 matched controls in the Cancer Prevention Study-II (CPS-II) Nutrition Cohort. The broad span of metabolites examined allowed us to simultaneously evaluate prior metabolic findings for replication and to search for novel metabolic risk factors that may not have been examined in studies preceding the advent of metabolomics.

2. Results

Study participants were an average (mean) 68 years of age and 98% were white (Table 1). The established risk factors for breast cancer of high BMI, early age at menarche, late age at menopause, late age at first live birth/parity, a history of benign breast disease, and family history of breast cancer were (as expected) more prevalent in cases than in controls.
Twenty-four of the 1275 metabolites were statistically significantly associated with breast cancer at a false discovery rate <0.20 in conditional logistic regression models (one below the Bonferroni threshold of 0.00004) (Figure 1; complete results in Supplementary Table S1). These metabolites were predominantly lipids (n = 14), followed by metabolites of unknown identity (n = 6), cofactors/vitamins (n = 2), amino acids (n = 1), and a xenobiotic (n = 1). No carbohydrates, energy, or peptide metabolites were significantly associated with breast cancer. Nine of the metabolites had a positive association with risk, with odds ratios (ORs) for a 90th versus 10th percentile comparison that ranged from 1.75 (95% confidence interval (CI): 1.29–2.36) for X-24293 to 1.45 (95% CI: 1.13–1.85) for 16alpha-hydroxy DHEA 3-sulfate, and 15 metabolites had an inverse association with risk, with ORs ranging from 0.59 (95% CI: 0.43–0.79) for 3,4-methyleneheptanoylcarnitine to 0.69 (95% CI: 0.55–0.87) for oxalate (Figure 2).
After confirming that results for unconditional models were like those of conditional models (the median absolute difference in OR between models was 0.01), we further examined associations between metabolites and breast cancer among non-hormone-users. We found that five metabolites were associated with breast cancer risk in this group at FDR < 0.2 (one below the Bonferroni threshold). The metabolite X-11795 was associated with breast cancer only in non-hormone-users, while the other four metabolites were associated with risk in both non-hormone-users and in all participants (X-24293, 1-palmitoyl-2-palmitoleoyl-GPC (16:0/16:1), Androstenediol (3beta,17beta) monosulfate (2), and 4-allylphenol sulfate) (Figure 3). The respective breast cancer ORs for these metabolites were 2.58 (95% CI: 1.70–3.90), 2.22 (95% CI: 1.45–3.40), 2.18 (95% CI: 1.40–3.39), 1.93 (95% CI: 1.32–2.82), and 0.45 (95% CI: 0.30–0.67). The heterogeneity by menopausal hormone use (current vs non-current use), was statistically significant for X-11795 (pheterogeneity = 7.34 × 10−6; qheterogeneity = 0.07). In our analysis of estrogen receptor-positive (ER+) cases, no associations met the FDR threshold for statistical significance. Results for the 20 metabolites most significantly associated with breast cancer among non-hormone-users and ER+ breast cancer are shown in Supplementary Tables S2 and S3. In exploratory analyses, associations did not vary by time to diagnosis (<5 years vs ≥5 years) or by cancer stage (localized vs regional/distant) (all q-values > 0.2).
Adjusting for additional covariates in our statistical models had little effect on the magnitude of metabolite-breast cancer associations, suggesting that associations are largely independent of known breast cancer risk factors (Table 2). Fifteen of the 24 associations changed by less than 5% when adjusting for these additional factors. The largest change was −9.7% for X-24293 (ORage-adjusted = 1.75 vs ORmultivariable = 1.58).
16alpha-hydroxy DHEA 3-sulfate and three other androgen metabolites are biochemically related to the hormone DHEA and, in our study, each of these metabolites was highly correlated with DHEA-sulfate, the body’s major reservoir of DHEA. The Pearson correlation of 16alpha-hydroxy DHEA 3-sulfate with DHEA-sulfate, for example, was 0.71. To test whether the associations of these metabolites with breast cancer risk were independent of DHEA, we added DHEA-sulfate to models but found little attenuation of associations. For example, the OR for 16alpha-hydroxy DHEA 3-sulfate was 1.54 (95% CI: 1.09–2.19) after adding DHEA-Sulfate, when compared with 1.45 (95% CI: 1.13–1.85) previously. In contrast, adding 16alpha-hydroxy DHEA 3-sulfate to models did affect the association of DHEA sulfate with breast cancer risk. The OR for DHEA sulfate was 0.91 (95% CI: 0.64–1.30) when 16alpha-hydroxy DHEA 3-sulfate was in the model, but 1.24 (95% CI: 0.97–1.60) when 16alpha-hydroxy DHEA 3-sulfate was not in the model.
In a metabolite correlation heatmap, we identified four distinct metabolite clusters encompassing 17 of the 24 breast cancer-associated metabolites (Figure 4). The first and largest cluster included four unknown metabolites (X-11478, X-21319, X-16944, X-18921) and three carnitines (3,4-methyleneheptanoylcarnitine, linolenoylcarnitine (C18:3), and linolenoylcarnitine (C18:2)). The second cluster comprised the four sex steroid hormones of androstenediol (3beta,17beta) monosulfate (2), androsteroid monosulfate (1), androstenediol (3beta,17beta) disulfate (1), and 16alpha-hydroxy DHEA 3-sulfate. The third cluster consisted of the three lipids of 1-(1-enyl-oleoyl)-GPE (P-18:1), 1-(1-enyl-stearoyl)-GPE (P-18:0), 1-(1-enyl-palmitoyl)-2-linoleoyl-GPC (P-16:0/18:2). The fourth and final cluster included the lipids of 1-palmitoyl-2-palmitoleoyl-GPC (16:0/16:1), 2-palmitoleoyl-GPC (16:1), 1-palmitoyl-2-oleoyl-GPC (16:0/18:1). Within each cluster, the correlations between metabolites ranged from 0.32 to 0.89.

3. Discussion

After controlling for multiple testing, we found that levels of 24 distinct metabolites were associated with the risk of postmenopausal breast cancer. Findings for three of the metabolites—16alpha-hydroxy DHEA 3-sulfate, androstenediol (3beta,17beta) monosulfate (2), and androstenediol (3beta,17beta) disulfate (1)—replicate results from recent studies [13,14], suggesting these associations are reliable and worth further exploring. The other findings consist of novel associations for carnitines, glycerophosphocholines, and glycerophosphoethanolamines that future studies should aim to replicate. Several associations, particularly those for hormonal metabolites, appeared stronger among women not using menopausal hormones at baseline, though the difference was statistically significant only for X-11795. Associations did not appear to vary by ER status. Overall, our study contributes to the growing data on metabolite-breast cancer associations by replicating three prior associations and by providing novel leads for future studies to replicate—findings that may help us to better understand breast cancer’s underlying mechanistic basis.
The three metabolites with replicating associations each share a common feature, namely DHEA as a precursor. While DHEA has previously been linked with breast cancer risk DHEA [16,17,18], the associations we identified were statistically independent of DHEA levels, much like in one of our prior analyses [14]. This suggests that the pathways these metabolites are part of could constitute novel targets for prevention efforts and that knowing levels of these metabolites could at least nominally improve breast cancer risk prediction models. Mechanistically, increased formation of 16alpha-hydroxy DHEA 3-sulfate, relative to other DHEA metabolites, may increase breast cancer risk because it can be further metabolized into 16alpha-hydroxyestrone and another 16-hydroxylation pathway estrogen metabolites [19], which studies suggest are especially carcinogenic [20,21]. The other two metabolites are derived specifically from 4-androstenediol, a DHEA metabolite that is more structurally similar to testosterone than other androstenediols [22]. Possibly, increased formation of 4-androstenediol and its metabolites, relative to other forms of androstenediol, may be particularly related to breast cancer risk. Of the seven 3β-androstenediol metabolites (e.g., 5alpha-androstan-3beta,17alpha-diol disulfate) we examined, none were associated with breast cancer risk, hinting at some specificity of association.
Eighteen other associations were novel to the literature and generally involved metabolites whose biological function is understood incompletely. One cluster included the metabolites of linoleoylcarnitine and linolenoylcarnitine, which are long-chain acyl fatty acid derivative esters of carnitines. Levels of such derivatives in serum can reflect shifts in mitochondrial metabolism of fats and proteins [23,24]. Two other clusters included glycerophosphocholines and glycerophosphoethanolamines, and circulating levels of these may indicate perturbations in the formation and maintenance of cellular membranes [25]. Additionally, cancer metabolism studies implicate perturbed choline metabolism as a potential metabolic hallmark of oncogenesis [26]. The last cluster included the sex steroid hormone metabolites, whose relevance was discussed above. If associations for these metabolites replicate, further targeted analyses will be needed to clarify the mechanistic basis of these relationships. An additional three associations--for threonate, oxalate, and 4-allyphenol sulfate—were examined in previous studies and were not statistically significant [13,14].
In total, five prospective studies on metabolomics and breast cancer preceded our own and they reported that an aggregate of 51 metabolites was associated with risk of breast cancer at the FDR < 0.2 level of statistical significance [10,11,12,13,14,15]. We evaluated 27 of these in relation to breast cancer, however, after control for multiple testing, only the three DHEA metabolites discussed above had replicated associations. Five further metabolites replicated at the nominal level of statistical significance (p < 0.05)—namely threonate, a marker for dietary supplement use [13]; two 4-androstenediol metabolites (androstenediol (3beta,17beta) monosulfate-1 and androstenediol (3beta,17beta) disulfate-2); 3-methyl-2-oxobutyrate, a branched-chain amino acid metabolite; and alpha-hydroxyisovalerate, an alcohol-related metabolite [13]. The remaining 19 metabolites had no discernable association with breast cancer risk in these data. If we include replication at the nominal level of significance, our results suggest an approximate 30% replication rate of prior findings, a less-than-ideal rate that indicates that large studies and/or meta-analyses may be needed in this research area [27,28,29].
Strengths of our study include a large number of cases and controls, prospective collection of serum, the detailed data on breast cancer risk factors, and the use of a reliable and robust metabolomics platform. The large number of metabolites we assessed made it possible to both identify novel associations and to attempt replication of previously reported findings.
Our study also has several limitations. We analyzed blood samples from a single point in time for each participant. Since metabolite levels can change over time, this may cause attenuation of odds ratios [28]. We did not use fasting samples, though our prior data suggest that few metabolites are strongly related (>10% of variance explained) to recent food intake [28]. The sample size among non-hormone-users was relatively small at 319 cases. Our study was observational and so unmeasured lifestyle or biological factors could confound associations. Additionally, though we adjusted for parity, we could not directly adjust for months of breastfeeding, which is an important breast cancer risk factor [4]. Nevertheless, our fully adjusted models differed little from age-adjusted models, suggesting that at least traditional breast cancer risk factors are not major confounders here. Our sample consisted predominately of non-Hispanic white women and results may not generalize to other women. Finally, the novel findings from our study require replication in other cohorts or in a consortium of cohorts, like the Consortium of Metabolomics Studies [29].
In summary, our metabolomics analysis delivered a strong replication of associations between DHEA metabolites and risk of breast cancer and identified novel associations that, if replicated, may provide important clues about breast cancer’s mechanistic basis. The study also provided some limited support for the notion that branched-chain amino acids and alcohol-related metabolites are associated with the risk of breast cancer. Future studies are needed to corroborate these findings and to further expand our understanding of the role of metabolism in breast cancer etiology.

4. Materials and Methods

4.1. Study Population and Design

The CPS-II Nutrition Cohort is a large, well-characterized prospective cohort established in 1992 by the American Cancer Society (ACS) that enrolled 184,185 men and women residing in 21 states with population-based state cancer registries [30]. The CPS-II Nutrition Cohort was a sub-cohort of CPS-II, a prospective study of cancer mortality among roughly 1.2 million adults begun in 1982. From June 1998 through May 2001, 21,963 of these women donated nonfasting blood samples [30]. Samples were aliquoted, and frozen in liquid nitrogen at approximately −130 °C for long-term storage. These women also completed a self-administered questionnaire at blood draw that assessed reproductive and medical history and lifestyle factors.
Our nested case-control analysis includes all 782 of the breast cancer cases (ICD-O C50) that occurred among the postmenopausal women who donated blood from 1998 to 2001. Cases were first self-reported on biennial questionnaires and then verified through medical records (90%) or through state cancer registries (10%) and they occurred a median of five years after blood donation. Stage data were available for 99.6% of cases and 172 cases were at a regional/distant stage. Estrogen receptor status data were available on ~94% of participants and 644 cases were estrogen receptor-positive (ER+). Controls with no history of breast cancer or other cancer (except for nonmelanoma skin cancer) were selected by incidence density sampling and 1:1 matched to cases by birth date (±6 months), race/ethnicity (Caucasian, African American, or other/unknown), and date of blood draw (±6 months). Controls were alive and had no history of cancer as of the date of diagnosis for the matched case. All aspects of the CPS-II Nutrition Cohort were approved by the Emory University (Atlanta, GA, USA) Institutional Review Board.

4.2. Metabolomics Assessment

Metabolomics analysis was performed by Metabolon (Metabolon Inc., Durham, NC, USA), whose platform and analytical process have been previously detailed [31,32]. Briefly, serum samples were treated with methanol and two minutes of vigorous shaking followed by centrifugation to precipitate proteins. The resulting extract was split into four fractions that were analyzed using ultra-performance liquid chromatography and a Q-Exactive high resolution/accurate mass spectrometer interfaced with a heated electrospray ionization source and an Orbitrap mass analyzer operated at 35,000 mass resolution. The fractions were examined under the following conditions: (1) acidic positive ion conditions chromatographically optimized for more hydrophilic compounds, (2) acidic positive ion conditions chromatographically optimized for more hydrophobic compounds, (3) basic negative ion optimized conditions, and (4) negative ionization following elution from a hydrophilic interaction liquid chromatography column. Metabolites were identified through comparison with Metabolon’s library of purified standards or recurrent unknown compounds based on retention time/index, mass to charge ratio, and chromatographic data.
Metabolon quantified the relative concentrations of 1385 metabolites in their analysis. Of these, 110 metabolites (primarily drug metabolites) had no detectable level in ≥90% of samples and we excluded them from analyses. Among the remaining 1275 metabolites, 61% of metabolites had fewer than 5% of their values below the limit of detection. Any values below the limit of detection were assigned the lowest value observed for that metabolite in subsequent analyses. To minimize batch variability and account for non-normal distributions, each metabolite’s value was divided by its daily median then log-transformed for statistical analyses.
Triplicates of 46 samples and duplicates of five samples were used as quality controls to assess the reproducibility of the platform. Across the 1275 metabolites, the median intraclass correlation coefficient (ICC) was 0.93 (25th and 75th percentile, respectively: 0.83, 0.98), indicating a high level of technical reliability comparable to that observed in previous analyses [13,28]. The median within-batch coefficient of variation was 0.13 and the median between-batch coefficient of variation was 0.18.

4.3. Statistical Analysis

For each metabolite, we estimated the age-adjusted odds ratio (OR) and 95% confidence intervals (CIs) of breast cancer risk using conditional logistic regression to account for matching. Odds ratios for metabolites were scaled so that they represent risk at the 90th percentile of metabolite values when compared with risk at the 10th percentile (OR = eβ(X90-X10) where β is the coefficient for the metabolite modeled continuously and X90 and X10 are metabolite values at the 90th and 10th percentiles). Associations were also separately estimated for cases among non-hormone-users and for ER+ cases, as prior studies suggest that metabolite associations may be evident predominantly in non-hormone users and/or for the ER+ cancer subtype [14,15]. For the analysis of non-hormone-users, it was necessary to break matched sets and analyze associations using unconditional logistic regression since cases and controls were not matched on hormone use. The threshold for statistical significance for all analyses was set at a false discovery rate [33,34] of less than 0.20, calculated separately for overall cancer, cancer among non-hormone users, and ER+ breast cancers. In exploratory analyses, we also evaluated whether associations varied by time to diagnosis (<5 years vs ≥5 years), cancer stage (localized vs regional/distant), and menopausal hormone use (current vs non-current use) using the Wald test for heterogeneity. Interactions were declared if p-values were less than the false discovery rate threshold of 0.2.
To clarify whether metabolite-breast cancer associations were related to, or confounded by, established breast cancer risk factors, we also separately modeled associations when further adjusting for body mass index, smoking status, alcohol intake, history of diabetes, menopausal hormone therapy use, age at menarche, type of menopause, age at menopause, age at first live birth and number of births, history of benign breast disease, first-degree family history of breast cancer, and moderate-vigorous intensity physical activity (all modeled as indicated in Table 1). We then calculated the differences in ORs between the original and fully adjusted models as an indicator of the degree to which covariates are implicated in metabolite-breast cancer associations. Models adjusted for fewer covariates were also examined in supplementary analyses.
For the associations of hormonal metabolites with breast cancer risk, we additionally considered whether associations were independent of DHEA by adjusting for DHEA-sulfate, the body’s primary reservoir of DHEA. DHEA is the parent compound for many hormonal metabolites and is known to be associated with breast cancer risk [16,17,18].
Finally, to evaluate whether metabolite associations were independent of one another, we examined metabolite-metabolite correlations in a heatmap, ordered using hierarchical clustering. All statistical analyses were conducted in SAS 9.4 (SAS Institute, Cary, NC, USA) and RStudio version 1.1.453 (RStudio Inc., Boston, MA, USA).

Supplementary Materials

The following are available online at https://www.mdpi.com/2218-1989/11/2/95/s1, Supplementary Table S1: Odds ratios, 95% confidence intervals, p-values, and q-values for postmenopausal breast cancer when comparing the 90th with the 10th percentile levels of metabolites in conditional logistic models adjusted for age at blood draw, Supplementary Table S2: Age-adjusted odds ratios (ORs) and 95% confidence intervals (CIs) among non-hormone users at baseline for postmenopausal breast cancer when comparing the 90th with the 10th percentile levels of metabolites (top 20 metabolites by statistical significance), Supplementary Table S3: Age-adjusted odds ratios (ORs) and 95% confidence intervals (CIs) for estrogen receptor-positive breast cancer when comparing the 90th with the 10th percentile levels of metabolites (top 20 metabolites by statistical significance).

Author Contributions

Conceptualization, S.C.M., J.N.S., V.L.S.; formal analysis, S.C.M., K.M.M., M.C.P.; writing—original draft preparation, S.C.M., K.M.M., V.L.S.; writing—review and editing, all authors. All authors have read and agreed to the published version of the manuscript.

Funding

The American Cancer Society funds the creation, maintenance, and updating of the Cancer Prevention Study-II cohort and the work herein was also supported by the Intramural Research Program of the National Cancer Institute, National Institutes of Health, Department of Health and Human Services.

Institutional Review Board Statement

All aspects of the CPS-II Nutrition Cohort were approved by the Emory University (Atlanta, GA, USA) Institutional Review Board.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Data described in the manuscript and analytic code are not available to protect participant confidentiality and in adherence with institutional policies.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Ferlay, J.; Soerjomataram, I.; Dikshit, R.; Eser, S.; Mathers, C.; Rebelo, M.; Parkin, D.M.; Forman, D.; Bray, F. Cancer incidence and mortality worldwide: Sources, methods and major patterns in GLOBOCAN 2012. Int. J. Cancer 2015, 136, E359–E386. [Google Scholar] [CrossRef] [PubMed]
  2. Ziegler, R.G.; Anderson, W.F.; Gail, M.H. Increasing breast cancer incidence in China: The numbers add up. J. Natl. Cancer Inst. 2008, 100, 1339–1341. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Gomez, S.L.; Noone, A.-M.; Lichtensztajn, D.Y.; Scoppa, S.; Gibson, J.T.; Liu, L.; Morris, C.; Kwong, S.; Fish, K.; Wilkens, L.R.; et al. Cancer incidence trends among Asian American populations in the United States, 1990–2008. J. Natl. Cancer Inst. 2013, 105, 1096–1110. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Collaborative Group on Hormonal Factors in Breast C. Breast cancer and breastfeeding: Collaborative reanalysis of individual data from 47 epidemiological studies in 30 countries, including 50,302 women with breast cancer and 96,973 women without the disease. Lancet 2002, 360, 187–195. [Google Scholar]
  5. Islami, F.; Sauer, A.G.; Miller, K.D.; Siegel, R.L.; Fedewa, S.A.; Jacobs, E.J.; McCullough, M.L.; Patel, A.V.; Ma, J.; Soerjomataram, I.; et al. Proportion and number of cancer cases and deaths attributable to potentially modifiable risk factors in the United States. CA Cancer J. Clin. 2018, 68, 31–54. [Google Scholar] [CrossRef]
  6. Patel, A.V.; Friedenreich, C.M.; Moore, S.C.; Hayes, S.C.; Silver, J.K.; Campbell, K.L.; Winters-Stone, K.; Gerber, L.H.; George, S.M.; Fulton, J.E.; et al. American College of Sports Medicine Roundtable Report on Physical Activity, Sedentary Behavior, and Cancer Prevention and Control. Med. Sci. Sports Exerc. 2019, 51, 2391–2402. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Lauby-Secretan, B.; Scoccianti, C.; Loomis, D.; Grosse, Y.; Bianchini, F.; Straif, K.; International Agency for Research on Cancer Handbook Working Group. Body Fatness and Cancer—Viewpoint of the IARC Working Group. New Engl. J. Med. 2016, 375, 794–798. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Hanahan, D.; Weinberg, R.A. Hallmarks of cancer: The next generation. Cell 2011, 144, 646–674. [Google Scholar] [CrossRef] [Green Version]
  9. Vander Heiden, M.G.; Cantley, L.C. Thompson CB. Understanding the Warburg effect: The metabolic requirements of cell proliferation. Science 2009, 324, 1029–1033. [Google Scholar] [CrossRef] [Green Version]
  10. Lécuyer, L.; Bala, A.V.; Deschasaux, M.; Bouchemal, N.; Triba, M.N.; Vasson, M.-P.; Rossary, A.; Demidem, A.; Galan, P.; Hercberg, S.; et al. NMR metabolomic signatures reveal predictive plasma metabolites associated with long-term risk of developing breast cancer. Int. J. Epidemiol. 2018, 47, 484–494. [Google Scholar] [CrossRef]
  11. Kühn, T.; Floegel, A.; Sookthai, D.; Johnson, T.; Rolle-Kampczyk, U.; Otto, W.; Von Bergen, M.; Boeing, H.; Kaaks, R. Higher plasma levels of lysophosphatidylcholine 18:0 are related to a lower risk of common cancers in a prospective metabolomics study. BMC Med. 2016, 14, 13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Tobias, D.K.; Hazra, A.; Lawler, P.R.; Chandler, P.D.; Chasman, D.I.; Buring, J.E.; Lee, I.-M.; Cheng, S.; Manson, J.E.; Mora, S. Circulating branched-chain amino acids and long-term risk of obesity-related cancers in women. Sci. Rep. 2020, 10, 16534. [Google Scholar] [CrossRef]
  13. Playdon, M.C.; Ziegler, R.G.; Sampson, J.N.; Stolzenberg-Solomon, R.; Thompson, H.J.; Irwin, M.L.; Mayne, S.T.; Hoover, R.N.; Moore, S.C. Nutritional metabolomics and breast cancer risk in a prospective study. Am. J. Clin. Nutr. 2017, 106, 637–649. [Google Scholar] [CrossRef] [PubMed]
  14. Moore, S.C.; Playdon, M.C.; Sampson, J.N.; Hoover, R.N.; Trabert, B.; E Matthews, C.; Ziegler, R.G. A Metabolomics Analysis of Body Mass Index and Postmenopausal Breast Cancer Risk. J. Natl. Cancer Inst. 2018, 110, 588–597. [Google Scholar] [CrossRef] [PubMed]
  15. His, M.; Viallon, V.; Dossus, L.; Gicquiau, A.; Achaintre, D.; Scalbert, A.; Ferrari, P.; Romieu, I.; Onland-Moret, N.C.; Weiderpass, E.; et al. Prospective analysis of circulating metabolites and breast cancer in EPIC. BMC Med. 2019, 17, 178. [Google Scholar] [CrossRef] [PubMed]
  16. Key, T.; Appleby, P.; Barnes, I.; Reeves, G.; Endogenous, H. Breast Cancer Collaborative G. Endogenous sex hormones and breast cancer in postmenopausal women: Reanalysis of nine prospective studies. J. Natl. Cancer Inst. 2002, 94, 606–616. [Google Scholar] [PubMed]
  17. Kaaks, R.; Rinaldi, S.; Key, T.J.; Berrino, F.; Peeters, P.H.M.; Biessy, C.; Dossus, L.; Lukanova, A.; Bingham, S.; Khaw, K.-T.; et al. Postmenopausal serum androgens, oestrogens and breast cancer risk: The European prospective investigation into cancer and nutrition. Endocr. Relat. Cancer 2005, 12, 1071–1082. [Google Scholar] [CrossRef]
  18. Zhang, X.; Tworoger, S.S.; Eliassen, A.H.; Hankinson, S.E. Postmenopausal plasma sex hormone levels and breast cancer risk over 20 years of follow-up. Breast Cancer Res. Treat. 2013, 137, 883–892. [Google Scholar] [CrossRef] [Green Version]
  19. Hampl, R.; Starka, L. Minireview: 16alpha-hydroxylated metabolites of dehydroepiandrosterone and their biological significance. Endocr. Regul. 2000, 34, 161–163. [Google Scholar]
  20. Ziegler, R.G.; Fuhrman, B.J.; Moore, S.C.; Matthews, C.E. Epidemiologic studies of estrogen metabolism and breast cancer. Steroids 2015, 99, 67–75. [Google Scholar] [CrossRef] [Green Version]
  21. Sampson, J.N.; Falk, R.T.; Schairer, C.; Moore, S.C.; Fuhrman, B.J.; Dallal, C.M.; Bauer, D.C.; Dorgan, J.F.; Shu, X.-O.; Zheng, W.; et al. Association of Estrogen Metabolism with Breast Cancer Risk in Different Cohorts of Postmenopausal Women. Cancer Res. 2017, 77, 918–925. [Google Scholar] [CrossRef] [Green Version]
  22. Chen, F.; Knecht, K.; Leu, C.; Rutledge, S.; Scafonas, A.; Gambone, C.; Vogel, R.; Zhang, H.; Kasparcova, V.; Bai, C.; et al. Partial agonist/antagonist properties of androstenedione and 4-androsten-3β,17β-diol. J. Steroid Biochem. Mol. Biol. 2004, 91, 247–257. [Google Scholar] [CrossRef]
  23. Domizio, S.; Romanelli, A.; Brindisino, P.; Puglielli, C.; Conte, E.; Sgarrella, M.C.; Sabatino, G. Glutaric aciduria type II: A case report. Int. J. Immunopathol. Pharmacol. 2005, 18, 805–808. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Violante, S.; Ijlst, L.; Brinke, H.T.; De Almeida, I.T.; Wanders, R.J.A.; Ventura, F.V.; Houten, S.M. Carnitine palmitoyltransferase 2 and carnitine/acylcarnitine translocase are involved in the mitochondrial synthesis and export of acylcarnitines. FASEB J. 2013, 27, 2039–2044. [Google Scholar] [CrossRef] [PubMed]
  25. Karanjia, R.N.; E Crossey, M.M.; Cox, I.J.; Fye, H.K.S.; Njie, R.; Goldin, R.D.; Taylor-Robinson, S.D. Hepatic steatosis and fibrosis: Non-invasive assessment. World J. Gastroenterol. 2016, 22, 9880–9897. [Google Scholar] [CrossRef] [Green Version]
  26. Glunde, K.; Bhujwalla, Z.M.; Ronen, S.M. Choline metabolism in malignant transformation. Nat. Rev. Cancer 2011, 11, 835–848. [Google Scholar] [CrossRef] [Green Version]
  27. Moore, S.C. Metabolomics and breast cancer: Scaling up for robust results. BMC Med. 2020, 18, 18. [Google Scholar] [CrossRef] [PubMed]
  28. Sampson, J.N.; Boca, S.M.; Shu, X.O.; Stolzenberg-Solomon, R.Z.; Matthews, C.E.; Hsing, A.W.; Tan, Y.T.; Ji, B.-T.; Chow, W.-H.; Cai, Q.; et al. Metabolomics in epidemiology: Sources of variability in metabolite measurements and implications. Cancer Epidemiol. Biomark. Prev. 2013, 22, 631–640. [Google Scholar] [CrossRef] [Green Version]
  29. Yu, B.; Zanetti, A.K.; Temprosa, M.; Albanes, D.; Appel, N.; Barrera, C.B.; Ben-Shlomo, Y.; Boerwinkle, E.; Casas, J.P.; Clish, C.; et al. The Consortium of Metabolomics Studies (COMETS): Metabolomics in 47 Prospective Cohort Studies. Am. J. Epidemiol. 2019, 188, 991–1012. [Google Scholar] [CrossRef] [Green Version]
  30. Calle, E.E.; Rodriguez, C.; Jacobs, E.J.; Almon, M.L.; Chao, A.; McCullough, M.L.; Feigelson, H.S.; Thun, M.J. The American Cancer Society Cancer Prevention Study II Nutrition Cohort: Rationale, study design, and baseline characteristics. Cancer 2002, 94, 500–511. [Google Scholar] [CrossRef]
  31. Evans, A.M.; DeHaven, C.D.; Barrett, T.; Mitchell, M.; Milgram, E. Integrated, nontargeted ultrahigh performance liquid chromatography/electrospray ionization tandem mass spectrometry platform for the identification and relative quantification of the small-molecule complement of biological systems. Anal. Chem. 2009, 81, 6656–6667. [Google Scholar] [CrossRef] [PubMed]
  32. Evans, A.M.B.B.; Liu, Q.; Mitchell, M.W.; Robinson, R.J.; Dai, H.; Stewart, S.J.; DeHaven, C.D.; Miller, L.A.D. High Resolution Mass Spectrometry Improves Data Quantity and Quality as Compared to Unit Mass Resolution Mass Spectrometry in High-Throughput Profiling Metabolomics. Metabolomics 2014, 4, 132. [Google Scholar]
  33. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. 1995, 57, 289–300. [Google Scholar] [CrossRef]
  34. Ganna, A.; Lee, D.; Ingelsson, E.; Pawitan, Y. Rediscovery rate estimation for assessing the validation of significant findings in high-throughput studies. Brief Bioinform. 2014, 16, 563–575. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Manhattan plot displaying the p-values for metabolite-breast cancer risk associations from conditional logistic models, separated by chemical class. For each chemical class, the key shows the total number of metabolites and the number with statistically significant associations (false discovery rate <0.20). Models were adjusted for age at blood draw.
Figure 1. Manhattan plot displaying the p-values for metabolite-breast cancer risk associations from conditional logistic models, separated by chemical class. For each chemical class, the key shows the total number of metabolites and the number with statistically significant associations (false discovery rate <0.20). Models were adjusted for age at blood draw.
Metabolites 11 00095 g001
Figure 2. Odds ratios, 95% confidence intervals, p-values, and q-values for postmenopausal breast cancer when comparing the 90th with the 10th percentile levels of metabolites in conditional logistic models adjusted for age at blood draw. Results presented only for metabolites below the false discovery rate threshold of 0.20. Markers indicate the odds ratio, and error bars, the 95% confidence intervals. The q-value is the estimated probability of a false discovery. 1 Abbreviations in order of appearance: 1-palmitoyl-2-palmitoleoyl-GPC (16:0/16:1), 1-palmitoyl-2-oleoyl-GPC (16:0/18:1), Androstenediol (3beta,17beta) monosulfate (2), Androstenediol (3beta,17beta) disulfate (1), 1-(1-enyl-palmitoyl)-2-linoleoyl-GPC (P-16:0/18:2).
Figure 2. Odds ratios, 95% confidence intervals, p-values, and q-values for postmenopausal breast cancer when comparing the 90th with the 10th percentile levels of metabolites in conditional logistic models adjusted for age at blood draw. Results presented only for metabolites below the false discovery rate threshold of 0.20. Markers indicate the odds ratio, and error bars, the 95% confidence intervals. The q-value is the estimated probability of a false discovery. 1 Abbreviations in order of appearance: 1-palmitoyl-2-palmitoleoyl-GPC (16:0/16:1), 1-palmitoyl-2-oleoyl-GPC (16:0/18:1), Androstenediol (3beta,17beta) monosulfate (2), Androstenediol (3beta,17beta) disulfate (1), 1-(1-enyl-palmitoyl)-2-linoleoyl-GPC (P-16:0/18:2).
Metabolites 11 00095 g002
Figure 3. Odds ratios, 95% confidence intervals, p-values, and Q-values among non-hormone users at baseline for postmenopausal breast cancer when comparing the 90th with the 10th percentile levels of metabolites (n = 319 cases). Results were from unconditional logistic models adjusted for age at blood draw and are presented only for metabolites below the false discovery rate threshold of 0.20. Markers indicate the odds ratio, and error bars, the 95% confidence intervals. The Q-value is the estimated probability of a false discovery.
Figure 3. Odds ratios, 95% confidence intervals, p-values, and Q-values among non-hormone users at baseline for postmenopausal breast cancer when comparing the 90th with the 10th percentile levels of metabolites (n = 319 cases). Results were from unconditional logistic models adjusted for age at blood draw and are presented only for metabolites below the false discovery rate threshold of 0.20. Markers indicate the odds ratio, and error bars, the 95% confidence intervals. The Q-value is the estimated probability of a false discovery.
Metabolites 11 00095 g003
Figure 4. Heat map depicting pairwise correlations between the 24 metabolites associated with breast cancer risk. Metabolites were ordered using hierarchical clustering analysis as shown by light grey lines on the top and left sides. The key shows how colors correspond to the different Pearson correlation coefficients. For reference, the light red square for linoleoylcarnitine (C18:2) and 1-(1-enyl-oleoyl)-GPE (P-18:1) represents a correlation of 0.25, the medium red square for 16alpha-hydroxy DHEA 3-sulfate and androstenediol (3beta,17beta) monosulfate (2) represents a correlation of 0.55, and the dark red square for X-21319 and X-18921 represents a correlation of 0.75. The diagonal dark red line is the identify line, i.e., the correlation of 1 when comparing a metabolite against itself. 1 Abbreviations in order of appearance (left to right): 1-palmitoyl-2-oleoyl-GPC (16:0/18:1), 1-palmitoyl-2-palmitoleoyl-GPC (16:0/16:1), 1-(1-enyl-palmitoyl)-2-linoleoyl-GPC (P-16:0/18:2), Androstenediol (3beta,17beta) monosulfate (2), Androstenediol (3beta,17beta) disulfate (1).
Figure 4. Heat map depicting pairwise correlations between the 24 metabolites associated with breast cancer risk. Metabolites were ordered using hierarchical clustering analysis as shown by light grey lines on the top and left sides. The key shows how colors correspond to the different Pearson correlation coefficients. For reference, the light red square for linoleoylcarnitine (C18:2) and 1-(1-enyl-oleoyl)-GPE (P-18:1) represents a correlation of 0.25, the medium red square for 16alpha-hydroxy DHEA 3-sulfate and androstenediol (3beta,17beta) monosulfate (2) represents a correlation of 0.55, and the dark red square for X-21319 and X-18921 represents a correlation of 0.75. The diagonal dark red line is the identify line, i.e., the correlation of 1 when comparing a metabolite against itself. 1 Abbreviations in order of appearance (left to right): 1-palmitoyl-2-oleoyl-GPC (16:0/18:1), 1-palmitoyl-2-palmitoleoyl-GPC (16:0/16:1), 1-(1-enyl-palmitoyl)-2-linoleoyl-GPC (P-16:0/18:2), Androstenediol (3beta,17beta) monosulfate (2), Androstenediol (3beta,17beta) disulfate (1).
Metabolites 11 00095 g004
Table 1. Participant characteristics at blood draw in a nested case-control study in the Cancer Prevention Study-II *,†.
Table 1. Participant characteristics at blood draw in a nested case-control study in the Cancer Prevention Study-II *,†.
Cases (n = 782)Controls (n = 782)p
Age (mean yrs ± SD)68 ± 668 ± 6Matched
Body Mass Index (kg/m2) 0.21
<25399 (51%)426 (55%)
25–29.9246 (32%)244 (31%)
30+132 (17%)109 (14%)
Race/ethnicity Matched
Non-Hispanic white764 (98%)768 (98%)
Other18 (2%)14 (2%)
Smoking status <0.0001
Never384 (49%)462 (59%)
Former363 (46%)281 (36%)
Current25 (3%)33 (4%)
Alcohol intake 0.005
Non-drinker254 (33%)299 (38%)
Current drinker513 (66%)449 (57%)
History of diabetes 0.41
No715 (91%)722 (92%)
Yes42 (5%)35 (5%)
Menopausal hormone therapy use 0.06
Never164 (21%)192 (25%)
Former155 (20%)170 (22%)
Current454 (58%)406 (52%)
Age at menarche (years) 0.61
≤1129 (4%)23 (3%)
12–13545 (70%)542 (69%)
14+196 (25%)207 (27%)
Type of menopause, age at menopause (years) 0.04
Natural and <4528 (4%)40 (5%)
Natural and 45–4991 (12%)111 (14%)
Natural and 50–54311 (40%)268 (34%)
Natural 55+85 (11%)73 (9%)
Oophorectomy or surgery226 (29%)255 (33%)
Drugs/treatment induced0 (0%)0 (0%)
Age at first live birth (years), number of live births 0.03
Nulliparous80 (10%)63 (8%)
<20 and 1+ children37 (5%)58 (7%)
20–29 and 1–2 children193 (25%)183 (23%)
20–29 and 3+ children382 (49%)409 (52%)
30+ and 1+ children86 (11%)66 (8%)
History of benign breast disease 0.30
No528 (68%)545 (70%)
Yes250 (32%)230 (30%)
Family history of breast cancer 0.11
No554 (71%)573 (73%)
Yes175 (22%)148 (19%)
Moderate-vigorous intensity physical activity (hours/week) 0.55
None/week398 (51%)386 (49%)
<1 h/week104 (13%)111 (14%)
1 h/week78 (10%)66 (8%)
2–3 h/week129 (17%)150 (19%)
4+ h/week68 (9%)66 (8%)
* Due to the incidence density sampling, 17 participants are included in both case and control columns. Percent estimates do not include missing data. The numbers of participants missing data were as follows: Body mass index: n = 8; Smoking status: n = 16; Alcohol intake: n = 49; History of diabetes; n = 50; Menopausal hormone therapy use: n = 50; Age at menarche: n = 22; Type of menopause and age at menopause: n = 76; Age at first live birth: n = 7; History of benign breast disease: n = 11; Family history of breast cancer: n = 114; Moderate-vigorous intensity physical activity: n = 8. p values derived from chi-square test.
Table 2. Comparison of odds ratios (OR) and 95% confidence intervals (CIs) of breast cancer risk from age-adjusted versus multivariate models.
Table 2. Comparison of odds ratios (OR) and 95% confidence intervals (CIs) of breast cancer risk from age-adjusted versus multivariate models.
MetaboliteAge-Adjusted OR (95% CI)Multivariate 2 OR (95% CI)Difference in OR, %
X-242931.75 (1.29–2.36)1.58 (1.13–2.22)−9.7%
Cysteine1.74 (1.25–2.42)1.74 (1.21–2.52)0.0%
1-palmitoyl-2-palmitoleoyl-GPC 11.62 (1.24–2.13)1.58 (1.16–2.16)−2.5%
2-palmitoleoyl-GPC (16:1)1.60 (1.17–2.20)1.52 (1.08–2.12)−5.0%
1-palmitoyl-2-oleoyl-GPC (16:0) 11.56 (1.19–2.03)1.52 (1.13–2.05)−2.6%
Androstenediol (3beta,17beta) mon 11.55 (1.21–1.99)1.45 (1.10–1.93)−6.5%
Androsteroid monosulfate (1)1.54 (1.18–2.02)1.50 (1.11–2.03)−2.6%
Androstenediol (3beta,17beta) dis 11.46 (1.13–1.88)1.36 (1.02–1.81)−6.9%
16alpha-hydroxy DHEA 3-sulfate1.45 (1.13–1.85)1.46 (1.11–1.93)0.7%
Oxalate (ethanedioate)0.69 (0.55–0.87)0.66 (0.51–0.86)−4.4%
3-hydroxybutyrylcarnitine (2)0.69 (0.54–0.87)0.73 (0.56–0.95)5.8%
X-137290.69 (0.53–0.88)0.72 (0.55–0.95)4.4%
1-(1-enyl-oleoyl)-GPE (P-18:1)0.67 (0.52–0.86)0.69 (0.51–0.92)3.0%
Threonate0.65 (0.50–0.84)0.61 (0.46–0.81)−6.2%
1-(1-enyl-palmitoyl)-2-linoleoy 10.65 (0.50–0.84)0.60 (0.45–0.81)−7.7%
X-213190.64 (0.49–0.83)0.64 (0.48–0.86)0.0%
Linolenoylcarnitine (C18:3)0.64 (0.49–0.83)0.68 (0.51–0.91)6.3%
1-(1-enyl-stearoyl)-GPE (P-18:0)0.64 (0.50–0.83)0.69 (0.51–0.92)7.8%
Linoleoylcarnitine (C18:2)0.64 (0.50–0.83)0.67 (0.50–0.90)4.7%
X-114780.63 (0.49–0.83)0.61 (0.46–0.82)−3.2%
X-169440.60 (0.46–0.77)0.58 (0.44–0.76)−3.3%
X-189210.60 (0.46–0.79)0.58 (0.43–0.77)−3.3%
4-allylphenol sulfate0.58 (0.45–0.75)0.60 (0.45–0.80)3.5%
3,4-methyleneheptanoylcarnitine0.58 (0.43–0.79)0.58 (0.42–0.81)0.0%
1: Abbreviations in order of appearance: 1-palmitoyl-2-palmitoleoyl-GPC (16:0/16:1), 1-palmitoyl-2-oleoyl-GPC (16:0/18:1), Androstenediol (3beta,17beta) monosulfate (2), Androstenediol (3beta,17beta) disulfate (1), 1-(1-enyl-palmitoyl)-2-linoleoyl-GPC (P-16:0/18:2). 2: Multivariate ORs were estimated with conditional logistic regression adjusted for age (years), body mass index (<25 kg/m2, 25.0–30 kg/m2, ≥30 kg/m2, unknown), smoking status (never, former, current), alcohol intake (non-drinker, current drinker, unknown), history of diabetes (no, yes, unknown), menopausal hormone therapy use (never, former, current, unknown), age at menarche (≤12 years, 12–13 years, ≥14 years, unknown), type of menopause and age at menopause (natural and <45 years, natural and 45–49 years, natural and 50–54 years, natural and ≥55 years, bilateral oophorectomy/surgery, drugs/radiation, hysterectomy, or unknown), age at first live birth and number of live births (nulliparous, age ≤19 years and ≥1 live births, age 20–29 years with 1 or 2 live births, age 20–29 with ≥3 live births, age 30+ with ≥1 live births, unknown), history of benign breast disease (no, yes, unknown), first degree family history of breast cancer (no, yes, unknown), and moderate-vigorous intensity physical activity (none, <1 h/week, 1 h/week, 2 h/week, 3 h/week, ≥4 h/week, unknown).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Moore, S.C.; Mazzilli, K.M.; Sampson, J.N.; Matthews, C.E.; Carter, B.D.; Playdon, M.C.; Wang, Y.; Stevens, V.L. A Metabolomics Analysis of Postmenopausal Breast Cancer Risk in the Cancer Prevention Study II. Metabolites 2021, 11, 95. https://doi.org/10.3390/metabo11020095

AMA Style

Moore SC, Mazzilli KM, Sampson JN, Matthews CE, Carter BD, Playdon MC, Wang Y, Stevens VL. A Metabolomics Analysis of Postmenopausal Breast Cancer Risk in the Cancer Prevention Study II. Metabolites. 2021; 11(2):95. https://doi.org/10.3390/metabo11020095

Chicago/Turabian Style

Moore, Steven C., Kaitlyn M. Mazzilli, Joshua N. Sampson, Charles E. Matthews, Brian D. Carter, Mary C. Playdon, Ying Wang, and Victoria L. Stevens. 2021. "A Metabolomics Analysis of Postmenopausal Breast Cancer Risk in the Cancer Prevention Study II" Metabolites 11, no. 2: 95. https://doi.org/10.3390/metabo11020095

APA Style

Moore, S. C., Mazzilli, K. M., Sampson, J. N., Matthews, C. E., Carter, B. D., Playdon, M. C., Wang, Y., & Stevens, V. L. (2021). A Metabolomics Analysis of Postmenopausal Breast Cancer Risk in the Cancer Prevention Study II. Metabolites, 11(2), 95. https://doi.org/10.3390/metabo11020095

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop