Cytotoxicity of binary mixtures of human pharmaceuticals in a fish cell line: approaches for non-monotonic concentration-response relationships

A Machine-Learning-Driven Pathophysiology-Based New Approach Method for the Dose-Dependent Assessment of Hazardous Chemical Mixtures and Experimental Validations

Environmental chemicals, such as PFAS, exist as mixtures and are frequently encountered at varying concentrations, which can lead to serious health effects, such as cancer. Therefore, understanding the dose-dependent toxicity of chemical mixtures is essential for health risk assessment. However, comprehensive methods to assess toxicity and identify the mechanisms of these harmful mixtures are currently absent. In this study, the dose-dependent toxicity assessments of chemical mixtures are performed in three methodologically distinct phases. In the first phase, we evaluated our machine-learning method (AI-HNN) and pathophysiology method (CPTM) for predicting toxicity. In the second phase, we integrated AI-HNN and CPTM to establish a comprehensive new approach method (NAM) framework called AI-CPTM that is targeted at refining prediction accuracy and providing a comprehensive understanding of toxicity mechanisms. The third phase involved experimental validations of the AI-CPTM predictions. Initially, we developed binary, multiclass classification, and regression models to predict binary, categorical toxicity, and toxic potencies using nearly a thousand experimental mixtures. This empirical dataset was expanded with assumption-based virtual mixtures, compensating for the lack of experimental data and broadening the scope of the dataset. For comparison, we also developed machine-learning models based on RF, Bagging, AdaBoost, SVR, GB, KR, DT, KN, and Consensus methods. The AI-HNN achieved overall accuracies of over 80%, with the AUC exceeding 90%. In the final phase, we demonstrated the superior performance and predictive capability of AI-CPTM, including for PFAS mixtures and their interaction effects, through rigorous literature and statistical validations, along with experimental dose-response zebrafish-embryo toxicity assays. Overall, the AI-CPTM approach significantly improves upon the limitations of standalone AI models, showing extensive enhancements in identifying toxic chemicals and mixtures and their mechanisms. This study is the first to develop a hybrid NAM that integrates AI with a pathophysiology method to comprehensively predict chemical-mixture toxicity, carcinogenicity, and mechanisms.

Metabolomics revealed disruptions in amino acid and antioxidant biochemistry in Daphnia magna exposed to industrial effluents associated with plastic and polymer production

Industrial wastewater effluents are a major source of chemicals in aquatic environments, and many of these chemicals may negatively impact aquatic life. In this study, the crustacean Daphnia magna, a common model organism in ecotoxicity studies, was exposed for 48 h to nine different industrial effluent samples from manufacturing facilities associated with the production of plastics, polymers, and coating products at a range of dilutions: 10, 25, 50, 100% (undiluted). A targeted metabolomic-based approach using liquid chromatography-tandem mass spectrometry (LC-MS/MS) was used to quantify polar metabolites from individual daphnids that survived the 48 h exposure. Multivariate analyses and metabolite changes revealed metabolic perturbations across all effluent samples studied, with non-monotonic responses and both up and downregulation relative to the unexposed control. Pathway analyses indicated the disruption of similar and distinct pathways, mostly connected to protein synthesis, amino acid metabolism, and antioxidant processes. Overall, we observed disruptions in Daphnia biochemistry that were similar across the effluent samples, but with unique features for each effluent sample. Additionally, non-monotonic heightened responses suggested additive and/or synergistic interactions between the chemicals within the industrial effluents. These findings demonstrate that targeted metabolomic approaches are a powerful tool for the biomonitoring of aquatic ecosystems in the context of complex mixtures, such as industrial wastewater effluents.

Combined cytotoxicity of phthalate esters on HepG2 cells: A comprehensive analysis of transcriptomics and metabolomics

Phthalate esters (PAEs), widely used as plasticizers, may pose a potential environmental and human hazard. The aim of this study was to compare the cytotoxicity of di(2-ethylhexyl) phthalates (DEHP) and dibutyl phthalate (DBP)) after their exposure to HepG2 cells alone or in combination. HepG2 cells treated with individual/combined DEHP and DBP at a dose of 10-2 M for 24 h were selected for metabolome and transcriptome analysis. The results demonstrated that exposure to the mixtures of DEHP and DBP caused enhanced or reduced toxic effects regarding 8 pathways with 1065 downregulated genes and 643 upregulated genes, in comparison with those of single chemicals. The combined toxicity of mixture revealed both synergistic and antagonistic interactions between DEHP and DBP. Besides, combined exposure to DEHP and DBP promoted TCA cycle, pyrimidine, and purine metabolism, while an antagonistic effect on fatty acid derangement should require further investigation. To summarize, our results suggest that DEHP exposed alone or combined with DBP caused a variety of metabolic disorders, and the type of combination effects varied among metabolic pathways.

