Research Article  Open Access
Lingyun Mo, Zongqiang Zhu, Yinian Zhu, Honghu Zeng, Yanhong Li, "Prediction and Evaluation of the Mixture Toxicity of Twelve Phenols and Ten Anilines to the Freshwater Photobacterium Vibrio qinghaiensis sp.Q67", Journal of Chemistry, vol. 2014, Article ID 728254, 9 pages, 2014. https://doi.org/10.1155/2014/728254
Prediction and Evaluation of the Mixture Toxicity of Twelve Phenols and Ten Anilines to the Freshwater Photobacterium Vibrio qinghaiensis sp.Q67
Abstract
Twelve substituted phenols and ten substituted anilines were chosen to investigate if the dose addition and independent action models can be used to evaluate the mixture toxicity of phenolic and aniline derivatives (PADs). Their photobacterium toxicity to the freshwater luminescent bacterium Vibrio qinghaiensis sp.Q67 showed that the twoparameter Weibull or Logit function could be effectively applied to describe the doseresponse relationships. The joint toxicity of three equivalenteffect concentration ratio (EECR) mixtures and twelve uniform design concentration ratio (UDCR) mixtures could be well evaluated using the dose addition (DA) or the independent action (IA) model within 95% confidence intervals.
1. Introduction
Phenol and aniline derivatives (PADs) are important chemical materials and intermediates in synthesis of other compounds, such as dyes, pharmaceuticals, polymers, and synthetic resins. However, they are also the major carcinogenic pollutants in aquatic environments and extremely toxic to organisms [1, 2]. Many of PADs are priority pollutants in the water [3].
Phenol and aniline derivatives usually occur as mixture contaminants in the aquatic environments [4]. However, most of the toxicity data obtained are information about pure and individual chemicals, not related to mixtures. Therefore, the reliability is questionable to evaluate the multicomponent mixture toxicity from the toxicity data about single chemicals. Several mixtures in different toxicity tests were studied by some researchers to get better results in the predictive evaluation [5].
In order to evaluate the combined toxicity, it is required firstly to assume the relationships between the toxicity of individual chemicals and their mixtures quantitatively. Essentially, two different basic models are available for this purpose, that is, the independent action (IA) and the dose addition (DA) [5]. Other early analysis methods of mixture toxicity included concepts such as toxic units, mixture toxicity index, and additive index. On these two different basic models, the similarity parameter was built [4, 5]. Joint toxicity studies on chemical mixtures have been slowly progressing over the past decades, with preliminary researches involving twocomponent mixtures in equieffect concentration (most of which are EC_{50}) [4, 6]. Since the year 2000, the nonlinear simulation of the doseeffect curve and the recurrent points have been used to estimate the low effects concentration, and the dose addition and independent action for multicomponent mixture predictions [5, 7, 8]. The mixture toxicity of nonpolar narcotic chemicals could be predicted by two QSAR approaches, that is, the forward stepwise multilinear regression (MLR) and the radial basis function neural networks (RBFNNs) from molecular descriptors which are calculated and defined as composite descriptors according to the fractional concentrations of the mixture components [9]. In order to enable the prioritization of ionic liquids with a favourable cytotoxicity profile, a Classification and Regression Tree (CART) classifier has been recently proposed due to its simplicity, reliability, good predictive capability, and transparent chemical interpretation based on structural molecular fragments [10, 11]. Toxicological evaluation of chemical mixtures is now considered to have truly entered into the field of multicomponent mixtures [4, 6, 12–15].
In this work, the joint toxicities of 22 species of phenol and aniline derivatives (PADs) and the toxicities of their mixtures to the freshwater photoluminescent bacterium Vibrio qinghaiensis sp.Q67 (Q67) were measured through a microplate toxicity test procedure [16–19]. To seek how the mixture toxicity changes with changing concentrations of individual chemicals in the mixture effectively, the compounds were mixed in equieffect concentrations ratio (EECR) or in uniform design concentration ratio (UDCR).
2. Materials and Methods
2.1. Materials
The PADs are common in water environment and all have a high solubility. Some physicalchemical properties of 22 tested PADs are listed in Table 1. Three EECR mixtures and twelve UDCR mixtures of the test agents were made by the stock solutions, which were prepared with pure agents in distilled water, kept in critically cleaned glass containers, and stored at 4°C.
 
