Transcriptional effects of binary combinations of PFAS in FaO cells.

Per- and polyfluoroalkyl substances (PFAS) represent a large class of structurally diverse chemicals of increasing public concern, mostly due to their chemical stability and undetermined toxicity profiles. In laboratory animals, adverse effects implicated for certain PFAS, perfluorooctanoic acid (PFOA) and perfluorooctane sulfonate (PFOS) in particular, include liver toxicity and the associated metabolic dysregulation, immune and thyroid alterations, reproductive toxicity, and selected tumors. The broad commercialization and environmental distribution of PFAS has drawn attention to the need for understanding risks associated with combined exposure to multiple PFAS in complex mixtures. The purpose of this investigation is to determine whether binary combinations of PFAS elicit a molecular response that is either greater than or less than the sum of the individual responses. Exposure of FaO rat hepatoma cells for 24 h to 25 µM-200 µM of the 4- and 8-carbon perfluorocarboxylic acids (PFBA and PFOA) or the 4, 6, and 8-carbon perfluorosulfonic acids (PFBS, PFHxS, and PFOS, respectively) individually caused a dose-dependent increase in PPARα-regulated expression of peroxisomal bifunctional enzyme (Ehhadh). Potency increased with carbon number, with the carboxylates eliciting a greater transcriptional response than the corresponding sulfonates. Combined exposure to PFOA and PFBA produced an effect that was significantly less than the sum of the individual responses. The response to the combination of PFOA and PFOS produced a summative effect at concentrations that were not cytotoxic. Combined exposures to PFOS and either PFBS or PFHxS at low noncytotoxic concentrations produced a transcriptional effect that was significantly less than the sum of the individual effects. The results demonstrate that among the five structurally related perfluoroalkyl acids included in this investigation, PPARα transcriptional activation in response to combined binary exposures is consistently at or below that predicted by the sum of the individual effects.

New models for the time dependent toxicity of individual and combined toxicants.

Time-dependent toxicity data of specific toxicants observed against A. fischeri were analysed using a single time-dependent Weibull cumulative distribution function (CDF). For binary mixtures the individual Weibull parameters provided the marginals for a bivariate copula model, characterised by a single extra parameter, the interaction exponent, θ. The copula model unites both the dose addition (DA, Loewe additivity) and dose independence (DI, Bliss independence) hypotheses of combinations into a single explicit equation and returns both hypotheses as special cases. The model predicts the linear isoboles from sham (like against like) experiments and the linear, concave, convex and mixed concave-convex isoboles from true binary mixtures. Systems are defined as being independent when θ = 1, additive when min(n 1, n 2) ≤ θ ≤ max(n 1, n 2), antagonistic when θ < min(n 1, n 2) and synergistic when θ > max(n 1, n 2), where n i are the individual dose (concentration) exponents. More complex mixtures were analysed by developing n-dimensional copulas: two ternary systems were analysed using a ternary copula that returned the three bivariate marginals and the three individual marginals of the mixture. The general model can also be used for time-independent studies by simply removing the time dependency. Research into combined effects has often assumed that if both the parameters for all the individual components in a mixture and the model for combinations were known then the additive effect of the whole could be predicted. This hypothesis has been shown to be false because without knowing how the components in a mixture interact the predictions from the standard DI or DA models provide only an initial, best guess, analysis.

Evaluation of time-dependent toxicity and combined effects for a series of mono-halogenated acetonitrile-containing binary mixtures.

Mixture and time-dependent toxicity (TDT) was assessed for a series of mono-halogenated acetonitrile-containing combinations. Inhibition of bioluminescence in Aliivibrio fischeri was measured after 15, 30 and 45-min of exposure. Concentration-response (x/y) curves were determined for each chemical alone at each timepoint, and used to develop predicted x/y curves for the dose-addition and independence models of combined effect. The x/y data for each binary mixture was then evaluated against the predicted mixture curves. Two metrics of mixture toxicity were calculated per combined effect model: (1) an EC50-based dose-addition (AQ) or independence (IQ) quotient and (2) the mixture/dose-addition (MX/DA) and mixture/independence (MX/I) metrics. For each single chemical and mixture tested, TDT was also calculated. After 45-min of exposure, 25 of 67 mixtures produced curves that were consistent with dose-addition using the MX/DA metric, with the other 42 being less toxic than predicted by MX/DA. Some mixtures had toxicity that was consistent with both dose-addition and independence. In general, those that were less toxic than predicted for dose-addition were also less toxic than predicted for independence. Of the 25 combinations that were consistent with dose-addition, 22 (88%) mixtures contained chemicals for which the individual TDT values were both >80%. In contrast, of the 42 non-dose-additive combinations, only 2 (4.8%) of the mixtures had both chemicals with individual TDT values >80%. The results support previous findings that TDT determinations can be useful for predicting chemical mixture toxicity.

