• The risk of RT-induced breast cancer after HL is strongly associated with a PRS for breast cancer in the general population.

  • A PRS, based on 9 SNPs interacting with RT in the occurrence of breast cancer after HL, also increased RT-induced breast cancer risk.

Female Hodgkin lymphoma (HL) patients treated with chest radiotherapy (RT) have a very high risk of breast cancer. The contribution of genetic factors to this risk is unclear. We therefore examined 211 155 germline single-nucleotide polymorphisms (SNPs) for gene-radiation interaction on breast cancer risk in a case-only analysis including 327 breast cancer patients after chest RT for HL and 4671 first primary breast cancer patients. Nine SNPs showed statistically significant interaction with RT on breast cancer risk (false discovery rate, <20%), of which 1 SNP in the PVT1 oncogene attained the Bonferroni threshold for statistical significance. A polygenic risk score (PRS) composed of these SNPs (RT-interaction-PRS) and a previously published breast cancer PRS (BC-PRS) derived in the general population were evaluated in a case-control analysis comprising the 327 chest-irradiated HL patients with breast cancer and 491 chest-irradiated HL patients without breast cancer. Patients in the highest tertile of the RT-interaction-PRS had a 1.6-fold higher breast cancer risk than those in the lowest tertile. Remarkably, we observed a fourfold increased RT-induced breast cancer risk in the highest compared with the lowest decile of the BC-PRS. On a continuous scale, breast cancer risk increased 1.4-fold per standard deviation of the BC-PRS, similar to the effect size found in the general population. This study demonstrates that genetic factors influence breast cancer risk after chest RT for HL. Given the high absolute breast cancer risk in radiation-exposed women, these results can have important implications for the management of current HL survivors and future patients.

Women who are treated at young ages with chest radiotherapy (RT) for Hodgkin lymphoma (HL) have a 5 to 20 times increased risk of breast cancer compared with the general population.1-11  The cumulative incidence of breast cancer up to 40 years after treatment with mantle-field RT is 30% to 40%,5,6,10  in the range of risks observed in BRCA1/2 mutation carriers.12  The risk of RT-induced breast cancer rises with increasing radiation dose and volume, but not all female HL survivors treated with high-dose, high-volume RT develop breast cancer. Some variation in risk is explained by age at RT exposure, which is inversely related with breast cancer risk, and premature menopause induced by concomitant alkylating-chemotherapy treatment, which reduces risk.13  However, variation in risk may also be due to genetic factors. The high risk of breast cancer in this population provides an excellent opportunity to investigate the genetic basis for differential sensitivity to radiation carcinogenesis. Although it is well known that ionizing radiation induces DNA damage, the molecular mechanisms underlying radiation-induced breast carcinogenesis are unclear. To date, there is no clear evidence that known high-risk breast cancer susceptibility genes contribute to RT-induced breast cancer risk in HL patients.14-17  However, there may be a more important role for common susceptibility variants, as suggested by genetic association studies in women exposed to low-dose radiation, albeit with conflicting results.18-26  The role of single-nucleotide polymorphisms (SNPs) in breast cancer risk after therapeutic high-dose radiation has been investigated in few studies: a small genome-wide association study (GWAS) on any second solid malignancy in childhood HL survivors27  and a GWAS on radiation-induced breast cancer in childhood cancer survivors.28  In addition, Ma et al investigated 14 SNPs previously associated with breast cancer in the general population in HL survivors.29 

In the current study, we used a 2-step design to investigate whether there are subgroups of women exposed to chest RT who are genetically more susceptible to radiation-induced breast cancer. We first used a case-only analysis to evaluate interactions between 211 155 SNPs and chest RT, by comparing patients with breast cancer after chest RT for HL with first primary breast cancer patients previously unexposed to RT. We then conducted a nested breast cancer case-control analysis among chest-irradiated HL survivors to evaluate a polygenic risk score (PRS) composed of RT-interacting SNPs from the case-only analysis (RT-interaction-PRS). As a separate aim, we studied the effect of a previously published breast cancer PRS (BC-PRS)30  in the general population on breast cancer risk among chest-irradiated HL survivors.

Study design

When studying interaction between RT and genetic variation on breast cancer after HL, a classical case-control study nested in a cohort of HL survivors would not be informative because, until recently, 90% of the patients with breast cancer after HL received RT, resulting in too few unexposed cases. Therefore, we used a 2-step design to identify susceptibility variants for radiation-induced breast cancer (Figure 1). First, we examined gene-radiation interaction for 211 155 SNPs in a case-only analysis comparing patients with breast cancer after chest RT for HL (further referred to as breast cancer after HL cases) and first primary breast cancer patients (further referred to as first primary breast cancer cases). For each SNP, we used logistic regression analysis to estimate the per-allele interaction odds ratio (IOR), a measure of departure from a multiplicative joint effect of the SNP and chest RT, for the risk of breast cancer, assuming independence between chest RT and the SNP in women from the general population.31 

Figure 1.

Study design. (A) First step: case-only study. (B) Second step: case-control study. IOR, interaction odds ratio. *See Mavaddat et al.30 

Figure 1.

Study design. (A) First step: case-only study. (B) Second step: case-control study. IOR, interaction odds ratio. *See Mavaddat et al.30 

Close modal

Second, we combined interacting SNPs in a PRS, that is, the sum of risk alleles weighted by their effect size (see supplemental Methods A [available on the Blood Web site] for details) and evaluated the association between this PRS and the risk of breast cancer after chest RT in a breast cancer case-control analysis among irradiated HL survivors, using an independent control group of chest-irradiated HL survivors without breast cancer as controls (further referred to as HL controls). We similarly evaluated a second PRS, which was previously reported to be associated with breast cancer in the general population (the BC-PRS).30 

