Effects of statin use on volumetric mammographic density: results from the KARMA study

Background Epidemiological data on statins and breast cancer risk have been inconclusive. The aim of this study was to clarify the role of statins in breast cancer risk by studying their effect on mammographic density. Methods The KARolinska MAmmography project for risk prediction of breast cancer (KARMA) includes 70,877 women who underwent either a screening or clinical mammography from January 2011 to December 2013. In total, 41,102 women responded to a web-based questionnaire, and had raw digital mammograms stored. Volumetric mammographic density was measured using Volpara™ and information on statin use was obtained through linkage with the Swedish National Prescription Register. Analysis of covariance was used to study the effect of statin use on mammographic density, adjusting for a large set of potential confounders. We also studied the effects of statin class and treatment duration and tested for potential effect modification by hormone replacement therapy (HRT). Results Statin use was recorded in 3,337 women (8.1 %) of the study population and lipophilic statins was the most commonly prescribed type (93.4 % of all statin users). After multivariable adjustment, percent dense volume was lower in statin users than in non-users (P < 0.001). This association was explained by a larger absolute non-dense volume in statin users (P < 0.001). Overall, no difference in absolute dense volume was detected, but interaction analyses revealed a larger dense volume among statin users who reported ever HRT use (P = 0.03). No differential effects were observed according to statin lipophilicity and treatment duration. Conclusions We observed no overall effect of statin use on mammographic density in terms of absolute dense volume, although a larger absolute dense volume was observed in statin users who reported ever HRT use, which requires further investigation.


Background
Statins are inhibitors of 3-hydroxy-3-methyl-glutaryl coenzyme A (HMG-CoA) reductase and are used worldwide as an efficient cholesterol-lowering medication with proven anti-inflammatory properties [1]. By reducing HMG-CoA reductase activity, levels of mevalonate decrease, which influences important down-stream players for cellular pathways that are crucial for cancer initiation, growth, and metastasis [2]. Pre-clinical studies and phase II trials have demonstrated anti-carcinogenic properties of statins [3][4][5]. Large-scale epidemiological observational studies have reported conflicting findings regarding the role of statins in breast cancer prevention [6][7][8][9][10][11]. On the other hand, available data on breast cancer prognosis are more consistent showing a lower risk of recurrence in lipophilic statin users [12][13][14]. These data suggest that statins may not prevent cancer occurrence, e.g. by affecting the malignant genotype, but rather inhibit the progress of an existing cancer by affecting the malignant phenotype [11]. Mammographic density (MD), which reflects the amount of fibroglandular or radio dense tissue on an X-ray of the breast (mammogram), is a strong determinant of breast cancer risk [15,16]. Women with extremely dense breasts (>75 % dense tissue on a mammogram) have a 4-6 fold higher risk of developing breast cancer compared to women having fatty or non-dense breasts (<5 % dense tissue) [15]. MD also shares many risk factors with breast cancer and is therefore seen as an intermediate in breast cancer etiology [17]. Several determinants of high MD have also been proven to cause cancer [18] and MD may be a useful marker in studies aiming at elucidating the role of potential risk factors for breast cancer [19,20]. To clarify the effect of statins on breast cancer risk, we decided to study the association between statin use and volumetric mammographic density in a large-screening based cohort with comprehensive information on medication use, lifestyle factors and volumetric mammographic density.

Study population
The KARolinska MAmmography project for the risk prediction of breast cancer (KARMA) is a prospective cohort study initiated in January 2011 and comprises 70,877 women who attended mammography screening or clinical mammography at four hospitals in Sweden [21]. Participants responded to a detailed web-based questionnaire, and permission was asked for storage of "for processing" (raw) full-field digital mammograms (FFDM) and linkage to Swedish national registers on inpatient care, prescriptions, cancer, and cause of death.
For the present study, we selected all women who attended the mammography screening program (40-74 years) with raw digital mammograms stored at baseline (N = 50,461). We excluded women with previous cancers other than non-melanoma skin cancer (N = 3,015), women who underwent breast enlargements/reductions/ surgery (N = 2,191), women lacking information on age and body mass index (BMI) (N = 3,908), and women pregnant 12 months prior to study entry (N = 40), leaving 41,307 women in the study. Of these, 41,102 women had medio-lateral oblique (MLO) mammograms and represented the final study population. The study was approved by the ethical review committee at the Karolinska Institute (Stockholm, Sweden), and all participants provided written informed consent.

