Evaluation of the Interaction-Based Hazard Index Formula Using Data on Four Trihalomethanes from U.S. EPA's Multiple-Purpose Design Study.

The interaction-based hazard index (HIINT), a mixtures approach to characterizing toxicologic interactions, is demonstrated and evaluated by statistically analyzing data on four regulated trihalomethanes (THMs). These THMs were the subject of a multipurpose toxicology study specifically designed to evaluate the HIINT formula. This HIINT evaluation uses single, binary and quaternary mixture THM data. While this research is considered preliminary, the results provide insights on the application of HIINT when toxicology mixture data are available and on improvements to the method. The results for relative liver weight show the HIINT was generally not conservative but did adjust the additive hazard index (HI) in the correct direction, predicting greater than dose-additivity, as seen in the mixture data. For the liver serum enzyme endpoint alanine aminotransferase, the results were mixed, with some indices giving an estimated effective dose lower than the observed mixture effective dose and others higher; in general, the HIINT adjusted the HI in the correct direction, predicting less than dose-additivity. In addition, a methodological improvement was made in the calculation of maximum interaction magnitude. Suggested refinements to the HIINT included mixture-specific replacements for default parameter values and approaches for supplementing the usual qualitative discussions of uncertainty with numerical descriptions.

Evaluation of a Proportional Response Addition Approach to Mixture Risk Assessment and Predictive Toxicology Using Data on Four Trihalomethanes from the U.S. EPA's Multiple-Purpose Design Study.

In this study, proportional response addition (Prop-RA), a model for predicting response from chemical mixture exposure, is demonstrated and evaluated by statistically analyzing data on all possible binary combinations of the four regulated trihalomethanes (THMs). These THMs were the subject of a multipurpose toxicology study specifically designed to evaluate Prop-RA. The experimental design used a set of doses common to all components and mixtures, providing hepatotoxicity data on the four single THMs and the binary combinations. In Prop-RA, the contribution of each component to mixture toxicity is proportional to its fraction in the mixture based on its response at the total mixture dose. The primary analysis consisted of 160 evaluations. Statistically significant departures from the Prop-RA prediction were found for seven evaluations, with three predications that were greater than and four that were less than the predicted response; interaction magnitudes (n-fold difference in response vs. prediction) ranged from 1.3 to 1.4 for the former and 2.6 to 3.8 for the latter. These predictions support the idea that Prop-RA works best with chemicals where the effective dose ranges overlap. Prop-RA does not assume the similarity of toxic action or independence, but it can be applied to a mixture of components that affect the same organ/system, with perhaps unknown toxic modes of action.

Environmental exposure to metals and the development of tauopathies, synucleinopathies, and TDP-43 proteinopathies: A systematic evidence map protocol.

BACKGROUND: Neurodegenerative diseases such as Alzheimer's disease, Parkinson's disease, and amyotrophic lateral sclerosis are incurable and expected to increase in prevalence in the upcoming decades. Environmental exposure to metals has been suggested as a contributing factor to the development of neurodegenerative disease. This systematic evidence map will identify and characterize the epidemiological and experimental data available on the intersection of eighteen metals of environmental concern (i.e., aluminum, antimony, arsenic, barium, beryllium, cadmium, chromium, cobalt, copper, lead, manganese, mercury, nickel, palladium, radium, silver, vanadium, and zinc) and three neurodegenerative disease clusters (i.e., tauopathies, synucleinopathies, and TDP-43 proteinopathies). We aim to describe the type and amount of evidence available (or lack thereof) for each metal and neurodegenerative disease combination and highlight important knowledge gaps and knowledge clusters for future research. METHODS: We will conduct a thorough search using two databases (MEDLINE and Web of Science Core Collection) and grey literature resources. Pre-defined criteria have been developed to identify studies which evaluate at least one of the selected metals and neurodegenerative disease-relevant outcomes (e.g., neuropathology, cognitive function, motor function, disease mortality). At each phase of review, studies will be evaluated by two reviewers. Studies determined to be relevant will be extracted for population, exposure, and outcome information. We will conduct a narrative review of the included studies, and the extracted data will be available in a database hosted on Tableau Public. CONCLUSION: This protocol documents the decisions made a priori to data collection regarding these objectives.

Genomic comparisons between hepatocarcinogenic and non-hepatocarcinogenic organophosphate insecticides in the mouse liver.

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.

Censoring Trace-Level Environmental Data: Statistical Analysis Considerations to Limit Bias.

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.

Analysis of proportional data in reproductive and developmental toxicity studies: Comparison of sensitivities of logit transformation, arcsine square root transformation, and nonparametric analysis.

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.

Global optimization of the Michaelis-Menten parameters using physiologically-based pharmacokinetic (PBPK) modeling and chloroform vapor uptake data in F344 rats.

Objective: To quantify metabolism, a physiologically based pharmacokinetic (PBPK) model for a volatile compound can be calibrated with the 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.Materials and methods: 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 64 bit, 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.Results and discussion: The final values estimated for Vmax (maximum metabolic rate) and Km (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 Vmax and Km 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 Km.Conclusion: 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.

Predictive Analysis Using Chemical-Gene Interaction Networks Consistent with Observed Endocrine Activity and Mutagenicity of U.S. Streams.

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.

Per- and polyfluoroalkyl substances in source and treated drinking waters of the United States.

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 <1 ng/L to 1102 ng/L. The median total PFAS concentration was 21.4 ng/L in the source water and 19.5 ng/L in the treated drinking water. Comparing the total PFAS concentration in source and treated water at each location, only five locations demonstrated statistically significant differences (i.e. P < 0.05) between the source and treated water. When the perfluorooctanoic acid (PFOA) and perfluorooctanesulfonic acid (PFOS) concentrations in the treated drinking water are compared to the existing US Environmental Protection Agency's PFOA and PFOS drinking water heath advisory of 70 ng/L for each chemical or their sum one DWTP exceeded the threshold. Six of the 25 DWTPs were along two large rivers. The DWTPs within each of the river systems had specific PFAS profiles, with the three DWTPs from one river being dominated by PFOA, while three DWTPs on the second river were dominated by perfluorobutyric acid (PFBA).

Chlorination of Source Water Containing Iodinated X-ray Contrast Media: Mutagenicity and Identification of New Iodinated Disinfection Byproducts.

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.

Evaluating In Vitro-In Vivo Extrapolation of Toxicokinetics.

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 under-estimated 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.

From the Cover: Genomic Effects of Androstenedione and Sex-Specific Liver Cancer Susceptibility in Mice.

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.

Method to assess component contribution to toxicity of complex mixtures: Assessment of puberty acquisition in rats exposed to disinfection byproducts.

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.

Human health screening and public health significance of contaminants of emerging concern detected in public water supplies.

The source water and treated drinking water from twenty five drinking water treatment plants (DWTPs) across the United States were sampled in 2010-2012. Samples were analyzed for 247 contaminants using 15 chemical and microbiological methods. Most of these contaminants are not regulated currently either in drinking water or in discharges to ambient water by the U. S. Environmental Protection Agency (USEPA) or other U.S. regulatory agencies. This analysis shows that there is little public health concern for most of the contaminants detected in treated water from the 25 DWTPs participating in this study. For vanadium, the calculated Margin of Exposure (MOE) was less than the screening MOE in two DWTPs. For silicon, the calculated MOE was less than the screening MOE in one DWTP. Additional study, for example a national survey may be needed to determine the number of people ingesting vanadium and silicon above a level of concern. In addition, the concentrations of lithium found in treated water from several DWTPs are within the range previous research has suggested to have a human health effect. Additional investigation of this issue is necessary. Finally, new toxicological data suggest that exposure to manganese at levels in public water supplies may present a public health concern which will require a robust assessment of this information.

Nationwide reconnaissance of contaminants of emerging concern in source and treated drinking waters of the United States.

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 from 29 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.

Dose and Effect Thresholds for Early Key Events in a PPARα-Mediated Mode of Action.

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.

Occurrence and Comparative Toxicity of Haloacetaldehyde Disinfection Byproducts in Drinking Water.

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 haloacetaldehydes (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.

Raising the bar for reproducible science at the U.S. Environmental Protection Agency Office of Research and Development.

Considerable concern has been raised regarding research reproducibility both within and outside the scientific community. Several factors possibly contribute to a lack of reproducibility, including a failure to adequately employ statistical considerations during study design, bias in sample selection or subject recruitment, errors in developing data inclusion/exclusion criteria, and flawed statistical analysis. To address some of these issues, several publishers have developed checklists that authors must complete. Others have either enhanced statistical expertise on existing editorial boards, or formed distinct statistics editorial boards. Although the U.S. Environmental Protection Agency, Office of Research and Development, already has a strong Quality Assurance Program, an initiative was undertaken to further strengthen statistics consideration and other factors in study design and also to ensure these same factors are evaluated during the review and approval of study protocols. To raise awareness of the importance of statistical issues and provide a forum for robust discussion, a Community of Practice for Statistics was formed in January 2014. In addition, three working groups were established to develop a series of questions or criteria that should be considered when designing or reviewing experimental, observational, or modeling focused research. This article describes the process used to develop these study design guidance documents, their contents, how they are being employed by the Agency's research enterprise, and expected benefits to Agency science. The process and guidance documents presented here may be of utility for any research enterprise interested in enhancing the reproducibility of its science.

Reproductive toxicity of a mixture of regulated drinking-water disinfection by-products in a multigenerational rat bioassay.

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.