Study population and genotyping

For the case-only analysis, we pooled 339 cases with breast cancer after HL from 3 breast cancer case-control studies29,32-34  nested in HL survivor cohorts: the Childhood Cancer Survivor study (CCSS),35  a British HL cohort,10  and the Dutch Hodgkin Lymphoma Cohort.6  Blood samples from these cases were genotyped using a custom Illumina iSelect Array comprising 211 155 SNPs, specially designed for the European Collaborative Oncological Gene-Environment Study (EU-COGS) project (referred to as iCOGS array).36  Extensive patient and HL treatment characteristics, as well as follow-up data, were available from medical records,4,29,35  through questionnaires sent to general practitioners and study participants, and from record linkages with national cancer registries.6,10,14,29,32-34  Female patients with breast cancer after HL were included in our study if they were diagnosed with primary breast cancer >8 years after chest RT for HL before the age of 41 years (see supplemental Methods B for definition of chest RT). Cases with breast cancer after HL were frequency matched (1:∼14) on age and year of breast cancer diagnosis (5-year intervals) and country, to 4673 first primary breast cancer cases of European origin not known to be exposed to chest RT. These were selected from 19 275 participants of 10 studies from The Netherlands, United Kingdom, and the United States within the Breast Cancer Association Consortium (BCAC)36  for whom iCOGS genotype data were available. When there were too few subjects in a specific age category, we oversampled in an adjacent age category in the same calendar year category.

For the case-control analysis, we included the 339 cases with breast cancer after HL mentioned previously and 508 HL survivors treated with chest RT who did not develop breast cancer until end of follow-up, available from the 3 breast cancer case-control studies described earlier in this section. For all HL controls without breast cancer, we collected similar data as described earlier in this section for cases with breast cancer after HL. In the published original case-control studies,13,32-34  which examined radiation dose-response, 1 to 4 controls were individually matched to each case. Controls had to have survived without breast cancer at least as long as the interval between HL and breast cancer for the corresponding case, and in case of the US study, had to have donated a blood sample. In addition, controls had to match the case on age at HL treatment (±3 years) and date of HL treatment (±5 years). Controls from the original case-control studies were excluded if they were not treated with chest RT, were treated at or after age 41 years, and/or did not donate a blood sample. In addition, controls were excluded if they developed breast cancer after the year of breast cancer diagnosis of the case to whom they had been matched. For the current study, we added recently diagnosed breast cancer after HL cases fulfilling the inclusion criteria. All separate studies involved in this collaboration were approved by the relevant institutional review boards, and all individuals gave written informed consent.

Quality control on genotype data

After quality control, 194 106 SNPs measured in 4671 first primary breast cancer cases, 327 cases with breast cancer after HL, and 491 HL controls without breast cancer remained for analyses. See supplemental Methods C for details on quality control.

Statistical analyses

In the case-only analysis, comparing breast cancer after HL cases and first primary breast cancer cases, we estimated the per-allele IOR by unconditional logistic regression analysis for all variants passing quality control, adjusting for the matching factors (age and year of breast cancer diagnosis, both continuous, and country) and the first principal component describing remaining genetic ethnic differences among European subjects (referred to as ethnicity) (see supplemental Methods D). P values for the IORs were calculated by the score test performed using the GenABEL package within R (see supplemental Methods E). Based on a conservative Bonferroni correction, SNPs with a P value < 2.6E-07 were considered statistically significant. Furthermore, we applied the false discovery rate (FDR) by Benjamini and Hochberg37  to identify SNPs among which the expected proportion of false positives is <20% (q value = 0.2). For significant SNPs in linkage disequilibrium (LD; r2 > 0.7), only the SNP with the lowest P value was included in the PRS.

Subsequently, for all subjects in the case-control analysis (breast cancer after HL cases and HL controls), we calculated the RT-interaction-PRS consisting of SNPs interacting with RT on breast cancer at 20% FDR and the 77 SNP BC-PRS. Missing genotypes were imputed by the mode among HL controls. Odds ratios (ORs) and corresponding 95% confidence intervals (CIs) for RT-induced breast cancer after HL were calculated by unconditional logistic regression per standard deviation (SD) increase in either the RT-interaction-PRS and/or the BC-PRS, adjusted for each other and for age at HL diagnosis (continuous), year of HL diagnosis (4 periods), country, and ethnicity. We also calculated ORs for breast cancer by categories of the PRSs (tertiles for the RT-interaction-PRS and deciles for the BC-PRS). P values for the ORs were based on Wald tests. Interaction between the RT-interaction-PRS and gonadotoxic treatment of HL (yes/no) and between the RT-interaction-PRS and age at HL treatment (≤20/>20 years) was tested by stratification on these factors. As a sensitivity analysis, we assessed the association of a PRS including only the SNPs attaining the Bonferroni threshold for statistical significance. For all Bonferroni-significant SNPs in the case-only analysis, we also tested their individual association with breast cancer after chest RT in the case-control analysis using logistic regression adjusted for age at and year of HL diagnosis, country, and ethnicity. All analyses were conducted with R software (www.r-project.org).

Data availability

Nonidentifiable data that support the findings of this study will be made available upon reasonable request. Access to the BCAC data are governed by the Data Access Coordinating Committee from BCAC. Data from the CCSS study can be retrieved from dbGAP using accession number phs001327.v1.p1.