Statin exposure
Information on statin use was obtained through linkage to the Swedish Prescribed Drug Register [22], which has had nationwide coverage since mid-2005. The register covers all drugs sold and dispensed by prescription at Swedish pharmacies to patients with a personal identification number. The following Anatomical Therapeutic Chemical (ATC) codes were extracted: 'C10AA01' (simvastatin), 'C10AA04' (fluvastatin), 'C10AA05' (atorvastatin) and 'C10AA08' for lipophilic statins, and 'C10AA03' (pravastatin) and 'C10AA07' (rosuvastain) for hydrophilic statins. Women were considered to be current users if they were dispensed statins during the year before mammography. We also calculated and categorized the cumulative statin duration in current users in three-year categories (<2 years, 2-5 years, >5 years).

Measurement of mammographic density
Mammographic density was measured from baseline mammograms using a fully automated volumetric method, Volpara, which has been validated against breast magnetic resonance imaging (MRI) data and appears to be robust to changes in imaging conditions [23]. We have previously shown that Volpara performs well in a high-throughput setting with both percent and absolute dense volumes being associated with established density determinants and breast cancer risk [24]. For our analyses, we used the average mammographic density from the left and right breast of the MLO view.
Technical details of the Volpara algorithm have been described elsewhere [23]. In brief, the algorithm computes the thickness of dense tissue at each pixel using the X-ray attenuation of an entirely fatty region as an internal reference. The absolute dense volume (cm 3 ) is measured by integrating the dense thickness at each pixel over the whole mammogram, and the total breast volume (cm 3 ) is derived by multiplying the breast area by the recorded breast thickness with an appropriate correction for the breast edge. The percent dense volume (%) is obtained from the ratio of these two measures, and the absolute non-dense volume is derived by subtracting the absolute dense volume from the total breast volume.

Covariate assessment
Information on potential confounders was derived from the web-based questionnaire. The questionnaire covers a wide range of breast cancer risk factors and other baseline characteristics such as anthropometrics, alcohol and smoking, education, reproductive and hormonal factors (age at menarche, number of births and age at first birth, menstruation status, use of oral contraceptives and hormone replacement therapy (HRT)), and family history of breast cancer. Body mass index (BMI) was calculated based on self-reported weight and height. Menopausal status was defined according to information on menstruation status, previous oophorectomy, and age at study entry. Women with no periods during the last year, a history of oophorectomy, or age 55 or older were defined as postmenopausal. Information on other medications was extracted from the Swedish Prescribed Drug Register using the following ATC codes: 'B01AC06' (low-dose aspirin) and ' A10BA02' (metformin). Current use was defined by dispensing within the year preceding mammography.

Statistical methods
Characteristics of statin users were compared with those of non-users using Chi-square tests for categorical variables and T-tests for continuous variables. All mammographic measures were log-transformed to approximate the normal distribution. We used analysis of covariance (ANCOVA) to examine associations between statin use and mammographic density, and geometric means were calculated for each category of statin use. We also studied the effect of statin class (lipophilic vs. hydrophilic) and treatment duration. Women using both lipophilic and hydrophilic statins were included in analyses of overall statin use and treatment duration but were excluded from analysis by statin class (N = 46).
All analyses were adjusted for potential confounders in a stepwise manner. We started with an age-adjusted model, after which BMI was added. In the third model, we further adjusted the analyses for other known determinants of mammographic density confounders, including menopausal status, HRT use, parity/age at first birth, age at menarche, education level, smoking, alcohol consumption, and history of benign breast disease. In the final multivariable adjusted model, we added comedication including low-dose aspirin and metformin use.
Since the effect of statins on breast cancer risk may depend on HRT use [25], we also tested for potential effect modification by HRT by adding a multiplicative interaction term between statin and any history of HRT use to the multivariable adjusted models. All statistical analyses were conducted using STATA version 12.0 (Stata Corp. College Station, TX, USA). Seventeen percent of the women had missing values on one or more covariates. On the basis of the missing data pattern, we assumed the missing values to be "missing at random" [26], and imputed these values using a multiple imputation technique (10 imputation sets) [27]. Table 1 shows the characteristics of statin users and nonusers. Statin use was recorded in 3,337 women (8.1 %) of the KARMA study population. As expected, statin users were more likely to be of older age, to be postmenopausal, and to have a higher BMI and a smoking history. Statin users also more often reported HRT use in the past and were more frequently on low-dose aspirin and metformin medication. Other variables associated with statin use were age at menarche, age at first birth, oral contraceptive use, and alcohol intake.