Time-dependence in mixture toxicity prediction.

The value of time-dependent toxicity (TDT) data in predicting mixture toxicity was examined. Single chemical (A and B) and mixture (A+B) toxicity tests using Microtox(®) were conducted with inhibition of bioluminescence (Vibrio fischeri) being quantified after 15, 30 and 45-min of exposure. Single chemical and mixture tests for 25 sham (A1:A2) and 125 true (A:B) combinations had a minimum of seven duplicated concentrations with a duplicated control treatment for each test. Concentration/response (x/y) data were fitted to sigmoid curves using the five-parameter logistic minus one parameter (5PL-1P) function, from which slope, EC25, EC50, EC75, asymmetry, maximum effect, and r(2) values were obtained for each chemical and mixture at each exposure duration. Toxicity data were used to calculate percentage-based TDT values for each individual chemical and mixture of each combination. Predicted TDT values for each mixture were calculated by averaging the TDT values of the individual components and regressed against the observed TDT values obtained in testing, resulting in strong correlations for both sham (r(2)=0.989, n=25) and true mixtures (r(2)=0.944, n=125). Additionally, regression analyses confirmed that observed mixture TDT values calculated for the 50% effect level were somewhat better correlated with predicted mixture TDT values than at the 25 and 75% effect levels. Single chemical and mixture TDT values were classified into five levels in order to discern trends. The results suggested that the ability to predict mixture TDT by averaging the TDT of the single agents was modestly reduced when one agent of the combination had a positive TDT value and the other had a minimal or negative TDT value.

Evaluation of an asymmetry parameter for curve-fitting in single-chemical and mixture toxicity assessment.

In mixture toxicity, concentration-effect data are often used to generate conclusions on combined effect. While models of combined effect are available for such assessments, proper fitting of the data is critical to obtaining accurate conclusions. In this study an asymmetry parameter (s) was evaluated for data-fitting and compared with our previous approach. Inhibition of bioluminescence was assessed with Vibrio fischeri at 15, 30 and 45-min of exposure with seven or eight concentrations and a control (each duplicated) for each single-chemical (A or B) and mixture (A:B). Concentration-effect data were fitted to sigmoid curves using the four-parameter logistic function (4PL) and the five-parameter logistic minus one-parameter (5PL-1P) function. For the 4PL, parameters included minimum effect, maximum effect, EC(50) and slope, while for the 5PL-1P the minimum effect parameter was removed and an asymmetry parameter was added. A total of 72 mixture toxicity data sets were evaluated, representing 432 single-chemical and 216 mixture curves. Mean coefficients of determination (r(2)) for all 648 curves showed that the 5PL-1P gave better fitting (0.9982 ± 0.0018) than the 4PL (0.9973 ± 0.0030). For both functions, the sum-of-squares of the residuals (SS-Res) was determined for each curve. The 5-parameter rational regression best described the relationship between the decrease in sum-of-squares of the residuals (i.e., 4PL: SS-Res - 5PL-1P: SS-Res) and log s, with fitting improved the most at low values of s (s<0.8). This held even when curves with r(2) values ≤ 0.9970 were removed from the analyses. Subsequent review of the combined effects obtained via the 4PL and the 5PL-1P functions resulted in a change in the interpretation of combined effect in 39/216 (18%) cases.

Time-dependence in mixture toxicity with soft-electrophiles: 2. Effects of relative reactivity level on time-dependent toxicity and combined effects for selected Michael acceptors.