Study populations of the case-only and case-control analysis

We included 327 breast cancer after HL cases from cohorts of female HL patients in The Netherlands, United Kingdom, and United States and 4671 frequency-matched first primary breast cancer cases previously unexposed to RT from the same countries in the case-only analysis. Furthermore, we included 491 HL controls in the case-control analysis (see Table 1 for the numbers of subjects by country). The median age at breast cancer diagnosis was 45 years (range, 24-76 years) for breast cancer after HL cases and 46 years (range, 22-84 years) for age-matched first primary breast cancer cases (Table 1). The median interval between HL and breast cancer diagnosis was 24 years (range, 9-46 years). For HL controls, median follow-up was 30 years (range, 9-49 years). Most HL cases and controls (87%) were treated with mantle-field irradiation, whereas 11% of the HL cases and controls received mediastinal radiotherapy without axillary node radiotherapy. Approximately one-half of the breast cancer after HL cases and almost 60% of HL controls were treated for HL with chemotherapy in addition to RT. Approximately 45% of breast cancer after HL cases and 57% of HL controls received gonadotoxic treatment (ie, alkylating chemotherapy and/or pelvic RT).

SNPs interacting with RT on breast cancer risk (case-only analysis)

We tested 194 106 SNPs that passed quality control for an interaction with chest RT in the case-only analysis of breast cancer patients (QQ plot is depicted in supplemental Figure 1). As shown in Table 2, 3 SNPs were statistically significantly associated at the Bonferroni threshold for multiple testing (P < 2.6E-07) and 7 additional SNPs met the 20% FDR threshold, of which 1 SNP was excluded because of strong LD (r2 = 0.9). The estimated per-allele IORs for these 9 SNPs ranged from 1.6 to 2.2. Most SNPs were quite common in the breast cancer after HL cases with minor allele frequencies (MAFs) between 2.8% and 43.7%.

PRS for RT-induced breast cancer (case-control analysis)

We constructed an RT-interaction-PRS of the 9 SNPs that showed a statistically significant (FDR, 20%) interaction with RT-induced breast cancer. The RT-interaction-PRS increased breast cancer risk after chest RT for HL with ORs of 1.2 (95% CI, 0.8-1.7; P = .348) and 1.6 (95% CI, 1.1-2.4; P = .007), respectively, for the middle and highest tertiles compared with the lowest tertile, adjusted for age and year of HL diagnosis, country, ethnicity, and the BC-PRS (Figure 2; supplemental Table 1). The OR per 1 SD of the RT-interaction-PRS was 1.3 (95% CI, 1.1-1.5; P = .002). Additional adjustment for gonadotoxic treatment did not affect the association of the RT-interaction-PRS with breast cancer risk (ORadjusted, 1.3; 95% CI, 1.1-1.5), suggesting that it is unlikely that chemotherapy has confounded our analyses. In addition, stratified analyses resulted in similar associations between the RT-interaction-PRS and breast cancer risk among women who received gonadotoxic treatment (alkylating chemotherapy and/or pelvic RT) and women who did not; we observed no statistically significant interaction between the RT-interaction-PRS and gonadotoxic treatment (P = .337; supplemental Table 1). Likewise, stratification by age at HL treatment (≤20, >20 years) did not result in different associations between the RT-interaction-PRS and breast cancer risk after RT for HL; there was no interaction between age at HL treatment and the RT-interaction PRS (P = .954).

Figure 2.

Risk of breast cancer after chest RT by tertiles of the RT-interaction-PRS among HL survivors.

Figure 2.

Risk of breast cancer after chest RT by tertiles of the RT-interaction-PRS among HL survivors.

Close modal

In a sensitivity analysis, we observed that a PRS containing only the 3 SNPs reaching the Bonferroni threshold for statistical significance also increased breast cancer risk with ORs of 1.4 (95% CI, 1.0-2.1; P = .070) and 1.6 (95% CI, 1.1-2.2; P = .018), respectively, for the middle and highest tertile compared with the lowest tertile which consisted of noncarriers. The OR per 1 SD of the 3 SNP RT-interaction-PRS was 1.2 (95% CI, 1.0-1.4; P = .014).

To confirm the observed associations, we also evaluated the individual effects of the 3 Bonferroni-significant SNPs on RT-induced breast cancer in the case-control analysis among chest-irradiated HL survivors (supplemental Table 2). Of these, an intronic variant in oncogene PVT1 (rs10505506) was associated with RT-induced breast cancer risk after HL with an OR of 1.3 (95% CI, 1.1-1.6; P = .007) per allele copy. Of note, rs10505506 is not in LD (r2 < 0.3 in Europeans from the 1000 Genomes Project38 ) with previously identified cancer risk variants in the PVT1 locus (supplemental Figure 2).

PRS based on known breast cancer SNPs (case-control analysis)

To evaluate the combined effect of known breast cancer SNPs, we studied a BC-PRS containing 76 SNPs that increase breast cancer risk in the general population30  in chest-irradiated HL survivors. The BC-PRS was associated with a 1.4-fold increased risk of RT-induced breast cancer (95% CI, 1.2-1.6; P = 9.1E-05) per SD increase in the BC-PRS. The ORs for developing breast cancer after chest RT for HL by deciles of the BC-PRS, compared with women in the middle quintile (40th to 60th percentile), are shown in Figure 3 and supplemental Table 3. The 10% of women with the lowest BC-PRS had an OR of 0.6 (95% CI, 0.3-1.1; P = .133) for developing RT-induced breast cancer compared with women in the middle quintile, whereas the OR for the 10% of women with the highest BC-PRS was 2.4 (95% CI, 1.4-4.2; P = .002), adjusted for age and year of HL diagnosis, country, ethnicity, and the RT-interaction-PRS (in tertiles). This results in a fourfold relative risk for the 10% of women with the highest compared with the lowest BC-PRS. There was no interaction between the RT-interaction-PRS and the BC-PRS (P = .645).

