Printer Friendly

Inhalable metal-rich air particles and histone H3K4 dimethylation and H3K9 acetylation in a cross-sectional study of steel workers.

BACKGROUND: Epidemiology investigations have linked exposure to ambient and occupational air particulate matter (PM) with increased risk of lung cancer. PM contains carcinogenic and toxic metals, including arsenic and nickel, which have been shown in in vitro studies to induce histone modifications that activate gene expression by inducing open-chromatin states. Whether inhalation of metal components of PM induces histone modifications in human subjects is undetermined.

OBJECTIVES: We investigated whether the metal components of PM determined activating histone modifications in 63 steel workers with well-characterized exposure to metal-rich PM.

METHODS: We determined histone 3 lysine 4 dimethylation (H3K4me2) and histone 3 lysine 9 acetylation (H3 K9ac) on histories from blood leukocytes. Exposure to inhalable metal components (aluminum, manganese, nickel, zinc, arsenic, lead, iron) and to total PM was estimated for each study subject.

RESULTS: Both H3K4me2 and H3K9ac increased in association with years of employment in the plant (p-trend = 0.04 and 0. [beta]6, respectively). H3K4me2 increased in association with air levels of nickel [p = 0.16; 95% confidence interval (CI), 0.03-0.3], arsenic [beta] = 0.16; 95% CI, 0.02-0.3), and iron [beta] = 0.14; 95% CI, 0.01-0.26). H3K9ac showed nonsignificant positive associations with air levels of nickel [beta] = 0.24; 95% CI, -0.02 to 0.51), arsenic [beta] = 0.21; 95% CI, -0.06 to 0.48), and iron ([beta] = 0.22; 95% CI, -0.03 to 0.47). Cumulative exposures to nickel and arsenic, defined as the product of years of employment by metal air levels, were positively correlated with both H3K4me2 (nickel: [beta] = 0.16; 95% CI, 0.01-0.3; arsenic: [beta] = 0.16; 95% CI, 0.03-0.29) and H3K9ac (nickel: [beta] = 0.27; 95% CI, 0.01-0.54; arsenic: [beta] = 0.28; 95% CI, 0.04-0.51).

CONCLUSIONS: Our results indicate histone modifications as a novel epigenetic mechanism induced in human subjects by long-term exposure to inhalable nickel and arsenic.

KEYWORDS: environmental carcinogens, epigenetics, histone modifications, metals, particulate matter. Environ Health Perspect 119:964-969 (2011). doi:10.1289/ehp. 1002955 [Online 8 March 2011]

Ambient and occupational exposure to inhalable particulate matter (PM) has been associated with increased risk of lung cancer (Dockery et al. 1993; Gibb et al. 2000; Kuo et al. 1999). Epidemiologic and in vivo studies suggest that the metal components of PM may be responsible for PM health effects, including lung cancer (Chang et al. 2005; Chen and Hwang 2005; Conroy et al. 2008; Corey et al. 2006; Coyle et al. 2006; Franklin et al. 2008; MacNee and Donaldson 2003; Roller 2009; Wild et al. 2009). Despite the well-recognized carcinogenic potentials of several toxic metals, the molecular mechanisms underlying their associations with cancer risk remain poorly understood. In particular, most carcinogenic metals are weak mutagens and do not induce DNA adduct formation, a key initiating event caused by other carcinogens (Salnikow and Zhitkovich 2008).

Growing evidence indicates that epigenetic dysregulation of gene expression plays a primary role in cancer etiology (Feinberg and Tycko 2004; Ke et al. 2006). Several toxic metals have been shown to bind more avidly to histone proteins than to other biopolymers such as DNA or RNA (Conroy et al. 2008; Zoroddu et al. 2000, 2002). Recent in vitro studies have shown that carcinogenic metals cause posttranslational epigenetic modifications of histone proteins, thus derailing the normal programming of gene expression (Ke et al. 2006; Yan and Boyd 2006; Yan et al. 2003; Zhou X et al. 2009). In lung epithelial cell lines and malignant transformation models, arsenic (Jensen et al. 2009; Zhou X et al. 2009), nickel (Zhou QX et al. 2009), and chromium (Sun et al. 2009; Zhou X et al. 2009) have been linked with specific activating histone modifications, such as H3K4 (histone 3 lysine 4) dimethylation (H3K4me2), that contribute to the formation of a relaxed or "open" chromatin structure permissive for gene transcription (Surralles et al. 1998; Yan et al. 2003). Metal-related induction of activating histone modifications have been suggested to contribute to metal carcinogenesis by causing the expression of tumor suppressor or other cancer-promoting genes (Wang et al. 2005; Zhou X et al. 2008). Inhaled PM pollutants have been previously shown to produce systemic changes in gene expression, which can be detected in peripheral blood of exposed individuals (Salnikow and Zhitkovich 2008). However, whether metals in inhalable PM induce alterations of histone modifications in human subjects has never been evaluated.

Although wide strata of the general population are exposed to low, background levels of toxic and carcinogenic inhalable compounds, industrial workers have been often exposed to significantly higher concentrations of potentially carcinogenic agents released from multiple sources (MacArthur et al. 2009). Foundry work is a specific condition of exposure to inhalable metal-rich PM that has been associated with increased risk of lung cancer in several early investigations (International Agency for Research on Cancer 1987). Even in modern foundry facilities, PM exhibits levels that are well above the concentrations found in ambient outdoor air and also have a larger proportion of toxic metal components (Alley et al. 2009; Fang et al. 2009; Roy et al. 2009; Vijay Bhaskar et al. 2009). In the present investigation in a steel furnace plant, we investigated whether metals in inhalable PM determined changes in H3K4me2 and histone 3 lysine 9 acetylation (H3K9ac), two histone modifications associated with open chromatin states, in peripheral blood leukocytes from workers with well-characterized exposure.

Address correspondence to A. Baccarelli, Department; of Environmental Health, Exposure, Epidemiology and Risk Program, Harvard School of Public Health, P.O. Box 15698, Boston, MA 02215 USA. Telephone: (617) 384-8742. Fax: (617) 384-8745. E-mail: abaccare@hsph.harvard.edu

This work was supported by Cassa di Risparmio delle Provincie Lombarde (CARIPLO) Foundation grant 2007-5469 [Epigenetic Studies of Toxic Hazards in Environmental Research (ESTHER)]; Italian Ministry of Scientific Research, Research Programs of National Interest (gram PRIN 2007-2S2HT8); and new investigator funding from the Harvard School of Public Health-National Institute of Environmental Health Sciences Center for Environmental Health (ES000002).

