|Home | About | Journals | Submit | Contact Us | Français|
The online version of this article has been published under an open access model. Users are entitled to use, reproduce, disseminate, or display the open access version of this article for non-commercial purposes provided that: the original authorship is properly and fully attributed; the Journal and Oxford University Press are attributed as the original place of publication with the correct citation details given; if an article is subsequently reproduced or disseminated not in its entirety but only in part or as a derivative work this must be clearly indicated. For commercial re-use, please contact email@example.com
Mechanistic evidence linking chromosomal aberration (CA) to early stages of cancer has been recently supported by the results of epidemiological studies that associated CA frequency in peripheral lymphocytes of healthy individuals to future cancer incidence. To overcome the limitations of single studies and to evaluate the strength of this association, a pooled analysis was carried out. The pooled database included 11 national cohorts and a total of 22 358 cancer-free individuals who underwent genetic screening with CA for biomonitoring purposes during 1965–2002 and were followed up for cancer incidence and/or mortality for an average of 10.1 years; 368 cancer deaths and 675 incident cancer cases were observed. Subjects were classified within each laboratory according to tertiles of CA frequency. The relative risk (RR) of cancer was increased for subjects in the medium [RR = 1.31, 95% confidence interval (CI) = 1.07–1.60] and in the high (RR = 1.41; 95% CI = 1.16–1.72) tertiles when compared with the low tertile. This increase was mostly driven by chromosome-type aberrations. The presence of ring chromosomes increased the RR to 2.22 (95% CI = 1.34–3.68). The strongest association was found for stomach cancer [RRmedium = 1.17 (95% CI = 0.37–3.70), RRhigh = 3.13 (95% CI = 1.17–8.39)]. Exposure to carcinogens did not modify the effect of CA levels on overall cancer risk. These results reinforce the evidence of a link between CA frequency and cancer risk and provide novel information on the role of aberration subclass and cancer type.
Lifestyle factors, environmental exposures and genetic background interact in the development of many cancers. The use of phenotypic markers that integrate the impact of environmental factors and genetic processes would be of great value in risk assessment and surveillance of populations occupationally or environmentally exposed to carcinogens. To date, the strongest mechanistic and epidemiological evidence for such markers is available for chromosome damage (1–3).
The first modern description of a specific chromosome abnormality in a human tumour, i.e. the Philadelphia chromosome in chronic myeloid leukaemia, was published almost 50 years ago (4), but only with the introduction of banding techniques, the major role of chromosomal rearrangements in carcinogenesis was fully recognized (5). Subsequent observations that many genes affected by chromosomal rearrangements were involved in critical stages in cell growth, development or survival focused the interest of researchers on how the rearrangements alter the function of these target genes. The model linking chromosome breakage and rearrangement to the activation/deactivation of genes relevant to early stages of carcinogenesis was further refined in the 1990s, when the use of the fluorescence in situ hybridization technique in irradiated cells contributed to understanding the role of double-strand breaks in the formation of chromosomal aberrations (CAs) and the implications of homology-dependent and homology-independent pathways of double-strand break repair (1,6).
Although mechanistic evidence supporting the role of chromosomal alterations in early stages of cancer has been available for a long time, the first epidemiological data showing that the frequency of CA in peripheral lymphocytes may predict cancer incidence in human populations were published in early 1990s when two small cohort studies were performed in healthy subjects scored for CA in the Nordic countries and Italy (7–9). Later analyses based on those cohorts suggested that the cancer risk predictivity of CA is not modified by occupational exposure to carcinogens or tobacco smoking, supporting the hypothesis that the association between CA and cancer might be independent of exposure (10). Other cohort studies published thereafter have generally confirmed the results of the first reports, with some discrepancy in the strength of the associations between CA and cancer risk and in the independence from exposure to carcinogens (11,12).
The aim of this pooled analysis, based on original data of 11 national cohort studies from Europe already published, was to provide a semi-quantitative estimate of the association between CA and cancer risk, a measure that could be used for risk assessment. Other goals were to estimate the association with specific cancer sites and major subclasses of CA and the role of exposure to carcinogenic agents as potential effect modifiers of the association between CA and cancer risk.
The study design of all cohort studies included in the pooled analysis was similar to the original approach used by the Nordic study group in the first studies published in the early 1990s (13). In brief, cancer-free subjects who underwent CA screening for biomonitoring purposes—mostly due to occupational exposure to carcinogens and mutagens—were followed up from the date of the first CA test until the date of death, cancer diagnosis, emigration, 85th birthday or the end of follow-up (2000–2003, depending on the country), whichever occurred first. To standardize for interlaboratory variation, all subjects were classified within each laboratory according to tertiles of CA frequency as follows: low (1–33 percentile), medium (34–66 percentile) or high (67–100 percentile). The approach based on tertiles was the most efficient according to the Akaike information criteria (score = 0.70 for total CA and all cancers) when compared in the pooled database with other methods of standardization such as z scores (score = 0.72) and deviation from the mean (score = 0.72) (14).
Information on cancer incidence and mortality was collected through linkage with national/local cancer and cause of death registries. In Hungary and Poland, an active system of follow-up for mortality and cancer incidence was set up (12).
CA frequency was measured in cultured peripheral blood lymphocytes using the classification criteria of Savage (15). Only individuals with at least 100 scored metaphases were included in the database (>100 metaphases were scored for 15% of subjects and >200 for 1.3%). Gaps were not included. The culture time was 48 or 72 h. Tertiles for CA data based on 48 h culture time (89.7% of the whole series) were evaluated separately from those with 72 h culture time. More detail about the study design and the quality procedures implemented can be found in the original papers (7–12). The Italian database included an additional unpublished cohort of 391 subjects, among whom 9 cancer cases occurred at the end of the follow-up.
After obtaining approval by local institutional review boards, individual cytogenetic and epidemiological data from the original cohorts were collected and entered in a pooled database. The information gathered in each cohort included demographic data of the subject, the laboratory performing the assay, information on exposure to genotoxic agents and smoking habits. CA data included the date of the test; culture time; the number of cells scored and the number of chromatid breaks, dicentrics, chromosome breaks, chromatid exchanges, ring chromosomes, marker chromosomes and aberrant cells. Total CA frequency was defined as the number of cells with aberrations, excluding gaps, per 100 cells. Chromosome-type aberrations included chromosome-type breaks, ring chromosomes, marker chromosomes and dicentrics, whereas chromatid-type aberrations included chromatid-type breaks and chromatid exchanges. The criteria for inclusion of individuals in the pooled analysis were valid demographic data, at least 15 years of age and no cancer history at the time of the CA test. Whenever a subject was analysed more than once, only the first test was considered for frequency classification.
The pooled databases included a total of 22 358 subjects from 11 countries and 45 laboratories screened for CA between 1965 and 2002 and followed up for a median period of 10.1 years. The Czech data accounted for 48.3% of the total persons-years. Details about the pooled cohort are reported in Table I.
Non-melanoma skin cancer was not systematically registered in the original cohorts, and therefore, subjects suffering from this cancer were excluded from the pooled analysis. This change contributed to small differences observed in the number of subjects (24 subjects overall) and in the national risk estimates with respect to published original data.
Information on occupational exposure at the time of the CA test was available for most subjects, although the quality of this information was heterogeneous among the cohorts. A common classification scheme based on the exposure groups identified by Tinnerberg et al. (16) was adopted. The most frequent exposure was to ionizing radiation (IR), with 2974 subjects, followed by exposure to cytostatic drugs (2521 subjects) and polycyclic aromatic hydrocarbon (2004 subjects). The group of non-exposed subjects included as controls in some of the studies consisted of 3662 subjects. No data on occupational exposure were available for cohorts from Sweden and Finland.
Information on smoking status at the time of the test was available for 18766 subjects (83.9% of the total), among whom, 9657 were classified as ever smokers (43.2%). The amount of cigarettes smoked per day was registered only in a few cohorts and the quality of these data was not comparable among countries; therefore, this information was not used in the statistical analysis.
The pooled effect of CA frequency on cancer rates was estimated with a negative binomial regression analysis (17). All estimates were adjusted for the concomitant effect of age (5 years age-class), gender, time since test (0–4, 5–9, 10–14, 15–19 and ≥20 years), smoking habit (ever/never/former) and occupational exposure to carcinogens (yes/no). Relative risks (RRs) were estimated according to different levels of CA frequency together with its 95% confidence intervals (CIs). The presence of a linear trend of RRs was tested with the χ2 trend test.
To take into account the different magnitude of cancer rates in each country, a random intercept model was fitted to the data to estimate the variance component due to the between-countries variability. A random slope model was performed to verify if the CA effect on cancer rates varied between countries (18). These models allowed us to estimate the amount of heterogeneity present in the data and accordingly provided suitable estimates of the standard errors of the combined effects. A sensitivity analysis was performed to examine the variations in the pooled effect estimates due to the exclusion of one country at a time. Similarly, other sensitivity analyses were performed to evaluate the potential impact of an unbalanced distribution of main cohort characteristics such as culture time or the number of cells scored.
The potential effect modification of exposure to carcinogens on the main association between CA frequency and cancer was explored by means of the likelihood ratio test, i.e. comparing two hierarchical regression models for each major exposure category identified by the exposure matrix, with and without the interaction term, respectively (18).
To test the possibility that pre-clinical stages of cancer may have influenced CA frequency, we repeated the analysis excluding the first 2 years of follow-up and also tested a possible effect modification due to time since test.
Additional analyses were performed by country, by type of CA including subclasses and for cancer sites with at least 10 observed cases.
Kaplan–Meier product-limit method was applied to estimate the elapsed time without cancer at different CA levels. The proportional hazard assumption was not rejected based on the approach suggested by Grambsch et al. (19) [χ2(2 df)=1.14, P=0.57].
An overall number of 675 incident cancers and 368 cancer deaths were registered at the end of the follow-up. The two end points were analysed independently. The most frequent cancer sites were lung (99 incident cases and 100 deaths), colon and rectum (88 and 53), bladder and kidney (76 and 28), breast (76 and 11) and the lymphohaematopoietic system (46 and 30).
The results on the association between overall cancer incidence and CA frequency in the whole database are reported in Table II. The RR was increased for subjects in the medium (RR = 1.31, 95% CI = 1.07–1.60) and in the high (RR = 1.41, 95% CI = 1.16–1.72) tertiles of the total CA distribution when compared with the low tertile (test for linear trend, P= 0.001). An increase in RR was seen within the subclass of chromosome-type aberrations, i.e. 1.29 (95% CI = 1.05–1.59) and 1.42 (95% CI = 1.17–1.71) for medium and high tertiles (test for linear trend, P< 0.001), respectively, but not for chromatid-type aberrations. Among the subclasses of CAs, the RR for the presence of ring chromosomes was 2.22 (95% CI = 1.34–3.68), whereas the non-significant RRs of other types of chromosomal abnormalities (dicentrics, chromatid exchanges and markers) ranged between 1.30 and 1.39.
A similar pattern of risk was found for mortality. Subjects in the medium tertile of total CA distribution experienced an RR of 1.36 (95% CI = 1.04–1.81, deaths = 136) and those in the high tertile an RR of 1.35 (95% CI = 1.03–1.80, deaths = 154). The RR for the two main classes of CA followed the same pattern observed for cancer incidence. The presence of marker chromosomes entailed a mortality RR of 1.71 (95% CI = 0.97–3.02).
Although increased risks were observed in all cohorts but Finland (with the highest RRs in Hungary and in Norway), country-specific results showed a certain degree of heterogeneity, and the relationship between RRs of medium and high CA tertiles was equivocal in some countries (Figure 1). However, fitting a random slope model to the data did not provide strong evidence of between-country heterogeneity in CA effect [χ2 = 1.92 (2 df), P= 0. 38]. Furthermore, to evaluate if the overall association between CA and cancer risk was driven by a specific country, a sensitivity analysis was performed, sequentially removing a country from the analysis. The most extreme difference in RRs was obtained by removing the Czech Republic (increase to RRmedium = 1.32, RRhigh = 1.64) and Hungary (decrease to RRmedium = 1.25, RRhigh = 1.29) from the pooled analyses. In all cases, RRs remained significantly increased.
No covariate tested as potential effect modifier, i.e. gender, smoking habit, age, time from test and occupational exposure, was found to significantly modify the effect of CA levels on overall cancer risk in the whole cohort. In the analysis excluding the first 2 years of observation, the results on the association between CA and cancer incidence were not modified [597 cancer cases, RRmedium = 1.33 (95% CI = 1.07–1.66), RRhigh = 1.46 (95% CI = 1.18–1.81)]. Similarly, the analyses stratified by median duration of follow-up did not yield different results in the two periods (data not shown).
In the analysis of specific types of cancer, the strongest association was found for stomach cancer [32 cases, RRmedium = 1.17 (95% CI = 0.37–3.70), RRhigh = 3.13 (95% CI = 1.17–8.39); test for linear trend, P = 0.009], whereas the results for other cancer sites were only suggestive of an association (Table III). Analysis based on cancer mortality showed comparable results except for neoplasms of the lymphohaematopoietic system, for which an association with CA frequency [30 cases, RRmedium = 2.39 (95% CI = 0.77–7.45), RRhigh = 2.61 (95% CI = 0.85–8.02)] was suggested. No country-specific effect was found for any cancer site, but this analysis was hampered by small numbers.
As a whole, subjects exposed to genotoxic agents did not present a significant increase of cancer risk when compared with non-exposed subjects [RR = 1.08 (95% CI = 0.86–1.36)]; however, some exposure groups showed increased risks for specific cancer sites, such as lung cancer in workers exposed to IR [RR = 2.71 (95% CI = 1.14–6.43)]. The association of cancer risk with CA frequency was particularly evident in workers exposed to IR [2974 subjects, RRmedium = 2.01 (95% CI = 1.15–3.51), RRhigh = 1.81 (95% CI = 1.03–3.20)] and at a lower extent in groups exposed to polyvinyl chloride and plastics [945 subjects, RRmedium = 1.62 (95% CI = 0.53–4.96), RRhigh = 1.73 (95% CI = 0.58–5.19)], metals [662 subjects, RRmedium = 3.10 (95% CI = 0.64–14.97), RRhigh = 2.75 (95% CI = 0.59–12.80)] and welding fumes [392 subjects, RRmedium = 2.01 (95% CI = 0.42–9.51), RRhigh = 2.63 (95% CI = 0.56–12.28)]. The RRs among subjects not exposed to occupational carcinogens were 1.39 (medium, 95% CI = 0.86–2.25) and 1.54 (high, 95% CI = 0.90–2.64). Each occupational category was tested as a possible effect modifier of the association between CA and overall cancer risk, and none of the interaction terms were significant. The presence of effect modification due to occupational exposure to carcinogens was tested also for subclasses of CA and for major groups of cancer. All interactions tested turned out to be non-significant.
A similar analysis was performed to test the effect modification of cigarette smoking. When compared with non-smokers, ever smokers experienced an increased risk of all cancers (RR = 1.53, 95% CI = 1.28–1.83) and especially lung cancer (RR = 6.65, 95% CI = 3.17–13.95). The association between CA and risk of cancer was stronger in ever smokers [316 cases, RRmedium = 1.42 (95% CI = 1.04–1.94), RRhigh = 1.54 (95% CI = 1.14–2.08); test for linear trend, P = 0.008] than in non-smokers [208 cases, RRmedium = 1.30 (95% CI = 0.93–1.82), RRhigh = 1.18 (95% CI = 0.83–1.68)], but the interaction was not statistically significant (P = 0.53).
The association between cancer risk and CA frequency was consistent over the follow-up time, as shown by the Kaplan–Meier curves of the probability to be cancer free (Figure 2), which indicated a significantly worse outcome in subjects in the medium and high CA tertiles (log-rank test P<.001). Survival curves based on cancer mortality showed similar results (data not shown).
The results of the present study, which included data from all CA cohorts published so far, confirm that CA frequency measured in peripheral blood lymphocytes of healthy subjects is associated with cancer risk, although the strength of this conclusion is somewhat weakened by the lack of a clear dose–response in most national cohorts. The independence of these findings from the time elapsed between cytogenetic testing and cancer diagnosis supports the assumption that increased CA levels in high-risk subjects are not attributable to clinically undetected cancers.
The pooled approach applied in this study gave several advantages over a meta-analysis of published data, including the possibility to perform interaction tests and subgroup analyses. Nevertheless, pooling data of biomarkers such as CA, which are characterized by variability of protocols, equipment and methods, may pose a serious question of heterogeneity. To address this issue, data from each laboratory were categorized in three frequency categories, an approach which according to the Akaike information criteria was more efficient than other methods of standardization.
A possible source of misclassification is the use of a single CA measure to classify subjects according to CA frequency. Due to the availability of a repeated measure for 11909 subjects in the Czech cohort (99.3%), we could evaluate the agreement between the two measures in this database. 76.1% of these subjects were classified in the same category for both tests, whereas for 11.6 and 12.2%, the second test would have resulted in a higher or lower class compared with the first test, respectively. To evaluate if the choice of the test may have influenced the strength or the direction of the association, we compared risk estimates based on tertiles calculated with the first and the second assay. The results were very similar, i.e. RRmedium = 1.27 and 1.24 and RRhigh = 1.17 and 1.21, respectively, showing that the classification based on a single measure of CA—at least in the Czech cohort—did not affect the estimate of the association between CA and cancer.
The interpretation of the study findings at country level is complex. Although a dose–response was present only in Hungary and Norway (test for linear trend, P = 0.002 and P = 0.030, respectively), the risk estimates were above unity in most cohorts, confirming the importance of a pooled design to reach the adequate statistical power to test this hypothesis.
The large size of the pooled database allowed us to investigate in greater detail hypotheses not properly tested in the original studies, in particular, whether the frequency of chromosome-type and chromatid-type aberrations has a different predictive value for cancer risk. The question is of great relevance since the induction of these aberrations reflects different mechanisms. In peripheral lymphocytes, chromosome-type CAs are induced by IR and other S-independent clastogens, causing DNA double-strand breaks, whereas chromatid-type CAs are produced by S-dependent agents such as many chemical mutagens. The published evidence is contrasting since earlier Nordic–Italian studies found that both subclasses are equally predictive of cancer risk (22), whereas the larger cohorts from the Czech Republic (11) and Central Eastern European countries (12) suggest that chromosome-type CA may show a stronger association. The results of the pooled estimates reflect the prevailing effect found in these latter cohorts. However, the limited improvement over the results based on total CA does not support the use of chromosome-type aberrations as the most predictive biomarker of risk.
The increased cancer risk observed in subjects bearing ring chromosomes (20 cases) and the borderline increase observed for marker chromosomes (limited to mortality) are not circumstantiated enough to be considered for future risk assessment protocols, and chance of error measurement remain plausible explanations. However, these findings should further be considered, given the mechanistic evidence linking chromosomal rearrangements to early stages of cancer.
Information about subclasses of chromosome- and chromatid-type aberrations was available only for a part of the cohort. To check if this subgroup was representative of the whole cohort, we re-evaluated cancer risks associated with total CA frequency in the same subset of subjects. The results were in agreement with a random selection of the population, showing risks very similar to those of the whole cohort [RRmedium = 1.40 (95% = 1.01–1.92), RRhigh = 1.67 (95% CI = 1.22–2.28)].
The hypothesis that CA frequency is associated with specific cancers has been investigated since the first cohort studies. A stronger association was initially described with lung cancer and haematological malignancies (8), whereas more recently, based on larger numbers, both the Czech cohort (11) and the Central and Eastern European country cohort (12) reported an increased risk of stomach cancer. The analysis of this pooled databases showed the association of CA frequency with stomach cancer as the only cancer site significantly associated with CA frequency (Table III). The hypothesis of a spurious association due to the number of comparisons cannot be ruled out; however, there is supportive mechanistic evidence. The presence of chromosome instability among the most relevant events in intestinal and diffuse stomach cancers has extensively been described (23,24) and linked to the metabolisms of agents involved in stomach carcinogenesis, such as folic acid and vitamin B12 (25,26). Interestingly, a recently published cohort study linking the frequency of micronuclei in lymphocytes of healthy subjects to the risk of cancer reported stomach cancer among the sites more specifically associated with micronuclei frequency (27). A recent report indicated that Helicobacter pylori infection is associated with an increased level of micronuclei in peripheral lymphocytes (28).
Another open issue is the role of exposures to genotoxic and carcinogenic agents experienced by study subjects at the time of CA testing. A nested case–control study carried out within the Nordic–Italian databases to specifically test the presence of interaction of CA frequency with major occupational exposure to carcinogens and with cigarette smoking did not find any difference in the risks (10). However, this independence was later put in doubt by the observation that in a subcohort of Czech miners exposed to radon, subjects with high levels of CA showed a much higher risk of cancer than those in other occupational groups of the same cohort (29). In the present study, the non-homogeneous quality of exposure assessment among studies limited the possibility to test by the presence of interaction. However, with an exploratory approach, study subjects were grouped by genotoxic exposure and risks were estimated within each subcohort. The small increase of cancer risk in the whole group of workers exposed to carcinogens when compared with non-exposed workers can be explained by exposure misclassification. The reliability of the design was confirmed by the high risks found for well-known associations between carcinogenic agents and tumour sites, such as those reported previously for lung cancer in exposure to IR. No clear evidence of effect modification was found for any exposure, although there was a suggestion of a higher risk associated with CA frequency for subjects exposed to IR. The estimate of cancer risk associated with smoking habit was more reliable, as supported by the close to expected RRs found in smokers for all cancers and lung cancer. However, also for this exposure, results of the interaction tests—although affected by a limited statistical power—are in agreement with the hypothesis that exposure is not an effect modifier of the association between CA and cancer.
The results of this study represent the synthesis of research conducted over two decades aiming at validating biomarkers of chromosomal damage as predictors of cancer occurrence in groups of healthy individuals. One of the indirect results of this effort was the establishment of an international database of over 20000 subjects screened for CA, which can be followed up to provide further and more reliable data regarding the use of this biomarker to identify high-risk populations.
European Union (Cytogenetic Biomarkers and Human Cancer Risk: QLK4-CT-2000-00628 and QLK4-CT-2002-02831) to S.B., H.N., U.S., A.C.-W., E.F., A.F., S.G., I.-L.H., L.E.K., J.L., P.R., R.J.S. and P.B.; Environmental Cancef Risk, Nutrition and Individual Susceptibility Network of Excellence (513943) to R.V. and P.B.; the Italian Association for Cancer Research to S.B.; the Italian Space Agency to S.B.; the National Cancer Institute of the USA to R.V., A.Z. and A.F.; the Finnish Work Environment Fund to H.N.; and the Czech Ministry of the Environment (VaV-SL/5/160/05) to P.R. and R.J.S.
The Cancer Registry of Norway is not responsible for the analysis or interpretation of the Norwegian cancer data presented. We wish to acknowledge the important contribution to this research by the late Lars Hagmar, Lund University, Sweden, who authored seminal cohort studies on the link between CA and cancer.
Conflict of Interest Statement: None declared.