Figure 3.

Risk of breast cancer after chest RT by deciles of the BC-PRS in the breast cancer after HL case-control analysis.

Figure 3.

Risk of breast cancer after chest RT by deciles of the BC-PRS in the breast cancer after HL case-control analysis.

Close modal

This study demonstrates that genetic factors influence the risk of breast cancer after chest RT for HL. We showed that a BC-PRS, consisting of 77 SNPs previously associated with breast cancer in the general population, also substantially increases the risk of breast cancer in chest-irradiated HL survivors. In addition, we identified 9 SNPs interacting with chest RT and the risk of breast cancer after HL and we showed a statistically significant association of a PRS composed of these interaction SNPs with breast cancer risk after chest RT for HL using an independent control group. These results imply that the absolute risk of breast cancer due to irradiation would be (even) larger among women at high genetic risk, which is relevant for clinical risk prediction.

Importantly, we validated the previously published BC-PRS in a high-risk population of female chest-irradiated HL survivors and found that there are large differences in risk between women with a low and high PRS. More specifically, we observed a fourfold increased relative risk between chest-irradiated HL survivors in the highest compared with the lowest decile of the BC-PRS. On a continuous scale, the effect size was very similar to that found in the general population (OR of 1.4 per SD in our study of HL survivors compared with ORs of 1.4 to 1.6 per SD in the general population).30,39  These results indicate that the effects of radiation exposure and common susceptibility variants, summarized in the PRS, combine approximately multiplicatively. Given the high absolute breast cancer risk in radiation-exposed women, these results have important implications for their management. The BC-PRS can be used to help guide treatment decisions in newly diagnosed HL patients as well as to help determine breast surveillance strategies for irradiated HL survivors. Annual breast cancer surveillance between the ages of 25 and 50 years is currently recommended by the International Late Effects of Childhood Cancer Guideline Harmonization Group for female survivors of childhood, adolescent, and young adult cancer who received ≥20 Gy chest radiation before age 30 years.40  Less clear is the evidence for surveillance in women treated at older ages, with lower dosages, or with different radiation volumes. Therefore, clinical prediction models for breast cancer that include both clinical and genetic factors can help to identify (additional) women who may benefit from breast cancer surveillance.

We chose to evaluate the 77 SNP BC-PRS by Mavaddat et al30  as this PRS has been associated with breast cancer risk in the general population and in high-risk groups such as BRCA1 and BRCA2 mutation carriers,39,41  allowing direct comparison of the reported effect sizes. Nevertheless, many more common susceptibility variants have recently been identified for breast cancer in the general population.42,43  Addition of these SNPs to the BC-PRS may further improve risk stratification for breast cancer in chest-irradiated HL survivors and in other high-risk groups. Inclusion of SNPs associated with hormone receptor-negative breast cancer may be of particular interest, as several studies have reported that HL survivors are more likely to develop hormone receptor-negative disease.44-46 

We applied an innovative design to examine the role of SNP-radiation interactions in breast cancer risk after HL. This is not feasible in a classical breast cancer case-control study in HL survivors, as, until recently, ∼90% of breast cancer cases after HL received chest RT. Therefore, we first performed a case-only analysis in breast cancer patients previously exposed and unexposed to chest RT, followed by a case-control analysis in HL survivors to evaluate the combined effect of the identified RT-interaction SNPs in a PRS. We used a 20% FDR as a cutoff to select SNPs interacting with RT for the RT-interaction-PRS, as it has been shown that the performance of a PRS improves when using more liberal thresholds than the conservative Bonferroni threshold.47,48  Although a PRS consisting of 3 SNPs statistically significant at the Bonferroni threshold showed a similar association with RT-induced breast cancer risk among HL survivors, the goodness of fit was better in the full PRS (data not shown).

The IORs that we estimated in the case-only analysis measure departure from a multiplicative joint effect of chest RT and the SNP, assuming independence between chest RT and the SNP in women from the general population.49  This assumption is likely to be justified except for SNPs associated with HL. SNPs associated with HL may also have shown a significant IOR in the case-only analysis. On the other hand, such SNPs may be associated with both HL and (radiation-induced) breast cancer and, therefore, we did not exclude SNPs previously associated with HL from inclusion in the RT-interaction-PRS. If they were only associated with HL they would have attenuated the association of the RT-interaction-PRS with breast cancer after chest RT in the case-control analysis. In the case-only analysis, we identified 1 SNP (rs9461776) interacting with radiation at 20% FDR significance located in the HLA region, which has extensively been reported to be associated with HL.50  rs9461776 showed no evidence of an association with breast cancer after chest RT (OR, 1.0; 95% CI, 0.8-1.4; P > .5) in the case-control analysis in HL survivors and may therefore have attenuated the association of the RT-interaction-PRS with the risk of breast cancer after chest RT.

