Skip to main content

Long-term exposure to insulin and volumetric mammographic density: observational and genetic associations in the Karma study

Abstract

Background

Long-term insulin exposure has been implicated in breast cancer etiology, but epidemiological evidence remains inconclusive. The aims of this study were to investigate the association of insulin therapy with mammographic density (MD) as an intermediate phenotype for breast cancer and to assess associations with long-term elevated circulating insulin levels using a genetic score comprising 18 insulin-associated variants.

Methods

We used data from the KARolinska MAmmography (Karma) project, a Swedish mammography screening cohort. Insulin-treated patients with type 1 (T1D, n = 122) and type 2 (T2D, n = 237) diabetes were identified through linkage with the Prescribed Drug Register and age-matched to 1771 women without diabetes. We assessed associations with treatment duration and insulin glargine use, and we further examined MD differences using non-insulin-treated T2D patients as an active comparator. MD was measured using a fully automated volumetric method, and analyses were adjusted for multiple potential confounders. Associations with the insulin genetic score were assessed in 9437 study participants without diabetes.

Results

Compared with age-matched women without diabetes, insulin-treated T1D patients had greater percent dense (8.7% vs. 11.4%) and absolute dense volumes (59.7 vs. 64.7 cm3), and a smaller absolute nondense volume (615 vs. 491 cm3). Similar associations were observed for insulin-treated T2D, and estimates were not materially different in analyses comparing insulin-treated T2D patients with T2D patients receiving noninsulin glucose-lowering medication. In both T1D and T2D, the magnitude of the association with the absolute dense volume was highest for long-term insulin therapy (≥ 5 years) and the long-acting insulin analog glargine. No consistent evidence of differential associations by insulin treatment duration or type was found for percent dense and absolute nondense volumes. Genetically predicted insulin levels were positively associated with percent dense and absolute dense volumes, but not with the absolute nondense volume (percentage difference [95% CI] per 1-SD increase in insulin genetic score = 0.8 [0.0; 1.6], 0.9 [0.1; 1.8], and 0.1 [− 0.8; 0.9], respectively).

Conclusions

The consistency in direction of association for insulin treatment and the insulin genetic score with the absolute dense volume suggest a causal influence of long-term increased insulin exposure on mammographic dense breast tissue.

Background

The role of insulin in breast cancer etiology has received growing attention in recent years [1]. Basic research suggests that long-term exposure to elevated exogenous and endogenous insulins promotes breast tumor growth, either directly by signaling mitogenic effects through the insulin receptor isoform A and the insulin-like growth factor 1 (IGF-1) receptor [2, 3] or indirectly by altering the levels of circulating estrogens [4]. The potential carcinogenic effect of insulin has been demonstrated in vitro in terms of increased proliferation in human breast epithelial cells and breast cancer cell lines [2, 3]. Whether these in vitro observations are relevant to humans and concerns surrounding risks of long-term exogenous insulin use and elevated circulating insulin levels are justified remains uncertain [5,6,7,8].

Thus far, observational studies [6, 9] and randomized clinical trials [10,11,12] have found no compelling evidence of an increased breast cancer incidence in insulin-treated diabetes patients, although recently some studies have suggested a possible elevated risk with long-term insulin glargine use [5, 13, 14], which might be related to the more constant pharmacokinetic profile and enhanced IGF-1 receptor affinity of this insulin analog [15, 16]. However, because of methodological limitations including short-term follow-up and confounding by indication, definitive conclusions cannot be drawn. Many agents affecting carcinogenesis have long latencies and require a minimum length of exposure. Most studies addressing recency or duration effects had limited control for confounding factors or were insufficiently powered to assess associations with long-term insulin therapy [6]. While randomized clinical trials are limited in size and follow-up to study cancer-specific outcomes, confounding by indication is a concern in observational studies, and the use of different comparators is necessary to disentangle treatment effects from their underlying indications. Studies investigating associations of circulating insulin levels with breast cancer risk have also yielded conflicting results, with either positive [17,18,19] or null associations [20]. Most of these findings, however, were based on small numbers of breast cancer patients and a single insulin measurement, which is not an ideal proxy for long-term insulin exposure 18].

Mammographic density (MD) refers to the amount of radiologically dense fibroglandular tissue in the breast, and high MD levels are a strong and independent predictor of breast cancer risk [21, 22]. Both traits also have several reproductive and lifestyle determinants in common, and MD is viewed as an intermediate phenotype in breast cancer etiology [23]. Because many insulin-treated diabetes patients are below the age at which breast cancer is usually diagnosed, and given the long latency of breast cancer, MD serves as an attractive intermediate endpoint for identifying potential carcinogenic effects.

In the present study, we aimed to investigate the association of long-term insulin exposure with MD in a mammography screening cohort using different methodological approaches. First, we assessed associations of insulin therapy with MD by comparing insulin-treated type 1 (T1D) and type 2 (T2D) diabetes patients with age-matched individuals without diabetes, overall and stratified by duration of insulin treatment and insulin glargine use. Potential confounding by indication was addressed in case-only analyses and additional analyses comparing insulin-treated diabetes patients with patients receiving other glucose-lowering medication. Second, as an alternative means of overcoming confounding, we explored associations with an insulin genetic score in nondiabetic women. This score comprising 18 insulin-associated variants represents genetic predisposition to elevated circulating insulin levels over the life course. Because genotypes sort randomly at conception, genetic association analyses are less likely affected by confounding and hence can provide additional evidence of the likelihood of a causal effect of long-term increased insulin exposure [24].

Methods

Study populations

