Skip directly to: content | left navigation | search

Induction of P-450 in workers exposed to dioxin

    Halperin, William; Kalow, Werner; Sweeney, Marie Haring; Tang, Bing K.; Fingerhut, Marilyn;
    Timpkins, Brent; Wille, Kathleen

    National Institute for Occupational Safety and Health, 4676 Columbia Parkway, Cincinnati, Ohio, USA (W Halperin, MH Sweeney, M Fingerhut, B Tompkins, K Wille).
    Department of Pharmacology, University of Toronto (W Kalow, B K Tang).
    Correspondence to: Dr William Halperin, Division of Surveillance, Hazard Evaluations, and Field Studies, National Institute for Occupational Safety and Health, Centers for Disease Control and Prevention, 4676 Columbia Parkway, Cincinnati, OH 45226, USA.
    Accepted 15 September 1994.

    Abstract

    Objectives: To examine the effects of occupational exposure to substances contaminated with 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) on cytochrome P-4501A2 activity in a cross sectional medical survey.

    Methods: The exposed workers had been employed at two chemical plants > 15 years earlier in the manufacture of 2,4, 5-trichlorophenol and its derivatives. The control group consisted of people with no occupational exposure to phenoxy herbicides and who lived within the communities of the exposed workers. A total of 58 workers and 125 unexposed controls participated in the analysis. Cytochrome P-450 activity was assessed with a test that measures caffeine metabolites in the urine. A ratio of metabolites of caffeine (CMR) constituted a measure of P-4501A2 activity.

    Results: Compared with the control group in multivariate logistic regression, raised non-significant associations were found for three of four categories of TCDD in exposed workers (TCDD < 20 pg/g, odds ratio (OR) 1.7, 95% confidence interval (95% CI) 0.6 to 5.0, TCDD 20–66, OR 0.3, 95% CI 0.0 to 1.7; TCDD 67–147, OR 2.3, 95% CI 0.6 to 8.8; TCDD greater than or equal to 148, OR 3.1, 95% CI 0.8 to 12.5). We found a strongly significant association of CMR and urinary cotinine, a measure of smoking, and urinary free ethanol. We found weak non-significant associations between P-4501A2 activity and increased serum TCDD among workers.

    Conclusions: The absence of an association between serum TCDD and cytochrome P-4501A2 may be due to the size of the study, insensitivity of the CMR to assess cytochrome P-4501A2 activity, or inadequate levels of exposure, although these were among the highest in human groups tested.

    (Occup Environ Med 1995;52:86–91)

    Keywords: dioxin; cytochrome P-450.

    Insights into the genetic basis for the toxicity of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD began with the discovery that TCDD was a potent inducer of cytochrome P-450 and that this induction was genetically controlled. This led to the recognition of a cellular macromolecule, the Ah receptor, which specifically recognises and binds with TCDD. After processing and transport across the nuclear membrane TCDD binds to these specific sites on the DNA where it controls synthesis of cellular proteins that include cytochrome P-450 enzymes. In animals TCDD induces synthesis of two members of the P-4501A family of cytochromes P-450, P-4501A1 (aryl hydrocarbon hydroxylase (AHH)) and P-4501A2 (phenacetin-O-deethylase). There is evidence for the existence of the Ah system in humans, which suggests that the mechanism involved in TCDD toxicity may be similar in humans and other susceptible species. To describe the functional relation between exposure to TCDD and enzyme induction, and to determine the sensitivity of the system, we studied a group of controls and a group of workers with known serum concentrations of TCDD from former exposure to TCDD in the manufacture of herbicides and their precursors.

    Grant et al described the conceptual basis for examination of profiles of the metabolites of caffeine as indicators of P-450 induction. Caffeine (1,3,7 trimethylxanthine (137X)) is metabolised to 1,7 dimethylxanthine (17X) by P-4501A2The 17X is then oxidised to 1,7 dimethyluric acid (17U) by a P-450 that is not substantially inducible, or is further 7-demethylated by P-4501A2 to three other metabolites, 5-acetylamino-6-formylamino-3-methyluracil (AFMU), 1-methylxanthine (1X), and a derivative of 1X, 1-methyl-uric acid (IU). The ratio of (AFMU + 1X + 1U)/17U in urine, caffeine metabolite ratio (CMR), has been proposed as an indicator of cytochrome P-450 activity inducible by polyaromatic hydrocarbons. The CMR has previously been shown to be increased after environmental exposure to polybrominated biphenyls (PBBs and to cigarette smokebut not by exposure to polyaromatic hydrocarbons in a foundry. The CMR seems to be minimally related to sex and is significantly reduced by pregnancy and exogenous oestrogens; physical training increases the CMR as does ingestion of broccoli

    Methods

    STUDY POPULATION

    The study population of exposed workers was recruited from current and former employees of two plants of the 12 included in the NIOSH mortality study of workers exposed to chemicals contaminated with TCDD. At the two plants, TCDD exposure resulted from the manufacture of 2,4,5-trichlorophenol (TCP) and its derivatives. Exposure to TCDD occurred between 1951 and 1969 at the New Jersey plant, and for about two years between 1968 and 1972 at the Missouri plant. A cross sectional medical study of workers from the two plants was conducted in 1987 and 1988. .

    A total of 586 workers, 490 from the New Jersey factory and 96 from the Missouri factory, were eligible for study. Of the 586 workers, 143 were dead (24%), and 43 (7.3%) were not traced, leaving 400 workers (68%) who were invited to participate in the study. A total of 357 of the eligible workers (89%) completed an occupational history; of those, 281 (70% of those eligible) also participated in the medical examination. A complete description of the study has been previously reported.Twenty office workers with no appreciable occupational exposure were excluded from this analysis. Contractual arrangements for collection of urine samples for caffeine metabolism were concluded only after the study was well under way. Urine samples for caffeine metabolites and serum for measurement of TCDD were collected from 58 exposed workers. This represents 100% of the workers who were examined after the start of collection of urine for the analysis of caffeine metabolites.

    Two hundred and sixty controls from the same neighbourhood, with no history of occupations in factories that manufactured chemicals contaminated with TCDD, participated in the medical examinations out of 325 controls who were interviewed. Age (within five years), race, and sex matched controls were identified by a door to door screening of each worker's residential census record. Seventy five per cent of the controls represented the first, second, or third person identified as a potential candidate.

    Of the 260 controls, serum TCDD was measured in only a sample of 79 people because only one of the sample controls had concentrations > 20 pg/g and assays were expensive. Of the 260 controls, all 125 examined after contractual arrangements were completed were tested for CMR. Of the 125 controls tested for CMR, 24 had also been tested for TCDD. For epidemiological analysis, the 101 controls not tested for serum TCDD were assigned the median serum TCDD concentration found in the controls tested, 6.08 pg/g. This allowed the 101 controls not tested for TCDD to stay in the analysis and provide information on the association of CMR with other potential confounders, such as cotinine. Also, we estimated TCDD in samples where serum concentrations fell below the limit of detection. The estimated concentration was calculated by dividing the limit of detection of each sample by radical 2.

    LABORATORY METHODS

    We have reported the methods of collection and preparation, analytical technique, and quality control standards of the laboratory analyses for clinical chemistry and TCDD analyses. For analysis of CMR, participants were asked to consume no less than one cup of coffee or tea or one can of soft drink and no more than six cups or cans on the day of urine collection. A nine hour urine collection began at 10 00 pm. At 7 00 am the collected urine was acidified to pH 3.5 and stored at -20 degrees Celsius. The set of frozen tubes was carried on dry ice by hand to the University of Toronto.

    High performance liquid chromatography was used to measure urates and 1X. The measurements of AFMU were performed according to Tang et al. Cotinine was determined by the same procedure without modification of the TSK-20 column.

    STATISTICAL ANALYSIS

    For purposes of data exploration and ease of interpretation we first used simple stratified analyses. Given the known association of smoking cigarettes with CMR, smoking was chosen as one variable, worker v control was the other. Participants with urinary cotinine concentrations > 0.0 were considered to be smokers. Variables considered included serum TCDD, age at the time of examination, alcohol-years (average number of alcoholic drinks consumed a day multiplied by the number of years alcohol was consumed), concentration of D-glucaric acid, pack-years (average number of cigarettes smoked a day divided by 20 multiplied by the number of years cigarettes were smoked), cotinine (/g creatinine), weight (kg), race, self reported history of smoking (smoked within the past year), and self reported alcohol consumption (drank alcohol within the past year). Statistical comparisons were done with Student's t test for continuous variables and chi squared test for categorical variables. All analyses were performed with SAS software. [24].

    We used multivariate linear regression analyses to examine the association of CMR with serum TCDD and controlled for confounding variables. Models were constructed iteratively by the addition of predictor variables based upon biological reasoning. The first model explored the expected relation between CMR, the dependent variable, and serum TCDD; urinary cotinine was then added to the model. Subsequent models included variables that were predicted from the start to have some contribution to CMR. In these analyses, the distributions of CMR, serum TCDD, cotinine, etc were examined for normality. Log transformation of both the independent and dependent variables did not improve the normality of the residuals—that is, the difference between the observed and predicted values for each predictor—or the explanatory power of the regression equations (R2). Non-normalised variables are presented for ease of interpretation.

    To model the relation of CMR with concentration of TCDD and to assess the contribution of other potential effect modifiers and confounders, such as cotinine, we conducted logistic regression analyses. For the logistic regression, CMR was dichotomised at the midpoint value (6.11) of the control population. Serum TCDD was included as a five level categorical variable and forced into all models regardless of significance. The controls consisted of 125 unexposed people, all of whom had TCDD concentrations < 20 pg/g. The TCDD for workers was stratified into four levels (< 20, 20–66, 67–147, and greater than or equal to 148 pg/g). The highest three categories of TCDD were determined by dividing the number of workers with concentrations of serum TCDD of greater than or equal to 20 pg/g into three groups nearly equal in size. After this, given the known association of cigarette smoking and enzyme induction, cotinine was entered as a categorical variable, but as the models did not converge, cotinine was entered into a substitute model as a continuous variable.

    To further delineate the relation of TCDD with CMR, other potential confounders and effect modifiers were examined by adding them one at a time to a base model containing TCDD, cotinine, and free ethanol. The additional variables were D-glucaric acid (a measure of phenobarbital induced enzymes), age at examination, cups of coffee consumed, weight, years since last exposed, sex, race, and specific cohort (New Jersey or Missouri). Values for D-glucaric acid and free ethanol were divided into quartiles and age and cups of coffee dichotomised based on the distribution in the control population. Tests for trends were conducted by entering TCDD as a continuous variable. An indicator variable was established to reflect whether TCDD concentrations in the controls were measured or estimated.

    Results

    This study population consisted of 58 exposed workers and 125 controls from the same neighbourhood. The age of the study population ranged from 31 to 77 years. Table 1 shows the distribution of smoking, drinking, race, age, weight, and serum TCDD, by category of serum TCDD. The upper three groups had TCDD concentrations that were substantially higher than the controls and other unexposed populations that have been examined and live in industrialised countries. As expected, mean TCDD was substantially higher in workers than in controls (157 v 6.7 pg/g; P < 0.001). Workers also weighed more than controls (84.9 v 80.1 kg; P < 0.05).

    Table 1. Characteristics of the study population by serum 2,3,7,8-TCDD category

    As smoking was a strong inducer of P-4501A2 and could mask an effect due to TCDD, we compared workers and controls by smoking (Table 2). Smokers, regardless of exposure, drank more coffee. The difference was significant only in workers. The CMR was also significantly higher in smokers among workers (10.9 v 6.8; P < 0.05) and controls (11.6 v 5.9; P < 0.001). The association of smoking and CMR served to validate the performance of the CMR assay in this study. When smokers among the workers and controls were compared, there were significant differences in weight.

    Table 2. Comparison of workers and controls by smoking

    The 24 controls who were tested for TCDD compared with the 101 who were estimated for TCDD (Table 3) had significantly higher concentrations of free ethanol and lower concentrations of D-glucaric acid (P < 0.001). As the controls for this test were chosen consecutively, the difference may represent either a chance occurrence or a temporal pattern in the characteristics of controls associated with their date of examination.

    Table 3. Mean (SD) concentrations for controls with measured and estimated TCDD

    No association of CMR with TCDD was evident from a scattergram (Figure 1) of only these variables.

    In multivariate linear regression analyses (Table 4), there was almost no relation of CMR with serum TCDD (beta coefficient 0.00012, P = 0.5). We found a significant positive association between urinary cotinine, an indication of current cigarette smoking and CMR (beta coefficient 0.194, P = 0.0001). Free ethanol, an indication of recent alcohol consumption, was positively associated with CMR (beta coefficient 0.0012, P = 0.003). Weight and CMR were negatively associated (beta coefficient -0.043, P = 0.03). A model that contained cotinine and serum TCDD predicted 28% of the variation in CMR, compared with only 0.25% for TCDD alone. The intercept (6.2 in model 2), a measure of the non-induced or constitutive enzyme activity in non-smokers was consistent with studies in other populations and in other laboratories.

    image

    Figure 1. Caffeine metabolite ratio (CMR) v concentrations of the dioxin TCDD in exposed workers (closed circle) and controls (open circle).

    Table 4. Predictors of CMR in models of linear regression analyses

    (Table 5) shows the results of analysis by logistic regression. In model 1, four groups of exposed workers were compared with the controls. Although odds ratios ORs were high for three of the four exposed groups, confidence intervals (CIs) of the ORs were broad and none exceeded expectancy, an OR of 1. With the known relation of CMR to smoking, in model 2 cotinine was added as a predictive variable. The OR (1.2; 95% CI 1.1 to 1.3) indicates a 20% increase in CMR for each unit rise in cotinine. The addition of the cotinine to the model did not substantially change the ORs for concentrations of TCDD, which suggests little confounding between these two variables. Urinary free ethanol was entered into model 3 because of variations in the concentration in free ethanol found in the stratified analyses. The OR increased with concentration of free ethanol, with the highest quartile significantly associated with CMR (OR 4.9; 95% CI 1.7 to 13.9). Other variables were then added one at a time (models 4–11; data not shown) to a model that included all values of serum TCDD, cotinine, and the highest concentration of free ethanol. These included D-glucaric acid, age at time of the examination, cups of coffee consumed before the examination, weight, years since last exposure to TCDD, sex, race, and whether a member of the New Jersey or Missouri cohort. None of these additional variables were significant predictors of CMR. In no model tested was the indicator variable for measured v estimated concentration of TCDD a significant predictor of CMR. To assess a trend in association between serum TCDD and CMR, TCDD was entered as a continuous variable in a model that included cotinine as a continuous variable, and values of urinary free ethanol. No trend was found; the OR for TCDD was 1.0010 (95% CI; 0.999 to 1.003).

    Table 5. Predictors of CMR* in logistic regression analysis model

    Discussion

    This study showed essentially no relation between serum TCDD and induction of P-450 enzyme as measured by the CMR. In comparison, and consistent with other reports, [25] a strong association was evident between induction of P-4501A2, as measured by CMR, and cigarette smoking as measured by urinary cotinine. The only indication of an association between CMR and TCDD was faint and was found in logistic regression. In logistic regression, for three of four categories of exposed workers defined by increasing concentrations of serum TCDD, the OR was high. The increases were not significant.

    There are several explanations for these findings. Size was a main limitation of this study. When the study population is divided for analysis into groups by exposure, there are small numbers of participants in the highly exposed groups. Although these chemical workers currently have among the highest serum concentrations of TCDD that have been found in human populations, and they have carried these burdens for > 15 years since last exposure, their numbers may be too small in this study to accurately reflect the association between TCDD and CMR. The limited power to detect an association between TCDD and CMR is particularly worrying because of the increased ORs in the workers of the two highest categories of exposure to TCDD (Table 5).

    This population showed a strong relation of CMR to cigarette smoking, to which TCDD, as an inducer, can be compared. The strength of association between CMR and cigarette smoking (R2 0.269) is lower than previously reported (R2 0.81), which may reflect the older age of our study population.

    Other biases that may have obscured a strong relation between CMR and serum TCDD include measures of exposure and measures of response. The measure of exposure among the workers was their serum TCDD. In our population there was a strong correlation (R 0.7, P < 0.001) between serum TCDD and duration of employment in an exposed job; hence substantial error in measurement of serum TCDD is unlikely. The measure of response in this population was the CMR, which has been used in several earlier studies that investigated the relation of enzyme induction and cigarette smoking,and exposure to polyaromatic hydrocarbons in foundries diet, exercise, and oestrogrens. These studies showed strong and consistent relations between smoking cigarettes and CMR, as was evident in our analysis. It is possible that the CMR, although an extremely practical test and clearly sufficiently sensitive to detect the strong relation of cigarette smoking with enzyme induction, may not be sufficiently sensitive and specific in our study population to measure P-4501A2 induction associated with TCDD. There has been substantial discussion about the enzymes actually reflected by the caffeine urine test compared with the caffeine breath test. There has been further discussion about the most appropriate caffeine metabolites to measure in the urine test. future studies might include comparisons of the sensitivity of alternative methods to assess cytochrome P-4501A2 induction.

    The stronger relation between cigarette smoking than TCDD on CMR suggests two conclusions.

    Cigarette smoking may be a stronger inducer of P-4501A2 than is TCDD in the highly exposed participants in this study, or CMR may measure another enzyme that is stimulated by cigarette smoking and that is distinct from P-4501A2. This possibility is very unlikely, based on in vitro studies with human P-450s derived from DNA which seemed to account for most or all of the caffeine metabolising capacity of microsomes of the human liver.

    Another possible explanation for the weak association between serum TCDD and CMR is that the dose of TCDD in this population, although among the highest found in humans, may not be sufficient in comparison with experimental doses given to animals to elicit an induction in P-450. It has been suggested that humans are among the least sensitive species for induction of hepatic enzyme.

    A final alternative is even more speculative. This study differs from experiments on laboratory animals in its time frame. The workers in this population had their last occupational exposure to TCDD 15 to 37 years before this study was performed and because of the long half life of TCDD, can be considered as having been continuously exposed during employment and subsequently. In comparison experimental studies usually measure enzyme induction comparatively soon after TCDD exposure. Might it be possible that TCDD induction of P-4501A2, in parallel to well understood hormones, is self regulating? Perhaps the induction of P-4501A2 leads to increasing unresponsiveness in vivo to the continued potency of TCDD as an inducer. Some evidence against this hypothesis is the lack of negative association between CMR and years since last exposure.

    In summary, this study showed the inability of serum TCDD in comparison with cigarette smoking, to raise CMR. Numerous methodological and mechanistic explanations are possible. These include questions about the sensitivity of the CMR, the responsiveness in humans of P-4501A2 induction by TCDD, the limited size of this study, or the inadequate exposure in this population to elicit a response, even in a group of the most highly occupationally exposed people.

    This study was funded, in part, by the Agency for Toxic Substances and Disease Registry (ATSDR).

    REFERENCES

      Silbergeld D, Gasiewicz T. Dioxins and the Ah receptor. Am J Ind Med 1989;16:455–74. [Medline Link] [Context Link]

      Poland A, Glover E. Chlorinated dibenzo-p-dioxins: potent inducers of delta-aminolevulinic acid synthetase and aryl hydrocarbon hydroxylase. II. A study of the structure-activity relationships. Mol Pharmacol 1973;9:736–47. [Medline Link] [Context Link]

      Nebert DW, Goujon FM, Gielen JE. Arlyhydrocarbon hydroxlase induction by polycyclic hydrocarbons: simple autosomal dominant trait in the mouse. Nature 1972;236:6836–42. [Context Link]

      Okey AB. Enzyme induction in the cytochrome P-450 system. Pharmacol Ther 1990;45:241–98. [Medline Link] [Context Link]

      Tritscher A, Clark, G, McCoy C, Greenlee W, Goldstein J, Lucier G. Dose response relationships for chronic exposure to 2,3,7,8-tetrachlorodibenzo-p-dioxin TCDD) in a rat tumor promotion model:2. Quantification and immunolocalization of cytochromes P450c1A1) and P450dlA2) in the liver [abstract]. Dioxin 91. 11th International Symposium on Chlorinated Dioxins and Related Compounds. Research Triangle Park, NC: 1991, 148. [Context Link]

      Kalow W. Genetic Variation in the human hepatic cytochrome p-450 system. Eur J Clin Pharmacol 1987;31:633–41. [Medline Link] [Context Link]

      Poland A, Knutson J. 2,3,7,8-tetrachlorodibenzo-p-dioxin and related halogenated aromatic hydrocarbons: examination of the mechanism of toxicity. Annual Review of Pharmacology and Toxicology 1982;22:517–54. [Context Link]

      Okey A, Roberts E, Harper P, Denison M. Induction of drug metabolizing enzymes: mechanism and consequences. Clin Biochem 1986;19:132–41. [Medline Link] [Context Link]

      Okey A, Riddick D, Harper P. The Ah receptor: mediator of the toxicity of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and related compounds. Toxicol Lett 1995;(in press). [Context Link]

      Grant DM, Tang BK, Kalow W. A simple test for acetylator phenotype using caffeine. Br J Clin Pharmacol 1984;17:459-66. [Medline Link] [Context Link]

      Vistisen K, Poulsen H, Loft S. Foreign compound metabolism capacity in man measured from metabolites of dietary caffeine. Carcinogenesis 1992;13:1561–8. [Medline Link] [Context Link]

      Butler MA, Iwasaki M, Guengerick FP, Kadlubar FF. Human cytochrome P-450pa P4501A2), phenacetin-O-deethylase, is primarily responsible for the hepatic d-demethylation of caffeine and N-oxidation of carcinogenic arylamines. Proc Natl Acad Sci USA 1989;86:7696–700. [Medline Link] [Context Link]

      Ikeya K, Jaiswal AK, Owens RA, Jones JE, Nebert DW, Kimura S. Human CYPIA2: sequence, gene structure, comparison with the mouse and rat homologous gene, and differences in liver mRNA expression. Mol Endocrinol 1989;9:1399–1408. [Medline Link] [Context Link]

      Campbell ME, Grant DM, Tadanobu I, Kalow W. Biotransformation of caffeine, parazanthine, theophylline, and theobromine by polycyclic aromatic hydrocarbon-inducible cytochromes) P-450 in human liver microsomes. Drug Metabolism and Disposition 1987;15:237–49. [Context Link]

      Lambert G, Schoeller D, Humphrey H, Kotake A, Lietz H, Campbell M, et al. The caffeine breath test and caffeine urinary metabolite ratios in the Michigan cohort exposed to polybrominated biphenyls: a preliminary study. Environ Health Perspect 1990;89:175–81. [Medline Link] [Context Link]

      Kalow W, Tang BK. Use of caffeine metabolite ratios to explore CYP1A2 and xanthine oxidase activities. Clin Pharmacol Ther 1991;50:508–19. [Medline Link] [Context Link]

      Sherson D, Sigagaard T, Overgaard E, Loft S, Poulsen HE, Jongeneelen FJ. Interaction of smoking, uptake of polycyclic aromatic hydrocarbons, and cytochrome P4501A2 activity among foundry workers. Br J Ind Med 1992;49:197–202. [Medline Link] [Context Link]

      Fingerhut MA, Halperin WE, Marlow D, Piacitelli L, Honchar P, Sweeney M, et al. Mortality among U S workers employed in the production of chemicals contaminated with 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). N Engl J Med 1991;324:212–8. [Medline Link] [Context Link]

      Sweeney M, Fingerhut M, Patterson D, Connally B, Piacitelli L, Morris J, et al. Comparison of serum levels of 2,3,7,8-tetrachlorodibenzo-p-dioxin 2,3,7,8-TCDD) in TCP production workers and in an unexposed comparison group. Chemosphere 1990;20:993–1000. [Context Link]

      Hornung RW, Reed LD. Estimation of average concentration in the presence of nondetectible values. App Occup Environ Hyg 1990;5:46–51. [Context Link]

      Patterson DG, Fingerhut MA, Roberts DW, Sweeney MH, Marlow D, Andrews J, Halperin WE. Levels of polychlorinated dibenzo-p-dioxin PCDD') and dibenzofurans DCDF's) in workers exposed to 2,3,7,8-tetrachlorobenzo-p-dioxin. Am J Ind Med 1989;16:135–46. [Medline Link] [Context Link]

      Campbell ME, Spielberg SP, Kalow W. A urinary metabolite ratio that reflects systemic caffeine clearance. Clin Pharmacol Ther 1987;42:157–65. [Medline Link] [Context Link]

      Tang BK, Kadar D, Kalow W. An alternate test for acetylator phenotyping with caffeine. Clin Pharmacol Ther 1987;42:509–13. [Medline Link] [Context Link]

      SAS, Version 6.03 ed. Cary, NC: SAS Institute, 1987. [Context Link]