Associations of Blood Lead, Dimercaptosuccinic Acid-Chelatable Lead, and Tibia Lead with Polymorphisms in the Vitamin D Receptor and [Delta]-Aminolevulinic Acid Dehydratase Genes.

A cross-sectional study was performed to evaluate the influence of polymorphisms in the [Delta]-aminolevulinic acid dehydratase (ALAD) and vitamin D receptor (VDR) genes on blood lead, tibia lead, and dimercaptosuccinic acid (DMSA)-chelatable lead levels in 798 lead workers and 135 controls without occupational lead exposure in the Republic of Korea. Tibia lead was assessed with a 30-min measurement by [sup.109]Cd-induced K-shell X-ray fluorescence, and DMSA-chelatable lead was estimated as 4-hr urinary lead excretion after oral administration of 10 mg/kg DMSA. The primary goals of the analysis were to examine blood lead, tibia lead, and DMSA-chelatable lead levels by ALAD and VDR genotypes, controlling for covariates; and to evaluate whether ALAD and VDR genotype modified relations among the different lead biomarkers. There was a wide range of blood lead (4-86 [micro]g/dL), tibia lead (-7-338 [micro]g Pb/g bone mineral), and DMSA-chelatable lead (4.8-2,103 [micro]g) levels among lead workers. Among lead workers, 9.9% (n = 79) were heterozygous for the [ALAD.sup.2] allele and there were no homozygotes. For VDR, 10.7% (n = 85) had the Bb genotype, and 0.5% (n = 4) had the BB genotype. Although the ALAD and VDR genes are located on different chromosomes, lead workers homozygous for the [ALAD.sup.1] allele were much less likely to have the VDR bb genotype (crude odds ratio = 0.29, 95% exact confidence interval = 0.06-0.91). In adjusted analyses, subjects with the [ALAD.sup.2] allele had higher blood lead levels (on average, 2.9 [micro]g/dL, p = 0.07) but no difference in tibia lead levels compared with subjects without the allele. In adjusted analyses, lead workers with the VDR B allele had significantly (p [is less than] 0.05) higher blood lead levels (on average, 4.2 [micro]g/dL), chelatable lead levels (on average, 37.3 [micro]g), and tibia lead levels (on average, 6.4 [micro]g/g) than did workers with the VDR bb genotype. The current data confirm past observations that the ALAD gene modifies the toxicokinetics of lead and also provides new evidence that the VDR gene does so as well. Key word: [Delta]-aminolevulinic acid dehydratase, bone lead, cross-sectional study, lead, polymorphisms, vitamin D receptor, X-ray fluorescence. Environ Health Perspect 108:949-954 (2000). [Online 31 August 2000]

http://ehpnet1.niehs.nih.gov/docs/2000 /108p949-954schwartz/abstract.html

An increasing body of evidence suggests that genetic factors modify the toxicokinetics of lead. The gene for the [Delta]-aminolevulinic acid dehydratase (ALAD) enzyme has been a focus of primary interest. Human ALAD is encoded by a single gene on chromosome 9p34 that has two alleles, [ALAD.sup.1] and [ALAD.sup.2], resulting in three isozymes, ALAD1-1, ALAD1-2, and ALAD2-2 (1,2). The prevalence of the [ALAD.sup.2] allele is approximately 10% in Asians and 20% in Caucasians (3-5). Subjects who have at least one copy of the [ALAD.sup.2] allele, compared to subjects with none, have higher blood lead levels (3-5); lower dimercaptosuccinic acid (DMSA)-chelatable lead levels (6); lower plasma aminolevulinic acid levels (7); a larger difference between trabecular and cortical bone lead levels (8); higher blood urea nitrogen and serum creatinine levels (8); less efficient uptake of lead into bone, especially trabecular bone (9); lower zinc protoporphyrin (ZPP) levels for given levels of blood lead (10); and lower urinary calcium and creatinine levels (11). ALAD has been identified as a principal lead-binding protein, and the proportion of lead bound to ALAD was greater for subjects with [ALAD.sup.2] (12).

Recent data suggest that polymorphisms in the vitamin D receptor (VDR) gene influence tibia lead levels (13). The VDR gene is located at chromosome 12cen-12 (14), and thus variant VDR alleles would not be expected to be linked with variant ALAD alleles. Most studies of the VDR gene have focused on the BsmI polymorphism; restriction enzyme digestion produces three genotypes commonly termed bb, Bb, and BB. The BB allele (defined by the absence of the polymorphic restriction site) has a prevalence of 7-32% in Caucasians (15). Study subjects (mainly women) with the BB genotype have bone mineral densities up to 10-15% lower than subjects with the bb genotype, with an overall difference across studies of 2-2.5% reported in a recent meta-analysis (15). As lead and calcium are known to behave similarly in biologic systems, these findings have motivated investigations of VDR genotype and bone lead levels. Subjects with the B allele had larger tibia lead concentrations with increasing age and lower tibia lead concentrations with increasing duration since last exposure to lead than did subjects without the B allele (13).

No prior studies have evaluated the joint influence of the VDR and ALAD genes on bone, blood, and chelatable lead levels. Here we report the largest such study to date in a cross-sectional analysis of 798 Korean lead workers and 135 controls without occupational lead exposure.

Materials and Methods

Study overview and design. The results presented here are a cross-sectional analysis of data from the first year of a 3-year longitudinal study of the health effects of occupational inorganic lead exposure (16). Enrollment began in October 1997 with the first of three annual evaluations for each study subject. The current report is focused on the influence of the ALAD and VDR polymorphisms on blood lead, tibia lead, and dimercaptosuccinic acid (DMSA)-chelatable lead levels measured during this first study evaluation and is based on the 798 lead workers and 135 controls without occupational lead exposure who were enrolled between 24 October 1997 and 19 August 1999. The study was reviewed and approved by institutional review boards at the Johns Hopkins School of Hygiene and Public Health, Baltimore, Maryland, and the Soonchunhyang University School of Medicine, Chonan, Korea.

Study population. Participation in the study was voluntary, and all participants provided written, informed consent. Subjects were paid approximately $30 for their participation. Lead workers were recruited from 24 different lead-using facilities, with participation in most facilities exceeding 80% (16). Retired workers from three facilities who had received medical surveillance services by Soonchunhyang University for several years were also recruited to participate in the study. Routine, government-mandated industrial hygiene sampling revealed that the study plants did not have significant amounts of other heavy metals such as cadmium. Controls without occupational lead exposure were recruited from an air conditioner assembly plant that did not use lead or other heavy metals and from hourly wage employees of Soonchunhyang University.

Data collection. Data collection methods have been previously reported (16). In brief, data were collected either at the Institute of Industrial Medicine at Soonchunhyang University in Chonan or on the premises of the lead-using facilities. The following data were collected or measured on all study subjects: a standardized interview for demographics, medical history, and occupational history; a neurobehavioral test battery consisting of examiner-administered tests; blood pressure; peripheral vibration threshold and pinch and grip strength; a 10-mL blood specimen by venipuncture that was stored at -70 [degrees] C as whole blood, plasma, and red blood cells; a spot urine sample; tibia lead concentration by X-ray fluorescence; and a 4-hour urine sample after oral administration of DMSA (in lead workers only).

Laboratory methods. Hemoglobin was assayed by the cyanmethemoglobin method (Beckman Coulter, Inc., Model Ac-T 8, Fullerton, CA), and hematocrit was measured by the capillary centrifugation method (17). Zinc protoporphyrin levels were measured with a portable hematofluorimeter (18). Urinary creatinine was measured using a Sigma kit (St. Louis, MO) (19). Blood lead levels were measured with a Zeeman background-corrected atomic absorption spectrophotometer (Hitachi Z-8100 model; Hitachi, Ltd., Tokyo, Japan) with the standard addition method of the National Institute of Occupational Safety and Health (20) at Soonchunhyang University Institute of Industrial Medicine, a certified reference laboratory for lead in Korea. Tibia lead was assessed, in units of micrograms lead per gram bone mineral, with a 30 min measurement at the left mid-tibial shaft using [sup.109]Cd-induced K-shell X-ray fluorescence (XRF), as previously described (21-23). XRF can provide negative point estimates of bone lead concentrations; however, all point estimates were retained in the statistical analyses, including negative values, because this method minimizes bias and does not require censoring of data (24).

We used 4-hour urinary lead excretion after oral administration of 10 mg/kg DMSA to measure DMSA-chelatable lead (25). Urine lead levels were measured in the laboratories of the Wadsworth Center at the New York State Department of Health, Albany, New York. Urinary lead concentrations were determined by electrothermal atomization atomic absorption spectrometry (Perkin-Elmer Model 4100ZL, Norwalk, CT) using previously published methods (26). Urinary lead excretion was highly correlated with lead excretion adjusted for differences, generally small, in urine collection times (Pearson's r = 0.98), so only the unadjusted data are presented.

ALAD and VDR genotyping. We completed ALAD and VDR genotyping on 798 and 795 subjects, respectively. VDR genotyping was completed using previously published methods (13). In brief, genomic DNA was extracted from whole blood by using the QIAamp Blood Kit (QIAGEN, Hilden, Germany), and the BsmI polymorphic site in intron 8 was amplified by polymerase chain reaction (PCR) using the primers originating in exon 7 (primer 1: 5'-CAACCAAGAC-TACAAGTACC-GCGTCAGTGA-3') and intron 8 (primer 2: 5'-AACCAGCGGGAAGAGGTCAAGGG-3'). Subjects homozygous for the presence of the BsmI restriction site are designated bb, heterozygotes are designated Bb, and those homozygous for the absence of the site are designated BB.

A modified PCR-based protocol was used for ALAD genotyping and has been previously described (3-5). In brief, the initial amplification, using 3' and 5' oligonucleotide primers (5'- AGACAGACATTAGCTCAGTA-3') and (5'-GGCAAAGAACACGTCCATTC-3'), generates a 916 base pair fragment. A second round of amplification uses a pair of nested primers (provided by J. Wetmur), sequences (5'-CAGAGCTGTTCCAACAGTGGA-3') and (5'-CCAGCACAATGTGGGAGTGA-3'), respectively, and generates an 887 base pair fragment. The amplified fragment was cleaved at the diagnostic Msp1 site, only present in the [ALAD.sup.2] allele, and three isozymes are observed, designated ALAD1-1, ALAD1-2, and ALAD2-2.

Statistical analysis. The primary goals of the analysis were to examine blood lead, tibia lead, and DMSA-chelatable lead levels by ALAD and VDR genotypes, controlling for covariates; and to evaluate whether ALAD and VDR genotype modified relations among the different lead biomarkers. Associations between ALAD and VDR genotype were evaluated in contingency tables using odds ratios (ORs) and 95% exact confidence intervals (CIs) calculated with Epi Info version 6.04b (Centers for Disease Control and Prevention, Atlanta, GA).

We used linear regression to separately model blood lead, tibia lead, and DMSA-chelatable lead levels, controlling for confounding variables, using statistical software programs of SAS Institute, Inc. (Cary, NC). In these regression models, only lead workers, not controls, were included. Covariates examined in linear regression models included age, sex, creatinine clearance (4-hr), hemoglobin, hematocrit, weight, height, body mass index, job duration, and tobacco and alcohol consumption (never, previous, and current use for each). Covariates were retained in the final regression models if they were either a significant predictor of blood lead, tibia lead, or DMSA-chelatable lead levels, or if they were a confounder of the relations between predictor variables and the lead biomarkers. The decisions regarding the variables in the final regression models were also made to be consistent with prior analyses of the data on the subjects presented here (23). Blood lead was modeled with and without adjustment for hematocrit; as this adjustment did not influence regression results, only unadjusted model results are presented.

DMSA-chelatable lead and tibia lead were log-transformed before regressing on covariates because of departures from normality. To estimate the mean adjusted differences between genotypes in the original scale of each lead measure, we exponentiated the predicted value from the regression, separately for each genotype, at the mean value of all continuous covariates and the reference value of all dichotomous covariates. To evaluate nonlinear relations, quadratic terms for continuous variables (i.e., age, job duration, weight, height) were evaluated. We evaluated effect modification by genotype by including cross-product terms between the genetic variables and relevant predictor variables (i.e., age, sex, tibia lead, creatinine clearance).

Results

Demographics and dose measures. Compared to controls without occupational lead exposure, lead-exposed subjects were older (40.5 vs. 34.5 years), had lower education levels (49.9% vs. 19.2% did not complete high school), and a lower proportion were male (79.4% vs. 91.9%; Table 1). The majority of both nonexposed and exposed subjects were current users of tobacco and alcohol products. There was a wide range of blood lead (4-86 [micro]g/dL), tibia lead (-7-338 [micro]g/g), and DMSA-chelatable lead (4.8-2,103 [micro]g) levels among lead workers (Table 1). The corresponding values among nonexposed control subjects were low. Among lead workers, tibia lead was moderately correlated with blood lead (Pearson's r = 0.42), DMSA-chelatable lead (r = 0.43), and job duration (r = 0.40) (all p-values [is less than] 0.01). The correlations of blood lead (r = 0.13) and DMSA-chelatable lead (r = 0.17) with job duration were much lower than were the correlations of these variables with tibia lead. Blood lead was highly correlated with DMSA-chelatable lead (r = 0.82).

Table 1. Description of study subjects, October 1997 to August 1999, Republic of Korea.

                                   Lead-exposed
Characteristic                     subjects (n = 798)

Age (years)                         40.5 [+ or -] 10.1 (17.8-64.8)
Lead work job duration (years)       8.2 [+ or -] 6.5 (0.1-36.2)
Height (cm)                        164.7 [+ or -] 8.1 (127.8-186.0)
Weight (kg)                         62.5 [+ or -] 9.1 (37.4-92.7)
Body mass index (kg/[cm.sup.2])     23.0 [+ or -] 3.0 (15.7-34.2)
Blood lead ([micro]g/dL)            32.0 [+ or -] 15.0 (4-86)
Tibia lead, ([micro]g Pb/g          37.2 [+ or -] 40.4 (-7-338)
  bone mineral)
DMSA-chelatable lead               186.0 [+ or -] 208.4 (4.8-2,103)
  ([micro]g)(b)
Hemoglobin (g/dL)                   14.2 [+ or -] 1.4 (6.5-17.9)
Creatinine clearance,              114.3 [+ or -] 33.9 (11.2-351.6)
  4-hr (mL/min)
Educational level(c)
  Lower school ([is less than        183 (23.0)
    or equal to] 6 years)
  Some middle school (7-8 years)      29 (3.6)
  Middle school graduate             155 (19.4)
    (9 years)
  Some high school (10-11 years)      31 (3.9)
  High school graduate               335 (42.0)
    (12 years)
  One or two years college            37 (4.6)
    (13-14 years)
  College graduate or more            27 (3.3)
  Missing                              1 (< 0.1)
Sex, male(c)                         634 (79.4)
Tobacco use(c)
  Never                              254 (31.9)
  Current use                        455 (57.1)
  Past use                            88 (11.0)
Alcohol use(c)
  Never                              231 (29.0)
  Current use                        518 (65.0)
  Past use                            48 (6.0)

Characteristic                     Controls (n = 135)

Age (years)                         34.5 [+ or -] 9.1 (22.0-60.2)
Lead work job duration (years)     NA(a)
Height (cm)                        167.9 [+ or -] 6.2 (148.0-183.4)
Weight (kg)                         66.9 [+ or -] 9.0 (48.0-93.5)
Body mass index (kg/[cm.sup.2])     23.7 [+ or -] 2.8 (18.5-30.1)
Blood lead ([micro]g/dL)             5.3 [+ or -] 1.8 (2-10)
Tibia lead, ([micro]g Pb/g           5.8 [+ or -] 7.0 (-11-27)
  bone mineral)
DMSA-chelatable lead               NA(a)
  ([micro]g)(b)
Hemoglobin (g/dL)                   15.3 [+ or -] 1.2 (11.1-18.2)
Creatinine clearance,              NA(a)
  4-hr (mL/min)
Educational level(c)
  Lower school ([is less than         10 (7.4)
    or equal to] 6 years)
  Some middle school (7-8 years)       3 (2.2)
  Middle school graduate              12 (8.9)
    (9 years)
  Some high school (10-11 years)       1 (0.7)
  High school graduate                93 (68.9)
    (12 years)
  One or two years college            11 (8.1)
    (13-14 years)
  College graduate or more             5 (3.7)
  Missing                              0 (0.0)
Sex, male(c)                         124 (91.9)
Tobacco use(c)
  Never                               35 (25.9)
  Current use                         87 (64.4)
  Past use                            13 (9.6)
Alcohol use(c)
  Never                               31 (23.0)
  Current use                         95 (70.4)
  Past use                             9 (6.7)

Values shown are mean [+ or -] SD except where indicated.

(a) The 4-hr urine collection was performed only in subjects who received DMSA.

(b) DMSA-chelatable lead ([micro]g) was estimated as 4-hr urinary lead excretion after oral administration of 10 mg/kg DMSA, in lead-exposed subjects only (784 subjects completed the urine collection).

(c) Values shown are number (%)

Prevalence and associations of genotypes. Among lead workers, 9.9% (n = 79) were heterozygous for the [ALAD.sup.2] allele, and there were no [ALAD.sup.2] homozygotes; 11.2% (n = 89) had at least one copy of the VDR B allele, and 0.5% (n = 4) had the BB genotype. The corresponding values for controls were 8.1% (n = 11) for the [ALAD.sup.2] allele and 8.9% (n = 12) and 0.7% (n = 1) for one and two copies of the VDR B allele, respectively. Because of the small number of subjects with the BB genotype, all subsequent analysis combined homozygous and heterozygous variant allele carriers.

In unadjusted (crude) analyses, there were no differences in age, job duration, hemoglobin, blood lead, tibia lead, or DMSA-chelatable lead by ALAD genotype (Table 2). In contrast, lead workers with the Bb or BB genotypes, compared to those with bb, were older and had higher DMSA-chelatable lead levels (both p-values [is less than] 0.05, Table 2).

Table 2. Selected demographic and lead biomarker variables (mean [+ or -] SD) by gene status in 798 lead-exposed subjects, October 1997 to August 1999, Republic of Korea.(a)

                                     ALAD genotype

Characteristic                            1-1

Number                              716
Age (years)                        40.5 [+ or -] 10.2
Job duration (years)                8.2 [+ or -] 6.6
Hemoglobin (g/dL)                  14.2 [+ or -] 1.4
Blood lead ([micro]g/dL)           31.7 [+ or -] 14.9
Tibia lead ([micro]g/g)            37.5 [+ or -] 40.6
DMSA-chelatable lead ([micro]g)   180.3 [+ or -] 181.2

                                     ALAD genotype

Characteristic                            1-2

Number                               79
Age (years)                        40.1 [+ or -] 9.7
Job duration (years)                8.2 [+ or -] 5.8
Hemoglobin (g/dL)                  14.2 [+ or -] 1.6
Blood lead ([micro]g/dL)           34.2 [+ or -] 15.9
Tibia lead ([micro]g/g)            31.4 [+ or -] 29.5
DMSA-chelatable lead ([micro]g)   161.7 [+ or -] 143.0

                                       VDR genotype

Characteristic                              bb

Number                              709
Age (years)                        40.2 [+ or -] 10.0(*)
Job duration (years)                8.4 [+ or -] 6.6
Hemoglobin (g/dL)                  14.2 [+ or -] 1.4
Blood lead ([micro]g/dL)           31.6 [+ or -] 14.8
Tibia lead ([micro]g/g)            37.1 [+ or -] 41.2
DMSA-chelatable lead ([micro]g)   173.5 [+ or -] 176.8(*)

                                       VDR genotype

Characteristic                           Bb or BB

Number                               89
Age (years)                        42.7 [+ or -] 10.3(*)
Job duration (years)                7.2 [+ or -] 5.6
Hemoglobin (g/dL)                  14.1 [+ or -] 1.4
Blood lead ([micro]g/dL)           34.8 [+ or -] 16.1
Tibia lead ([micro]g/g)            38.1 [+ or -] 33.5
DMSA-chelatable lead ([micro]g)   217.2 [+ or -] 179.7(*)

(a) ALAD and VDR genotyping were completed on 795 and 798 lead workers, respectively.

(*) p < 0.05.

The ALAD gene is located on chromosome 9, and the VDR gene is on chromosome 12. Nonetheless, an association was observed between the ALAD and VDR genotypes (Table 3). Among lead workers, subjects homozygous for the [ALAD.sup.1] allele were much less likely to have the VDR bb genotype (crude OR = 0.28; 95% CI, 0.06-0.89). In contrast, among controls, subjects homozygous for the ALAD1 allele were more likely to have the VDR bb genotype (crude OR = 2.53; 95% CI, 0.23-14.84). Although there were only two controls with ALAD1-2 and VDR Bb or BB, in the stratified analysis, the ORs between the two genes among lead workers and controls were significantly different (test for homogeneity of stratum-specific ORs, p = 0.04).

Table 3. Association of VDR genotype status by ALAD genotype status.

                                     VDR genotype

                                               Bb or BB,
ALAD                            bb, n (%)(a)   n (%)(a)

All study participants
  1-1                           743 (80.0)      96 (10.3)
  1-2                            85 (9.2)        5 (0.5)
  Total                         828 (89.1)     101 (10.9)

Lead workers
  1-1                           629 (79.2)      86 (10.8)
  1-2                            76 (9.6)        3 (0.4)
  Total                         705 (88.8)      89 (11.2)

Controls without occupational
lead exposure
  1-1                           114 (84.4)      10 (7.4)
  1-2                             9 (6.7)        2 (1.5)
  Total                         123 (91.1)      12 (8.9)

ALAD                               Total      OR (95% CI)(b)

All study participants
  1-1                           839 (90.3)   0.46 (0.14-1.15)
  1-2                            90 (9.7)
  Total                         929 (100)

Lead workers
  1-1                           715 (90.1)   0.29 (0.06-0.91)
  1-2                            79 (10.0)
  Total                         794 (100)

Controls without occupational
lead exposure
  1-1                           124 (91.9)   2.53 (0.23-14.84)
  1-2                            11 (8.1)
  Total                         135 (100)

(a) Percentage of table totals.

(b) Test for homogeneity of stratum-specific ORs, p = 0.04, indicating that the association of the two genotypes in lead workers and controls was different.

Predictors of blood lead levels in lead-exposed subjects. After adjustment for age, sex, and current tobacco use using linear regression, lead workers with the VDR B polymorphism had higher blood lead levels (p [is less than] 0.01), in models with or without control for ALAD genotype (model with ALAD genotype; Table 4). Subjects with the [ALAD.sup.2] allele also had higher blood lead levels (p [is less than] 0.05) in models without and with control for VDR genotype (model with VDR genotype; Table 4). On average, after controlling for the two genotypes, lead workers with the VDR B polymorphism had blood lead levels 4.2 [micro]g/dL higher than subjects with bb, and lead workers with the [ALAD.sup.2] allele had blood lead levels 3.6 [micro]g/dL higher than subjects without the allele. There was no evidence of gene-gene interaction in these models (evaluated by inclusion of an ALAD-VDR cross-product term). There was also no evidence of effect odification by ALAD or VDR genotype on the relations of the predictor variables with blood lead levels. The final linear regression model (Table 4) accounted for 35% of the variance in blood lead levels.

Table 4. Linear regression modeling of blood lead, Korean lead workers, 1997-1999.(a)

Independent       Units of [Beta]          [Beta]
variable          coefficient              coefficient

Age               [micro]g/dL/year             0.286
Female            [micro]g/dL                -13.782
Current smoker    [micro]g/dL                  3.406
Tibia lead        [micro]g/dL/[micro]g/g       0.131
VDR, Bx vs. bb    [micro]g/dL                  4.183
ALAD, 12 vs. 11   [micro]g/dL                  3.627

Independent       [Beta]
variable            SE     p-Value

Age               0.049    < 0.001
Female            1.382    < 0.001
Current smoker    1.081      0.002
Tibia lead        0.011    < 0.00
VDR, Bx vs. bb    1.376      0.002
ALAD, 12 vs. 11   1.445      0.010

(a) Model [r.sup.2] = 0.35.

Predictors of DMSA-chelatable lead levels in lead-exposed subjects. After adjustment for covariates (age, sex, current tobacco use, body mass index, and 4-hr creatinine clearance), subjects with the VDR B polymorphism had higher DMSA-chelatable lead levels (on average, 32%, or 37.3 [micro]g higher than subjects with VDR bb, p [is less than] 0.01). In contrast, the [ALAD.sup.2] allele was not significantly associated with chelatable lead levels (p = 0.69). The relation between creatinine clearance and DMSA-chelatable lead was modified by ALAD genotype (Table 5, model 3). The intercept for lead workers with the [ALAD.sup.2] allele was 59% lower than for lead workers with the [ALAD.sup.1] allele (p = 0.05). Among lead workers with only the [ALAD.sup.1] allele, chelatable lead levels increased 5.8 [micro]g for each increase of 10 mL/min in creatinine clearance near its mean value (p [is less than] 0.01); in contrast, among lead workers with the [ALAD.sup.2] allele, chelatable lead levels increased 15.5 [micro]g for each increase of 10 mL/min in creatinine clearance near its mean value (difference in two slopes, p = 0.04). The final linear regression models accounted for 25-26% of the variance in DMSA-chelatable lead levels. Addition of blood lead to the models of DMSA-chelatable lead increased the model [r.sup.2] to 79-80%. There were no interactions in these models between blood lead and either of the two genes.

Table 5. Linear regression modeling of DMSA-chelatable lead, Korean lead workers, 1997-1999.

                                   Units of
Independent                        [Beta]
variable                           coefficient

Model 1
  Age                              [micro]g/year
  Female                           [micro]g
  Current smoker                   [micro]g
  Body mass index                  [micro]g/kg/[cm.sup.2]
  Creatinine clearance             [micro]g/mL/min
  VDR, Bx vs. bb                   [micro]g

Model 2
  Age                              [micro]g/year
  Female                           [micro]g
  Current smoker                   [micro]g
  Body mass index                  [micro]g/kg/[cm.sup.2]
  Creatinine clearance             [micro]g/mL/min
  ALAD, 12 vs. 11                  [micro]g

Model 3
  Age                              [micro]g/year
  Female                           [micro]g
  Current smoker                   [micro]g
  Body mass index                  [micro]g/kg/[cm.sup.2]
  Creatinine clearance             [micro]g/mL/min
  ALAD, 12 vs. 11                  [micro]g
  ALAD x creatinine clearance(a)   [micro]g/mL/min

Independent                          [Beta]      [Beta]
variable                           coefficient     SE

Model 1
  Age                                  0.033     0.004
  Female                              -0.958     0.105
  Current smoker                       0.299     0.080
  Body mass index                      0.010     0.012
  Creatinine clearance                 0.006     0.001
  VDR, Bx vs. bb                       0.282     0.103

Model 2
  Age                                  0.034     0.004
  Female                              -0.946     0.106
  Current smoker                       0.310     0.081
  Body mass index                      0.013     0.012
  Creatinine clearance                 0.006     0.001
  ALAD, 12 vs. 11                      0.044     0.108

Model 3
  Age                                  0.034     0.004
  Female                              -0.937     0.106
  Current smoker                       0.307     0.081
  Body mass index                      0.012     0.012
  Creatinine clearance                 0.005     0.001
  ALAD, 12 vs. 11                     -0.890     0.458
  ALAD x creatinine clearance(a)       0.008     0.004

Independent                                    Model
variable                           p-Value   [r.sup.2]

Model 1                                        0.26
  Age                              < 0.001
  Female                           < 0.001
  Current smoker                   < 0.001
  Body mass index                    0.40
  Creatinine clearance             < 0.001
  VDR, Bx vs. bb                     0.006

Model 2                                        0.25
  Age                              < 0.001
  Female                           < 0.001
  Current smoker                   < 0.001
  Body mass index                    0.26
  Creatinine clearance             < 0.001
  ALAD, 12 vs. 11                    0.69

Model 3                                        0.26
  Age                              < 0.001
  Female                           < 0.001
  Current smoker                   < 0.001
  Body mass index                    0.29
  Creatinine clearance             < 0.001
  ALAD, 12 vs. 11                    0.05
  ALAD x creatinine clearance(a)     0.04

(a) DMSA-chelatable lead was log-transformed for these regressions because of departure from normality.

Predictors of tibia lead levels in lead-exposed subjects. After adjustment for age (linear and quadratic terms), sex, job duration, and body mass index, ALAD genotype was not associated with tibia lead levels (Table 6, model 1, p = 0.73), but VDR genotype was associated (Table 6, model 2, p = 0.03). On average, subjects with the VDR B allele had tibia lead levels that were 29%, or 6.4 [micro]g/g, higher than did subjects without the allele. The final regression models accounted for 15% of the variance in tibia lead levels.

Table 6. Linear regression modeling of tibia lead, Korean lead workers, 1997-1999.(a)

                     Units of
Independent          [Beta]                     [Beta]
variable             coefficient                coefficient

Model 1(b)
  Age                [micro]g/g/year                0.014
  [Age.sup.2]        [micro]g/g/[year.sup.2]        0.001
  Female             [micro]g/g                    -0.407
  Job duration       [micro]g/g/year                0.048
  Body mass index    [micro]g/g/kg/[cm.sup.2]       0.033
  ALAD (12 vs. 11)   [micro]g/g                     0.042

Model 2(b)
  Age                [micro]g/g/year                0.013
  [Age.sup.2]        [micro]g/g/[year.sup.2]        0.001
  Female             [micro]g/g                    -0.412
  Job duration       [micro]g/g/year                0.050
  Body mass index    [micro]g/g/kg/[cm.sup.2]       0.030
  VDR (Bx vs. bb)    [micro]g/g                     0.254

Independent          [Beta]
variable               SE     p-Value

Model 1(b)
  Age                0.005      0.002
  [Age.sup.2]        0.0003     0.003
  Female             0.104    < 0.001
  Job duration       0.007    < 0.001
  Body mass index    0.013      0.01
  ALAD (12 vs. 11)   0.122      0.73

Model 2(b)
  Age                0.005      0.006
  [Age.sup.2]        0.0003     0.002
  Female             0.103    < 0.001
  Job duration       0.007    < 0.001
  Body mass index    0.013      0.02
  VDR (Bx vs. bb)    0.117      0.03

(a) Tibia lead was log-transformed for these regressions because of departure from normality.

(b) Model [r.sup.2] = 0.15.

VDR genotype reportedly influences bone mineral density, so we examined the influence of VDR genotype on age- and sex-associated differences in tibia lead levels. There were no interactions observed between VDR genotype and age on tibia lead levels in all lead workers or in specific subgroups (i.e., males, females, older subjects), and no interaction between VDR genotype and sex on tibia lead levels.

Predictors of lead biomarkers in controls without occupational lead exposure. In linear regression models including only controls, none of these genetic associations with blood lead, DMSA-chelatable lead, or tibia lead levels were observed.

Discussion

The current study confirms past observations that the ALAD gene modifies the toxicokinetics of lead and also provides new evidence that the VDR gene does so as well. In fact, the influence of the VDR B allele on blood lead levels was larger than was the influence of the [ALAD.sup.2] allele. The mechanism by which these genes influence blood lead levels may differ. [ALAD.sup.2] and VDR B were associated with higher blood lead levels; however, only VDR B was associated with higher tibia lead levels (p = 0.03).

In adjusted analyses, subjects with the VDR B allele had significantly (p [is less than] 0.05) higher blood lead levels (on average, 4.2 [micro]g/dL), chelatable lead levels (on average, 37.3 [micro]g), and tibia lead levels (on average, 6.4 [micro]g/g) than did subjects with the VDb genotype. VDR genotype did not modify relations between such factors as age, sex, and renal function and any of the lead dose measures. In part, these observations may be explained by the greater intestinal absorption of lead, or greater uptake and subsequent release of lead from bone, in individuals with VDR B (13,27-32). Vitamin D, after binding to the VDR receptor, increases intestinal absorption of calcium and lead. The VDR B allele has been associated with lower bone mineral densities and higher tibia lead levels (13,15), but the mechanism underlying these observations is not currently known. Interpretation of the observation that lead workers with the VDR B allele have higher tibia lead levels than do workers with VDR bb is complicated by the fact that VDR genotype is likely to influence the content of both calcium and lead in bone and tibia lead concentration as measured by XRF is standardized to bone mineral content. Thus, higher tibia lead concentrations can be due to higher lead content, lower calcium content, or both.

In adjusted analyses, subjects with the [ALAD.sup.2] allele had higher blood lead levels (on average, 3.6 [micro]g/dL; p = 0.01) but no differences in tibia or chelatable lead levels compared to subjects without the allele. Creatinine clearance was an important predictor of chelatable lead levels and ALAD genotpe modified the relation between creatinine clearance and chelatable lead. Subjects with the [ALAD.sup.2] allele had larger increases in chelatable lead levels with increasing creatinine clearance than did subjects without the allele. We previously reported that subjects with [ALAD.sup.2] had lower DMSA-chelatable lead levels than did lead workers with [ALAD.sup.1] (6). It is important to note that in the previous study, DMSA was administered at 5 mg/kg, and mean DMSA-chelatable lead levels were approximately half of those in the current study, in which workers were administered 10 mg/kg DMSA. These data are consistent with earlier observations that the [ALAD.sup.2] allele increases erythrocytic binding of lead (3-12). This increase in intraerythrocytic lead may decrease the relative deposition of lead in critical target organs and thus protect against the toxicity of lead.

An unexpected observation was that ALAD and VDR genotypes were associated. These genes are located on different chromosomes. Lead workers with the ALAD1-1 genotype were much less likely to have the VDR bb genotype (OR = 0.29, p [is less than] 0.05). Although this observation has to be interpreted with caution because only three lead workers had the ALAD1-2 and VDR Bb or BB genotypes, the exact confidence interval did not include 1.0. The data also suggested that the association between the two genotypes in the controls without occupational lead exposure was different, in that controls with the ALAD1-1 genotype were more likely to have the VDR bb genotype (OR = 2.5), and the stratum-specific ORs were significantly different. This is not a stable estimate, however, due to small cell sizes, and requires confirmation. One study reported that differential selection in the lead industry may occur by ALAD genotype (3), and the current data are further evidence that genetic factors may influence the duration of work in the lead industry.

To date, data would suggest that the [ALAD.sup.1] allele is more likely to confer health risks associated with lead exposure. For the VDR genotype, data are insufficient to determine whether the polymorphisms are likely to modify health risks due to lead and which allele is the allele of risk. We speculate that the alleles for either ALAD or VDR that confer health risk should become less prevalent with increasing duration of occupational exposure to lead (3); this could occur, for example, if the at-risk alleles are associated with the development of acute symptoms that increase the probability of quitting jobs with lead exposure. We have no information in former lead workers on either symptoms or work duration by ALAD or VDR genotype, and serial blood lead measurements from the start of employment are not available for the majority of lead workers. These limitations weaken the inferences that we can make at this time. However, compared to the controls, it appears that lead workers have higher prevalences of both the [ALAD.sup.2] and VDR B alleles (10.3% vs. 8.1% and 11.4% vs. 8.9%, respectively). This observation would support the inference that the [ALAD.sup.2] allele is "protective," as is the VDR B allele, and that there may be selection by genotype among lead workers, but this speculation requires further study.

REFERENCES AND NOTES

(1.) Potluri VR, Astrin KH, Wetmur JK, Bishop DF, Desnick RJ. Human [Delta]-aminolevulate dehydratase: chromosomal localization to 9q34 by in situ hybridization. Hum Genet 76:236-239 (1987).

(2.) Battistuzzi G, Petrucci R, Silvagni L, Urbani FR, Caiola S. [Delta]-Aminolevulinate dehydratase: a new genetic polymorphism in man. Ann Hum Genet 45:223-229 (1981).

(3.) Schwartz BS, Lee B-K, Stewart W, Ahn K-D, Springer K, Kelsey K. Associations of [Delta]-aminolevulinic acid dehydratase genotype with plant, exposure duration, and blood lead and zinc protoporphyrin levels in Korean lead workers. Am J Epidemiol 142:738-745 (1995).

(4.) Wetmur JG, Lehnert G, Desnick RJ. The [Delta]-aminolevulinate dehydratase polymorphism: higher blood lead levels in lead workers and environmentally exposed children with the 1-2 and 2-2 isozymes. Environ Res 56:109-119 (1991).

(5.) Ziemsen B, Angerer J, Lehnert G, Benkmann HG, Goedde HW. Polymorphism of [Delta]-aminolevulinic acid dehydratase in lead-exposed workers. Int Arch Occup Environ Health 58:245-247 (1986).

(6.) Schwartz BS, Lee B-K, Stewart W, Ahn K-D, Kelsey KT. [Delta]-Aminolevulinic acid dehydratase genotype modifies 4-hour urinary lead excretion after oral administration of dimercaptosuccinic acid. Occup Environ Med 54:241-240 (1997).

(7.) Sithisarankul P, Schwartz BS, Lee B-K, Kelsey KT, Strickland PT. Aminolevulinic acid dehydratase genotype mediates plasma levels of the neurotoxin, 5-aminolevulinic acid, in lead-exposed workers. Am J Ind Med 32:15-20 (1997).

(8.) Smith CM, Wang X, Hu H, Kelsey KT. A polymorphism in the [Delta]-aminolevulinic acid dehydratase gene may modify the pharmacokinetics and toxicity of lead. Environ Health Perspect 103:248-253 (1995).

(9.) Fleming DEB, Chettle DR, Wetmur JG, Desnick RJ, Robin J-P, Boulay D, Richard NS, Gordon CL, Webber CE. Effect of the [Delta]-aminolevulinate dehydratase polymorphism on the accumulation of lead in bone and blood in lead smelter workers. Environ Res 77:49-61 (1998).

(10.) Alexander BH, Checkoway H, Costa-Mallen P, Faustman EM, Woods JS, Kelsey KT, van Netten C, Costa LG. Interaction of blood lead and [Delta]-aminolevulinic acid dehydratase genotype on markers of heme synthesis and sperm production in lead smelter workers. Environ Health Perspect 106:213-216 (1998).

(11.) Bergdahl IA, Gerhardsson L, Schutz A, Desnick RJ, Wetmur JG, Skerfving S. Delta-aminolevulinic acid dehydratase polymorphism: influence on lead levels and kidney function in humans. Arch Environ Health 52:91-96 (1997).

(12.) Bergdahl IA, Grubb A, Schutz A, Desnick RJ, Wetmur JG, Sassa S, Skerfving S. Lead binding to [Delta]-aminolevulinic acid dehydratase (ALAD) in human erythrocytes. Pharmacol Toxicol 81:153-158 (1997).

(13.) Schwartz BS, Stewart WF, Kelsey KT, Simon D, Park S, Links JM, Todd AC. Associations of tibia lead levels with BsmI polymorphisms in the vitamin D receptor in former organolead manufacturing workers. Environ Health Perspect 108:199-203 (2000).

(14.) Taymans SE, Pack S, Pak E, Orban Z, Barsony J, Zhuang Z, Stratakis CA. The human vitamin D receptor gene (VDR) is localized to region 12cen-q12 by fluorescent in situ hybridization and radiation hybrid mapping: genetic and physical VDR map. J Bone Miner Res 14:1163-1166 (1999).

(15.) Cooper GS, Umbach DM. Are vitamin D receptor polymorphisms associated with bone mineral density? A meta-analysis. J Bone Miner Res 11:1841-1849 (1996).

(16.) Schwartz BS, Lee B-K, Lee G-S, Stewart WF, Lee S-S, Hwang K-Y, Ahn K-D, Kim Y-B, Bolla KI, Simon D, Parsons PJ, Todd AC. Associations of blood lead, DMSA-chelatable lead, tibia lead, and job duration with neurobehavioral test scores in Korean lead workers. Am J Epidemiol (in press).

(17.) Thomas WJ, Collins TM. Comparison of venipuncture blood counts with microcapillary measurements in screening for anemia in one-year-old infant. J Pediatr 101:32-35 (1982).

(18.) Blumberg WE, Eisinger J, Lamola AA, Zuckerman DM. Zinc protoporphyrin level in blood determination by a portable hematofluorometer: a screening device for lead poisoning. J Lab Clin Med 89:712-723 (1977).

(19.) Heinegard D, Tiderstrom G. Determination of serum creatinine by a direct colorimetric method. Clin Chem Acta 43:305 (1973).

(20.) Kneip TJ, Crable JV. Methods for biological monitoring: a manual for assessing human exposure to hazardous substances. Washington, DC:American Public Health Association, 1988;199-201.

(21.) Todd AC, McNeill FE. In Vivo Measurements of Lead in Bone Using a Cd Spot Source. In: Human Body Composition Studies. New York:Plenum Press, 1993;299-302.

(22.) Todd AC. Contamination of in vivo bone-lead measurements. Phys Meal Biol 45:229-240 (2000).

(23.) Todd AC, Lee B-K, Lee G-S, Schwartz BS. Predictors of blood lead, tibia lead, and DMSA-chelatable lead in 802 Korean lead workers. Occup Environ Med, submitted.

(24.) Kim R, Aro A, Rotnitzky A, Amarasiriwardena C, Hu H. K X-ray fluorescence measurements of bone lead concentration: the analysis of low-level data. Phys Med Biol 40:1475-1485 (1995).

(25.) Lee B-K, Schwartz BS, Stewart W, Ahn K-D. Urinary lead excretion after DMSA and EDTA: evidence for differential access to lead storage sites. Occup Environ Med 52:13-19 (1995).

(26.) Parsons PJ, Slavin W. Electrothermal atomization atomic absorption spectrometry for the determination of lead in urine: results of an interlaboratory study. Spectrochim Acta Part B 54:853-864 (1999).

(27.) Fullmer CS. Intestinal lead and calcium absorption: effect of 1,25-dihydroxycholecalciferol and lead status. Proc Soc Exp Biol Med 194:258-264 (1990).

(28.) Fullmer CS. Intestinal interactions of lead and calcium. Neurotoxicology 13:799-807 (1992).

(29.) Smith CM, DeLuca HF, Tanaka Y, Mahaffey KR. Effect of lead ingestion on functions of vitamin D and its metabolites. J Nutr 111:1321-1329 (1981).

(30.) Morrison NA, Yeoman R, Kelly PJ, Eisman JA. Contribution of trans-acting factor alleles to normal physiological variability: vitamin D receptor gene polymorphism and osteocalcin. Proc Natl Acad Sci USA 89:6665-6669 (1992).

(31.) Barger-Lux MJ, Heaney RP, Hayes J, DeLuca JF, Johnson ML, Gong G. Vitamin D receptor gene polymorphism, bone mass, body size, and vitamin D receptor density. Calcif Tissue Int 57:161-162 (1995).

(32.) Morrison NA, Qi JC, Tokita A, Kelly PJ, Crofts L, Nguyen TV, Sambrook PN, Eisman JA. Prediction of bone density from vitamin D receptor alleles. Nature 367:284-287 (1994).

Brian S. Schwartz,(1,2,3) Byung-Kook Lee,(4) Gap-Soo Lee,(4) Walter F. Stewart,(1,3) David Simon,(3) Karl Kelsey,(5) and Andrew C. Todd(6)

(1) Department of Environmental Health Sciences, Johns Hopkins School of Hygiene and Public Health, Baltimore, Maryland, USA; (2) Department of Medicine, Johns Hopkins School of Medicine, Baltimore, Maryland, USA; (3) Department of Epidemiology, Johns Hopkins School of Hygiene and Public Health, Baltimore, Maryland, USA; (4) Institute of Industrial Medicine, Soonchunhyang University, Chonan, Korea; (5) Department of Cancer Cell Biology, Harvard School of Public Health, Boston, Massachusetts, USA; (6) Department of Community and Preventive Medicine, Mount Sinai Medical Center, New York, New York, USA.

Address correspondence to B.S. Schwartz, Division of Occupational and Environmental Health, Johns Hopkins School of Hygiene and Public Health, Room 7041, 615 Wolfe Street, Baltimore, MD 21205 USA. Telephone: [phone]. Fax: [phone]. E-mail: [e-mail]

We thank P.J. Parsons for performing the urine lead measurements and Y-B. Kim, K-Y. Hwang, S-S. Lee, and K-D. Ahn for assisting in data collection in Korea.

This research was supported by grant ES07198 (to B.S.S.) from the National Institute of Environmental Health Sciences.

Received 17 April 2000; accepted 30 May 2000.