Combined cytotoxicity of ZnO nanoparticles and chlorpyrifos in the rainbow trout, Oncorhynchus mikyss, gonadal cell line RTG-2

The toxicity of ZnO nanoparticles (ZnO NPs) in aquatic organisms has been extensively studied, but little information is available on the effects associated with their interaction with other contaminants. In this context, the in vitro effects of co-exposure of chlorpyrifos (CPF) and ZnO NPs on fish-derived cells were investigated. A selection of concentrations was tested in single and binary exposures: CPF (0.312 - 75 mg/L) and ZnO NPs (10 - 100 mg/L). Cytotoxicity was measured using commonly used cellular endpoints: Alamar Blue/CFDA-AM for viability and plasma membrane integrity, NRU for lysosomal disruption and MTT for mitochondrial function. In addition, specific mechanisms of toxicity for CPF and ZnO NPs were tested: acetylcholinesterase (AChE) activity and ROS generation, respectively. AChE was by far the most sensitive assay for single exposure to CPF. There was no concentration-response relationship for ROS after single exposure to ZnO NPs, but 10 mg/L produced significant effects only for this cellular endpoint. Co-exposure of CPF with 10 m/L of ZnO NPs produced significant effects in almost all endpoints tested, which were enhanced by co-exposure with 100 mg/L of ZnO NPs. AChE testing of additional co-exposures with bulk ZnO, together with the application of the Independent Action (IA) prediction model, which allowed us to draw more in-depth conclusions on the toxicological behavior of the mixture. Synergism was observed at 0.625 mg/L CPF concentration and antagonism at 5 mg/L CPF in mixtures containing 100 mg/L of both ZnO NPs and bulk ZnO. However, more cases of synergism between CPF and ZnO NPs occurred at intermediate CPF concentrations, demonstrating that nano-sized particles have a more toxic interaction with CPF than bulk ZnO. Therefore it can be argued that in vitro assays allow the identification of interaction profiles of NP-containing mixtures by achieving multiple endpoints with a large number of concentration combinations.

Evaluation of Rare Earth Element-Associated Hormetic Effects in Candidate Fertilizers and Livestock Feed Additives

Rare earth elements (REEs) are recognized as emerging contaminants with implications in human and environmental health. Apart from their adverse effects, REEs have been reported as having positive effects when amended to fertilizers and livestock feed additives, thus suggesting a hormetic trend, implying a concentration-related shift from stimulation to inhibition and toxicity, with analogous trends that have been assessed for a number of xenobiotics. In view of optimizing the success of REE mixtures in stimulating crop yield and/or livestock growth or egg production, one should foresee the comparative concentration-related effects of individual REEs (e.g., Ce and La) vs. their mixtures, which may display distinct trends. The results might prompt further explorations on the use of REE mixtures vs. single REEs aimed at optimizing the preparation of fertilizers and feed additives, in view of the potential recognition of their use in agronomy and zootechny.

Statistical modeling of hormesis quantities in inverted U-shaped dose-response relationships by reparameterization of a bilogistic model

Statistical procedures that allow quantitative determination of hormesis features are required for quantitative characterization of hormesis to provide information on the biphasic dose-response phenomenon and its variability. Only a direct estimate of individual effective doses in hormetic dose-response relationships is possible using prior extensions of the bilogistic model of Beckon and coworkers. This study presented further extensions of the model to determine the toxic potency and hormetic dose zone by estimating two effective doses simultaneously. In addition, the extended models allow for partitioning the hormetic dose zone through the dose of maximum stimulation. This study demonstrated a 4-step statistical modeling approach to quantify 20 hormesis quantities. The applicability and challenges of the mathematical procedures are discussed based on a few examples of hormetic dose-response relationships. The syntaxes for the analyses were provided as Appendix to demonstrate its implementation in SAS^® statistical software. Given the variability of hormetic dose-responses generated from toxicological studies in many disciplines, the proposed approach cannot apply to all dose-response patterns. However, we hope the proposed extensions could provide versatile statistical tools for quantitatively exploring a variety of biphasic dose-response curves.

Transcriptomic Points of Departure Calculated from Rainbow Trout Gill, Liver, and Gut Cell Lines Exposed to Methylmercury and Fluoxetine

Ethical and resource limitation concerns are pushing chemicals management to develop alternatives to animal testing strategies. The objective of our study was to determine whether transcriptomic point of departure (tPOD) values could be derived from studies that followed Organisation for Economic Co-operation and Development (OECD) Test No. 249 (rainbow trout gill cell line), as well as from studies on trout liver and gut cells. Gill, liver, and gut cell lines were exposed to methylmercury and fluoxetine. Concentrations causing 50% cytotoxicity (LC50) were derived, the whole transcriptome was sequenced, and gene tPOD and pathway benchmark dose (BMD) values were derived from transcriptomic dose-response analysis. Differences in LC50 and transcriptomic responses across the cell lines were noted. For methylmercury, the tPOD_mode values were 14.5, 20.5, and 17.8 ppb for the gill, liver, and gut cells, respectively. The most sensitive pathway (pathway BMDs in parentheses) was ferroptosis in the gill (3.1 ppb) and liver (3.5 ppb), and glutathione metabolism in the gut (6.6 ppb). For fluoxetine, the tPOD_mode values were 109.4, 108.4, and 97.4 ppb for the gill, liver, and gut cells, respectively. The most sensitive pathway was neurotrophin signaling in the gill (147 ppb) and dopaminergic signaling in the gut (86.3 ppb). For both chemicals, the gene tPOD and pathway BMD values were lower than cytotoxic concentrations in vitro, and within 10-fold below the in vivo LC50s. By bringing together transcriptomics and dose-response analysis with an OECD test method in three cell lines, the results help to establish an in vitro method yielding tPOD values that are hypothesized to be protective of in vivo concentrations associated with adverse outcomes, and also give insights into mechanisms of action. Environ Toxicol Chem 2022;41:1982-1992. © 2022 SETAC.

Assessment of the Potential Ecotoxicological Effects of Pharmaceuticals in the World's Rivers

During their production, use, and disposal, active pharmaceutical ingredients (APIs) are released into aquatic systems. Because they are biologically active molecules, APIs have the potential to adversely affect nontarget organisms. We used the results of a global monitoring study of 61 APIs alongside available ecotoxicological and pharmacological data to assess the potential ecotoxicological effects of APIs in rivers across the world. Approximately 43.5% (461 sites) of the 1052 sampling locations monitored across 104 countries in a recent global study had concentrations of APIs of concern based on apical, nonapical, and mode of action-related endpoints. Approximately 34.1% of the 137 sampling campaigns had at least one location where concentrations were of ecotoxicological concern. Twenty-three APIs occurred at concentrations exceeding "safe" concentrations, including substances from the antidepressant, antimicrobial, antihistamine, β-blocker, anticonvulsant, antihyperglycemic, antimalarial, antifungal, calcium channel blocker, benzodiazepine, painkiller, progestin, and lifestyle compound classes. At the most polluted sites, effects are predicted on different trophic levels and on different endpoint types. Overall, the results show that API pollution is a global problem that is likely negatively affecting the health of the world's rivers. To meet the United Nations' Sustainable Development Goals, work is urgently needed to tackle the problem and bring concentrations down to an acceptable level. Environ Toxicol Chem 2022;41:2008-2020. © 2022 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Occurrence, detection and ecotoxicity studies of selected pharmaceuticals in aqueous ecosystems- a systematic appraisal

Pharmaceutical compounds (PCs) have globally emerged as a significant group of environmental contaminants due to the constant detection of their residues in the environment. The main scope of this review is to fill the void of information on the knowledge on the African occurrence of selected PCs in environmental matrices in comparison with those outside Africa and their respective toxic actions on both aquatic and non-aquatic biota through ecotoxicity bioassays. To achieve this objective, the study focused on commonly used and detected pharmaceutical drugs (residues). Based on the conducted literature survey, Africa has the highest levels of ciprofloxacin, sulfamethoxazole, lamivudine, acetaminophen, and diclofenac while Europe has the lowest of all these PC residues in her physical environments. For ecotoxicity bioassays, the few data available are mostly on individual groups of pharmaceuticals whereas there is sparsely available data on their combined forms.

Embryonic toxicity of 3,4-dichloroaniline (3,4-DCA) on Javanese medaka (Oryzias javanicus Bleeker, 1854)

•

The 96 h LC₅₀ of 3,4-dichloroaniline in Javanese medaka embryo is 32.87 mg.L⁻¹.

•

3,4-DCA lowers heart rate of developing Javanese medaka embryos.

•

The sublethal concentration of 3,4-DCA delays hatching in Javanese medaka embryo.

•

The LOEC for deformities in embryos of Javanese medaka was 0.5 mg.L⁻¹.

Early-life exposure to toxic chemicals causes irreversible morphological and physiological abnormalities that may last for a lifetime. The present study aimed to determine the toxicity effect of 3,4-Dichloroaniline (3,4-DCA) on Javanese medaka (Oryzias javanicus) embryos. Healthy embryos were exposed to various 3,4-DCA concentrations for acute toxicity (5, 10, 25, 50, and 100 mg.L⁻¹) and sublethal toxicity (0.10, 0.50, 1.25, 2.50, and 5.00 mg.L⁻¹) for 96 h and 20 days respectively. Acute toxicity test revealed that the median lethal concentration (96h-LC₅₀) was 32.87 mg.L⁻¹ (95 % CI = 27.90–38.74, R² = 0.95). Sublethal exposure revealed that 1.25 mg.L^-1 at 3 days post-exposure (3 dpe) has a significant lower heartrate (120 ± 12.3 beats/min., p < 0.01), while at 7 dpe those exposed to 5 mg.L⁻¹ (141.8 ± 8.3 beats/min) had significantly (p < 0.01) lower heart rate compared to other treatments. Likewise, at 13 dpe, 5.00 mg.L⁻¹ (110.4 ± 17.3 beats/min) and 2.5 mg.L^-1 (130.4 ± 8.3 beats/min) were significantly lower (p < 0.001) compared to control. None of the embryos in 5.00 mg.L⁻¹ and 2.50 mg.L^-1 treatment groups survived at the end of the experiment. The results indicated a concentration-dependent response. The lowest observed effect concentration (LOEC) that exerted developmental deformities was 0.5 mg.L⁻¹. Javanese medaka embryo have low sensitivity to acute toxicity of 3,4-DCA, but developmental abnormalities at sublethal concentrations were observed.

Cell-based assays as an alternative for the study of aquatic toxicity of pharmaceuticals

An increasing number and amount of pharmaceuticals for human and veterinary use currently reach the aquatic environment, and the determination of their effects on aquatic organisms becomes of major importance. The 96-h fish lethal test is one of the conventional assays required for environmental hazardous assessment, but it is extremely time-consuming and costly, and it raises ethical concerns. In a broad study, we compared the ability of cell-based assays to detect, in absolute terms, lethal toxicity in fish due to pharmaceuticals in order to select sensitive cell lines to be posteriorly used as an alternative to fish testing. This study also explored the sensitivity of the rat cardiomyoblast H9c2(2-1) cell line and the suitability of the sulforhodamine B colorimetric assay regarding 15 pharmaceuticals belonging to 9 different therapeutic classes. The relation between in vivo and in vitro data was expressed as LC_50,96h/EC₅₀ ratios, and 66% of concordant data were attained. Accordingly, it was possible to conclude that cell-based assays could be considered a suitable alternative to fish lethal testing for pharmaceuticals, which, after validation, may dramatically reduce the number of fish required for environmental hazardous assessment. Several cell lines were selected as promising alternatives, but H9c2(2-1), HepG2, PLHC-1, and RTG-2 could be considered suitable starting cell types for further studies, as relevant results were obtained with low exposure times.

Behavior and histopathology as biomarkers for evaluation of the effects of paracetamol and propranolol in the neotropical fish species Phalloceros harpagos

Pharmaceutical drugs in the aquatic environment can induce adverse effects on nontarget organisms. This study aimed to assess the short-term effects of sublethal concentrations of both paracetamol and propranolol on the fish Phalloceros harpagos, specifically light/dark preference, swimming patterns, skin pigmentation, histopathology, and liver glycogen levels. Fish were acutely exposed to sublethal concentrations of both paracetamol (0.008, 0.08, 0.8, 8, 80 mg L^-1) and propranolol (0.0001, 0.001, 0.01, 0.1, 1 mg L^-1) under controlled conditions. For scototaxis, a significant preference for the dark compartment was observed for the group exposed to the highest concentration of paracetamol (80 mg L^-1). Propranolol exposure significantly altered the swimming pattern, especially in fish exposed to the 0.001 mg L^-1 concentration. Pigmentation was reduced in propranolol-exposed fish (0.1, 1 mg L^-1). The lowest concentration of propranolol (0.0001 mg L^-1) induced a decrease of histochemical reaction for hepatic glycogen. These data demonstrate that pharmaceuticals can induce sublethal effects in nontarget organisms, even at low concentrations, compromising specific functions of the individual with ecological relevance, such as energy balance and behavior.

In vitro nuclear receptor inhibition and cytotoxicity of hydraulic fracturing chemicals and their binary mixtures

The widespread use of hydraulic fracturing (HF) in oil and gas extraction operations has led to concern over environmental risks posed by chemicals used in HF fluids. Here we employed a suite of stable luciferase reporter gene assays to investigate the potential for selected HF chemicals or geogenics to activate or antagonise nuclear receptor signalling. We screened three biocides (bronopol [BP], glutaraldehyde [GA], and tetrakis(hydroxymethyl)phosphonium sulfate [THPS]), a surfactant (2-butoxyethanol), a friction reducer (polyacrylamide), and a coal seam geogenic (o-cresol) for their potential to act as agonists or antagonists of the estrogen receptor, androgen receptor, progesterone receptor (PR), glucocorticoid receptor or peroxisome proliferator-activated receptor gamma (PPARγ). None of the chemicals induced luciferase activity in any of assays used in the study. In antagonistic mode, BP, GA and THPS caused reductions in luciferase activity in the reporter assays at higher concentrations (50-100 μM), while at low concentrations (2-10 μM) GA and THPS enhanced luciferase activity in some assays relative to controls. None of the other tested chemicals exhibited antagonism in the selected assays. In most cases, altered receptor signalling only occurred at concentrations exhibiting cytotoxicity. However, PPARγ activity, and to a lesser extent PR activity, were inhibited by THPS at sub-cytotoxic concentrations. The majority of binary combinations tested exhibited significantly less-than-additive cytotoxicity, and none of the combinations exhibited synergistic cytotoxicity. In summary, the results of the present study indicate that the selected chemicals are not likely to function as direct agonists of the nuclear receptors tested, and only one chemical, THPS was an apparent partial antagonist of two nuclear receptors.

Nonlinear responses to waterborne cadmium exposure in zebrafish. An in vivo study

Cadmium (Cd) has proved to be associated with numerous toxic effects in aquatic organisms via waterborne exposure. With a view to investigate Cd toxicity along a broad spectrum of exposures reaching from environmental to toxic, we employed adult zebrafish (Danio rerio) for an in vivo study. A number of 10 fish per tank were placed in 40L tanks and were exposed for 30 days to 0.0, 5.0, 25, 50, 75, 100 and 1000μgCd per liter. There were 2 tanks for each Cd exposure (duplicate experiment). Mortality was recorded daily, dead fish were collected and tissue samples were obtained for histologic observation, whereas remaining tissues were stored for Cd burden determination. Surviving fish were collected at the end of the experiment. Median overall survival (OS) in days was found to be 9.0, 11.0, 8.0 and 7.0 for 25μg/L, 50μg/L, 75μg/L and 100μg/L respectively, with all of them showing mortality greater than 50%. Remarkably, fish exposed to the highest Cd concentration (1000μg/L) survived the longest exhibiting a mean OS of 29.2 days. Cd determination in fish tissue was conducted with an in house ICP-MS method and levels ranged from 3.1 to 29.1ng/mg. Log Cd tissue levels were significantly correlated with the log Cd exposure levels (r = 0.535, p < 0.001). The highest Cd burden was determined for fish exposed to 1000μgCd /L (mean = 12.2ng/mg). Histopathology supported these results. Our findings disclose a deviation in toxic responses through the range of Cd concentrations, leading to nonlinear responses. These differentiated responses, could be linked to hormesis phenomena.

Predicting biphasic responses in binary mixtures: Pelargonic acid versus glyphosate

Predicting hormesis in mixtures is challenging, but essential considering that chemical exposures often occur in mixtures and at low doses. This study investigated mixture effects with two herbicides prone to induce hormesis and to interact, namely pelargonic acid versus glyphosate. Five independent mixture experiments were conducted in vitro to assess effects on root growth of lettuce. Mixture effects on the dose were analyzed using classical joint-action models in terms of deviation from the reference model of concentration addition. For effects on the hormetic magnitude (ymax), a linear reference model was utilized. Hormesis was inconsistent across rays, so that effects on inhibitory doses and ymax could be evaluated, but not effects on hormetic doses. Mixture effects on the dose were additive at lower doses changing to strong high-dose synergism. Mixture effects on ymax followed a linear change with mixture ratio or significantly deviated from linearity with a one-sided trend across rays in two experiments. The trend was antipodal between experiments, but well described by a curved ymax model based on single dose-response relationships. Atypical ymax deviations were associated with strong synergism at ED50, suggesting that the linearity model applies for chemicals showing no/minor interaction at ED50, while for strongly interacting chemicals ymax predictions seem more critical. The study unambiguously proved synergism on the dose for pelargonic acid versus glyphosate and indicated an impact of these joint effects on ymax. The study confirms the predictability of hormesis in mixtures and provides a further methodological step towards an incorporation of hormesis into mixture-toxicity evaluations.

An automated fitting procedure and software for dose-response curves with multiphasic features

In cancer pharmacology (and many other areas), most dose-response curves are satisfactorily described by a classical Hill equation (i.e. 4 parameters logistical). Nevertheless, there are instances where the marked presence of more than one point of inflection, or the presence of combined agonist and antagonist effects, prevents straight-forward modelling of the data via a standard Hill equation. Here we propose a modified model and automated fitting procedure to describe dose-response curves with multiphasic features. The resulting general model enables interpreting each phase of the dose-response as an independent dose-dependent process. We developed an algorithm which automatically generates and ranks dose-response models with varying degrees of multiphasic features. The algorithm was implemented in new freely available Dr Fit software (sourceforge.net/projects/drfit/). We show how our approach is successful in describing dose-response curves with multiphasic features. Additionally, we analysed a large cancer cell viability screen involving 11650 dose-response curves. Based on our algorithm, we found that 28% of cases were better described by a multiphasic model than by the Hill model. We thus provide a robust approach to fit dose-response curves with various degrees of complexity, which, together with the provided software implementation, should enable a wide audience to easily process their own data.

Genotoxic and immunotoxic potential effects of selected psychotropic drugs and antibiotics on blue mussel (Mytilus edulis) hemocytes

The potential toxicity of pharmaceuticals towards aquatic invertebrates is still poorly understood and sometimes controversial. This study aims to document the in vitro genotoxicity and immunotoxicity of psychotropic drugs and antibiotics on Mytilus edulis. Mussel hemocytes were exposed to fluoxetine, paroxetine, venlafaxine, carbamazepine, sulfamethoxazole, trimethoprim and erythromycin, at concentrations ranging from μg/L to mg/L. Paroxetine at 1.5 μg/L led to DNA damage while the same concentration of venlafaxine caused immunomodulation. Fluoxetine exposure resulted in genotoxicity, immunotoxicity and cytotoxicity. In the case of antibiotics, trimethoprim was genotoxic at 200 μg/L and immunotoxic at 20 mg/L whereas erythromycin elicited same detrimental effects at higher concentrations. DNA metabolism seems to be a highly sensitive target for psychotropic drugs and antibiotics. Furthermore, these compounds affect the immune system of bivalves, with varying intensity. This attests the relevance of these endpoints to assess the toxic mode of action of pharmaceuticals in the aquatic environment.

Toxicological effects of clofibric acid and diclofenac on plasma thyroid hormones of an Indian major carp, Cirrhinus mrigala during short and long-term exposures

In the present investigation, the toxicity of most commonly detected pharmaceuticals in the aquatic environment namely clofibric acid (CA) and diclofenac (DCF) was investigated in an Indian major carp Cirrhinus mrigala. Fingerlings of C. mrigala were exposed to different concentrations (1, 10 and 100μgL(-1)) of CA and DCF for a period of 96h (short term) and 35 days (long term). The toxic effects of CA and DCF on thyroid hormones (THs) such as thyroid stimulating hormone (TSH), thyroxine (T4) and triiodothyronine (T3) levels were evaluated. During the short and long-term exposure period TSH level was found to be decreased at all concentrations of CA (except at the end of 14(th) day in 1 and 10μgL(-l) and 21(st) day in 1μgL(-l)) whereas in DCF exposed fish TSH level was found to be increased when compared to control groups. T4 level was found to be decreased at 1 and 100μgL(-l) of CA exposure at the end of 96h. However, T4 level was decreased at all concentrations of CA and DCF during long-term (35 days) exposure period. Fish exposed to all concentrations of CA and DCF had lower level of T3 in both the treatments. These results suggest that both CA and DCF drugs induced significant changes (P<0.01 and P<0.05) on thyroid hormonal levels of C. mrigala. The alterations of these hormonal levels can be used as potential biomarkers in monitoring of pharmaceutical drugs in aquatic organisms.