The authors declare they have no actual or potential competing financial interests.

Received 9 September 2010; accepted 8 March 2011.

Materials and Methods

Study subjects. We recruited 63 male healthy workers, free of cancer and cardiopulmonary disease, who had been working in a steel production plant in Brescia, Italy, for at least 1 year. Individual written informed consent and approval from the local institutional review board were obtained before the study. All of the study subjects had a rotating weekly schedule based on 4 consecutive working days of 8 hr each, followed by 2 days of rest, The study subjects worked in 11 different areas within the plant, which were selected to provide a wide contrast of exposures between the study subjects. The exposure of each of the study subjects in the plant was monitored for the first 3 working days of a work week. On the fourth day, each subject donated a 20-mL peripheral blood sample. We used ethylene diamine tetraacetic acid (EDTA) tubes to collect 7 mL whole blood that was promptly (within 30 min from the blood drawing) centrifuged on site at 2,500 rpm for 15 min. The bully coat (400 [micro] L) was separated and transferred in a cryovial, immediately snap frozen in vapor phase of liquid nitrogen, shipped in nitrogen dry shippers to the laboratory, and kept stored in vapor phase of liquid nitrogen until histone extraction. A self-administered questionnaire was used to collect detailed information on lifestyle, drug use, medical conditions, body mass index (BMI), education, and residential history. Records from the factory administrative files were used to extract information on occupational history.

Exposure assessment. PM metal components (aluminum, manganese, nickel, zinc, arsenic, lead, iron) and PM mass [PM with aerodynamic diameters [less than or equal to] 10 um and [less than or equal to] 1 [micro] m ([PM.sub.10] and [PM.sub.1]], respectively)] were measured in each of the 11 work areas of the steel production plant. We measured air concentrations of individual metal PM components in [PM.sub.10], through multielement analysis performed by means of inductively coupled-plasma mass spectrometer (ELAN DRC 11; PerkinElmer, Waltham, MA, USA) using the total quant method. External calibration was performed using calibration standard 3, stock multielement (10 [micro] g/mL; PerkinElmer). [PM.sub.10], and [PM.sub.1] were measured using a Grimm 1100 light-scattering dust analyzer (Grimm Technologies, Inc., Douglasville, GA, USA).

Study subjects recorded in a personal log the time they spent in each of the work areas. Personal exposure was calculated as the average of work area levels weighted by the time spent in each area. All metal and PM levels were expressed in micrograms per cubic meter, in the plant, exposure levels have shown very little variability over time, as measures repeated over 3 years in a subset of the study population showed very high correlations ([r.sup.2] Greater than 0.90). Therefore, the time-weighted levels of metals and PM represented, in addition to the exposure during the week of the study, also a measure of the usual exposure of the study subjects (Tarantini et al. 2009). We estimated cumulative exposures as the product of the time-weighted levels of metals and PM during the study by the years of employment in the plant.

Total histone extraction and histone modification analysis. We obtained buffy coat from peripheral blood collected in EDTA tubes centrifuged at room temperature (22-25 [degrees] C) for 15 min at 1,500 x g. Red blood cell lysis solution (s.r.l. cod.A7933; Promega, Madison, WI, USA) was added to the buffy coat to wash out red blood cells. After 10 min at room temperature, the mixture was centrifuged at 2,500 x g for 15 min, and the supernatant was discarded. Remaining monolayer cells were processed according to protocol used by Chen ec al. (2006). Briefly, cells were lysed in 1 mL ice-cold radioimmunoprecipitation assay buffer (Santa Cruz Biotechnology, Santa Cruz, CA, USA) supplemented with a protease inhibitor mixture (Roche Applied Sciences, Indianapolis, IN, USA) for 10 min. The sample was then collected and centrifuged at 10,000 x g for 10 min. After discarding supernatant, remaining pellet was resuspended in 0.4 N [H.sub.2] [SO.sub.4]. After incubation on ice for 90 min, the sample was centrifuged at 14,000 x g for 15 min. The supernatant was mixed with cold acetone and kept at--20 [degrees] C overnight. The histones were collected by centrifugation at 14,000 x g for 15 min. After one wash with acetone, the histones were air dried and resuspended in 500 [micro] L water. We measured total proteins in each sample by means of the Bradford assay according to manufacturer's instructions (protein assay kit 500-0002; Bio-Rad Laboratories, Milan, Italy). We used equal amounts of proteins (4 [micro] g) to normalize results of the subsequent analysis on histones.

We used a solid-phase sandwich enzyme-linked immunosorbent assay (ELISA), using monoclonal antibodies, to detect endogenous levels of H3K4me2 and H3K9ac (PathScan Sandwich ELISA Kits; Cell Signaling Technology, Beverly, MA, USA) according to the manufacturer's protocol. The assays we used in our study are the PathScan dimethyL-histone H3 (Lys4) Sandwich ELISA Kit 7124 and PathScan acetyl-histone H3 (Lys9) Sandwich ELISA Kit 7121. The assays use dimethyl-histone H3 (Lys4) (C64G9) rabbit monoclonal antibodv 9725 and acetyl-histone H3 (Lys9) antibody 9671, respectively, which have been shown by the manufacturer to be specific for the histone modifications of concern. Sample measurements were performed in duplicate. We used a Synergy HT-BioTek spectrophotometer (Winooski, VT, USA) to read 450 nm absorbance. The absorbance values at 450 nm directly reflected the concentration of modified histones (Deligezer et al. 2010). According to the Beer-Lambert law, optical density (OD; absorbance) is used for colorimetric analysis so that readings relate directly to concentration. The coefficient of variation in replicate samples of the assays was 0.30% for H3K4me2 and 0.42% for H3K9ac.

Statistical analysis. We performed trend tests using the continuous variable in the regression model and presenting the corresponding two-sided p-values. We evaluated the association of PM metal components and PM mass levels with histone modifications using simple linear regression models, as well as multivariable models adjusting for age, BMI, pack-years, and percent granulocytes in the differential blood count as continuous variables and education as categorical variable. The independent variables used in muitivariable models were selected a priori and included general characteristics potentially associated with cancer risks or other carcinogenic exposures. In addition, we adjusted for percent granulocytes to account for possible shifts in the proportion of leukocytes subtypes associated with the exposures. As sensitivity analyses, we fitted in the models as independent variables data from differential white blood counts (i.e., percent lymphocytes, monocytes, eosinophils, or basophils) or duration of sample storage and found no major changes in the results. Regression diagnostics were performed separately for each model. We examined whether the exposure--response relationships were linear through graphical inspection. Furthermore, we fitted a polynomial regression by including a quadratic term for exposure and compared these models with the linear model using the likelihood ratio test. Neither graphical inspection nor the likelihood ratio tests suggested any departure from linearity.

