Abstract
Background:
A positive association between socio-economic position and breast cancer has been widely observed, but not hitherto within twin pairs, where shared familial factors were adjusted for.
Methods:
We used data on education and other factors from the Danish Twin Registry, The Danish Cancer Registry, and official registers on a total of 16 310 twins. Unpaired and intrapair Cox regression analyses were compared.
Results:
In the unpaired analysis, an educational gradient in breast cancer risk was found. Similar results were seen in the intrapair analyses of all twins, although no longer statistically significant. When intrapair analyses were stratified on zygosity, the effect of education was attenuated in the monozygotic twins.
Conclusion:
The main findings support an effect of education beyond shared familial factors.
Similar content being viewed by others
Main
In Denmark, breast cancer (BC) incidence has increased during recent decades, and it is now the most frequent cancer in women (Ferlay et al, 2007). In contrast to most health outcomes, a positive social gradient has been observed in BC, showing a higher risk in higher socio-economic position (SEP) groups (Kongevinas et al, 1997).
Several BC risk factors are well established: high age at first birth (AFB), low parity (Ma et al, 2006), hormonal replacement therapy (Porch et al, 2002; Tjonneland et al, 2004), high body mass index (van den Brandt et al, 2000; Ahlgren et al, 2004), and high alcohol consumption (Boyle and Boffetta, 2009). In recent years, there has also been interest in in utero- and early-life factors, including high birth weight (Silva et al, 2008; Oberg et al, 2009). Moreover, evidence suggests that child growth and childhood diet may convey the risk (Frazier et al, 2003; Ahlgren et al 2004; De Stavola et al, 2004; Michels et al, 2006). The BC aetiology seems to be influenced by different exposures that exert their influence throughout the life course (Lynch and Smith, 2005). Although modest, a genetic component of BC has also been demonstrated. In a large Nordic twin study, non-shared environment was found to influence the risk of BC the most (explaining 67% of the variance), while shared environmental factors seemed to have a negligible role (Lichtenstein et al, 2000). This may suggest that adult factors have a much larger effect than early-life factors. However, the power to disentangle genetic and shared environmental effects in heritability studies is limited.
The observed social differentials in BC may be driven by socially patterned adult factors such as AFB and parity, as better-educated women are more likely to postpone their pregnancies and have fewer children. Alternatively, the differentials may be a result of selection processes in which underlying factors, such as genes or childhood environment, influence SEP and the risk of BC, thus creating a spurious association (Smith et al, 1994). Previous studies on SEP and cancer have mostly dealt with mortality as opposed to incidence, which may complicate the picture if the social patterning of incidence and survival from cancer differs. However, a recent study found no effect of SEP in childhood on BC incidence when adult SEP was adjusted for (de Kok et al, 2008).
Twins provide a unique model to disentangle contributions from adult, childhood, and genetic factors as twins share their intrauterine- and early environment and are partly (dizygotic (DZ) twins) or fully (monozygotic (MZ) twins) matched on genetic setup. We aimed, by using education as an indicator of SEP, to separate the effect of SEP on BC from that of shared familial factors.
Materials and methods
In the Danish Twin Registry, all female twins born between 1921 and 1974 were considered for the study. After restrictions to sex-liked intact twin pairs of known zygosity, alive, and resident in Denmark at the beginning of follow-up (1 January 1980 or the date the participants turned 30 years if younger than 30 years in 1980) with no previous BC diagnosis and with full information on education, the final study population consisted of 16 310 twins, including 6268 MZ and 10 042 dizygotic same sex twins (DZSS).
Information on education came from Statistics Denmark and was operationalised as educational status according to the International Standard Classification of Education (Institute For Statistics, 2006) and length of education. Educational status was categorised into primary/lower secondary (primary) and upper secondary/tertiary education (secondary/tertiary) while length of education was included as a quantitative variable measured in years.
As parity and AFB are known to be strong determinants of BC and strongly correlated with education, we included these variables in a subanalysis to investigate the degree to which these factors might be driving the association between education and BC. Parity was constructed as a time-dependent variable (0, 1, 2+ births), while AFB was categorised as follows: no children, <26 years, 26–30 years, 31–35 years, >35 years. Owing to collinearity between the two variables, they were included in the model as a combined variable.
The Danish Cancer Registry contains information on all primary incident cancers in Denmark since 1943, except non-melanoma skin cancer. Reporting is mandatory and registered information comes from a cross-checking of different sources. Cases of BC were indentified through ICD-7 (1980–1993) and ICD-10 (1994–2006) codes (170 and C50, respectively).
Data analysis
Using Cox regression, person years of follow-up were accumulated from baseline (age at 1 January 1980 or age 30 years) to the age at the end of follow-up (age at 31 December 2006), diagnosis of BC, death, or emigration, whichever was earliest. The analyses were adjusted for calendar time. Unpaired analyses were compared with intrapair analyses. The intrapair analyses were carried out by stratifying on ‘pair-number’ (Holt and Prentice, 1974), which permits the baseline hazard to vary freely between pairs while at the same time fixing it within pairs. The hazard function for twin j in pair i is given by hij(t,z)=h0i(t) exp(βzij), where h0i(t) is the pair-specific baseline hazard and β is the effect of education that applies to all subjects.
Differences in genetic relatedness between MZ (genetically identical) and DZ twins (share on average 50% of segregating genes) could further be exploited to make inferences about genetic and environmental confounding, respectively; assuming an effect of education in the unpaired analyses, genetic confounding would be indicated if the intrapair analyses showed a partial attenuation of the association in DZSS twins and a full attenuation in the MZ twins. Similarly, confounding from shared environment would be likely to account for the association if a full attenuation were observed in both DZSS and MZ twins. Finally, persistence of the association would be expected if education had an independent effect on BC beyond the effect of shared familial factors. Hazard ratios (HRs) were calculated and 95% confidence limits were corrected for interdependence of observations within twin pairs. In the intrapair analyses, a strata option was included (StataCorp. 2009, Stata Statistical Software: Release 11, College Station, TX, USA and StataCorp LP, College Station, TX, USA). Potential interaction between education and zygosity was examined.
Results
In the cohort of 16 310 female twins, 518 developed BC during the follow-up period 1980–2006. As expected, more DZSS than MZ twins were discordant on education (Table 1). More twins were discordant by the quantitative definition, which expresses the difference in years of education within a twin pair. A larger mean difference in length of education was also observed among the DZSS twins. Additionally, a higher proportion of DZSS (44%) twins had a primary education compared with MZ twins (36%). The unpaired analysis of risk according to educational status (Table 2) shows a statistically significant increased risk associated with a secondary/tertiary education compared with primary education (HR=1.39, 95% confidence interval (CI)=1.15–1.67) for all twins. Similar results were seen when DZSS and MZ twins were analysed separately, and an analysis of potential interaction showed no evidence of modification of educational effect by zygosity (P=0.66). For the length of education, we observed a 4% increased risk of BC associated with one additional year of education (HR=1.04, 95% CI=1.01–1.07). The same pattern was seen when stratified on zygosity. The intrapair analysis of all twins showed an association with educational status of the same magnitude as in the unpaired analyses (HR=1.38), although no longer statistically significant (95% CI=0.95–2.01). When the analysis was split on zygosity, the association was attenuated in the MZ twins (HR=0.82, 95% CI=0.41–1.67). However, CIs were wide due to the limited number of events. For DZSS twins, the association persisted (HR=1.70, 95% CI=1.08–2.69). Similar trends were suggested in the analyses by length of education. No statistically significant interaction between education and zygosity was found in the intrapair analysis (P=0.35). When the analysis was adjusted for parity and AFB, the effect of education changed only negligibly (educational status: HR=1.37, 95% CI=1.14–1.65 and educational length: HR=1.04, 95% CI=1.01–1.07) and the effect of the combined covariate was statistically insignificant (P=0.61).
Discussion
In this study, we examined the association between education and BC risk in Danish female twins with and without adjusting for familial background factors.
In the unpaired analyses we found the expected increased risk associated with a secondary/tertiary compared with a primary educational status of similar magnitude to that previously reported (Kongevinas et al, 1997; Carlsen et al, 2008). This association was replicated in the intrapair analyses for all twins, suggesting that the effect of education is not due to confounding from familial factors. This accords with the Nordic twin study (Lichtenstein et al, 2000) that suggests that non-shared factors are most important risk factors for BC as in the study that showed no effect of childhood SEP on incidence, when adult factors were adjusted for (de Kok et al, 2008).
When MZ twins were analysed separately, the association was completely attenuated. However, as only few MZ twin pairs were available for the intrapair analysis, this is likely a chance finding. If not, the results are suggestive of underlying genetic factors associated with both a high education and a high risk of BC. The fact that no partial attenuation was observed in DZSS twins, contrary to expectation, might be explained by the presence of complex gene–gene interactions where several genes need to be present simultaneously in order for a phenotype to occur. Alternatively, if MZ twins also share their environment more closely than DZSS twins (Horwitz et al, 2003), for example fertility patterns, this could account for the observed differences in attenuation of effect between MZ and DZSS twins. However, this was not supported in our data set where parity and AFB did not have a statistically significant impact on the risk of BC.
This study is the first to use the discordant twin-pair design to examine social inequality in BC. Several previous studies have used this approach to address causality in the social patterning of health and disease, but the population size available for this study allowed us to investigate a hard end point such as BC, which presents an interesting hypothesis because of its curious positive association with SEP. The prospective design and almost complete information in the registers ensure valid information about educational status and BC in this largely unselected twin population.
A key assumption in the discordant twin-pair logic is that twins who have grown up together are matched on early environment. Survey data indicate that this was the case for about 98% of the population (Madsen et al, 2010), so we do not expect this to be an important source of bias. Selection is another relevant issue. Twins who are discordant on education may comprise a selected group as they have obtained different educational statuses in spite of their common background and may not be representative of the twin population at large. Nevertheless, the discordant twin-pair method offers an exceptional opportunity to enhance validity in non-experimental designs.
In conclusion, the finding of an effect of education in the combined intrapair analysis suggests that education has an effect on BC risk beyond shared familial factors; parity and AFB did not seem to be driving this association. The finding of an attenuation of effect in MZ twins, suggesting genetic confounding, may serve as a hypothesis for testing in larger data sets.
Change history
29 March 2012
This paper was modified 12 months after initial publication to switch to Creative Commons licence terms, as noted at publication
References
Ahlgren M, Melbye M, Wohlfahrt J, Sørensen TI (2004) Growth patterns and the risk of breast cancer in women. N Engl J Med 351: 1619–1626
Boyle P, Boffetta P (2009) Alcohol consumption and breast cancer risk. Breast Cancer Res 11 (Suppl 3): S3
Carlsen K, Høybye MT, Dalton SO, Tjønneland A (2008) Social inequality and incidence of and survival from breast cancer in a population-based study in Denmark, 1994–2003. Eur J Cancer 44: 1996–2002
de Kok IM, van Lenthe FJ, Avendano M, Louwman M, Coebergh JW, Mackenbach JP (2008) Childhood social class and cancer incidence: results of the globe study. Soc Sci Med 66: 1131–1139
De Stavola BL, dos Santos Silva I, McCormack V, Hardy RJ, Kuh DJ, Wadsworth ME (2004) Childhood growth and breast cancer. Am J Epidemiol 159: 671–682
Ferlay J, Autier P, Boniol M, Heanue M, Colombet M, Boyle P (2007) Estimates of the cancer incidence and mortality in Europe in 2006. Ann Oncol 18: 581–592
Frazier AL, Ryan CT, Rockett H, Willett WC, Colditz GA (2003) Adolescent diet and risk of breast cancer. Breast Cancer Res 5: R59–R64
Holt JD, Prentice RL (1974) Survival analyses in twin studies and matched pair experiments. Biometrika 61: 17–30
Horwitz AV, Videon TM, Schmitz MF, Davis D (2003) Rethinking twins and environments: possible social sources for assumed genetic influences in twin research. J Health Soc Behav 44: 111–129
Institute For Statistics (2006) ISCED 1997 International Standard Classification of Education. UNESCO: Montreal, Canada
Kongevinas M, Pearce N, Susser M, Boffetta P (1997) IARC Scientific Publications. International Agency for Research on Cancer: Lyon
Lichtenstein P, Holm NV, Verkasalo PK, Iliadou A, Kaprio J, Koskenvuo M, Pukkala E, Skytthe A, Hemminki K (2000) Environmental and heritable factors in the causation of cancer—analyses of cohorts of twins from Sweden, Denmark, and Finland. N Engl J Med 343: 78–85
Lynch J, Smith GD (2005) A life course approach to chronic disease epidemiology. Annu Rev Public Health 26: 1–35
Ma H, Bernstein L, Pike MC, Ursin G (2006) Reproductive factors and breast cancer risk according to joint estrogen and progesterone receptor status: a meta-analysis of epidemiological studies. Breast Cancer Res 8: R43
Madsen M, Andersen AM, Christensen K, Andersen PK, Osler M (2010) Does educational status impact adult mortality in Denmark? A twin approach. Am J Epidemiol 172 (2): 225–234
Michels KB, Rosner BA, Chumlea WC, Colditz GA, Willett WC (2006) Preschool diet and adult risk of breast cancer. Int J Cancer 118: 749–754
Oberg S, Cnattingius S, Sandin S, Lichtenstein P, Iliadou A (2009) Birth weight-breast cancer revisited: is the association confounded by familial factors? Cancer Epidemiol Biomarkers Prev 18: 2447–2452
Porch JV, Lee IM, Cook NR, Rexrode KM, Burin JE (2002) Estrogen-progestin replacement therapy and breast cancer risk: the Women's Health Study (United States). Cancer Causes Control 13: 847–854
Silva IS, De Stavola B, McCormack V (2008) Birth size and breast cancer risk: re-analysis of individual participant data from 32 studies. PLoS Med 5: e193
Smith GD, Blane D, Bartley M (1994) Explanations for socio-economic differentials in mortality. Evidence from Britain and elsewhere. Eur J Public Health 4: 131–144
Tjonneland A, Christensen J, Thomsen BL, Olsen A, Overvad K, Ewertz M, Mellemkjaer L (2004) Hormone replacement therapy in relation to breast carcinoma incidence rate ratios: a prospective Danish cohort study. Cancer 100: 2328–2337
van den Brandt PA, Spiegelman D, Yaun SS, Adami HO, Beeson L, Folsom AR, Fraser G, Goldbohm RA, Graham S, Kushi L, Marshall JR, Miller AB, Rohan T, Smith-Warner SA, Speizer FE, Willett WC, Wolk A, Hunter DJ (2000) Pooled analysis of prospective cohort studies on height, weight, and breast cancer risk. Am J Epidemiol 152: 514–527
Acknowledgements
The Danish Aging Research Center is supported by a grant from the VELUX Foundation (95-103-11419) and from the Danish Cancer Society (SU08001). We thank Dorthe Almind Pedersen, MSc for her valuable help with the management of the data.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Competing interests
The authors declare no conflict of interest.
Rights and permissions
From twelve months after its original publication, this work is licensed under the Creative Commons Attribution-NonCommercial-Share Alike 3.0 Unported License. To view a copy of this license, visit http://creativecommons.org/licenses/by-nc-sa/3.0/
About this article
Cite this article
Madsen, M., Andersen, P., Gerster, M. et al. Does the association of education with breast cancer replicate within twin pairs? A register-based study on Danish female twins. Br J Cancer 104, 520–523 (2011). https://doi.org/10.1038/sj.bjc.6606090
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/sj.bjc.6606090
Keywords
This article is cited by
-
Matched survival data in a co-twin control design
Lifetime Data Analysis (2014)