A review of mammalian carcinogenicity study design and potential effects of alternate test procedures on the safety evaluation of food ingredients.

Extensive experience in conducting long term cancer bioassays has been gained over the past 50 years of animal testing on drugs, pesticides, industrial chemicals, food additives and consumer products. Testing protocols for the conduct of carcinogenicity studies in rodents have been developed in Guidelines promulgated by regulatory agencies, including the US EPA (Environmental Protection Agency), the US FDA (Food and Drug Administration), the OECD (Organization for Economic Co-operation and Development) for the EU member states and the MAFF (Ministries of Agriculture, Forestries and Fisheries) and MHW (Ministry of Health and Welfare) in Japan. The basis of critical elements of the study design that lead to an accepted identification of the carcinogenic hazard of substances in food and beverages is the focus of this review. The approaches used by entities well-known for carcinogenicity testing and/or guideline development are discussed. Particular focus is placed on comparison of testing programs used by the US National Toxicology Program (NTP) and advocated in OECD guidelines to the testing programs of the European Ramazzini Foundation (ERF), an organization with numerous published carcinogenicity studies. This focus allows for a good comparison of differences in approaches to carcinogenicity testing and allows for a critical consideration of elements important to appropriate carcinogenicity study designs and practices. OECD protocols serve as good standard models for carcinogenicity testing protocol design. Additionally, the detailed design of any protocol should include attention to the rationale for inclusion of particular elements, including the impact of those elements on study interpretations. Appropriate interpretation of study results is dependent on rigorous evaluation of the study design and conduct, including differences from standard practices. Important considerations are differences in the strain of animal used, diet and housing practices, rigorousness of test procedures, dose selection, histopathology procedures, application of historical control data, statistical evaluations and whether statistical extrapolations are supported by, or are beyond the limits of, the data generated. Without due consideration, there can be result conflicting data interpretations and uncertainty about the relevance of a study's results to human risk. This paper discusses the critical elements of rodent (rat) carcinogenicity studies, particularly with respect to the study of food ingredients. It also highlights study practices and procedures that can detract from the appropriate evaluation of human relevance of results, indicating the importance of adherence to international consensus protocols, such as those detailed by OECD.

Association between cage shelf level and spontaneous and induced neoplasms in mice.

Data from a chronic feeding study with 2-acetylaminofluorene [(2-AAF) CAS: 53-96-3; N-fluoren-2-yl-acetamide] done on 20,880 female BALB/c mice were analyzed for associations between cage shelf level and occurrence of induced and spontaneous neoplasms. Each cage was maintained on a rack at a given shelf level throughout the experiment, allowing analysis of data by shelf level. Differences in the crude incidence of 2-AAF-induced liver and bladder neoplasms appeared to be shelf related, but these differences were small and disappeared when shelf-level analyses of time-to-tumor onset distributions were performed. There was evidence for a shelf-level influence on 5 of the 6 major spontaneous neoplasms noted. Time to onset of uterine polyps and reticulum cell sarcomas was significantly delayed on the top shelf of five of six animal rooms. Also, there was a significant delay in onset of lymphomas, adrenocortical adenomas, and lung alveolar cell tumors on the top shelf when data were combined from all six animal rooms, but these delays on the top shelf were significant in no more than two of six animal rooms when rooms were analyzed separately. There was no indication of any shelf-level influence on the development of harderian gland adenomas. In conclusion, shelf level is an environmental factor that should be considered in the design and analysis of carcinogenesis studies.

Primary harderian gland neoplasms in mice.

Five hundred and seventy-four spontaneous harderian gland neoplasms were observed in 3,123 male and 9,024 female untreated experimental mice of various ages and of 6 inbred strains and substrains and 3 hybrid stocks. Adenomas occurred in 86 males and 458 females, and adenocarcinomas were present in 3 males and 27 females. The highest tumor incidence occurred in mice after 24 months of age; no tumors were observed in any mice prior to 6 months of age. With the exception of 1 age group, both sexes of the BALB/c mice had the highest incidence, and also this strain accounted for 463 of the neoplasms. The adenomas were classified into one of four different histologic types: papillary, cystic papillary, acinar, and cystic. The papillary type was the most frequently observed in both sexes. The cystic papillary type was the second most frequently observed and was the type most frequently associated with gross lesions. This type may represent a developmental variant of the papillary type. The incidence of the acinar type was similar in both sexes. The cystic type had the lowest incidence. All except 3 adenocarcinomas appeared to evolve from adenomas, with evolvement occurring primarily in the acinar type. Five of the adenocarcinomas either metastasized to the lung or invaded periorbital tissues. The remaining 25 were confined to the harderian gland.

Influence of cause-of-death assignment on time-to-tumor analyses in animal carcinogenesis studies.

The importance of cause-of-death determination in an animal carcinogenesis study with respect to estimation of time-to-tumor distributions of internally occurring (occult) tumors is discussed. A nontechnical description of time-to-tumor estimation is presented. The information obtained from time-to-tumor estimation when cause-of-death designation was used is illustrated for liver tumors in female mice of the inbred strain BALB/cStCrlfC3Hf/Nctr from the ED01 study with N-2-fluorenylacetamide done at the National Center for Toxicological Research. A time-to-tumor analysis of reticulum cell sarcoma data from the same study has provided insight into some difficulties involved in routine case-by-case determination of cause of death. A more flexible system for assigning of cause of death to dead animals and cause of morbidity to moribund animals is described as a way to improve cause-of-death assignment.

Correlation of sister chromatid exchange frequency with mutation induction and cell survival in Chinese hamster ovary cells exposed to N-oxidized arylamines.

In this study the relationships between sister chromatid exchange (SCE) frequency, mutation induction at the hypoxanthineguanine phosphoribosyl transferase locus, and cell survival were established in Chinese hamster ovary cells exposed to one of the three N-oxidized arylamines. The toxicants used were N-hydroxy-2-aminofluorene, N-hydroxy-N'-acetylbenzidine, and 1-nitrosopyrene. Mutation induction as measured by resistance to 6-thioguanine related poorly to reduced cloning ability and SCE frequency. Although SCEs were formed at concentrations of toxicants which produced no measurable reduction in cell survival, a strong correlation was observed between these two biological responses. This relationship was strengthened when the results obtained in this study were combined with those from a previous study which examined the relationship between SCE induction and cell survival in Chinese hamster ovary cells exposed to simple alkylating agents. These results support the contention that a common step is involved in the induction of SCE and cellular toxicity in Chinese hamster ovary cells exposed to certain classes of chemical toxicants.

Classification ensembles for unbalanced class sizes in predictive toxicology.

This paper investigates the effects of the ratio of positive-to-negative samples on the sensitivity, specificity, and concordance. When the class sizes in the training samples are not equal, the classification rule derived will favor the majority class and result in a low sensitivity on the minority class prediction. We propose an ensemble classification approach to adjust for differential class sizes in a binary classifier system. An ensemble classifier consists of a set of base classifiers; its prediction rule is based on a summary measure of individual classifications by the base classifiers. Two re-sampling methods, augmentation and abatement, are proposed to generate different bootstrap samples of equal class size to build the base classifiers. The augmentation method balances the two class sizes by bootstrapping additional samples from the minority class, whereas the abatement method balances the two class sizes by sampling only a subset of samples from the majority class. The proposed procedure is applied to a data set to predict estrogen receptor binding activity and to a data set to predict animal liver carcinogenicity using SAR (structure-activity relationship) models as base classifiers. The abatement method appears to perform well in balancing sensitivity and specificity.

Two-stage and Weibull models for carcinogenesis applied to the ED01 discontinued dosing data.

The two-stage clonal expansion model for a single, less-than-lifetime period of dosing is formulated and applied to the liver and bladder tumor data from the ED01 study. The model successfully predicts liver tumor incidence for time points beyond termination of dosing with 2-acetylaminofluorene, but it is unsuccessful for bladder tumor incidence. A discontinued dosing version of the Weibull model is proposed and is shown to predict successfully both liver and bladder tumor incidences for time points after termination of dosing.

Combining uncertainty factors in deriving human exposure levels of noncarcinogenic toxicants.

Acceptable levels of human exposure to noncarcinogenic toxicants in environmental and occupational settings generally are derived by reducing experimental no-observed-adverse-effect levels (NOAELs) or benchmark doses (BDs) by a product of uncertainty factors (Barnes and Dourson, Ref. 1). These factors are presumed to ensure safety by accounting for uncertainty in dose extrapolation, uncertainty in duration extrapolation, differential sensitivity between humans and animals, and differential sensitivity among humans. The common default value for each uncertainty factor is 10. This paper shows how estimates of means and standard deviations of the approximately log-normal distributions of individual uncertainty factors can be used to estimate percentiles of the distribution of the product of uncertainty factors. An appropriately selected upper percentile, for example, 95th or 99th, of the distribution of the product can be used as a combined uncertainty factor to replace the conventional product of default factors.

Dose-response trend tests for tumorigenesis, adjusted for body weight.

Several studies have demonstrated a relationship between rodent body weight and tumor incidence for some tissue/organ sites. It is not uncommon for a chemical tested for carcinogenicity to also affect body weight. In such cases, comparisons of tumor incidence may be biased by body-weight differences across dose groups. A simple procedure was investigated for reducing this bias. This procedure divides the animals into a few groups based on body weight. Body weight at 12 months was used, before the appearance of a tumor was likely to affect body weight. Statistics for dose-response trend tests are calculated within body weight strata and pooled to obtain an overall dose-response trend test. This procedure is analogous to that currently used, of stratifying animals, based on their age at the time of removal from a study. Age stratification is used to account for differences in animal age across dose groups, which can affect comparisons of tumor incidence. Several examples were investigated where the high-dose group had reduced body weights and associated reductions in tumor incidence. When the data were analyzed by body-weight strata, some positive dose-response trends for tumor incidence were demonstrated. In one case, the weight-adjusted analysis indicated that a negative dose-response trend in tumor incidence was a real effect, in addition to a body weight reduction. These examples indicate that it is important to consider the effects of body weight changes as low as 10%, and perhaps below, that were caused by chemicals in 2-year bioassays for carcinogenesis. The simple procedure of analyzing tumor incidence within body-weight strata can reduce the bias introduced by weight differences across dose groups.

Dose-response trend tests for tumorigenesis adjusted for differences in survival and body weight across doses.

A relationship between rodent body weight and tumor incidence for some tissue/organ sites has been demonstrated in many studies. It is not uncommon for a chemical tested for carcinogenicity to also affect body weight due to toxicity and/or food consumption. In such cases, comparisons of tumor incidence may be biased by body weight differences across dose groups. A simple procedure was investigated for reducing this bias. This procedure divides the animals into a few groups on the basis of body weight. Body weight at 12 months was used, before the appearance of a tumor was likely to affect body weight. Statistics for dose-response trend tests are calculated within body weight strata and pooled to obtain an overall dose-response trend test. This procedure is analogous to stratifying animals on the basis of age at the time of removal from a study to account for differences in ages of animals across dose groups that can affect comparisons of tumor incidence. In this paper, differences in survival times of animals were adjusted by the Poly-3 technique used by the National Toxicology Program. This technique does not require the assignment of cause of death. Several examples from rodent chronic bioassays were investigated, where the high dose group had reduced body weights and associated reductions in tumor incidence. When we analyzed the data by body weight strata, some positive dose-response trends for tumor incidence were demonstrated. In one case, the body weight adjusted analysis indicated that a negative dose-response trend in tumor incidence was a real effect in addition to a body weight reduction. These examples indicate that it is important to consider the effects of body weight changes as low as 10%, and perhaps less, as possibly being caused by chemicals in 2-year bioassays for carcinogenesis. The simple procedure of analyzing tumor incidence within body weight strata can reduce the bias introduced by body weight differences across dose groups.

Bioassays of shortened duration for drugs: statistical implications.

Declining survival rates in rodent carcinogenesis bioassays have raised a concern that continuing the practice of terminating such studies at 24 months could result in too few live animals at termination for adequate pathological evaluation. One option for ensuring sufficient numbers of animals at the terminal sacrifice is to shorten the duration of the bioassay, but this approach is accompanied by a reduction in statistical power for detecting carcinogenic potential. The present study was conducted to evaluate the loss of power associated with early termination. Data from drug studies in rats were used to formulate biologically based dose-response models of carcinogenesis using the 2-stage clonal expansion model as a context. These dose-response models, which were chosen to represent 6 variations of the initiation-promotion-completion cancer model, were employed to generate a large number of representative bioassay data sets using Monte Carlo simulation techniques. For a variety of tumor dose-response trends, tumor lethality, and competing risk-survival rates, the power of age-adjusted statistical tests to assess the significance of carcinogenic potential was evaluated at 18 and 21 months, and compared to the power at the normal 24-month stopping time. The results showed that stopping at 18 months would reduce power to an unacceptable level for all 6 submodels of the 2-stage clonal expansion model, with the pure-promoter and pure-completer models being most adversely affected. For the 21-month stopping time, the results showed that, unless pure promotion can be ruled out a priori as a potential carcinogenic mode of action, the loss of power is too great to warrant early stopping.

Statistical issues in the analysis of low-dose endocrine disruptor data.

UNLABELLED: The National Institute of Environmental Health Sciences (NIEHS) and the U.S. Environmental Protection Agency (U.S. EPA) recently cosponsored the Endocrine Disruptors Low-Dose Peer REVIEW: The purpose of this meeting was to examine data supporting the presence or absence of low-dose effects of endocrine disruptors in specific studies and then to evaluate the likelihood and significance of these and/or other potential low-dose effects for humans. All invited speakers agreed to provide their raw data in advance of the meeting to a Statistics Subpanel, which was asked to reevaluate the authors' experimental design, data analysis, and interpretation of experimental results. The purpose of this statistical reevaluation was to provide an independent assessment of the experimental design and data analysis used in each of the studies and to identify key statistical issues relevant to the evaluation and interpretation of the data. This paper presents a summary of the Statistics Subpanel's evaluation. Specific examples are presented to illustrate problems that arose in the experimental design and data analysis of certain studies. The statistical principles and issues that are discussed in this paper are not unique to endocrine disruptor studies and should provide important guidelines regarding appropriate experimental design and statistical analysis for other types of laboratory investigations.

Flow cytometric evaluation of cell-cycle progression in ethyl methanesulfonate and methyl methanesulfonate-exposed P3 cells: relationship to the induction of sister-chromatid exchanges and cellular toxicity.

In order to determine the relationships among the reduction in relative cloning efficiency (RCE), sister-chromatid exchange (SCE) formation, and interference with progression through the cell-cycle, human teratocarcinoma-derived (P3) cells were exposed to either ethyl methanesulfonate or to methyl methanesulfonate. The relationship between SCE and toxicity was quantified, the progression through the cell-cycle was evaluated with flow cytometric methods, and the effects of these chemicals on cell growth and average generation time (AGT) were determined. A strong correlation existed between RCE and SCE (r = -0.978, p less than .001) which was accompanied by an inhibition of growth as evidenced by a significant (p less than .0001) negative linear effect of concentration on the relative cell count from 24 to 72 hours after exposure and by a concentration-dependent increase (p less than .0001) in the AGT. Delays in the transit through S-phase were evident 4 hours after exposure to toxic concentrations of either carcinogen and by 8 to 12 hours post-exposure at the lower concentrations. Increases in the percentage of nuclei in G2 + M, indicative of G2 arrest, occurred from 12 to 24 hours after exposure. One interpretation of these results is that those effects of EMS and MMS exposure which result in S-phase delay and G2 arrest may be those elements common to the induction of SCE and cellular toxicity.

Detection of TEM-induced reciprocal translocations in F1 sons of CD-1 male mice: comparison of sequential fertility evaluation and cytogenetic analysis.

To determine the positive and negative classification error rates associated with the HTA in our laboratory, F1 sons of TEM-exposed CD-1 male mice were evaluated by the sequential fertility method with subsequent cytogenetic analysis. Males who sired three litters of size 10 or less when mated to primiparous females from either the B6C3F1 or the BCF1 strain were classified as partial steriles. When meiotic chromosome analyses revealed the presence of at least two cells containing multivalent figures, males were classified as translocation heterozygotes. When the fertility evaluation and the cytogenetic analysis were compared, normal fertility was observed on 5 of 83 (6.02%) translocation-bearing F1 males mated to B6C3F1 tester females and on 3 of 83 (3.61%) F1 males mated to BCF1 tester females. Thus, the false-negative error rates were 6.02% and 3.61% with these two tester strains. Multivalent figures were not observed in the meiotic chromosomes of 410 F1 males. Of these, 12 (2.93%) had reduced fertility when mated to the B6C3F1 tester strain as did 7 (1.71%) mated to the BCF1 strain. Thus, the false-positive error rates with these two tester strains were 2.93% for the B6C3F1 strain and 1.71% for the BCF1 strain. Our results indicate that non-zero error rates, both false-positive and false-negative, are associated with the sequential mating method HTA. In addition, the magnitude of these error rates was influenced not only by the tester female strain but also by the genotype of the F1 male.

Flow cytometric analysis of bromodeoxyuridine-induced inhibition of cell proliferation in the human teratocarcinoma-derived cell line, P3.

We have developed methods in our laboratory whereby the effects of toxicant exposure on cell proliferation can be evaluated flow cytometrically. We sought to relate the flow cytometric analyses to other biological response measurements. Thus, we exposed P3 cells to increasing concentrations of bromodeoxyuridine (BRdU) and measured sister-chromatid exchange (SCE) frequency, average generation time (AGT), and relative cloning ability. Each of these is well documented (see introduction) to respond to BRdU exposure in a concentration-dependent manner. In this study, SCE frequency remained constant between the concentrations of 2.5 and 10 microM of BRdU. However, a small, but significant, increase in SCE frequency was observed between the concentrations of 10 microM and 50 microM BRdU. A significant increase in AGT was noted in 50 microM BRdU-exposed cells. Relative cloning efficiency decreased in a concentration-dependent manner when cells were cultured for 24, 48, or 72 hours with BRdU. When cell proliferation was assessed by flow cytometric analysis in cells exposed to 0, 10, or 50 microM BRdU, a statistically significant delay in the cell-cycle was observed in BRdU-exposed cells. These results may be interpreted to mean that inhibition of cell proliferation is detected by this type of analysis at toxicant concentrations that induce other biological endpoints. The inclusion of flow cytometric analysis in a test battery to evaluate toxicant effects is warranted.

Effect of interleukin-2 on cell proliferation, sister-chromatid exchange induction, and nuclear stress protein phosphorylation in PHA-stimulated Fischer 344 rat spleen lymphocytes: modulation by 2-mercaptoethanol.

The effect of interleukin-2 (IL-2) on cell proliferation, sister-chromatid exchange (SCE) frequency, and the phosphorylation of nuclear stress proteins was evaluated in phytohemagglutinin (PHA)-stimulated spleen lymphocytes isolated from Fischer 344 rats. In addition, the ability of 2-mercaptoethanol (2-ME) to modulate the induction of these biological responses was characterized. Cell proliferation, as measured by the mitotic index, increased significantly (P less than .003) from a range of 3-4% in PHA-stimulated cultures to a range of 8-11% in PHA-stimulated cultures exposed to IL-2. The average generation time (AGT) did not respond to IL-2 in a concentration-dependent manner and decreased significantly (P less than .05) when 20 microM 2-ME was included with IL-2 in the culture medium. The number of SCE increased significantly (P less than .004) from control frequencies, which ranged from 13.1 to 15.6 SCE per cell, to frequencies of 18.5 to 21.5 SCE per cell as the concentration of IL-2 in the culture medium increased to 50 half-maximal units per ml. A reduction in SCE frequency was observed when cells were cultured with 20 microM 2-ME and IL-2 compared to IL-2 alone. Three nuclear proteins, with relative molecular masses of approximately 13,000-18,000, 20,000, and 80,000, were phosphorylated in IL-2-exposed G1-phase nuclei. Elicitation of these nuclear proteins in IL-2-exposed cells was not affected by exposure to 2-ME.

Upper bound risk estimates for mixtures of carcinogens.

The excess cancer risk that might result from exposure to a mixture of chemical carcinogens usually is estimated with data from experiments conducted on individual chemicals. An upper bound on the total excess risk is estimated commonly by summing individual upper bound risk estimates. The degree to which this approach might overstate the true risk associated with the mixture has not been evaluated previously. This paper reports the results of a Monte Carlo simulation study on the degree of reduction in conservation that might be achieved using alternative methods for calculating mixture upper bounds. An unexpected finding is that for chemicals that exhibit strongly linear dose-response relationships, the summing of multistage-model-based upper bounds on excess risk can be anti-conservative, that is, it can provide less than the nominal 100(1-alpha)% coverage.

The effects of ascorbic acid on the induction of urothelial lesions in mice by 2-acetylaminofluorene.

288 BALB/c male mice were allocated to: Group 1-control; Group 2-500 ppm of 2-acetylaminofluorene (2-AAF) in the food, with control water; Group 4-500 ppm of 2-AAF in the food and 250 mg/100 ml of vitamin C in the water; and Group 3-control food and 250 mg% of vitamin C in the water. Major histopathologic changes at the end of 28 days included both inflammation of the lamina propria and hyperplasia of the transitional epithelium of the urinary bladder in mice receiving 2-AAF. The most severe lesions were seen in the mice administered the combination of 2-AAF and vitamin C. It was postulated that either decreased water consumption and concentration of the urine and/or reduction in urinary pH may have contributed to the severity of the lesions.

Dominant lethal test response with IMS and TEM using different combinations of male and female stocks of mice.

We conducted dominant lethal studies with chemical mutagens IMS and TEM to investigate the influence which different treated male stocks might have upon individual female response. In both studies, treated males from a CD-1 random bred stock and each of two F1 hybrid stocks (C57BL/6N X C3H; C57BL/6J X BALB/c) were mated to untreated females from the two F1 hybrid stocks. We found that variation in response due to the treatment effect was modulated by the specific male stock/female stock combination involved. The male-dependent variation observed prevented any meaningful evaluation of relative female stock sensitivity but it was obvious that factors other than differences in repair capacity of female germ cells contributed to the response differences observed.

Correlation between specific DNA-methylation products and mutation induction at the HGPRT locus in Chinese hamster ovary cells.

Suspension cultures of Chinese hamster ovary (CHO) cells were exposed to methyl methanesulfonate (MMS) or methylnitrosourea (MNU) and assayed for mutation induction (6-thioguanine resistance) and for specific DNA adducts. DNA methylation at the 1-, 3- and 7-positions of adenine, the 3-, O6- and 7-positions of guanine, and phosphate was detected in cultures exposed to MMS, while MNU produced 3- and 7-methyladenine, 3-methylcytosine, 3-, O6- and 7-methylguanine, O4-methylthymidine and methylated phosphodiesters. When mutations induced by MMS and MNU were compared by linear correlation analysis with levels of each of these adducts, only O6-methylguanine displayed a strong correlation with mutations (r = 0.879, p less than 0.001). The relationship between O6-methylguanine and induced mutations in CHO cells is similar to that previously reported in CHO cells for O6-ethylguanine and mutations (Heflich et al., 1982) and indicates that alkylation-induced mutations at the HGPRT locus in CHO cells are primarily associated with O6-alkylguanine formation.