Skip to main content
Download PDF

Lipidome of mammographic breast density in premenopausal women

Breast Cancer Research volume 25, Article number: 121 (2023) Cite this article

Abstract

Background

High mammographic breast density (MBD) is a strong risk factor for breast cancer development, but the biological mechanisms underlying MBD are unclear. Lipids play important roles in cell differentiation, and perturbations in lipid metabolism are implicated in cancer development. Nevertheless, no study has applied untargeted lipidomics to profile the lipidome of MBD. Through this study, our goal is to characterize the lipidome of MBD in premenopausal women.

Methods

Premenopausal women were recruited during their annual screening mammogram at the Washington University School of Medicine in St. Louis, MO. Untargeted lipidomic profiling for 982 lipid species was performed at Metabolon (Durham, NC®), and volumetric measures of MBD (volumetric percent density (VPD), dense volume (DV), and non-dense volume (NDV)) was assessed using Volpara 1.5 (Volpara Health®). We performed multivariable linear regression models to investigate the associations of lipid species with MBD and calculated the covariate-adjusted least square mean of MBD by quartiles of lipid species. MBD measures were log10 transformed, and lipid species were standardized. Linear coefficients of MBD were back-transformed and considered significant if the Bonferroni corrected p-value was < 0.05.

Results

Of the 705 premenopausal women, 72% were non-Hispanic white, and 23% were non-Hispanic black. Mean age, and BMI were 46 years and 30 kg/m2, respectively. Fifty-six lipid species were significantly associated with VPD (52 inversely and 4 positively). The lipid species with positive associations were phosphatidylcholine (PC)(18:1/18:1), lysophosphatidylcholine (LPC)(18:1), lactosylceramide (LCER)(14:0), and phosphatidylinositol (PI)(18:1/18:1). VPD increased across quartiles of PI(18:1/18:1): (Q1 = 7.5%, Q2 = 7.7%, Q3 = 8.4%, Q4 = 9.4%, Bonferroni p-trend = 0.02). The lipid species that were inversely associated with VPD were mostly from the triacylglycerol (N = 43) and diacylglycerol (N = 7) sub-pathways. Lipid species explained some of the variation in VPD. The inclusion of lipid species increased the adjusted R2 from 0.45, for a model that includes known determinants of VPD, to 0.59.

Conclusions

We report novel lipid species that are associated with MBD in premenopausal women. Studies are needed to validate our results and the translational potential.

Background

High mammographic breast density (MBD) is a strong risk factor for breast cancer development [1]. Women with extremely dense breasts have a 4–sixfold increased risk of developing breast cancer compared to women with scattered areas of fibroglanduar density [1, 2]. It is estimated that ~ 29% of breast cancers in premenopausal women could be potentially prevented if women with dense or heterogeneously dense breasts reduced their density to scattered areas of fibroglandular density [3]. Further, MBD is an intermediate phenotype [4] for breast cancer development; hence, understanding the mechanisms underlying high MBD may open up new avenues for breast cancer prevention.

Adiposity is strongly associated with MBD measures. We and others have demonstrated that childhood, early adulthood, and adulthood adiposity, as well as changes in adiposity over the life course, are associated with MBD in premenopausal women [5,6,7]. Adiposity explains the greatest percentage of variability in MBD, approximately 22% in non-Hispanic black women and 26% in non-Hispanic white women [8]. Studies on the associations of adiposity biomarkers with MBD yielded inconsistent results [9,10,11,12]. Thus, there is an urgent need to investigate the associations of novel biomarkers of adiposity with MBD.

The lipidome is the complete lipid profile within a tissue [13]. Lipidomics, the comprehensive analysis of lipid molecules, is essential to understanding lipid biology. Recent technological advances have allowed the quantification of lipid species unbiasedly with increasing accuracy. Lipids play essential roles in cell structure, generation of cell membranes, signaling, and fuel/storage [14]. Lipids also function as second messengers and as hormones and are responsible for various cell functions. Changes in lipid metabolism have been shown to impact cell growth, proliferation, differentiation, and motility [15], and perturbations in lipid metabolism are implicated in metabolic diseases [16, 17]. Applying lipidomics to MBD should, therefore, provide meaningful insights into the biologic mechanisms underlying MBD, but to the best of our knowledge, no study has characterized the lipidome of MBD. Our study aimed to characterize, for the first time, the comprehensive lipidome of MBD in premenopausal women.

Methods

Study population

This study population is comprised of 705 premenopausal women who were recruited during annual screening mammogram at the Joanne Knight Breast Health Center at the Siteman Cancer Center at Washington University School of Medicine (WU) in St. Louis, MO. Women were eligible to participate in the study if they were premenopausal, not pregnant, and were able to comply with study procedures. Women with a history of cancer, breast augmentation (implants or reduction), and currently use or have used selective estrogen receptor modulators in the past six months were excluded [6]. Women were considered premenopausal if they had a regular menstrual cycle in the past 12 months and did not have a history of hormone replacement therapy or bilateral oophorectomy [6]. All participants provided written informed consent, and the study was performed in accordance with the Declaration of Helsinki. We received approval for this study from the WU Institutional Review Board.

