Use of skincare products and risk of cancer of the breast and endometrium: a prospective cohort study
Environmental Health volume 18, Article number: 105 (2019)
Concerns have been raised that extensive use of personal care products that contain endocrine disrupting compounds increase the risk of hormone sensitive cancers.
To assess the effect of skincare product use on the risk of pre- and postmenopausal breast cancer, estrogen receptor positive (ER+) and negative (ER-) breast cancer and cancer of the endometrium.
We used data from 106,978 participants in the population-based Norwegian Women and Cancer cohort. Participants were categorized into non-, light, moderate, frequent and heavy users of skincare products based on self-reported use of hand and facial cream and body lotion. Cancer incidence information from the Cancer Registry of Norway was linked to individual data through the unique identity number of Norwegian citizens. Multivariable Cox proportional hazard regression was used to assess the effect of skincare product use on the risk of cancer of the breast and endometrium. We used multiple imputation by chained equations to evaluate the effect of missing data on observed associations.
We found no associations between use of skincare products and incidence of premenopausal breast cancer (frequent/heavy versus non−/light use: hazard ratio [HR] =1.10, 95% confidence interval [CI]: 0.92–1.32), postmenopausal breast cancer (heavy versus light use: HR = 0.87, 95% CI: 0.65–1.18, frequent versus light use: HR = 0.97, 95% CI: 0.88, 1.07) or endometrial cancer (frequent/heavy versus non−/light use: HR = 0.97, 95% CI: 0.79–1.20). Use of skincare products did not increase the risk of ER+ or ER- breast cancer and there was no difference in effect across ER status (0.58 ≤ pheterogeneity ≤ 0.99). The magnitude and direction of the effect estimates based on complete case analyses and multiple imputation were similar.
Heavy use of skincare products, i.e. creaming the body up to two times per day during mid-life, did not increase the risk of cancer of the breast or endometrium.
Use of skincare products such as body lotion, facial cream and hand cream is common among Norwegian women . Several components of these personal care products (PCPs) are classified as known or suspected endocrine disruptors (EDs), i.e., compounds that are able to interfere with the endocrine function in humans and wildlife . EDs are easily absorbed by the skin into the central blood circulation after dermal application, and have been detected in various concentrations in human blood, urine and breastmilk . Concerns have therefore been raised whether extensive use of cosmetics and skincare products could increase the risk of hormone related cancer, for instance breast cancer .
Among frequently used EDs in PCPs are phthalates, ultraviolet (UV) filters, triclosan and parabens . Phthalates are also commonly used as softeners in consumer products such as food packaging material, children’s toys and building material  and have been suggested to interfere with the testosterone production or action . Recently, phthalates have shown anti-estrogenic effect in breast cancer cell lines  and induced cell proliferation in normal breast cells . However, recent epidemiological studies reported no significant association between phthalate exposure and breast cancer [9, 10].
UV-filters are a large group of compounds used as constituents in sunscreen as they are able to absorb UV-radiation. They are also included in other PCPs to increase shelf-life , and have been detected in human urine [11, 12], breastmilk  and breast tissue . Many UV-filters exert estrogenic activity in in vivo/in vitro experiments , however epidemiological studies of endocrine disrupting effects of UV-filters are rare. Triclosan is primarily used as an antimicrobial agent in soap, toothpaste, cosmetics and pharmaceuticals and has shown endocrine disrupting properties in experimental settings . Prenatal triclosan concentrations was recently reported associated with earlier menarche in American girls, whereas there was no effect of prenatal or peripubertal concentrations of the UV filter Benzophenone-3 .
Parabens are alkyl esters of p-hydroxybenzoic acid which is naturally occurring in several plants and berries, such as blueberries, strawberries, red and white currants [17, 18]. For decades, methyl-, ethyl, propyl- and butyl-paraben were the most frequently used preservatives in skincare products due to their antimicrobial properties, low toxicity, cost and weak sensitizing properties. From the beginning of this century, restrictions on the use of parabens in consumer products have been implemented within the European Union as a result of scientific reporting of weak estrogenic activity of parabens. Several in vitro studies have shown that parabens are able to bind to the estrogen receptor and stimulate proliferation within human breast cancer cell lines [19, 20], also in concentrations similar to what have been detected in human breast cancer cells  and in human breast tissue . Subcutaneous exposure to butyl and isobutyl-paraben has also been linked to uterus enlargement in rodents as a result of estrogenic activity  and peripubertal concentrations of methyl- and propoyl-paraben have been associated with measures of pubertal timing in girls . In our previous study, we showed that women who cream their body once daily or more had elevated plasma concentrations of methyl-, ethyl- and propyl- parabens . Other studies have found similar results , and there are reason to believe that women who use skincare products frequently experience higher body burdens of other EDs as well, for instance phthalates, UV filters and triclosan [12, 26].
Epidemiological evidence of the effects of cosmetics, skincare products or constituents of PCPs on hormone sensitive cancers, such as cancer of the breast and endometrium, is sparse. A recent cohort study reported 13% increased risk of breast cancer by frequent use of skincare products and 15% increased risk by frequent use of beauty products . Parada et al. found positive associations between paraben concentrations in urine (methyl- propyl and ∑parabens) and prevalent breast cancer, and negative, non-significant associations with breast cancer mortality . Out of four case-control studies of antiperspirant use and breast cancer risk, two reported no effect [29, 30] and one found an increased risk of breast cancer by use of underarm cosmetics . The forth study reported an association between earlier age at breast cancer diagnosis and longer duration of deodorant use and underarm shaving . Use of skin lighteners has recently been reported not being associated with breast cancer . Thus, given the widespread use of skincare products and the potential estrogenic effects of product constituents, there is a clear lack of epidemiological studies addressing the effect of skincare product use on hormone sensitive cancers. In this population-based study, we aimed to investigate the effect of skincare product use on the risk of pre- and postmenopausal breast cancer, estrogen receptor positive (ER+) and negative (ER-) breast cancer and cancer of the endometrium.
Study design, participants and sub-samples
The Norwegian Women and Cancer study (NOWAC) is a national representative cohort study initiated in 1991 with the aim of exploring associations between lifestyle and cancer among Norwegian women . Women aged 30–70 years were randomly selected from the National Registry and invited to participate in the study through a mailed invitation letter to their home address that also included a detailed questionnaire. Women that agreed to participate have been followed-up regularly with consecutive questionnaires. Since the initiation of the study, the cohort has been expanded several times and includes today about 172,000 women. All participants have answered between one and four questionnaires regarding their current health status and lifestyle factors. The questionnaires have been distributed in waves and the level of details has varied between the questionnaires. The external validity of NOWAC has been studied extensively and been found satisfactory and without selection bias .
Questionnaires that included questions about usage frequencies of hand cream, facial cream and body lotion were mailed to 192,648 women, of which 114,794 responded. The questions regarding usage of skincare products were the same in all questionnaires. We excluded 1431 women that belonged to a certain wave of recruitment that was not randomly sampled and 5 women that died or emigrated before their baseline questionnaire were registered. Additionally 6286 women were excluded as they were diagnosed with cancer (except non-melanoma skin cancer, International Statistical Classification of Diseases, Injuries and Causes of Death Revision 10 code [ICD-10]: C44) prior to answering the questionnaire. Women reporting implausible values on height, weight, age at menarche, age at first full term pregnancy and age at menopause were also excluded from the study (n = 94). Thus, the final study sample used for studying the associations between skincare product use and breast cancer subtypes included 106,978 women (Fig. 1).
After excluding 85,721 women that were either postmenopausal (n = 80,245), used menopausal hormone therapy (MHT, n = 3543) or were ≥ 53 years at baseline (n = 1933), the study sample for premenopausal breast cancer included 21,257 women (Fig. 1). Women were considered postmenopausal if they stated that their menstrual bleeding had stopped or reported use of MHT if they were ≥ 53 years. This cut-off is based on the definition used in the Million Women Study  and has been validated against plasma concentrations of sex hormones in NOWAC .
The study sample for postmenopausal breast cancer included 106,328 women (Fig. 1). Women that were not categorized as postmenopausal at baseline were included in the study from the age of 53 years or the age of reported menopause. Due to this left truncation, we excluded 567 women that were diagnosed with cancer and 83 that died or emigrated before the start of follow-up. The study sample for endometrial cancer included 95,577 women. Here, we excluded women that reported hysterectomy (n = 11,401).
Measures of skincare products use and included covariates
In the NOWAC questionnaires, women were asked to record how often they used skincare products such as body lotion, hand cream and facial cream (never/seldom, 1–3 times/month, 1 time/week, 2–4 times/week, 5–6 times/week, 1 time/day and ≥ 2 times/day). The recorded frequencies were converted to percentage body surface covered with cream per day. Use of hand cream, facial cream and body lotion once per day corresponded to 100% of body surface covered in cream per day. Calculated percentages were later categorized into five groups; non-users (0–0.001%), light users (0.002- < 35.0%), moderate users (35.0- < 65.0%), frequent users (65.0- < 115.0%) and heavy users (115.0–200%), where 200% corresponded to using hand cream, facial cream and body lotion twice per day. The conversion is described in detail in Aniansson et al.  and has also been used previously to assess the correlation between skincare product use and plasma concentrations of parabens .
Information on the women’s age was extracted from the National Registry and the other covariates from the baseline NOWAC questionnaire. Education was categorized based on years of completed schooling corresponding to secondary school (< 10 years), high school (10–12 years) and higher education (> 12 years). Body mass index (BMI) was calculated from self-reported body weight (kg) and height (m), and categorized as under−/normal weight (< 25.0 kg/m2), overweight (25.0–29.9 kg/m2) or obese (≥30 kg/m2). The women’s menopausal status (premenopausal, perimenopausal, postmenopausal, unknown) was determined by reported regular menstrual bleeding or not, use of MHT and the women’s age . Age at menarche was categorized into three groups (≤12 years, 13–14 years, ≥15 years) and age at first full term pregnancy and parity were combined into one variable (nullipara, < 30 years at first full term pregnancy and unipara, ≥30 years and unipara, < 30 years and multipara, ≥30 years and multipara) in the statistical analysis. Use of oral contraceptives (OCs) and MHT were categorized in two ways; never/former/current was used for the breast cancer analysis, and the never/ever categorization for endometrial cancer. Use of intrauterine device (IUD) was also assessed as never/ever use. Smoking status was categorized as never, former and current smoker. Physical activity was recorded on an ordinal 10-point scale and categorized as low (1–4), moderate (5–6) and high (≥7). Alcohol consumption was recorded by a food frequency questionnaire and used in continuous scale.
Cancer incidence, tumor receptor status, death and emigration
The NOWAC study receives annual updates from the Cancer Registry of Norway in order to identify study participants diagnosed with cancer during the preceding year. Information on date of diagnosis and hormone receptor status for the breast cancer cases are also included. In the present study, women diagnosed with a first primary invasive malignant neoplasm of the breast (ICD-10: C50) and uterus (ICD-10: C54) were identified. Endometrial cancer cases were identified from morphological codes 8380, 8382, 8480, 8481, 8560, 8570, 8020, 8041, 8045, 8255, 8310, 8441, 8460 or 8323 from the International Classification of Disease for Oncology ICD-O-3. Information about deaths and emigration was extracted from the Causes of Death Registry and the National Registry. End of follow-up was December 31, 2016.
We followed the official Norwegian thresholds for classification of hormone receptor status . From February 2012 and onwards, a tumor was classified as ER- if it displayed < 1% reactivity, whereas prior to February 2012, the threshold for ER- was < 10%. Accordingly, tumors classified as ER+ displayed ≥10% reactivity prior to February 2012 and ≥ 1% after. The change in threshold was due to changes in clinical practice at Norwegian hospitals.
We investigated the associations between skincare product use and the risk of breast and endometrial cancer using Cox proportional hazard regression with age as the time scale. Entry time was age at answering the questions on skincare products, and exit time was age at cancer diagnosis, death, emigration or end of follow-up, whichever occurred first. In the analysis of premenopausal breast cancer, women were censored at the age of 53 or age at reported menopause if that occurred prior to cancer diagnosis, death, emigration or end of follow-up. Due to few cancer cases among the non- and heavy users of skincare products, we used three exposure groups (non−/light, moderate, frequent/heavy) in the analysis of premenopausal breast cancer, endometrial cancer and the subtype-specific analysis of breast cancer. Light users (or non−light users) were used as the reference group in all analysis. We used the “purposeful selection method” described by Hosmer and Lemeshow  to evaluate which covariates to include in the final multivariable models. Covariates included in each regression model are listed in the table footnotes. To test for linear trend, we replaced the group identifier with the median use of % skin covered in cream per day per group, and included that variable in the multivariable models. Additionally, to model the relationship between cancer risk and use of skincare products in continuous scale and to allow for non-linear effects, we fitted regression models with natural cubic spline transformations (4 knots) of the exposure variable (% skin creamed per day). The knots were placed at equally spaced percentiles. We evaluated non-linearity by testing the null hypothesis of the second and third spline coefficients jointly equalled zero. We assessed effect modification by MHT use for postmenopausal breast cancer, but did not evaluate any other interactions, due lack of statistical power. Departures from the proportional hazard assumption were assessed by inspection of the Schoenfeldts residuals. Participants with missing values of included covariates were excluded from the complete case analysis.
To assess differences in sub-type specific breast cancer risk, we used Cox proportional hazard regression for ER+ and ER- separately, where women who were diagnosed with another breast cancer subtype, were censored at the time of diagnosis . We tested for heterogeneity in associations between subtypes by a chi-squared (contrast) test .
In the complete case analysis of premenopausal breast cancer, 92% of the observations were included. Corresponding proportions for postmenopausal breast cancer, ER+, ER- breast cancer and endometrial cancer were 83, 84, 90 and 69%, respectively. To evaluate the effect of missing information on the observed results, we used multiple imputation by chained equations to obtain 20 imputed datasets with complete observations for each outcome. The hazard ratios (HRs) estimated on the imputed data sets were pooled together using Rubin’s rule to obtain valid statistical inferences .
Several additional analyses were performed. We assessed the associations between recorded usage frequencies of body lotion (categorized as never/seldom, 1–4 times/month, 2–6 times/week and 1- ≥ 2 times/day) and breast and endometrial cancer incidence. Furthermore, we summarized the usage frequencies of skincare products (scores from 0 [no use of any of skincare products] to 18 [use of hand cream, facial cream and body lotion ≥2 times/day]; categorized as low [0–5], moderate low [6–10], moderate high [11–14] and high [15–18]) and studied the associations with risk of breast and endometrial cancer.
All P-values were two-sided and a 5% level of significance was used. The statistical analysis was conducted using Stata, version 15.1 (StataCorp, College Station, TX, USA).
In the total study sample of 106,978 women, mean (standard deviation [SD]) age at enrollment was 54.7 (4.8) years. During the mean follow-up time of 10.7 (2.6) years, 3408 women were diagnosed with breast cancer and 681 with endometrial cancer. Mean age at breast and endometrial cancer diagnosis were 60.1 (5.5) and 61.0 (5.2) years, respectively.
Age at enrollment was similar in all categories of skincare product use. Compared to light users, a larger proportion of women with under−/normal weight and postmenopausal status were found among the heavy users. Further, heavy users were slightly younger when they had their first child and a larger proportions had < 3 children compared to the light users. Heavy users were also more frequent current and former users of OC, MHT and cigarettes, more physically active and reported higher intake of alcohol than the light users (Table 1).
Frequent/heavy users of skincare products did not experience increased risk of premenopausal breast cancer (HR = 1.10, 95% confidence interval [CI]: 0.92–1.32), postmenopausal breast cancer (heavy versus light use: HR = 0.87, 95% CI: 0.65–1.18, frequent versus light use: HR = 0.97, 95% CI: 0.88, 1.07) or endometrial cancer (HR = 0.97, 95% CI: 0.79–1.20) compared to the non−/light users. Nor did the moderate users, and there was no linear trend across effect estimates (0.27 ≤ ptrend ≤ 0.63) In fact, moderate use of skincare products was associated with 23% decreased risk of endometrial cancer (HR = 0.77, 95% CI: 0.61, 0.98) (Table 2). The regression models with natural cubic spline transformations of %skin covered in cream per day, displayed no increased risk of pre- or postmenopausal breast cancer or endometrial cancer by increasing use of skincare products. However, an inverse association between skincare product use and postmenopausal breast cancer risk was suggested (Fig. 2 and Additional file 1). Further, there was no evidence of an interaction between skincare product use and MHT use in relation to postmenopausal breast cancer risk (p interaction(Wald test) = 0.44).
The subtype-specific breast cancer analyses suggested no increased risk of ER+ or ER- breast cancer by moderate or frequent/heavy use of skincare products and no difference in effect by estrogen receptor status (0.58 ≤ pheterogeneity ≤ 0.99, Table 3).
The magnitude and direction of the effect estimates based on complete case analyses and multiple imputation were similar, except for the lack of inverse association between moderate skincare product use and risk of endometrial cancer (Additional files 2 and 3). Thus, missing values did not bias the observed associations in the complete case analyses.
In the additional analyses, we found no significant associations between usage frequencies of body lotion and the risk of pre- or postmenopausal breast cancer, endometrial cancer, ER+ or ER- breast cancer (body lotion: 0.30 ≤ ptrend ≤ 0.55) (Additional file 4). However, there was a significant inverse trend for usage frequencies of skincare products combined and postmenopausal breast cancer and ER+ breast cancer (0.02 ≤ ptrend ≤ 0.046), although we did not find significantly reduced HRs for any of the categories (postmenopausal breast cancer: moderate low versus low use: HR = 1.01, 95% CI: 0.88,1.17; moderate high versus low use: HR = 0.94, 95% CI: 0.81,1.08; high versus low use: HR = 0.89, 95% CI: 0.77,1.03; ER+ breast cancer: moderate low versus low use: HR = 1.07, 95% CI: 0.91,1.24; moderate high versus low use: HR = 0.95, 95% CI: 0.82,1.11; high versus low use: HR = 0.93, 95% CI: 0.80,1.09) (Additional file 5). Endometrial cancer and ER- breast cancer were not associated with usage frequencies of skincare products combined (0.20 ≤ ptrend ≤ 0.45) (Additional file 5).
In this large and national representative cohort of women in Norway, we found no evidence of increased risk of pre- or postmenopausal breast cancer by use of skincare products. Nor did we observe evidence of increased risk of ER+ or ER- breast cancer or endometrial cancer. We have previously reported that frequent and heavy users of skincare products in the NOWAC study experience elevated plasma concentrations of methyl-, ethyl and propyl-parabens, compounds that exert endocrine disrupting properties and have been linked to breast cancer in in vivo and in vitro studies, as well in a recent case control study [28, 44, 45]. We can now conclude that these women did not experience increased risk of hormone sensitive cancers. In fact, an inverse relationship between skincare product use modelled in continuous scale with restricted cubic splines and postmenopausal breast cancer risk was observed, which could be a chance finding or a result of residual confounding.
To this date, our study is the largest prospective cohort study assessing the effect of skincare product use on cancer risk. Our results may seem to contrast Taylor et al. who recently reported increased risk of breast cancer by frequent use of skincare products among non-Hispanic white women from the United States . However, that study included both in situ and invasive breast cancer and stratified analyses displayed no significant association between frequent use of skincare products and invasive breast cancer, which is in correspondence with our results. Further, the authors found no significant associations between skincare product use and ER+ or ER- breast cancer which is in line with our findings. Notably, Taylor and colleagues included nine different skincare products (face cream, cleansing cream, anti-aging cream, foot cream, body lotion, hand lotion, petroleum jelly, and two different talcum powders) in comparison to three in the present study, which could explain differences in results. In the analyses stratified by breast cancer type, Taylor et al. observed a borderline significant increased risk of in situ breast cancer among frequent users of skincare products compared to infrequent users, which may suggest that the association with skincare product use may be different in in situ and invasive breast cancer. Additionally, skincare product constituents in the United States and in Europe may vary due to differences in regulations, which could also affect the results.
Recently, Parada et al. reported positive associations between urinary concentrations of parabens and breast cancer status that may seem to contrast our results . It is, however, challenging to compare two such different studies of which one rely on self-reported skincare product use collected years prior to breast cancer diagnosis and the other on paraben measurements in urine samples collected after breast cancer diagnosis. Clearly, Parada et al. showed that women with breast cancer experienced slightly higher concentrations of parabens compared to healthy controls, however reverse causation cannot be fully ruled out as the samples were donated after disease diagnosis.
Our results further suggest that use of skincare products does not increase the risk of endometrial cancer, an estrogen sensitive cancer type that have increased in incidence over the last 50 years in Norway . In fact, compared to non−/light users of skincare products, moderate users experienced 23% decreased risk of endometrial cancer. However, due to the lack of linear trend across effect estimates and the lack of inverse association when modelling skincare product use in continuous scale and after multiple imputation, we consider this a random finding.
Use of PCPs is considered an important source of paraben exposure for humans [5, 12]. Several in vivo and in vitro studies have also demonstrated that parabens have the ability to interfere with the endocrine system and may increase the risk of breast cancer through various mechanisms [21, 22, 47, 48]. As a consequence, the use of parabens in consumer products within the European Union has been restricted and other biocides such as isothiazolinones have been introduced as replacements compounds. Lately, concern have been raised by dermatologists as these compounds poses much stronger sensitizing properties than parabens and an increasing incidence of allergic contact dermatitis is therefore expected . As there is a lack of epidemiological evidence of the effects of parabens on human health, one may question whether replacing parabens with a more potent sensitizer is wise? There are however, several possible explanations for the disagreement between in vivo/in vitro studies and epidemiological studies that include i) the exposure classification in large scale epidemiological studies is challenging, as exposure usually is estimated from questionnaires. This limits the individual exposure assessment as detailed plasma concentration of EDs for each participant is not available. The questionnaire based approach has however been used for decades in nutritional epidemiology and if a proper validation study is performed, as in the our pervious study , the questionnaires will be able to rank individuals into broader categories based on their exposure to EDs from PCPs. ii) It is possible that humans have a different susceptibility to EDs compared to cell lines and experimental animals and that effects seen in in vivo and/or in vitro studies are not as pronounced among humans. Animal models will continue to provide important insight into the human biology, however there are also important differences between mice and men that needs to be addressed . iii) The exposure measure in epidemiological studies may not necessarily cover exposure during etiological relevant periods. A large challenge when studying EDs in relations to hormone sensitive cancers is the limited knowledge about which period in life humans are most susceptible for ED exposure in relation to the outcome of interest. Recently, Harley et al. reported that urinary concentrations of parabens at age 9 were associated with earlier breast development, menarche and development of pubic hair . They also found that prenatal triclosan concentrations were associated with earlier menarche. Thus, this study suggest that prenatal and peripubertal exposure to EDs from PCPs may modify important risk factors for breast and endometrial cancer, such as age at menarche. Many questionnaire-based epidemiological studies assess exposure during the previous year and assume that it reflect past exposure. However, as lifestyle and consumers products often change from puberty to adulthood, the exposure assessed later in life does not necessarily reflect exposure during etiologically sensitive time windows. Thus, due to the widespread use of cosmetics and skincare products among young women today, the precautionary principle of phasing out potentially harmful EDs may be favorable even though they may be replaced by compounds with stronger sensitizing properties. However, there is clearly a lack of epidemiological studies that assess the effect of ED exposure from PCPs during various life-phases and the subsequent risk of hormone related cancer.
The strengths of this study includes the large and national representative sample, the linkage to national registries that ensures complete information about death, emigration and cancer diagnosis and the detailed information about sociodemographic, reproductive and lifestyle variables that we included in our analyses. Our previous study of plasma concentrations of parabens in relation to reported use of skincare products confirmed that individuals classified as frequent or heavy users of skincare products experience elevated circulatory concentrations of methyl-, ethyl- and propyl parabens compared to the non-users . Philippat et al. also confirmed a positive association between the number of PCPs used and increasing urinary paraben concentrations. They further suggested that questionnaires can be used to measure exposure to parabens . However, skincare products are not only sources of parabens; they are sources of other EDs with different properties and mechanisms of action, individually and as mixtures. As we have not measured other EDs than parabens in plasma from skincare products users in NOWAC, we cannot be certain that those classified as frequent or heavy users also experienced the highest concentrations of other EDs such as phthalates, UV filters and triclosan that individually or combined may exert estrogenic activity.
We have previously assessed the change in use of skincare products over an eight year period. We found moderate agreement (weighted kappa: 0.52) and that almost 92% of women were classified into the same user group ±1 category . Thus, the use of skincare products among these middle-aged women is a well-established habit that is relatively stable over time. Nevertheless, we cannot assume that the reported use of skincare products reflect exposure during childhood/puberty/adolescence. In fact, given the advanced age of the study participants, it is likely that their use of skincare products during childhood and adolescence were lower than what is common among children and adolescents today as the general Norwegian household economy was much weaker in the 1950’s and 60’s than nowadays and the availability to skincare products were much more limited. Caution should therefore be made if generalizing these results to younger birth cohorts (born 1960 and onwards).
The NOWAC questionnaires did not include questions about which part of the body that was covered in body lotion; therefore we assumed that one application of body lotion covered nearly the whole body (91%). This may be seen upon as a limitation as many women only use body lotion on arms and legs, suggesting that our measure “% skin covered in cream per day” may overestimate the exposure. However, we have previously shown that “% skin covered in cream per day” was strongly associated with plasma concentrations of parabens within the NOWAC cohort , and we believe that this measure better reflect blood concentrations of EDs from skincare product use than the frequency of use. Nevertheless, we conducted several additional analyses of reported usage frequencies of skincare products (Additional files 4 and 5). These analyses confirmed our main results, i.e. no increased risk of breast or endometrial cancer by frequent use of skincare products. Still, it is possible that some NOWAC participants over- or under-reported their usage frequencies of skincare products. As all information was collected before cancer diagnosis, this measurement error is likely non-differential. In our age-adjusted models, assuming non-differential, nonsystematic errors, misclassification would attenuate the HR of the higher skincare category, but the test for trend would be valid . As we have several exposure categories and many confounding variables included in our multivariable models, non-differential misclassification can bias the estimates both towards and away from the null .
Our study was also unable to address whether the location of body lotion application was associated with breast or endometrial cancer risk. This may be important if dermal absorption of EDs on specific areas is more relevant for breast/endometrial carcinogenesis than continuously elevated systemic concentrations. Finally, we did not have information on antiperspirant and cosmetic use, which could be important sources of paraben exposure. Nevertheless, previous studies have called for epidemiological investigations in order to confirm or refute in vivo/in vitro findings regarding paraben exposure and the risk of breast cancer. Our finding that use of skincare products during mid-life was not associated with breast and endometrial cancer risk adds important knowledge to this research field.
Our population-based prospective cohort study provides evidence that heavy use of skincare products, i.e. creaming the body up to two time per day during mid-life, does not increase the risk of pre- or postmenopausal breast cancer, ER+ or ER- breast cancer or endometrial cancer.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Estrogen receptor positive
Estrogen receptor negative
Personal care product
the Norwegian Women And Cancer study
International Statistical Classification of Diseases, Injuries and Causes of Death Revision 10
Menopausal hormone therapy
Body mass index
Aniansson B, Veierød MB, Rylander C, Lund E, Sandanger TM. Characterization of heavy users of skin care products among Norwegian women from 2003 to 2011. Arch Public Health. 2016;74(1):53.
Witorsch RJ, Thomas JA. Personal care products and endocrine disruption: a critical review of the literature. Crit Rev Toxicol. 2010;40(Suppl 3):1–30.
Azzouz A, Rascon AJ, Ballesteros E. Determination of free and conjugated forms of endocrine-disrupting chemicals in human biological fluids by GC-MS. Bioanalysis. 2016;8(11):1145–58.
Darbre PD. Environmental oestrogens, cosmetics and breast cancer. Best Pract Res Clin Endocrinol Metab. 2006;20(1):121–43.
Guo Y, Kannan K. A survey of phthalates and parabens in personal care products from the United States and its implications for human exposure. Environ Sci Technol. 2013;47(24):14442–9.
Schettler T. Human exposure to phthalates via consumer products. Int J Androl. 2006;29(1):134–9.
Okazaki H, Takeda S, Matsuo S, Matsumoto M, Furuta E, Kohro-Ikeda E, et al. Inhibitory modulation of human estrogen receptor alpha and beta activities by dicyclohexyl phthalate in human breast cancer cell lines. J Toxicol Sci. 2017;42(4):417–25.
Chen FP, Chien MH, Chen HY, Ng YT. Effects of phthalates on normal human breast cells co-cultured with different fibroblasts. PLoS ONE. 2018;13(6):e0199596.
Fu Z, Zhao F, Chen K, Xu J, Li P, Xia D, et al. Association between urinary phthalate metabolites and risk of breast cancer and uterine leiomyoma. Reprod Toxicol. 2017;74:134–42.
Parada H Jr, Gammon MD, Chen J, Calafat AM, Neugut AI, Santella RM, et al. Urinary phthalate metabolite concentrations and breast Cancer incidence and survival following breast Cancer: the Long Island breast Cancer study project. Environl Health Perspect. 2018;126(4):047013.
Calafat AM, Wong LY, Ye X, Reidy JA, Needham LL. Concentrations of the sunscreen agent benzophenone-3 in residents of the United States: National Health and Nutrition Examination Survey 2003--2004. Environ Health Perspect. 2008;116(7):893–7.
Philippat C, Bennett D, Calafat AM, Picciotto IH. Exposure to select phthalates and phenols through use of personal care products among Californian adults and their children. Environ Res. 2015;140:369–76.
Rodriguez-Gomez R, Zafra-Gomez A, Dorival-Garcia N, Ballesteros O, Navalon A. Determination of benzophenone-UV filters in human milk samples using ultrasound-assisted extraction and clean-up with dispersive sorbents followed by UHPLC-MS/MS analysis. Talanta. 2015;134:657–64.
Barr L, Alamer M, Darbre PD. Measurement of concentrations of four chemical ultraviolet filters in human breast tissue at serial locations across the breast. J Appl Toxicol. 2018;38(8):1112–20.
Alamer M, Darbre PD. Effects of exposure to six chemical ultraviolet filters commonly used in personal care products on motility of MCF-7 and MDA-MB-231 human breast cancer cells in vitro. J Appl Toxicol. 2018;38(2):148–59.
Harley KG, Berger KP, Kogut K, Parra K, Lustig RH, Greenspan LC, et al. Association of phthalates, parabens and phenols found in personal care products with pubertal timing in girls and boys. Hum Reprod. 2019;34(1):109–17.
Maatta K, Kamal-Eldin A, Torronen R. Phenolic compounds in berries of black, red, green, and white currants (Ribes sp.). Antioxid Redox Signal. 2001;3(6):981–93.
Sellappan S, Akoh CC, Krewer G. Phenolic compounds and antioxidant capacity of Georgia-grown blueberries and blackberries. J Agric Food Chem. 2002;50(8):2432–8.
Okubo T, Yokoyama Y, Kano K, Kano I. ER-dependent estrogenic activity of parabens assessed by proliferation of human breast cancer MCF-7 cells and expression of ERalpha and PR. Food Chem Toxicol. 2001;39(12):1225–32.
Byford JR, Shaw LE, Drew MG, Pope GS, Sauer MJ, Darbre PD. Oestrogenic activity of parabens in MCF7 human breast cancer cells. J Steroid Biochem Mol Biol. 2002;80(1):49–60.
Charles AK, Darbre PD. Combinations of parabens at concentrations measured in human breast tissue can increase proliferation of MCF-7 human breast cancer cells. J Appl Toxicol. 2013;33(5):390–8.
Pan S, Yuan C, Tagmount A, Rudel RA, Ackerman JM, Yaswen P, et al. Parabens and human epidermal growth factor receptor ligand cross-talk in breast Cancer cells. Environ Health Perspect. 2016;124(5):563–9.
Routledge EJ, Parker J, Odum J, Ashby J, Sumpter JP. Some alkyl hydroxy benzoate preservatives (parabens) are estrogenic. Toxicol Appl Pharmacol. 1998;153(1):12–9.
Sandanger TM, Huber S, Moe MK, Braathen T, Leknes H, Lund E. Plasma concentrations of parabens in postmenopausal women and self-reported use of personal care products: the NOWAC postgenome study. J Expo Sci Environ Epidemiol. 2011;21(6):595–600.
Larsson K, Ljung Bjorklund K, Palm B, Wennberg M, Kaj L, Lindh CH, et al. Exposure determinants of phthalates, parabens, bisphenol a and triclosan in Swedish mothers and their children. Environ Int. 2014;73:323–33.
Sakhi AK, Sabaredzovic A, Cequier E, Thomsen C. Phthalate metabolites in Norwegian mothers and children: Levels, diurnal variation and use of personal care products. Sci Total Environ. 2017;599–600:1984–92.
Taylor KW, Troester MA, Herring AH, Engel LS, Nichols HB, Sandler DP, et al. Associations between personal care product use patterns and breast Cancer risk among white and black women in the sister study. Environ Health Perspect. 2018;126(2):027011.
Parada H Jr, Gammon MD, Ettore HL, Chen J, Calafat AM, Neugut AI, et al. Urinary concentrations of environmental phenols and their associations with breast cancer incidence and mortality following breast cancer. Environ Int. 2019;130:104890.
Mirick DK, Davis S, Thomas DB. Antiperspirant use and the risk of breast cancer. J Natl Cancer Inst. 2002;94(20):1578–80.
Fakri S, Al-Azzawi A, Al-Tawil N. Antiperspirant use as a risk factor for breast cancer in Iraq. East Mediterr Health J. 2006;12(3–4):478–82.
Linhart C, Talasz H, Morandi EM, Exley C, Lindner HH, Taucher S, et al. Use of underarm cosmetic products in relation to risk of breast Cancer: a case-control study. EBioMedicine. 2017;21:79–85.
McGrath KG. An earlier age of breast cancer diagnosis related to more frequent use of antiperspirants/deodorants and underarm shaving. Eur J Cancer Prev. 2003;12(6):479–85.
Brinton LA, Figueroa JD, Ansong D, Nyarko KM, Wiafe S, Yarney J, et al. Skin lighteners and hair relaxers as risk factors for breast cancer: results from the Ghana breast health study. Carcinogenesis. 2018;39(4):571–9.
Lund E, Dumeaux V, Braaten T, Hjartaker A, Engeset D, Skeie G, et al. Cohort profile: the Norwegian women and cancer study - NOWAC - Kvinner og kreft. Int J Epidemiol. 2008;37(1):36–41.
Lund E, Kumle M, Braaten T, Hjartåker A, Bakken K, Eggen E, et al. External validity in a population-based national prospective study - the Norwegian women and Cancer study (NOWAC). Cancer Causes Control. 2003;14(10):1001–8.
Beral V. Breast cancer and hormone-replacement therapy in the million women study. Lancet. 2003;362(9382):419–27.
Waaseth M, Bakken K, Dumeaux V, Olsen KS, Rylander C, Figenschau Y, et al. Hormone replacement therapy use and plasma levels of sex hormones in the Norwegian Women and Cancer postgenome cohort - a cross-sectional analysis. BMC Womens Health. 2008;8:1.
The Directorate of Health ( Helsedirektoratet). Nasjonalt handlingsprogram med retningslinjer for diagnostikk, behandling og oppfølging av pasienter med brystkreft. 2011. https://helsedirektoratet.no/Lists/Publikasjoner/Attachments/920/Brystkreft-handlingsprogram-med-retningslinjer-IS-2316.pdf. Accessed 11 May 2019.
Hosmer D, Lemeshow S. Logistic regression. In: Veierod MB, editor. Medical statistics in clinical and epidemiological research. Oslo, Norway: Gyldendahl Akademisk; 2012.
Prentice RL, Kalbfleisch JD, Peterson AV Jr, Flournoy N, Farewell VT, Breslow NE. The analysis of failure times in the presence of competing risks. Biometrics. 1978;34(4):541–54.
Wang M, Spiegelman D, Kuchiba A, Lochhead P, Kim S, Chan AT, et al. Statistical methods for studying disease subtype heterogeneity. Stat Med. 2016;35(5):782–800.
Rubin DB, Schenker N. Multiple imputation for interval estimation from simple random samples with ignorable nonresponse. J Am Stat Assoc. 1986;81(394):366–74.
Darbre PD, Charles AK. Environmental oestrogens and breast cancer: evidence for combined involvement of dietary, household and cosmetic xenoestrogens. Anticancer Res. 2010;30(3):815–27.
Darbre PD, Harvey PW. Paraben esters: review of recent studies of endocrine toxicity, absorption, esterase and human exposure, and discussion of potential human health risks. J Appl Toxicol. 2008;28(5):561–78.
Cancer Registry of Norway. Cancer in Norway 2017. 2018. https://www.kreftregisteret.no/globalassets/cancer-in-norway/2017/cin-2017.pdf. Accessed 11 May 2019.
Gopalakrishnan K, Teitelbaum SL, Lambertini L, Wetmur J, Manservisi F, Falcioni L, et al. Changes in mammary histology and transcriptome profiles by low-dose exposure to environmental phenols at critical windows of development. Environ Res. 2017;152:233–43.
Sun L, Yu T, Guo J, Zhang Z, Hu Y, Xiao X, et al. The estrogenicity of methylparaben and ethylparaben at doses close to the acceptable daily intake in immature Sprague-Dawley rats. Sci Rep. 2016;6:25173.
Sasseville D, Alfalah M, Lacroix JP. Parabenoia" debunked, or "Who's afraid of Parabens? Dermatitis. 2015;26(6):254–9.
Mestas J, Hughes CC. Of mice and not men: differences between mouse and human immunology. J Immunol. 2004;172(5):2731–8.
Buonaccorsi JP. Measurement error: models, methods, and applications: Boca Raton, FL, USA; Chapman and Hall/CRC; 2010.
Jurek AM, Greenland S, Maldonado G, Church TR. Proper interpretation of non-differential misclassification effects: expectations vs observations. Int J Epidemiol. 2005;34(3):680–7.
We would like to thank Dr. Génevièr in France who has shown great enthusiasm and interest in this topic. We are also grateful to all women who have participated in the NOWAC study.
The publication charges for this article have been funded by a grant from the publication fund of UiT The Arctic University of Norway.
Ethics approval and consent to participate
NOWAC has been approved by the Norwegian Data Inspectorate and the Regional Committee for Medical Research Ethics in Northern Norway. All participants have provided written informed consent.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Hazard ratios (HRs) and 95% confidence intervals (CIs) for the association between skincare product use and risk of cancer. Skincare product use modelled in continuous scale using restricted cubic spline transformations of “% skin covered in cream per day” with 4 knots.
Hazard ratios (HRs) and 95% confidence intervals (CIs) for the associations between skincare product use and risk of pre- and postmenopausal breast cancer and endometrial cancer after multiple imputation by chained equations of missing values of included covariates.
Hazard ratios (HRs) and 95% confidence intervals (CIs) for the associations between skincare product use and risk of estrogen receptor positive (ER+) and negative (ER-) breast cancer after multiple imputation by chained equations of missing values of included covariates.
Hazard ratios (HRs) and 95% confidence intervals (CIs) for the associations between usage frequencies of body lotion and risk of pre- and postmenopausal breast cancer, endometrial cancer, ER+ and ER- breast cancer.
Hazard ratios (HRs) and 95% confidence intervals (CIs) for the associations between usage frequencies of body lotion, hand cream and facial cream combined and risk of pre- and postmenopausal breast cancer, endometrial cancer, ER+ and ER- breast cancer.
About this article
Cite this article
Rylander, C., Veierød, M.B., Weiderpass, E. et al. Use of skincare products and risk of cancer of the breast and endometrium: a prospective cohort study. Environ Health 18, 105 (2019). https://doi.org/10.1186/s12940-019-0547-6
- Skin care
- Body lotion
- Hand cream
- Facial cream
- Personal care products