Of the 9 SNPs (MAF > 1%) interacting with RT on breast cancer risk at 20% FDR, 1 attained the genome-wide level (P < 5 × 10−8) of statistical significance. This SNP (rs10505506) was also associated with breast cancer risk in chest-irradiated HL survivors (OR, 1.3; 95% CI, 1.1-1.6; P = .007). SNP rs10505506 is located in the intronic region of PVT1, which is a known oncogene regulated by tumor-suppressor p53 encoding a long noncoding RNA and several microRNAs.51,52 PVT1 has been shown to interact with the adjacent proto-oncogene MYC and translocations in this locus have been associated with Burkitt lymphoma. In addition, overexpression of PVT1 is associated with several types of cancers including breast cancer, acute myeloid leukemia, and HL. Likewise, GWAS studies have identified several conditionally independent SNPs in this locus associated with cancer, including breast cancer and HL,53,54  but none of these are in LD with rs10505506. A potential link with radiation has recently been suggested in a mouse model after whole-body irradiation.55 

The association of the RT-interaction-PRS with breast cancer risk after HL was not weakened in “low-risk” groups of women irradiated at older age (ie, 20 years or older) or women treated with gonadotoxic treatment. In addition, we did not observe interaction between the RT-interaction-PRS and either gonadotoxic treatment or age at HL treatment. This is in line with the notion that gonadotoxic treatment and age at HL treatment are independent risk factors for breast cancer risk after HL. This suggests that age and treatment-related risk factors for breast cancer after HL and the genetic risk scores (both the RT-interaction-PRS and the BC-PRS) combine multiplicatively as has previously been shown for several reproductive risk factors and the 77 SNP BC-PRS in the general population.56 

A limitation of this study is that the study populations for the construction and evaluation of the RT-interaction-PRS were not independent, as the breast cancer after HL cases were included in both analyses. External validation of the RT-interaction-PRS in an independent study is therefore needed to confirm our findings. In addition, we excluded SNPs with a low MAF (<1%) from our analyses, as these low-frequency SNPs are more prone to genotyping errors. However, Morton et al recently reported 2 suggestive associations for low-frequency variants at 11q23 and 1q32.3, both not present on the iCOGS array, with breast cancer risk after childhood cancer,28  suggesting a potential role for low-frequency SNPs in RT-induced breast cancer. Inclusion of these SNPs to the RT-interaction-PRS might strengthen its association with RT-induced breast cancer. Likewise, additional SNPs interacting with RT on breast cancer may be identified when assessing SNP data from denser genotyping chips imputed to a reference panel. However, in this first analysis, we focused on high-quality SNPs specifically selected for the iCOGS array.

In conclusion, we showed that a BC-PRS previously developed in the general population also applies in a high-risk breast cancer population of chest-irradiated HL survivors. In addition, we developed an RT-interaction-PRS composed of 9 SNPs interacting with radiation that was associated with raised breast cancer risk after chest RT for HL. Although our RT-interaction-PRS needs validation in an independent sample, the BC-PRS can already be applied in clinical practice. This can benefit treatment decision-making in future HL patients as well as identification of high-risk survivors eligible for breast cancer surveillance.

The online version of this article contains a data supplement.

The publication costs of this article were defrayed in part by page charge payment. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734.

For the Dutch Hodgkin Lymphoma cohort, the authors thank Ph. M. P. Poortmans from Elisabeth Hospital Tilburg at the time of the study, M. L. M. Lybeert from Catherina Hospital Eindhoven, G. W. Imhoff from Groningen University Medical Center, J. M. M. Raemaekers from Rijnstate Hospital/Radboud University Medical Center Nijmegen at the time of the study, J. M. Zijlstra from VU Medical Center Amsterdam, and J. H. Borger from Maastricht University Medical Center for inclusion of HL patients. The UK Hodgkin Lymphoma cohort thanks the study participants, study staff, and clinicians who participated in the study as listed in the authorship and acknowledgments in Cooke et al.32 

This work was supported by Dutch Cancer Society grant 2010-4720 (F.E.v.L.). The UK Hodgkin Lymphoma cohort was supported by Breast Cancer Now and the European Commission. The Institute of Cancer Research acknowledges funding from the National Institute for Health Research to the Biomedical Research Centre. The Childhood Cancer Survivor study was supported by the National Institutes of Health, National Cancer Institute (CA55727 [principal investigator: G. T. Armstrong]). Support to St. Jude Children’s Research Hospital was also provided by the Cancer Center Support (CORE) grant (CA21765 [principal investigator: C. Roberts]) and the American Lebanese-Syrian Associated Charities (ALSAC).

All sponsors of this work are public or nonprofit organizations that support science in general. They had no role in gathering, analyzing, or interpreting the data.

Contribution: A.W.J.O.-v.W., M.H., M.K.S., H.G.d.H., and F.E.v.L. performed the analyses, interpreted the data, and drafted the manuscript; A.W.J.O.-v.W., H.G.d.H., F.H.v.d.B., M.H., and M.K.S. performed statistical analysis; F.E.v.L. N.S.R., C.P.M.J., A.D.G.K., J.D., H.A.-C., C.A.H, E.J.S., A.C., P.D., M.J.H., J.P., F.J.C., P.P., N.O., D.F.E., B.M.P.A., L.C.S., S.B., R.C., L.L.R., and A.J.S. contributed in the inclusion of patients and data collection; H.G.d.H., M.L.D.B., and A.M.v.E. provided administrative support and coordinated the work; M.H., M.K.S., A.B., and F.E.v.L. supervised the study; and all authors critically revised the manuscript.

Conflict-of-interest disclosure: The authors declare no competing financial interests.