Participants completed a questionnaire on behavioral, reproductive, demographic, and clinical characteristics. They also had their height, weight, and body fat percentage measured on the day of their screening mammogram visit [6]. Women were asked to fast on the day of their mammogram prior to providing a blood sample [6]. Blood samples were sent to the Tissue Procurement Core at WU Siteman Cancer Center within 30 min of collection and stored at −80 °C [18].

Lipidomics profiling

Blood samples from each woman were sent to Metabolon (Durham, NC®) for comprehensive quantitative lipidomic profiling. Lipidomic profiling quantified 982 lipid species in 3 lipid super pathways and 14 sub-pathways: phospholipids (phosphatidylcholines (PC), lysophosphatidylcholines (LPC), phosphatidylethanolamines (PE), lysophosphatidylethanolamines (LPE), and phosphatidylinositols (PI)), sphingolipids (ceramides (CER), dihydroceramides (DCER), hexosylceramides (HCER), lactosylceramides (LCER), and sphingomyelins (SM)) and neutral complex lipids (cholesteryl esters (CE), diacylglycerols (DAG), triacylglycerols (TAG), and monoacylglycerols (MAG). Individual lipid species were quantified by taking the ratio of the signal intensity of each target compound to that of its assigned internal standard, then multiplying by the concentration of internal standard added to the sample [19]. Quality control samples, a large pool of human plasma maintained by Metabolon, were included in the assay runs and median relative standard deviation (RSD) percent values were calculated for each run. The overall median RSD across runs was 8%. Lipid sub-pathway concentrations were calculated from the sum of all molecular species within a sub-pathway [19].

Mammographic breast density assessment

We assessed volumetric measures of MBD (volumetric percent density (%) (VPD), dense volume (cm3) (DV), and non-dense volume (cm3) (NDV)) using Volpara 1.5 (Volpara Health®). VPD was calculated by using the maximum fibroglandular volume between the left and right breasts. VPD was calculated by dividing DV by total breast volume and multiplying by 100. VPD was also categorized into: < 3.5%, 3.5% ~ 7.5%, > 7.5% ~ 15.5%, and > 15.5%, corresponding to the MBD groups of almost entirely fatty, scattered areas of fibroglandular density, heterogeneously dense, and extremely dense.

Statistical analysis

Distributions of demographic, reproductive, and clinical factors were summarized by the four VPD categories into means and standard deviations for continuous variables and counts and percentages for categorical variables. One hundred twenty-five lipid species missing in 300 or more of the 705 samples were excluded from the analyses, leaving 857 lipid species. Lipid species with less than 300 missing samples were imputed using the 10-nearest neighbor methods using the R package impute.[20] We excluded 5 women with missing MBD measures, leaving 700 in the analytic sample. Spearman correlation coefficients were calculated across the lipid sub-pathways and lipid species with MBD measures.

We performed covariate-adjusted multivariable linear regression models fitting each MBD outcome with the concentration of lipid sub-pathways and lipid species, in quantitative scale and in quartiles, to determine their associations. MBD measures were log10 transformed for normality. In the quantitative scale, each lipid sub-pathway/species was standardized to have a zero mean and a unit standard deviation. The following covariates were accounted for in covariate-adjusted analyses: age (continuous), age at menarche (continuous), body shape at age 10 (based on Stunkard pictogram), body fat % (continuous), race (non-Hispanic white, non-Hispanic black, other), family history of breast cancer (yes, no), oral contraceptive use (never, less than 1 year, 1–4 years, 5–9 years, more than 10 years), alcohol consumption (never, < 1 drink/week, 1–2 drinks/week, 3–5 drinks/week, and 6 + drinks/week), and parity/age at first birth (nulliparous, 1–2 children & < 25 years, 1–2 children & 25–29 years, 1–2 children & ≥ 30 years, ≥ 3 children & < 25 years, ≥ 3 children & ≥ 25 years). Body fat % was used instead of body mass index (BMI) because it explained a slightly greater proportion of variation in VPD (R2 = 0.45) than BMI (R2 = 0.43). BMI and body fat % were highly correlated (r = 0.88). Some of the covariates included in the analyses had > 1% missingness, including body shape at age 10 (N = 40, 5.7%), and body fat % (N = 25, 3.5%); hence, missing values in all covariates were imputed using multivariate imputation by chain equations method via the R package mice.[21] We report the least square means (LSM) of VPD, DV, and NDV by quartile of lipid sub-pathway and lipid species. We calculated the p-value for trend in the adjusted LSM along quartiles by setting all values within each quartile to the median of that quartile range and operationalizing it as an ordinal variable. Linear coefficients were back-transformed. Residuals from the linear regression models were graphically examined for model goodness-of-fit. The proportional odds ratio model was applied to the four VPD categories with each lipid class as a predictor with the inclusion of the covariates to estimate the adjusted odds ratio with 95% CI. We corrected for multiple testing and the family-wise error rate using false discovery rate (FDR), and Bonferroni method, respectively.

To identify variables that influenced VPD the most, we bootstrapped the original dataset 200 times and applied least absolute shrinkage and selection operator (LASSO) penalized multivariable linear regression to each bootstrapped dataset to fit VPD with consideration of all available variables (the lipid sub-pathways, lipid species, and the covariates). We calculated the frequency of each variable being retained in the model fitted from the 200 bootstrapped datasets as an importance measure of the variables in predicting VPD. Model goodness-of-fit was evaluated based on the adjusted R2 value after sequentially adding lipid species to the model.

Results

The mean age, BMI, and body fat % of study participants were 46 years, 30 kg/m2, and 40.4%, respectively; approximately 71.8% of participants were non-Hispanic white, and 23.1% were non-Hispanic black (Table 1). Almost half of the women in the study (44.9%) had scattered areas of fibroglandular density (VPD of 3.5% ~ 7.5%), and 28.1% had heterogeneously dense breasts.

Table 1 Characteristics of Premenopausal Women Recruited During Annual Screening Mammogram across Volumetric Percent Density Categoriesa
Full size table

Correlations between lipid sub-pathways/species and MBD

Lipid species in the TAG, LCER, LPC, and PC sub-pathways were most strongly correlated with VPD, NDV, and DV (Additional file 1: Fig. 1, and Additional file 2: Table 1). The strongest positive correlations observed with VPD were with LPC(18:1) (r = 0.39, p-value = 4.4E−27); with NDV was TAG54:6-FA20:4 (r = 0.51, p-value = 3.3E−47); and with DV was TAG58:10-FA20:4 (r = 0.14, p-value = 0.0003). The strongest inverse correlations for VPD were TAG54:6-FA20:4 (r = − 0.47, p−value = 8.6E−40); with NDV was LPC(18:1) (r = − 0.44, p-value = 1.9E−34); and with DV was PC(18:1/18:2) (r = −0.10, p-value = 0.008), (Additional file 2: Table 1).

Multivariable linear regression between lipid sub-pathways/species and MBD

At the sub-pathway level, 4 of the 9 lipid sub-pathways (DCER, LCER, DAG, and TAG) were significantly associated with VPD (Bonferroni p value < 0.05) in the covariate-adjusted linear regression analyses (Fig. 1A and Additional file 1: Table 2). Fifty-six lipid species across seven lipid sub-pathways were significantly associated with VPD: TAG (N = 43), DAG (N = 7), PC (N = 2), PI (N = 1), LPC (N = 1), LCER (N = 1), and CE (N = 1). Four of the 56 lipid species were positively associated: LCER(14:0), LPC(18:1), PC(18:1/18:1), and PI(18:1/18:1), while the remaining 52 were inversely associated with VPD (Fig. 2A). Of the 52 lipid species that were inversely associated with VPD, 43 were from the TAG sub-pathway. The TAG species with similar chain lengths that were significantly associated with VPD displayed strong positive correlations ranging from 0.72–0.99 for TAG50 species, 0.58–0.99 for TAG52, and 0.69–0.97 for TAG54. One standard deviation increase in TAG54:6-FA20:4 was associated with a > 10% decrease in VPD, while one standard deviation increase in LCER(14:0) and PI(18:1/18:1) was associated with a > 10% increase in VPD. Three lipid sub-pathways (DCER, DAG, and TAG) were significantly associated with NDV (Fig. 1B and Additional file 1: Table 2). One hundred and seventeen lipid species were significantly associated with NDV, with 113 lipid species from the DAG (N = 14) and TAG (N = 99) sub-pathways showing positive associations, and 4 species (LCER(14:0), LPC(18:1), PC(18:1/18:1), and PI(18:1/18:1)) showing inverse associations (Fig. 2B). Three lipid sub-pathways and 48 lipid species (38.4%) were significantly associated with both VPD and NDV (Fig. 1C and Fig. 2C), and in opposite directions as expected. Of note, one standard deviation in TAG54:6-FA20:4 was associated with a > 10% change in both VPD (decrease) and NDV (increase). No lipid sub-pathway or species was significantly associated with DV.

Fig. 1
figure 1

Covariate-adjusted Associations between Lipid Sub-pathways with VPD and NDV A volcano plot of lipid sub pathways for VPD, B volcano plot lipid sub-pathways for NDV, C Venn diagram of lipid sub-pathways for VPD, DV and NDV Abbreviations: volumetric percent density (VPD), non-dense volume (NDV) dihydroceramide (DCER), lactosylceramide (LCER), diacylglycerol (DAG), triacylglycerol (TAG)

Full size image
Fig. 2
figure 2

Covariate-adjusted Associations between Lipid Species with VPD and NDV A volcano plot of lipid species for VPD, B volcano plot of lipid species for NDV, C Venn diagram of lipid species for VPD, DV and NDV. Abbreviations: volumetric percent density (VPD), non-dense volume (NDV) lysophosphatidylcholine (LPC), phosphatidylinositol (PI), lactosylceramide (LCER), diacylglycerol (DAG), phosphatidylcholine (PC), triacylglycerol (TAG)

Full size image

Multivariable covariate-adjusted least square means of volumetric percent density by quartiles of lipid sub-pathways/species

We categorized the lipid sub-pathways and species into quartiles and then calculated the least squared means (LSM) of VPD and NDV across the quartiles for those that were significant in the multiple linear regression analysis (Figs. 1&2). The four lipid sub-pathways, DCER, LCER, DAG, and TAG, were still associated with VPD after both FDR and Bonferroni correction (Table 2). All the 56 lipid species that were associated with VPD in continuous scale (Fig. 2A) were still significantly associated with VPD in quartiles at an FDR p-value < 0.05 and 48 at a more stringent Bonferroni p-value < 0.05 (Table 3). PI (18:1/18:1) was the lipid species with the strongest positive association with VPD. VPD increased across quartiles of PI(18:1/18:1): (Q1 = 7.5%, Q2 = 7.7%, Q3 = 8.4%, Q4 = 9.4%,FDR p-trend = 8.8E-05; Bonferroni p-trend = 0.02), (Table 3). All the 117 lipid species that were significantly associated with NDV in continuous scale (Fig. 2B) were still significantly associated with NDV at an FDR p-trend < 0.05, but only 101 were associated with NDV at a Bonferroni p-trend < 0.05, (Additional file 1: Table 3). Quartiles of DCER, LCER, DAG, and TAG were associated with NDV (Additional file 1: Table 4). There were no significant associations between lipid sub-pathways/lipid species and DV (Additional file 1: Table 4).

Table 2 Covariate-adjusted Least Square Means of Volumetric Percent Density by Quartiles of the 14 Lipid Sub-pathways
Full size table
Table 3 Covariate-adjusted Least Square Means of Volumetric Percent Density (VPD) by Quartiles of Lipid Species that were Significantly Associated with VPD at a Bonferroni P-value < 0.05 a,b,c
Full size table

Proportional odds model

Additional file 1: Fig. 2 presents the results from the covariate-adjusted proportional odds model investigating the associations between lipid sub-pathways and the four VPD categories. LPC and LCER were significantly positively associated with VPD categories, while DCER, DAG, and TAG were inversely associated.

LASSO regression and bootstrapping

Body fat %, body shape at age 10, parity/age at first birth, and family history of breast cancer were most strongly associated with VPD and were retained in the final lasso-penalized multivariable linear regression model in 100% of the 200 bootstrapped datasets. Some lipid species (e.g., PI(18:1/18:1)—99%, LCER (14:0), and PE(O-16:0/22:6)—98.5%) were selected at higher frequencies than classic covariates, such as age, race and age at menarche, that are associated with VPD. For VPD, the adjusted R2 (a model goodness-of-fit measure) derived from the model containing only the classic covariates was 0.45. A reasonably selected (based on the scree plot of the adjusted R2) model, which included all the covariates and 57 lipid species, rendered an improved adjusted R2 = 0.59 (results not shown).

Discussion

Using untargeted comprehensive lipidomics profiling, we demonstrate, for the first time, the associations of lipid species with MBD in premenopausal women. We observed 56 lipid species that were associated with VPD and 117 lipid species with NDV. Most of these lipid species were from the large TAG sub-pathway. Some of the lipid species explain greater variation in VPD than well-established determinants of MBD and appear to improve model goodness-of-fit based on the adjusted R2 as the covariate-only model rendered an adjusted R2 = 0.45 while additionally including 57 lipid species increased the adjusted R2 to 0.59.

To the best of our knowledge, this is the first study to characterize the lipidome of MBD. Several lipid species were associated with MBD after adjustment for adiposity. Studies have investigated the associations of biomarkers of adiposity, such as adipokines and insulin-like growth factors with, MBD, but the findings are not consistent. The association of most of these biomarkers with MBD, particularly the adipokines, is influenced by BMI. Several studies found a null association between these biomarkers and MBD after adjusting for measures of adiposity [10, 11, 22,23,24,25], although some studies reported an inverse association between leptin and percent density [10, 11, 26]. The association of C-peptide with MBD is similarly attenuated once adjusted for BMI [9, 27]. Similarly, the associations of insulin-like growth factors with MBD are not consistent [12, 28,29,30,31]. Our study, therefore, provides important novel insights on the biological mechanisms underlying MBD.

Although no population-based studies have comprehensively characterized the lipidome of breast cancer, preclinical studies suggest that breast cancer subtypes have unique lipidomic profiles [32,33,34], and a few studies have evaluated the utility of a limited set of lipid species for their diagnostic performance in breast cancer development. One study profiled 110 lipid species in 121 breast cancer cases and 45 healthy controls;19 lipid species distinguished women with triple-negative breast cancer (TNBC) from non-cancer controls, and 5 lipid species distinguished women with TNBC from women with other types of breast cancer [35]. Another study compared lipid profiles in breast tissue of 42 women with breast cancer to 19 healthy controls, reporting 48 significantly differential lipid species. They identified higher levels of lipid species in the LPC sub-pathway but lower levels of lipids in ceramides, DAG, PC, and PE sub-pathways in breast tissues of women with breast cancer compared to healthy controls [36]. They found higher levels of LPC(18:1) in cancer tissue compared to normal tissue [36]. We found a significant positive association between LPC(18:1) and VPD; hence, LPC(18:1) may be a novel biomarker of VPD and breast cancer development. LPC(18:1) has been found to be abundant in MCF-7 cell lines (ER + /PR +) compared to MCF10A cell lines (normal) [32].

Lipids are generated in cells through de novo synthesis or by intake through exogenous sources [15]. Lipid species are a diverse group of compounds and are categorized by their primary function, which is determined by their polar head and are further differentiated by the addition of a hydrocarbon chain [37]. Lipids are functionally important in maintaining cell membranes and fueling the cell [15]. Lipids within the phospholipid and sphingolipid super pathways are a key component in cell membrane development, while neutral complex lipids, particularly DAG and TAG provide energy to fuel the cell [34].

Three of the 4 lipid species (PC(18:1/18:1), LPC(18:1), PI(18:1/18:1)) that were positively associated with VPD were from the phospholipid super pathway and the associations are biologically plausible. They contain at least one monounsaturated fatty acid (MUFA) chain. Higher levels of MUFAs have been reported in breast tumor tissues compared to normal tissues from the same woman [38], and some, but not others, have suggested associations of MUFAs with breast cancer risk [39]. PI(18:1/18:1) promotes cell survival and is responsive to diet and stress [40]. LPCs have an impact on innate immunity and have been shown to possess both anti-and pro-inflammatory properties [41]. Two studies exploring the relationship between circulating triglycerides and MBD observed null or inverse findings.[42, 43]. Our findings on the inverse associations of several TAG species with VPD, therefore, provide new information on triglycerides and MBD. Because of the very strong correlations of TAG species with the same chain lengths, studies are needed to determine which of the TAG species are more amenable to targeting breast cancer prevention.

Strengths/limitations

Our study population is large and well annotated, with detailed information on demographic, reproductive, and anthropometric measures. This allowed us to perform robust analyses with control for confounders. The study population is diverse and mirrors the underlying population attending screening mammogram at our institution; hence, our results are generalizable to the source population which our study participants were recruited from or other populations of similar characteristics [8]. We performed lipidomic profiling on fasting blood samples, which should provide a more reliable measure of the lipid species than using non-fasted samples [44].

There are also some limitations. The study is cross-sectional; therefore, we cannot establish longitudinal trajectories of lipids. Due to the large number of lipid species evaluated, some findings may be due to chance. Nevertheless, we applied the stringent Bonferonni correction; hence, the potential of chance findings are considerably reduced. Our study is limited to premenopausal women; therefore, our results cannot be generalized to postmenopausal women.

Conclusion

Our study offers new insights into the biological mechanisms underlying high MBD in premenopausal women. Additional studies are needed to validate our findings.

Availability of data and materials

De-identified data are available from the corresponding author upon reasonable request.

Abbreviations

MBD:

Mammographic breast density

WU:

Washington university school of medicine

PC:

Phosphatidylcholines

LPC:

Lysophosphatidylcholines

PE:

Phosphatidylethanolamines

LPE:

Lysophosphatidylethanolamines

PI:

Phosphatidylinositols

CER:

Ceramides

DCER:

Dihydroceramides

HCER:

Hexosylceramides

LCER:

Lactosylceramides

SM:

Sphingomyelins

CE:

Cholesteryl esters

DAG:

Diacylglycerols

TAG:

Triacylglycerols

MAG:

Monoacylglycerols

VPD:

Volumetric percent density

DV:

Dense volume

NDV:

Non-dense volume

BMI:

Body mass index

LSM:

Least square means

FDR:

False discovery rate

LASSO:

Least absolute shrinkage and selection operator

TNBC:

Triple-negative breast cancer

MUFA:

Monounsaturated fatty acid

SD:

Standard deviation

CI:

Confidence interval

Q1-4:

Quartiles 1–4

RSD:

Relative standard deviation

References

  1. Boyd NF, Guo H, Martin LJ, Sun L, Stone J, Fishell E, Jong RA, Hislop G, Chiarelli A, Minkin S, et al. Mammographic density and the risk and detection of breast cancer. N Engl J Med. 2007;356(3):227–36.

    Article  CAS  PubMed  Google Scholar 

  2. McCormack VA, dos Santos SI. Breast density and parenchymal patterns as markers of breast cancer risk: a meta-analysis. Cancer Epidemiol Biomarkers Prev. 2006;15(6):1159–69.

    Article  PubMed  Google Scholar 

  3. Engmann NJ, Golmakani MK, Miglioretti DL, Sprague BL, Kerlikowske K. Breast cancer surveillance C: population-attributable risk proportion of clinical risk factors for breast cancer. JAMA Oncol. 2017;3(9):1228–36.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Boyd NF, Rommens JM, Vogt K, Lee V, Hopper JL, Yaffe MJ, Paterson AD. Mammographic breast density as an intermediate phenotype for breast cancer. Lancet Oncol. 2005;6(10):798–808.

    Article  PubMed  Google Scholar 

  5. Alimujiang A, Imm KR, Appleton CM, Colditz GA, Berkey CS, Toriola AT. Adiposity at age 10 and mammographic density among premenopausal women. Cancer Prev Res (Phila). 2018;11(5):287–94.

    Article  PubMed  Google Scholar 

  6. Alimujiang A, Appleton C, Colditz GA, Toriola AT. Adiposity during early adulthood, changes in adiposity during adulthood, attained adiposity, and mammographic density among premenopausal women. Breast Cancer Res Treat. 2017;166(1):197–206.

    Article  PubMed  Google Scholar 

  7. Soguel L, Diorio C. Anthropometric factors, adult weight gain, and mammographic features. Cancer Causes Control. 2016;27(3):333–40.

    Article  PubMed  Google Scholar 

  8. Moore JX, Han Y, Appleton C, Colditz G, Toriola AT: Determinants of Mammographic Breast Density by Race Among a Large Screening Population. JNCI Cancer Spectr 2020, 4(2):pkaa010.

  9. Ahern TP, Hankinson SE, Willett WC, Pollak MN, Eliassen AH, Tamimi RM. Plasma C-peptide, mammographic breast density, and risk of invasive breast cancer. Cancer Epidemiol Biomarkers Prev. 2013;22(10):1786–96.

    Article  CAS  PubMed  Google Scholar 

  10. Woolcott CG, Courneya KS, Boyd NF, Yaffe MJ, McTiernan A, Brant R, Jones CA, Stanczyk FZ, Terry T, Cook LS, et al. Association between sex hormones, glucose homeostasis, adipokines, and inflammatory markers and mammographic density among postmenopausal women. Breast Cancer Res Treat. 2013;139(1):255–65.

    Article  CAS  PubMed  Google Scholar 

  11. Dossus L, Rinaldi S, Biessy C, Hernandez M, Lajous M, Monge A, Ortiz-Panozo E, Yunes E, Lopez-Ridaura R, Torres-Mejia G, et al. Circulating leptin and adiponectin, and breast density in premenopausal Mexican women: the Mexican Teachers’ Cohort. Cancer Causes Control. 2017;28(9):939–46.

    Article  CAS  PubMed  Google Scholar 

  12. Rice MS, Tworoger SS, Rosner BA, Pollak MN, Hankinson SE, Tamimi RM. Insulin-like growth factor-1, insulin-like growth factor-binding protein-3, growth hormone, and mammographic density in the Nurses’ Health Studies. Breast Cancer Res Treat. 2012;136(3):805–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Fahy E, Cotter D, Sud M, Subramaniam S. Lipid classification, structures and tools. Biochim Biophys Acta. 2011;1811(11):637–47.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. van Meer G, Voelker DR, Feigenson GW. Membrane lipids: where they are and how they behave. Nat Rev Mol Cell Biol. 2008;9(2):112–24.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Santos CR, Schulze A. Lipid metabolism in cancer. FEBS J. 2012;279(15):2610–23.

    Article  CAS  PubMed  Google Scholar 

  16. Quehenberger O, Dennis EA. The human plasma lipidome. N Engl J Med. 2011;365(19):1812–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Wenk MR. The emerging field of lipidomics. Nat Rev Drug Discov. 2005;4(7):594–610.

    Article  CAS  PubMed  Google Scholar 

  18. Toriola AT, Appleton CM, Zong X, Luo J, Weilbaecher K, Tamimi RM, Colditz GA. Circulating receptor activator of nuclear factor-κB (RANK), RANK ligand (RANKL), and mammographic density in premenopausal women. Cancer Prev Res (Phila). 2018;11(12):789–96.

    Article  CAS  PubMed  Google Scholar 

  19. Complex Lipids Targeted Panel [https://www.metabolon.com/lp/complex-lipid-panel/]

  20. Hastie T, Tibshirani R, Narasimhan B, Chu G: impute: impute: Imputation for microarray data. In. R package version 1.70.0.; 2022.

  21. van Buuren S, Groothuis-Oudshoorn K. Mice: multivariate imputation by chained equations in R. J Stat Softw. 2011;45(3):1–67.

    Article  Google Scholar 

  22. Maskarinec G, Woolcott C, Steude JS, Franke AA, Cooney RV. The relation of leptin and adiponectin with breast density among premenopausal women. Eur J Cancer Prev. 2010;19(1):55–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Stuedal A, Ursin G, Veierod MB, Bremnes Y, Reseland JE, Drevon CA, Gram IT. Plasma levels of leptin and mammographic density among postmenopausal women: a cross-sectional study. Breast Cancer Res. 2006;8(5):R55.

    Article  PubMed  PubMed Central  Google Scholar 

  24. McCormack VA, Dowsett M, Folkerd E, Johnson N, Palles C, Coupland B, Holly JM, Vinnicombe SJ, Perry NM, dos Santos SI. Sex steroids, growth factors and mammographic density: a cross-sectional study of UK postmenopausal Caucasian and Afro-Caribbean women. Breast Cancer Res. 2009;11(3):R38.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Hanna M, Dumas I, Orain M, Jacob S, Tetu B, Sanschagrin F, Bureau A, Poirier B, Diorio C. Association between expression of inflammatory markers in normal breast tissue and mammographic density among premenopausal and postmenopausal women. Menopause. 2017;24(5):524–35.

    Article  PubMed  Google Scholar 

  26. Dallal C, Garte S, Ragin C, Chen J, Lloyd S, Modugno F, Weissfeld J, Taioli E. Plasma leptin levels, LEPR Q223R polymorphism and mammographic breast density: a cross-sectional study. Int J Biol Markers. 2013;28(2):161–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Diorio C, Pollak M, Byrne C, Masse B, Hebert-Croteau N, Yaffe M, Cote G, Berube S, Brisson J. Levels of C-peptide and mammographic breast density. Cancer Epidemiol Biomarkers Prev. 2005;14(11 Pt 1):2661–4.

    Article  CAS  PubMed  Google Scholar 

  28. Diorio C, Pollak M, Byrne C, Masse B, Hebert-Croteau N, Yaffe M, Cote G, Berube S, Morin C, Brisson J. Insulin-like growth factor-I, IGF-binding protein-3, and mammographic breast density. Cancer Epidemiol Biomarkers Prev. 2005;14(5):1065–73.

    Article  CAS  PubMed  Google Scholar 

  29. Maskarinec G, Williams AE, Kaaks R. A cross-sectional investigation of breast density and insulin-like growth factor I. Int J Cancer. 2003;107(6):991–6.

    Article  CAS  PubMed  Google Scholar 

  30. Lai JH, Vesprini D, Zhang W, Yaffe MJ, Pollak M, Narod SA. A polymorphic locus in the promoter region of the IGFBP3 gene is related to mammographic breast density. Cancer Epidemiol Biomarkers Prev. 2004;13(4):573–82.

    Article  CAS  PubMed  Google Scholar 

  31. dos Santos SI, Johnson N, De Stavola B, Torres-Mejia G, Fletcher O, Allen DS, Allen NE, Key TJ, Fentiman IS, Holly JM, et al. The insulin-like growth factor system and mammographic features in premenopausal and postmenopausal women. Cancer Epidemiol Biomarkers Prev. 2006;15(3):449–55.

    Article  Google Scholar 

  32. Eiriksson FF, Nohr MK, Costa M, Bodvarsdottir SK, Ogmundsdottir HM, Thorsteinsdottir M. Lipidomic study of cell lines reveals differences between breast cancer subtypes. PLoS ONE. 2020;15(4): e0231289.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Giudetti AM, De Domenico S, Ragusa A, Lunetti P, Gaballo A, Franck J, Simeone P, Nicolardi G, De Nuccio F, Santino A, et al. A specific lipid metabolic profile is associated with the epithelial mesenchymal transition program. Biochim Biophys Acta Mol Cell Biol Lipids. 2019;1864(3):344–57.

    Article  CAS  PubMed  Google Scholar 

  34. Ward AV, Anderson SM, Sartorius CA. Advances in analyzing the breast cancer lipidome and its relevance to disease progression and treatment. J Mammary Gland Biol Neoplasia. 2021;26(4):399–417.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Eghlimi R, Shi X, Hrovat J, Xi B, Gu H. Triple negative breast cancer detection using LC-MS/MS Lipidomic profiling. J Proteome Res. 2020;19(6):2367–78.

    Article  CAS  PubMed  Google Scholar 

  36. Giallourou N, Urbaniak C, Puebla-Barragan S, Vorkas PA, Swann JR, Reid G. Characterizing the breast cancer lipidome and its interaction with the tissue microbiota. Commun Biol. 2021;4(1):1229.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Kimura T, Jennings W, Epand RM. Roles of specific lipid species in the cell and their molecular mechanism. Prog Lipid Res. 2016;62:75–92.

    Article  CAS  PubMed  Google Scholar 

  38. Mohammadzadeh F, Mosayebi G, Montazeri V, Darabi M, Fayezi S, Shaaker M, Rahmati M, Baradaran B, Mehdizadeh A, Darabi M. Fatty acid composition of tissue cultured breast carcinoma and the effect of Stearoyl-CoA desaturase 1 inhibition. J Breast Cancer. 2014;17(2):136–42.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Lodi M, Kiehl A, Qu FL, Gabriele V, Tomasetto C, Mathelin C. Lipid intake and breast cancer risk: is there a link? A new focus and meta-analysis. Eur J Breast Health. 2022;18(2):108–26.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Thurmer M, Gollowitzer A, Pein H, Neukirch K, Gelmez E, Waltl L, Wielsch N, Winkler R, Loser K, Grander J, et al. PI(18:1/18:1) is a SCD1-derived lipokine that limits stress signaling. Nat Commun. 2022;13(1):2982.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Knuplez E, Marsche G: An Updated Review of Pro- and Anti-Inflammatory Properties of Plasma Lysophosphatidylcholines in the Vascular System. Int J Mol Sci 2020, 21(12).

  42. Lucht SA, Eliassen AH, Bertrand KA, Ahern TP, Borgquist S, Rosner B, Hankinson SE, Tamimi RM. Circulating lipids, mammographic density, and risk of breast cancer in the Nurses’ Health Study and Nurses’ Health Study II. Cancer Causes Control. 2019;30(9):943–53.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Boyd NF, Connelly P, Byng J, Yaffe M, Draper H, Little L, Jones D, Martin LJ, Lockwood G, Tritchler D. Plasma lipids, lipoproteins, and mammographic densities. Cancer Epidemiol Biomarkers Prev. 1995;4(7):727–33.

    CAS  PubMed  Google Scholar 

  44. Carayol M, Licaj I, Achaintre D, Sacerdote C, Vineis P, Key TJ, Onland Moret NC, Scalbert A, Rinaldi S, Ferrari P. Reliability of serum metabolites over a two-year period: a targeted metabolomic approach in fasting and non-fasting samples from EPIC. PLoS ONE. 2015;10(8): e0135437.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We would like to thank all of the women who participated in this study.

Funding

This work was supported by funding from NIH/NCI (R01CA246592 to A.T. Toriola), and K.R. Getz was supported by T32CA190194. The content is solely the responsibility of the authors and does not represent the official view of the NIH.

Author information

Authors and Affiliations

  1. Division of Public Health Sciences, Department of Surgery, School of Medicine, Washington University, 660 South Euclid Avenue, Box 8100, St. Louis, MO, 63110, USA

    Kayla R. Getz, Myung Sik Jeon, Chongliang Luo, Jingqin Luo & Adetunji T. Toriola

  2. Siteman Cancer Center Biostatistics Shared Resource, Division of Public Health Sciences, Department of Surgery, Washington University School of Medicine, St. Louis, MO, USA

    Myung Sik Jeon, Chongliang Luo & Jingqin Luo

  3. Siteman Cancer Center, Washington University School of Medicine, St. Louis, MO, USA

    Adetunji T. Toriola

Authors
  1. Kayla R. Getz
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Myung Sik Jeon
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Chongliang Luo
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Jingqin Luo
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Adetunji T. Toriola
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

ATT conceptualized the study design and acquired the data. JL and MSJ analyzed the data. ATT, KRG and JL interpreted the data. ATT and KG wrote the first draft of the manuscript and all authors provided critical revision of the manuscript.

Corresponding author

Correspondence to Adetunji T. Toriola.

Ethics declarations

Ethics approval consent to participate:

The study was performed in accordance with the Declaration of Helsinki and was approved by the Washington University Institutional Review Board. All participants provided written informed consent.

Consent to publish

Not applicable.

Competing Interests

The authors have no relevant financial or non-financial interests to disclose.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1. Supplementary Figure 1.

Spearman Correlation Heatmap Plot across Lipid Sub-pathways and Mammographic Breast Density Measures. Supplementary Table 2. Multivariable Covariate-adjusted Linear Regression between Lipid Sub-pathways and Mammographic Breast Density Measures. Supplementary Table 3. Covariate-adjusted Least Square Means of Non-dense Volume (NDV) by Quartiles of Lipid Species that were Significantly Associated with NDV at a Bonferroni P-value<0.05. Supplementary Table 4: Covariate-adjusted Least Square Mean of Dense Volume and Non-dense Volume by Quartile of Lipid Sub-pathways. Supplementary Figure 2: Covariate-adjusted Proportional Odds Model between the 14 Lipid Sub-pathways and Volumetric Percent Density Categories.

Additional file 2. Supplementary Table 1.

Spearman Correlations across Lipid Species and Mammographic Breast Density Measures (Excel file).

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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) 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

Getz, K.R., Jeon, M.S., Luo, C. et al. Lipidome of mammographic breast density in premenopausal women. Breast Cancer Res 25, 121 (2023). https://doi.org/10.1186/s13058-023-01725-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13058-023-01725-1

Keywords

Breast Cancer Research

ISSN: 1465-542X

Contact us