Results
The distribution of statin use by lipophilicity, type, and duration is summarized in Table 2. Most statin users were prescribed a lipophilic statin (93.4 %), with simvastatin being the most commonly prescribed type. Among current users, the percentages of participants using statins <2 years, 2-5 years, and >5 years were 15.1 %, 38.5 % and 46.4 % respectively. Table 3 shows the geometric means of mammographic density by statin use, lipophilicity, and statin duration.
Overall, percent dense volume was lower in statin users than in non-users (P < 0.001). The largest difference in percent dense volume was observed in the age-adjusted model, while the effect of statin use was attenuated after adjustment for BMI, resulting in smaller difference in geometric means despite being significant. Inspection of the individual density components revealed that the inverse association with percent dense volume was attributed to a larger absolute non-dense volume in statin users (P < 0.001). No effect of statin use on the absolute dense volume was found after multivariable adjustment, although a trend towards a higher dense volume was observed after adjusting for concomitant use of low-dose aspirin and metformin (P = 0.06). There was also no evidence of differential effects by lipophilicity and treatment duration according to overlapping confidence intervals. Both lipophilic and hydrophilic statin users had a lower percent dense volume and higher absolute non-dense volume compared to non-users. The differences in density measures by statin duration was driven by the strong overall effect of statin use, without a clear trend by length of exposure in current users as demonstrated by the U-shaped associations.
Interaction analyses showed evidence of effect modification by HRT. A positive association between statin use and the absolute dense volume was found in those who had used HRT, while no such association was found in never HRT users (P interaction = 0.03) ( Table 4).