Outliers were excluded from regression analysis by dropping observations with studentized residuals that exceeded Plus 3 or Hyphen or Minus 3. The number of outliers removed ranged from a minimum of 0 to a maximum of 3.

Regression coefficients were computed with ordinary least squares estimators. To compare the magnitude of the associations of H3K4mc2 and H3K9ac with different exposures, we calculated standardized regression coefficients and 95% confidence intervals (CIs) that express the change in histone modifications associated with an increase in exposure equal to the difference between the 90th and 10th percentile of the exposure distribution. We checked regression assumptions by performing diagnostic tests for each model, which included the Shapiro-Wilk test to verify normality of residuals and the White test to verify the homogeneity of variance of the residuals. A two-sided p - value Less than 0.05 was considered statistically significant. Statistical analyses were performed in SAS (version 9.1.3; SAS Institute Inc., Gary, NC, USA) and R (R Foundation for Statistical Computing Vienna, Austria).

Results

Subjects' characteristics and exposure levels. The mean age of the study subjects was 44.00 years (range, 27-55 years). Twenty-five subjects (40%) were current smokers, who reported a median number of 15 cigarettes smoked every' day (range, 0-30 cigarettes/day). The BMI of the study participants showed a median of 26 kg/ [m.sup.2] (range, 20-33 kg/ [m.sup.2]). Table 1 shows the average levels of inhalable air metal concentrations and PM mass estimated. For both metal levels and PM mass, the study subjects showed wide ranges of exposures. The subject with the maximum individual exposure level was at least 17 times more exposed than the subject with the minimum exposure level (e.g., lor PM10: maximum, 1220.17 [micro] g/ [m sup.3]; minimum, 73.72 [micro]g/[m sup.3]). For some of the exposures (aluminum, manganese, zinc, lead), the maximum individual exposure was Greater than 200 times higher than the minimum individual exposure (Table 1).
Table 1. Levels of personal exposure to metal components of PM
and total PM mass.

Exposure (a)           n      Mean [+ or -] SD  Minimum    25th
([micro]g/[m.sup.3])

Aluminum              63    8.50 [+or -] 18.07     0.40    1.48

Manganese             63   11.26 [+or -] 30.41     0.11    1.20

Nickel                63     0.30 [+or -] 0.18     0.02    0.23

Zinc                  63  18.85 [+ or -] 26.37     0.26    1.47

Arsenic               63     0.10 [+ or -] 0.1     0.01    0.02

Lead                  63    7.53 [+ or -]17.46     0.13    0.63

Iron                  63  32.02 [+ or -] 22.08     0.96   18.00

[PM.sub.10]           63  233.42 [+ or-]214.56    73.72  152.23

[PM.sub.1]            63     8.48 [+ or-] 6.18     1.71    3.51

                         Percentile

Exposure (a)                50th    75th  Maximum
([micro]g/[m.sup.3])

Aluminum                    2.05    7.41    84.07

Manganese                   4.63   10.77   174.79

Nickel                      0.25    0.46     0.72

Zinc                        8.45   32.28   129.06

Arsenic                     0.07    0.17     0.31

Lead                        2.87    9.52    99.90

Iron                       25.64   48.69    88.43

[PM.sub.10]               179.45  222.86  1220.17

[PM.sub.1]                  9.01   11.35    30.49

(a)Metal components were measured on the [PM.sub.10] fraction of PM
mass. Coarse PM levels were calculated as the difference between
[PM.sub.10] and PM|.


The correlation matrix between exposure levels showed correlations of different strengths between the exposures (Table 2). We found high correlations [Pearson's correlation coefficient (r) = 0.7-1] for aluminum with manganese, lead, and [PM.sub.10]; for manganese with lead and PM1o; for nickel with arsenic and iron; and between arsenic and iron, [PM.sub.10] and lead, and [PM.sub.1] and [PM.sub.10] Correlations were moderate (r = 0.4-0.7) for [PM.sub.1] with aluminum, manganese, and lead, between nickel and aluminum, and between iron and zinc. Except for PM1 and arsenic (r = -0.20), all other correlations were positive with r Less than 0.40.
Table 2. Matrix of correlations [(r).sup.a] among individual exposures
(PM metal components and PM mass measures).

Exposure              Aluminum     Manganese       Nickel     Zinc

Aluminum                     1

Manganese    0.75     < 0.0001             1

Nickel             0.46 0.0001   0.39 0.0017             1

Zinc               0.18 0.1625   0.21 0.1034   0.32 0.0105       1

Arsenic                   0.15          0.31          0.84    0.25

                        0.2378        0.0143        0.0000  0.0483

Lead                      0.75          0.99          0.38    0.35

                      < 0.0001      < 0.0001        0.0018  0.0047

Iron                      0.15          0.28          0.76    0.48

                        0.2451        0.0261      < 0.0001  0.0001

[PM.sub.10]               0.81          0.82          0.34    0.34

                      < 0.0001      < 0.0001        0.0056  0.0070

[PM.sub.1]                0.63          0.60          0.22    0.28

                      < 0.0001      < 0.0001        0.0897  0.0277

Exposure       Arsenic         Lead       Iron   [PM.sub.10]  [PM.sub.1]

Aluminum

Manganese

Nickel

Zinc

Arsenic                 1

Lead                 0.28             1

                   0.0265

Iron                 0.70          0.32       1

                 < 0.0001        0.0113

[PM.sub.10]          0.04          0.85    0.27             1

                   0.7446      < 0.0001  0.0294

[PM.sub.1]          -0.20          0.64    0.26          0.90         1

                   0.1195      < 0.0001  0.0408      < 0.0001

(a)Pearson's product-moment correlation coefficient and corresponding
p-values.


Associations of histone modifications with subjects' characteristics and years of employment. The range of values of the histone modifications was 0.26-1.20 OD for H3K4me2 (median, 0.99 OD) and -0.50 to 1.06 OD for H3K9ac (median, 0.37 OD).

H3K4me2 was moderately but significantly correlated with H3K9ac (Pearson's r = 0.45, p Less than 0.001). H3K4me2 and H3K9ac were not associated with age, BM, smoking, number of cigarettes/days, smoking duration, pack-years, area of residence (city center, suburbs, rural), self-reported traffic intensity near home, or percent granulocytes (Table 3). H3K9ac was positively associated with education (p = 0.04; Table 3). Both H3K4me2 and H3K9ac increased in association with the years of employment of the study subjects in the steel plant (Table 3). We confirmed the association between years of employment and histone modifications in age-adjusted regression models.
Table 3. Association of the subjects' characteristics with H3K4me2and
H3K9ac (in OD units) measured on the fourth day of a workweek.

                                        H3K4me2

Variable                    n           Mean OD      p-Value (a)
                                       (95% CI)
Aqe (years)
<39                        22  0.90 (0.79-1.02)
39-47                      20  1.00 (0.96-1.05)
>47                        21  0.94 (0.88-1.00)             0.47
BMI (kg/[m.sup.2])
<25                        21   0.96 (0.9-1.02)
25-27.5                    21  0.89 (0.78-1.01)
>27.5                      21  0.99 (0.95-1.03)             0.57
Smoking

Nonsmoker                  24  0.96 (0.88-1.04)
Former smoker              14  0.87 (0.76-0.97)
Current smoker             25  0.98 (0.92-1.05)           0.14 *
No. of cigarettes/day
0                          38  0.93 (0.87-0.99)
1-10                       10  0.98 (0.92-1.05)
>10                            0.98 (0.87-1.09)             0.29
Smoking
duration (years)
0                          23  0.97 (0.88-1.05)
0-19                       19   0.90 (0.8-1.00)
>19                        21  0.97 (0.91-1.03)             0.94
Pack-years
0                          24  0.96 (0.88-1.04)
0-4                        17  0.92 (0.84-0.99)
>4                         20  0.95 (0.86-1.04)             0.79

                                 H3K9ac

Variable                        Mean OD      p-Value (a)
                               (95% CI)
Aqe (years)
<39                    0.41 (0.22-0.59)
39-47                  0.43 (0.26-0.59)
>47                    0.43 (0.27-0.58)             0.88
BMI (kg/[m.sup.2])
<25                    0.48 (0.33-0.62)
25-27.5                0.20 (0.04-0.35)
>27.5                  0.59 (0.45-0.74)             0.36
Smoking

Nonsmoker              0.44 (0.28-0.61)
Former smoker          0.35 (0.12-0.59)
Current smoker         0.43 (0.30-0.57)             0.76
No. of cigarettes/day
0                      0.41 (0.28-0.54)
1-10                   0.39 (0.23-0.55)
>10                    0.46 (0.24-0.68)             0.69
Smoking
duration (years)
0                       0.45(0.28-0.62)
0-19                    0.35(0.15-0.54)
>19                    0.45 (0.32-0.58)             0.98
Pack-years
0                      0.44 (0.28-0.61)
0-4                    0.39 (0.22-0.57)
<4                     0.43 (0.25-0.61)             0.90

                                                           H3K4me2

Variable                            Variable   n           Mean OD
                                                          (95% CI)
Aqe (years)                        Education
<39                           Primary school  12  0.96 (0.83-1.09)
39-47                          Middle school  37   0.94(0.89-0.99)
>47                              High school  14   0.96 0.81-1.10)
BMI (kg/[m.sup.2])         Area of residence
<25                              City center   8   0.89(0.54-1.23)
25-27.5                              Suburbs  41   0.97(0.92-1.01)
>27.5                                  Rural  12   0.91(0.79-1.03)
Smoking                Self-reported traffic
                         intensity near home
Nonsmoker                               High   5   0.86(0.43-1.29)
Former smoker                         Medium  38   0.94(0.89-1.00)
Current smoker                           Low  18   0.97(0.91-1.03)
No. of cigarettes/day        Granulocyte (%)
0                                      45-55  21   0.95(0.88-1.01)
1-10                                   55-61  21   0.94(0.82-1.06)
>10                                    61-76  21   0.96(0.91-1.02)
Smoking                             Years of
duration (years)                  employment
0                                         <9  20   0.88(0.75-1.00)
0-19                                    9-21  17   0.96(0.91-1.02)
>19                                      >21  22  0.99 (0.95-1.04)
Pack-years
0
0-4
>4

Variable                   p-Value (a)           Mean OD  p-Value (a)
                                                (95% CI)
Aqe (years)
<39                                     0.25 (0.05-0.45)
39-47                                   0.43 (0.31-0.56)
>47                               0.98  0.55 (0.36-0.74)         0.04
BMI (kg/[m.sup.2])
<25                                     0.30 (0.05-0.64)
25-27.5                                 0.49 (0.38-0.60)
>27.5                             0.94  0.26 (0.02-0.50)         0.43
Smoking

Nonsmoker                               0.39 (0.15-0.92)
Former smoker                           0.46 (0.33-0.581
Current smoker                    0.26  0.36 (0.18-0.53)         0.57
No. of cigarettes/day
0                                        0.43 (0.3-0.57)
1-10                                    0.42 (0.19-0.65)
>10                               0.76  0.41 (0.28-0.54)         0.86
Smoking
duration (years)
0                                       0.22 (0.04-0.39)
0-19                                    0.43 (0.28-0.58)        0.006
>19                               0.04  0.51 (0.37-0.66)
Pack-years
0
0-4
>4

(a)p-Value test for trend. (b)p-Value, one-way analysis of variance.


Association of levels of exposure to inbalable metals with histone modifications. We evaluated whether levels of H3K4me2 and H3K9ac were associated with the levels of personal exposure to metals in inhalable PM, as well as to PM mass, in both simple regression models and multivariable models adjusted for age, BMI, education, pack-years, and percent granulocytes. Results from unadjusted and adjusted models showed similar results (Table 4). In adjusted models, H3K4me2 increased in association with nickel, arsenic, and iron. H3K4me2 was not associated with aluminum, manganese, zinc, lead, PM10, or PM 1 exposure levels in unadjusted or adjusted regression models (Table 4). H3K9ac was positively but not significantly associated with nickel and iron levels. H3K9ac was not associated with the levels of the remaining metals, [PM.sub.10] or [PM.sub.1] (Table4).
Table 4. Association of personal level of exposure to PM metal
components and total PM mass with H3K4me2 and H3K9ac.

                       Unadjusted                   Adjusted (a)
Exposure   p-Coefficient (95% CI)  p-Vaiue  (3-Coefficient (95% CI)
Aluminium    0.01 (-0.02 to 0.04)     0.36     0.02 (-0.01 to 0.05)
Manganese    0.02 (-0.01 to 0.04)     0.21     0.02 (-0.01 to 0.05)
Nickel        0.15 (0.03 to 0.28)     0.02       0.16 (0.03 to 0.3)
Zinc          0.051-0.01 to 0.11)     0.08      0.05(-0.01 to 0.12)
Arsenic       0.16 (0.03 to 0.28)     0.02       0.16 (0.02 to 0.3)
Lead         0.02 (-0.01 to 0.04)     0.19     0.02 (-0.01 to 0.05)
Iron           0.12(0.01 to 0.24)     0.04       0.14(0.01 to 0.26)
PM,0         0.03 (-0.03 to 0.09)     0.34     0.04 (-0.03 to 0.11)
PM1           0.02 (-0.06 to 0.1)     0.64     0.03 (-0.06 to 0.13)

                                Unadjusted
Exposure   p-Value  p-Coefficient (95% CI)  p-Value
Aluminium     0.24   -0.04 (-0.11 to 0.02)     0.15
Manganese     0.17   -0.03 (-0.08 to 0.02)     0.25
Nickel        0.02    0.22 (-0.04 to 0.48)     0.10
Zinc          0.12   -0.05 (-0.18 to 0.07)     0.41
Arsenic       0.02    0.18 (-0.08 to 0.45)     0.17
Lead          0.16   -0.04 (-0.09 to 0.02)     0.16
Iron          0.03    0.21 (-0.02 to 0.45)     0.08
PM,0          0.23    -0.08 (-0.2 to 0.04)     0.17
PM1           0.45    -0.04 (-0.2 to 0.13)     0.65

                   Adjusted (a)
Exposure   p-Coefficient (95% CI)  p-Value
Aluminium    -0.03 (-0.1 to 0.03)     0.33
Manganese   -0.01 (-0.06 to 0.04)     0.73
Nickel       0.24 (-0.02 to 0.51)     0.07
Zinc        -0.06 (-0.19 to 0.07)     0.38
Arsenic      0.21 (-0.06 to 0.48)     0.13
Lead        -0.02 (-0.08 to 0.04)     0.52
Iron         0.22 (-0.03 to 0.47)     0.08
PM,0        -0.04 (-0.17 to 0.09)     0.51
PM1          0.00 (-0.18 to 0.18)     0.97

Data are regression coefficient [beta] and 95% CI expressing the change
in histone modifications (OD units) associated with an increase in
exposure equal to the difference between the 90th and 10th percentile of
the exposure distribution. (a)Mu Invariable regression models adjusted
for age, BMI, education, pack-years, and percent granulocytes.


Association of cumulative exposure to inbalable metals with histone modifications. Table 5 shows the associations of H3K4me2 and H3K9ac with cumulative exposures, estimated as the product of level of exposure to PM metal components or total PM mass and years of employment. Again, unadjusted and adjusted models showed similar results. In adjusted models, H3K4me2 increased in association with cumulative exposure to nickel and arsenic. H3K9ac showed significant associations with cumulative exposure to nickel and arsenic. Cumulative exposure to iron was positively but not significantly associated with both H3K4me2 and H3K9ac.
Table 5. Association of cumulative level of exposure (a) to PM metal
components and total PM mass with H3K4me2 and H3K9ac.

                                                   Adjusted (b)
Exposure      [beta]-Coefficient   p-Value   [beta]-Coefficient
                   (95% CI)                         (95% CI)

Aluminium    0.01 (-0.03 to 0.04)     0.61  0.02 (-0.02 to 0.06)
Manganese    0.01 (-0.02 to 0.04)     0.33  0.02 (-0.02 to 0.05)
Nickel       0.11 (-0.01 to 0.23)     0.07    0.16 (0.01 to 0.3)
Zinc         0.06 (-0.03 to 0.15)     0.21  0.06 (-0.05 to 0.17)
Arsenic       0.14 (0.03 to 0.26)     0.01   0.16 (0.03 to 0.29)
Lead         0.02 (-0.02 to 0.06)     0.35  0.02 (-0.03 to 0.06)
Iron         0.09 (-0.03 to 0.21)     0.14  0.12 (-0.02 to 0.27)
[PM.sub.10]  0.03 (-0.04 to 0.11)     0.41  0.05 (-0.04 to 0.14)
[PM.sub.1]   0.03 (-0.07 to 0.13)     0.53  0.07 (-0.06 to 0.21)

                              Unadjusted

Exposure     p-Value    [beta]-Coefficient     p-Value
                                    (95%CI)

Aluminium       0.44  -0.02 (-0.09 to 0.04)     0.48

Manganese       0.35  -0.01 (-0.07 to 0.04)     0.63

Nickel          0.03    0.26 (0.04 to 0.49)     0.02

Zinc            0.30  -0.02 (-0.21 to 0.16)     0.80

Arsenic         0.02    0.30 (0.09 to 0.51)     0.01

Lead            0.38   -0.031-0.11 to 0.05)     0.42

Iron            0.09    0.25 (0.03 to 0.48)     0.03

[PM.sub.10]     0.26   0.00 (-0.14 to 0.15)     0.99

[PM.sub.1]      0.29    0.081-0.11 to 0.28)     0.40

                   Adjusted (b)
Exposure       [beta]-oefficient    p-Value
                           (95%CI)

Aluminium     -0.03 (-0.1 to 0.05)     0.48
Manganese     0.00 (-0.06 to 0.05)     0.86
Nickel         0.27 (0.01 to 0.54)     0.04
Zinc         -0.03 (-0.24 to 0.18)     0.78
Arsenic        0.28 (0.04 to 0.51)     0.02
Lead          -0.02 (-0.1 to 0.06)     0.65
Iron           0.24 (-0.03 to 0.5)     0.08
[PM.sub.10]   0.00 (-0.17 to 0.17)     0.97
[PM.sub.1]    0.08 (-0.17 to 0.32)     0.54

Data are regression coefficient (|[beta]) and 95% CI expressing the
change in histone modifications (OD units) associated with an increase
in exposure equal to the difference between the 90th and 10th
percentile of the exposure distribution.
(a)Cumulative levels of exposure were estimated as the product of
personal level of exposure to PM metal components or total PM mass
and years of employment in the job. (b)Multivariable regression models
adjusted for age, BMI, education, pack-years, and percent granulocytes.


Discussion

The present study, based on a healthy worker population from a steel plant near Brescia, Italy, showed that exposure to some metal components of PM was associated with increased activating histone modifications measured in blood leukocyte samples. In particular, we found that both H3K4me2 and H3K9ac were higher in individuals with more years of employment in the plant and higher estimated cumulative exposures to arsenic and nickel.

To the best of our knowledge, this is the first study showing associations between metal exposures, such as nickel, arsenic, and iron, and histone modifications in human subjects. Our study was based on measures of histone modifications in blood leukocyte samples from healthy subjects, thus suggesting that exposure-related alterations of histone modifications may occur in normal tissues and possibly anticipate the onset of disease. Estimated effects of metal exposures on H3K4me2 in blood leukocytes were consistent with previous in vitro toxicology studies that showed that carcinogenic metals increased H3K4me2 in A549 human lung carcinoma cells (Sun et al. 2009; Zhou X et al. 2009). In particular, nickel and arsenic, which were most consistently associated with increased H3K4me2 in our analyses, were also previously found to increase H3K4me2 in vitro (Zhou X et al. 2009)- In our study, we found that H3K4me2 was correlated with the levels of exposure to nickel, arsenic, and iron, but only nickel and arsenic showed significant associations with H3K4me2 when we evaluated cumulative exposures. Also, we found significant associations of cumulative exposures to nickel and arsenic with increased H3K9ac.

The changes observed in previous in vitro studies on lung carcinoma cells (Sun et al. 2009; Zhou X et al. 2009) and in our study in blood leukocyte samples suggest that induction of gene-activating histone modifications such as H3K4me2 and H3K9ac might be a systemic process detectable across different tissues. However, whether H3K4me2 and H3K9ac are induced in human lung tissues exposed in vivo to carcinogenic metals needs to be confirmed in future investigations. Metal components of inhaled PM have been shown to induce oxidative stress and inflammatory processes, which are known to affect histone modifications (Donaldson et al. 2003; Gilmour et al. 2003) and might specifically affect blood leukocyte measures. A series of experimental studies using ambient PM10 collected in Utah Valley (USA) near a local steel plant demonstrated that anthropogenic PM containing bioavailable transition metals have heightened acute inflammatory effects, and that PM oxidant generation ability was enhanced in PM with higher metal content (Dye et al. 2001; Frampton et al. 1999; Ghio and Devlin 2001). Previous studies have shown that exposures to PM or related airborne pollutants are associated with changes in proinflammatory and cancer-related gene expression in blood leukocytes (Feinberg and Tycko 2004; Franklin et al. 2008; MacNee and Donaldson 2003; Salnikow and Zhitkovich 2008) and with markers of gene expression control such as DNA methylation (Baccarelli et al. 2009; Boleti et al.2007; Pavanello et al. 2009; Tarantini et al. 2009).

We investigated a population with well-characterized exposure that allowed us to compare subjects over a wide range of different exposure levels. We controlled for several potential confounders by fitting multivariable models that included several individual characteristics as independent variables. However, we cannot exclude the possibility that other unmeasured exposures that are present in foundry facilities, such as heat, carbon monoxide, and nonionizing radiation (Gomes et al. 2002; Lewis et al. 1992), might have influenced H3K4me2 and H3K9ac. In univariate analyses, we found that H3K9ac was positively associated with the education level of the study subjects. Although our study did not provide information to evaluate the biological basis for this association, our results showed that adjusted and unadjusted estimates for metal effects on histone modifications were remarkably similar, so confounding from education and other variables included in our models is unlikely.

Because of the limited number of study subjects, it is possible that the associations observed were attributable to chance. However, the occupational exposure and relatively controlled environment of a foundry provide a good setting for evaluating these mechanistic questions and limit bias and chance findings. Our study was based on subjects working in several work areas of the same factory and did not include a different population of subjects without a specific condition of exposure to inhaled pollutants. Limiting our investigation to individuals who have all been worKing in the same work facility avoided potential concerns related to the selection of external referents who might have differed from the exposed population in terms of socioeconomic factors and other characteristics determining hiring into the plant (Pearce et al. 2007). Nonetheless, the differences in the personal levels of exposure in our study group were large, providing sufficient contrast for identifying exposure-related changes in histone modifications. For example, the lowest level of [PM.sub.10] observed in our study population (73.72 [micro] g/ [m.sup.3]) was only marginally higher than ambient [PM.sub.10] levels measured in the geographic area in which the plant is located (average annual ambient PMjo levels between 41 and 57 [micro] g/[m.sup.3] were recorded in the year of the study by different ambient monitoring stations in Brescia area) (Anselmi and Patelli 2006), whereas the highest level was 1220.17 ag/m3. The subjects with the highest exposures to arsenic and nickel had exposure levels that were 36 and 31 times higher, respectively, than those for the least-exposed subjects. Although we based our study on a group of steel workers with average exposures higher than that of the general population, the levels of exposure to metals in our study were all lower than the commonly accepted threshold limits for industrial settings (American Conference of Governmental Industrial Hygienists 2009).

Because we determined differences in exposures within our study population by.the different tasks routinely performed by each of the study subjects, the personal exposure levels we measured in the week of the study reflected the usual exposure of the study subjects. We confirmed this by the high correlation ([r.sup.2] Greater than 0.90) between exposure measures repeated over 3 years in a subset of the study population. Employment records showed that all the subjects included in the present study performed the same tasks for all the years they had been employed in the plant, suggesting that contrasts of exposure within the study population might have remained stable over time. Therefore, we estimated long-term cumulative exposure as the product of the levels of exposure by years of employment in the plant. However, because most of the study subjects had worked for more than a decade in the plant, it is possible that changes in production or exposure protection regulations might have occurred over such an extended period of time. Therefore, we recognize that the estimated cumulative exposures we used in the present study are prone to exposure misclassification, and the results based on these metrics of exposure should be interpreted with caution.

Conclusions

Our results suggest that exposure to some metal components of PM, including nickel and arsenic, increased the two genes activating H3K4me2 and H3K9ac in blood leukocytes collected from healthy steel workers. Changes in the genomic levels of histone modifications may affect gene expression and contribute to the carcinogenic properties of inhalable nickel and arsenic. Further studies are required to directly link these changes with exposure-related increases in the risk of cancer, as well as to identify specific genes and pathways that are affected by the exposure-related changes in histone modifications.

References

Alley D, Langley-Turnbaugh S, Gordon N, Wise J, Van Epps G, Jalbert A. 2009. The effect of [PM.sub.10] on human lung fibroblasts. Toxicol Ind Health 25(2): 111-120.

American Conference of Governmental Industrial Hygienists. 2009. 2009 Threshold Limit Values and Biological Exposure Indices. Cincinnati: American Conference of Governmental Industrial Hygienists.

Anselmi U, Patelli R. 2006. Rapporto sulla qualita dell'aria di Brescia e provincia [in Italian]. Milan: ARPA Lombardia.

Baccarelli A, Wright RO, Bollati V, Tarantini L, Litonjua AA, Suh HH, et al. 2009. Rapid DNA methylation changes after exposure to traffic particles. Am J Respir Crit Care Med 179(7): 572-578.

Bollati V, Baccarelli A, Hou L, Bonzini M, Fustinoni S, Cavallo D, et al. 2007. Changes in DNA methylation patterns in subjects exposed to low-dose benzene. Cancer Res 67(3): 876-880.

Chang CC, Hwang JS, Chan CC, Wang PY, Hu TH, Cheng TJ. 2005. Effects of concentrated ambient particles on heart rate variability in spontaneously hypertensive rats. J Occup Health 47(6): 471-480.

Chen H, Ke Q, Kluz T, Yan Y, Costa M. 2006. Nickel ions increase histone H3 lysine 9 dimethylation and induce transgene silencing. Mol Cell Biol 26(10): 3728-3737.

Chen LCr Hwang JS. 2005. Effects of subchronic exposures to concentrated ambient particles (CAPs) in mice. IV. Characterization of acute and chronic effects of ambient air fine particulate matter exposures on heart-rate variability. Inhal Toxicol 17(4-5): 209-216.

Conroy J, Byrne SJ, Gun'ko YK, Rakovich YP, Donegan JF, Davies A, et al. 2008. CdTe nanoparticles display tropism to core histones and histone-rich cell organelles. Small 4(11):2006-2015.

Corey LM, Baker C, Luchtel DL 2006. Heart-rate variability in the apolipoprotein E knockout transgenic mouse following exposure to Seattle particulate matter. J Toxicol Environ Health A 69(10) 953-965.

Coyle YM, Minahjuddin AT, Hynan LS, Minna JD. 2006. An ecological study of the association of metal air pollutants with lung cancer incidence in Texas. J Thorac Oncol 1(7): 654-661.

Deligezer U, Yaman F, Darendeliler E, Dizdar Y, Holdenrieder S, Kovancilar M, et al. 2010. Post-treatment circulating plasma BMP6 mRNA and H3K27 methylation levels discriminate metastatic prostate cancer from localized disease. Clin Chim Acta 411(19-20): 1452-1456.

Dockery DW, Pope CA III, Xu X, Spengler JD, Ware JH, Fay ME, et al. 1993. An association between air pollution and mortality in six U.S. cities. N Engl J Med 329(24): 1753-1759.

Donaldson K, Stone V, Borm PJ, Jimenez LA, Gilmour PS, Schins RP, et al. 2003. Oxidative stress and calcium signaling in the adverse effects of environmental particles (PM10). Free Radic Biol Med 34(11): 1369-1382.

Dye JA, Lehmann JR, McGee JK, Winsett DW, Ledbetter AD, Everitt Jl, et al. 2001. Acute pulmonary toxicity of particulate matter filter extracts in rats: coherence with epidemiologic studies in Utah Valley residents. Environ Health Perspect 109(suppl 3):395-403.

Fang GC, Lin SJ, Lee JF, Chang CC. 2009. A study of particulates and metallic element concentrations in temple. Toxicol Ind Health 25(2): 93-100.

Feinberg AP, Tycko B. 2004. The history of cancer epigenetics. Nat Rev Cancer 4(2): 143-153.

Frampton MW, Ghio AJ, Samet JM, Carson JL, Carter JD, Devlin RB. 1999. Effects of aqueous extracts of PM10 filters from the Utah Valley on human airway epithelial cells. Am J Physiol 277(5 pt1): L960-L967.

Franklin M, Koutrakis P, Schwartz P. 2008. The role of particle composition on the association between PM2.5 and mortality. Epidemiology 19(5): 680-689.

Ghio AJ, Devlin RB. 2001. Inflammatory lung injury after bronchial instillation of air pollution particles. Am J Respir Crit Care Med 164(4): 704-708.

Gibb HJ, Lees PS, Pinsky PF, Rooney BC. 2000. Lung cancer among workers in chromium chemical production. Am J Ind Med 38(2): 115-126.

Gilmour PS, Rahman I, Donaldson K, MacNee W. 2003. Histone acetylation regulates epithelial IL-8 release mediated by oxidative stress from environmental particles. Am J Physiol Lung Cell Mol Physiol 284(3): L533-L540.

Gomes J, Lloyd 0, Norman N. 2002. The health of the workers in a rapidly developing country: effects of occupational exposure to noise and heat. Occup Med (Lond) 52(3): 121-128.

International Agency for Research on Cancer. 1987. Overall evaluations of carcinogenicity: an updating of IARC Monographs volumes 1 to 42. Monogr Eval Carcinog Risks Hum Suppl 7: 1-440.

Jensen TJ, Wozniak RJ, Eblin KE, Wnek SM, Gandolfi AJ, Futscher BW. 2009. Epigenetic mediated transcriptional activation of WIMT5A participates in arsenical-associated malignant transformation. Toxicol Appal Pharmacol 235(1): 39-46.

Ke Q, Davidson T, Chen H, Kluz T, Costa M. 2006. Alterations of histone modifications and transgene silencing by nickel chloride. Carcinogenesis 27(7): 1481-1488.

Kuo HW, Chang CL, Liang WM, Chung BC. 1999. Respiratory abnormalities among male foundry workers in central Taiwan. Occup Med (Lond) 49(8): 499-505.

Lewis S, Mason C, Srna J. 1992. Carbon monoxide exposure in blastfurnace workers. Aust J Public Health 16(3): 262-268.