This study was nested within the KARolinska MAmmography Project for Risk Prediction of Breast Cancer (Karma), a prospective cohort of 70,877 women attending mammography screening or clinical mammography at one of four mammography units in Sweden between January 2011 and March 2013 [25]. All participants responded to a web-based questionnaire and ~ 55,000 women without a history of cancer, breast enlargement, reduction, or surgery had raw digital mammograms collected and stored at study entry, representing the study base for the present analysis. Information on diabetes diagnoses and insulin prescriptions was retrieved through linkage with the Swedish Patient Register [26] and Prescribed Drug Register [27]. The Patient Register has nationwide coverage and includes inpatient hospitalizations since 1987 and outpatient physician visits since 2001. The Prescribed Drug Register covers all drugs sold and dispensed by prescription since July 1, 2005.

Associations between insulin therapy and MD were analyzed in a matched cohort design including all insulin-treated T1D and T2D patients and a sample of age-matched individuals without diabetes. Insulin-treated diabetes patients were identified through the Prescribed Drug Register, with current insulin use defined as at least one dispensed prescription for insulin or insulin analogs (Anatomical Therapeutic Chemical Classification System [ATC] code A10A, including A10AE04 for insulin glargine) in the year prior to study entry (see Additional file 1: Table S1 for all identified insulin prescriptions). Diabetes diagnoses were retrieved from the Karma questionnaire and Patient Register. All register-based diagnoses were based on International Classification of Diseases (ICD) code 250 (ICD-8 and ICD-9) until 1996 and unique codes for T1D and T2D as introduced from 1997 onward (ICD-10 codes E10 and E11, respectively). Because no distinction between T1D and T2D was made in the Karma questionnaire and earlier ICD versions (ICD-8 and ICD-9), T1D and T2D patients were differentiated on the basis of previously established cutoffs of diagnosis age (T1D, ≤ 30 years; T2D, ≥ 40 years) [28, 29] when diabetes-specific codes were missing. Using these prescription and diagnostic criteria, 359 insulin-treated diabetes patients were identified, including 122 T1D and 237 T2D patients (of whom 97 T1D [79.5%] and 76 T2D [32.1%] patients were differentiated on the basis of diagnosis age). For each patient, we randomly sampled up to five individuals without a diabetes diagnosis from the study base, matched on birth year. In total, 21 T2D patients could not be matched to a maximum of 5 individuals, leaving 1771 age-matched individuals without diabetes for analyses (i.e., 610 and 1161 for insulin-treated T1D and T2D, respectively).

We further evaluated insulin treatment effects in an analysis comparing T2D patients treated with insulin only (n = 112) with T2D patients receiving other noninsulin glucose-lowering medication (n = 407). This active comparator group comprised all T2D patients with at least one dispensed prescription for glucose-lowering medication (ATC code A10B), excluding insulins in the year prior to study entry (see Additional file 1: Table S2 for all identified noninsulin glucose-lowering prescriptions). Because insulin is the mainstay of therapy for T1D, we were unable to do a similar analysis for this group of patients.

Finally, we examined associations with an insulin genetic score comprising 18 insulin-associated variants (see details below) as an instrument to proxy long-term exposure to circulating insulin levels. Associations with the insulin genetic score were assessed in a subcohort of 9437 participants with available genotyping data. All women in the subcohort had no history of cancer or diabetes at the time of study entry when blood samples were obtained. The Karma study was approved by the ethical review committee at Karolinska Institutet, and all participants provided written informed consent.

Volumetric mammographic density

MD was estimated from raw digital mammograms collected at study entry using Volpara ™ version 1.5.0 (Volpara Solutions, Wellington, New Zealand) [30]. Volpara volumetric MD measures show good agreement with breast magnetic resonance imaging data [31] and have been validated as being predictive of breast cancer risk [30, 32]. The Volpara algorithm estimates 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 (cm3) is computed by integrating the dense thickness at each pixel over the whole mammogram, and the total breast volume (cm3) is derived by multiplying the breast area by the recorded breast thickness, with an appropriate correction for the breast edge. From these measures, the absolute nondense volume (cm3) and percentage of the breast covered by dense tissue (%) can be obtained. The average measurement of the left and right breasts of the mediolateral oblique view was taken for all analyses.

Covariates

The following potential confounders known to be associated with MD were extracted from the baseline questionnaire: education level, body mass index (BMI, based on self-reported height and weight), lifestyle measures (smoking, alcohol intake, and physical activity), reproductive and hormonal factors (age at menarche, number of births/age at first birth, menopausal status, use of oral contraceptives and hormone replacement therapy), and personal history of previous benign breast disease and breast cancer heredity. We also extracted information on prescriptions of comedication through linkage with the Prescribed Drug Register using the ATC coding system (including low-dose aspirin [ATC code B01AC06], statins [ATC codes C10AA01, C10AA03, C10AA04, C10AA05, C10AA07, C10AA08], and metformin [ATC code A10BA02]) and summarized data on comorbid conditions derived from the Patient Register into the Charlson comorbidity index score [33].

Insulin genetic score

Associations with the insulin genetic score were assessed in the Karma subcohort with genotyping data. Whole-blood samples of 9437 study participants were genotyped using Illumina iSelect arrays (iCOGS [n = 3909] and OncoArray [n = 5528]), details of which have been described elsewhere [34, 35], and missing genotypes for common variants across the genome were imputed using the 1000 Genomes Project March 2012 release as a reference. The insulin genetic score was constructed using 18 independent single-nucleotide polymorphisms (SNPs) shown to be robustly associated with circulating log insulin levels at P < 5.0 × 10− 8 [36] (Additional file 1: Table S3). The score was calculated on the basis of a weighted method according to each SNP’s effect size (β) obtained from the literature [24]:

$$ \mathrm{Insulin}\ \mathrm{genetic}\ \mathrm{score}={\beta}_1{x}_1+{\beta}_2{x}_2+.\dots \kern0.5em {\beta}_k{x}_k+{\beta}_n{x}_n, $$

where β is the per-allele beta value for log-transformed fasting insulin levels associated with the effect allele for SNP k and xk is the number of alleles for the same SNP (0, 1, 2) and n is the total number of SNPs included in the score (n = 18).

Statistical analyses

To approximate the normal distribution, all mammographic measures were log-transformed prior to analyses, and geometric means and percentage (%) differences were calculated [37]. Differences in MD by insulin therapy were analyzed using generalized linear models, including a basic model with adjustment for age and BMI and a multivariable model with inclusion of other potential confounders (education level, age at menarche, parity and age at first birth, oral contraceptives, menopausal status, hormone replacement therapy, alcohol intake, smoking status, statins, low-dose aspirin, Charlson comorbidity index, history of benign breast disease, and family history of breast cancer). Analyses were additionally adjusted for metformin comedication to account for potential antiproliferative effects of this diabetes drug [6]. Differences in MD were assessed overall and by insulin glargine use and treatment duration. In T2D patients, treatment duration was defined from the first dispensed insulin prescription encountered in the Prescribed Drug Register. Because age at T1D onset is well below the age at which women undergo mammography screening, and with prescription data being available only from July 2005 onward, insulin treatment duration in T1D patients was calculated from the age at T1D diagnosis to study entry.

To address possible residual confounding by underlying disease, we also examined MD differences comparing insulin-treated with non-insulin-treated T2D patients. This active comparator analysis was adjusted for the same covariates as listed above and additionally for diabetes duration to account for differences in disease onset between the two patient groups.

Associations with the insulin genetic score were assessed using linear regression, with beta values representing percentage differences in MD per 1-SD increment in the score. Because there was no evidence of heterogeneity by genotyping array (iCOGS vs. OncoArray), a one-sample approach was undertaken. Genetic score analyses were adjusted for age, BMI, menopausal status, genotyping array, and six principal components to account for population stratification. To investigate the independence of genetic effects of other potential confounders, we also examined associations of the score with covariates entered in the insulin treatment analysis.

All analyses were undertaken using STATA version 14 (StataCorp, College Station, TX, USA) and PLINK version 1.9 [38]. Missing values on covariates were imputed using multivariate multiple imputation with chained equations, and ten imputed datasets were generated [39]. Imputation models included the outcome (MD), exposure (insulin treatment/insulin genetic score), and all covariates included in any of the analysis models (see Additional file 2: Supplementary Methods).

Results

Descriptive characteristics of the study populations are summarized in Table 1. Mean ages at diagnosis were 20 and 55 years for the insulin-treated T1D and T2D patients, respectively. Insulin glargine was prescribed to 53% of T1D patients and 31% of T2D patients, and combined therapy with metformin was given in 4% and 53% of T1D and T2D patients, respectively. Compared with age-matched individuals without diabetes, insulin-treated T1D and T2D patients were younger at first child’s birth, less often used oral contraceptives in the past, more frequently reported a family history of breast cancer, and were more often alcohol abstainers and less physically active at study entry. As expected, insulin-treated T1D and T2D patients also presented with more comorbid conditions and were more often on statin and low-dose aspirin medications. Univariate associations with other participant characteristics were comparable for the two patient groups, except for BMI and age at menarche. Study participants of the subcohort for genetic analysis had characteristics similar to those of the age-matched individuals without T2D, except for a larger proportion being premenopausal. Descriptive characteristics of insulin-treated and non-insulin-treated T2D patients are summarized in Additional file 1: Table S4. Compared with T2D patients receiving other glucose-lowering medication, insulin-treated T2D patients were younger, had a lower BMI, and were more likely to have comorbid conditions and a history of benign breast disease. They also had a longer disease duration (8.3 vs. 4.8 years), as reflected by the younger age at diagnosis.

Table 1 Descriptive characteristics of the study population

A summary plot of MD percentage differences estimated in each observational analysis is presented in Additional file 3: Figure S1 to facilitate comparison across the different analyses. Geometric means of MD in insulin-treated diabetics and age-matched individuals without diabetes are listed in Table 2 together with corresponding percentage differences. Compared with age-matched women without diabetes, T1D patients had a greater percent dense (11.3% vs. 8.7%) and absolute dense (65.7 vs. 59.6 cm3) volume and a smaller absolute nondense volume (510 vs. 610 cm3). These associations were not materially different after multivariable adjustment (percent dense volume [11.4 vs. 8.7%], absolute dense volume [64.7 vs. 59.7 cm3], and absolute nondense volume [491 vs. 615 cm3]). Compared with women without diabetes, the largest difference in absolute dense volume was found for current use of insulin glargine, whereas for percent dense and absolute nondense volumes, no notable difference in magnitude of association was observed by insulin type. Overall, similar associations were found for insulin-treated T2D. Compared with age-matched individuals without diabetes, insulin-treated T2D patients had greater percent dense (7.8% vs. 6.5%) and absolute dense (58.0 vs. 52.6 cm3) volumes, as well as a smaller absolute nondense volume (754 vs. 679 cm3). As for T1D, insulin glargine use was associated with the largest difference in absolute dense volume in T2D patients.

Table 2 Geometric means and percentage differences of volumetric mammographic density comparing insulin-treated T1D and T2D patients with age-matched individuals without diabetes