Correspondence: Flora E. van Leeuwen, Department of Epidemiology, The Netherlands Cancer Institute, Plesmanlaan 121, 1066 CX Amsterdam, The Netherlands; e-mail: f.v.leeuwen@nki.nl.

1.
Aisenberg
AC
,
Finkelstein
DM
,
Doppke
KP
,
Koerner
FC
,
Boivin
JF
,
Willett
CG
.
High risk of breast carcinoma after irradiation of young women with Hodgkin’s disease
.
Cancer
.
1997
;
79
(
6
):
1203
-
1210
.
2.
Bhatia
S
,
Yasui
Y
,
Robison
LL
, et al
;
Late Effects Study Group
.
High risk of subsequent neoplasms continues with extended follow-up of childhood Hodgkin’s disease: report from the Late Effects Study Group
.
J Clin Oncol
.
2003
;
21
(
23
):
4386
-
4394
.
3.
De Bruin
ML
,
Sparidans
J
,
van’t Veer
MB
, et al
.
Breast cancer risk in female survivors of Hodgkin’s lymphoma: lower risk after smaller radiation volumes
.
J Clin Oncol
.
2009
;
27
(
26
):
4239
-
4246
.
4.
Hodgson
DC
,
Gilbert
ES
,
Dores
GM
, et al
.
Long-term solid cancer risk among 5-year survivors of Hodgkin’s lymphoma
.
J Clin Oncol
.
2007
;
25
(
12
):
1489
-
1497
.
5.
Moskowitz
CS
,
Chou
JF
,
Wolden
SL
, et al
.
Breast cancer after chest radiation therapy for childhood cancer
.
J Clin Oncol
.
2014
;
32
(
21
):
2217
-
2223
.
6.
Schaapveld
M
,
Aleman
BM
,
van Eggermond
AM
, et al
.
Second cancer risk up to 40 years after treatment for Hodgkin’s lymphoma
.
N Engl J Med
.
2015
;
373
(
26
):
2499
-
2511
.
7.
Travis
LB
,
Hill
D
,
Dores
GM
, et al
.
Cumulative absolute breast cancer risk for young women treated for Hodgkin lymphoma
.
J Natl Cancer Inst
.
2005
;
97
(
19
):
1428
-
1437
.
8.
van Leeuwen
FE
,
Klokman
WJ
,
Veer
MB
, et al
.
Long-term risk of second malignancy in survivors of Hodgkin’s disease treated during adolescence or young adulthood
.
J Clin Oncol
.
2000
;
18
(
3
):
487
-
497
.
9.
Hancock
SL
,
Tucker
MA
,
Hoppe
RT
.
Breast cancer after treatment of Hodgkin’s disease
.
J Natl Cancer Inst
.
1993
;
85
(
1
):
25
-
31
.
10.
Swerdlow
AJ
,
Cooke
R
,
Bates
A
, et al
.
Breast cancer risk after supradiaphragmatic radiotherapy for Hodgkin’s lymphoma in England and Wales: a national cohort study
.
J Clin Oncol
.
2012
;
30
(
22
):
2745
-
2752
.
11.
Ng
AK
,
Bernardo
MV
,
Weller
E
, et al
.
Second malignancy after Hodgkin disease treated with radiation therapy with or without chemotherapy: long-term risks and risk factors
.
Blood
.
2002
;
100
(
6
):
1989
-
1996
.
12.
Kuchenbaecker
KB
,
Hopper
JL
,
Barnes
DR
, et al
;
BRCA1 and BRCA2 Cohort Consortium
.
Risks of breast, ovarian, and contralateral breast cancer for BRCA1 and BRCA2 mutation carriers
.
JAMA
.
2017
;
317
(
23
):
2402
-
2416
.
13.
Krul
IM
,
Opstal-van Winden
AWJ
,
Aleman
BMP
, et al
.
Breast cancer risk after radiation therapy for hodgkin lymphoma: influence of gonadal hormone exposure
.
Int J Radiat Oncol Biol Phys
.
2017
;
99
(
4
):
843
-
853
.
14.
Broeks
A
,
Russell
NS
,
Floore
AN
, et al
.
Increased risk of breast cancer following irradiation for Hodgkin’s disease is not a result of ATM germline mutations
.
Int J Radiat Biol
.
2000
;
76
(
5
):
693
-
698
.
15.
Gaffney
DK
,
Hemmersmeier
J
,
Holden
J
, et al
.
Breast cancer after mantle irradiation for Hodgkin’s disease: correlation of clinical, pathologic, and molecular features including loss of heterozygosity at BRCA1 and BRCA2
.
Int J Radiat Oncol Biol Phys
.
2001
;
49
(
2
):
539
-
546
.
16.
Nichols
KE
,
Heath
JA
,
Friedman
D
, et al
.
TP53, BRCA1, and BRCA2 tumor suppressor genes are not commonly mutated in survivors of Hodgkin’s disease with second primary neoplasms
.
J Clin Oncol
.
2003
;
21
(
24
):
4505
-
4509
.
17.
Nichols
KE
,
Levitz
S
,
Shannon
KE
, et al
.
Heterozygous germline ATM mutations do not contribute to radiation-associated malignancies after Hodgkin’s disease
.
J Clin Oncol
.
1999
;
17
(
4
):
1259
.
18.
Bhatti
P
,
Doody
MM
,
Alexander
BH
, et al
.
Breast cancer risk polymorphisms and interaction with ionizing radiation among U.S. radiologic technologists
.
Cancer Epidemiol Biomarkers Prev
.
2008
;
17
(
8
):
2007
-
2011
.
19.
Bhatti
P
,
Doody
MM
,
Rajaraman
P
, et al
.
Novel breast cancer risk alleles and interaction with ionizing radiation among U.S. radiologic technologists
.
Radiat Res
.
2010
;
173
(
2
):
214
-
224
.
20.
Bhatti
P
,
Struewing
JP
,
Alexander
BH
, et al
.
Polymorphisms in DNA repair genes, ionizing radiation exposure and risk of breast cancer in U.S. Radiologic technologists
.
Int J Cancer
.
2008
;
122
(
1
):
177
-
182
.
21.
Rajaraman
P
,
Bhatti
P
,
Doody
MM
, et al
.
Nucleotide excision repair polymorphisms may modify ionizing radiation-related breast cancer risk in US radiologic technologists
.
Int J Cancer
.
2008
;
123
(
11
):
2713
-
2716
.
22.
Schonfeld
SJ
,
Bhatti
P
,
Brown
EE
, et al
.
Polymorphisms in oxidative stress and inflammation pathway genes, low-dose ionizing radiation, and the risk of breast cancer among US radiologic technologists
.
Cancer Causes Control
.
2010
;
21
(
11
):
1857
-
1866
.
23.
Sigurdson
AJ
,
Bhatti
P
,
Chang
SC
, et al
.
Polymorphisms in estrogen biosynthesis and metabolism-related genes, ionizing radiation exposure, and risk of breast cancer among US radiologic technologists
.
Breast Cancer Res Treat
.
2009
;
118
(
1
):
177
-
184
.
24.
Sigurdson
AJ
,
Bhatti
P
,
Doody
MM
, et al
.
Polymorphisms in apoptosis- and proliferation-related genes, ionizing radiation exposure, and risk of breast cancer among U.S. Radiologic Technologists
.
Cancer Epidemiol Biomarkers Prev
.
2007
;
16
(
10
):
2000
-
2007
.
25.
van Heemst
D
,
den Reijer
PM
,
Westendorp
RG
.
Ageing or cancer: a review on the role of caretakers and gatekeepers
.
Eur J Cancer
.
2007
;
43
(
15
):
2144
-
2152
.
26.
Varma
G
,
Varma
R
,
Huang
H
, et al
.
Array comparative genomic hybridisation (aCGH) analysis of premenopausal breast cancers from a nuclear fallout area and matched cases from western New York
.
Br J Cancer
.
2005
;
93
(
6
):
699
-
708
.
27.
Best
T
,
Li
D
,
Skol
AD
, et al
.
Variants at 6q21 implicate PRDM1 in the etiology of therapy-induced second malignancies after Hodgkin’s lymphoma
.
Nat Med
.
2011
;
17
(
8
):
941
-
943
.
28.
Morton
LM
,
Sampson
JN
,
Armstrong
GT
, et al
.
Genome-wide association study to identify susceptibility loci that modify radiation-related risk for breast cancer after childhood cancer
.
J Natl Cancer Inst
.
2017
;
109
(
11
).
29.
Ma
YP
,
van Leeuwen
FE
,
Cooke
R
, et al
.
FGFR2 genotype and risk of radiation-associated breast cancer in Hodgkin lymphoma
.
Blood
.
2012
;
119
(
4
):
1029
-
1031
.
30.
Mavaddat
N
,
Pharoah
PD
,
Michailidou
K
, et al
.
Prediction of breast cancer risk based on profiling with common genetic variants
.
J Natl Cancer Inst
.
2015
;
107
(
5
).
31.
Piegorsch
WW
,
Weinberg
CR
,
Taylor
JA
.
Non-hierarchical logistic models and case-only designs for assessing susceptibility in population-based case-control studies
.
Stat Med
.
1994
;
13
(
2
):
153
-
162
.
32.
Cooke
R
,
Jones
ME
,
Cunningham
D
, et al
;
England and Wales Hodgkin Lymphoma Follow-up Group
.
Breast cancer risk following Hodgkin lymphoma radiotherapy in relation to menstrual and reproductive factors
.
Br J Cancer
.
2013
;
108
(
11
):
2399
-
2406
.
33.
Inskip
PD
,
Robison
LL
,
Stovall
M
, et al
.
Radiation dose and breast cancer risk in the childhood cancer survivor study
.
J Clin Oncol
.
2009
;
27
(
24
):
3901
-
3907
.
34.
van Leeuwen
FE
,
Klokman
WJ
,
Stovall
M
, et al
.
Roles of radiation dose, chemotherapy, and hormonal factors in breast cancer following Hodgkin’s disease
.
J Natl Cancer Inst
.
2003
;
95
(
13
):
971
-
980
.
35.
Kenney
LB
,
Yasui
Y
,
Inskip
PD
, et al
.
Breast cancer after childhood cancer: a report from the Childhood Cancer Survivor Study
.
Ann Intern Med
.
2004
;
141
(
8
):
590
-
597
.
36.
Michailidou
K
,
Hall
P
,
Gonzalez-Neira
A
, et al
.
Large-scale genotyping identifies 41 new loci associated with breast cancer risk
.
Nat Genet
.
2013
;45(4):353-361.
37.
Benjamini
Y
,
Hochberg
Y
.
Controlling the false discovery rate: a practical and powerful approach to multiple testing
.
J R Stat Soc B
.
1995
;
57
(
1
):
289
-
300
.
38.
Auton
A
,
Brooks
LD
,
Durbin
RM
, et al
;
1000 Genomes Project Consortium
.
A global reference for human genetic variation
.
Nature
.
2015
;
526
(
7571
):
68
-
74
.
39.
Muranen
TA
,
Mavaddat
N
,
Khan
S
, et al
.
Polygenic risk score is associated with increased disease risk in 52 Finnish breast cancer families
.
Breast Cancer Res Treat
.
2016
;
158
(
3
):
463
-
469
.
40.
Mulder
RL
,
Kremer
LC
,
Hudson
MM
, et al
;
International Late Effects of Childhood Cancer Guideline Harmonization Group
.
Recommendations for breast cancer surveillance for female survivors of childhood, adolescent, and young adult cancer given chest radiation: a report from the International Late Effects of Childhood Cancer Guideline Harmonization Group
.
Lancet Oncol
.
2013
;
14
(
13
):
e621
-
e629
.
41.
Kuchenbaecker
KB
,
McGuffog
L
,
Barrowdale
D
, et al
.
Evaluation of polygenic risk scores for breast and ovarian cancer risk prediction in BRCA1 and BRCA2 mutation carriers
.
J Natl Cancer Inst
.
2017
;
109
(
7
).
42.
Michailidou
K
,
Lindström
S
,
Dennis
J
, et al
;
ConFab/AOCS Investigators
.
Association analysis identifies 65 new breast cancer risk loci
.
Nature
.
2017
;
551
(
7678
):
92
-
94
.
43.
Milne
RL
,
Kuchenbaecker
KB
,
Michailidou
K
, et al
;
NBSC Collaborators
.
Identification of ten variants associated with risk of estrogen-receptor-negative breast cancer
.
Nat Genet
.
2017
;
49
(
12
):
1767
-
1778
.
44.
Dores
GM
,
Anderson
WF
,
Beane Freeman
LE
,
Fraumeni
JF
Jr
,
Curtis
RE
.
Risk of breast cancer according to clinicopathologic features among long-term survivors of Hodgkin’s lymphoma treated with radiotherapy
.
Br J Cancer
.
2010
;
103
(
7
):
1081
-
1084
.
45.
Horst
KC
,
Hancock
SL
,
Ognibene
G
, et al
.
Histologic subtypes of breast cancer following radiotherapy for Hodgkin lymphoma
.
Ann Oncol
.
2014
;
25
(
4
):
848
-
851
.
46.
Milano
MT
,
Li
H
,
Gail
MH
,
Constine
LS
,
Travis
LB
.
Long-term survival among patients with Hodgkin’s lymphoma who developed breast cancer: a population-based study
.
J Clin Oncol
.
2010
;
28
(
34
):
5088
-
5096
.
47.
Dudbridge
F
.
Power and predictive accuracy of polygenic risk scores [published correction appears in PLoS Genet. 2013;9(4)]
.
PLoS Genet
.
2013
;
9
(
3
):
e1003348
.
48.
Purcell
SM
,
Wray
NR
,
Stone
JL
, et al
;
International Schizophrenia Consortium
.
Common polygenic variation contributes to risk of schizophrenia and bipolar disorder
.
Nature
.
2009
;
460
(
7256
):
748
-
752
.
49.
Khoury
MJ
,
Millikan
R
,
Gwinn
M
.
Genetic and molecular epidemiology
. In:
Rothman
KJ
,
Greenland
S
,
Lash
TL
, eds.
Modern Epidemiology
. 3rd ed.
Philadelphia, PA
:
Lippincott Williams & Wilkins
;
2008
:
573
.
50.
Moutsianas
L
,
Enciso-Mora
V
,
Ma
YP
, et al
.
Multiple Hodgkin lymphoma-associated loci within the HLA region at chromosome 6p21.3
.
Blood
.
2011
;
118
(
3
):
670
-
674
.
51.
Barsotti
AM
,
Beckerman
R
,
Laptenko
O
,
Huppi
K
,
Caplen
NJ
,
Prives
C
.
p53-dependent induction of PVT1 and miR-1204
.
J Biol Chem
.
2012
;
287
(
4
):
2509
-
2519
.
52.
Huppi
K
,
Volfovsky
N
,
Runfola
T
, et al
.
The identification of microRNAs in a genomically unstable region of human chromosome 8q24
.
Mol Cancer Res
.
2008
;
6
(
2
):
212
-
221
.
53.
Enciso-Mora
V
,
Broderick
P
,
Ma
Y
, et al
.
A genome-wide association study of Hodgkin’s lymphoma identifies new susceptibility loci at 2p16.1 (REL), 8q24.21 and 10p14 (GATA3)
.
Nat Genet
.
2010
;
42
(
12
):
1126
-
1130
.
54.
Shi
J
,
Zhang
Y
,
Zheng
W
, et al
;
kConFab Investigators
.
Fine-scale mapping of 8q24 locus identifies multiple independent risk variants for breast cancer
.
Int J Cancer
.
2016
;
139
(
6
):
1303
-
1317
.
55.
Aryankalayil
MJ
,
Chopra
S
,
Levin
J
, et al
.
Radiation-induced long noncoding RNAs in a mouse model after whole-body irradiation
.
Radiat Res
.
2018
;
189
(
3
):
251
-
263
.
56.
Rudolph
A
,
Song
M
,
Brook
MN
, et al
.
Joint associations of a polygenic risk score and environmental risk factors for breast cancer in the Breast Cancer Association Consortium
.
Int J Epidemiol
.
2018
;
47
(
2
):
526
-
536
.

Supplemental data

Sign in via your Institution