MacArthur AC, Le ND, Fang R, Band PR. 2009. Identification of occupational cancer risk in British Columbia: a population-based case-control study of 2,998 lung cancers by histopathologica subtype. Am J Ind Med 52 (3): 221-232.

MacNee W, Donaldson K. 2003. Mechanism of lung injury caused by PM10 and ultrafine particles with special reference to COPD. Eur Respir J Suppl 40: 47s-51s.

Pavanello S, Bollati V, Pesatori AC, Kapka L, Bolognesi C, Bertazzi PA, et al. 2009. Global and gene-specific promoter methylation changes are related to anti-B[a|PDE - DNA adduct levels and influence micronuclei levels in poly-cyclic aromatic hydrocarbon-exposed individuals. Int J Cancer 125(7): 1692-1697.

Pearce N, Checkoway H, Kriebel D. 2007. Bias in occupational epidemiology studies. Occup Environ Med 64(8): 562-568.

Roller M. 2009. Carcinogenicity of inhaled nanoparticles. Inhal Toxicol 21(suppl 1): 144-157.

Roy AA, Baxla SP, Gupta T, Bandyopadhyaya R, Tripathi SN. 2009. Particles emitted from indoor combustion sources: size distribution measurement and chemical analysis. Inhal Toxicol 21(10): 837-848.

Salnikow K, Zhitkovich A. 2008. Genetic and epigenetic mechanisms in metal carcinogenesis and cocarcinogen-esis: nickel, arsenic, and chromium. Chem Res Toxicol 21(1): 28-44.

Sun H, Zhou X, Chen H, Li Q, Costa M. 2009. Modulation of histone methylation and MLH1 gene silencing by hexavalent chromium. Toxicol Appl Pharmacol 237(3): 258-266.

Surralles J, Puerto S, Ramirez MJ, Creus A, Marcos Rf Mullenders LH, et al. 1998. Links between chromatin structure, DNA repair and chromosome fragility. Mutat Res 404(1-2): 39-44.

Tarantini L, Bonzini M, Apostoli P, Pegoraro V, Bollati V, Marinelli B, et al. 2009. Effects of particulate matter on genomic DNA methylation content and iNOS promoter methylation. Environ Health Perspect 117: 217-222.

Vijay Bhaskar B, Jeba Rajasekhar RV, Muthusubramanian P, Kesarkar AP. 2009. Ionic and heavy metal composition of respirable particulate in Madurai, India. Environ Monit Assess 164(l-4): 323-336.

Wang Z, Neuburg D, Li C, Su L, Kim JY, Chen JC, et al. 2005. Global gene expression profiling in whole-blood samples from individuals exposed to metal fumes. Environ Health Perspect 113: 233-241.

Wild P, Bourgkard E, Paris C. 2009. Lung cancer and exposure to metals; the epidemiological evidence. Methods Mol Biol 472: 139-167.

Yan C, Boyd DD. 2006. Histone H3 acetylation and H3 K4 methylation define distinct chromatin regions permissive for transgene expression. Mol Cell Biol 26(17): 6357-6371.

Yan Y, Kluz T, Zhang P, Chen HB, Costa M. 2003. Analysis of specific lysine histone H3 and H4 acetylation and methyia-tion status in clones of cells with a gene silenced by nickel exposure. Toxicol Appl Pharmacol 190(3):272-277.

Zhou QX, Zhao XN, Xiao JP. 2009. Preconcentration of nickel and cadmium by Ti02 nanotubes as solid-phase extraction adsorbents coupled with flame atomic absorption spectrometry. Talanta 77(5): 1774-1777.

Zhou X, Li Q, Arita A, Sun H, Costa M. 2009. Effects of nickel, chromate, and arsenite on histone 3 lysine methylation. Toxicol Appl Pharmacol 236(1): 78-84.

Zhou X, Sun H, Ellen TP, Chen H, Costa M. 2008. Arsenite alters global histone H3 methylation. Carcinogenesis 29(91: 1831-1836.

Zoroddu MA, Kowalik-Jankowska T, Kozlowski H, Molinari H, Salnikow K, Broday L, et al. 2000. Interaction of Ni(ll) and Cu(ll) with a metal binding sequence of histone H4: AKRHRK, a model of the H4 tail. Biochim Biophys Acta 1475(2): 163-168.

Zoroddu MA, Schinocca L, Kowalik-Jankowska T, Kozlowski H, Salnikow K, Costa M. 2002. Molecular mechanisms in nickel carcinogenesis: modeling Ni(ll) binding site in histone H4. Environ Health Perspect 110(suppl 5): 719-723.

Laura Cantone, (1) Francesco Nordio, (1), (2) Lifang Hou, (3) Pietro Apostoli, (4) Matteo Bonzini, (5) Letizia Tarantini, (1) Laura Angelici, (1) Valentina Bollati, (1) Antonella Zanobetti, (6) Joel Schwartz, (6) Pier A, Bertazzi,1 and Andrea Baccarelli (6)

(1)Department of Environmental and Occupational Health, Universita di Milano and Istituto Di Ricovero e Cura a Carattere Scientifico, and Maggiore Hospital, Mangiagalli and Regina Elena Foundation, Milan, Italy; (2)Department of Clinical Medicine, Nephrology and Health Sciences, University of Parma Medical School, Parma, Italy; (3)Department of Preventive Medicine, Feinberg School of Medicine, Northwestern University, Chicago, Illinois, USA; (4)Department of Experimental and Applied Medicine, Occupational Medicine and Industrial Hygiene, University of Brescia, Brescia, Italy; (5)Department of Clinical and Biological Sciences, University of Insubria, Varese, Italy; (6)Laboratory of Environmental Epigenetics, Exposure Epidemiology and Risk Program, Department of Environmental Health, Harvard School of Public Health, Boston, Massachusetts, USA
COPYRIGHT 2011 National Institute of Environmental Health Sciences
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2011 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Title Annotation:Research
Author:Cantone, Laura; Nordio, Francesco; Hou, Lifang; Apostoli, Pietro; Bonzini, Matteo; Tarantini, Letizi
Publication:Environmental Health Perspectives
Article Type:Report
Geographic Code:4EUIT
Date:Jul 1, 2011
Words:7815
Previous Article:Shorter anogenital distance predicts poorer semen quality in young men in Rochester, New York.
Next Article:Determinants of agricultural pesticide concentrations in carpet dust.
Topics:

Terms of use | Privacy policy | Copyright © 2019 Farlex, Inc. | Feedback | For webmasters