Toxicity assessments for organic chemical mixtures are often described as being approximately additive. Recent mixture studies with soft electrophiles have suggested that agents with less-than fully time-dependent toxicity (TDT) may actually induce toxicity by more than one mode of toxic action within the same series of concentrations. To evaluate this concept further, four Michael acceptor electrophiles, each with a different rate of in chemico reactivity and different level of TDT, were tested with each other and in sham combinations (a single chemical tested as if it were a binary mixture) using the Microtox system. For each binary combination, each agent was tested alone and in a mixture, with toxicity assessed as inhibition of bioluminescence at 15-, 30- and 45-min of exposure. Each single agent and mixture test included seven duplicated concentrations and a duplicated control treatment. To evaluate relative reactivity, each agent was also tested with the model nucleophile glutathione (GSH). Agents with greater in chemico reactivity (mean RC(50) mM) showed greater toxicity (mean 45-min EC(50) - mM) but these were inversely related to the TDT levels of the agents. Combined effects for the sham combinations, as quantified by additivity quotient values for the EC(50) of the mixture, tended to be close to 1.00 (i.e., the dose-addition EC(50)-AQ). For true binary combinations (i.e., two chemicals tested together), the EC(50)-AQ tended to be increasingly above 1.00 when TDT levels of the agents in the mixture were more disparate. The results of this study with Michael acceptors suggested that: (i) when reactivity was fast, there was most likely a single prominent mode of toxic action, i.e., electro(nucleo)philic reactivity, leading to time-dependent toxicity at the full or high levels, (ii) when the reaction rate for a chemical was slower, two modes of action, electro(nucleo)philic reactivity and narcosis, were apparent such that the time-dependent toxicity level was lower as well, (iii) mixtures of the former agents show a combined effect that was strictly dose-additive, whereas (iv) mixtures which included one (or more) agent with a lower reaction rate had a combined effect that was approximately additive rather than strictly dose-additive.

Chemical mixture toxicity testing with Vibrio fischeri: combined effects of binary mixtures for ten soft electrophiles.

The toxicity of 30 binary combinations of 10 soft electrophiles was examined in Microtox using dose-response curve (DRC) analysis. Chemicals from three groups of soft electrophiles-vinyl Michael acceptors (I--react with a thiol group), dicarbonyl reactive agents (II--react with a primary amine), and alpha-haloactivation compounds (III--react with a thiol group)--were selected for testing to evaluate the relationship between molecular site of chemical action and combined toxic effect. For each combination tested, each single agent was tested alone at six duplicated concentrations and three 1:1 mixtures of the agents were also tested, each at six duplicated concentrations. Exposure duration was 15 min for each single agent and mixture test. Sigmoid DRCs for each single chemical and mixture were constructed and the single chemical curves were used to develop a theoretical dose-addition DRC for the combination. Additivity quotient (AQ) values for slope and EC50 were calculated by dividing the actual mixture slope or EC50 for a given combination by the predicted slope or EC50, respectively, from the theoretical dose-addition DRC. Three criteria were selected for value in determining the combined effect obtained for each combination: (1) slope AQ 95% confidence interval (CI) overlap with 1.0 (1.0=dose addition), (2) EC50 AQ 95% CI overlap with 1.0, and (3) mean mixture data point 95% and 99% CI overlap with the theoretical dose-addition DRC. Each of three sham combinations showed combined effects consistent with dose addition for each criterion. Dose addition was expected for 15 nonsham combinations (nine within-group combinations and six group I:III combinations) and a nondose-additive effect was expected for 12 combinations (all I:II and II:III combinations). Actual combined effects obtained by incorporating all three criteria (noted above) showed only six instances of dose addition. Therefore, time-dependent toxicity (TDT) tests of each soft electrophile alone and for three nonpolar narcotic chemicals alone were conducted, using 15-, 30-, and 45-min exposure durations, to assess the time-dependent nature of the toxicity. Results of the TDT tests suggested that five had fully (or nearly fully) TDT (interpreted as an irreversible effect representing one molecular site of action), five of the soft electrophiles had partially TDT (i.e., representing two or more molecular sites of action for the agents, one irreversible and one reversible), and the three nonpolar narcotics had no TDT (i.e., a fully reversible toxic effect). With this TDT information, the combined effects for 25 of the 27 mixtures, although rather complex, could be explained. It is noteworthy that all combined effects obtained, whether concluded to be dose-additive or not, were close to dose-additive for hazard assessment purposes.

Evaluation of dose-response curve analysis in delineating shared or different molecular sites of action for osteolathyrogens.

Single-chemical and mixture concentration-response curves generated using a frog embryo model were examined for value in assessing whether chemicals exert toxic effects at the same or at different molecular sites of action. Toxicity tests were conducted on a series of osteolathyrogens, i.e. chemicals that inhibit cross-linking of developing connective tissue fibers. Induction of osteolathyrism, which manifests as lesions in the notochord of exposed tadpoles, has several possible molecular sites of action, including agent-cofactor reactivity during the enzyme-mediated cross-linking process. UV-VIS spectrophotometry of osteolathyrogen-cofactor reactivity (i.e. in vitro analysis) was coupled with the 96-h frog embryo mixture toxicity assay (i.e. in vivo toxicity) to compare molecular sites of action for several osteolathyrogens with the combined osteolathyritic effects of the agents. Single-chemical concentration-response curves were used to calculate theoretical curves for the dose-addition model of combined effect. Slope and EC(50) values for both theoretical and experimental mixture curves were then generated to statistically examine the hypothesis that agents with shared sites of action have dose-response curve (DRC) slopes that are similar when given alone and in combination, and slope and EC(50) values that, when administered together, are consistent with those calculated for dose-addition. For combinations of cofactor-binding agents (semicarbazide, thiosemicarbazide, aminoacetonitrile), slope values were generally similar with additivity quotients near 1.0 (1.0=dose-additive) and combined osteolathyritic effects that were consistent with dose-addition. None of these were true for combinations that included agents that did not show rapid cofactor binding (ß-aminopropionitrile, methyleneaminoacetonitrile). The results suggest that DRC analysis could be a useful tool for delineating common or different molecular sites of toxic action and that the approaches used warrant further study for evaluating the mechanistic basis for combined effects of toxicants.

Housing and husbandry of Xenopus for oocyte production.

Despite their importance as a research model, particularly in developmental toxicology investigations, there are few established standards for maintaining Xenopus spp. frogs in the laboratory. The authors review the literature on handling, housing, nutrition, and breeding of Xenopus spp. for optimal oocyte production.

Biochemical and toxicological evaluation of agent-cofactor reactivity as a mechanism of action for osteolathyrism.

In vitro reactivity for each of four osteolathyrogens with a model compound for the lysyl oxidase (LO) cofactor was evaluated and coupled with mixture toxicity testing to evaluate agent-cofactor reactivity as a potential mechanism of action for osteolathyrism. Reactivity of the model cofactor (mLTQ: 4-butylamino-5-methyl-o-quinone), with each of two ureides, semicarbazide (SC) and thiosemicarbazide (TSC), and each of two aminonitriles, aminoacetonitrile (AAN) and beta-aminopropionitrile (betaAPN), was assessed using UV-vis spectrophotometry; both in the absence and presence of Cu(II)-bipyridine (bipy) complex. Two sets of mixture toxicity experiments were conducted using a frog embryo assay that assessed the incidence of osteolathyrism in the notochord of tadpoles after 96-h exposure. The resulting concentration-response curves for each set were evaluated (chi(2) goodness-of-fit test) against theoretical curves for two combined effects models: dose-addition and independence, to determine the combined effect of each osteolathyrogen combination. The agents SC, TSC and AAN each showed rapid, irreversible reactivity with mLTQ, both in the absence and presence of Cu(II)-bipy complex, as indicated by bleaching of the mLTQ peak (504 nm) and formation of an adduct at 350 nm. betaAPN showed no apparent reactivity in the absence of prolonged incubation with mLTQ, whether Cu(II)-bipy complex was present or not. After prolonged incubation (24-144 h) a new peak formed at 350 nm, suggesting that betaAPN reacts weakly with the cofactor, but in a manner different from the other agents examined. The toxicity tests indicated a dose-additive combined effect for the SC:TSC, AAN:SC and AAN:SC:TSC mixtures (0.1