REVIEW article

Front. Microbiol., 31 July 2020

Sec. Microbiotechnology

Volume 11 - 2020 | https://doi.org/10.3389/fmicb.2020.01651

Mathematical Modeling Approaches for Assessing the Joint Toxicity of Chemical Mixtures Based on Luminescent Bacteria: A Systematic Review

  • 1. Department of Environmental Science and Engineering, Xi'an Jiaotong University, Shaanxi, China

  • 2. Department of Structures and Environmental Engineering, University of Agriculture, Faisalabad, Pakistan

Abstract

Developments in industrial applications inevitably accelerate the discharge of enormous substances into the environment, whereas multi-component mixtures commonly cause joint toxicity which is distinct from the simple sum of independent effect. Thus, ecotoxicological assessment, by luminescent bioassays has recently brought increasing attention to overcome the environmental risks. Based on the above viewpoint, this review included a brief introduction to the occurrence and characteristics of toxic bioassay based on the luminescent bacteria. In order to assess the environmental risk of mixtures, a series of models for the prediction of the joint effect of multi-component mixtures have been summarized and discussed in-depth. Among them, Quantitative Structure-Activity Relationship (QSAR) method which was widely applied in silico has been described in detail. Furthermore, the reported potential mechanisms of joint toxicity on the luminescent bacteria were also overviewed, including the Trojan-horse type mechanism, funnel hypothesis, and fishing hypothesis. The future perspectives toward the development and application of toxicity assessment based on luminescent bacteria were proposed.

Introduction

With the intensive industrial emissions and human activities, natural compounds and anthropogenic pollutants discharged into the environment accumulate in ecosystems, which contributed a lot to the environmental pollution. The toxicity of the pollutants, especially multi-component mixtures in a single sample are of more considerable concern due to their potential negative impact on human and other living organisms, such as modifying the stability of genetic information and inhibition of growth. The multi-component contaminants may have an adverse effect on organism even the concentration of single pollutant is less than No Observable Effect Concentration (NOEC) (Hageman et al., 2019). Furthermore, when two or more component mixtures act on a living organism simultaneously, the joint toxicity may significantly differ from the single-substance toxicity.

The sustained attention on the toxic effects of pollutants in the environment has brought forth the rapid development of sensitive and controllable monitoring techniques. In comparison with the physical and chemical analysis, that could only determine the compositions and concentration of the pollutants in the environmental samples, the biological assay could intuitively reflect the detrimental effects of contaminants toward target organisms which is of great importance and contribute to revealing the mechanisms of the toxic effects (Li et al., 2013a,b). Among various bioassay types, acute toxicity bioassays have been rapidly developing as a preferred tool for environmental risk assessment in recent decades which can detect the adverse effects of a single component or multi-component mixtures in a short period and determine the relationship between dose and effect of toxic substances. Acute toxicity bioassays generally can be conducted with both prokaryote and eukaryote, containing bacteria (Ma et al., 2018; Rao et al., 2019; Feng et al., 2020), algae (Costa et al., 2018; Staveley et al., 2018), invertebrates (Chen et al., 2020; Sackey et al., 2020; Vilela et al., 2020), vertebrates (Clemow and Wilkie, 2015; Thornton et al., 2017; Barron et al., 2018; Ivey et al., 2019), and so on. Inevitably, the application of multicellular organisms exhibits the undesirable limitations of relatively longer operation time and higher cost. Therefore, microbial bioassays, especially with luminescent bacteria increased widely interests in the acute toxicity bioassays due to its operational simplicity, less time and cost consuming, higher sensitivity, specificity, and reproducibility (Conforti et al., 2008; Skotti et al., 2014; Yang et al., 2016).

With the growing number of disparate predictions available to perform multi-component mixtures biotoxicity analyses with luminescent bacteria, the choice of methods presents a major challenge to researchers as there is no review of comparative methods currently available. In this review, we summarize the occurrence of bioluminescent assays and the prevalence of determination of joint toxicity of multi-component mixtures and describes the various types of models involved in joint toxicity prediction. The representative methods for the prediction and evaluation of acute joint toxicity are overviewed including graphical method, joint effect indexes, and integrative models. Furthermore, the QSAR method as a hotspot in joint toxicity researches, are also presented in detail.

The Physiological Mechanism and Development of Toxic Bioassay by Luminescent Bacteria

Currently, luminescent bacteria have been applied worldwide in the bioassays of toxicity assessment, including Aliivibrio fisheri (A. fischeri, previously identified as Vibrio fisheri or Photobacterium fisheri), Photobacterium phosphoreum, Vibrio qinghaiensis, Vibrio harveyi and so on (Urbanczyk et al., 2007). Luminescent bacteria that could emit visible bioluminescence under certain physiological conditions were primarily isolated from the ocean environment except for V. cholerea and V. qinghaiensis sp. Nov (Zhu et al., 1994; Ma et al., 2014).

Physiological Mechanism of Bioluminescence From Luminescent Bacteria

The bioluminescence mechanism of luminescent bacteria has been fully studied and characterized (Inouye, 1994; Moore and James, 1995; Bourgois et al., 2001; Vetrova et al., 2007; Nijvipakul et al., 2010; Bergner et al., 2015; Brodl et al., 2018; Tanet et al., 2019). Generally, luciferase, reduced flavin mononucleotide (FMNH2), oxygen (O2), and long-chain fatty aldehyde (RCHO) are prerequisite for the luminescence emitted by luminescent bacteria in which flavin mononucleotide (FMN) was reduced to FMNH2 by the catalyzation of NAD(P)H: FMN oxidoreductase. Then FMNH2 and molecular oxygen will bind with luciferase which has flavin reductase activity to form the intermediate. Subsequently, the intermediate will decompose in the presence of aldehyde, FMNH2 are oxidized into FMN and H2O by the catalyzation of luciferase enzyme, along with the emitting of continuous blue-green light (490 nm) (Figure 1) (Inouye, 1994).

Figure 1