Common name; CAS RN: Chemical Abstracts Services Registry Number; MW: molecular weight. 
The freshwater photoluminescent bacterium Vibrio qinghaiensis sp.Q67 (Q67) was used as the test organism that was kindly provided by East China Normal University, China. The culture condition and media have been described in detail in some of our earlier works [16–19]. To prepare the culture medium, the following chemicals were dissolved into 1000 mL distilled water: 5.0 g tryptone, 5.0 g yeast extract, 3.0 g glycerin, 35.8 mg Na_{2}HPO_{4}·12H_{2}O, 13.6 mg KH_{2}PO_{4}, 33.0 mg CaCl_{2}, 1.54 g NaCl, 0.61 g MgCl_{2}·6H_{2}O, 1.34 g NaHCO_{3}, 0.25 g MgSO_{4}·7H_{2}O. The solution was adjusted to pH 8.8–9.0 and turned into a solid culture medium after adding 2% agar. 50 mL culture medium was put into a 100 mL conical flask, which was then obturated by using a brown paper and sterilized at 121°C for 20 min through high pressure steam. For each experiment, the photoluminescent bacteria Q67 were inoculated from a stock culture medium into a fresh culture medium, which was then maintained and cultured for 24 h at 22 ± 1°C. Finally, the bacteria Q67 were further cultured in a liquid medium at 22 ± 1°C by shaking (120 r/min) for 18 h [16, 18–20].
2.2. Methods
A Veritas luminometer with a 96well microplate was bought from the Turner BioSystems Inc., USA, and used to examine the toxicity of PADs and their mixtures. Regarding our earlier microplate toxicity test procedure, twelve concentration gradients were planned for each of the twelve phenols and ten anilines [16, 17]. 100 μL MilliQ water was injected into twelve wells of the first row of the microplate, which were regarded as twelve controls (without any toxicant). Twelve different toxicant solutions were injected into twelve wells of the second row of the microplate and then filled with MilliQ water to the total volume of 100 μL. To obtain a chemiluminescence inhibition (response) from 1% to the maximum, the twelve toxicant solutions were prepared after a proper dilution factor. Similarly, some other toxicant solutions were prepared and injected into twelve wells of the third or fourth row of the microplate. The fifth row of the microplate was used as a test duplication of the second row. Finally, each well of the microplate was filled with 100 μL bacterial Q67 suspension. All microplate tests were triplicated. Considering the sum of all component concentrations as the total concentration of a PADs mixture, the microplate toxicity test with the PADs mixture was conducted in the same manner as the procedure of a single PAD toxicant.
In each test, the photoluminescent bacteria Q67 were exposed to different treatments and controls, and the relative light unit (RLU) was measured. An inhibition ratio of the luminous intensity of the bacterium ( or ) was applied to express the PADs toxicity to the photoluminescent bacteria and can be calculated by the following equation: where is a mean RLUs value of a treatment and a mean RLUs value of all untreated controls.
2.3. Mixture Design
The test mixtures were prepared by both the EECR and UDCR methods. According to the measured results for the 22 PAD toxic pollutants, the best nonlinear simulation was performed and the given effect (%) was expressed as a concentration . For the three EECR mixtures, the EC_{1}, EC_{10}, and EC_{50} percentages of the total mixture concentration () were designated as the concentration ratios of each PAD in the mixtures, that is, EE2201, EE2210, and EE2250, respectively. The concentrations percentages of different PAD compounds in the mixtures were showed in Table 2.

To avoid the significant range limitation of the concentration percentages in EECR mixtures, 12 UDCR mixtures (denoted as UD2201, UD2202, UD2203, UD2204, UD2205, UD2206, UD2207, UD2208, UD2209, UD2210, UD2211, and UD2212) were designed using the former 11 lines in the uniform design table (Utable), U_{13}(13^{22}) (Table 2). The uniform design (UD) experiment is planned in terms of the uniform design table (Utable). Like orthogonal design, UD offers lots of experimental tables to the users for convenient utilization. Utables of form (), where the subscript refers to the number of the mixture experiments and the superscript to the number of pesticides (factors) in the experiment (mixture), are purposely chosen to mimic the tables of orthogonal designs, (), except that the number of levels equals the number of experiments () [16]. In light of U_{13}(13^{22}) given in Table 2, twelve 22component mixtures and a control group were designed. For each mixture design, proper concentrations were selected based on the results of some preexperiments and no uniform dilution factors were applied. As an effective experimental design method, the uniform design (UD) can save time and effort in investigating the changes of concentration in threedimension space. When the factors in the uniform design (the numbers of mixture components) and the levels of the factors (the concentration levels) are large, it is very necessary and efficient to use the uniform design (UD) to investigate the joint effect of a multiplecomponent mixture because the UD can reduce significantly the mixture experiment efforts.
By using the computer program APTox, the doseresponse curves (DRCs) were predicted based on both of the DA and IA models within a confidence interval of 95%; the root mean square error (RMSE_{x}) andthe correlation coefficient () of the models and the mixture values were also calculated [16, 20]. Mathematically, the DA model is expressed as where is the number of mixture components, the concentration of the th component that provokes % effect when applied singly, and the concentration of the th component in the mixture [16, 20]. The IA model can be formulated mathematically as where and are the total concentration and total effect of the mixture, respectively, and is the effect of the th component with the concentration of in the mixture [16, 20].
3. Results and Discussion
3.1. Photobacterial Toxicity of Single PADs
Based on the nonlinear leastsquares fitting the experimental concentrationresponse data of the 22 species of PADs, the best doseresponse curves (DRCs) were obtained. The DRCs of the 22 PADs to the photobacterium could be well described by the twoparameter function Weibull or Logit. According to the experimental results of the 12 concentration points designed with the microplate method and the experience rules that concentration points should be 56 times the parameter values of the DRCs model, the twoparameter model is markedly more effective than the threeparameter model in fitting the experimental data. Therefore, all DRCs simulations in the present paper have applied only two parameters with the Weibull and Logit model in the data analysis.
The Weibull and Logit functions can be expressed in (2) and (3), respectively,
From the nonlinear leastsquares fitting to the experimental data, the best simulated DRCs of the 22 PADs were obtained (Table 3). The results indicated that the doseresponse curves (DRCs) of the 22 PADs on the luminescent bacterium Q67 could be described by the twoparameter Weibull or Logit models. The RMSE, , and parameters of the fitting models are also calculated and given in Table 3. By substituting the and values into the Weibull or Logit functions, the values of the concentration () that inhibited the chemiluminescence by (%) at each effect level were calculated. In Table 3, the pEC_{50}, that is, the negative logarithm of the EC_{50} (−logEC_{50}), for 22 PADs are also given in comparison with classical toxic indices. Both of the Weibull and Logit functions could well describe the concentrationtoxicity data of 22 PADs with and , which indicated the high stability and good calibration ability of the fitting of the two models [19].

A believable concentrationresponse examination of individual chemicals is indispensable to evaluate mixture toxicities. As a typical example for the concentrationinhibition ratio relationship of the single substance, the experimental and fitting results for omethylaniline (PAD16) are illustrated in Figure 1(a). The photobacterial toxicity relationship for PAD16 was a nonsymmetrical and sigmoidal curve, which was sharper in the high inhibition section than in the low inhibition section, which was a general character of the 22 PADs tested. Additionally, Figure 1(a) reveals the relation between EC_{1} and EC_{50} in the low inhibition section. For the estimations, was chosen as a lower limit of the effect level, since the obtained test data commonly tolerated the interpolative approximation of inhibitions down to this level.
(a)
(b)
The nonlinear fitting with the regression model was depicted in Figure 1(b). The obtained concentrationresponse curves showed significant discrepancy in position and shape. In the strict mathematical term, however, the intersections among the curves obviously indicated that the curves did not occur in parallel with each other. It was very obvious, as the ratio between EC_{1} and EC_{50} was used as a gradient indicator in the low effect sections of the doseresponse curves. The ratio was diverse for the substance examined and varied from 0.00357 of pchlorophenol (PAD03) to 0.04138 of ochloroaniline (PAD20). In Table 3, for the extreme case PAD03, the EC_{50} was more than 280 times greater than the EC_{1}, but the EC_{50} was only 24 times higher than the EC_{1} for PAD20. The toxic dose range of PAD03 was the widest among the 22 compounds. The most toxic PAD was pnitroaniline (PAD17, ) and the least toxic PAD was aniline (PAD13, ). The toxicity order of the 22 PADs was pnitroaniline (PAD17, ) > 2,4dichlorophenol (PAD05, 3.66) > onitroaniline (PAD18, 3.65) > pchlorophenol (PAD03, 3.43) > pnitrophenol (PAD11, 3.29) > onitrophenol (PAD12, 3.20) > mnitrophenol (PAD10, 3.16) > 2,3dimethylphenol (PAD02, 3.08) > pchloroaniline (PAD22, 2.93) > ocresol (PAD08, 3.00) > pcresol (PAD09, 2.99) > mnitroaniline (PAD19, 2.93) > mchloroaniline (PAD21, 2.89) > mcresol (PAD07, 2.83) > ochlorophenol (PAD04, 2.81) > ochloroaniline (PAD20, 2.78) > 3,5dihydroxytoluene (PAD01, 2.77) > pmethylaniline (PAD15, 2.52) > phenol (PAD06, 2.47) > mmethylaniline (PAD14, 2.42) > omethylaniline (PAD16, 2.36) > aniline (PAD13, 2.04).
3.2. Photobacterial Toxicity of PADs Mixtures
The concentrationresponse relationships of three EECR mixtures and twelve UDCR mixtures were determined on the Veritas luminometer in the same procedure as described above. The three EECR mixtures (EE2201, EE2210, and EE2250) and 12 UDCR mixtures (UD2201UD2212) consisted of all 22 PADs. The experimental data of the mixture toxicities and the predictions of the DA and IA models are illustrated in Figures 2 and 3. The dependence of the inhibition of photobacterial luminosity on the total concentration of the mixture constituent () can be expressed by the concentrationresponse equations.
The fitting parameters of the Weibull or Logit models, , , RSME, and _{,} were calculated for the toxicity prediction of the 15 PADs mixtures (Table 4). The Weibull or Logit functions could well fit the concentrationtoxicity data of the 15 PADs mixtures with and , which indicated the high stability and good calibration ability of the fitting of the two models. The values of the 15 PADs mixtures varied from 2.64 for UD2203 to 2.94 for UD2209, which were within the range from the least toxic PAD (2.04 for PAD13) to the most toxic PAD (3.68 for PAD17), indicating the absence of strong synergistic or antagonistic interactions within the mixture [13].

A mixture concentration can be reflected only within a certain space of a direction for the equieffect concentrations ratio (EECR) method; that is, it cannot describe the whole situation and can only represent some alike directions of space concentrations because of the existing relationship among the effect concentrations (), even though different EECR mixtures were measured. Instead, the PAD mixtures can be designed by the uniform design concentration ratio (UDCR) method to investigate the toxicity rule at the various space concentrations of PADs [16, 19]. To study the toxicity interaction among different PADs in the mixtures, the joint toxicities of all 15 PADs mixtures were also evaluated according to the DA and IA models (Table 4). Based upon the optimal nonlinear Weibull or Logit models for the 22 individual PADs, the 20 points varying between 1% and 99% for the inhibition concentrations of the DA model () and the inhibition concentrations of the IA model () were determined by the APTox computer program.
In Figures 2 and 3, the comparisons of the experimental data with the results calculated using the DA and IA models are illustrated. The DA and IA models could well predict the toxicity effects of all 15 PADs mixtures within 95% confidences of the experimental DRCs. The calculated DRC model parameters () indicated that the difference between the predicted results of these two models was insignificant in the present experiment.
All of the predicted curves are very similar to each other, but the inhibition effects evaluated by DA model were slightly lower than those predicted by IA model in all cases except for the EE2201 mixture (Figures 2 and 3), which might be related to a very small synergistic or antagonistic effect [15, 19]. With regard to the existing quantitative structureactivity relationship (QSAR) method for the prediction of mixture toxicity, the independent action (IA) and the dose addition (DA) are two ideal and prominent reference models of additive behavior for the analysis of joint activity, but the two approaches do not consider the interactions of constituents in a mixture, which can result in complex and substantial changes in the apparent properties of its components, causing synergistic or antagonistic effects [9]. The result of the present work corresponded with the pharmacological presumption that the combined toxicities of the rigorously similaracting chemicals could be well evaluated using the DA model [6, 13]. This finding was well in accordance with the result of a previous experimental study on phenylurea mixtures; that is, the identical evaluations of the DA and IA models showed that the DA model can be used to predict the toxicities of the similaracting chemicals exactly [6]. Assuming the inherent variance of natural periphyton and epipsammon communities, the joint inhibitions of these 22 rigorously congeneric phenylureas having an analogous toxicological action mechanism could be well evaluated by the DA and IA models. It was found that the RMSE_{DA} and RMSE_{IA}, the and were almost identical for both the DA and IA models (Table 4). As shown in Figures 2 and 3, both concepts created nearly the same concentrationresponse curves in all cases, the obtained data points of three EECR and twelve UDCR mixtures scattered around the solid lines of the DA model and the dotted lines of the IA model. Thus, the effects of 15 PADs mixtures could be well predicted using both DA and IA models. For these 15 various PADs mixtures, the values for the toxic effect calculations using IA model were lower than or identical to the values for the toxic effect calculations using DA model (Table 4).
The polar and narcotic PADs show a fairly broad range of toxicities for numerous aquatic organisms that have been comprehensively studied [21–24]. Their toxicity was related to the property of being waterrepellent, weak H bonding acceptor ability, and strong H bonding donor ability [23]. The mixture toxicities of two aniline compounds to Photobacterium phosphoreum could be principally described by the simple summing of the twocomponent toxicities [25]. The quantitative structureactivity relationship could be well used to express the PADs toxicities to Scenedesmus obliquus. Additionally, the toxicity of mixtures could be successfully predicted by the ideal based on structural factors of individual chemicals in twocomponent mixtures [26]. The data of the present study indicated that the toxicity inhibitions of the PADs mixtures were greater than that of the most toxic PAD compound. The DA model can be well utilized to evaluate the toxic effects of multicomponent mixtures of benzene derivatives, regardless of the inhibition levels and the PADs concentration ratios [19]. The PADs toxic evaluation in water systems is no more limited to individual pure chemicals. Both of the DA and IA approaches may be used to evaluate the toxic effects of PADs mixtures.
4. Conclusions
Assessment research of the individual toxicity of the 22 PADs confirmed that these compounds elicited the common effect on the inhibition luminescence of the photobacterium Q67 at concentrations below those that were lethal to the organisms. Based on the nonlinear leastsquares fitting the experimental doseresponse data of the 22 species of PADs, the best doseresponse curves (DRCs) obtained for the 22 PADs to the photobacterium can be accurately expressed using the twoparameter Weibull or Logit model. The toxic inhibitions of multicomponent PADs mixtures can be well evaluated using the DA and IA models based on the concentrationresponse curves of individual components. Both of these models were proved to be suitable methods for the toxic evaluation of PADs mixtures in the aquatic ecosystem, when the models for the evaluation of multicomponent PADs mixtures were used with the EECR or UDCR mixture design.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The paper has greatly benefited from insightful comments by the editor and anonymous reviewer. This research was financially assisted by National Natural Science Foundation of China (NSFC21207024, NSFC51268008) and the Guangxi Provincial Natural Science Foundation (2011GXNSFA018059).
References
 H. Gao, S. H. Zhang, D. Q. Xiong, N. Liu, W. M. Gong, and Q. Wang, “Study on actue toxicities of phenol and aniline to two marine organisms,” Marine Environmental Science, vol. 25, p. 33, 2006. View at: Google Scholar
 K. W. Thomulka, C. G. Abbas, D. A. Young, and J. H. Lange, “Evaluating median effective concentrations of chemicals with bioluminescent bacteria,” Bulletin of Environmental Contamination and Toxicology, vol. 56, no. 3, pp. 446–452, 1996. View at: Publisher Site  Google Scholar
 W. M. Zhou, D. Q. Fu, and Z. G. Sun, “List for priority pollutants in waters,” Environmental Monitoring in China, vol. 6, no. 4, pp. 1–3, 1990. View at: Google Scholar
 R. Altenburger, H. Schmitt, and G. Schüürmann, “Algal toxicity of nitrobenzenes: combined effect analysis as a pharmacological probe for similar modes of interaction,” Environmental Toxicology and Chemistry, vol. 24, no. 2, pp. 324–333, 2005. View at: Publisher Site  Google Scholar
 T. Backhaus, R. Altenburger, W. Boedeker, M. Faust, M. Scholze, and L. H. Grimme, “Predictability of the toxicity of a multiple mixture of dissimilarly acting chemicals to Vibrio fischeri,” Environmental Toxicology and Chemistry, vol. 19, no. 9, pp. 2348–2356, 2000. View at: Publisher Site  Google Scholar
 T. Backhaus, Å. Arrhenius, and H. Blanck, “Toxicity of a mixture of dissimilarly acting substances to natural algal communities: predictive power and limitations of independent action and concentration addition,” Environmental Science and Technology, vol. 38, no. 23, pp. 6363–6370, 2004. View at: Publisher Site  Google Scholar
 M. Scholze, W. Boedeker, M. Faust, T. Backhaus, R. Altenburger, and L. H. Grimme, “A general bestfit method for concentrationresponse curves and the estimation of loweffect concentrations,” Environmental Toxicology and Chemistry, vol. 20, no. 2, pp. 448–457, 2001. View at: Publisher Site  Google Scholar
 R. Altenburger, M. Nendza, and G. Schüürmann, “Mixture toxicity and its modeling by quantitative structureactivity relationships,” Environmental Toxicology and Chemistry, vol. 22, no. 8, pp. 1900–1915, 2003. View at: Publisher Site  Google Scholar
 F. Luan, X. Xu, H. Liu, and M. N. D. S. Cordeiro, “Prediction of the baseline toxicity of nonpolar narcotic chemical mixtures by QSAR approach,” Chemosphere, vol. 90, no. 6, pp. 1980–1986, 2013. View at: Publisher Site  Google Scholar
 M. CruzMonteagudo, E. AncedeGallardo, M. Jorge, and M. N. Dias Soeiro Cordeiro, “Chemoinformatics profiling of ionic liquids—automatic and chemically interpretable cytotoxicity profiling, virtual screening, and cytotoxicophore identification,” Toxicological Sciences, vol. 136, no. 2, pp. 548–565, 2013. View at: Publisher Site  Google Scholar
 M. CruzMonteagudo and M. N. D. S. Cordeiro, “Chemoinformatics profiling of ionic liquids–uncovering structurecytotoxicity relationships with networklike similarity graphs,” Toxicological Sciences, vol. 138, no. 1, pp. 191–204, 2014. View at: Publisher Site  Google Scholar
 M. Faust, R. Altenburger, T. Backhaus et al., “Joint algal toxicity of 16 dissimilarly acting chemicals is predictable by the concept of independent action,” Aquatic Toxicology, vol. 63, no. 1, pp. 43–63, 2003. View at: Publisher Site  Google Scholar
 Å. Arrhenius, F. Grönvall, M. Scholze, T. Backhaus, and H. Blanck, “Predictability of the mixture toxicity of 12 similarly acting congeneric inhibitors of photosystem II in marine periphyton and epipsammon communities,” Aquatic Toxicology, vol. 68, no. 4, pp. 351–367, 2004. View at: Publisher Site  Google Scholar
 M. Junghans, T. Backhaus, M. Faust, M. Scholze, and L. H. Grimme, “Application and validation of approaches for the predictive hazard assessment of realistic pesticide mixtures,” Aquatic Toxicology, vol. 76, no. 2, pp. 93–110, 2006. View at: Publisher Site  Google Scholar
 W. Y. Huang, F. Liu, S. S. Liu, H. L. Ge, and H. H. Chen, “Predicting mixture toxicity of seven phenolic compounds with similar and dissimilar action mechanisms to Vibrio qinghaiensis sp.nov.Q67,” Ecotoxicology and Environmental Safety, vol. 74, no. 6, pp. 1600–1606, 2011. View at: Publisher Site  Google Scholar
 S. S. Liu, X. Q. Song, H. L. Liu, Y. H. Zhang, and J. Zhang, “Combined photobacterium toxicity of herbicide mixtures containing one insecticide,” Chemosphere, vol. 75, no. 3, pp. 381–388, 2009. View at: Publisher Site  Google Scholar
 J. Zhang, S. S. Liu, and H. L. Liu, “Effect of ionic liquid on the toxicity of pesticide to Vibrioqinghaiensis sp.Q67,” Journal of Hazardous Materials, vol. 170, no. 23, pp. 920–927, 2009. View at: Publisher Site  Google Scholar
 J. Zhang, S. S. Liu, R. N. Dou, and H. L. Liu, “Evaluation on the toxicity of ionic liquid mixture with antagonism and synergism to Vibrio qinghaiensis sp.Q67,” Chemosphere, vol. 82, no. 7, pp. 1024–1029, 2011. View at: Publisher Site  Google Scholar
 L. Y. Mo, S. S. Liu, Y. N. Zhu, H. L. Liu, H. Y. Liu, and Z. S. Yi, “Combined toxicity of the mixtures of phenol and aniline derivatives to Vibrio qinghaiensis sp.Q67,” Bulletin of Environmental Contamination and Toxicology, vol. 87, no. 4, pp. 473–479, 2011. View at: Publisher Site  Google Scholar
 S. S. Liu, J. Zhang, Y. H. Zhang, and L. T. Qin, “Aptox: Assessment and prediction on toxicity of chemical mixtures,” Acta Chimica Sinica, vol. 70, no. 14, Article ID 15111517, pp. 1511–1517, 2012. View at: Publisher Site  Google Scholar
 H. Könemann and A. Musch, “Quantitative structureactivity relationships in fish toxicity studies—part 2: the influence of pH on the QSAR of chlorophenols,” Toxicology, vol. 19, no. 3, pp. 223–228, 1981. View at: Publisher Site  Google Scholar
 M. T. D. Cronin and T. W. Schultz, “Structuretoxicity relationships for phenols to Tetrahymena pyriformis,” Chemosphere, vol. 32, no. 8, pp. 1453–1468, 1996. View at: Publisher Site  Google Scholar
 E. U. Ramos, W. H. J. Vaes, P. Mayer, and J. L. M. Hermens, “Algal growth inhibition of Chlorella pyrenoidosa by polar narcotic pollutants: toxic cell concentrations and QSAR modeling,” Aquatic Toxicology, vol. 46, no. 1, pp. 1–10, 1999. View at: Publisher Site  Google Scholar
 G. H. Lu, C. Wang, and X. L. Guo, “Prediction of toxicity of phenols and anilines to algae by quantitative structureactivity relationship,” Biomedical and Environmental Sciences, vol. 21, no. 3, pp. 193–196, 2008. View at: Publisher Site  Google Scholar
 X. Yuan, L. M. Su, and L. L. Song, “Evaluation and prediction of joint toxicity of aniline and its derivatives to photobacterium phosphoreum,” Chemical Research in Chinese Universities, vol. 20, p. 32, 2004. View at: Google Scholar
 C. Wang, G. H. Lu, Z. Y. Tang, and X. L. Guo, “Quantitative structureactivity relationships for joint toxicity of substituted phenols and anilines to Scenedesmus obliquus,” Journal of Environmental Sciences, vol. 20, no. 1, pp. 115–119, 2008. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Lingyun Mo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.