Discussion
In this study including 41,102 women from a large screening-based cohort, we found no evidence of an overall effect of statin use on absolute dense volume. Statin users had a significantly lower percent dense volume, but this difference was mainly attributable to a larger non-dense volume in statin users. While a large number of studies have investigated the effect of statin use on breast cancer risk and prognosis [6][7][8][9][10][11][12][13][14], only few studies have examined the association with mammographic density. Our results are in line with a longitudinal study showing no effect of statins on change in area-based mammographic density [28]. No differential effects were observed by statin class or treatment duration [13,29] but interaction analyses revealed potential effect modification by HRT with a larger absolute dense volume among statin users who also reported HRT use.
Mammographic density can be expressed in both absolute and relative terms. The preferred MD measure for breast cancer risk prediction is debated, although recent studies show that the proportion of dense breast tissue may have a stronger predictive value than the absolute amount of dense tissue, thus suggesting non-dense tissue as a contributing risk factor [30,31]. Percent dense volume mirrors the relationship between dense and non-dense parenchyma of the breast. Because the proportion of dense volume is strongly influenced by the amount of non-dense tissue in the breast, any factor that is associated with a difference in non-dense volume can lead to a substantial difference in percent dense volume, even if there is no association with the absolute dense volume.
Since most statin users have a medical history of hypercholesterolemia and are more likely to have larger fat deposits than non-users, their lower percent dense volume, as observed in this study, is most likely due to the underlying indication, as reflected by their larger absolute nondense volume. We therefore believe that the absolute Missingness on all variables < 5 %, except for age at first birth (5.6 %) and hormone replacement therapy (6.9 %) dense volume is a more optimal measure for studying the effect of statins on mammographic density, as this measure is less influenced by this type of confounding [32]. The role of statins in primary breast cancer prevention is controversial, and several factors need to be acknowledged when interpreting the available evidence. First of all, our study highlights the issue of confounding by indication, an effect that cannot easily be controlled for in observational studies [33]. Any protective effect of statins may be counteracted in observational analyses due to the underlying indication of dyslipidemia, as high total cholesterol and low HDL cholesterol levels have been associated with an increase in breast cancer risk in postmenopausal women [34,35], although a recent publication report opposite findings with positive associations between HDL and breast cancer risk [36]. The notion that statin use and body composition are strongly interrelated, another effect not easily adjusted for in association analyses, may also partly explain the conflicting observational data regarding the effect of statins on breast cancer risk [6][7][8]13].
Another shortcoming in many studies may be the lack of separate analyses for lipophilic and hydrophilic statins [7]. Theoretically, lipophilic statins, in contrast to hydrophilic statins, are able to penetrate the cell membrane and affect cell proliferation and survival, hence having anticarcinogenic effects [37]. Regarding uptake in the hepatocyte cell, lipophilic statins are passively diffused through membranes, whereas hydrophilic statins have an active carrier-mediated uptake [38]. The different mechanisms of uptake generate less hepato-selectivity in lipophilic statins [38]. Animal studies have shown greater distribution in extrahepatic tissue for lipophilic statins [39]. Thus, lipophilic statins, in comparison to hydrophilic statins, have superior influence on extrahepatic tissue. However, in this study, analyses stratified for statin lipophilicity did not show a significant difference in absolute dense volume. Finally, bioavailability differs among statins, and the distribution to different tissues is dose-dependent [2]. The effectiveness of statins, which were initially designed as hepatoselective drugs [40], has been questioned as a cancerpreventive drug considering its low systemic availability [10]. Recently, we conducted a window-of-opportunity trial, prescribing a high-dose lipophilic statin (atorvastatin) to newly diagnosed breast cancer patients two weeks prior to primary breast cancer surgery. The study showed significantly reduced proliferation (Ki67) in the subgroup of patients where tumors expressed the target for statins, HMG-CoA reductase [4]. These results imply that statins may exert systemic effects including tumor-targeted effects in breast cancer, despite its hepatoselective properties. Similar antiproliferative effect of lipophilic statins have been demonstrated in another window trial [5]. Therefore, statins and their role as anti-cancerous drug, requires further investigation, in particular long-term use of high-dose lipophilic statins [41].
Despite the lack of an overall association with statin use, the higher absolute dense volume found in statin users who also reported HRT is of interest. HRT use is associated with an increase in MD [42] and breast cancer risk [43,44], and the effect of HRT seems to depend on the pretreatment level of endogenous estrogens. Data from the Women's Health Initiative trial show that the adverse effect of HRT on breast cancer risk is largest in those having the lowest estrogen levels [45]. Data on the effect of statins on endogenous estrogen levels are scarce, but a small biomarker study reported lower estrone sulfate levels in women using simvastatin [28]. Given the larger HRT effect in a low endogenous estrogen environment, this may explain the larger absolute dense volume observed with statins in ever HRT users. Moreover, an interaction between statin and HRT use has previously been reported for breast cancer [25], in a direction similar to that observed for the absolute dense volume. An alternative explanation for the observed interaction is the underlying indication. A metabolite of cholesterol has been associated with estrogen positive breast cancers [46], supporting a cross talk between estrogen and cholesterol in breast cancer. Thus, the larger absolute dense volume in statin users who ever used HRT could also reflect the combined effect of HRT and underlying disease of hypercholesterolemia (rather than the actual statin prescription)  on mammographic density. Since mammographic density represents one of the strongest risk factors for breast cancer, potential adverse effects of statins and/or hypercholesterolemia in women using HRT require further exploration.
Our study has several strengths, including the large screening-based cohort with detailed information on lifestyle factors as well as registry-based information on statin use and other medications. We used a fully-automated method for measuring mammographic density, which in contrast to other methods is not prone to subjective measurement error. Moreover, volumetric approaches yield more accurate density measures, as they account for interindividual differences in breast thickness. Both dense and non-dense tissues were evaluated in this study, which strengthens the interpretation of our findings and highlighted the issue of confounding by indication. Our study also has some limitations. The cross-sectional nature limits the ability to explore the temporal association between statins and mammographic density in terms of initiation and discontinuation. Also, we could not study the effect of long-term treatment (>10 years), because the population-based prescription register has only reached full coverage since 2005. Finally, we did not have data available on the dosages prescribed.

Conclusions
In conclusion, we found no strong evidence for an overall effect of statin use on volumetric mammographic density in terms of absolute dense volume. The observed interaction with HRT with a larger absolute dense volume in statin users, who also reported HRT use, requires further investigation. All mammographic density measures were log transformed prior to analysis and the values shown are on a back-transformed scale. b In postmenopausal women only. P value refers to F-test for between-statin group differences in geometric means. P interaction refers to the P value of the multiplicative interaction term Model 4: Adjusted for age (years), body mass index (kg/m 2 ), menopausal status (pre, peri/unknown, post), HRT use (never, former, current), parity/age at first birth (nulliparous, 1 child age at first birth < 25 years, 1 child age at first birth ≥ 25 years, 2 children age at first birth < 25 years, 2 children age at first birth ≥ 25 years, ≥ 3 children age at first birth < 25 years, ≥ 3 children age at first birth ≥ 25 years), age at menarche (years), education level (compulsory, gymnasium, university), smoking (never, former, current), alcohol consumption (nondrinker, 0,1-1,9 g/day, 2.0-4.9 g/day, 5.0-9.9 g/day, ≥10 g/day) and benign breast disease (yes vs. no), low-dose aspirin use (yes vs. no) and metformin use (yes vs. no) Submit your next manuscript to BioMed Central and take full advantage of: