Search

In response to Timothy Darvill's article, ‘Mythical rings?’ (this issue), which argues for an alternative interpretation of Waun Mawn circle and its relationship with Stonehenge, Parker Pearson and colleagues report new evidence from the Welsh site and elaborate on aspects of their original argument. The discovery of a hearth at the centre of the circle, as well as further features around its circumference, reinforces the authors’ original interpretation. The authors explore the evidence for the construction sequence, which was abandoned before the completion of the monument. Contesting Darvill's argument that the Aubrey Holes at Stonehenge originally held posts, the authors reassert their interpretation of this circle of cut features as Bluestone settings.

Trace-level environmental data typically include values near or below detection and quantitation thresholds where health effects may result from low-concentration exposures to one chemical over time or to multiple chemicals. In a cook stove case study, bias in dibenzo[a,h]anthracene concentration means and standard deviations (SDs) was assessed following censoring at thresholds for selected analysis approaches: substituting threshold/2, maximum likelihood estimation, robust regression on order statistics, Kaplan-Meier, and omitting censored observations. Means and SDs for gas chromatography-mass spectrometry-determined concentrations were calculated after censoring at detection and calibration thresholds, 17% and 55% of the data, respectively. Threshold/2 substitution was the least biased. Measurement values were subsequently simulated from two log-normal distributions at two sample sizes. Means and SDs were calculated for 30%, 50%, and 80% censoring levels and compared to known distribution counterparts. Simulation results illustrated (1) threshold/2 substitution to be inferior to modern after-censoring statistical approaches and (2) all after-censoring approaches to be inferior to including all measurement data in analysis. Additionally, differences in stove-specific group means were tested for uncensored samples and after censoring. Group differences of means tests varied depending on censoring and distributional decisions. Investigators should guard against censoring-related bias from (explicit or implicit) distributional and analysis approach decisions.

PBPK models are well-established frameworks used to describe absorption, distribution, metabolism, and excretion of xenobiotics. To quantify metabolism, a PBPK model for a volatile compound can be calibrated with closed chamber (i.e., vapor uptake) inhalation data. Here, we introduce global optimization as a novel component of the predictive process and use it to illustrate a procedure for metabolic parameter estimation. Male F344 rats were exposed in vapor uptake chambers to initial concentrations of 100, 500, 1000, and 3000 ppm chloroform. Chamber time-course data from these experiments, in combination with optimization using a chemical specific PBPK model, were used to estimate Michaelis-Menten metabolic constants. Matlab® simulation software was used to integrate the mass balance equations and to perform the global optimizations using MEIGO (MEtaheuristics for systems biology and bIoinformatics Global Optimization – Version 64bit, R2016A), a toolbox written for Matlab®). The cost function used the chamber time-course data and least squares to minimize the difference between data and simulation values. The final values estimated for V(max) (maximum metabolic rate) and K(m) (affinity constant) were 1.2 mg/h and a range between 0.0005 and 0.6 mg/L, respectively. Also, cost function plots were used to analyze the dose-dependent capacity to estimate V(max) and K(m) within the experimental range used. Sensitivity analysis was used to assess identifiability for both parameters and show these kinetic data may not be sufficient to identify K(m). The importance of an accurate estimation of metabolism using vapor uptake data is discussed. In summary, this work should help toxicologists interested in optimization techniques understand the overall process employed when calibrating metabolic parameters in a PBPK model with inhalation data.

Trace-level environmental data typically include values near or below detection and quantitation thresholds where health effects may result from low-concentration exposures to one chemical over time or to multiple chemicals. In a cook stove case study, bias in dibenzo[a,h]anthracene concentration means and standard deviations (SDs) was assessed following censoring at thresholds for selected analysis approaches: substituting threshold/2, maximum likelihood estimation, robust regression on order statistics, Kaplan-Meier, and omitting censored observations. Means and SDs for gas chromatography-mass spectrometry-determined concentrations were calculated after censoring at detection and calibration thresholds, 17% and 55% of the data, respectively. Threshold/2 substitution was the least biased. Measurement values were subsequently simulated from two log-normal distributions at two sample sizes. Means and SDs were calculated for 30%, 50%, and 80% censoring levels and compared to known distribution counterparts. Simulation results illustrated (1) threshold/2 substitution to be inferior to modern after-censoring statistical approaches and (2) all after-censoring approaches to be inferior to including all measurement data in analysis. Additionally, differences in stove-specific group means were tested for uncensored samples and after censoring. Group differences of means tests varied depending on censoring and distributional decisions. Investigators should guard against censoring-related bias from (explicit or implicit) distributional and analysis approach decisions.

In a recent U.S. Geological Survey/U.S. Environmental Protection Agency study assessing more than 700 organic compounds in 38 streams, in vitro assays indicated generally low estrogen, androgen, and glucocorticoid receptor activities, with 13 surface waters with 17β-estradiol-equivalent (E2Eq) activities greater than a 1-ng/L estimated effects-based trigger value for estrogenic effects in male fish. Among the 36 samples assayed for mutagenicity in the Salmonella bioassay (reported here), 25% had low mutagenic activity and 75% were not mutagenic. Endocrine and mutagenic activities of the water samples were well correlated with each other and with the total number and cumulative concentrations of detected chemical contaminants. To test the predictive utility of knowledge-base-leveraging approaches, site-specific predicted chemical-gene (pCGA) and predicted analogous pathway-linked (pPLA) association networks identified in the Comparative Toxicogenomics Database were compared with observed endocrine/mutagenic bioactivities. We evaluated pCGA/pPLA patterns among sites by cluster analysis and principal component analysis and grouped the pPLA into broad mode-of-action classes. Measured E2eq and mutagenic activities correlated well with predicted pathways. The pPLA analysis also revealed correlations with signaling, metabolic, and regulatory groups, suggesting that other effects pathways may be associated with chemical contaminants in these waters and indicating the need for broader bioassay coverage to assess potential adverse impacts.

Geologists and archaeologists have long known that the bluestones of Stonehenge came from the Preseli Hills of west Wales, 230km away, but only recently have some of their exact geological sources been identified. Two of these quarries—Carn Goedog and Craig Rhos-y-felin—have now been excavated to reveal evidence of megalith quarrying around 3000 BC—the same period as the first stage of the construction of Stonehenge. The authors present evidence for the extraction of the stone pillars and consider how they were transported, including the possibility that they were erected in a temporary monument close to the quarries, before completing their journey to Stonehenge.

Contaminants of emerging concern (CECs), including per- and polyfluoroalkyl substances (PFAS), are of interest to regulators, water treatment utilities, the general public and scientists. This study measured 17 PFAS in source and treated water from 25 drinking water treatment plants (DWTPs) as part of a broader study of CECs in drinking water across the United States. PFAS were quantitatively detected in all 50 samples, with summed concentrations of the 17 PFAS ranging from, GRAPHICAL ABSTRACT

The discovery of a dismantled stone circle—close to Stonehenge's bluestone quarries in west Wales—raises the possibility that a 900-year-old legend about Stonehenge being built from an earlier stone circle contains a grain of truth. Radiocarbon and OSL dating of Waun Mawn indicate construction c. 3000 BC, shortly before the initial construction of Stonehenge. The identical diameters of Waun Mawn and the enclosing ditch of Stonehenge, and their orientations on the midsummer solstice sunrise, suggest that at least part of the Waun Mawn circle was brought from west Wales to Salisbury Plain. This interpretation complements recent isotope work that supports a hypothesis of migration of both people and animals from Wales to Stonehenge.

BACKGROUND: In developmental and reproductive toxicity studies, analysis of litter-based binary endpoints (e.g., incidence of malformed fetuses) is complex in that littermates often are not entirely independent of one another. It is well established that the litter, not the individual fetus, is the proper independent experimental unit in statistical analysis. Accordingly, analysis is often based on the proportion affected per litter and the litter proportions are analyzed as continuous data. Because these proportional data generally do not meet assumptions of symmetry or normality, data are typically analyzed by nonparametric methods, arcsine square root transformation, or logit transformation. METHODS: We conducted power calculations to compare different approaches (nonparametric, arcsine square root-transformed, logit-transformed, untransformed) for analyzing litter-based proportional data. A reproductive toxicity study with a control and one treated group provided data for two endpoints: prenatal loss, and fertility by in utero insemination (IUI). Type 1 error and power were estimated by 10,000 simulations based on two-sample one-tailed t-tests with varying numbers of litters per group. To further compare the different approaches, we conducted additional analyses with shifted mean proportions to produce illustrative scenarios. RESULTS: Analyses based on logit-transformed proportions had greater power than those based on untransformed or arcsine square root-transformed proportions, or nonparametric procedures. CONCLUSION: The logit transformation is preferred to the other approaches considered when making inferences concerning litter-based proportional endpoints, particularly with skewed distributions. The improved performance of the logit transformation becomes increasingly pronounced as the response proportions are increasingly close to the boundaries of the parameter space.

Short-term biomarkers of toxicity have an increasingly important role in the screening and prioritization of new chemicals. In this study, we examined early indicators of liver toxicity for three reference organophosphate (OP) chemicals, which are among the most widely used insecticides in the world. The OP methidathion was previously shown to increase the incidence of liver toxicity, including hepatocellular tumors, in male mice. To provide insights into the adverse outcome pathway (AOP) that underlies these tumors, effects of methidathion in the male mouse liver were examined after 7 and 28 day exposures and compared to those of two other OPs that either do not increase (fenthion) or possibly suppress liver cancer (parathion) in mice. None of the chemicals caused increases in liver weight/body weight or histopathological changes in the liver. Parathion decreased liver cell proliferation after 7 and 28 days while the other chemicals had no effects. There was no evidence for hepatotoxicity in any of the treatment groups. Full-genome microarray analysis of the livers from the 7 and 28 day treatments demonstrated that methidathion and fenthion regulated a large number of overlapping genes, while parathion regulated a unique set of genes. Examination of cytochrome P450 enzyme activities and use of predictive gene expression biomarkers found no consistent evidence for activation of AhR, CAR, PXR, or PPARα. Parathion suppressed the male-specific gene expression pattern through STAT5b, similar to genetic and dietary conditions that decrease liver tumor incidence in mice. Overall, these findings indicate that methidathion causes liver cancer by a mechanism that does not involve common mechanisms of liver cancer induction.

Current strategies for predicting carcinogenic mode of action for nongenotoxic chemicals are based on identification of early key events in toxicity pathways. The goal of this study was to evaluate short-term key event indicators resulting from exposure to androstenedione (A4), an androgen receptor agonist and known liver carcinogen in mice. Liver cancer is more prevalent in men compared with women, but androgen-related pathways underlying this sex difference have not been clearly identified. Short-term hepatic effects of A4 were compared with reference agonists of the estrogen receptor (ethinyl estradiol, EE) and glucocorticoid receptor (prednisone, PRED). Male B6C3F1 mice were exposed for 7 or 28 days to A4, EE, or PRED. EE increased and PRED suppressed hepatocyte proliferation, while A4 had no detectable effects. In a microarray analysis, EE and PRED altered >3000 and >670 genes, respectively, in a dose-dependent manner, whereas A4 did not significantly alter any genes. Gene expression was subsequently examined in archival liver samples from male and female B6C3F1 mice exposed to A4 for 90 days. A4 altered more genes in females than males and did not alter expression of genes linked to activation of the mitogenic xenobiotic receptors AhR, CAR, and PPARα in either sex. A gene expression biomarker was used to show that in female mice, the high dose of A4 activated the growth hormone-regulated transcription factor STAT5b, which controls sexually dimorphic gene expression in the liver. These findings suggest that A4 induces subtle age-related effects on STAT5b signaling that may contribute to the higher risk of liver cancer in males compared with females.

A method based on regression modeling was developed to discern the contribution of component chemicals to the toxicity of highly complex, environmentally realistic mixtures of disinfection byproducts (DBPs). Chemical disinfection of drinking water forms DBP mixtures. Because of concerns about possible reproductive and developmental toxicity, a whole mixture (WM) of DBPs produced by chlorination of a water concentrate was administered as drinking water to Sprague–Dawley (S–D) rats in a multigenerational study. Age of puberty acquisition, i.e., preputial separation (PPS) and vaginal opening (VO), was examined in male and female offspring, respectively. When compared to controls, a slight, but statistically significant delay in puberty acquisition was observed in females but not in males. WM-induced differences in the age at puberty acquisition were compared to those reported in S–D rats administered either a defined mixture (DM) of nine regulated DBPs or individual DBPs. Regression models were developed using individual animal data on age at PPS or VO from the DM study. Puberty acquisition data reported in the WM and individual DBP studies were then compared with the DM models. The delay in puberty acquisition observed in the WM-treated female rats could not be distinguished from delays predicted by the DM regression model, suggesting that the nine regulated DBPs in the DM might account for much of the delay observed in the WM. This method is applicable to mixtures of other types of chemicals and other endpoints.

When chemical or microbial contaminants are assessed for potential effect or possible regulation in ambient and drinking waters, a critical first step is determining if the contaminants occur and if they are at concentrations that may cause human or ecological health concerns. To this end, source and treated drinking water samples from29 drinking water treatment plants (DWTPs) were analyzed as part of a two-phase study to determine whether chemical and microbial constituents, many of which are considered contaminants of emerging concern, were detectable in the waters. Of the 84 chemicals monitored in the 9 Phase I DWTPs, 27 were detected at least once in the source water, and 21 were detected at least once in treated drinking water. In Phase II, which was a broader and more comprehensive assessment, 247 chemical and microbial analytes were measured in 25 DWTPs, with 148 detected at least once in the source water, and 121 detected at least once in the treated drinking water. The frequency of detection was often related to the analyte’s contaminant class, as pharmaceuticals and anthropogenic waste indicators tended to be infrequently detected and more easily removed during treatment, while per and polyfluoroalkyl substances and inorganic constituents were both more frequently detected and, overall, more resistant to treatment. The data collected as part of this project will be used to help inform evaluation of unregulated contaminants in surface water, groundwater, and drinking water., GRAPHICAL ABSTRACT

There is increasing public concern regarding health risks from multiple chemical exposures and combined exposure to chemical and nonchemical stressors (cumulative risk). These exposures may involve more than one exposure pathway and typically vary over time with respect to the total exposure level and the relative proportions of component chemicals. Research on toxicology and risk assessment of multichemical exposures is complicated further because toxicologic interactions vary with exposure sequence, duration and route, as well as with total dose, relative chemical proportions, and endpoint.

The U.S. Environmental Protection Agency (EPA) is faced with the challenge of efficiently and credibly evaluating chemical safety often with limited or no available toxicity data. The expanding number of chemicals found in commerce and the environment, coupled with time and resource requirements for traditional toxicity testing and exposure characterization, continue to underscore the need for new approaches. In 2005, EPA charted a new course to address this challenge by embracing computational toxicology (CompTox) and investing in the technologies and capabilities to push the field forward. The return on this investment has been demonstrated through results and applications across a range of human and environmental health problems, as well as initial application to regulatory decision-making within programs such as the EPA’s Endocrine Disruptor Screening Program. The CompTox initiative at EPA is more than a decade old. This manuscript presents a blueprint to guide the strategic and operational direction over the next 5 years. The primary goal is to obtain broader acceptance of the CompTox approaches for application to higher tier regulatory decisions, such as chemical assessments. To achieve this goal, the blueprint expands and refines the use of high-throughput and computational modeling approaches to transform the components in chemical risk assessment, while systematically addressing key challenges that have hindered progress. In addition, the blueprint outlines additional investments in cross-cutting efforts to characterize uncertainty and variability, develop software and information technology tools, provide outreach and training, and establish scientific confidence for application to different public health and environmental regulatory decisions.

Iodinated contrast media (ICM) are nonmutagenic agents administered for X-ray imaging of soft tissues. ICM can reach μg/L levels in surface waters because they are administered in high doses, excreted largely unmetabolized, and poorly removed by wastewater treatment. Iodinated disinfection byproducts (I-DBPs) are highly genotoxic and have been reported in disinfected waters containing ICM. We assessed the mutagenicity in Salmonella of extracts of chlorinated source water containing one of four ICM (iopamidol, iopromide, iohexol, and diatrizoate). We quantified 21 regulated and nonregulated DBPs and 11 target I-DBPs and conducted a nontarget, comprehensive broad-screen identification of I-DBPs. We detected one new iodomethane (trichloroiodomethane), three new iodoacids (dichloroiodoacetic acid, chlorodiiodoacetic acid, bromochloroiodoacetic acid), and two new nitrogenous I-DBPs (iodoacetonitrile and chloroiodoacetonitrile). Their formation depended on the presence of iopamidol as the iodine source; identities were confirmed with authentic standards when available. This is the first identification in simulated drinking water of chloroiodoacetonitrile and iodoacetonitrile, the latter of which is highly cytotoxic and genotoxic in mammalian cells. Iopamidol (5 μM) altered the concentrations and relative distribution of several DBP classes, increasing total haloacetonitriles by >10-fold. Chlorination of ICM-containing source water increased I-DBP concentrations but not mutagenicity, indicating that such I-DBPs were either not mutagenic or at concentrations too low to affect mutagenicity.

Prioritizing the risk posed by thousands of chemicals potentially present in the environment requires exposure, toxicity, and toxicokinetic (TK) data, which are often unavailable. Relatively high throughput, in vitro TK (HTTK) assays and in vitro-to-in vivo extrapolation (IVIVE) methods have been developed to predict TK, but most of the in vivo TK data available to benchmark these methods are from pharmaceuticals. Here we report on new, in vivo rat TK experiments for 26 non-pharmaceutical chemicals with environmental relevance. Both intravenous and oral dosing were used to calculate bioavailability. These chemicals, and an additional 19 chemicals (including some pharmaceuticals) from previously published in vivo rat studies, were systematically analyzed to estimate in vivo TK parameters (e.g., volume of distribution [Vd], elimination rate). For each of the chemicals, rat-specific HTTK data were available and key TK predictions were examined: oral bioavailability, clearance, Vd, and uncertainty. For the non-pharmaceutical chemicals, predictions for bioavailability were not effective. While no pharmaceutical was absorbed at less than 10%, the fraction bioavailable for non-pharmaceutical chemicals was as low as 0.3%. Total clearance was generally more underestimated for nonpharmaceuticals and Vd methods calibrated to pharmaceuticals may not be appropriate for other chemicals. However, the steady-state, peak, and time-integrated plasma concentrations of nonpharmaceuticals were predicted with reasonable accuracy. The plasma concentration predictions improved when experimental measurements of bioavailability were incorporated. In summary, HTTK and IVIVE methods are adequately robust to be applied to high throughput in vitro toxicity screening data of environmentally relevant chemicals for prioritizing based on human health risks. © The Author(s) 2018. Published by Oxford University Press on behalf of the Society of Toxicology. All rights reserved.

Current strategies for predicting adverse health outcomes of environmental chemicals are centered on early key events in toxicity pathways. However, quantitative relationships between early molecular changes in a given pathway and later health effects are often poorly defined. The goal of this study was to evaluate short-term key event indicators using qualitative and quantitative methods in an established pathway of mouse liver tumorigenesis mediated by peroxisome proliferator-activated receptor alpha (PPARα). Male B6C3F1 mice were exposed for 7 days to di (2-ethylhexyl) phthalate (DEHP), di-n-octyl phthalate (DNOP), and n-butyl benzyl phthalate (BBP), which vary in PPARα activity and liver tumorigenicity. Each phthalate increased expression of select PPARα target genes at 7 days, while only DEHP significantly increased liver cell proliferation labeling index (LI). Transcriptional benchmark dose (BMDT) estimates for dose-related genomic markers stratified phthalates according to hypothetical tumorigenic potencies, unlike BMDs for non-genomic endpoints (relative liver weights or proliferation). The 7-day BMDT values for Acot1 as a surrogate measure for PPARα activation were 29, 370, and 676 mg/kg/day for DEHP, DNOP, and BBP, respectively, distinguishing DEHP (liver tumor BMD of 35 mg/kg/day) from non-tumorigenic DNOP and BBP. Effect thresholds were generated using linear regression of DEHP effects at 7 days and 2-year tumor incidence values to anchor early response molecular indicators and a later phenotypic outcome. Thresholds varied widely by marker, from 2-fold (Pdk4 and proliferation LI) to 30-fold (Acot1) induction to reach hypothetical tumorigenic expression levels. These findings highlight key issues in defining thresholds for biological adversity based on molecular changes.

The discovery of Neolithic houses at Durrington Walls that are contemporary with the main construction phase of Stonehenge raised questions as to their interrelationship. Was Durrington Walls the residence of the builders of Stonehenge? Were the activities there more significant than simply domestic subsistence? Using lipid residue analysis, this paper identifies the preferential use of certain pottery types for the preparation of particular food groups and differential consumption of dairy and meat products between monumental and domestic areas of the site. Supported by the analysis of faunal remains, the results suggest seasonal feasting and perhaps organised culinary unification of a diverse community.

The introduction of drinking water disinfection greatly reduced waterborne diseases. However, the reaction between disinfectants and natural organic matter in the source water leads to an unintended consequence, the formation of drinking water disinfection byproducts (DBPs). The haloace-taldehydes (HALs) are the third largest group by weight of identified DBPs in drinking water. The primary objective of this study was to analyze the occurrence and comparative toxicity of the emerging HAL DBPs. A new HAL DBP, iodoacetaldehyde (IAL) was identified. This study provided the first systematic, quantitative comparison of HAL toxicity in Chinese hamster ovary cells. The rank order of HAL cytotoxicity is tribromoacetaldehyde (TBAL) ≈ chloroacetaldehyde (CAL) > dibromoacetaldehyde (DBAL) ≈ bromochloroacetaldehyde (BCAL) ≈ dibromochloroacetaldehyde (DBCAL) > IAL > bromoacetaldehyde (BAL) ≈ bromodichloroacetaldehyde (BDCAL) > dichloroacetaldehyde (DCAL) > trichloroacetaldehyde (TCAL). The HALs were highly cytotoxic compared to other DBP chemical classes. The rank order of HAL genotoxicity is DBAL > CAL ≈ DBCAL > TBAL ≈ BAL > BDCAL > BCAL ≈ DCAL > IAL. TCAL was not genotoxic. Because of their toxicity and abundance, further research is needed to investigate their mode of action to protect the public health and the environment.

Background Trihalomethanes (THMs) and haloacetic acids (HAAs) are regulated disinfection by-products (DBPs); their joint reproductive toxicity in drinking water is unknown. Objective We aimed to evaluate a drinking water mixture of the four regulated THMs and five regulated HAAs in a multigenerational reproductive toxicity bioassay. Methods Sprague-Dawley rats were exposed (parental, F1, and F2 generations) from gestation day 0 of the parental generation to postnatal day (PND) 6 of the F2 generation to a realistically proportioned mixture of THMs and HAAs at 0, 500×, 1,000×, or 2,000× of the U.S. Environmental Protection Agency’s maximum contaminant levels (MCLs). Results Maternal water consumption was reduced at ≥ 1,000×; body weights were reduced at 2,000×. Prenatal and postnatal survival were unaffected. F1 pup weights were unaffected at birth but reduced at 2,000× on PND6 and at ≥ 1,000× on PND21. Postweaning F1 body weights were reduced at 2,000×, and water consumption was reduced at ≥ 500×. Males at 2,000× had a small but significantly increased incidence of retained nipples and compromised sperm motility. Onset of puberty was delayed at 1,000× and 2,000×. F1 estrous cycles and fertility were unaffected, and F2 litters showed no effects on pup weight or survival. Histologically, P0 (parental) dams had nephropathy and adrenal cortical pathology at 2,000×. Conclusions A mixture of regulated DBPs at up to 2,000× the MCLs had no adverse effects on fertility, pregnancy maintenance, prenatal survival, postnatal survival, or birth weights. Delayed puberty at ≥ 1,000× may have been secondary to reduced water consumption. Male nipple retention and compromised sperm motility at 2,000× may have been secondary to reduced body weights. Citation Narotsky MG, Klinefelter GR, Goldman JM, DeAngelo AB, Best DS, McDonald A, Strader LF, Murr AS, Suarez JD, George MH, Hunter ES III, Simmons JE. 2015. Reproductive toxicity of a mixture of regulated drinking-water disinfection by-products in a multigenerational rat bioassay. Environ Health Perspect 123:564–570; http://dx.doi.org/10.1289/ehp.1408579

There is a general consensus, based on a number of surveys and analytic efforts, that mixture study designs have historically been lacking. Although there has been considerable progress over the past decades, further improvement is necessary both in the development and application of experimental designs to yield data suitable for quantitative analytic methods and in the implementation of appropriate statistical analyses. This chapter reviews the state of the science with regard to the experimental and statistical quality of mixture studies. The importance of properly powering mixture experiments is emphasized; in particular when the focus is on the low-dose/low-effect region. Issues with powering defined mixture and complex mixture experiments are explored. Some designs that have proven useful in mixture experimentation are reviewed, including full and fractional factorial designs and statistical mixture designs such as the isobologram and the fixed ratio ray. General considerations are provided that will aid in the development of both experimental design and analysis strategies that address the question(s) being asked.

Researchers in numerous fields (e.g., pharmacology, entomology, toxicology, and epidemiology) have attempted to model the joint action of chemicals using simple formulas based only on knowledge of individual chemical toxicity or pharmacological effect (i.e., dose-response relationships). Collectively, these formulas are referred to as “additivity models,” and they are based on concepts of additivity that include dose addition, independent action, integrated addition, and effect summation. In toxicology, additivity-based predictions are often compared to observed mixture data to assess the presence and magnitude of interactions (greater-than-additive or less-than-additive) among chemicals. These models can also be used to estimate the toxicity of a defined mixture for comparison to the observed toxicity of a related, but more complex, mixture. Alternatively, additivity models have been used to explore mechanisms of joint action. In general, the steps for investigating joint toxicity using additivity models include (1) deciding on which additivity model(s) to apply (e.g., dose addition, independent action, or both), (2) collecting dose-response data on individual chemicals and the mixture, (3) incorporating individual chemical data in an additivity model to generate predictions, and (4) comparing predicted to observed mixture responses. Many of the additivity models have a long and sometimes controversial history. This chapter provides background on several of the common additivity models, illustrates their application with examples, and discusses their advantages and limitations.

Due to the ever-increasing number of chemicals coming to market, and the cost of performing traditional in vivo studies, there has been a shift toward the use of less costly alternative techniques. The adverse outcome pathway (AOP) concept has emerged as a scaffold for organizing mechanistic information from these methods. Two main elements – key events (KEs) and key event relationships (KERs) – are utilized to describe the underlying mechanism outlined by the AOP. Each KE depicts the measureable changes in the state of the biological system at each level of organization that are essential for the progression along the pathway. The KERs, meanwhile, contain the biological information that connects each of the KEs. This chapter covers some of the potential applications for AOPs when performing risk assessment of chemical mixtures. The structure of the AOP provides much more precision when considering mechanistic data in a mixtures assessment. The use of this concept provides a means to allow more specificity when deciding whether to use dose addition, independent action or integrated addition risk assessment methodologies. Furthermore, AOPs enable novel approaches for determining chemical groups and how they may be utilized within mixtures risk assessment.

As the use of alternative drinking water treatment increases, it is important to understand potential public health implications associated with these processes. The objective of this study was to evaluate the formation of disinfection byproducts (DBPs) and cytotoxicity of natural organic matter (NOM) concentrates treated with chlorine, chloramine, and medium pressure ultraviolet (UV) irradiation followed by chlorine or chloramine, with and without nitrate or iodide spiking. The use of concentrated NOM conserved volatile DBPs and allowed for direct analysis of the treated water. Treatment with UV prior to chlorine in ambient (unspiked) samples did not affect cytotoxicity as measured using an in vitro normal human colon cell (NCM460) assay, compared to chlorination alone when toxicity is expressed on the basis of dissolved organic carbon (DOC). Nitrate-spiked UV+chlorine treatment produced greater cytotoxicity than nitrate-spiked chlorine alone or ambient UV+chlorine samples, on both a DOC and total organic halogen basis. Samples treated with UV+chloramine were more cytotoxic than those treated with only chloramine using either dose metric. This study demonstrated the combination of cytotoxicity and DBP measurements for process evaluation in drinking water treatment. The results highlight the importance of dose metric when considering the relative toxicity of complex DBP mixtures formed under different disinfection scenarios.

The extensification of agricultural systems into marginal lands is a common response to environmental, economic, and political pressures for more cultivable land. Yet the course that extensification takes in particular instances is unpredictable given the choices available to producers. This article investigates an instance of extensification during the late second millennium BCE on the semi-arid Eastern Karak Plateau in west-central Jordan. Architectural, faunal, and archaeobotanical evidence is presented from Khirbat al-Mudayna al-'Aliya, one of several communities that participated in an extensified settlement system on the edge of the Wadi al-Mujib and its tributaries. Producers practiced agriculture and pastoralism in a low-intensity subsistence economy that supported a nucleated settlement of households. Faunal analysis determined goats were kept, and wild animals supplemented diets. Archaeobotanical analysis of charred plant remains from storage bins in a building destroyed by fire indicated that barley was stored in a semi-processed state and that harvesting by uprooting was practiced, thus resulting in the maximization of the straw harvest. The riparian zone beneath the settlement was a key venue for subsistence activities. This Early Iron Age example contrasts with later episodes of extensification whose settlement systems were more dispersed and agro-pastoralist regime more integrated.

Previously, we reported that acute treatment with propoxur or carbaryl decreased the duration of the photic after discharge (PhAD) of flash evoked potentials (FEPs). In the current studies, we compared the effects of acute or repeated exposure to a mixture of carbaryl and propoxur (1:1.45 ratio; propoxur:carbaryl) on the duration of the PhAD and brain ChE activity in Long Evans rats. Animals were exposed (po) either to a single dose (0, 3, 10, 45 or 75 mg/kg), or 14 daily dosages (0, 3, 10, 30, 45 mg/kg), of the mixture. Acute and repeated treatment with 3mg/kg (or greater) of the mixture produced dose-related inhibition of brain ChE activity. Compared to controls, the PhAD duration decreased after acute administration of 75 mg/kg or repeated treatment with 30 mg/kg of the mixture. The linear relationship between the percent of control brain ChE activity and the PhAD duration was similar for both exposure paradigms. Dose-response models for the acute and repeated exposure data did not differ for brain ChE activity or the duration of the PhAD. Repeated treatment with the mixture resulted in slightly less (13-22%) erythrocyte ChE inhibition than acute exposure. Both acute and repeated treatment resulted in dose-additive results for the PhAD duration and less than dose-additive responses (6-16%) for brain ChE activity for the middle range of dosages. Acute treatment resulted in greater than dose-additive erythrocyte ChE inhibition (15-18%) at the highest dosages. In contrast, repeated treatment resulted in less than dose-additive erythrocyte ChE inhibition (16-22%) at the middle dosages. Brain and plasma levels of propoxur and carbaryl did not differ between the acute and repeated dosing paradigms. In summary, a physiological measure of central nervous system function and brain ChE activity had similar responses after acute or repeated treatment with the carbamate mixture, and brain ChE showed only small deviations from dose-additivity. Erythrocyte ChE activity had larger differences between the acute and repeated treatment paradigms, and showed slightly greater deviations from dose-additivity. Because these treatments utilized larger dosages than anticipated environmental exposures, concern for non-additive effects in humans is minimized. The small magnitude of the deviations from dose-additivity also suggest that in the absence of repeated exposure data, results from an acute study of readily reversible carbamate toxicity can be used to estimate the response to repeated daily exposures.

Statistical design and environmental relevance are important aspects of studies of chemical mixtures, such as pesticides. We used a dose-additivity model to test experimentally the default assumptions of dose additivity for two mixtures of seven N-methylcarbamates (carbaryl, carbofuran, formetanate, methomyl, methiocarb, oxamyl, and propoxur). The best-fitting models were selected for the single-chemical dose-response data and used to develop a combined prediction model, which was then compared with the experimental mixture data. We evaluated behavioral (motor activity) and cholinesterase (ChE)-inhibitory (brain, red blood cells) outcomes at the time of peak acute effects following oral gavage in adult and preweanling (17 days old) Long-Evans male rats. The mixtures varied only in their mixing ratios. In the relative potency mixture, proportions of each carbamate were set at equitoxic component doses. A California environmental mixture was based on the 2005 sales of each carbamate in California. In adult rats, the relative potency mixture showed dose additivity for red blood cell ChE and motor activity, and brain ChE inhibition showed a modest greater-than additive (synergistic) response, but only at a middle dose. In rat pups, the relative potency mixture was either dose-additive (brain ChE inhibition, motor activity) or slightly less-than additive (red blood cell ChE inhibition). On the other hand, at both ages, the environmental mixture showed greater-than additive responses on all three endpoints, with significant deviations from predicted at most to all doses tested. Thus, we observed different interactive properties for different mixing ratios of these chemicals. These approaches for studying pesticide mixtures can improve evaluations of potential toxicity under varying experimental conditions that may mimic human exposures.

A developmental toxicity bioassay was used in three experiments to evaluate water concentrates for suitability in multigenerational studies. First, chlorinated water was concentrated 135-fold by reverse osmosis; select lost disinfection by-products were spiked back. Concentrate was provided as drinking water to Sprague-Dawley and F344 rats from gestation day 6 to postnatal day 6. Maternal serum levels of luteinizing hormone on gestation day 10 were unaffected by treatment for both strains. Treated dams had increased water consumption, and increased incidences of polyuria, diarrhea, and (in Sprague-Dawley rats) red perinasal staining. Pup weights were reduced. An increased incidence of eye defects was seen in F344 litters. Chemical analysis of the concentrate revealed high sodium (6.6 g/l) and sulfate (10.4 g/l) levels. To confirm that these chemicals caused polyuria and osmotic diarrhea, respectively, Na₂SO₄ (5-20 g/l) or NaCl (16.5 g/l) was provided to rats in drinking water. Water consumption was increased at 5- and 10-g Na₂SO₄/l and with NaCl. Pup weights were reduced at 20-g Na₂SO₄/l. Dose-related incidences and severity of polyuria and diarrhea occurred in Na₂SO₄-treated rats; perinasal staining was seen at 20 g/l. NaCl caused polyuria and perinasal staining, but not diarrhea. Subsequently, water was concentrated ∼120-fold and sulfate levels were reduced by barium hydroxide before chlorination, yielding lower sodium (≤1.5 g/l) and sulfate (≤2.1 g/l) levels. Treatment resulted in increased water consumption, but pup weight and survival were unaffected. There were no treatment-related clinical findings, indicating that mixtures produced by the second method are suitable for multigenerational testing.

In rodent bioassays where chemicals are administered in the drinking water, water consumption data for individual animals are needed to estimate chemical exposures accurately. If multiple animals share a common water source, as occurs in some studies, only the total amount of drinking water consumed by all animals utilizing the common source is directly measurable, and water consumption rates for individual animals are not available. In the Four Lab Study of the US Environmental Protection Agency, which included a multigenerational rodent bioassay, a complex mixture of drinking water disinfection by-products was delivered to multiple Sprague-Dawley rats from a common drinking water container. To estimate disinfection by-product mixture exposure for each animal, authors developed four log-linear regression models to allocate water consumption among rats sharing a common water container. The four models represented three animal lifestages: Gestation, Lactation, and Postweaning, with separate Postweaning models for male and female. Authors used data from six Sprague-Dawley rat bioassays to develop these models from available individual cage data for the Postweaning models, and available individual animal data for the Gestation and Lactation models. The r(2) values for the model fits were good, ranging from 0.67 to 0.92. The Gestation and Lactation models were generally quite accurate in predicting average daily water consumption whereas the Postweaning models were less robust. These models can be generalized for use in other reproductive and developmental bioassays where common water sources are used and data on the explanatory variables are available.

In complex mixture toxicology, there is growing emphasis on testing environmentally representative doses that improve the relevance of results for health risk assessment, but are typically much lower than those used in traditional toxicology studies. Traditional experimental designs with typical sample sizes may have insufficient statistical power to detect effects caused by environmentally relevant doses. Proper study design, with adequate statistical power, is critical to ensuring that experimental results are useful for environmental health risk assessment. Studies with environmentally realistic complex mixtures have practical constraints on sample concentration factor and sample volume as well as the number of animals that can be accommodated. This article describes methodology for calculation of statistical power for non-independent observations for a multigenerational rodent reproductive/developmental bioassay. The use of the methodology is illustrated using the U.S. EPA’s Four Lab study in which rodents were exposed to chlorinated water concentrates containing complex mixtures of drinking water disinfection by-products. Possible experimental designs included two single-block designs and a two-block design. Considering the possible study designs and constraints, a design of two blocks of 100 females with a 40:60 ratio of control:treated animals and a significance level of 0.05 yielded maximum prospective power (~90%) to detect pup weight decreases, while providing the most power to detect increased prenatal loss.

The U.S. Environmental Protection Agency's "Four Lab Study" involved participation of researchers from four national Laboratories and Centers of the Office of Research and Development along with collaborators from the water industry and academia. The study evaluated toxicological effects of complex disinfection byproduct (DBP) mixtures, with an emphasis on reproductive and developmental effects that have been associated with DBP exposures in some human epidemiologic studies. This paper describes a new procedure for producing chlorinated drinking water concentrate for animal toxicology experiments, comprehensive identification of100 DBPs, and quantification of 75 priority and regulated DBPs. In the research reported herein, complex mixtures of DBPs were produced by concentrating a natural source water with reverse osmosis membranes, followed by addition of bromide and treatment with chlorine. By concentrating natural organic matter in the source water first and disinfecting with chlorine afterward, DBPs (including volatiles and semivolatiles) were formed and maintained in a water matrix suitable for animal studies. DBP levels in the chlorinated concentrate compared well to those from EPA's Information Collection Rule (ICR) and a nationwide study of priority unregulated DBPs when normalized by total organic carbon (TOC). DBPs were relatively stable over the course of the animal studies (125 days) with multiple chlorination events (every 5-14 days), and a significant portion of total organic halogen was accounted for through a comprehensive identification approach. DBPs quantified included regulated DBPs, priority unregulated DBPs, and additional DBPs targeted by the ICR. Many DBPs are reported for the first time, including previously undetected and unreported haloacids and haloamides. The new concentration procedure not only produced a concentrated drinking water suitable for animal experiments, but also provided a greater TOC concentration factor (136×), enhancing the detection of trace DBPs that are often below detection using conventional approaches.

For evaluation of the adverse health effects associated with exposures to complex chemical mixtures in the environment, the U.S. Environmental Protection Agency (EPA) (2000) states, "if no data are available on the mixture of concern, but health effects data are available on a similar mixture ... a decision must be made whether the mixture on which health effects are available is 'sufficiently' similar to the mixture of concern to permit a risk assessment." This article provides a detailed discussion of statistical considerations for evaluation of the similarity of mixtures. Multivariate statistical procedures are suggested to determine whether individual samples of drinking-water disinfection by-products (DBPs) vary significantly from a group of samples that are considered to be similar. The application of principal components analysis to (1) reduce the dimensionality of the vectors of water samples and (2) permit visualization and statistical comparisons in lower dimensional space is suggested. Formal analysis of variance tests of homogeneity are illustrated. These multivariate statistical procedures are applied to a data set describing samples from multiple water treatment plants. Essential data required for carrying out sensitive analyses include (1) identification and measurement of toxicologically sensitive process input and output characteristics, and (2) estimates of variability within the data to construct statistically efficient estimates and tests.

This article describes the disinfection by-product (DBP) characterization portion of a series of experiments designed for comprehensive chemical and toxicological evaluation of two drinking-water concentrates containing highly complex mixtures of DBPs. This project, called the Four Lab Study, involved the participation of scientists from four laboratories and centers of the U.S. Environmental Protection Agency (EPA) Office of Research and Development, along with collaborators from the water industry and academia, and addressed toxicologic effects of complex DBP mixtures, with an emphasis on reproductive and developmental effects that are associated with DBP exposures in epidemiologic studies. Complex mixtures of DBPs from two different disinfection schemes (chlorination and ozonation/postchlorination) were concentrated successfully, while maintaining a water matrix suitable for animal studies. An array of chlorinated/brominated/iodinated DBPs was created. The DBPs were relatively stable over the course of the animal experiments, and a significant portion of the halogenated DBPs formed in the drinking water was accounted for through a comprehensive qualitative and quantitative identification approach. DBPs quantified included priority DBPs that are not regulated but have been predicted to produce adverse health effects, as well as those currently regulated in the United States and those targeted during implementation of the Information Collection Rule. New by-products were also reported for the first time. These included previously undetected and unreported bromo- and chloroacids, iodinated compounds, bromo- and iodophenols, and bromoalkyltins.

This article describes disinfection of the same source water by two commonly used disinfection treatment scenarios for purposes of subsequent concentration, chemical analysis, and toxicological evaluation. Accompanying articles in this issue of the Journal of Toxicology and Environmental Health describe concentration of these finished waters by reverse osmosis techniques, chemical characterization of the resulting disinfection by-product (DBP) concentrates, in vivo and in vitro toxicological results, and risk assessment methods developed to analyze data from this project. This project, called the "Four Lab Study," involved participation of scientists from four laboratories/centers of the U.S. Environmental Protection Agency Office of Research and Development as well as extramural collaborators from the water industry and academia. One of the two finished waters was prepared by conventional treatment and disinfected by chlorination. The other finished water was also prepared by conventional treatment and disinfected by ozonation followed by chlorination (ozonation/postchlorination). Chlorination conditions of dose, time and temperature were similar for both treatment scenarios, allowing for a comparison. Both finished waters had acceptably low levels of particulates and bacteria, representative pH and chlorine levels, and contained numerous DBP. Known effects of ozonation were observed in that, relative to the water that was chlorinated only, the ozonated/postchlorinated water had lower concentrations of total organic halogen, trihalomethanes (THM), haloacetic acids (HAA), and higher concentrations of bromate, and aldehydes.

Large-scale differential gene expression analysis was used to examine the biological effects of disinfected surface waters on cultured rat hepatocytes. Source water from East Fork Lake (Harsha Lake), a reservoir on the Little Miami River in Ohio, was spiked with iodide and bromide and disinfected by chlorination or ozonation/postchlorination. The chlorinated and ozonated/postchlorinated waters were concentrated, respectively, 136- and 124-fold (full strength) by reverse-osmosis membrane techniques. Volatile disinfection by-products (DBP) lost during concentration were restored to the extent possible. Primary rat hepatocytes were exposed to either full-strength or 1:10 or 1:20 dilutions of the concentrates for 24 h and assayed for cytotoxicity and gene expression alterations. The full-strength concentrates were cytotoxic, whereas the diluted samples exhibited no detectable cytotoxicity. Differential gene expression analysis provided evidence for the underlying causes of the severe cytotoxicity observed in rat hepatocytes treated with the full-strength ozonation/postchlorination concentrate (e.g., cell cycle arrest, metabolic stasis, oxidative stress). Many gene expression responses were shared among the hepatocyte cultures treated with dilutions of the ozonation/ postchlorination and chlorination concentrates. The shift in the character of the response between the full-strength concentrates and the diluted samples indicated a threshold for toxicity. A small subset of gene expression changes was identified that was observed in the response of hepatocytes to peroxisome proliferators, phthalate esters, and haloacetic acids, suggesting a peroxisome proliferative response.

Chemical disinfection of drinking water is a major public health triumph of the 20th century, resulting in significant decreases in morbidity and mortality from waterborne diseases. Disinfection by-products (DBP) are chemicals formed by the reaction of oxidizing disinfectants with inorganic and organic materials in the source water. To address potential health concerns that cannot be answered directly by toxicological research on individual DBPs or defined DBP mixtures, scientists residing within the various organizations of the U.S. Environmental Protection Agency's Office of Research and Development (the National Health and Environmental Effects Research Laboratory, the National Risk Management Research Laboratory, the National Exposure Research Laboratory, and the National Center for Environmental Assessment) engaged in joint investigation of environmentally realistic complex mixtures of DBP. Research on complex mixtures of DBP is motivated by three factors: (a) DBP exposure is ubiquitous to all segments of the population; (b) some positive epidemiologic studies are suggestive of potential developmental, reproductive, or carcinogenic health effects in humans exposed to DBP; and (c) significant amounts of the material that makes up the total organic halide portion of the DBP have not been identified. The goal of the Integrated Disinfection Byproducts Mixtures Research Project (the 4Lab Study) is provision of sound, defensible, experimental data on environmentally relevant mixtures of DBP and an improved estimation of the potential health risks associated with exposure to the mixtures of DBP formed during disinfection of drinking water. A phased research plan was developed and implemented. The present series of articles provides the results from the first series of experiments.

Fixed-ratio ray designs have been used for detecting and characterizing interactions of large numbers of chemicals in combination. Single-chemical dose-response data are used to predict an “additivity curve” along an environmentally relevant ray. A “mixture curve” is estimated from the mixture dose-response data along the ray. A test of additivity is equivalent to a test of coincidence of these two curves, which is based on the traditional hypothesis testing framework that assumes additivity in the null hypothesis and rejects with evidence of interaction. However, failure to reject may be due to lack of statistical power, making the claim of additivity problematic. As a solution we have developed rigorous methodology to test for additivity using statistical equivalence testing logic in which additivity is claimed based on pre-specified biologically important additivity margins, if the data support such a claim. Using the principle of confidence interval inclusion, a confidence region about the difference of meaningful functions of model parameters from the mixture model and that predicted under additivity is computed. When the confidence region is completely contained within the additivity margins then additivity is claimed with a Type I error rate chosen a priori to be some acceptably small value. The method is illustrated using an environmentally relevant fixed-ratio mixture of nine haloacetic acids where cytotoxic response is measured.

As part of the DOE RIA R&D effort we investigated the possibilities and problems of beam strippers in the different heavy ion accelerator components of a possible Rare Isotope Accelerator (RIA) facility. We focused on two beam stripping positions in the RIA heavy ion driver where benchmark currents of up to 5 particle μA 238-U were projected at energies of 10.5 MeV/u and 85 MeV/u respectively. In order to select feasible stripper materials, data from experiments with Uranium beams at Texas A&M and GSI were evaluated. Based on these results thermal estimates for a possible design were calculated and cooling simulations with commercially available software performed. Additionally, we performed simulations with the GEANT4 code on evaluating the radiation environment for our beam stripping solution at the 85 MeV/u position in the RIA driver.

When an interaction has been detected among the chemicals in a mixture, it may be of interest to predict the interaction threshold. A method is presented for estimation of an interaction threshold along a mixture ray which allows differences in the shapes of the dose-response curves of the individual components (e.g., mixtures of full and partial agonists with differing response maxima). A point estimate and confidence interval for the interaction threshold may be estimated. The methods are illustrated with data from a study of a mixture of 18 polyhalogenated aromatic hydrocarbons (PHAHs) in rats exposed by oral gavage for four consecutive days. Serum total thyroxine (T4) was the response variable. Previous analysis of these data demonstrated a dose-dependent interaction among the 18 chemicals in the mixture, with additivity suggested in the lower portion of the dose-response curve and synergy (greater than additive response) in the higher portion of the dose-response curve. The present work builds on this analysis by construction of an interaction threshold model along the mixture ray. This interaction threshold model has two components: an implicit additivity region and an explicit region that describes the departure from additivity; the interaction threshold is the boundary between the two regions. Estimation of the interaction threshold within the observed experimental region suggested evidence of additivity in the low dose region. Total doses of the mixture that exceed the upper limit of the confidence interval on the interaction threshold were associated with a greater-than-additive interaction.