Primary DNA damage and genetic polymorphisms for CYP1A1, EPHX and GSTM1 in workers at a graphite electrode manufacturing plant
© Moretti et al; licensee BioMed Central Ltd. 2007
Received: 18 June 2007
Accepted: 01 October 2007
Published: 01 October 2007
The results of a cross-sectional study aimed to evaluate whether genetic polymorphisms (biomarkers of susceptibility) for CYP1A1, EPHX and GSTM1 genes that affect polycyclic aromatic hydrocarbons (PAH) activation and detoxification might influence the extent of primary DNA damage (biomarker of biologically effective dose) in PAH exposed workers are presented. PAH-exposure of the study populations was assessed by determining the concentration of 1-hydroxypyrene (1OHP) in urine samples (biomarker of exposure dose).
The exposed group consisted of workers (n = 109) at a graphite electrode manufacturing plant, occupationally exposed to PAH. Urinary 1OHP was measured by HPLC. Primary DNA damage was evaluated by the alkaline comet assay in peripheral blood leukocytes. Genetic polymorphisms for CYP1A1, EPHX and GSTM1 were determined by PCR or PCR/RFLP analysis.
1OHP and primary DNA damage were significantly higher in electrode workers compared to reference subjects. Moreover, categorization of subjects as normal or outlier highlighted an increased genotoxic risk OR = 2.59 (CI95% 1.32–5.05) associated to exposure to PAH. Polymorphisms in EPHX exons 3 and 4 was associated to higher urinary concentrations of 1OHP, whereas none of the genotypes analyzed (CYP1A1, EPHX, and GSTM1) had any significant influence on primary DNA damage as evaluated by the comet assay.
The outcomes of the present study show that molecular epidemiology approaches (i.e. cross-sectional studies of genotoxicity biomarkers) can play a role in identifying common genetic risk factors, also attempting to associate the effects with measured exposure data. Moreover, categorization of subjects as normal or outlier allowed the evaluation of the association between occupational exposure to PAH and DNA damage highlighting an increased genotoxic risk.
Graphite electrodes are used in the metal industry, mainly in electric arc furnace steel production, to refine steel in ladle furnaces, and in other smelting processes . The production process of graphite electrodes involves the use of coal tar, coal tar pitch and petroleum coke with workers possibly exposed to polycyclic aromatic hydrocarbons (PAH) by inhalation (of PAH both volatile and bound to respirable particulate matter) and dermal contact . Occupational exposure of workers to PAH-containing coal tar pitch volatiles, pitch and coke occurs during the manufacturing process particularly during the heating of raw materials [3–6]. Occupational exposure to PAH was associated with increased risk of developing lung, skin, bladder and prostate cancer among graphite electrode manufacturing workers [2, 7–10].
Several hundred PAH have been characterized for their chemistry and many individual PAH (e.g. benzo [a]pyrene, benzo [a]anthracene, dibenzo [a, h]anthracene) are regarded as probably or possibly carcinogenic to humans by the International Agency for Research on Cancer (group 2A or 2B) . The carcinogenic activity of PAH is related to the DNA damaging activity of some of their metabolites which can covalently bind to nucleophilic residues of DNA bases. The PAH are activated to the corresponding electrophilic diolepoxides (ultimate carcinogens) via cytochrome P450 (CYP) enzyme system . The inactivation of diolepoxide metabolites occurs mainly through conjugation with reduced glutathione (GSH) by glutathione S-transferase (GST) enzymes .
Many genes encoding carcinogens metabolising enzymes have been found to be polymorphic in humans, with relevance to the individual response to carcinogens, probably acting as modifiers of exposure biomarkers (susceptibility markers) [14–16].
A critical polymorphic gene which contributes to the bioactivation of many PAH is the CYP1A1 encoding the cytochrome P4501A1 enzyme (CYP1A1), an inducible enzyme with aryl hydrocarbon hydroxylase activity. Among several polymorphisms identified in the CYP1A1 gene, two closely linked mutations have been extensively studied in relation to cancer risk. The CYP1A1 Ile/Val (m2) mutation in the heme-binding region doubles the microsomal enzyme activity and it is in linkage disequilibrium in Caucasians with the CYP1A1 MspI (m1) mutation, that has also been associated experimentally with increased catalytic activity . Positive associations between the presence of these variant alleles and increased PAH-DNA adducts have been reported [18–20].
A further critical enzyme in PAH metabolism is microsomal epoxide hydrolase (mEH) which catalyses the hydrolysis of epoxides into dihydrodiols, and as such plays an important role in the detoxification of toxic, highly reactive, intermediates formed by cytochrome P450-mediated reactions . Low mEH activity has been associated with adverse drug responses or diseases states . Two polymorphic sites within the mEH gene (EPHX) have been identified. A substitution at codon 113 (Tyr → His) in exon 3 is associated with decreased mEH activity, whereas a substitution at codon 139 (His → Arg) in exon 4 is associated with increased mEH activity. These mutations affect mEH enzyme activity by altering protein stability without affecting the specific activity .
The GSTM1 gene encodes for the cytosolic enzyme glutathione S-transferase μ1 (GSTM1) that detoxify activated forms of chemical carcinogens such as polyaromatic hydrocarbon epoxides. This gene is deleted in about 50% of Caucasians, with a reported variation of from 38 to 65% . The inherited absence of the GSTM1 gene (the GSTM1 null genotype) is therefore theoretically associated at a higher risk to the toxic effects of chemicals and its influence on various biomarkers of exposure has been widely studied [18, 25–27].
In the present paper we reported the results of a cross-sectional study aimed to evaluate whether genetic polymorphisms (biomarkers of susceptibility) for CYP1A1, EPHX and GSTM1 genes that affect PAH activation and detoxification might influence the extent of primary DNA damage (biomarker of biologically effective dose) in PAH exposed workers (n = 109) and in unexposed controls (n = 82). PAH-exposure of the study populations was assessed by determining the concentration of 1-hydroxypyrene (1OHP) in urine samples (biomarker of exposure dose). The urinary concentration of 1OHP, a metabolite of the non carcinogenic PAH pyrene, is a well validated marker for PAH-exposure  giving an accurate assessment of total PAH exposure from all exposure routes . The extent of primary DNA damage was evaluated with the comet assay in peripheral blood leukocytes (PBL) . The study subjects were genotyped for polymorphisms in the CYP1A1 gene (mutations associated with increased catalytic activity) encoding the cytochrome P4501A1 enzyme, in the EPHX gene (mutations influencing the mEH enzyme activity) encoding for the mEH enzyme, and the GSTM1 gene (inherited absence) encoding for the cytosolic glutathione S-transferase μ1 (GSTM1).
Study Population and Samples Collection
A total of 191 healthy men living in the same area of Central Italy were enrolled in the study. The exposed group was composed of 109 workers at a graphite electrode producing plant, occupationally exposed to PAH. The electrode plant produces large electrodes (length 2.5÷3 m, diameter up to 1 m) and small bricks (diameter less than 0.1 m) of electrode paste for steel industry . As a control group, 82 reference subjects (not occupationally exposed to PAH) were recruited from the technical and maintenance staff of the University of Perugia.
Workers and reference subjects were interviewed to obtain personal data and information on current job, smoking habit, alcohol consumption, diet, and current and past health status (including information on medicine intake, X-ray examinations, and viral infections). Following informed consent was obtained from all individuals enrolled in the study, workers and reference subjects provided a urine sample for 1OHP determination and a peripheral venous blood sample (using heparinized vacutainer tubes) for the comet assay and genotyping. Biological samples were collected from workers at the end of the shift after at least 4 consecutive days at work. The samples were coded and immediately transferred to the laboratory in refrigerated boxes.
The research protocol was approved by the Ethics Committee of the University of Perugia.
CYP1A1, EPHX, and GSTM1 Genotyping
Total DNA was extracted from peripheral blood cells using standard techniques. All genotypes were determined after gene amplification using polymerase chain reactions (PCR) . The T → C mutation (m1) in the 3'-flanking region of CYP1A1 gene was detected by PCR followed by restriction fragment length polymorphism (RFLP) analysis using the restriction enzyme MspI . The genotype w1/w1 (lacking the MspI site) forms an uncleaved 340 bp band, while the genotype m1/m1 (homozygous for the allele carrying the mutation with the MspI site) generates two bands of 200 and 140 bp. The heterozygous genotype m1/w1 corresponds to three bands of 340, 200, and 140 bp.
The CYP1A1 Ile/Val replacement (m2) was detected by BsrDI-RFLP analysis . The Ile/Val polymorphism arises from a A → G base change resulting in the replacement of isoleucine by valine at residue 462 in the heme binding region of the enzyme. The Val allele variant shows an almost two-fold higher catalytic enzyme activity than Ile form.
Amplicons of exons 3 and 4 of EPHX gene (162, and 381 bp, respectively) were obtained by PCR, RFLP digestions were then performed to determine the exon 3 (Tyr113His) and exon 4 (His139Arg) genotypes, using the restriction enzymes EcoRV and RsaI, respectively . On the basis of the polymorphisms at codon 113 (exon 3) and 139 (exon 4) of EPHX gene, the subjects were classified according to expected mEH enzyme activity (low mEH, intermediate mEH, or high mEH activity) .
GSTM1 genotyping for gene deletions was carried out by detecting the presence or the absence of the intact gene . The absence of GSTM1 specific amplification products revealed the corresponding null genotype (homozygous deletion of the GSTM1 gene, resulting in deficiency of GSTM1 activity). The GSTM1 positive genotype, detected by the presence of GSTM1 specific band of 215 bp, contained wild-type homozygotes and heterozygotes for the deletion (not differentiated in the analysis), both expressing GSTM1 enzyme. Co-amplification of β-globin gene was used as an internal control (presence of amplifiable DNA in the sample).
Analysis of 1-Hydroxypyrene in Urine
Urinary concentrations of 1OHP were determined by HPLC in enzymatically hydrolyzed urine samples . Urine samples, adjusted to pH 5.0, were treated overnight at 37°C with β-glucuronidase and aryl sulfatase and then purified with solid phase extraction with Sep-Pack C18 cartridges primed with methanol. The cartridges were then washed with high purity water and 1OHP was eluted with methanol. The eluate was gently evaporated to dryness under nitrogen and reconstituted in methanol. Of the reconstituted eluate, 15 μl were injected into an HPLC and 1OHP, eluting at a retention time of 8 min, was detected with excitation and emission wavelengths of 347 and 388 nm, respectively.
Analysis of Primary DNA Damage (Comet Assay) in Leukocytes
PBL were obtained from whole blood by lysis of erythrocytes . Viability of cells after isolation was determined by the fluorochrome-mediated (simultaneous staining with fluorescein diacetate and propidium iodide) viability test . Isolated PBL were processed in the comet assay following the standard alkaline protocol , with minor modification [38, 39].
The cells (2 × 105) were mixed with 0.7% low melting temperature agarose (total volume 75 μl/slide) and sandwiched between a layer of 0.5% normal melting temperature agarose (75 μl) and a top layer of 0.7% low melting temperature agarose (65 μl) onto conventional microscope slides. Lysis of cellular and nuclear membranes of the embedded cells was performed by immersing the slides for 60 min, at 4°C in the dark, in ice-cold freshly prepared lysis solution (10 mM Tris-HCl, 1% sodium N-lauroylsarcosinate, 2.5 M NaCI, 100 mM Na2EDTA, 1% triton X-100, and 10% DMSO; pH 10). The slides were removed from the lysis solution and then placed on a horizontal electrophoresis box. The unit was filled with freshly made alkaline buffer (300 mM NaOH, 1 mM Na2EDTA; pH > 13) to a level of 0.25 cm over the slides. To allow DNA unwinding and expression of alkali labile damage, the embedded cells were exposed to alkali for 20 min, then the electrophoresis was performed in the same buffer for 20 min by applying an electric field of 25 V (1 V/cm) and adjusting the current to 300 mA.
To control the assay conditions, particularly slides preparation procedure and electrophoresis efficiency, negative and positive internal controls (Jürkat cells, human lymphoblastoid T-cells) were processed in parallel with whole blood samples. Jürkat cells were untreated (negative control) or incubated for 1 h with 1 μg/ml 4-nitroquinoline-N-oxide (positive control). Electrophoresis runs were considered valid only if the internal controls yielded the expected results.
After electrophoresis, the slides were first washed gently with 0.4 M Tris-HCl buffer (pH 7.5) to neutralize the alkali, and the DNA was then stained by adding 100 μl of ethidium bromide (2 μg/ml). The slides were kept in a humidified sealed box to prevent drying of the gel and analyzed within 48–72 hours.
Comets in each gel were analyzed (blind) at 500× magnification using an epi-fluorescent microscope (excitation filter, 515–560 nm; barrier filter, 590 nm) equipped with a high sensitivity black and white CCD camera. For each subject, the average tail length, tail intensity and tail moment values were determined scoring 150 comets (50 comets/slide, from at least three replicate slides). Imaging was performed using a specialized analysis system ("Comet Assay III", Perceptive Instruments).
Statistical Analysis of Data
The analyses were carried out using the SPSS 10.0 statistical software package (SPSS Inc., IL, USA). For each subject, the averaged migration extents (i.e. tail length, tail intensity, and tail moment) of the 150 cells analyzed were used as a data point in data analysis. Cells were also classified as either "undamaged" or "damaged" by considering threshold levels indicating the cells with abnormal size tail (AST) (i.e. the 95th percentile of the distribution of the tail parameters among controls) . Cells with tail parameters values below the cut-off (18.87 μm, 16.73%, and 1.90 for tail length, tail intensity, and tail moment, respectively) were classified as "undamaged", and those with higher values as "damaged". Based on these binary outcomes, each subject was classified as "normal" or "outlier" by considering the upper bound of the expected number of AST which was calculated for each subject from the binomial distribution . The latter outcome variable (i.e. each subject defined as normal or outlier) is again binary, becoming the subject the statistical unit.
Group differences in concentration of 1OHP and extent of primary DNA damage (individual averaged tail parameters) were tested with the non-parametric Mann-Whitney U-test (two-tailed). The Hardy-Weinberg equilibrium test for genotype distribution was performed using a χ2 test with 1 degree of freedom. The GSTM1 genotype was coded as positive (wild-type homozygotes and heterozygotes for the deletion) or null (homozygous deletion), making direct calculation of Hardy-Weinberg equilibrium impossible. The Pearson-χ2 test was used to determine significant differences in the distribution of outliers subjects and of allele frequencies of the considered genotypes among the groups. Outlier subjects prevalence odds ratio (OR) with 95% confidence intervals (95% CI) were calculated by means of cross tabulation. The level of statistical significance was set at p < 0.05.
Characteristics of the Study population
Demographic characteristics of the study population grouped according to exposure status.
43.32 ± 0.56
42.99 ± 0.73
≤ 40 years
> 40 years
Duration of employmentb
14.34 ± 8.05
≤ 10 years
> 10 years
18.85 ± 9.14
14.92 ± 8.11
Distribution of CYP1A1, EPHX, and GSTM1 genotypes in exposed and control subjects. Data are reported as the number of subjects (percentage between brackets).
CYP1A1 (MspI)b w1/w1
w1/m1 + m1/m1
CYP1A1 (Ile/Val)c w2/w2
w2/m2 + m2/m2
EPHX d Low
Medium + High
Influence of Exposure on the Biological End-Points
Statistically significant correlations were found between the three different parameters measuring the extent of DNA damage in leukocytes (i.e. tail length, tail intensity and tail moment), thus data presentation will be limited to tail intensity, also taking into account that relative tail intensity (percentage of DNA migrated in the comet tail) is considered to be the most useful parameter, as it bears a linear relationship to strand-breaks frequency and is relatively unaffected by threshold settings in the computerized analysis system .
Urinary concentration of 1-hydroxypyrene (1OHP) and extent of primary DNA damage in peripheral blood leukocytes in exposed workers and control subjects. Data reported as the group mean values (± SEM) of individual: urinary concentration of 1OHP (expressed as μg 1OHP/g creatinine), averaged tail intensity values (% of DNA migrated in the comet tail evaluated in 150 cells) and number of AST.
2.64 ± 0.29*
5,28 ± 0,21*
11,75 ± 0,81*
≤ 40 years
2,96 ± 0,43*
4,63 ± 0,36
9,59 ± 1,35
> 40 years
2,49 ± 0,37*
5,58 ± 0,26*
12,73 ± 0,99*
Employment ≤ 10 years
2,23 ± 0,34
4,90 ± 0,29
10,14 ± 1,16
Employment > 10 years
2,84 ± 0,40
5,47 ± 0,28
12,55 ± 1,05
2.85 ± 0.47*
5.47 ± 0.30*
12,81 ± 1,14*
2.43 ± 0.34*
5.11 ± 0.31
10,71 ± 1,14*
0.15 ± 0.02
4,33 ± 0,22
7,50 ± 0,98
≤ 40 years
0.17 ± 0.03
4.34 ± 0.35
6,83 ± 1,51
> 40 years
0,14 ± 0,02
4,33 ± 0,28
7,88 ± 1,29
0.10 ± 0.02
4.23 ± 0.31
8,36 ± 1,60
0.20 ± 0.02§
4.44 ± 0.32
6,60 ± 1,12
Proportions of outlier subjects in exposed and control groups. Data are referred to tail intensity (% of DNA migrated in the comet tail).
OR (95% CI)
≤ 40 years
> 40 years
Employment ≤ 10 years
Employment > 10 years
Urinary 1OHP was significantly higher in electrode workers than in reference subjects, with 1OHP concentrations about 18-fold higher in exposed workers than in controls (p < 0.001). Smoking resulted in a statistically significant increase of the 1OHP levels in controls (p < 0.001), but not in the exposed group. Moreover, urinary 1OHP was significantly lower (p = 0.046) in electrode workers older than 40 respect to the youngest subjects. Urinary 1OHP concentrations found in the present study correlate (r = 0.732) with the exposure to total PAH (data not shown).
Factory workers showed a statistical significant increase in averaged DNA damage over that for controls (p = 0.002). Smoking habit did not increase the extent of DNA damage, both in the exposed and in the reference groups. Primary DNA damage was significantly higher (p = 0.036) in electrode workers older than 40 respect to the youngest subjects.
The group mean value of damaged cells (AST) is significantly higher in the exposed than in control subjects (p < 0.001). Moreover, the number of 'outlier' subjects (subjects showing a high number of AST) is significantly higher in the exposed (38.5%) than in the control group (19.5%), with an OR = 2.59 (CI95% 1.32–5.05). In the group of exposed workers, the higher frequency of AST and of outliers is mainly accounted by the subjects with an age over 40 years, with an OR = 2.63 (CI95% 1.17–5.90).
Influence of Genotypes on the Biological End-Points
Urinary concentration of 1-hydroxypyrene (1OHP) and extent of primary DNA damage in peripheral blood leukocytes in exposed workers and control subjects with respect to metabolic genotypes (CYP1A1, EPHX, and GSTM1). Data reported as the group mean values (± SEM) of individual: urinary concentration of 1OHP (expressed as μg 1OHP/g creatinine), averaged tail intensity values (% of DNA migrated in the comet tail evaluated in 150 cells) and number of AST.
CYP1A1 (MspI)b w1/w1
2.37 ± 0.29*
5.08 ± 2.45*
11.23 ± 0.92*
w1/m1 + m1/m1
3.33 ± 0.71*
5.81 ± 0.43
13.13 ± 1.65*
CYP1A1 (Ile/Val)c w2/w2
2.71 ± 0.31*
5.22 ± 0.22*
11.51 ± 0.84*
w2/m2 + m2/m2
1.96 ± 0.61*
5.85 ± 0.75
13.91 ± 2.95
EPHX d Low
2.17 ± 0.37*,§
5.22 ± 0.30
12.63 ± 1.12*
Medium + High
3.19 ± 0.45*
5.00 ± 0.30*
10.72 ± 1.16*
2.53 ± 0.38*
5.46 ± 0.29*
12.14 ± 1,09*
2.79 ± 0.45*
5.05 ± 0.32
11.22 ± 1,21*
CYP1A1 (MspI)b w1/w1
0.14 ± 0.02
4.28 ± 0.25
7.24 ± 1.10
w1/m1 + m1/m1
0.21 ± 0.05
4.87 ± 0.53
6.54 ± 1.30
CYP1A1 (Ile/Val)c w2/w2
0.14 ± 0.02
4.38 ± 0.24
6.77 ± 0.91
w2/m2 + m2/m2
0.20 ± 0.07
4.25 ± 0.30
10.86 ± 5.08
EPHX d Low
0.13 ± 0.02
4.58 ± 0.33
7.34 ± 1.13
Medium + High
0.16 ± 0.02
4.21 ± 0.30
6.96 ± 1.43
0.15 ± 0.02
4.49 ± 0.32
8.32 ± 1,48
0.15 ± 0.03
4.23 ± 0.31
5.67 ± 1,02
This study was designed to evaluate the occupational exposure to PAH, as assessed by biological monitoring techniques, in workers employed at a graphite electrode manufacturing plant and in reference subjects with presumed lower exposure to PAH. Current exposure to PAH was assessed by determining the urinary concentration of 1OHP (biomarker of internal dose), whereas biological effect monitoring focused on the evaluation of primary DNA damage extent as evaluated by the comet assay in PBL (biomarker of biologically effective dose). Moreover, in this study we analysed genetic polymorphisms in the CYP1A1 (MspI and Ile/Val sites), EPHX (exons 3 and 4), and GSTM1 genes to evaluate the impact of these metabolic genotypes (biomarker of individual susceptibility) on the levels of urinary 1OHP and the extent of primary DNA damage.
Urinary 1OHP levels were significantly higher in exposed workers than in matched controls. This finding is consistent with the results of other studies aimed to determine the effects of occupational exposure to PAH on urinary 1OHP concentrations [3, 43–45]. The statistically significant increase in urinary 1OHP levels observed in control smokers as compared to control non smokers was not confirmed in the group of exposed workers. This aspect could be explained in terms of saturating dose, probably due to the large work-related effects, at which no further effect can be seen at higher doses.
In the comet assay, DNA strand breakage is quantified from geometric (e.g. migration distance) and fluorescence (e.g. per cent of fluorescence migrated in the tail) measurements both by eye (i.e. comet score) or computerized image analysis. In the presence of DNA damage, the distribution of tail parameters results strongly skewed to the right and the two sides of the distribution have different spreads. In this situation the mean is strongly influenced by extreme observations and, as the resulting large standard deviation, appears to be rather uninformative. In this situation, skewed not symmetrical distributions could be defined more accurately by the median value (i.e. 50th percentile) or the values corresponding to the 75th/95th percentile . This statistical approach in the comet assay was adopted in several researches [47–49].
In epidemiological studies, enrolled subjects should be characterized for a given aspect in a simple manner, preferably with a single number. In this molecular epidemiology approach we have chosen to describe the extent of DNA damage only in terms of tail intensity, a parameter describing the percentage of DNA migrated in the tail which is relatively unaffected by threshold settings in the computerized imaging .
As regard the extent of primary DNA damage, the comet assay has been employed only in few studies aimed to evaluate the PAH genotoxic effects in exposed workers [45, 50–55], with both positive and negative findings. DNA single-strand breakage did not differ between 99 potroom workers at an aluminum reduction plant and 55 unexposed referents . Occupational exposure to PAH did not result in increased DNA strand breaks coke oven workers as compared to unexposed subjects . The analysis of DNA damage did not show significant differences between 42 primary aluminium industry workers and 16 local residents with no occupational exposure to PAH . No effect of occupational exposure was observed in 50 coke oven workers as compared to 50 control workers not exposed to PAH in the extent of DNA damage . Exposure to PAH caused a significantly higher single strand DNA breakage in lymphocytes and granulocytes of 24 workers from automobile emission inspection companies and 28 workers from a waste incinerating company as compared to 43 matched unexposed subjects . The extent of primary DNA damage evaluated with the comet assay was found to be 3.13 times higher for graphite-electrode-producing plant workers (n = 29) when compared with controls (n = 32) . Thus, the findings of the present work are in line with the results of the unique study considering graphite-electrode-producing plant workers  and support the evidence that occupational exposure to PAH during graphite electrode manufacturing can result in primary DNA damage (strand breakage as evaluated with the comet assay).
Smoking habit did not increase any of the DNA damage parameters, either in the exposed than in the reference group. The effect of smoking as a potential confounder in occupational studies has been recently evaluated in a meta-analysis study of the available, conflicting results obtained with the comet assay . The authors concluded that an effect of smoking could not be formally demonstrated when the evaluation of DNA damage was based on image analysis.
Homozygous variant carriers of the CYP1A1 polymorphisms (MspI and Ile/Val) are extremely rare in Caucasians . The frequencies (all subjects) that we found (2.65 and 0.53%, for MspI and Ile/Val, respectively) agree with these observations. The EPHX allele frequencies found in this study are similar to those reported in non-Hispanic whites . The frequency of the GSTM1 null genotype (43.4%) among the studied population (all subjects) agrees with the frequencies reported in the literature for the Caucasian population indicating that 40÷50% of the considered subjects lack of GSTM1 activity [59, 60].
Our results indicate that polymorphisms in EPHX exons 3 and 4 is associated to higher urinary concentrations of 1OHP, whereas the presence of variants in the CYP1A1 gene (i.e. MspI and Ile/Val mutations) as well as the presence of the GSTM1 null genotype showed to have no effect on urinary 1OHP excretion. The increased urinary excretion of 1OHP in subjects having a high EPHX activity agrees with previously published results . Whereas, the absence of a relationship between urinary levels of 1OHP and the presence of the GSTM1 null genotype is not unexpected, as 1OHP is mainly excreted as glucuronide conjugate .
None of the genotypes analyzed (CYP1A1, EPHX, and GSTM1) had any significant influence on primary DNA damage as evaluated by the comet assay. However, it was reported that DNA damage by benzo(a)pyrene (i.e. benzo [a]pyrene diolepoxide-DNA adducts) in PAH-exposed coke oven workers is influenced by smoking habits and GSTM1 polymorphisms [63, 64]. Thus, it could be of interest in the future, in workers exposed to high concentration of PAH, to compare damage caused by benzo(a)pyrene (such as BPDE-DNA adducts) and DNA strand breakage (as evaluated with the comet assay) also in relation to genetic polymorphisms.
No significant correlations were observed between urinary levels of 1OHP and the extent of DNA strand breakage (Spearman's correlation coefficients: r = 0.338, r = 0.157, and r = 0.205, for tail length, tail intensity, and tail moment, respectively) in this study. The absence of significant correlations could be explained in terms of different persistence for the considered biomarkers. In fact, urinary metabolites (i.e. 1OHP) mirror the exposure during the last workshift and some days before , whereas the alkaline comet assay measures temporary strand breaks that happen before the DNA repair systems or cell turnover occur .
The results of a cohort study among workers in a graphite electrode production plant in Italy showed an excess of mortality for cancer in these workers with a standardised mortality ratio of 1.27 (CI95% 1.07–1.50) . The main strategy of primary cancer prevention is minimize exposures to recognized genotoxic/carcinogenic risk factors. The outcomes of the present study, together with the results previously published by Marczynski et al. , show that molecular epidemiology approaches (i.e. cross-sectional studies of genotoxicity biomarkers) can play a role in identifying common genetic risk factors, also attempting to associate the effects with measured exposure data. Moreover, categorization of subjects as normal or outlier, as performed in this study, allowed the evaluation of the association between exposure to genotoxins in this occupational branch and DNA damage, highlighting an increased genotoxic risk with a statistically significant OR = 2.59 (CI95% 1.32–5.05).
The Authors wish to thank Prof. Giuseppina Scassellati-Sforzolini, Emeritus Professor of Hygiene, University of Perugia, for her helpful assistance during this research and constructive comments on the manuscript.
- Boffetta P, Jourenkova N, Gustavsson P: Cancer risk from occupational and environmental exposure to polycyclic aromatic hydrocarbons. Cancer Causes Control. 1997, 8 (3): 444-472. 10.1023/A:1018465507029.View ArticlePubMedGoogle Scholar
- Moulin JJ, Wild P, Mur JM, Lafontaine M, Lefer M, Mercier-Gallay M, Villemot P, Whebi V, Coulon JP: Risk of lung, larynx, pharynx and buccal cavity cancers among carbon electrode manufacturing workers. Scand J Work Environ Health. 1989, 15 (1): 30-37.View ArticlePubMedGoogle Scholar
- Dell'Omo M, Muzi G, Marchionna G, Latini L, Carrieri P, Paolemili P, Abbritti G: Preventive measures reduce exposure to polycyclic aromatic hydrocarbons at a graphite electrode plant. Occup Environ Med. 1998, 55 (6): 401-406.View ArticlePubMedPubMed CentralGoogle Scholar
- Monarca S, Pasquini R, Sforzolini GS, Viola V, Fagioli F: Application of the Salmonella mutagenicity assay and determination of polycyclic aromatic hydrocarbons in workplaces exposed to petroleum pitch and petroleum coke. Int Arch Occup Environ Health. 1982, 49 (3–4): 223-239. 10.1007/BF00377932.View ArticlePubMedGoogle Scholar
- Pasquini R, Monarca S, Sforzolini GS, Conti R, Fagioli F: Mutagens in urine of carbon electrode workers. Int Arch Occup Environ Health. 1982, 50 (4): 387-395. 10.1007/BF00377835.View ArticlePubMedGoogle Scholar
- Petry T, Schmid P, Schlatter C: Airborne exposure to polycyclic aromatic hydrocarbons (PAHs) and urinary excretion of 1-hydroxypyrene of carbon anode plant workers. Ann Occup Hyg. 1996, 40 (3): 345-357.View ArticlePubMedGoogle Scholar
- Teta MJ, Ott MG, Schnatter AR: Population based mortality surveillance in carbon products manufacturing plants. Br JInd Med. 1987, 44 (5): 344-350.Google Scholar
- Mori I: Cancer mortality among man-made graphite electrode manufacturing workers: results of a 38 year follow up. Occup Environ Med. 2002, 59 (7): 473-480. 10.1136/oem.59.7.473.View ArticlePubMedPubMed CentralGoogle Scholar
- Gustavsson P, Bellander T, Johansson L, Salmonsson S: Surveillance of mortality and cancer incidence among Swedish graphite electrode workers. Environ Res. 1995, 70 (1): 7-10. 10.1006/enrs.1995.1039.View ArticlePubMedGoogle Scholar
- Donato F, Monarca S, Marchionna G, Rossi A, Cicioni C, Chiesa R, Colin D, Boffetta P: Mortality from cancer and chronic respiratory diseases among workers who manufacture carbon electrodes. Occup Environ Med. 2000, 57 (7): 484-487. 10.1136/oem.57.7.484.View ArticlePubMedPubMed CentralGoogle Scholar
- IARC: Polynuclear aromatic compounds, Part 1, Chemical, environmental and experimental data. IARC Monograph on the Evaluation of the Carcinogenic Risk of Chemicals to Humans. 1983, Lyon (France): International Agency for Research on Cancer, 32:Google Scholar
- Crofts F, Taioli E, Trachman J, Cosma GN, Currie D, Toniolo P, Garte SJ: Functional significance of different human CYP1A1 genotypes. Carcinogenesis. 1994, 15 (12): 2961-2963. 10.1093/carcin/15.12.2961.View ArticlePubMedGoogle Scholar
- Coles B, Ketterer B: The role of glutathione and glutathione transferases in chemical carcinogenesis. Crit Rev BiochemMol Biol. 1990, 25 (1): 47-70.View ArticleGoogle Scholar
- Ingelman-Sundberg M: Genetic variability in susceptibility and response to toxicants. Toxicol Lett. 2001, 120 (1–3): 259-268. 10.1016/S0378-4274(01)00278-8.View ArticlePubMedGoogle Scholar
- Pavanello S, Clonfero E: Biological indicators of genotoxic risk and metabolic polymorphisms. Mutat Res. 2000, 463 (3): 285-308. 10.1016/S1383-5742(00)00051-X.View ArticlePubMedGoogle Scholar
- Schoket B, Papp G, Levay K, Mrackova G, Kadlubar FF, Vincze I: Impact of metabolic genotypes on levels of biomarkers of genotoxic exposure. Mutat Res. 2001, 482 (1–2): 57-69.View ArticlePubMedGoogle Scholar
- Bartsch H, Nair U, Risch A, Rojas M, Wikman H, Alexandrov K: Genetic polymorphism of CYP genes, alone or in combination, as a risk modifier of tobacco-related cancers. Cancer Epidemiol Biomarkers Prev. 2000, 9 (1): 3-28.PubMedGoogle Scholar
- Pastorelli R, Guanci M, Cerri A, Negri E, La Vecchia C, Fumagalli F, Mezzetti M, Cappelli R, Panigalli T, Fanelli R, et al: Impact of inherited polymorphisms in glutathione S-transferase M1, microsomal epoxide hydrolase, cytochrome P450 enzymes on DNA, and blood protein adducts of benzo(a)pyrene-diolepoxide. Cancer Epidemiol Biomarkers Prev. 1998, 7 (8): 703-709.PubMedGoogle Scholar
- Rojas M, Alexandrov K, Cascorbi I, Brockmoller J, Likhachev A, Pozharisski K, Bouvier G, Auburtin G, Mayer L, Kopp-Schneider A, et al: High benzo[a]pyrene diol-epoxide DNA adduct levels in lung and blood cells from individuals with combined CYP1A1 MspI/Msp-GSTM1*0/*0 genotypes. Pharmacogenetics. 1998, 8 (2): 109-118. 10.1097/00008571-199804000-00003.View ArticlePubMedGoogle Scholar
- Teixeira JP, Gaspar J, Martinho G, Silva S, Rodrigues S, Mayan O, Martin E, Farmer PB, Rueff J: Aromatic DNA adduct levels in coke oven workers: correlation with polymorphisms in genes GSTP1, GSTM1, GSTT1 and CYP1A1. Mutat Res. 2002, 517 (1–2): 147-155.View ArticlePubMedGoogle Scholar
- Au WW, Oh HY, Grady J, Salama SA, Heo MY: Usefulness of genetic susceptibility and biomarkers for evaluation of environmental health risk. Environ Mol Mutagen. 2001, 37 (3): 215-225. 10.1002/em.1030.View ArticlePubMedGoogle Scholar
- Autrup H: Genetic polymorphisms in human xenobiotica metabolizing enzymes as susceptibility factors in toxic response. Mutat Res. 2000, 464 (1): 65-76.View ArticlePubMedGoogle Scholar
- Raaka S, Hassett C, Omiencinski CJ: Human microsomal epoxide hydrolase: 5'-flanking region genetic polymorphisms. Carcinogenesis. 1998, 19 (3): 387-393. 10.1093/carcin/19.3.387.View ArticlePubMedGoogle Scholar
- Board P, Coggan M, Johnston P, Ross V, Suzuki T, Webb G: Genetic heterogeneity of the human glutathione transferases: a complex of gene families. Pharmacol Ther. 1990, 48 (3): 357-369. 10.1016/0163-7258(90)90054-6.View ArticlePubMedGoogle Scholar
- Binkova B, Topinka J, Mrackova G, Gajdosova D, Vidova P, Stavkova Z, Peterka V, Pilcik T, Rimar V, Dobias L, et al: Coke oven workers study: the effect of exposure and GSTM1 and NAT2 genotypes on DNA adduct levels in white blood cells and lymphocytes as determined by 32P-postlabelling. Mutat Res. 1998, 416 (1–2): 67-84.View ArticlePubMedGoogle Scholar
- Alexandrie AK, Warholm M, Carstensen U, Axmon A, Hagmar L, Levin JO, Ostman C, Rannug A: CYP1A1 and GSTM1 polymorphisms affect urinary 1-hydroxypyrene levels after PAH exposure. Carcinogenesis. 2000, 21 (4): 669-676. 10.1093/carcin/21.4.669.View ArticlePubMedGoogle Scholar
- Sram RJ, Binkova B: Molecular epidemiology studies on occupational and environmental exposure to mutagens and carcinogens, 1997–1999. Environ Health Perspect. 2000, 108 (Suppl 1): 57-70. 10.2307/3454632.View ArticlePubMedPubMed CentralGoogle Scholar
- Jongeneelen FJ, Anzion RB, Henderson PT: Determination of hydroxylated metabolites of polycyclic aromatic hydrocarbons in urine. J Chromatogr. 1987, 413: 227-232. 10.1016/0378-4347(87)80230-X.View ArticlePubMedGoogle Scholar
- VanRooij JG, Bodelier-Bade MM, Jongeneelen FJ: Estimation of individual dermal and respiratory uptake of polycyclic aromatic hydrocarbons in 12 coke oven workers. Br J Ind Med. 1993, 50 (7): 623-632.PubMedPubMed CentralGoogle Scholar
- Tice RR, Agurell E, Anderson D, Burlinson B, Hartmann A, Kobayashi H, Miyamae Y, Rojas E, Ryu JC, Sasaki YF: Single cell gel/comet assay: guidelines for in vitro and in vivo genetic toxicology testing. Environ Mol Mutagen. 2000, 35 (3): 206-221. 10.1002/(SICI)1098-2280(2000)35:3<206::AID-EM8>3.0.CO;2-J.View ArticlePubMedGoogle Scholar
- Pastorelli R, Cerri A, Rozio M, Dell'Omo M, Muzi G, Abbritti G, Fanelli R, Airoldi L: Benzo (a) pyrene diolepoxide adducts to albumin in workers exposed to polycyclic aromatic hydrocarbons: association with specific CYP1A1, GSTM1, GSTP1 and EHPX genotypes. Biomarkers. 2001, 6 (5): 357-374. 10.1080/13547500110044267.View ArticlePubMedGoogle Scholar
- Hayashi S, Watanabe J, Kawajiri K: Genetic polymorphisms in the 5'-flanking region change transcriptional regulation of the human cytochrome P450IIE1 gene. J Biochem (Tokyo). 1991, 110 (4): 559-565.Google Scholar
- Cascorbi I, Brockmoller J, Roots I: A C4887A polymorphism in exon 7 of human CYP1A1: population frequency, mutation linkages, and impact on lung cancer susceptibility. Cancer Res. 1996, 56 (21): 4965-4969.PubMedGoogle Scholar
- Bell DA, Taylor JA, Paulson DF, Robertson CN, Mohler JL, Lucier GW: Genetic risk and carcinogen exposure: a common inherited defect of the carcinogen-metabolism gene glutathione S-transferase M1 (GSTM1) that increases susceptibility to bladder cancer. J Natl Cancer Inst. 1993, 85 (14): 1159-1164. 10.1093/jnci/85.14.1159.View ArticlePubMedGoogle Scholar
- Brendler-Schwaab SY, Schmezer P, Liegibel U, Weber S, Michalek K, Tompa A, Pool-Zobel BL: Cells of different tissues for in vitro and in vivo studies in toxicology: Compilation of isolation methods. Toxicology In Vitro. 1994, 8 (6): 1285-1302. 10.1016/0887-2333(94)90122-8.View ArticlePubMedGoogle Scholar
- Jones KH, Senft JA: An improved method to determine cell viability by simultaneous staining with fluorescein diacetate-propidium iodide. J Histochem Cytochem. 1985, 33 (1): 77-79.View ArticlePubMedGoogle Scholar
- Singh NP, McCoy MT, Tice RR, Schneider EL: A simple technique for quantitation of low levels of DNA damage in individual cells. Exp Cell Res. 1988, 175 (1): 184-191. 10.1016/0014-4827(88)90265-0.View ArticlePubMedGoogle Scholar
- Villarini M, Moretti M, Pasquini R, Scassellati-Sforzolini G, Fatigoni C, Marcarelli M, Monarca S, Rodriguez AV: In vitro genotoxic effects of the insecticide deltamethrin in human peripheral blood leukocytes: DNA damage ('comet' assay) in relation to the induction of sister-chromatid exchanges and micronuclei. Toxicology. 1998, 130 (2–3): 129-139. 10.1016/S0300-483X(98)00097-3.View ArticlePubMedGoogle Scholar
- Villarini M, Scassellati-Sforzolini G, Moretti M, Pasquini R: In vitro genotoxicity of terbutryn evaluated by the alkaline single-cell microgel-electrophoresis "comet" assay. Cell Biol Toxicol. 2000, 16 (5): 285-292. 10.1023/A:1026794213308.View ArticlePubMedGoogle Scholar
- Moretti M, Villarini M, Scassellati-Sforzolini G, Monarca S, Salucci A, Vicent Rodriguez A: Application of the single-cell gel-electrophoresis ('comet') assay to the detection of primary DNA damage in workers of the rubber industry: comparison of manual and computerized analysis. Toxicological and Environmental Chemistry. 1999, 13-24.Google Scholar
- Bonassi S, Fontana V, Ceppi M, Barale R, Biggeri A: Analysis of correlated data in human biomonitoring studies. The case of high sister chromatid exchange frequency cells. Mutat Res. 1999, 438 (1): 13-21.View ArticlePubMedGoogle Scholar
- Collins AR: The comet assay for DNA damage and repair: principles, applications, and limitations. Mol Biotechnol. 2004, 26 (3): 249-261. 10.1385/MB:26:3:249.View ArticlePubMedGoogle Scholar
- Nan HM, Kim H, Lim HS, Choi JK, Kawamoto T, Kang JW, Lee CH, Kim YD, Kwon EH: Effects of occupation, lifestyle and genetic polymorphisms of CYP1A1, CYP2E1, GSTM1 and GSTT1 on urinary 1-hydroxypyrene and 2-naphthol concentrations. Carcinogenesis. 2001, 22 (5): 787-793. 10.1093/carcin/22.5.787.View ArticlePubMedGoogle Scholar
- Ovrebo S, Haugen A, Fjeldstad PE, Hemminki K, Szyfter K: Biological monitoring of exposure to polycyclic aromatic hydrocarbon in an electrode paste plant. J Occup Med. 1994, 36 (3): 303-310. 10.1097/00043764-199403000-00007.View ArticlePubMedGoogle Scholar
- Siwinska E, Mielzynska D, Kapka L: Association between urinary 1-hydroxypyrene and genotoxic effects in coke oven workers. Occup Environ Med. 2004, 61 (3): e10-10.1136/oem.2002.006643.View ArticlePubMedPubMed CentralGoogle Scholar
- Duez P, Dehon G, Kumps A, Dubois J: Statistics of the Comet assay: a key to discriminate between genotoxic effects. Mutagenesis. 2003, 18 (2): 159-166. 10.1093/mutage/18.2.159.View ArticlePubMedGoogle Scholar
- Buschini A, De Palma G, Poli P, Martino A, Rossi C, Mozzoni P, Scotti E, Buzio L, Bergamaschi E, Mutti A: Genetic polymorphism of drug-metabolizing enzymes and styrene-induced DNA damage. Environ Mol Mutagen. 2003, 41 (4): 243-252. 10.1002/em.10150.View ArticlePubMedGoogle Scholar
- Lu Y, Morimoto K, Nakayama K: Health practices and leukocyte DNA damage in Japanese hard-metal workers. Prev Med. 2006, 43 (2): 140-144. 10.1016/j.ypmed.2006.03.008.View ArticlePubMedGoogle Scholar
- Sigurdson AJ, Hauptmann M, Alexander BH, Doody MM, Thomas CB, Struewing JP, Jones IM: DNA damage among thyroid cancer and multiple cancer cases, controls, and long-lived individuals. Mutat Res. 2005, 586 (2): 173-188.View ArticlePubMedGoogle Scholar
- Carstensen U, Hou SM, Alexandrie AK, Hogstedt B, Tagesson C, Warholm M, Rannug A, Lambert B, Axmon A, Hagmar L: Influence of genetic polymorphisms of biotransformation enzymes on gene mutations, strand breaks of deoxyribonucleic acid, and micronuclei in mononuclear blood cells and urinary 8-hydroxydeoxyguanosine in potroom workers exposed to polyaromatic hydrocarbons. Scand J Work Environ Health. 1999, 25 (4): 351-360.View ArticlePubMedGoogle Scholar
- Crebelli R, Carta P, Andreoli C, Aru G, Dobrowolny G, Rossi S, Zijno A: Biomonitoring of primary aluminium industry workers: detection of micronuclei and repairable DNA lesions by alkaline SCGE. Mutat Res. 2002, 516 (1–2): 63-70.View ArticlePubMedGoogle Scholar
- Marczynski B, Rihs HP, Rossbach B, Holzer J, Angerer J, Scherenberg M, Hoffmann G, Bruning T, Wilhelm M: Analysis of 8-oxo-7,8-dihydro-2'-deoxyguanosine and DNA strand breaks in white blood cells of occupationally exposed workers: comparison with ambient monitoring, urinary metabolites and enzyme polymorphisms. Carcinogenesis. 2002, 23 (2): 273-281. 10.1093/carcin/23.2.273.View ArticlePubMedGoogle Scholar
- Moretti M, Villarini M, Scassellati-Sforzolini G, Monarca S, Libraro M, Fatigoni C, Donato F, Leonardis C, Perego L: Biological monitoring of genotoxic hazard in workers of the rubber industry. Environ Health Perspect. 1996, 104 (Suppl 3): 543-545. 10.2307/3432820.View ArticlePubMedPubMed CentralGoogle Scholar
- Sul D, Oh E, Im H, Yang M, Kim CW, Lee E: DNA damage in T- and B-lymphocytes and granulocytes in emission inspection and incineration workers exposed to polycyclic aromatic hydrocarbons. Mutat Res. 2003, 538 (1–2): 109-119.View ArticlePubMedGoogle Scholar
- van Delft JH, Steenwinkel MS, van Asten JG, de Vogel N, Bruijntjes-Rozier TC, Schouten T, Cramers P, Maas L, van Herwijnen MH, van Schooten F, et al: Biological monitoring the exposure to polycyclic aromatic hydrocarbons of coke oven workers in relation to smoking and genetic polymorphisms for GSTM1 and GSTT1. Ann Occup Hyg. 2001, 45 (5): 395-408.View ArticlePubMedGoogle Scholar
- Hoffmann H, Hogel J, Speit G: The effect of smoking on DNA effects in the comet assay: a meta-analysis. Mutagenesis. 2005, 20 (6): 455-466. 10.1093/mutage/gei064.View ArticlePubMedGoogle Scholar
- Cosma G, Crofts F, Currie D, Wirgin I, Toniolo P, Garte SJ: Racial differences in restriction fragment length polymorphisms and messenger RNA inducibility of the human CYP1A1 gene. Cancer Epidemiol Biomarkers Prev. 1993, 2 (1): 53-57.PubMedGoogle Scholar
- Zhao H, Spitz MR, Gwyn KM, Wu X: Microsomal epoxide hydrolase polymorphisms and lung cancer risk in non-Hispanic whites. Mol Carcinog. 2002, 33 (2): 99-104. 10.1002/mc.10023.View ArticlePubMedGoogle Scholar
- Hirvonen A, Husgafvel-Pursiainen K, Anttila S, Vainio H: The GSTM1 null genotype as a potential risk modifier for squamous cell carcinoma of the lung. Carcinogenesis. 1993, 14 (7): 1479-1481. 10.1093/carcin/14.7.1479.View ArticlePubMedGoogle Scholar
- Somorovska M, Jahnova E, Tulinska J, Zamecnikova M, Sarmanova J, Terenova A, Vodickova L, Liskova A, Vallova B, Soucek P, et al: Biomonitoring of occupational exposure to styrene in a plastics lamination plant. Mutat Res. 1999, 428 (1–2): 255-269.View ArticlePubMedGoogle Scholar
- Kuljukka-Rabb T, Nylund L, Vaaranrinta R, Savela K, Mutanen P, Veidebaum T, Sorsa M, Rannug A, Peltonen K: The effect of relevant genotypes on PAH exposure-related biomarkers. J Expo Anal Environ Epidemiol. 2002, 12 (1): 81-91. 10.1038/sj.jea.7500204.View ArticlePubMedGoogle Scholar
- Apostoli P, Neri G, Lucas D, Manno M, Berthou F: Influence of genetic polymorphisms of CYP1A1 and GSTM1 on the urinary levels of 1-hydroxypyrene. Toxicol Lett. 2003, 144 (1): 27-34. 10.1016/S0378-4274(03)00227-3.View ArticlePubMedGoogle Scholar
- Pavanello S, Pulliero A, Siwinska E, Mielzynska D, Clonfero E: Reduced nucleotide excision repair and GSTM1-null genotypes influence anti-B[a]PDE-DNA adduct levels in mononuclear white blood cells of highly PAH-exposed coke oven workers. Carcinogenesis. 2005, 26 (1): 169-175. 10.1093/carcin/bgh303.View ArticlePubMedGoogle Scholar
- Rojas M, Cascorbi I, Alexandrov K, Kriek E, Auburtin G, Mayer L, Kopp-Schneider A, Roots I, Bartsch H: Modulation of benzo[a]pyrene diolepoxide-DNA adduct levels in human white blood cells by CYP1A1, GSTM1 and GSTT1 polymorphism. Carcinogenesis. 2000, 21 (1): 35-41. 10.1093/carcin/21.1.35.View ArticlePubMedGoogle Scholar
- Cotelle S, Ferard JF: Comet assay in genetic ecotoxicology: a review. Environ Mol Mutagen. 1999, 34 (4): 246-255. 10.1002/(SICI)1098-2280(1999)34:4<246::AID-EM4>3.0.CO;2-V.View ArticlePubMedGoogle Scholar
- Merlo DF, Garattini S, Gelatti U, Simonati C, Covolo L, Ceppi M, Donato F: A mortality cohort study among workers in a graphite electrode production plant in Italy. Occup Environ Med. 2004, 61 (2): e9-10.1136/oem.2003.009357.View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2458/7/270/prepub