The enzyme involved in bioluminescence in luminescent bacteria is commonly encoded by lux genes which are usually clustered together as lux operon (luxCDABE) (Matheson and Lee, 1981). Luciferase possesses α and β subunits, that are encoded by luxA and luxB genes, respectively. Moreover, the α subunit which exhibits major catalytic activity is similar to the β subunit in amino acid sequence (around 30%) (Sharifian et al., 2017). In general, the organization of lux genes and the amino acid sequences from different luminescent bacteria are similar which indicates their origin from common ancestor and conservation in evolution. Fatty acid reductase complex is encoded by luxC, luxD, and luxE, which is responsible for the transformation between long-chain fatty acids and fatty aldehydes (Figure 1). In certain luminescent bacteria such as Aliivibrio species, the presence of regulatory genes luxI and luxR are detected which link bioluminescence with quorum sensing (QS) (Miyashiro and Ruby, 2012). The gene luxI encodes the enzyme required for the synthesis of the autoinducer molecule, and LuxR is the autoinducer-dependent transcription factor, the binding of autoinducer-LuxR to the promoter will initiate the transcription of lux operon (Engebrecht and Silverman, 1984; Fuqua et al., 1994). Some additional lux genes including luxG, luxF and luxH have also been found in some luminescent lux operons, but there is no evidence to confirm that these genes were directly related to bioluminescence (O'Grady and Wimpee, 2008; Bergner et al., 2015; Deeva et al., 2019).

Development Application of Toxicity Assessment by Luminescent Bacteria

The toxicity bioassay with luminescent bacteria has been successfully employed in the comprehensive ecotoxicity assessments. In 1978, Beckman Company of the United States first commercialized the toxic bioassay as “Microtox®” with A. fischeri, the sensitivity of which was comparable to the 96 h acute toxicity test of fish (Flokstra et al., 2008). Since then, the toxicity assessment by luminescent bacteria has been developing rapidly toward higher sensitivity, accuracy, reproducibility, anti-interference, and extensive application. As a comprehensive indicator, luminescent bacteria have been successfully applied for monitoring the ecotoxicity of the natural environment, including wastewater, river water, sewage sludge, landfill leachate, and so on (Jarque et al., 2016; Pino-Otín et al., 2019). Zhang et al., examined the toxicity of surface water from the Huangpu River of China by V. qinghaiensis sp.-Q67, and found the seasonal variation of toxicity (Zhang et al., 2015). Likewise, Jarque et al. held the view that V. fishieri luminescence assay to evaluate the toxicity associated with sediments from the Morava River and its tributary Drevnice River followed seasonal patterns (Jarque et al., 2016). Luminescent bacteria toxicity bioassay could also be carried out to evaluate the ecological risk of wastewater before and after diverse treatment processes, including ozonation, Fenton, UV, electrochemical degradation, biological process, and so on (De Schepper et al., 2012; Hamdi El Najjar et al., 2014; Lütke Eversloh et al., 2014; Wei et al., 2015; Nebout et al., 2016; Xu et al., 2018).

No confined to aqueous bodies, luminescent bacteria bioassay are also conducted in another environmental medium. Aammi et al. evaluated the acute toxicity and genotoxicity levels of atmospheric particulate matter (PM) samples from three districts in Istanbul, Turkey. V.fischeri was used to test the biotoxicity of the collected ambient air samples. The study identified significant locational and seasonal differences in the toxicity level of PM2.5−10 extracts. In terms of seasonal differences, more toxic and anthropogenic particle samples present in winter. Most of the toxic PM samples were identified in central urban, it can be attributed mainly to the inner-city and transit traffic emissions (Aammi et al., 2017). Moreover, Microtox® system along with V.fischeri has also been applied to determine the ecotoxicity of PM emitted from three different light-duty vehicles, the results demonstrated that PM from both gasoline and diesel vehicle equipped with the silicon-carbide catalyzed diesel particle filter (CDPF) appeared higher toxicity on a per mass basis than the corresponding vehicle with conventional exhaust, which indicated the decrease of PM emission level was not equated to the reduction of toxicity (Vouitsis et al., 2009). To reveal the influence factors on the formation of environmental persistent free radicals (EPFRs) in soil and the associated biotoxicity, Zhang et al. mimicked the formation of the EPFRs in the catechol-contaminated soil, they concluded that significantly inhibited the bioluminescence of P. phosphoreum when the concentration of EPFRs reached the maximum at pyrolysis temperature of 300. It could be speculated that the production of ·OH or reactive oxygen species (ROS) may be the dominant contributor to its acute toxicity to luminescent bacteria (Zhang et al., 2019b).

The conventional biotoxicity evaluation methods based on luminescent bacteria suffered from only being able to detect the overall toxicity without selectivity. A great deal of previous research has focused on specific toxins toxicity evaluation of target substances to overcome the drawback. Surveys such as that conducted by Huang et al. (2019a) have shown that dual detection by a combination of luminescent bacteria and aggregation aggregation-induced emission (AIE) probe could detect bioaccumulated Hg2+. The toxic effect of Hg2+ on P. phosphoreum give rise to quenching its bioluminescence inside the bacteria. Moreover, the AIE-active probe 2-AFN-I entered the damaged bacteria and selectively bound to Hg2+, which would turn-on the fluorescence and complete the overall dual detection (Huang et al., 2019a). What's more, microbial biosensors as one of the most significant applications of luminescent bacteria has been developing gradually with the incomparable advantages of selective identification and quantification of specific substances directly in air and water. Microbial biosensors are usually composed of immobilized microorganisms (recognition element) in the transducer to produce a measurable signal such as current, potential, or optical a signal, which can be further amplified and analyzed. Since 1997, Boyd and Mowat confirmed the change of luminescence emitted by living microorganisms in the presence of a target analyte (Boyd et al., 1997), bioluminescence genes (lux) have been recognized as an effective and preferable reporter system of biosensors for toxicity evaluation. The luxCDABE genes have been transferred into specific bioreporter bacteria under the control of a particular promoter to integrate into the recognition element of a biosensor. These recombinant luminescent bacteria can be divided into two categories according to expression pattern (Figure 2). If the lux genes are transferred into the constitutive expression system, the luminescence will be emitted continuously and decrease with the increasing quantity of toxic substances (“lights off”), which is non-specific and suitable for the total toxicity assessment of the mixture. Most natural bioluminescent bacteria belong to this group (Belkin, 2003; Gu et al., 2004). In the other group, lux genes are fused with inducible promoters which can be activated under the exposure of specific chemicals. Therefore, the produced luminescence will increase with the increasing amount of specific toxic chemicals (“lights on”). Most of the microbial biosensors based on luminescent bacteria could be classified as this specific group (Belkin, 2003). Nowadays, biosensors are rapidly developing toward higher sensitivity, specificity, reproducibility, automation and high throughput (Conforti et al., 2008; Skotti et al., 2014; Yang et al., 2016). Furthermore, some specific biosensors with recombinant bioluminescent bacteria have also been designed to detect cell damage, like DNA damage, protein damage, cell membrane damage and stress response (Ahn et al., 2009; Sazykin et al., 2016; Jiang et al., 2017).

Figure 2

Evaluation and Prediction of Joint Toxicity of Mixed Substances

The multi-component mixtures could cause combined toxicity which is distinct from independent effect, due to the apparent differences between the toxic mode of actions (TMOA) and the interactions among constituents. Even if the characteristics of each component are acquainted, the integrative effect of the mixture could not be directly inferred. Therefore, the evaluation and prediction of the joint effect of mixed contaminants have been attracting increasing interests. In the research field of environmental toxicology, the most reported studies on joint toxicity focused on the evaluation of the combined effects of components in the mixed system (Monosson, 2005). According to the different interaction modes of the constituents in a mixture, combined toxicity can be categorized into diverse effects, including antagonism (ANT), additive (ADD), independent (IND), and synergistic (SYN) effects (Plackett and Hewlett, 1952). With the development of toxicology and mathematics, the existing qualitative evaluation methods neither meet the requirements of determining the toxicity contribution of a single component to mixed pollutants nor meet the requirements of predicting the toxicity of a mixture of specific known components. Therefore, more efforts have been made on the quantitative assessments and predictions of joint toxicity of multi-component mixture (Altenburger et al., 2003). Generally, the combined toxicity of the multi-component was determined using a fixed ratio design. Always, the EC50 value of each component within the mixture is constituted a fixed ratio design to evaluate the joint effect. The EC50 value is defined as the concentration that produces a 50% decrease in light emission, obtained from the concentration-response relationship, which could calculate from regression models (Table 1). Next, the joint toxicity of multi-component mixtures would be predicted using disparate models when mastering some knowledge of EC50 of an individual component.

Table 1

NameFunction regression expression
Probit
Logit
Morgan-mericier
Weibull/GomertzP(Con) = 1−exp(−exp(θ1 + θ2lg(Con)))
Generalized logit 1
Generalisied logit 2
Aranda-ordaz
Logit with box-cox-transformation
Weibull with box-cox- transformation
Probit with box-cox- transformation

Classical regression models for calculating the concentration-response relationships.

Con denotes the concentration of the tested substance, P (Con) is the mean effect of the model and θ1, θ2, θ3 are estimated model parameters. Where Φ (y) is the cumulative normal (Gauss) distribution, meaning that the probability that a standard normal random variable is less than y.

Graphical Methods

Isobologram

Isobologram was originally developed by Loewe and Muischnek (1926). This method can be used to determine whether the multi-component mixtures behaved in similar modes of action. The straight-line presents the theoretical additive effect of the two components connect two individual effects caused by the fixed concentration. The accuracy of the estimated intercepts of the theoretical isobole with the axis affects the degree of deviation of observed data and experimental data. In many cases, the deviation may result in inadequacy predication. Bernbaumhas put forward interaction index (CI) to directly present the effects on non-interactive multi-compounds mixtures from the dose-effect curve of the individual compound (Berenbaum, 1981). The graphical representation is shown in Figure 3A. In the binary mixture, the equation of the CI is showed in Equation (1). In this equation, a1 and b1 are the dose of components in the mixture, A and B are the dosed of the individual agents giving rise to the same effects as the mixtures. The combination index (CI) determinates the joint effect of the combinations, that shows in Figure 6.

Figure 3

Applied with the 95% confidence interval to simplify the process of the judgment, Ding proposed an improvement approach to characterize the binary toxicities of mixtures and applied this judgment to evaluate the joint toxicity of antibiotics successfully (Figure 3B) (Ding et al., 2015). The EC50 of compounds A and B are presented by a2 and b2, respectively. (a1, a3) is the 95% confidence limit of a2, while (b1, b3) is the 95% confidence limit of b2. The values a0 and b0 represent the equivalent concentration of EC50 for mixture A and B. According to the location of P (a0, b0), the joint toxicity effect of mixture A and B can be judged. If P is located between line a1b1 and a3b3, the joint effect is additive, if P is underline a1b1, the joint toxicity is synergy if P is above line a3b3 and in the triangle a3b3Q, the joint toxicity is independent and if P is located outside the rectangle a3Ob3Q, the joint effect is antagonistic. Much of the available literature on isobologram deals with the question of mixtures toxicity of multi-component. Using the experimental data of poly and perfluoroalkyl substances (PFAS) based on an amphibian fibroblast cell line, Hoover et al. employed the isobolograms to ascertain the interaction type of the binary mixtures of PFAS, the additive effects were shown in most combination tests (Hoover et al., 2019). Mori et al. in his research conducted isobologram analysis which based on toxic unite (TU) to explore the synergism of the mixture of Tetramethylammonium hydroxide (TMAH) and potassium iodide (KI) to Daphnia magna, using four mixture ratios of TMAH and KI (4:1, 3:2, 1:1, 1:4), as shown in Figure 3C (Mori et al., 2015). Similarly, the bioluminescent inhibition of V. fisheri by exposure to binary mixtures of fullerene soot and nonylphenol (Figure 3E), triclosan (Figure 3F), malathion (Figure 3G), and glyphosate (Figure 3H) using isobologram were discussed, indicating antagonism (Sanchis et al., 2016). Moreover, if three substances are used together, the equivalent surface is shown in Figure 3D. The method of judging the toxicity effect of the three substances is the same as that of a binary substance (Huang et al., 2019b). One should be cautious, isobologram is a graphical method so that variability of date can-not be reflected in this approach and visualization of the effect of mixtures is different, and the graphics draw require large data sets to produce sufficiently reliable results (Berenbaum, 1981; Cassee et al., 1998; Meadows et al., 2002). It is worth mentioning isobologram is unconcerned about chemical structures or modes of action of the individual chemical.

Response Surface

Response or effect surface methodology is based on mathematical relationships, which is obtained by multiple linear regression. The zero interaction responses surface for dose additivity method is depicted in the Equation (2).

where dA, dB represents the dose of components A and B in the mixture. E is the zero interaction combinations effect. The variable α and μ mean the EC50 and slope, respectively.

In the case of significant, the response surface for binary components mixtures can be displayed in three dimensions. The concentration of the two substances is expressed as X-axis and Y-axis, respectively. The combined effect of the two substances is showed as Z-axis. The response surface is displayed by different mathematical model fitting. Neither the synergistic effect nor antagonistic effect of the response surface is flat, because of the two components non-existence of interaction. Ren et al. analyzed the binary effect of copper and lead toward luminescent bacteria based on the response surface, as shown in Figure 4A. They found the transition between different effects with the concentration of metal in the mixtures (Ren et al., 2004). Similarity, Tsai et ai. used an isobologram analysis to test combinations of two drugs etoposide and cisplatin against NCI-H226 cell lines, the results obtained from the preliminary analysis are shown that there was no synergism between etoposide and cisplatin at the 50% level, the sketch modified from Groten et al. (Figure 4B) (Tsai et al., 1989; Groten et al., 2001). It should be noted if predictions for doses extend intended range may result in deviations. Moreover, the approach is related to the relatively large number of parameters that must be determined by non-linear regression techniques, which requires a lot of experiments.

Figure 4

Methods Based on Toxic Unite

Sprague and Ramsay first proposed Toxic Unit (TU) concept to describe the contribution of the individual component to the overall effect (Sprague and Ramsay, 1965). In the luminescent bacteria toxicity assay, the relative luminescence intensity value (RLU) decreased with the toxic substance mass concentration increased. While TU is derived from scaling a measured compound concentration to its inherent EC50 (Fulladosa et al., 2004). The TU concept has been applied to weigh the effect of single chemicals for the assessment of joint toxicity, can be calculated by Equations (3–5):

where ci represents the concentration of component i. When the mixture is at its EC50, TUi is the toxic unit of i, M is the sum of the toxic units, (TUi)max is the maximum toxic unit of the mixture, and M0 is the ratio of M and (TUi)max.

In such a way, M quantified the joint effect that showed in Figure 5 (Chen and Huang, 1996). This approach is relatively simple and can roughly reflect the relationship between the reaction mixtures, the central premise of this concept is that each of these components alone cannot produce any effect.

Figure 5

Based on the concept of TU, some other joint effect indexes have been proposed to perfect the evaluation system (Table 2). The evaluation expressions are concise and the model expectations are replaced by a value to quantify the degree of deviation. However, there are still some inevitable limitations: (1) Only EC50 value is provided in the indexes, the evaluation about different levels are unavoidable lost. (2) Similar to the isobologram method, experimental variability also consists of these joint effect indexes method, it is unlikely to meet calculated responses precisely (Altenburger et al., 2003). In some cases, the combination index associated with the isobologram was put forward to address the issue of combination ecotoxicological assessments, effectively (Chou, 2006; Wang et al., 2015d; Godoy et al., 2019).

Table 2

Type of index and type of interactionJoint effect index
Additivity index (AI)Similarity parameter (λ)Mixture toxicity index (MTI)
Mathematical modelThe expression depicted as Equation (4)
if M <1, AI = 1/M−1,
if M > 1, AI = 1/M + 1,
if M = 1, AI = M−1
SYNAI > 0λ > 1MTI > 1
Simple ADDAI = 0λ = 1MTI = 1
INDMTI = 0
ANTAI <00 < λ <1MTI <0
Partly ADD0 < MTI <1
ReferenceL, 1977Prakash et al., 1996Konemann, 1981

Criteria for the judgment of joint effects of different evaluation indexes.

the effect is inexistence when used this type.

Apart from the fore-mentioned joint effect indexes, other mathematical models can predict and analyze the joint effect of the binary chemical mixtures. The models always analyze the difference between the calculated value from the model and the observed value from the actual experiments. In this review, two popular evaluation approaches are recommended in detail.

The first approach (referred to statistical approach 1) was proposed by Ribo and Rogers (1990), which assumed that two in a mixture have the same TMOA do not react with each other. The interactive toxicity in that binary mixtures was studied by calculating the difference (EC50,Diff) between the observed toxicities (EC50,Obs) of the binary mixtures with the calculated toxicities (EC50,Calc). The computed mixture toxicities were calculated as Equations (6) and (7).

where CM is the sum of the concentrations of the substances in the binary mixture (Cx and Cy) and the EC50,x, and EC50,y is the effective concentrations representing the observed toxicity of each substance separately. The obtained EC50,Diff value implies the type of interaction between the two substances in the mixture (Figure 5). ADD effect is expressed when the difference is considered statistically non-significant (at a 95% confidence level), the definition of ADD effect is also applied to the following statistical approach 2. With this approach, Fulladosa et al. predicted the toxicity of the possible binary equitoxic mixtures of Co, Cd, Cu, Zn, and Pb to A. fischeri. The research showed the most combined effect of the metals was antagonistic for Co-Cd, Cd-Zn, Cd-Pb and Cu-Pb, while the synergistic effect for Co-Cu and Zn-Pb and merely additive in some cases (Fulladosa et al., 2005).

The second popular approach (referred to statistical approach 2) postulated that there may be interactions between two chemicals that act independently on biological models through different mechanisms. For that reason, the computed mixture toxicities can be expressed as new toxicity units (TU), defined as TU= 100/EC50, the definition is different from the above classical TU. The results can be calculated by Equations (8) and (9).

where TUX and TUY are the toxicity units of X and Y, TUHo (X+Y) is the toxicity units representing the a priori expected toxicity of a mixture of two components X and Y. The interaction effect was indicated by TUDiff value (Figure 5). Fulladosa et al. studied binary As (V) species mixtures based on the statistical analysis of TUObs and TUHo, they noted that As (V) species mixture interaction types were ANT at the pH of 6 and 7 (Fulladosa et al., 2004).

Based on the derivative models associated with TU, several studies have been conducted for ecotoxicity assay of diverse pollutants toward luminescent bacteria, which are categorically reviewed as shown in Table 3.

Table 3

CategoriesTest substancesTest organismMixture ratiopHIndicatorsEC50 (Experimental)SignificanceJoint effectReferences
Heavy metals mixturesCd (II) + Cu (II)P. phosphoreumEC50: EC505.5–6MTI/TU8.581 mg/L0.969ANT/ANTZeb et al., 2016
Cd (II) + Pb (II)P. phosphoreumEC50: EC505.5–6MTI/TU2.076 mg/L0.9036ANT/ANTZeb et al., 2016
Cu (II) + Pb (II)P. phosphoreumEC50: EC505.5–6MTI/TU1.781 mg/L0.9266ADE*/ADE*Zeb et al., 2016
Cd (II) + Cu (II) + Pb (II)P. phosphoreumEC50: EC505.5–6MTI/TU2.443 mg/L0.9548ADE*/ADE*Zeb et al., 2016
Pb (II) + Cd (II)V. qinghaiensisEC50: EC50TU1.2 TU0.9719ANTXiong et al., 2007
Pb (II) + Cr (VI)V. qinghaiensisEC50: EC50TU2.47 TU0.93ANTXiong et al., 2007
Cd (II) + Cr (VI)V. qinghaiensisEC50: EC50TU0.67 TU0.9554SYNXiong et al., 2007
Co (II) + Cd (II)A. fischerEC50: EC506A1/A20.755 mmol/L0.001/0.031ANT/ANTFulladosa et al., 2005
Co (II) + Zn (II)A. fischerEC50: EC506A1/A20.666 mmol/L0.098/0.069ADD/ADDFulladosa et al., 2005
Co (II) + Cu (II)A. fischerEC50: EC506A1/A20.258 mmol/L1.2E−6/6.4E−6SYN/SYNFulladosa et al., 2005
Co (II) + Pb (II)A. fischerEC50: EC506A1/A20.359 mmol/L0.284/0.12ADD/ADDFulladosa et al., 2005
Cd (II) + Zn (II)A. fischerEC50: EC506A1/A20.126 mmol/L0.002/0.004ANT/ANTFulladosa et al., 2005
Cd (II) + Cu (II)A. fischerEC50: EC506A1/A20.078 mmol/L0.181/0.113ADD/ADDFulladosa et al., 2005
Cd (II) + Pb (II)A. fischerEC50: EC506A1/A20.152 mmol/L0.001/2.8E−6ANT/ANTFulladosa et al., 2005
Zn (II) + Cu (II)A. fischerEC50: EC506A1/A20.012 mmol/L0.735/0.440ADD/ADDFulladosa et al., 2005
Zn (II) + Pb (II)A. fischerEC50: EC506A1/A20.006 mmol/L0.023/0.001SYN/SYNFulladosa et al., 2005
Cu (II) + Pb (II)A. fischerEC50: EC506A1/A20.007 mmol/L0.002/3.1E−6ANT/ANTFulladosa et al., 2005
Cu (II) + Zn (II)E. coli HB101pUCD6070.48: 1 (mg/L: mg/L)7TU/AI/MTI8.54 mg/LANT/ANT/ANTZhou et al., 2018
Cu (II) + Zn (II)E. coli HB101pUCD6071.44: 1 (mg/L: mg/L)7TU/AI/MTI8.14 mg/LADD*/ANT/ADD*Zhou et al., 2018
Cu (II) + Zn (II)E. coli HB101pUCD6074.33: 1 (mg/L: mg/L)7TU/AI/MTI7.89 mg/LANT/ANT/ANTZhou et al., 2018
Zn (II) + Hg (II)E. coli HB101pUCD60712.73: 1 (mg/L: mg/L)7TU/AI/MTI9.51 mg/LANT/ANT/ANTZhou et al., 2018
Zn (II) + Hg (II)E. coli HB101pUCD6074.24: 1 (mg/L: mg/L)7TU/AI/MTI7.42 mg/LADD*/ANT/ADD*Zhou et al., 2018
Zn (II) + Hg (II)E. coli HB101pUCD6071.41: 1 (mg/L: mg/L)7TU/AI/MTI5.05 mg/LANT/ANT/ANTZhou et al., 2018
Hg (II) + Cu (II)E. coli HB101pUCD6070.49: 1 (mg/L: mg/L)7TU/AI/MTI4.88 mg/LADD*/ANT/ADD*Zhou et al., 2018
Hg (II) + Cu (II)E. coli HB101pUCD6070.16: 1 (mg/L: mg/L)7TU/AI/MTI5.53 mg/LADD*/ANT/ADD*Zhou et al., 2018
Hg (II) + Cu (II)E. coli HB101pUCD6070.05: 1 (mg/L: mg/L)7TU/AI/MTI0.72 mg/LANT/ANT/ANTZhou et al., 2018
Dibutyl phthalate and antibioticsChlorotetracycline hydrochloride + Dibutyl phthalateA. fischerEC50: EC50TU49.82 mg/L0.9976SYNWei et al., 2018
Sulfamethazine + Dibutyl phthalateA. fischerEC50: EC50TU46.51 mg/L0.9625SYNWei et al., 2018
Sulfamethazine + Dibutyl phthalateA. fischerEC50: EC50TU55.82 mg/L0.9989SYNWei et al., 2018
Sulfadiazine + Dibutyl phthalateA. fischerEC50: EC50TU26.309 mg/L0.9642SYNWei et al., 2018
Antibiotics mixturesChlortetracycline hydrochloride + SineptinaP. phosphoreumEC50: EC507 ± 0.5TU0–2 TU0.9610SYNRen et al., 2011
Chlortetracycline hydrochloride + SalinomycinP. phosphoreumEC50: EC507 ± 0.5TU0–2 TU0.9613ANTRen et al., 2011
Chlortetracycline hydrochloride + FlavomycinP. phosphoreumEC50: EC507 ± 0.5TU0–2 TU0.9554ANTRen et al., 2011
Sineptina + SalinomycinP. phosphoreumEC50: EC507 ± 0.5TU0–2 TU0.9701SYNRen et al., 2011
Sineptina + FlavomycinP. phosphoreumEC50: EC507 ± 0.5TU0-2 TU0.9794SYNRen et al., 2011
Salinomycin + FlavomycinP. phosphoreumEC50: EC507 ± 0.5TU0–2 TU0.9872ANTRen et al., 2011
Herbicidal mixturesFenuron + SimetrynP. phosphoreumEC50: EC50TU1.18 TUADDGe et al., 2014
Fenuron + AtrazineP. phosphoreumEC50: EC50TU1.44 TUANTGe et al., 2014
Fenuron + PrometonP. phosphoreumEC50: EC50TU2.28 TUANTGe et al., 2014
Fenuron + PrometrynP. phosphoreumEC50: EC50TU2.52 TUANTGe et al., 2014
Monuron + SimetrynP. phosphoreumEC50: EC50TU0.71 TUSYNGe et al., 2014
Monuron + AtrazineP. phosphoreumEC50: EC50TU2.84 TUANTGe et al., 2014
Monuron + PrometonP. phosphoreumEC50: EC50TU1.09 TUADDGe et al., 2014
Monuron + PrometrynP. phosphoreumEC50: EC50TU1.82 TUANTGe et al., 2014
Monolinuron + SimetrynP. phosphoreumEC50: EC50TU0.90 TUSYNGe et al., 2014
Monolinuron + AtrazineP. phosphoreumEC50: EC50TU1.11 TUADDGe et al., 2014
Monolinuron + PrometonP. phosphoreumEC50: EC50TU2.33 TUANTGe et al., 2014
Monolinuron + PrometrynP. phosphoreumEC50: EC50TU1.44 TUANTGe et al., 2014
Diuron + SimetrynP. phosphoreumEC50: EC50TU0.68 TUSYNGe et al., 2014
Diuron + AtrazineP. phosphoreumEC50: EC50TU1.22 TUANTGe et al., 2014
Diuron + PrometonP. phosphoreumEC50: EC50TU2.59 TUANTGe et al., 2014
Diuron + PrometrynP. phosphoreumEC50: EC50TU1.40 TUANTGe et al., 2014
Pesticides mixturesAcifluorfen + Dimethoate + ProchlorazA. fischeriEC50: EC50: EC50TU0.66 TUSYNCedergreen et al., 2012
Acifluorfen + Diquat + ProchlorazA. fischeriEC50: EC50: EC50TU0.51 TUSYNCedergreen et al., 2012
Dimethoate + Prochloraz + ChlorfenvinphosA. fischeriEC50: EC50: EC50TU0.81 TUSYNCedergreen et al., 2012
Heavy metal and other compoundsCu + nitroaromatic compoundsP. phosphoreum0.2 (Cu):1 (EC50: EC50) 0.5 (Cu):1 (EC50: EC50) 0.8 (Cu):1 (EC50: EC50)TU/AIAt low Cu concentration, most binary mixtures show simple addition. Three binary mixtures are antagonistic. At medium Cu concentration, all effect types are existed. At high Cu concentration, the binary combined mixtures most showed antagonism, simple addition is only found for Cu combined with o-dinitrobenzene, p-nitrobromobenzene, p-nitrobenzoic acid and 2,4-dinitropheolSu et al., 2012
Pesticides and MTBEMTBE + Dichlofluanid/TBT/Linuron Diuron /Sea nine/IrgarolA. fischeri0.1 mg/L (MTBE): EC50 of individual pesticideMTIThe combinations of MTBE and pure Diuron, Dichlofluanid, TBT, and Linuron are showed synergy. The addition of Sea nine and Irgarol did not increase the toxicity of the MTBEHernando et al., 2003

Characterization of the acute toxicity of multi-component mixtures assessments associated with TU by luminescent bacteria.

*

means the partly effect.

– means this item not mentioned in the article.

A1 method is the popular evaluation method which is mentioned as Equations (6) and (7), and A2 methods is showed as Equations (8) and (9).

All of the methods mentioned above are not considered the toxic mode of actions. The isobologram and response surface method are always addressed the binary toxicity of two-components mixtures. The toxic unit and corresponding methods except the statistical approach 2 (referred to TUDiff), assumed that no interaction among the components in the mixture. It is decisive to take advantages and disadvantages into consideration and choose the appropriate method for joint toxicity evaluation.

Methods of Concentration Addition and Independence Action

Concentration Addition Method

As reported by the EPA (EPA, 2000), the additivity is defined as “when the effect of the combination is estimated by the sum of the exposure levels or the effect of the individual chemicals,” In the context of the above definitions, the definitions about “effect” and “sum” draw forth two classical mathematical models. One is Concentration addition (CA, also termed Dose Addition or Loewe additivity), which assumes the constituents in a mixture has the same TMOA, differing in their efficacy. The joint effect of compound mixtures that have the same TMOA can be calculated using the CA model based on the concentration-response relationship of single substances (Loewe and Muischnek, 1926). The mathematical expression of CA by Equation (10):

where ECx,mix is the effective concentration of the mixture eliciting x% effect, ECx,i denotes the concentration of the ith component when exists individually and elicits the same effect (x%) as the mixture, pi is the molar concentration ratio of the ith component in the mixture.

CA implies that two (or more) components act on acceptors with common TMOA indiscriminately. That is to say, components can be substituted for each other in the mixtures (Figure 6A). The concept of CA hypothesized the combined effect may develop in the presence of multi-components. This hypothesis is significant because different categories of organisms in the ecosystem are exposed to a mixture of pollutants.

Figure 6

Independent Action Method

In another case, the resistance of the organism to one component can not be replaced by another component in the mixture. That is to say, the response of two substances is statistically independent, the action mode of which belongs to different TMOAs (Figure 6B). The model IA model (or response addition) is another approach was hereby proposed to assess the mixture toxicity of components with dissimilar TMOAs, which means multi-component interaction with different target sites (Fraser, 1872; Bliss, 1939). IA model is commonly defined as:

where E(cmix) is the overall effect caused by the total effect of the mixture, cmix is the total concentration of a mixture; E(ci) is the value that denotes the single effect of the ith component corresponds to the concentration of ciin the mixture. The IA model is apt to predict the joint effect of mixtures with definite compositions.

Which of the Two Models, CA or IA, Predict the Date More Precise?

From IA expression, we can conclude that toxicity assessment based on IA relies on the knowledge of the individual effect, which required from the concentration-response curves of all individual component. Moreover, E(ci) is lower than E(cmix) commonly (Junghans et al., 2006; Backhaus and Faust, 2012). That is, the limitation of IA may provoke the accuracy deviation of some predictions. Considering the simplicity ease to utilization, TU based on the CA model has been used in the context of unknown toxic effects of a single component.

There are partly researches were implemented to investigate the feasibility of CA and IA (Zhang et al., 2008; Liu et al., 2009; Zhou et al., 2010; Mo et al., 2014; Villa et al., 2014; Tong et al., 2015). It is not to say, CA (or IA) is only propitious to the similar (dissimilar) toxicology modes of action. The accuracy of CA even for the prediction of the toxicity of multi-component mixture composed with different TMOAs is investigated. Design uniform equivalent-effect concentration experimental method, the combined toxicity prediction of six organophosphorus pesticides (Fenitrothion, Malathion, Dicapthon, Chlormephos, Methyl parathion and Famphur) based on V. qinghaiensis sp.-Q67 can be precisely assessed by CA (Zhang et al., 2008). Others (Junghans et al., 2006; Zhou et al., 2010) also validated that CA may be a reasonable assumption for the toxicity of the multiple-component mixtures of pesticides regardless of the similarity or dissimilarity of their TMOA. Simultaneously, CA and IA are feasible for predicting the toxicological effect of some specific mixtures. Villa et al. tested the toxicity of complex mixtures of chemicals including narcotics, polar narcotics, herbicides, insecticides and fungicides toward A. fisheri. In this experiment, the set of chemicals was attributed to different groups on account of dissimilar toxicology modes of action. The finding demonstrated that CA and IA models are both feasible for predicting the toxicological effect of mixtures, and the difference between the two models was never higher than a factor of four (Villa et al., 2014). However, some researchers believe that the applicability of the two models remains a matter of debate in terms of accuracy in predicting the combined effects of mixtures. The potential risk management options should be studied more carefully (Carbajo et al., 2015; Godoy and Kummrow, 2017; Yang et al., 2017).

Integrative Model of CA and IA

Some toxicants are composed of the components with similar or completely different toxic modes of action in the ecological environment. That is to say, some of which may have the same mode of action, while others may differ. Therefore, both the CA and the IA model may be not suitable for evaluating the combined toxicity of these mixtures. To work this out, the integrative models were developed to predict the toxicity of those complex mixtures (Ra et al., 2006; Qin et al., 2011; Kim et al., 2014; Mo et al., 2017). Among them, the TSP (two-stage prediction) model is one of the most commonly used approaches based on the integrative model concept (Mo et al., 2017). The basic principle of the TSP model is to apply the CA model and IA model that divided in stages to predict the joint toxicity of mixtures as following: In the first stage, the components in the mixture with similar TMOA are assigned into the same group, and joint toxicity of each group was calculated using CA model (Equation 10). In the second stage, the IA model is used to predict the mixture toxicity of all the groups with different TMOAs (Equation 11). In this way, the prediction of the combined toxicity of the mixture with similar and different TMOAs can be predicted.

Mo et al. selected six phenolic compounds (Methyl phenol, 2-Nitrophenol, 4-Nitrophenol, 2.4-Dichlorophenol, Phenol and 2-chlorophenol) and heavy metals (Cr, Cu, Ni, Cd, Ag, and Hg) with different TMOAs to compose the mixtures with V. qinghaiensis sp.-Q67, a series of mixtures were designed with equivalent-effect concentration and fixed concentration ration. By comparing the prediction error of CA, IA and TSP models, the TSP model exhibited better performance on predicting the overall effect of a mixture containing compounds with different TMOAs than CA and IA models (Mo et al., 2017). The predicted combined toxicity of mixtures of imidazolium and pyridinium ionic liquids in the ratios of their EC50, EC10, and NOEC on luciferase also showed that the TSP method preceded the CA and IA methods effectively (Ge et al., 2014).

The above methods are both attributed to traditional experimental methods, with the aid of the individual concentration-effect relationship so as the final judgment of the multi-component mixture can be confirmed. Moreover, experimental toxicology studies usually used experimental animals, tissues, bacteria, cell. It is no doubt these requirements have improved the complexity of the operations. Hence, researchers are devoted to developing alternative methods to conquer these limitations of traditional experimental methods. QSAR (Roy et al., 2015a,b), Read-across (RA) (Jeliazkova et al., 2010), Molecular docking (Yao et al., 2013) and Expert systems (Roy and Kar, 2016), are deemed faster, cheaper and can manifest more information (Mwense et al., 2004). In these developing alternative methods, the QSAR model is one of the most widely used methods in silico.

QSAR Assisted Toxicity Prediction

QSAR, as a calculated method has been widely applied in toxicology. QSAR is an acronym for Quantitative Structure-Activity Relationship. The concept was delivered in 1872, and the modern age of QSAR analysis originated from the works of Hansch et al. (1962). By mathematical function, QSAR is a statistical approach which express a relationship between the magnitude of biological effect (BA) and changes in a molecular structure, as Equation (12) showed:

where i denotes a specific chemical of a series, this series may be of homogeneous, or of heterogeneous substances.

In the study of constructing a QSAR model, it is decisive to acquire and screen molecular structure descriptors (also called indicators). At present, more than 3,000 molecular descriptors have been defined and applied to QSAR models to predict group biotoxicity. Choose distinct indicators when adopting the QSAR model that will affect the validity of the predicted results, therefore the preparation of the indicators is the key to the construction (Khan et al., 2020). Among abundant kinds of indicators, the partition coefficient of a chemical between n-octanol and water is supposed to be the most effective indicator (Tichy et al., 2008). However, it is difficult to use Kow to determine the combined toxicity of mixtures, since the available date was obtained by UV spectrophotometry or HPLC is only suitable for determining the Kow of single chemicals. Furthermore, Verhaar extended the C18-EmporeTM disks/water partition coefficient (KMD) to predict the bioconcentration of mixtures, which was found to have a close relationship with logKow (Verhaar et al., 1995). The C18-EmporeTM disks/water partition coefficient (KMD) can be calculated from Equation (13), the unmeasurable problem was successfully solved.

where W is the volume of solution, V is the volume of hydrophobic phase, is the initial amount of chemical i in water, n is the total number of individual chemicals in the mixture, and KSDi is the partition coefficient of individual chemical i. The value of W/V was suggested to be 6.8 × 105.

Not limited to the partition coefficient of chemical between n-octanol and water, studies concerned other kind indicators to employ QSAR to assess ecotoxicity evaluations of diverse pollutants toward familiar luminescent bacteria, A. fisheri, P. phosphoreum and V. qinghaiensis sp.-Q67, which are categorically reviewed as shown in Table 4.

Table 4

IndicatorsComponentsModelsTest organismReferences
logKowNarcoticslg1/Cnarc(mmol/L) = 0.94logKow − 2.61A. fischeriVighi et al., 2009
n = 23, R2 = 0.92
Polar narcotics compoundslg1/Cpolar(mmol/L) = 0.502logKow + 0.294A. fischeriVighi et al., 2009
n = 25, R2 = 0.82
KMDHalogenated benzeneslgEC50 = 0.928logKMD + 0.224P. phosphoreumLin et al., 2002
n = 74, R2 = 0.973, SE = 0.113, F = 287.785, P < 0.001
EbindSulfonamides (SA) –SA potentiators (SAP)[[Mathtype-mtef1-eqn-40.mtf]]A. fischeriWang et al., 2018
n = 18, R2 = 0.699, RMSE = 0.161, F = 17.432, P = 0.000, Q2loo = 0.600, RMSEloo = 0.174, , RMSEP = 0.067
EbindSulfonamides (SA) –tetracyclines (TC)A. fischeriWang et al., 2018
n = 15, R2 = 0.813, RMSE = 0.068, F = 26.138, P = 0.000, Q2loo = 0.714, RMSEloo = 0.075, = 0.601, RMSEP = 0.039
EbindSA potentiators (SAP) –tetracyclines (TC)A. fischeriWang et al., 2018
n = 8, R2 = 0.829, RMSE = 0.092, F = 12.109, P = 0.012, Q2loo = 0.603, RMSEloo = 0.111,
= 0.760, RMSEP = 0.069
EHOMO, ELUMO2,4 Dinitrotoluene –aromatic compounds−lg(1/EC50) = 3.672−0.882ELUMOA. fischeriYuan et al., 2002
n = 8, r = 0.925, s = 0.406
Ebind, pKaSulfonamide antibiotics mixturesP. phosphoreumZou et al., 2012
n = 7, R2 = 0.954, SE = 0.028, F = 20.659, P = 0.017
EHOMO, ELUMO(0.2 × EC50) Cd –chlorinated anilineslgEC50 = −28.212(EHOMOELUMO)−4.187P. phosphoreumJin et al., 2014
n = 9, R2 = 0.832, SE = 0.0603, F = 34.769, P = 0.001, , RMSCV = 0.064
EHOMO, ELUMO(0.5 × EC50) Cd –chlorinated anilineslgEC50 = −42.770(EHOMOELUMO)−7.243P. phosphoreumJin et al., 2014
n = 9, R2 = 0.837, SE = 0.0900, F = 35.87, P = 0.001, , RMSCV = 0.094
EHOMO, ELUMO(0.8 × EC50) Cd –chlorinated anilineslgEC50 = −55.638(EHOMOELUMO)−9.791P. phosphoreumJin et al., 2014
n = 9, R2 = 0.903, SE = 0.0869, F = 65.18, P = 0.000, = 0.8204, RMSCV = 0.107
CSEV, S(0.2 × EC50) Cu –nitroaromatic compoundslg1/EC50 = −1.459+0.043CSEV+0.684SP. phosphoreumSu et al., 2012
n = 11, R2 = 0.828, SE = 0.251, F = 19.2, P = 0.001
CAA(0.5 × EC50) Cu –nitroaromatic compoundslg1/EC50 = −8.039+0.043CAAP. phosphoreumSu et al., 2012
n = 11, R2 = 0.727, SE = 0.411, F = 25.0, P = 0.001
CAA(0.8 × EC50) Cu –nitroaromatic compoundslg1/EC50 = −8.754+0.046CAAP. phosphoreumSu et al., 2012
n = 11, R2 = 0.732, SE = 0.432, F = 24.6, P = 0.001
VE2_B(p) TIC3(0.2 × EC50) Zn–nitro-substituted benzeneslg1/EC50 = 20.540−53.948VE2_B(p)−0.033TIC3P. phosphoreumZhang et al., 2019a
n = 11, R2 = 0.933, SE = 0.163, F = 55.554, P < 0.001
Eig06_AEA(dm)(0.5 × EC50) Zn–nitro-substituted benzeneslg1/EC50 = 3.760+1.100Eig06_AEA(dm)P. phosphoreumZhang et al., 2019a
n = 11, R2 = 0.856, SE = 0.232, F = 53.318, P < 0.001
Eig06_AEA(dm)(0.8 × EC50) Zn –nitro-substituted benzeneslg1/EC50 = 3.908+1.511Eig06_AEA(dm)P. phosphoreum(Zhang et al., 2019a
n = 11, R2 = 0.937, SE = 0.201, F = 13.351, P < 0.001
Benzene—its derivativesA. fischeriChang et al., 2016
n = 40, R2 = 0.827, , SD = 0.269, , F = 110.935, P < 0.0001.
σp, C*Aldehyde –cyanogenic toxicantP. phosphoreumLin et al., 2005
R2 = 0.868, SE = 0.232, F = 121.769, P = 0.000, M is the sum of toxic units
Oaldehyd, CcyanogenicAldehyde –cyanogenic toxicantTU = 1.00 ± 0.20 when Oaldehyde- Ccyanogenic >-0.125P. phosphoreumTian et al., 2012a
TU = −27.6 × Oaldehyde − 5.22 × Ccyanogenic − 6.97 when Oaldehyde- Ccyanogenic ≤ -0.125
n = 40, r = 0.887, SE = 0.195, F = 140, p < 0.001, = 0.748
log Kmbw, log KmcwNonpolar/polar—narcotic—chemical mixtures, P. phosphoreumLin et al., 2003
n = 84, R2 = 0.948, SE = 0.166, F = 745.201, p = 0.000
RDF035m HATSs H-047Antibiotic and pesticidespEC50,mix = (2.2780 ± 0.2173)−(0.0546 ± 0.0050) × (RDF035m)mix
−(0.0158 ± 0.0024) × (HATSs)mix − (0.3124 ± 0.0214) × (H−047)mix
A. fischeriQin et al., 2018
n = 31, R2 = 0.9366, RMSE = 0.1345, F = 132.89, = 0.9087, = 0.9087
Mor13s, L/Bw, Eig08_EA(ed) Eig09_EA(ed)Aromatic brominated –chlorinated disinfection byproductspEC50,mix = 3.7370(±0.0324)−0.1710(±0.0135)Mor13s−0.0421(±0.0135)L/Bw
+0.0866(±0.0093)Eig08_EA(ed)−0.0085(±0.0017)Eig09_EA(ed)
V. qinghaiensis sp.-Q67Chen et al., 2019
n = 31, = 0.7804, = 0.7782, = 0.8057, CCC = 0.8778, = 0.7700
R7_TPSA, R2_pKa,Natural deep eutectic solventslogEC50 = −5.596R7_TPSA+0.145R2_pKaA. fischeriGiner et al., 2020
R2 = 0.858, SE = 0.166, F = 103, P < 0.0001

Characterizations of various types of QSAR models for prediction of joint toxicity.

Researches about QSAR are still of significant interest in the development of innovative models in environmental toxicity prediction. The emergence of the QSAR models fills the gap in predicting the combined toxicity of organic compounds and heavy metals since the information in this field is still scarce. Jin et al. developed three QSARs to determine the individual EC50 of Cd and nine chlorinated anilines (o-chloroaniline, m-chloroaniline, p-chloroaniline, 2,3-dichloroaniline, 2,4-dichloroaniline, 2,5-dichloroaniline, 2,6-dichloroaniline, 3,4-dichloroaniline, and 2,4,5-trichloroaniline) with P. phosphoreum, setting three different levels of Cd concentrations (low, medium and high levels) to mix with chlorinated anilines and the results showed that the number of chlorinated anilines manifesting synergy with Cd is decreasing as the concentration of Cd increases. The robustness of the models was confirmed by comparing the experimental and predicted values, and all the relative error values remain within 16% (Jin et al., 2014). Likewise, the joint toxicities of Cu (low, medium and high levels) with 11 nitroaromatic compounds (nitrobenzene, o-dinitrobenzene, m-nitrobromobenzene, p-nitrobromobenzene, o-nitroaniline, p-nitroaniline, p-nitrobenzoic acid, o-nitrophenol, m-nitrophenol, p-nitrophenol and 2,4-dinitrophenol) were studied by P. phosphoreum with developed QSAR analysis, there is a good agreement between the predicted values and experimental with R2 = 0.764, P = 0.000, (Su et al., 2012). By drawing on the study of Su's, Zhang's, Su, Zhang, Li, Qin and Zhang (2019a) work on the acute toxic effect of mixtures between metal Zn and above-mentioned 11 nitroaromatic compounds established robust QSAR models to predict the joint interaction when combined with Zn at low, medium, and high concentrations (Zhang et al., 2019a). Not limited to heavy metals and organics, the joint toxicity of mixed organics can also be predicted by QSAR models. Qin et al. developed a generalized QSAR model for predicting the additive and non-additive toxicities of multi-component mixtures, the experiment tested the joint toxicity of 45 multi-component mixtures composed of two antibiotics (Tetracycline hydrochloric and Chloramphenicol) and four pesticides (Metribuzine, Trichlorfon, dichlorvos, Linuron) to A. fisheri. Compared with classical CA and IA models, the result demonstrated that the QSAR model exhibited high predictive capability for predicting joint toxicity (Qin et al., 2018).

The Potential Mechanisms of Joint Toxicity

Though more and more efforts have been devoted to investigating mechanisms with mixture toxicity, joint toxicity mechanisms, especially for specific pollutants, are still too complex to figure out. In recent years, there has been growing evidence regarding the experimental factors including mixture ratio setting and complexity of toxic compounds play a vital role in shaping the joint effect, which in turn impacts the determination of toxicological mechanism. At the same time, considering the environmental factors on the joint effect study, the ecological environmental media should be researched more adequately. Because of the above-mentioned reasons, the studies concentrated on the mechanisms of joint toxicity toward luminescent bacteria are merely few. Analogous to other test organisms, impair mechanisms of the luminescent bacteria are speculated in this section.

Trojan-Horse Effect and Reverse Trojan-Horse Effect

The cellular membrane acts as a significant site for interactions of chemicals. A wealth of reports showed that physiological activities in the inner of an organism are related to the cellular membranes. The functions of the cell membrane are considered in many aspects, such as transmembrane transport of small molecules, energy exchange, cell recognition, cell-mediated immunity, nerve conduction and metabolism regulation and so on (Oxender and Fox, 1978; Alberts et al., 1989). The impact of NPs on the membrane or wall integrity of algal cells (Angel et al., 2015; Sendra et al., 2017; Sousa et al., 2019) and bacteria cells (Kaweeteerawat et al., 2015; Martín-de-Lucía et al., 2017; Pulido-Reyes et al., 2017) (Figure 7A), has been evaluated extensively. Indications of damaged cell membrane make it interesting to study the potential biotoxicity effect of multi-component mixtures. A principal effect hypothesis concerning the cellular membrane damage which is described as the “Trojan-horse effect” has garnered considerable attention (Limbach et al., 2007). This hypothesis was first proposed and adopted gradually in nanomaterial. Owing to the large surface area, the ability of nanomaterials to act as carriers for a wide range of heavy metals and organic pollutants (Liu et al., 2018). CNTs and C60 are confirmed to carry heavy metals and dichlorodiphenyldichloroethylene (p, p'-DDE) entry into the inner cell via damaging cell membranes, respectively. Except for membrane disruptions, the entry process can impair other intracellular structures (receptors, ion channels, transporters, glycoproteins) (De La Torre-Roche et al., 2012; Wang et al., 2015a). The interaction toxicity of the Pb-nanotube also showed a synergistic effect, increasing over five times compared with the toxicity of single Pb (Martinez et al., 2013). Carbon adsorption of pollutants by nanomaterials was common in other organisms, like Japanese medaka (Su et al., 2013), Daphnia magna (Simon et al., 2015), Cyprinus carpio (Zhang et al., 2007). Based on the Trojan-horse effect, the schematic illustration of adsorption of metals ions on nanodiamonds (NDs) and the AgNP (Silver Nanoparticles) exposure, were shown in the Figures 7B,C (Quadros and Marr, 2012; Zhu et al., 2015). That mechanism also applied to luminescent bacteria, which elucidated that the synergistic effect of sulfonamides (SAs) with metal oxide nanoparticles (Meo-NPs) toward A. fischer (Figure 7D; Wang et al., 2016). As absorbent and carrier, metal oxide nanoparticles facilitated the entry of SAs into the inner of membrane completely and finally upon reaching the target site and release the encapsulated SAs molecules by making use of local physiological stimuli present here.

Figure 7

At the opposite extreme from the “Trojan-horse effect,” the term “Reverse Trojan-horse effect” is a relatively new name for referring to describing these phenomenons that antagonism is the predominant effect of multi-component mixtures (Barranger et al., 2019). Using A.fisheri, Sanchiss' study of the combined ecotoxicity of fullerene-soot and co-contaminants as malathion, diuron, triclosan, and nonylphenol indicated the joint effect of antagonism in all of the cases (Sanchis et al., 2016). In their analysis of this result, Sanchis et al. mentioned the special situation of aromatic rings as a tentative substance of allowing the production of micelles and changing the state of aggregation of fullerene aggregates. Therefore, the simple surface sorption process cannot clarify the antagonism between the fullerene-soot and co-contaminants. Interestingly antagonism appeared at low effect levels in wastewater–nanoparticle mixtures using recombinant bioluminescent cyanobacterium Anabaena sp. PCC 7120 strain CPB4337 (Martín-de-Lucía et al., 2017). Furthermore, the other four categories were added by Naasz et al. to discriminate different categories to recommend the joint toxicity upon mixture exposure, the additional mechanisms that are Surface enrichment, Retention, Inertism and Coalism (Naasz et al., 2018).

Competition for the Active Site

In the plasma membrane, some transmembrane proteins serve as receptors to detect and transduce chemical signals in the cell environment. According to the receptor theory, toxic targets are mainly receptors on the cell surface or in the nucleus and cytoplasm (Wu et al., 2016). The receptor theory represented that one pollutant with better absorption may replace those with poor adsorption when pollutants with similar chemical properties compete on the cell surface that will cause the specific biological accumulation and biotoxicity (Elliott et al., 1986; Nirmalakhandan et al., 1994; Posthuma et al., 1997). A study on the mixture toxicity of heavy metals demonstrated that the majority of sites were already immobilized by ions to the extent that new chemicals hardly have a chance to bind, such a strategy may fluctuate the joint effect along with competition for toxic sites among coexisting cations (Zeng et al., 2015).

Other Factors Affect Toxicity Prediction

It has been mentioned, the joint effect of chemical mixtures varies due to its type of compounds and toxic ratios (at their equitoxic ratios or non-equitoxic ratios). The multi-chemical mixture may show dissimilar toxicity response with a different ratio of the same chemicals (Kar and Leszczynski, 2019). For narcotic toxicants, funnel hypothesis was put forward to interpret the fact that as the number of components in a mixture increase, the range of deviation from toxic additivity decreases (Rayburn et al., 1991; Warne and Hawker, 1995) (Figures 8A,B). With the number of components increases, the joint effect tends to show additivity effect and toxicity reduction can be regarded as a result of dissolution and dilution (Altenburger et al., 2000; Backhaus et al., 2000; Faust et al., 2001). The equitoxic mixtures of the reactive chemicals also have the same tendency, mixtures effect is close to additive as the increased number of the components (Tian et al., 2012b). It is necessary to note that there are inconsistent with funnel hypothesis as some ternary mixtures effect showed more synergistic or antagonistic than binary mixtures effect (Cedergreen et al., 2012; Chen et al., 2014, 2015; Wang et al., 2015b,c). Besides, the fishing hypothesis was proposed to explain the variation rules of the joint effects of cyanogenic toxicants and aldehydes (Figure 8C; Li et al., 2014). The hypothesis elucidated why the interactions based on the equitoxic ratios of the different chemicals were stronger than the effects at non-equitoxic ratios. In this hypothesis, the CN acted as fishhooks and aldehydes acted as fish. A hook only caught one fish when fishing. There is a competition among different types of fish for one hook during the process of fishing, while no competitions among different types of fish. However, there is no uniform explanation for the prediction of different types of interactions in the multi-component mixtures, which remain a rough challenge for ecotoxicology.

Figure 8

Conclusion and Future Perspectives

Compared with other biological assay methods, luminescent bacteria toxicity assay has the advantages of simple operation, short test time and high repeatability. The results of numerous researches have illustrated the joint effect of multi-component mixtures and described the various types of models involved in joint toxicity prediction based on bacteria toxicity assay. Other studies also attempt to evaluate the toxicity of specific toxins by combining luminescent bacteria with other means such as dual detections by a combination of luminescent bacteria and probe. However, the effect of major factors such as the number of mixed components, the dominating components, and the toxic ratio of individual toxicants, on the joint effect of multi-component mixtures should be considered necessary. To address variable co-exposure scenarios, further research is required to decipher more related indicators to construct a fitting model for the assessment of multi-component mixtures. It would contribute to promoting thus further use of effect-based prediction of ecological environmental quality. Furthermore, the current luminescent bacteria toxicity assay is mostly used for acute toxicity testing, the determination of chronic toxicity exposure based on luminescent bacteria should be refined within the larger regulatory network. The differences between chronic and acute toxicity require further elucidation.

Statements

Author contributions

All authors wrote the manuscript jointly. DW, SW, and LB contributed to reviewing the literature and the write up of the manuscript. MN, SL, and WY contributed to the editing of the manuscript.

Funding

This work was supported by the National Key R&D Program of China (Grant number 2019YFB2103003), the National Science Foundation of China (Grant number 31670512), and Natural Science Basic Research Plan in Shaanxi Province of China (Grant number 2018JM3039).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  • 1

    AammiS.KaracaF.PetekM. (2017). A toxicological and genotoxicological indexing study of ambient aerosols (PM2.5-10) using in vitro bioassays. Chemosphere174, 490498. 10.1016/j.chemosphere.2017.01.141

  • 2

    AhnJ. M.HwangE. T.YounC. H.BanuD. L.KimB. C.NiaziJ. H.et al. (2009). Prediction and classification of the modes of genotoxic actions using bacterial biosensors specific for DNA damages. Biosens. Bioelectron.25, 767772. 10.1016/j.bios.2009.08.025

  • 3

    AlbertsB.JohnsonA.LewisJ. (1989). Molecular Biology of the Cell.New York, NY: Garland Science Press.

  • 4

    AltenburgerR.BackhausT.BoedekerW.FaustM.ScholzeM.GrimmeL. H. (2000). Predictability of the toxicity of multiple chemical mixtures to Vibrio fischeri: mixtures composed of similarly acting chemicals. Environ. Toxicol. Chem.19, 23412347. 10.1002/etc.5620190926

  • 5

    AltenburgerR.NendzaM.SchuurmannG. (2003). Mixture toxicity and its modeling by quantitative structure-activity relationships. Environ. Toxicol. Chem.22, 19001915. 10.1897/01-386

  • 6

    AngelB. M.VallottonP.ApteS. C. (2015). On the mechanism of nanoparticulate CeO2 toxicity to freshwater algae. Aquat. Toxicol.168, 9097. 10.1016/j.aquatox.2015.09.015

  • 7

    BackhausT.FaustM. (2012). Predictive environmental risk assessment of chemical mixtures: a conceptual framework. Environ. Sci. Technol.46, 25642573. 10.1021/es2034125

  • 8

    BackhausT.ScholzeM.GrimmeL. H. (2000). The single substance and mixture toxicity of quinolones to the bioluminescent bacterium Vibrio fischeri. Aquat. Toxicol.49, 4961. 10.1016/S0166-445X(99)00069-7

  • 9

    BarrangerA.RanceG. A.AminotY.DallasL. J.SforziniS.WestonN. J.et al. (2019). An integrated approach to determine interactive genotoxic and global gene expression effects of multiwalled carbon nanotubes (MWCNTs) and benzo[a]pyrene (BaP) on marine mussels: evidence of reverse ‘Trojan Horse’ effects. Nanotoxicology13, 13241343. 10.1080/17435390.2019.1654003

  • 10

    BarronM. G.ConmyR. N.HolderE. L.MeyerP.WilsonG. J.PrincipeV. E.et al. (2018). Toxicity of cold lake blend and Western Canadian select dilbits to standard aquatic test species. Chemosphere191, 16. 10.1016/j.chemosphere.2017.10.014

  • 11

    BelkinS. (2003). Microbial whole-cell sensing systems of environmental pollutants. Curr. Opin. Microbiol.6, 206212. 10.1016/S1369-5274(03)00059-6

  • 12

    BerenbaumM. C. (1981). Criteria for analyzing interactions between biologically-active agents. Adv. Cancer Res.35, 269335. 10.1016/S0065-230X(08)60912-4

  • 13

    BergnerT.TabibC. R.WinklerA.StipsitsS.KayerH.LeeJ.et al. (2015). Structural and biochemical properties of LuxF from Photobacterium leiognathi. Biochim. Biophys. Acta. 1854(10 Pt A), 14661475. 10.1016/j.bbapap.2015.07.008

  • 14

    BlissC. I. (1939). The toxicity of poisons applied jointly. Ann. Appl. Biol.26, 585615. 10.1111/j.1744-7348.1939.tb06990.x

  • 15

    BourgoisJ. J.SluseF. E.BaguetF.MallefetJ. (2001). Kinetics of Light Emission and Oxygen Consumption by Bioluminescent Bacteria. J. Bioenerg. Biomembr.33, 353363. 10.1023/A:1010615508916

  • 16

    BoydE. M.KillhamK.WrightJ.RumfordS.HetheridgeM.CummingR.et al. (1997). Toxicity assessment of xenobiotic contaminated groundwater using lux modified Pseudomonas fluorescens. Chemosphere35, 19671985. 10.1016/S0045-6535(97)00271-3

  • 17

    BrodlE.WinklerA.MacherouxP. (2018). Molecular mechanisms of bacterial bioluminescence. Comput. Struct. Biotechnol. J. 16, 551564. 10.1016/j.csbj.2018.11.003

  • 18

    CarbajoJ. B.Perdigon-MelonJ. A.PetreA. L.RosalR.LetonP.Garcia-CalvoE. (2015). Personal care product preservatives: risk assessment and mixture toxicities with an industrial wastewater. Water Res.72, 174185. 10.1016/j.watres.2014.12.040

  • 19

    CasseeF. R.GrotenJ. P.van BladerenP. J.FeronV. J. (1998). Toxicological evaluation and risk assessment of chemical mixtures. Crit. Rev. Toxicol.28, 73101. 10.1080/10408449891344164

  • 20

    CedergreenN.SørensenH.SvendsenC. (2012). Can the joint effect of ternary mixtures be predicted from binary mixture toxicity results?Sci. Total Environ. 427–428, 229237. 10.1016/j.scitotenv.2012.03.086

  • 21

    ChangC. M.OuY. H.LiuT. C.LuS. Y.WangM. K. (2016). A quantitative structure-activity relationship approach for assessing toxicity of mixture of organic compounds. SAR QSAR Environ. Res.27, 441453. 10.1080/1062936X.2016.1207204

  • 22

    ChenC.WangY.QianY.ZhaoX.WangQ. (2015). The synergistic toxicity of the multiple chemical mixtures: implications for risk assessment in the terrestrial environment. Environ. Int.77, 95105. 10.1016/j.envint.2015.01.014

  • 23

    ChenC.WangY.ZhaoX.WangQ.QianY. (2014). Comparative and combined acute toxicity of butachlor, imidacloprid and chlorpyrifos on earthworm, Eisenia fetida. Chemosphere100, 111115. 10.1016/j.chemosphere.2013.12.023

  • 24

    ChenC. Y.HuangC. F. (1996). Toxicity of organic mixtures containing cyanogenic toxicants. Environ. Toxicol. Chem.15, 14641469. 10.1002/etc.5620150906

  • 25

    ChenQ.LiY.LiB. (2020). Is color a matter of concern during microplastic exposure to Scenedesmus obliquus and Daphnia magna?J. Hazard. Mater.383:121224. 10.1016/j.jhazmat.2019.121224

  • 26

    ChenY. H.QinL. T.MoL. Y.ZhaoD. N.ZengH. H.LiangY. P. (2019). Synergetic effects of novel aromatic brominated and chlorinated disinfection byproducts on Vibrio qinghaiensis sp.-Q67. Environ. Pollut.250, 375385. 10.1016/j.envpol.2019.04.009

  • 27

    ChouT. C. (2006). Theoretical basis, experimental design, and computerized simulation of synergism and antagonism in drug combination studies. Pharmacol. Rev.58, 621681. 10.1124/pr.58.3.10

  • 28

    ClemowY. H.WilkieM. P. (2015). Effects of Pb plus Cd mixtures on toxicity, and internal electrolyte and osmotic balance in the rainbow trout (Oncorhynchus mykiss). Aquat. Toxicol.161, 176188. 10.1016/j.aquatox.2015.01.032

  • 29

    ConfortiF.IoeleG.StattiG. A.MarrelliM.RagnoG.MenichiniF. (2008). Antiproliferative activity against human tumor cell lines and toxicity test on Mediterranean dietary plants. Food Chem. Toxicol.46, 33253332. 10.1016/j.fct.2008.08.004

  • 30

    CostaM. B.TavaresF. V.MartinezC. B.ColaresI. G.MartinsC. M. G. (2018). Accumulation and effects of copper on aquatic macrophytes Potamogeton pectinatus L.: potential application to environmental monitoring and phytoremediation. Ecotoxicol. Environ. Saf.155, 117124. 10.1016/j.ecoenv.2018.01.062

  • 31

    De La Torre-RocheR.HawthorneJ.DengY.XingB.CaiW.NewmanL. A.et al. (2012). Fullerene-enhanced accumulation of p,p '-DDE in agricultural crop species. Environ. Sci. Technol.46, 93159323. 10.1021/es301982w

  • 32

    De SchepperW.DriesJ.GeuensL.BlustR. (2012). Sequential partial ozonation for the treatment of wastewater concentrate: practical implications from a conventional and (Eco) toxicological perspective. Ozone Sci. Eng.34, 163173. 10.1080/01919512.2012.663328

  • 33

    DeevaA. A.ZykovaE. A.NemtsevaE. V.KratasyukV. A. (2019). Functional divergence between evolutionary-related LuxG and Fre oxidoreductases of luminous bacteria. Proteins87, 723729. 10.1002/prot.25696

  • 34

    DingS.WuJ.ZhangM.LuH.MahmoodQ.ZhengP. (2015). Acute toxicity assessment of ANAMMOX substrates and antibiotics by luminescent bacteria test. Chemosphere140, 174183. 10.1016/j.chemosphere.2015.03.057

  • 35

    ElliottN. G.SwainR.RitzD. A. (1986). Metal interaction during accumulation by the mussel mytilus-edulis-planulatus. Mar. Biol.93, 395399. 10.1007/BF00401107

  • 36

    EngebrechtJ.SilvermanM. (1984). Identification of genes and gene-products necessary for bacterial bioluminescence. Proc. Natl. Acad. Sci. U.S.A.81, 41544158. 10.1073/pnas.81.13.4154

  • 37

    EPA (2000). Supplementary Guidance for Conducting Health Risk Assessment of Chemical Mixtures. Washington, DC: US Environmental Protection Agency.

  • 38

    FaustM.AltenburgerR.BackhausT.BlanckH.BoedekerW.GramaticaP.et al. (2001). Predicting the joint algal toxicity of multi-component s-triazine mixtures at low-effect concentrations of individual toxicants. Aquat. Toxicol.56, 1332. 10.1016/S0166-445X(01)00187-4

  • 39

    FengJ. N.GuoX. P.ChenY. R.LuD. P.NiuZ. S.TouF. Y.et al. (2020). Time-dependent effects of ZnO nanoparticles on bacteria in an estuarine aquatic environment. Sci. Total Environ.698:134298. 10.1016/j.scitotenv.2019.134298

  • 40

    FlokstraB. R.AkenB. V.SchnoorJ. L. (2008). Microtox® toxicity test: detoxification of TNT and RDX contaminated solutions by poplar tissue cultures. Chemosphere71, 19701976. 10.1016/j.chemosphere.2007.12.020

  • 41

    FraserT. R. (1872). Lecture on the antagonism between the actions of active substances. Br. Med. J.2, 485487. 10.1136/bmj.2.618.485

  • 42

    FulladosaE.MuratJ. C.MartinezM.VillaescusaI. (2004). Effect of pH on arsenate and arsenite toxicity to luminescent bacteria (Vibrio fischeri). Arch. Environ. Contam. Toxicol.46, 176182. 10.1007/s00244-003-2291-7

  • 43

    FulladosaE.MuratJ. C.VillaescusaI. (2005). Study on the toxicity of binary equitoxic mixtures of metals using the luminescent bacteria Vibrio fischeri as a biological target. Chemosphere58, 551557. 10.1016/j.chemosphere.2004.08.007

  • 44

    FuquaW. C.WinansS. C.GreenbergE. P. (1994). Quorum sensing in bacteria - the luxr-luxi family of cell density-responsive transcriptional regulators. J. Bacteriol.176, 269275. 10.1128/JB.176.2.269-275.1994

  • 45

    GeH.LinZ.YaoZ.GaoY.CongY.YuH. (2014). Balance between herbicidal activity and toxicity effect: a case study of the joint effects of triazine and phenylurea herbicides on Selenastrum capricornutum and Photobacterium phosphoreum. Aquat. Toxicol.150, 165174. 10.1016/j.aquatox.2014.03.007

  • 46

    GinerB.LafuenteC.LapenaD.ErrazquinD.LombaL. (2020). QSAR study for predicting the ecotoxicity of NADES towards Aliivibrio fischeri. Exploring the use of mixing rules. Ecotoxicol. Environ. Saf.191:110004. 10.1016/j.ecoenv.2019.110004

  • 47

    GodoyA. A.KummrowF. (2017). What do we know about the ecotoxicology of pharmaceutical and personal care product mixtures? A critical review. Crit. Rev. Environ. Sci. Technol.47, 14531496. 10.1080/10643389.2017.1370991

  • 48

    GodoyA. A.OliveiraA. C.SilvaJ. G. M.AzevedoC. C. J.DominguesI.NogueiraA. J. A.et al. (2019). Single and mixture toxicity of four pharmaceuticals of environmental concern to aquatic organisms, including a behavioral assessment. Chemosphere235, 373382. 10.1016/j.chemosphere.2019.06.200

  • 49

    GrotenJ. P.FeronV. J.SuhnelJ. (2001). Toxicology of simple and complex mixtures. Trends Pharmacol. Sci.22, 316322. 10.1016/S0165-6147(00)01720-X

  • 50

    GuM. B.MitchellR. J.KimB. C. (2004). Whole-cell-based biosensors for environmental biomonitoring and application. Adv. Biochem. Eng. Biotechnol.87, 269305. 10.1007/b13533

  • 51

    HagemanK. J.AebigC. H. F.LuongK. H.KaserzonS. L.WongC. S.ReeksT.et al. (2019). Current-use pesticides in New Zealand streams: comparing results from grab samples and three types of passive samplers. Environ. Pollut.254(Pt A), 112973. 10.1016/j.envpol.2019.112973

  • 52

    Hamdi El NajjarN.TouffetA.DebordeM.JournelR.Karpel Vel LeitnerN. (2014). Kinetics of paracetamol oxidation by ozone and hydroxyl radicals, formation of transformation products and toxicity. Sep. Purif. Technol. 136, 137143. 10.1016/j.seppur.2014.09.004

  • 53

    HanschC.MaloneyP. P.FujitaT. (1962). Correlation of biological activity of phenoxyacetic acids with hammett substituent constants and partition coefficients. Nature194:178. 10.1038/194178b0

  • 54

    HernandoM. D.EjerhoonM.Fernández-AlbaA. R.ChistiY. (2003). Combined toxicity effects of MTBE and pesticides measured with Vibrio fischeri and Daphnia magna bioassays. Water Res.37, 40914098. 10.1016/S0043-1354(03)00348-8

  • 55

    HooverG.KarS.GuffeyS.LeszczynskiJ.SepulvedaM. S. (2019). In vitro and in silico modeling of perfluoroalkyl substances mixture toxicity in an amphibian fibroblast cell line. Chemosphere233, 2533. 10.1016/j.chemosphere.2019.05.065

  • 56

    HuangL.LiS.LingX.ZhangJ.QinA.ZhuangJ.et al. (2019a). Dual detection of bioaccumulated Hg(2+) based on luminescent bacteria and aggregation-induced emission. Chem. Commun.55, 74587461. 10.1039/C9CC02782C

  • 57

    HuangR. Y.PeiL.LiuQ.ChenS.DouH.ShuG.et al. (2019b). Isobologram analysis: a comprehensive review of methodology and current research. Front. Pharmacol.10:1222. 10.3389/fphar.2019.01222

  • 58

    InouyeS. (1994). NAD(P)H-flavin oxidoreductase from the bioluminescent bacteria, Vibrio-fischeri ATCC-7744, is a flavoprotein. FEBS Lett.347, 163168. 10.1016/0014-5793(94)00528-1

  • 59

    IveyC. D.BesserJ. M.SteevensJ. A.WaltherM. J.MeltonV. D. (2019). Influence of dissolved organic carbon on the acute toxicity of copper and zinc to white sturgeon (Acipenser transmontanus) and a cladoceran (Ceriodaphnia dubia). Environ. Toxicol. Chem.38, 26822687. 10.1002/etc.4592

  • 60

    JarqueS.MasnerP.KlanovaJ.ProkesR.BlahaL. (2016). Bioluminescent vibrio fischeri assays in the assessment of seasonal and spatial patterns in toxicity of contaminated river sediments. Front. Microbiol.7:1738. 10.3389/fmicb.2016.01738

  • 61

    JeliazkovaN.JaworskaJ.WorthA. P. (2010). Open source tools for read-across and category formation. Silico Toxicol.17, 408445. 10.1039/9781849732093-00408

  • 62

    JiangB.LiG.XingY.ZhangD.JiaJ.CuiZ.et al. (2017). A whole-cell bioreporter assay for quantitative genotoxicity evaluation of environmental samples. Chemosphere184, 384392. 10.1016/j.chemosphere.2017.05.159

  • 63

    JinH.WangC.ShiJ.ChenL. (2014). Evaluation on joint toxicity of chlorinated anilines and cadmium to 0Photobacterium phosphoreum and QSAR analysis. J. Hazard. Mater.279, 156162. 10.1016/j.jhazmat.2014.06.068

  • 64

    JunghansM.BackhausT.FaustM.ScholzeM.GrimmeL. (2006). Application and validation of approaches for the predictive hazard assessment of realistic pesticide mixtures. Aquat. Toxicol.76, 93110. 10.1016/j.aquatox.2005.10.001

  • 65

    KarS.LeszczynskiJ. (2019). Exploration of computational approaches to predict the toxicity of chemical mixtures. Toxics7:15. 10.3390/toxics7010015

  • 66

    KaweeteerawatC.ChangC. H.RoyK. R.LiuR.LiR.TosoD.et al. (2015). Cu Nanoparticles have different impacts in escherichia coli and lactobacillus brevis than their microsized and ionic analogues. Acs Nano9, 72157225. 10.1021/acsnano.5b02021

  • 67

    KhanP. M.KarS.RoyK. (2020). Ecotoxicological QSARs of mixtures, in Ecotoxicological QSARs, ed RoyK.. (Humana Imprint Press), 437475. 10.1007/978-1-0716-0150-1_19

  • 68

    KimJ.KimS.SchaumannG. E. (2014). Development of a partial least squares-based integrated addition model for predicting mixture toxicity. Hum. Ecol. Risk Assess.20, 174200. 10.1080/10807039.2012.754312

  • 69

    KonemannH. (1981). Fish toxicity tests with mixtures of more than 2 chemicals - a proposal for a quantitative approach and experimental resultS. Toxicology19, 229238. 10.1016/0300-483X(81)90132-3

  • 70

    LM. L. (1977). Aquatic toxicology and hazard evaluation. Baltim. Am. Soc. Test. Mater.99108.

  • 71

    LiS.AndersonT. A.MaulJ. D.ShresthaB.GreenM. J.Canas-CarrellJ. E. (2013a). Comparative studies of multi-walled carbon nanotubes (MWNTs) and octadecyl (C18) as sorbents in passive sampling devices for biomimetic uptake of polycyclic aromatic hydrocarbons (PAHs) from soils. Sci. Total Environ.461–462, 560567. 10.1016/j.scitotenv.2013.05.048

  • 72

    LiS.IrinF.AtoreF. O.GreenM. J.Canas-CarrellJ. E. (2013b). Determination of multi-walled carbon nanotube bioaccumulation in earthworms measured by a microwave-based detection technique. Sci. Total Environ.445–446, 913. 10.1016/j.scitotenv.2012.12.037

  • 73

    LiW.TianD.LinZ.WangD.YuH. (2014). Study on the variation rules of the joint effects for multicomponent mixtures containing cyanogenic toxicants and aldehydes based on the transition state theory. J. Hazard. Mater.267, 98108. 10.1016/j.jhazmat.2013.12.035

  • 74

    LimbachL. K.WickP.ManserP.GrassR. N.BruininkA.StarkW. J. (2007). Exposure of engineered nanoparticles to human lung epithelial cells: influence of chemical composition and catalytic activity on oxidative stress. Environ. Sci. Technol.41, 41584163. 10.1021/es062629t

  • 75

    LinZ. F.NiuX. J.WuC.YinK. D.CaiZ. P. (2005). Prediction of the toxicological joint effects between cyanogenic toxicants and aldehydes to Photobacterium phosphoreum. QSAR Comb. Sci.24, 354363. 10.1002/qsar.200430882

  • 76

    LinZ. F.YuH. X.WeiD. B.WangG. H.FengJ. F.WangL. S. (2002). Prediction of mixture toxicity with its total hydrophobicity. Chemosphere46, 305310. 10.1016/S0045-6535(01)00083-2

  • 77

    LinZ. F.ZhongP.YinK. D.WangL. S.YuH. X. (2003). Quantification of joint effect for hydrogen bond and development of QSARs for predicting mixture toxicity. Chemosphere52, 11991208. 10.1016/S0045-6535(03)00329-1

  • 78

    LiuS. S.SongX. Q.LiuH. L.ZhangY. H.ZhangJ. (2009). Combined photobacterium toxicity of herbicide mixtures containing one insecticide. Chemosphere75, 381388. 10.1016/j.chemosphere.2008.12.026

  • 79

    LiuY.NieY.WangJ.WangJ.WangX.ChenS.et al. (2018). Mechanisms involved in the impact of engineered nanomaterials on the joint toxicity with environmental pollutants. Ecotoxicol. Environ. Saf.162, 92102. 10.1016/j.ecoenv.2018.06.079

  • 80

    LoeweS.MuischnekH. (1926). Combinated effects I Announcement - Implements to the problem. Naunyn Schmiedebergs Arch. Exp. Pathol. Pharmakol.114, 313326. 10.1007/BF01952257

  • 81

    Lütke EverslohC.HenningN.SchulzM.TernesT. A. (2014). Electrochemical treatment of iopromide under conditions of reverse osmosis concentrates – Elucidation of the degradation pathway. Water Res.48, 237246. 10.1016/j.watres.2013.09.035

  • 82

    MaW.HanY.XuC.HanH.ZhuH.LiK.et al. (2018). Biotoxicity assessment and toxicity mechanism on coal gasification wastewater (CGW): a comparative analysis of effluent from different treatment processes. Sci. Total Environ.637–638, 18. 10.1016/j.scitotenv.2018.04.404

  • 83

    MaX. Y.WangX. C.NgoH. H.GuoW.WuM. N.WangN. (2014). Bioassay based luminescent bacteria: interferences, improvements, and applications. Sci. Total Environ.468–469, 111. 10.1016/j.scitotenv.2013.08.028

  • 84

    Martín-de-LucíaI.Campos-MañasM. C.AgüeraA.Rodea-PalomaresI.Pulido-ReyesG.LeganésF.et al. (2017). Reverse Trojan-horse effect decreased wastewater toxicity in the presence of inorganic nanoparticles. Environ. Sci. Nano.4, 12731282. 10.1039/C6EN00708B

  • 85

    MartinezD. S. T.AlvesO. L.BarbieriE.MartinezD. S. T.AlvesO. L.BarbieriE. (2013). Carbon nanotubes enhanced the lead toxicity on the freshwater fish. J. Phys. Conf. Ser.429:012043. 10.1088/1742-6596/429/1/012043

  • 86

    MathesonI. B. C.LeeJ. (1981). An efficient bacterial bioluminescence with reduced lumichromE. Biochem. Biophys. Res. Commun.100, 532536. 10.1016/S0006-291X(81)80209-4

  • 87

    MeadowsS. L.GenningsC.CarterW. H.BaeD. S. (2002). Experimental designs for mixtures of chemicals along fixed ratio rays. Environ. Health Perspect.110, 979983. 10.1289/ehp.02110s6979

  • 88

    MiyashiroT.RubyE. G. (2012). Shedding light on bioluminescence regulation in Vibrio fischeri. Mol. Microbiol.84, 795806. 10.1111/j.1365-2958.2012.08065.x

  • 89

    MoL. Y.LiuJ.QinL. T.ZengH.-H.LiangY. P. (2017). Two-stage prediction on effects of mixtures containing phenolic compounds and heavy metals on vibrio qinghaiensis sp Q67. Bull. Environ. Contam. Toxicol.99, 1722. 10.1007/s00128-017-2099-1

  • 90

    MoL.ZhuZ.ZhuY.ZengH.LiY. (2014). Prediction and evaluation of the mixture toxicity of twelve phenols and ten anilines to the freshwater photobacterium vibrio qinghaiensis sp.-Q67. J. Chem.9:728254. 10.1155/2014/728254

  • 91

    MonossonE. (2005). Chemical mixtures: considering the evolution of toxicology and chemical assessment. Environ. Health Perspect.113, 383390. 10.1289/ehp.6987

  • 92

    MooreS. A.JamesM. N. G. (1995). Structural refinement of the Nonfluorescent flavoprotein from photobacterium-leiognathi AT 1.60 angstrom resolution. J. Mol. Biol.249, 195214. 10.1006/jmbi.1995.0289

  • 93

    MoriI. C.Arias-BarreiroC. R.KoutsaftisA.OgoA.KawanoT.YoshizukaK.et al. (2015). Toxicity of tetramethylammonium hydroxide to aquatic organisms and its synergistic action with potassium iodide. Chemosphere120, 299304. 10.1016/j.chemosphere.2014.07.011

  • 94

    MwenseM.WangX. Z.BuontempoF. V.HoranN.YoungA.OsbornD. (2004). Prediction of noninteractive mixture toxicity of organic compounds based on a fuzzy set method. J. Chem. Inf. Comput. Sci.44, 17631773. 10.1021/ci0499368

  • 95

    NaaszS.AltenburgerR.KuhnelD. (2018). Environmental mixtures of nanomaterials and chemicals: the Trojan-horse phenomenon and its relevance for ecotoxicity. Sci. Total Environ.635, 11701181. 10.1016/j.scitotenv.2018.04.180

  • 96

    NeboutP.CagnonB.DelpeuxS.Di GiustoA.ChedevilleO. (2016). Comparison of the efficiency of adsorption, ozonation, and ozone/activated carbon coupling for the removal of pharmaceuticals from water. J. Environ. Eng.142:42. 10.1061/(ASCE)EE.1943-7870.0001042

  • 97

    NijvipakulS.BallouD. P.ChaiyenP. (2010). Reduction kinetics of a flavin oxidoreductase LuxG from Photobacterium leiognathi (TH1): half-sites reactivity. Biochemistry49, 92419248. 10.1021/bi1009985

  • 98

    NirmalakhandanN.ArulgnanendranV.MohsinM.SunB.CadenaF. (1994). Toxicity of mixtures of organic-chemicals to microorganisms. Water Res.28, 543551. 10.1016/0043-1354(94)90005-1

  • 99

    O'GradyE. AWimpeeC. F. (2008). Mutations in the lux operon of natural dark mutants in the genus Vibrio. Appl. Environ. Microbiol.74, 6166. 10.1128/AEM.01199-07

  • 100

    OxenderD.FoxC. F. (1978). Molecular aspects of membrane transport. J. Supramol. Struct.9, 120167. 10.1002/jss.400090505

  • 101

    Pino-OtínM. R.BallesteroD.NavarroE.González-ColomaA.ValJ.MainarA. M. (2019). Ecotoxicity of a novel biopesticide from Artemisia absinthium on non-target aquatic organisms. Chemosphere216, 131146. 10.1016/j.chemosphere.2018.09.071

  • 102

    PlackettR. L.HewlettP. S. (1952). Quantal responses to mixtures of poisons. J. R. Statist. Soc. B, 14, 14163. 10.2307/2983865

  • 103

    PosthumaL.BaerselmanR.Van VeenR. P. M.Dirven-Van BreemenE. M. (1997). Single and joint toxic effects of copper and zinc on reproduction of Enchytraeus crypticus in relation to sorption of metals in soils. Ecotoxicol. Environ. Saf.38, 108121. 10.1006/eesa.1997.1568

  • 104

    PrakashJ.NirmalakhandanN.SunB.PeaceJ. (1996). Toxicity of binary mixtures of organic chemicals to microorganisms. Water Res.30, 14591463. 10.1016/0043-1354(95)00246-4

  • 105

    Pulido-ReyesG.LeganesF.Fernandez-PinasF.RosalR. (2017). Bio-nano interface and environment: a critical review. Environ. Toxicol. Chem.36, 31813193. 10.1002/etc.3924

  • 106

    QinL. T.ChenY. H.ZhangX.MoL. Y.ZengH.-H.LiangY.-P. (2018). QSAR prediction of additive and non-additive mixture toxicities of antibiotics and pesticide. Chemosphere198, 122129. 10.1016/j.chemosphere.2018.01.142

  • 107

    QinL. T.LiuS. S.ZhangJ.XiaoQ. F. (2011). A novel model integrated concentration addition with independent action for the prediction of toxicity of multi-component mixture. Toxicology280, 164172. 10.1016/j.tox.2010.12.007

  • 108

    QuadrosM. E.MarrL. C. (2012). Environmental and human health risks of aerosolized silver nanoparticles. J. Air Waste Manage. Assoc.60, 770781. 10.3155/1047-3289.60.7.770

  • 109

    RaJ. S.LeeB. C.ChangN. I.KimS. D. (2006). Estimating the combined toxicity by two-step prediction model on the complicated chemical mixtures from wastewater treatment plant effluents. Environ. Toxicol. Chem.25, 21072113. 10.1897/05-484R.1

  • 110

    RaoY.HanF.ChenQ.WangD.XueD.WangH.et al. (2019). Efficient degradation of diclofenac by LaFeO3-Catalyzed peroxymonosulfate oxidation—kinetics and toxicity assessment. Chemosphere218, 299307. 10.1016/j.chemosphere.2018.11.105

  • 111

    RayburnJ. R.FortD. J.McNewR.BantleJ. A. (1991). synergism and antagonism induced by 3 carrier solvents with t-retinoic acid and 6-aminonicotinamide using fetax. Bull. Environ. Contam. Toxicol.46, 625632. 10.1007/BF01688209

  • 112

    RenH.WangJ.ChenX.LiuJ. (2011). Primary research on single and combined toxicity of antibiotics in feeds by luminescent bacteria. J. Anim. Sci.21, 4953.

  • 113

    RenS.MeeR. W.FrymierP. D. (2004). Using factorial experiments to study the toxicity of metal mixtures. Ecotoxicol. Environ. Saf.59, 3843. 10.1016/S0147-6513(03)00099-X

  • 114

    RiboJ. M.RogersF. (1990). Toxicity of mixtures of aquatic contaminants using the luminescent bacteria bioassaY. Toxic. Assess.5, 135152. 10.1002/tox.2540050203

  • 115

    RoyK.KarS. (2016). In Silico Methods for Predicting Drug Toxicity, Methods in Molecular Biology. Berlin: Springer. Vol. 1425.

  • 116

    RoyK.KarS.DasR. N. (2015a). Understanding the Basics of QSAR for Applications in Pharmaceutical Sciences and Risk Assessment. Cambridge, MA: Academic Press.

  • 117

    RoyK.KarS.DasR. N. (2015b). A Primer on QSAR/QSPR Modeling. Springer Briefs in Molecular Science. (New York, NY: Springer Press). 10.1007/978-3-319-17281-1

  • 118

    SackeyL. N. A.KociV.van GestelC. A. M. (2020). Ecotoxicological effects on Lemna minor and Daphnia magna of leachates from differently aged landfills of Ghana. Sci. Total Environ.698:134295. 10.1016/j.scitotenv.2019.134295

  • 119

    SanchisJ.OlmosM.VincentP.FarreM.BarceloD. (2016). New insights on the influence of organic co-contaminants on the aquatic toxicology of carbon nanomaterials. Environ. Sci. Technol.50, 961969. 10.1021/acs.est.5b03966

  • 120

    SazykinI. S.SazykinaM. A.KhmelevtsovaL. E.MirinaE. A.KudeevskayaE. M.RogulinE. A.et al. (2016). Biosensor-based comparison of the ecotoxicological contamination of the wastewaters of Southern Russia and Southern Germany. Int. J. Environ. Sci. Technol. (Tehran)13, 945954. 10.1007/s13762-016-0936-0

  • 121

    SendraM.YesteP. M.Moreno-GarridoI.GaticaJ. M.BlascoJ. (2017). CeO2 NPs, toxic or protective to phytoplankton? Charge of nanoparticles and cell wall as factors which cause changes in cell complexity. Sci. Total Environ.590–591, 304315. 10.1016/j.scitotenv.2017.03.007

  • 122

    SharifianS.HomaeiA.HemmatiR.KhajehK. (2017). Light emission miracle in the sea and preeminent applications of bioluminescence in recent new biotechnology. J. Photochem. Photobiol. B172, 115128. 10.1016/j.jphotobiol.2017.05.021

  • 123

    SimonA.PreussT. G.SchaefferA.HollertH.MaesH. M. (2015). Population level effects of multiwalled carbon nanotubes in Daphnia magna exposed to pulses of triclocarban. Ecotoxicology24, 11991212. 10.1007/s10646-015-1479-8

  • 124

    SkottiE.AnastasakiE.KanellouG.PolissiouM.TarantilisP. A. (2014). Total phenolic content, antioxidant activity and toxicity of aqueous extracts from selected Greek medicinal and aromatic plants. Ind. Crops Prod.53, 4654. 10.1016/j.indcrop.2013.12.013

  • 125

    SousaC. A.SoaresH.SoaresE. V. (2019). Chronic exposure of the freshwater alga Pseudokirchneriella subcapitata to five oxide nanoparticles: hazard assessment and cytotoxicity mechanisms. Aquat. Toxicol.214:105265. 10.1016/j.aquatox.2019.105265

  • 126

    SpragueJ. B.RamsayB. A. (1965). Lethal levels mixed copper-zinc solutions for juvenile salmon. J. Fish. Res. Board Can.22:425. 10.1139/f65-042

  • 127

    StaveleyJ. P.GreenJ. W.NuszJ.EdwardsD.HenryK.KernM. (2018). Variability in nontarget terrestrial plant studies should inform endpoint selection. Integr. Environ. Assess. Manag.14, 639648. 10.1002/ieam.4055

  • 128

    SuL.ZhangX.YuanX.ZhaoY.ZhangD.QinW. (2012). Evaluation of joint toxicity of nitroaromatic compounds and copper to Photobacterium phosphoreum and QSAR analysis. J. Hazard. Mater.241, 450455. 10.1016/j.jhazmat.2012.09.065

  • 129

    SuY.YanX.PuY.XiaoF.WangD.YangM. (2013). Risks of single-walled carbon nanotubes acting as contaminants-carriers: potential release of phenanthrene in japanese medaka (Oryzias latipes). Environ. Sci. Technol.47, 47044710. 10.1021/es304479w

  • 130

    TanetL.TamburiniC.BaumasC.GarelM.SimonG.CasalotL. (2019). Bacterial bioluminescence: light emission in photobacterium phosphoreum is not under quorum-sensing control. Front. Microbiol.10:365. 10.3389/fmicb.2019.00365

  • 131

    ThorntonL. M.LeSueurM. C.YostA. T.StephensD. A.OrisJ. T.Sellin JeffriesM. K. (2017). Characterization of basic immune function parameters in the fathead minnow (Pimephales promelas), a common model in environmental toxicity testing. Fish Shellfish Immunol.61, 163172. 10.1016/j.fsi.2016.12.033

  • 132

    TianD.LinZ.YinD.ZhangY.KongD. (2012a). Atomic charges of individual reactive chemicals in binary mixtures determine their joint effects: an example of cyanogenic toxicants and aldehydes. Environ. Toxicol. Chem.31, 270278. 10.1002/etc.1701

  • 133

    TianD.LinZ.YuJ.YinD. (2012b). Influence factors of multicomponent mixtures containing reactive chemicals and their joint effects. Chemosphere88, 9941000. 10.1016/j.chemosphere.2012.03.043

  • 134

    TichyM.HanzlikovaI.RuckiM.PokornaA.UzlovaR.TumovaJ. (2008). Acute toxicity of binary mixtures: alternative methods, QSAR and mechanisms. Interdiscip. Toxicol.1, 1517. 10.2478/v10102-010-0025-z

  • 135

    TongF.ZhaoY.GuX.GuC.LeeC. C. (2015). Joint toxicity of tetracycline with copper(II) and cadmium(II) to Vibrio fischeri: effect of complexation reaction. Ecotoxicology24, 346355. 10.1007/s10646-014-1383-7

  • 136

    TsaiC. M.GazdarA. F.VenzonD. J.SteinbergS. M.DedrickR. L.MulshineJ. L.et al. (1989). Lack of invitro synergy between etoposide and cis-diammine-dichloroplatinum(II). Cancer Res.49, 23902397.

  • 137

    UrbanczykH.AstJ. C.HigginsM. J.CarsonJ.DunlapP. V. (2007). Reclassification of Vibrio fischeri, Vibrio logei, Vibrio salmonicida and vibrio wodanis as Aliivibrio fischeri gen. nov., comb. nov., Aliivibrio logei comb. nov., Aliivibrio salmonicida comb. nov. and Aliivibrio wodanis comb. nov. Int. J. Syst. Evol. Microbiol.57(Pt 12), 28232829. 10.1099/ijs.0.65081-0

  • 138

    VerhaarH. J. M.BusserF. J. M.HermensJ. L. M. (1995). Surrogate parameter for the base-line toxicity content of contaminated water - simulating the bioconcentration of mixtures of pollutants and counting molecules. Environ. Sci. Technol.29, 726734. 10.1021/es00003a021

  • 139

    VetrovaE. V.KudryashevaN. S.KratasyukV. A. (2007). Redox compounds influence on the NAD(P)H:FMN-oxidoreductase-luciferase bioluminescent system. Photochem. Photobiol. Sci. 6, 3540. 10.1039/B608152E

  • 140

    VighiM.MiglioratiS.MontiG. S. (2009). Toxicity on the luminescent bacterium Vibrio fischeri (Beijerinck). I: QSAR equation for narcotics and polar narcotics. Ecotoxicol. Environ. Saf.72, 154161. 10.1016/j.ecoenv.2008.05.008

  • 141

    VilelaP.JacomeG.KimS. Y.NamK.YooC. (2020). Population response modeling and habitat suitability of Cobitis choii fish species in South Korea for climate change adaptation. Ecotoxicol. Environ. Saf.189:109949. 10.1016/j.ecoenv.2019.109949

  • 142

    VillaS.VighiM.FinizioA. (2014). Experimental and predicted acute toxicity of antibacterial compounds and their mixtures using the luminescent bacterium Vibrio fischeri. Chemosphere108, 239244. 10.1016/j.chemosphere.2014.01.046

  • 143

    VouitsisE.NtziachristosL.PistikopoulosP.SamarasZ.ChrysikouL.SamaraC. (2009). An investigation on the physical, chemical and ecotoxicological characteristics of particulate matter emitted from light-duty vehicles. Environ. Pollut.157, 23202327. 10.1016/j.envpol.2009.03.028

  • 144

    WangD.WuX.LinZ.DingY. (2018). A comparative study on the binary and ternary mixture toxicity of antibiotics towards three bacteria based on QSAR investigation. Environ. Res.162, 127134. 10.1016/j.envres.2017.12.015

  • 145

    WangF.YaoJ.LiuH.LiuR.ChenH.YiZ.et al. (2015a). Cu and Cr enhanced the effect of various carbon nanotubes on microbial communities in an aquatic environment. J. Hazard. Mater.292, 137145. 10.1016/j.jhazmat.2015.03.032

  • 146

    WangN.WangX. C.MaX. (2015b). Characteristics of concentration–inhibition curves of individual chemicals and applicability of the concentration addition model for mixture toxicity prediction. Ecotoxicol. Environ. Saf.113, 176182. 10.1016/j.ecoenv.2014.12.008

  • 147

    WangY.ChenC.QianY.ZhaoX.WangQ. (2015c). Ternary toxicological interactions of insecticides, herbicides, and a heavy metal on the earthworm Eisenia fetida. J. Hazard. Mater.284, 233240. 10.1016/j.jhazmat.2014.11.017

  • 148

    WangY.ChenC.QianY.ZhaoX.WangQ.KongX. (2015d). Toxicity of mixtures of lambda-cyhalothrin, imidacloprid and cadmium on the earthworm Eisenia fetida by combination index (CI)-isobologram method. Ecotoxicol. Environ. Saf.111, 242247. 10.1016/j.ecoenv.2014.10.015

  • 149

    WangD.GeH.LinZ.LiS. (2016). Acute joint effects and mechanisms of binary toxicity of sulfonamides and metal oxide nanoparticles to vibrio fischeri. Asian J. Ecotoxicol.11, 6570. 10.7524/AJE.1673-5897.20160526001

  • 150

    WarneM. S. J.HawkerD. W. (1995). The number of components in a mixture determines whether synergistic and antagonistic or additive toxicity predominate - the funnel hypothesis. Ecotoxicol. Environ. Saf.31, 2328. 10.1006/eesa.1995.1039

  • 151

    WeiJ.SongY.MengX.PicJ.-S. (2015). Combination of Fenton oxidation and sequencing batch membrane bioreactor for treatment of dry-spun acrylic fiber wastewater. Environ. Earth Sci.73, 49114921. 10.1007/s12665-015-4091-6

  • 152

    WeiS.WangF.ChenY.LanT.ZhangS. (2018). The joint toxicity effect of five antibiotics and dibutyl phthalate to luminescent bacteria (Vibrio fischeri). Environ. Sci. Pollut. Res.25, 2650426511. 10.1007/s11356-018-2720-9

  • 153

    WuX.CobbinaS. J.MaoG.XuH.ZhangZ.YangL. (2016). A review of toxicity and mechanisms of individual and mixtures of heavy metals in the environment. Environ. Sci. Pollut. Res.23, 82448259. 10.1007/s11356-016-6333-x

  • 154

    XiongW.WuS.XuY.ZhuW. (2007). The joint toxicity effects of cadmium, chromium and lead on luminescent bacteria. Ecol. Environ.4, 10851087. 10.16258/j.cnki.1674-5906.2007.04.003

  • 155

    XuX.ChenJ.WangS.GeJ.QuR.FengM.et al. (2018). Degradation kinetics and transformation products of chlorophene by aqueous permanganate. Water Res.138, 293300. 10.1016/j.watres.2018.03.057

  • 156

    YangG.ChenC.WangY.PengQ.ZhaoH.GuoD.et al. (2017). Mixture toxicity of four commonly used pesticides at different effect levels to the epigeic earthworm, Eisenia fetida. Ecotoxicol. Environ. Saf.142, 2939. 10.1016/j.ecoenv.2017.03.037

  • 157

    YangX.JiY.WangF.XuJ.LiuX.MaK.et al. (2016). Comparison of organics and heavy metals acute toxicities to Vibrio fischeri. J. Serb. Chem. Soc.81, 697705. 10.2298/JSC151124011Y

  • 158

    YaoZ.LinZ.WangT.TianD.ZouX.GaoY.et al. (2013). Using molecular docking-based binding energy to predict toxicity of binary mixture with different binding sites. Chemosphere92, 11691176. 10.1016/j.chemosphere.2013.01.081

  • 159

    YuanX.LuG. H.ZhaoJ. S. (2002). QSAR study on the joint toxicity of 2,4-dinitrotoluene with aromatic compounds to Vibrio fischeri. J. Environ. Sci. Health, Part A: Toxic/Hazard. Subst. Environ. Eng.37, 573578. 10.1081/ESE-120003238

  • 160

    ZebB.PingZ.MahmoodQ.LinQ.PervezA.IrshadM.et al. (2016). Assessment of combined toxicity of heavy metals from industrial wastewaters on Photobacterium phosphoreum T3S. Appl. Water Sci.7, 20432050. 10.1007/s13201-016-0385-4

  • 161

    ZengY.WangL.JiangL.CaiX.LiY. (2015). Joint Toxicity of lead, chromium, cobalt and nickel to photobacterium phosphoreum at no observed effect concentration. Bull. Environ. Contam. Toxicol.95, 260264. 10.1007/s00128-015-1568-7

  • 162

    ZhangL.LiQ.ChenL.ZhangA.HeJ.WenZ.et al. (2015). Toxicity of surface water from Huangpu River to luminous bacteria (Vibrio qinghaiensis SP. Q67) and zebrafish (Danio rerio) embryos. Ecotoxicol. Environ. Saf.112, 137143. 10.1016/j.ecoenv.2014.10.037

  • 163

    ZhangS.SuL.ZhangX.LiC.QinW.ZhangD.et al. (2019a). Combined toxicity of nitro-substituted benzenes and zinc to photobacterium phosphoreum: evaluation and QSAR analysis. Int. J. Environ. Res. Public Health16:6. 10.3390/ijerph16061041

  • 164

    ZhangX.SunH.ZhangZ.NiuQ.ChenY.CrittendenJ. C. (2007). Enhanced bioaccumulation of cadmium in carp in the presence of titanium dioxide nanoparticles. Chemosphere67, 160166. 10.1016/j.chemosphere.2006.09.003

  • 165

    ZhangY.GuoX.SiX.YangR.ZhouJ.QuanX. (2019b). Environmentally persistent free radical generation on contaminated soil and their potential biotoxicity to luminous bacteria. Sci. Total Environ.687, 348354. 10.1016/j.scitotenv.2019.06.137

  • 166

    ZhangY. H.LiuS. S.SongX. Q.GeH. L. (2008). Prediction for the mixture toxicity of six organophosphorus pesticides to the luminescent bacterium Q67. Ecotoxicol. Environ. Saf.71, 880888. 10.1016/j.ecoenv.2008.01.014

  • 167

    ZhouX.SangW.LiuS.ZhangY.GeH. (2010). Modeling and prediction for the acute toxicity of pesticide mixtures to the freshwater luminescent bacterium Vibrio qinghaiensis sp.-Q67. J. Environ. Sci.22, 433440. 10.1016/S1001-0742(09)60126-1

  • 168

    ZhouS.ZhouJ.HuangG.YingG. (2018). Toxicity and joint effects of heavy metals zn, cu and hg to recombinant luminescent bacterium. J. South China Normal Univers. (Nat. Sci. Edn.)50, 3337. 10.6054/j.jscnun.2018008

  • 169

    ZhuW.WangJ.ChenX.ZhaxiC.YangY.SongY. (1994). A new species of luminescent bacterium: vibrio qinghaiensis sp.nov. Oceanol. Limnol. Sin.25, 273279.

  • 170

    ZhuY.ZhangY.ShiG.YangJ.ZhangJ.LiW.et al. (2015). Nanodiamonds act as Trojan horse for intracellular delivery of metal ions to trigger cytotoxicity. Part Fibre Toxicol. 12:2. 10.1186/s12989-014-0075-z

  • 171

    ZouX.LinZ.DengZ.YinD.ZhangY. (2012). The joint effects of sulfonamides and their potentiator on Photobacterium phosphoreum: differences between the acute and chronic mixture toxicity mechanisms. Chemosphere86, 3035. 10.1016/j.chemosphere.2011.08.046

Summary

Keywords

joint toxicity, luminescent bacteria, model, QSAR, mechanism

Citation

Wang D, Wang S, Bai L, Nasir MS, Li S and Yan W (2020) Mathematical Modeling Approaches for Assessing the Joint Toxicity of Chemical Mixtures Based on Luminescent Bacteria: A Systematic Review. Front. Microbiol. 11:1651. doi: 10.3389/fmicb.2020.01651

Received

01 April 2020

Accepted

25 June 2020

Published

31 July 2020

Volume

11 - 2020

Edited by

Gérald Thouand, Université de Nantes, France

Reviewed by

Uwe Strotmann, Westfälische Hochschule, Germany; Supratik Kar, Jackson State University, United States

Updates

Copyright

*Correspondence: Shanshan Li

This article was submitted to Microbiotechnology, a section of the journal Frontiers in Microbiology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics