BioMed Research International
Volume 2016 (2016), Article ID 4158451, 6 pages
http://dx.doi.org/10.1155/2016/4158451

Research Article

Quantitative Characterization of the Toxicities of Cd-Ni and Cd-Cr Binary Mixtures Using Combination Index Method

Lingyun Mo,1,2 Minyi Zheng,1,2 Meng Qin,1,2 Xin Zhang,1,2 Jie Liu,1,2 Litang Qin,1,2 Honghu Zeng,1,2 and Yanpeng Liang1,2

¹Guangxi Key Laboratory of Environmental Pollution Control Theory and Technology, Guilin University of Technology, Guilin 541004, China
²Collaborative Innovation Center for Water Pollution Control and Water Safety in Karst Area, Guilin University of Technology, Guilin 541004, China

Received 24 July 2016; Revised 23 October 2016; Accepted 1 November 2016

Academic Editor: Zongming Ren

Copyright © 2016 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.

Abstract.

Direct equipartition ray design was used to construct Cd-Ni and Cd-Cr binary mixtures. Microplate toxicity analysis was used to evaluate the toxicity of individual substance and the Cd-Ni and Cd-Cr mixtures on Chlorella pyrenoidosa and Selenastrum capricornutum. The interacting toxicity of the mixture was analyzed with concentration addition (CA) model. In addition, combination index method (CI) was proposed and used to quantitatively characterize the toxicity of the binary mixtures of Cd-Ni and Cd-Cr observed in experiment and find the degree of deviation from the predicted outcome of the CA model, that is, the intensity of interacting toxicity. Results indicate that most of the 20 binary mixtures exhibit enhancing and synergistic effect, and only Cd-Cr-R4 and Cd-Cr-R5 mixtures have relatively high antagonistic effects against C. pyrenoidosa. Based on confidence interval, CI can compare the intensities of interaction of the mixtures under varying levels of effect. The characterization methods are applicable for analyzing binary mixture with complex interaction.

1. Introduction

The evaluation and prediction of combined effect of chemical pollutants in the environment have become an interesting research topic in environmental chemistry [1, 2]. Several research on mixtures indicate that the dependency of concentration ratio and effect level with interaction is a common phenomenon [3–7] that requires a reasonable and easily applicable method to represent the toxic interaction of the mixture so as to ensure the reliability of the evaluation results. At present, two additional reference models, namely, concentration addition (CA) and independence action (IA), are widely used to evaluate mixture toxicity. CA model is applicable for mixtures with similar mode of action, whereas the IA model is applicable for mixtures with different modes of action [1, 8–11]. However, the modes of action of many substances are still unknown, and selecting different standards is still difficult, so CA or IA cannot be selected [12, 13]. Through theoretical research and experimental verification, Chou [14–16] proposed to evaluate the toxicity interaction of mixtures using the combination index (CI) based on 50% effective equation without depending on the mode of action. CI has been widely used in assessing the toxicity of drug combinations [17–19]. In recent years, CI has attracted the attention of scholars of environmental sciences [20–22].

All toxicity tests contain experimental error. Different toxicants or pollutants have different concentration-response features, and concentration-response curve (CRC) fitting also contains error of fitting. All these affect the uncertainty of CI index. This research considers the influence of experimental error/different concentration-effect feature and error of fitting on CI. The toxicity interaction is also analyzed by taking the binary mixture of Cd-Ni and Cd-Cr as examples.

2. Materials and Methods

2.1. Type of Algae, Key Instruments, and Drugs

Chlorella pyrenoidosa and Selenastrum capricornutum with serial number FACHB-5 and FACHB-271 are purchased from the FACHB-collection of typical culture collection committee of Chinese Academy of Science. BG11 culture medium, which is placed into manual climatic box, is used for cultivation at 22°C, illumination intensity of 3000 lux, and light-dark cycle of 12 h : 12 h. The algae were collected during logarithmic phase for experiment.

Other materials used were as follows: TECAN infinite 200 microplate reader (Switzerland), illumination incubator with constant temperature and artificial climate (Jintan Yitong Electronics Co. Ltd.), 96-hole transparent polystyrene microplate with flat base (Greiner), YXQ-LS-70A vertical pressure steam sterilization, SW-CF-IFD clean bench, Milli-Q ultrapure water system, Haier drug freezer, and Eppendorf pipette.

CdCl₂·2.5H₂O[Cd] (Damas-Beta), Ni(NO₃)₂·6H₂O[Ni] (Guangdong Guanghua Technology Stock Co. Ltd.), and K₂Cr₂O₇[Cr] (Damas-Beta) were also used. The prepared solution of the compounds was stored in a refrigerator at 4°C to prepare for the test. See the basic properties of heavy metal compounds in Table 1.

Table 1: Physical properties of three heavy metals.

Compound
CAS
Molecular weight
Purity (%)

CdCl₂⋅2.5H₂O
7790-78-5
228.36
≥98

Ni(NO₃)₂⋅6H₂O
13478-00-7
290.79
≥98

K₂Cr₂O₇
7778-50-9
294.18
≥99

2.2. Experimental Design and Toxicity Test

To explore the toxicity interaction of Cd-Ni and Cd-Cr mixtures within the whole space range of concentration, the concentration design of the binary mixtures was created through direct equipartition ray design [4]. The points of 50% effective concentration (EC₅₀) of each component on the 2D coordinate planes of concentration comprising two components of the mixture were selected. The two points of EC₅₀ were then connected, and the segment was divided into six equal parts. Five radials were made (R1, R2, R3, R4, and R5) from the origin through the point of division. Appropriate dilution factors were selected and 12 points of concentration on each radial were designed. The concentration of each radial can be seen in Tables 2 and 3.

Table 2: Concentration ratio of Cd-Ni and EC₅₀ value of 10 binary mixture rays.

Binary mixture
rays

EC₅₀

Cd
Ni

Cd-Ni toxicity to C. pyrenoidosa
R1
0.1034
0.8966

R2
0.2239
0.7761

R3
0.3658
0.6342

R4
0.5357
0.4643

R5
0.7425
0.2575

Cd-Ni toxicity to S. capricornutum
R1
0.1009
0.8991

R2
0.2192
0.7808

R3
0.3595
0.6405

R4
0.5289
0.4711

R5
0.7373
0.2627

Table 3: Concentration ratio of Cd-Cr and EC₅₀ value of 10 binary mixture rays.

Binary mixture
Rays

EC₅₀

Cd
Cr

Cd-Cr toxicity to C. pyrenoidosa
R1
0.9341
0.0659

R2
0.8500
0.1500

R3
0.7392
0.2608

R4
0.5863
0.4137

R5
0.3618
0.6382

Cd-Cr toxicity to S. capricornutum
R1
0.0645
0.9355

R2
0.1470
0.8530

R3
0.2563
0.7437

R4
0.4081
0.5919

R5
0.6328
0.3672

The toxic effects of individual heavy metal and the Cd-Ni and Cd-Cr mixtures on C. pyrenoidosa and S. capricornutum were evaluated through microplate toxicity analysis (MTA) [23, 24]. Water was added around the 96-hole microplate. On the plate, all the 60 holes on the 6 × 10 array in the middle are used as test holes in toxicity test of algae. There are 24 holes in columns 2, 6, 7, and 11. A total of 100 µL ultrapure water was added in each hole to act as the blank control group. In the remaining 36 holes, we performed three parallel toxicity tests on 12 pollutants with different concentrations designed based on geometric progression, with total volume of the test solution at 100 μL. Finally, 100 μL algae solution was added at logarithmic phase in each of the 60 test holes for toxicity test to come up with a total volume of 200 μL.

Microplate reader tests the absorbancy of each hole on the microplate, and wave length is set as 682 nm. The inhibition of the pollutant on C. pyrenoidosa at different exposure duration was calculated using the following equation:

where indicates the inhibition of individual compound or mixture with certain concentration on C. pyrenoidosa or S. capricornutum at different exposure duration; refers to the point of exposure duration of the pollutant, , and 4 corresponding to the exposure duration of 0 h, 24 h, 36 h, 72 h, and 96 h, respectively; is the value of algae solution in pollutant treatment group at ; is the value of algae solution in blank control group at ; is the value of algae solution in pollutant treatment group at ; and is the value of algae solution in blank control group at .

2.3. Analysis on Toxicity Interaction of Mixture

The synergism, additive effect, or antagonism quantification of binary mixtures was simulated for CI values using CompuSyn [12, 14–16, 25]. The general equation for -drug combination at inhibition is shown as follows [12]:

where is the combination index for drugs at inhibition, CI < 1, CI = 1, and CI > 1 indicate synergism, additive effect, and antagonism, respectively, is the doses of drug alone that inhibit , is the portion of drug in combination that also inhibits , is the sum of the dose of drugs that exerts inhibition in combination, is the proportionality of the dose of each of drugs that exerts inhibition in combination, and is the dose of each drug alone that exerts inhibition, where is the median-effect dose (antilog of the -intercept of the median-effect plot), is the fractional inhibition at inhibition, and is the slope of the median-effect plot, which depicts the shape of the dose-effect curve.

3. Results and Discussion

3.1. Toxicity of Individual Substance on C. pyrenoidosa

Through MTA, the concentration-response data points (Figures 1 and 2) of Cd, Ni, and Cr on C. pyrenoidosa and S. capricornutum were tested. The median-effect dose (effect concentration, EC₅₀) that inhibits the pyrenoidosa and S. capricornutum under study by 50% was calculated by CompuSyn software [26]. The (EC₅₀) value was calculated by the following:

Figure 1: Concentration-response curves of three heavy metals on C. pyrenoidosa (CP).

Figure 2: Concentration-response curves of three heavy metals on S. capricornutum (SP).

The order of toxicities of the three heavy metals on C. pyrenoidosa is Cr (EC₅₀ = ) > Cd (EC₅₀ = ) > Ni (EC₅₀ = ), and on S. capricornutum, the order of toxicities is Cd (EC₅₀ = ) > Ni (EC₅₀ = ) > Cr (EC₅₀ = ). By comparing the toxicities of three heavy metals on C. pyrenoidosa and S. capricornutum, different green algae were observed to have different tolerance against heavy metals. The tolerance from large to small is Ni (C. pyrenoidosa > S. capricornutum), Cd (C. pyrenoidosa > S. capricornutum), and Cr (S. capricornutum > C. pyrenoidosa). The order of tolerance of S. capricornutum and C. pyrenoidosa against Cr is consistent with the research results of Chen et al. [27].

3.2. Toxicity Analysis of Binary Mixture of Heavy Metal on C. pyrenoidosa and S. capricornutum

3.2.1. Mixture Characterization

The CI values of all 20 mixture rays in the binary mixture system for two groups (Cd-Ni and Cd-Cr) were calculated using formula (2). Figures 3 and 4 present the change curve of CI with the effect level (). The concentration ratio of individual in mixtures and EC₅₀ for two groups mixtures were listed in Tables 2 and 3.

(a)

(b)

Figure 3: Plot of the combination index versus the effect expressed by Cd-Ni mixed system of binary mixtures. ((a) C. pyrenoidosa; (b) S. capricornutum).

(a)

(b)

Figure 4: Plot of the combination index versus the effect expressed by Cd-Cr mixed system of binary mixtures. ((a) C. pyrenoidosa; (b) S. capricornutum).

3.2.2. Cd-Ni Mixture System

Figure 3 indicates the CI of Cd-Ni mixture system on C. pyrenoidosa (CP) and S. capricornutum (SP). As can be seen from Figure 3(a), the toxicity of Cd-Ni on C. pyrenoidosa in low-effect area (lower than 10%) presents an antagonism and synergistic effect in high-effect area (higher than 20%). Mixtures R1–R5 present a synergistic effect in the area with inhibition greater than 30%.

As shown in Figure 3(b), the toxicity of Cd-Ni on S. capricornutum in the whole-effect area presents synergism, additive effect, and antagonism. Mixtures R1–R4 present synergistic effect in the area with inhibition greater than 30%, and mixture R5 presents a synergistic effect in the area with inhibition greater than 40%.

In general, in R1–R5 of Cd-Ni mixture, the concentration proportion of Ni decreases with the increase of the concentration proportion of Cd, and the area of synergistic effect had an increasing trend.

3.2.3. Cd-Cr Mixture System

Figure 4 indicates the relationship between the CI and inhibition of Cd-Cr mixture system on C. pyrenoidosa and S. capricornutum for 96 h.

As seen in Figure 4(a), in Cd-Cr mixture system of R1 toxicity to C. pyrenoidosa presents synergistic effect in the high-effect area (higher than 35%), while the effect lower than 35% presents antaganism. Mixture R2 presents synergistic effect in the effect higher than 10% and mixtures R3 and R4 show synergistic effect for the whole effect. Mixture R5 present synergistic effect in the effect lower than 70% and presents antaganism in high-effect area (higher than 70%).

As seen in Figure 4(b), in Cd-Cr mixtures of R1–R5 toxicity to S. capricornutum presents synergistic effect in the whole effect. For the toxicity of Cd-Cr system on S. capricornutum, Cr concentration gradually increased as the Cd concentration gradually decreased, and the synergistic effect presents in the whole-effect area.

3.2.4. Characterization Results

As seen from Figures 3 and 4, without considering confidence interval, the CI values of radials of many mixtures all deviated from 1 in certain range of effect and interaction exists. For example, the effect of Cd-Ni (R3) on C. pyrenoidosa upon 96 h exposure is , that is, synergistic effect, and the effect of Cd-Cr (R2) on C. pyrenoidosa 96 h exposure is , that is, antagonism. However, accidental error exists in the toxicity test and error of fitting in nonlinear fitting, so the calculated CI has some degree of uncertainty (indicated by 95% confidence interval in this research).

The results above indicate that, for most mixtures, the same mixture presents synergistic effect or addition action and antagonistic effect on C. pyrenoidosa or S. capricornutum in the whole effect area. In addition, all the interactions between the five rays of the same mixture system have different changes, mainly because of the difference in mixing ratio. Hence, the concentration proportion of the mixture has greater influence on the type of interaction of the mixtures.

4. Conclusions

(1)The toxicity interaction of Cd-Ni and Cd-Cr mixtures on C. pyrenoidosa or S. capricornutum mainly contains addition action and synergistic effect; however, the toxicity effect of Cd-Cr on C. pyrenoidosa is addition action in low-effect area. In high-effect area, the effect is antagonistic at effect more than 60%. The concentration range containing antagonism and the intensity of antagonism are affected by several factors, including the components of the mixture, concentration ratio of the components, and the effect level.(2)CI can quantitatively characterize the toxicity effect of the binary mixture from the angle of effect. CI can also compare the intensities of the interaction between the mixtures under any effect level based on confidence interval. The characterization method is applicable for analyzing the binary mixture with complicated interaction.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

The authors are especially grateful for the financial support from the National Natural Science Foundation of China (21667013 and 21407032), the Provincial Natural Science Foundation of Guangxi (2014GXNSFBA118233), Science Research and Technology Development Project of Guangxi (Guikehe1599005-2-2), the Project of High Level Innovation Team and Outstanding Scholar in Guangxi Colleges and Universities (002401013001), and Guangxi Scientific Experiment Center of Mining, Metallurgy, and Environment.

References

D. J. Spurgeon, O. A. H. Jones, J.-L. C. M. Dorne, C. Svendsen, S. Swain, and S. R. Stürzenbaum, “Systems toxicology approaches for understanding the joint effects of environmental chemical mixtures,” Science of the Total Environment, vol. 408, no. 18, pp. 3725–3734, 2010.

T. Backhaus and M. Faust, “Predictive environmental risk assessment of chemical mixtures: a conceptual framework,” Environmental Science and Technology, vol. 46, no. 5, pp. 2564–2573, 2012.

A. Koutsaftis and I. Aoyama, “Toxicity of four antifouling biocides and their mixtures on the brine shrimp Artemia salina,” Science of the Total Environment, vol. 387, no. 1–3, pp. 166–174, 2007.

R.-N. Dou, S.-S. Liu, L.-Y. Mo, H.-L. Liu, and F.-C. Deng, “A novel direct equipartition ray design (EquRay) procedure for toxicity interaction between ionic liquid and dichlorvos,” Environmental Science and Pollution Research, vol. 18, no. 5, pp. 734–742, 2011.

J. Zhang, S.-S. Liu, R.-N. Dou, H.-L. Liu, and J. Zhang, “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.

J. Zhang, S.-S. Liu, Z.-Y. Yu, and J. Zhang, “Time-dependent hormetic effects of 1-alkyl-3-methylimidazolium bromide on Vibrio qinghaiensis sp.-Q67: luminescence, redox reactants and antioxidases,” Chemosphere, vol. 91, no. 4, pp. 462–467, 2013.

J. Zhang and S. S. Liu, “Toxicity interaction between ionic liquids and organophosphate pesticides,” Asian Journal of Ecotoxicolog, vol. 7, no. 4, pp. 408–414, 2012 (Chinese).

R. Altenburger, T. Backhaus, W. Boedeker, M. Faust, M. Scholze, and L. H. Grimme, “Predictability of the toxicity of multiple chemical mixtures to Vibrio fischeri: mixtures composed of similarly acting chemicals,” Environmental Toxicology and Chemistry, vol. 19, no. 9, pp. 2341–2347, 2000.

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.

M. Faust, R. Altenburger, T. Backhaus et al., “Predicting the joint algal toxicity of multi-component s-triazine mixtures at low-effect concentrations of individual toxicants,” Aquatic Toxicology, vol. 56, no. 1, pp. 13–32, 2001.

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.

T.-C. Chou, “Theoretical basis, experimental design, and computerized simulation of synergism and antagonism in drug combination studies,” Pharmacological Reviews, vol. 58, no. 3, pp. 621–681, 2006.

M. Cleuvers, “Aquatic ecotoxicity of pharmaceuticals including the assessment of combination effects,” Toxicology Letters, vol. 142, no. 3, pp. 185–194, 2003.

T.-C. Chou, “The mass-action law based algorithms for quantitative econo-green bio-research,” Integrative Biology, vol. 3, no. 5, pp. 548–559, 2011.

T.-C. Chou, “Drug combination studies and their synergy quantification using the chou-talalay method,” Cancer Research, vol. 70, no. 2, pp. 440–446, 2010.

T.-C. Chou, “Preclinical versus clinical drug combination studies,” Leukemia & Lymphoma, vol. 49, no. 11, pp. 2059–2080, 2008.

H. Kostková, T. Etrych, B. Říhová et al., “HPMA copolymer conjugates of DOX and mitomycin C for combination therapy: physicochemical characterization, cytotoxic effects, combination index analysis, and anti-tumor efficacy,” Macromolecular Bioscience, vol. 13, no. 12, pp. 1648–1660, 2013.

F.-Y. Liu, Z.-J. Zhao, P. Li, X. Ding, and C.-F. Sun, “Effect of rapamycin combined with cisplatin on head and neck squamous cancer cells regulated by CCL19,” Chinese Journal of Stomatology, vol. 46, no. 4, pp. 197–200, 2011.

Y. Yamamoto, I. Hyodo, M. Takigahira et al., “Effect of combined treatment with the epirubicin-incorporating micelles (NC-6300) and 1,2-diaminocyclohexane platinum (II)-incorporating micelles (NC-4016) on a human gastric cancer model,” International Journal of Cancer, vol. 135, no. 1, pp. 214–223, 2014.

I. Rodea-Palomares, A. L. Petre, K. Boltes et al., “Application of the combination index (CI)-isobologram equation to study the toxicological interactions of lipid regulators in two aquatic bioluminescent organisms,” Water Research, vol. 44, no. 2, pp. 427–438, 2010.

I. González-Sánchez, A. Lira-Rocha, A. Navarrete et al., “Synergistic anticancer activity of thiazolo[5,4-b]quinoline derivative D3CLP in combination with cisplatin in human cervical cancer cells,” Anticancer Research, vol. 32, no. 12, pp. 5159–5165, 2012.

I. Rodea-Palomares, F. Leganés, R. Rosal, and F. Fernández-Piñas, “Toxicological interactions of perfluorooctane sulfonic acid (PFOS) and perfluorooctanoic acid (PFOA) with selected pollutants,” Journal of Hazardous Materials, vol. 201-202, pp. 209–218, 2012.

H. L. Ge, S. S. Liu, X. W. Zhu, and L. J. Wang, “Microplate algae toxicity of 9 pesticides to Scenedesmus obliquue based on absorbance,” Asian Journal of Ecotoxicolog, vol. 3, no. 6, pp. 606–612, 2008 (Chinese).

Y.-H. Zhang, S.-S. Liu, X.-Q. Song, and H.-L. Ge, “Prediction for the mixture toxicity of six organophosphorus pesticides to the luminescent bacterium Q67,” Ecotoxicology and Environmental Safety, vol. 71, no. 3, pp. 880–888, 2008.

T.-C. Chou and P. Talalay, “Quantitative analysis of dose-effect relationships: the combined effects of multiple drugs or enzyme inhibitors,” Advances in Enzyme Regulation, vol. 22, pp. 27–55, 1984.

T. C. Chou and N. Martin, CompuSyn for Drug Combinations: PC Software and User's Guide: A Computer Program for Quantitation of Synergism and Antagonism in Drug Combinations, and the Determination of IC50 and ED50 and LD50 Values, ComboSyn Inc, Paramus, NJ, USA, 2005.

H. Chen, G. Pan, H. Yan, and Y. Qin, “Toxic effects of hexavalent chromium on the growth of blue-green microalgae,” Huanjing Kexue/Environmental Science, vol. 24, no. 2, pp. 13–18, 2003 (Chinese).