A positive association of insulin glargine use with the absolute dense volume was also observed in case-only analyses including insulin-treated T1D and T2D patients only (Additional file 1: Table S5), despite the lower level of statistical significance due to smaller sample size. The association of insulin glargine with the absolute nondense volume, however, was directionally inconsistent with the association observed in analyses using age-matched nondiabetics as a comparator (Additional file 1: Table S5, Additional file 3: Figure S1). In case-only analyses, glargine insulin users had a greater (T2D) or similar-sized (T1D) absolute nondense volume (and consequently similar-sized (T2D) or greater (T1D) percent dense volume) compared with nonglargine insulin users. Associations between insulin therapy and MD were not very different in analyses comparing insulin-treated T2D patients with T2D patients receiving noninsulin glucose-lowering medication; that is, insulin-treated T2D patients had a greater percent and absolute dense volume and a smaller absolute nondense volume than non-insulin-treated T2D patients, with the absolute dense volume association being specific for insulin glargine use (Additional file 1: Table S6).

Results of analyses by insulin treatment duration are summarized in Fig. 1 and Additional file 1: Table S7. Compared with age-matched individuals without diabetes, a gradual increase in absolute dense volume was found with duration of insulin treatment in T1D patients. In T2D patients, a similar increase in absolute dense volume was observed with treatment duration. Case-only analyses in insulin-treated T1D patients supported the presence of a treatment duration effect on the absolute dense volume (P trend = 0.04), though numbers were small, and the association with long-term therapy (beyond 28 years) attenuated following multivariable adjustment (P trend = 0.23) (Additional file 1: Table S8). No consistent evidence of a treatment duration effect was found for percent dense and absolute nondense volumes, overall and in case-only analyses (Additional file 3: Figure S1, Additional file 1: Tables S7 and S8).

Fig. 1
figure 1

Associations of duration of insulin therapy with volumetric mammographic density in type 1 diabetes (T1D) and type 2 diabetes (T2D) patients. Geometric means and 95% CIs of volumetric mammographic density by duration of insulin therapy in T1D and T2D patients. Model 1 (open circles): adjusted for age and body mass index. Model 2 (closed circles): adjusted for age, body mass index, education level, age at menarche, parity, age at first birth, menopausal status, oral contraceptives use, hormone replacement therapy, alcohol intake, physical activity, smoking status, statins, low-dose aspirin, Charlson comorbidity index, benign breast disease, family history of breast cancer, and metformin therapy

Associations of the insulin genetic score with MD are shown in Fig. 2, and a summary of individual SNP estimates is provided in Additional file 1: Table S9. Genetically predicted insulin levels were positively associated with both percent dense and absolute dense volume (Fig. 2) (beta [95% CI] per 1-SD increase in genetic score = 0.80 [0.00–1.59] and 0.93 [0.05–1.81], respectively), but not with the absolute nondense volume (beta [95% CI] per 1-SD increase in genetic score = 0.05 [− 0.78; 0.87]). There was no evidence of associations being driven by potential confounders (Additional file 4: Figure S2), except for genetically elevated insulin levels being weakly associated with a lower odds of previous benign breast disease. Given the positive association between benign breast disease and MD, analyses conditioning on this variable strengthened the effect estimates (Additional file 5: Figure S3).

Fig. 2
figure 2

Association of the insulin genetic score with volumetric mammographic density. Association of 18-single-nucleotide polymorphism insulin genetic score with volumetric mammographic density in nondiabetic women, overall and stratified by genotyping array. Associations with volumetric mammographic density were analyzed in a linear regression model, adjusted for age, body mass index, menopausal status, and six principal components. Analyses in the total study population were additionally adjusted for genotyping array. All volumetric mammographic density measures were log-transformed prior to analyses, with beta values representing percentage differences in volumetric mammographic density per 1-SD increase in insulin genetic score

Discussion

In this mammography screening cohort, we found that insulin-treated T1D and T2D patients had higher MD levels as compared with age-matched individuals without diabetes and diabetes patients receiving other noninsulin glucose-lowering medication. Associations with the absolute dense volume were driven mainly by long-term insulin use and the long-acting insulin analog glargine, whereas no associations with treatment duration or insulin type were found for the absolute nondense volume. We further observed positive associations of genetically predicted insulin levels with percent dense and absolute dense volumes, but not with absolute nondense volume. The consistency in direction of association for insulin treatment and the insulin genetic score with the absolute dense volume strengthen the evidence of a potential causal effect of long-term increased insulin exposure on mammographic dense breast tissue.

Observational studies and randomized clinical trials are limited in their ability to investigate associations of long-term insulin exposure with relatively rare outcomes such as breast cancer. As a continuous intermediate trait, MD represents an attractive endpoint for identifying potential carcinogenic effects [40]. Because MD was routinely collected in all study participants regardless of screening outcome, this outcome is also not susceptible to ascertainment bias. To date, only one study assessed associations between exogenous insulin use and MD. This study [41] reported a suggestive association of insulin therapy with a higher prevalence of mixed/dense breast density patterns, but it included only 20 diabetes patients and could not examine associations with treatment duration or insulin glargine use, nor could it distinguish associations based on fibroglandular dense and adipose nondense tissues. Though previous observational data on breast cancer risk have been somewhat conflicting, recent findings [14] of a positive association between long-term insulin glargine use and breast cancer risk cohere with the results for the absolute dense volume reported herein. Altogether, these findings suggest that the carcinogenic potential of exogenous insulins might be greatest for insulin glargine, possibly because of its unique characteristics in terms of receptor affinity and pharmacokinetic profile with prolonged duration of action.

Because our study is observational in nature, we aimed to integrate evidence from different methodological approaches to assess the likelihood of a potential causal effect of long-term insulin exposure. First, we investigated the influence of exogenous insulin use in diabetes patients. To address confounding by indication, we assessed associations with insulin-treated T1D and T2D separately. Because T1D and T2D differ in pathophysiology and underlying risk factors (with T1D being an autoimmune disease, whereas insulin resistance, driven mainly by obesity, is the hallmark of T2D), consistent results for insulin-treated T1D and T2D are suggestive of an insulin therapy effect independent of underlying disease etiology. This analysis approach using nondiabetics as a comparator, however, does not rule out confounding by the indication itself (T1D or T2D). Hence, to address possible residual confounding, we also assessed MD differences by duration of insulin treatment and insulin glargine use in case-only analyses, and we performed additional analyses using non-insulin-treated diabetes patients as an active comparator. Although none of these observational assessments may be completely free of bias, sources and directions of bias in each of these are likely to be different. Therefore, consistency of direction of association across the different approaches can be interpreted as evidence for a potential causal association. To further investigate the likelihood of a long-term insulin effect, we also explored associations with an insulin genetic score [24] as an instrument for long-term exposure to elevated circulating insulin levels. Although effect sizes for this endogenous genetic proxy are not directly comparable to those observed for exogenous insulins, because of differences in measurement scale and magnitudes of anticipated physiological effects [13], the consistency in direction of associations observed for the absolute dense volume strengthen the evidence of a causal influence of long-term elevated insulin levels on mammographic dense tissue. On the other hand, direction of association for insulin treatment and the insulin genetic score were not consistent for the absolute nondense volume, arguing against a causal influence of insulin on adipose breast tissue.

To our knowledge, this is the first study investigating associations of long-term insulin exposure with MD combining data from observational and genetic analyses. Although absolute mammographic dense tissue is a well-known risk factor and intermediate phenotype of breast cancer, further studies investigating the effect of  long-term insulin exposure on breast cancer risk are warranted, ideally with MD measurement to address the extent to which breast cancer risk associations are mediated by MD. Because long-term exogenous insulin use in diabetes patients may have a greater impact on breast tissue than genetic predisposition to long-term elevations in endogenous insulin, it will also be relevant to assess potential differential associations by exposure type. Moreover, because insulin is the mainstay of treatment for T1D and uncontrolled T2D, and because diabetes patients tend to participate less in mammography screening programs [42, 43], increasing screening awareness and participation among insulin-treated patients with diabetes may be a first step in reducing insulin-associated adverse effects from a clinical perspective.

Some limitations of the present study are noteworthy. Historical information on insulin treatment duration was not complete, because the Prescribed Drug Register has had nationwide coverage since July 2005. This limited the analysis contrasting short- vs. long-term effects of insulin use, especially in insulin-treated T2D patients, where no assumption regarding treatment initiation could be made. Also, insulin glargine and treatment duration analyses were limited by small numbers of patients, which resulted in some uncertainty in effect estimates and low statistical power when analyses were restricted to diabetes patients only. We also cannot rule out misclassification of T1D and T2D patients in instances where diagnostic age cutoffs were used. T2D incidence below the age of 30 years was low in Sweden during the study period, and most diabetes patients diagnosed before 30 years of age are likely to be true T1D cases. Finally, our study was insufficiently powered to study effects of insulin analogs other than insulin glargine. Strengths of the current study are the screening-based setting, extensive information on potential confounders, and the use of a fully automated method for MD measurement that is not prone to subjective measurement error. The unique study design further allowed us to address the likelihood of a causal association by making relevant patient comparisons and by using a genetic score robustly associated with circulating insulin levels independent of other metabolic markers and confounding factors [24].

Conclusions

Our study provides evidence of observational and genetic associations of long-term increased insulin exposure with the absolute dense volume. Apart from identifying a potential causal effect of long-term increased insulin exposure on mammographic dense breast tissue, these findings support efforts to improve screening awareness and participation among insulin-treated patients with diabetes.

Abbreviations

ATC:

Anatomical Therapeutic Chemical Classification System

BMI:

Body mass index

ICD:

International Classification of Diseases

iCOGS:

Illumina iSelect array

IGF-1:

Insulin like-growth factor 1

Karma:

KARolinska MAmmography Project for Risk Prediction of Breast Cancer

MD:

Mammographic density

SNP:

Single-nucleotide polymorphism

T1D:

Type 1 diabetes

T2D:

Type 2 diabetes

References

  1. Giovannucci E, Harlan DM, Archer MC, Bergenstal RM, Gapstur SM, Habel LA, Pollak M, Regensteiner JG, Yee D. Diabetes and cancer: a consensus report. Diabetes Care. 2010;33(7):1674–85.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Chappell J, Leitner JW, Solomon S, Golovchenko I, Goalstone ML, Draznin B. Effect of insulin on cell cycle progression in MCF-7 breast cancer cells: direct and potentiating influence. J Biol Chem. 2001;276(41):38023–8.

    PubMed  CAS  Google Scholar 

  3. Ish-Shalom D, Christoffersen CT, Vorwerk P, Sacerdoti-Sierra N, Shymko RM, Naor D, De Meyts P. Mitogenic properties of insulin and insulin analogues mediated by the insulin receptor. Diabetologia. 1997;40(Suppl 2):S25–31.

    Article  PubMed  CAS  Google Scholar 

  4. Rose DP, Vona-Davis L. The cellular and molecular mechanisms by which insulin influences breast cancer risk and progression. Endocr Relat Cancer. 2012;19(6):R225–41.

    Article  PubMed  CAS  Google Scholar 

  5. Wu JW, Filion KB, Azoulay L, Doll MK, Suissa S. Effect of long-acting insulin analogs on the risk of cancer: a systematic review of observational studies. Diabetes Care. 2016;39(3):486–94.

    Article  PubMed  CAS  Google Scholar 

  6. Bronsveld HK, ter Braak B, Karlstad O, Vestergaard P, Starup-Linde J, Bazelier MT, De Bruin ML, de Boer A, Siezen CL, van de Water B, et al. Treatment with insulin (analogues) and breast cancer risk in diabetics; a systematic review and meta-analysis of in vitro, animal and human evidence. Breast Cancer Res. 2015;17:100.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  7. Formica V, Tesauro M, Cardillo C, Roselli M. Insulinemia and the risk of breast cancer and its relapse. Diabetes Obes Metab. 2012;14(12):1073–80.

    Article  PubMed  CAS  Google Scholar 

  8. Pisani P. Hyper-insulinaemia and cancer, meta-analyses of epidemiological studies. Arch Physiol Biochem. 2008;114(1):63–70.

    Article  PubMed  CAS  Google Scholar 

  9. Larsson SC, Mantzoros CS, Wolk A. Diabetes mellitus and risk of breast cancer: a meta-analysis. Int J Cancer. 2007;121(4):856–62.

    Article  PubMed  CAS  Google Scholar 

  10. Bordeleau L, Yakubovich N, Dagenais GR, Rosenstock J, Probstfield J, Chang Yu P, Ryden LE, Pirags V, Spinas GA, Birkeland KI, et al. The association of basal insulin glargine and/or n-3 fatty acids with incident cancers in patients with dysglycemia. Diabetes Care. 2014;37(5):1360–6.

    Article  PubMed  CAS  Google Scholar 

  11. Rosenstock J, Fonseca V, McGill JB, Riddle M, Halle JP, Hramiak I, Johnston P, Davis M. Similar risk of malignancy with insulin glargine and neutral protamine Hagedorn (NPH) insulin in patients with type 2 diabetes: findings from a 5 year randomised, open-label study. Diabetologia. 2009;52(9):1971–3.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  12. Sturmer T, Marquis MA, Zhou H, Meigs JB, Lim S, Blonde L, Macdonald E, Wang R, Lavange LM, Pate V, et al. Cancer incidence among those initiating insulin therapy with glargine versus human NPH insulin. Diabetes Care. 2013;36(11):3517–25.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. Suissa S, Azoulay L, Dell’Aniello S, Evans M, Vora J, Pollak M. Long-term effects of insulin glargine on the risk of breast cancer. Diabetologia. 2011;54(9):2254–62.

    Article  PubMed  CAS  Google Scholar 

  14. Wu JW, Azoulay L, Majdan A, Boivin JF, Pollak M, Suissa S. Long-term use of long-acting insulin analogs and breast cancer incidence in women with type 2 diabetes. J Clin Oncol. 2017;35(32):3647–53.

    Article  PubMed  Google Scholar 

  15. Kurtzhals P, Schaffer L, Sorensen A, Kristensen C, Jonassen I, Schmid C, Trub T. Correlations of receptor binding and metabolic and mitogenic potencies of insulin analogs designed for clinical use. Diabetes. 2000;49(6):999–1005.

    Article  PubMed  CAS  Google Scholar 

  16. Yehezkel E, Weinstein D, Simon M, Sarfstein R, Laron Z, Werner H. Long-acting insulin analogues elicit atypical signalling events mediated by the insulin receptor and insulin-like growth factor-I receptor. Diabetologia. 2010;53(12):2667–75.

    Article  PubMed  CAS  Google Scholar 

  17. Gunter MJ, Hoover DR, Yu H, Wassertheil-Smoller S, Rohan TE, Manson JE, Li J, Ho GY, Xue X, Anderson GL, et al. Insulin, insulin-like growth factor-I, and risk of breast cancer in postmenopausal women. J Natl Cancer Inst. 2009;101(1):48–60.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  18. Kabat GC, Kim M, Caan BJ, Chlebowski RT, Gunter MJ, Ho GY, Rodriguez BL, Shikany JM, Strickler HD, Vitolins MZ, et al. Repeated measures of serum glucose and insulin in relation to postmenopausal breast cancer. Int J Cancer. 2009;125(11):2704–10.

    Article  PubMed  CAS  Google Scholar 

  19. Lawlor DA, Smith GD, Ebrahim S. Hyperinsulinaemia and increased risk of breast cancer: findings from the British Women’s Heart and Health study. Cancer Causes Control. 2004;15(3):267–75.

    Article  PubMed  Google Scholar 

  20. Mink PJ, Shahar E, Rosamond WD, Alberg AJ, Folsom AR. Serum insulin and glucose levels and breast cancer incidence: the atherosclerosis risk in communities study. Am J Epidemiol. 2002;156(4):349–52.

    Article  PubMed  Google Scholar 

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

    Article  Google Scholar 

  22. Wanders JO, Holland K, Veldhuis WB, Mann RM, Pijnappel RM, Peeters PH, van Gils CH, Karssemeijer N. Volumetric breast density affects performance of digital screening mammography. Breast Cancer Res Treat. 2017;162(1):95–103.

    Article  PubMed  Google Scholar 

  23. 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 

  24. Nead KT, Sharp SJ, Thompson DJ, Painter JN, Savage DB, Semple RK, Barker A, Australian National Endometrial Cancer Study Group (ANECS), Perry JR, Attia J, et al. Evidence of a causal association between insulinemia and endometrial cancer: a Mendelian randomization analysis. J Natl Cancer Inst. 2015;107(9):djv178.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. Gabrielson M, Eriksson M, Hammarström M, Borgquist S, Leifland K, Czene K, Hall P. Cohort profile: the Karolinska Mammography Project for Risk Prediction of Breast Cancer (KARMA). Int J Epidemiol. 2017;46(6):1740–1g.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Ludvigsson JF, Andersson E, Ekbom A, Feychting M, Kim JL, Reuterwall C, Heurgren M, Olausson PO. External review and validation of the Swedish national inpatient register. BMC Public Health. 2011;11:450.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Wettermark B, Hammar N, Fored CM, Leimanis A, Otterblad Olausson P, Bergman U, Persson I, Sundstrom A, Westerholm B, Rosen M. The new Swedish prescribed drug register—opportunities for pharmacoepidemiological research and experience from the first six months. Pharmacoepidemiol Drug Saf. 2007;16(7):726–35.

    Article  PubMed  Google Scholar 

  28. Zendehdel K, Nyren O, Ostenson CG, Adami HO, Ekbom A, Ye W. Cancer incidence in patients with type 1 diabetes mellitus: a population-based cohort study in Sweden. J Natl Cancer Inst. 2003;95(23):1797–800.

    Article  PubMed  Google Scholar 

  29. Hemminki K, Liu X, Forsti A, Sundquist J, Sundquist K, Ji J. Subsequent type 2 diabetes in patients with autoimmune disease. Sci Rep. 2015;5:13871.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Brand JS, Czene K, Shepherd JA, Leifland K, Heddson B, Sundbom A, Eriksson M, Li J, Humphreys K, Hall P. Automated measurement of volumetric mammographic density: a tool for widespread breast cancer risk assessment. Cancer Epidemiol Biomark Prev. 2014;23(9):1764–72.

    Article  Google Scholar 

  31. Highnam R, Brady M, Yaffe MJ, Karssemeijer N, Harvey J. Robust breast composition measurement - VolparaTM. Lect Notes Comput Sci. 2010;6136:342–9.

    Article  Google Scholar 

  32. Eng A, Gallant Z, Shepherd J, McCormack V, Li J, Dowsett M, Vinnicombe S, Allen S, dos -Santos-Silva I. Digital mammographic density and breast cancer risk: a case-control study of six alternative density assessment methods. Breast Cancer Res. 2014;16(5):439.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Charlson ME, Pompei P, Ales KL, MacKenzie CR. A new method of classifying prognostic comorbidity in longitudinal studies: development and validation. J Chronic Dis. 1987;40(5):373–83.

    Article  PubMed  CAS  Google Scholar 

  34. Michailidou K, Hall P, Gonzalez-Neira A, Ghoussaini M, Dennis J, Milne RL, Schmidt MK, Chang-Claude J, Bojesen SE, Bolla MK, et al. Large-scale genotyping identifies 41 new loci associated with breast cancer risk. Nat Genet. 2013;45(4):353–61. e2

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  35. Amos CI, Dennis J, Wang Z, Byun J, Schumacher FR, Gayther SA, Casey G, Hunter DJ, Sellers TA, Gruber SB, et al. The OncoArray Consortium: a network for understanding the genetic architecture of common cancers. Cancer Epidemiol Biomark Prev. 2017;26(1):126–35.

    Article  Google Scholar 

  36. Scott RA, Lagou V, Welch RP, Wheeler E, Montasser ME, Luan J, Magi R, Strawbridge RJ, Rehnberg E, Gustafsson S, et al. Large-scale association analyses identify new loci influencing glycemic traits and provide insight into the underlying biological pathways. Nat Genet. 2012;44(9):991–1005.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  37. Cole TJ, Altman DG. Statistics notes: percentage differences, symmetry, and natural logarithms. BMJ. 2017;358:j3683.

    Article  PubMed  Google Scholar 

  38. Chang CC, Chow CC, Tellier LC, Vattikuti S, Purcell SM, Lee JJ. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience. 2015;4:7.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  39. White IR, Royston P, Wood AM. Multiple imputation using chained equations: issues and guidance for practice. Stat Med. 2011;30(4):377–99.

    Article  PubMed  Google Scholar 

  40. Boyd NF, Martin LJ, Yaffe MJ, Minkin S. Mammographic density and breast cancer risk: current understanding and future prospects. Breast Cancer Res. 2011;13(6):223.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Buschard K, Thomassen K, Lynge E, Vejborg I, Tjonneland A, von Euler-Chelpin M, Andersen ZJ. Diabetes, diabetes treatment, and mammographic density in Danish diet, Cancer, and Health cohort. Cancer Causes Control. 2017;28(1):13–21.

    Article  PubMed  Google Scholar 

  42. Jimenez-Garcia R, Hernandez-Barrera V, Carrasco-Garrido P, Gil A. Prevalence and predictors of breast and cervical cancer screening among Spanish women with diabetes. Diabetes Care. 2009;32(8):1470–2.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  43. Lipscombe LL, Hux JE, Booth GL. Reduced screening mammography among women with diabetes. Arch Intern Med. 2005;165(18):2090–5.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

The authors thank all participants of the Karma study, the study personnel for their devoted work during data collection, and Dr. Ralph Highnam and colleagues for their technical support with the Volpara software.

Funding

This work was financed by the Swedish Research Council (grant no. 2014-2271), the Swedish Cancer Society (grant no. CAN 2013/469), and the Cancer Society in Stockholm (grant no. 141092). The KARMA study is supported by the Märit and Hans Rausing Initiative Against Breast Cancer and the Cancer and Risk Prediction Center (CRisP), a Linneus Centre (grant no. 70867902) financed by the Swedish Research Council.

Genotyping of the OncoArray was principally funded from three sources: the PERSPECTIVE project, funded from the Government of Canada through Genome Canada and the Canadian Institutes of Health Research, the Ministère de l’Économie, de la Science et de l’Innovation du Québec through Genome Québec, and the Quebec Breast Cancer Foundation; the National Cancer Institute Genetic Associations and Mechanisms in Oncology (GAME-ON) initiative and Discovery, Biology and Risk of Inherited Variants in Breast Cancer (DRIVE) project (National Institutes of Health grants U19 CA148065 and X01HG007492); and Cancer Research UK (C1287/A10118 and C1287/A16563). The Breast Cancer Association Consortium (BCAC) is funded by Cancer Research UK (C1287/A16563), by the European Community’s Seventh Framework Programme under grant agreement 223175 (HEALTH-F2-2009-223175) (COGS), and by the European Union’s Horizon 2020 Research and Innovation Programme under grant agreements 633784 (B-CAST) and 634935 (BRIDGES). Genotyping of the iCOGS array was funded by the European Union (HEALTH-F2-2009-223175), Cancer Research UK (C1287/A10710), the Canadian Institutes of Health Research for the “CIHR Team in Familial Risks of Breast Cancer” program, and the Ministry of Economic Development, Innovation and Export Trade of Quebec (grant no. PSR-SIIRI-701. The funders had no role in the design of the study; the collection, analysis, or interpretation of the data; the writing of the manuscript; or the decision to submit the manuscript for publication.

Availability of data and materials

The datasets analyzed in the present study are available from the corresponding author upon reasonable request.

Author information

Authors and Affiliations

Authors

Contributions

SB and JSB conceived of and designed the study with input from AHR, NBP, and MD. JSB performed the statistical analyses. SB and JSB drafted the manuscript. All authors contributed to the interpretation of the data and critically reviewed the manuscript for publication. All authors read and approved the final version submitted for publication and are accountable for the accuracy and integrity of the work.

Corresponding author

Correspondence to Signe Borgquist.

Ethics declarations

Ethics approval and consent to participate

All participants signed informed consent forms, and the ethical review board at Karolinska Institutet approved the study (2010/958-31/1).

Consent for publication

All authors approved the manuscript and consented to its publication.

Competing interests

MD declares being employed by Novo Nordisk (Denmark). All other authors declare that they have no competing interests. Novo Nordisk had no role in the design of the study; the collection, analysis, or interpretation of the data; the writing of the manuscript; or the decision to submit the manuscript for publication.

Publisher’s Note

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

Additional files

Additional file 1:

Table S1. Prescriptions of insulin and insulin analogs dispensed in the year prior to study entry in T1D and T2D patients. Table S2. Prescriptions of glucose-lowering medication dispensed in the year prior to study entry in T2D patients not receiving insulin therapy. Table S3. Single-nucleotide polymorphisms included in the insulin genetic score. Table S4. Descriptive characteristics of insulin-treated T2D patients and non-insulin-treated T2D patients. Table S5. Geometric means and percentage differences of volumetric mammographic density comparing glargine insulin users to non-glargine insulin users (case only analyes). Table S6. Geometric means and percentage differences of volumetric mammographic density comparing insulin-treated T2D patients to non-insulin treated T2D patients. Table S7. Geometric means and percentage differences of volumetric mammographic density comparing insulin-treated T1D and T2D patients to age-matched individuals without diabetes by insulin treatment duration. Table S8. Geometric means and percentage differences of volumetric mammographic density by treatment duration in insulin-treated T1D patients (case only analyses). Table S9. Associations of fasting insulin single nucleotide polymorphisms with volumetric mammographic density. (DOCX 88 kb)

Additional file 2:

Supplementary methods: handing of missing covariate data. (DOCX 13 kb)

Additional file 3:

Figure S1. Summary of percentage differences in volumetric mammography density observed by insulin therapy across the different observational analyses. Association of insulin therapy with volumetric mammographic density in different observational analyses with age and BMI adjusted estimates (open circles) and multivariable adjusted estimates (closed circles). Betas represent % differences in percent dense, absolute dense and absolute non-dense volume for differences in insulin exposure. (PDF 476 kb)

Additional file 4:

Figure S2. Association of the insulin genetic score with potential confounders. Association of 18-SNP insulin genetic score with potential confounders in Karma sub-cohort of non-diabetic women with genotyping data. Associations were tested using linear regression and (multinomial) logistic regression, adjusting for age, six principal components and genotyping array. Betas represent differences in covariate level per 1-standard deviation increment in insulin genetic score. (PDF 225 kb)

Additional file 5:

Figure S3. Association of the insulin genetic score with volumetric mammographic density measures after additional adjustment for benign breast disease. Association of 18-SNP insulin genetic score with volumetric mammographic density in Karma sub-cohort of non-diabetic women with genotyping data, overall and stratified by genotyping array and with additional adjustment for benign breast disease. Associations with volumetric mammographic density were tested by linear regression, adjusting for age, body mass index, menopausal status, six principal components and benign breast disease. Analyses in total sub-cohort were additionally adjusted for genotyping array. All volumetric mammographic density measures were log-transformed prior to analyses, with betas representing % differences in volumetric mammographic density per 1-standard deviation increment in insulin genetic score. (PDF 146 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Borgquist, S., Rosendahl, A.H., Czene, K. et al. Long-term exposure to insulin and volumetric mammographic density: observational and genetic associations in the Karma study. Breast Cancer Res 20, 93 (2018). https://doi.org/10.1186/s13058-018-1026-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13058-018-1026-7

Keywords