Screening the critical protein subnetwork to delineate potential mechanisms and protective agents associated with arsenic-induced cutaneous squamous cell carcinoma: A toxicogenomic study

Recent studies show that complex mechanisms are involved in arsenic-induced malignant transformation of cells. This study aimed to decipher molecular mechanisms associated with arsenic-induced cutaneous squamous cell carcinoma (cSCC) and suggest potential protective factors. RNA-seq-based differentially expressed genes between arsenic-exposed human keratinocytes (HaCaT) and controls were used to construct a protein-protein interaction (PPI) network and discover critical subnetwork-based mechanisms. Protective compounds against arsenic toxicity were determined and their target interactions in the core sub-network were identified by the comparative toxicogenomic database (CTD). The binding affinity between the effective factor and target was calculated by molecular docking. A total of 15 key proteins were screened out as critical arsenic-responsive subnetwork (FN1, IL-1A, CCN2, PECAM1, FGF5, EDN1, FGF1, PXDN, DNAJB9, XBP1, ERN1, PDIA4, DNAJB11, FOS, PDIA6) and 7 effective protective agents were identified (folic acid, quercetin, zinc, acetylcysteine, methionine, catechin, selenium). The GeneMANIA predicted detailed interactions of the subnetwork and revealed terms related to unfolded protein response as the main processes. FN1, IL1A and CCN2, as top significant genes, had good docking affinity with folic acid and quercetin, as selected key compounds. Integration of gene expression and protein-protein interaction related to arsenic exposure in cSCC explored the potential mechanisms and protective agents.

Sublethal effects of DBE-DBCH diastereomers on physiology, behavior, and gene expression of Daphnia magna

1,2-dibromo-4-(1,2-dibromoethyl)-cyclohexane (DBE-DBCH) is a brominated flame retardant used in commercial and industrial applications. The use of DBE-DBCH containing products has resulted in an increased release into the environment. However, limited information is available on the long-term effects of DBE-DBCH and its effects in aquatic invertebrates. Thus, the present study was aimed at determining how DBE-DBCH diastereomers (αβ and γδ) affects aquatic invertebrates using Daphnia magna as a model organism. Survival, reproduction, feeding, swimming behavior and toxicogenomic responses to environmental relevant concentrations of DBE-DBCH were analyzed. Chronic exposure to DBE-DBCH resulted in decreased lifespan, and reduced fecundity. Expression of genes involved in reproductive processes, vtg1 and jhe, were also inhibited. DBE-DBCH also induced hypoxia by inhibiting the transcription of genes involved in heme biosynthesis and oxygen transport. Furthermore, DBE-DBCH also inhibited feeding resulting in emptiness of the alimentary canal. Increased expression of the stress response biomarkers was observed following DBE-DBCH exposure. In addition, DBE-DBCH diastereomers also altered the swimming behavior of Daphnia magna. The present study demonstrates that DBE-DBCH cause multiple deleterious effects on Daphnia magna, including effects on reproduction and hormonal systems. These endocrine disrupting effects are in agreement with effects observed on vertebrates. Furthermore, as is the case in vertebrates, DBE-DBCH γδ exerted stronger effects than DBE-DBCH αβ on Daphnia magna. This indicate that DBE-DBCH γδ has properties making it more toxic to all so far studied animals than DBE-DBCH αβ.

Deregulation of whole-transcriptome gene expression in zebrafish (Danio rerio) after chronic exposure to low doses of imatinib mesylate in a complete life cycle study

Imatinib mesylate (IM) is an anticancer drug that belongs to tyrosine kinase inhibitors. We report the results of the first investigation of the chronic exposure of zebrafish (Danio rerio) to IM. The exposure to IM (0.01, 1 and 100 μg/L) was initiated in adult fish and continued through hatching and the offspring generation for seven months. In addition to standard toxicological endpoints, induction of genotoxic effects and whole-genome transcriptome of liver samples of offspring generation of zebrafish were analysed. Exposure to IM did not affect the survival and growth of zebrafish, did not cause any histopathological changes, but it induced a marginal increase in the chromosomal damage in blood cells. The whole-genome transcriptome analyses demonstrated dose-dependent increase in the number of differentially expressed genes with a significantly higher number of deregulated genes in female fish compared to male. Differentially expressed genes included genes involved in response to DNA damage, cell cycle control and regulation of circadian rhythm. Based on the low genotoxic activity and the pattern of the changes in DNA damage responsive genes we consider that at current environmental exposure levels, IM represents low risk for genotoxic effects in aquatic organisms. Exposure to IM also induced deregulation of the expression of genes associated with steroidogenesis and hormone metabolism and function, which indicates hormone-disrupting activity of IM that has not been studied so far. The study provide new information on the potential consequences of chronic exposure to the residues of tyrosine kinase inhibitors, which remain to be further explored.

De novo transcriptome assembly and differential gene expression analysis of the calanoid copepod Acartia tonsa exposed to nickel nanoparticles

The calanoid copepod Acartia tonsa is a reference species in standardized ecotoxicology bioassay. Despite this interest, there is a lack of knowledge on molecular responses of A. tonsa to contaminants. We generated a de novo assembled transcriptome of A. tonsa exposed 4 days to 8.5 and 17 mg/L nickel nanoparticles (NiNPs), which have been shown to reduce egg hatching success and larval survival but had no effects on the adults. Aims of our study were to 1) improve the knowledge on the molecular responses of A. tonsa copepod and 2) increase the genomic resources of this copepod for further identification of potential biomarkers of NP exposure. The de novo assembled transcriptome of A. tonsa consisted of 53,619 unigenes, which were further annotated to nr, GO, KOG and KEGG databases. In particular, most unigenes were assigned to Metabolic and Cellular processes (34-45%) GO terms, and to Human disease (28%) and Organismal systems (23%) KEGG categories. Comparison among treatments showed that 373 unigenes were differentially expressed in A. tonsa exposed to NiNPs at 8.5 and 17 mg/L, with respect to control. Most of these genes were downregulated and took part in ribosome biogenesis, translation and protein turnover, thus suggesting that NiNPs could affect the copepod ribosome synthesis machinery and functioning. Overall, our study highlights the potential of toxicogenomic approach in gaining more mechanistic and functional information about the mode of action of emerging compounds on marine organisms, for biomarker discovering in crustaceans.

Toxicogenomic and bioinformatics platforms to identify key molecular mechanisms of a curcumin-analogue DM-1 toxicity in melanoma cells

Melanoma is a highly invasive and metastatic cancer with high mortality rates and chemoresistance. Around 50% of melanomas are driven by activating mutations in BRAF that has led to the development of potent anti-BRAF inhibitors. However resistance to anti-BRAF therapy usually develops within a few months and consequently there is a need to identify alternative therapies that will bypass BRAF inhibitor resistance. The curcumin analogue DM-1 (sodium 4-[5-(4-hydroxy-3-methoxy-phenyl)-3-oxo-penta-1,4-dienyl]-2-methoxy-phenolate) has substantial anti-tumor activity in melanoma, but its mechanism of action remains unclear. Here we use a synthetic lethal genetic screen in Saccharomyces cerevisiae to identify 211 genes implicated in sensitivity to DM-1 toxicity. From these 211 genes, 74 had close human orthologues implicated in oxidative phosphorylation, insulin signaling and iron and RNA metabolism. Further analysis identified 7 target genes (ADK, ATP6V0B, PEMT, TOP1, ZFP36, ZFP36L1, ZFP36L2) with differential expression during melanoma progression implicated in regulation of tumor progression, cell differentiation, and epithelial-mesenchymal transition. Of these TOP1 and ADK were regulated by DM-1 in treatment-naïve and vemurafenib-resistant melanoma cells respectively. These data reveal that the anticancer effect of curcumin analogues is likely to be mediated via multiple targets and identify several genes that represent candidates for combinatorial targeting in melanoma.

Linking genomic responses of gonads with reproductive impairment in marine medaka (Oryzias melastigma) exposed chronically to the chemopreventive and antifouling agent, 3,3'-diindolylmethane (DIM)

3,3'-Diindolylmethane (DIM) has been promoted as an effective chemopreventive and antifouling additive. However, the concurrent risks or side effects of DIM are not fully understood, especially on tissues responsive to estrogen. Therefore, this study employed marine medaka (Oryzias melastigma) as a test model to evaluate relative safety and explore mechanisms of toxic action of DIM on development and function of gonad after chronic (28days) aqueous exposure to relatively low doses (0μg/L or 8.5μg/L). Integration of comprehensive toxicogenomic analysis at the transcriptome and proteome levels with apical endpoints, such as production of eggs and swimming performance of larvae, elucidated the molecular linkage in gonad from bottom up along the reproductive adverse outcome pathway. A series of sequential changes at the transcript and protein levels were linked to lesser fecundity and viability of larvae exposed to DIM. Anomalous production of vitellogenin (VTG) and eggshell proteins in testis confirmed the estrogenic potency of DIM. In the ovary, although storage of VTG was greater, lesser expressions of cathepsin enzymes blocked cleavage and incorporation of VTG into oocytes as yolk, which acted together with lower eggshell proteins to inhibit maturation of primary oocyte and thus contributed to impairment of fecundity. Overall, this study demonstrated that exposure to DIM impaired reproductive fitness. Diverse molecular initiating changes in gonads were linked to apical endpoints that could be used in assessment of risks posed by DIM on gametogenesis. In combination with chemical stability and potent endocrine disruption, the results of this study can inform decisions about the use of DIM either as chemopreventive or antifouling agent.

Comparison of toxicogenomic responses to phthalate ester exposure in an organotypic testis co-culture model and responses observed in vivo

We have developed a three-dimensional testicular co-culture system (3D-TCS) which mimics in vivo testes. In this study, transcriptomic responses to phthalate esters (PE's) were compared in the 3D-TCS with responses in rat testes in vivo. Microarray data from the 3D-TCS and from in vivo testes were used to compare changes in gene expression patterns after exposure to developmentally toxic (DTPE) or developmentally non-toxic (DNTPE) phthalate esters. DTPE treatments clustered separately from DNTPE treatments based on principle components analysis both in vitro and in vivo. Pathway analysis using GO-Elite software showed that terms relating to steroid metabolism or reproductive development were enriched both in vitro and in vivo after DTPE exposure. Processes such as cell cycle, cell proliferation and apoptosis were enriched for DTPE treatments in vitro, but not in vivo. Based on these analyses we concluded that transcriptomic responses in the 3D-TCS reflect key aspects of in vivo phthalate toxicity.

Assessment of toxicity and genotoxicity of low doses of 5-fluorouracil in zebrafish (Danio rerio) two-generation study

Residues of anti-neoplastic drugs represent new and emerging pollutants in aquatic environments. Many of these drugs are genotoxic, and it has been postulated that they can cause adverse effects in aquatic ecosystems. 5-Fluorouracil (5-FU) is one of the most extensively used anti-neoplastic drugs in cancer therapy, and this article describes the results of the first investigation using a two-generation toxicity study design with zebrafish (Danio rerio). Exposure of zebrafish to 5-FU (0.01, 1.0 and 100 μg/L) was initiated with adult zebrafish (F0 generation) and continued through the hatchings and adults of the F1 generation, and the hatchings of the F2 generation, to day 33 post-fertilisation. The exposure did not affect survival, growth and reproduction of the zebrafish; however, histopathological changes were observed in the liver and kidney, along with genotoxic effects, at all 5-FU concentrations. Increases in DNA damage determined using the comet assay were significant in the liver and blood cells, but not in the gills and gonads. In erythrocytes, a significant, dose-dependent increase in frequency of micronuclei was observed at all 5-FU concentrations. Whole genome transcriptomic analysis of liver samples of F1 generation zebrafish exposed to 0.01 μg/L and 1 μg/L 5-FU revealed dose-dependent increases in the number of differentially expressed genes, including up-regulation of several DNA-damage-responsive genes and oncogenes (i.e., jun, myca). Although this chronic exposure to environmentally relevant concentrations of 5-FU did not affect the reproduction of the exposed zebrafish, it cannot be excluded that 5-FU can lead to degenerative changes, including cancers, which over long-term exposure of several generations might affect fish populations. The data from this study contribute to a better understanding of the potential consequences of chronic exposure of fish to low concentrations of anti-neoplastic drugs, and they demonstrate that further studies into multi-generation toxicity are needed.

Relationship between head and neck cancer therapy and some genetic endpoints

Head and neck cancer (HNC) is the sixth most common human malignancy worldwide. The main forms of treatment for HNC are surgery, radiotherapy (RT) and chemotherapy (CT). However, the choice of therapy depends on the tumor staging and approaches, which are aimed at organ preservation. Because of systemic RT and CT genotoxicity, one of the important side effects is a secondary cancer that can result from the activity of radiation and antineoplastic drugs on healthy cells. Ionizing radiation can affect the DNA, causing single and double-strand breaks, DNA-protein crosslinks and oxidative damage. The severity of radiotoxicity can be directly associated with the radiation dosimetry and the dose-volume differences. Regarding CT, cisplatin is still the standard protocol for the treatment of squamous cell carcinoma, the most common cancer located in the oral cavity. However, simultaneous treatment with cisplatin, bleomycin and 5-fluorouracil or treatment with paclitaxel and cisplatin are also used. These drugs can interact with the DNA, causing DNA crosslinks, double and single-strand breaks and changes in gene expression. Currently, the late effects of therapy have become a recurring problem, mainly due to the increased survival of HNC patients. Herein, we present an update of the systemic activity of RT and CT for HNC, with a focus on their toxicogenetic and toxicogenomic effects.

An approach for integrating toxicogenomic data in risk assessment: the dibutyl phthalate case study

An approach for evaluating and integrating genomic data in chemical risk assessment was developed based on the lessons learned from performing a case study for the chemical dibutyl phthalate. A case study prototype approach was first developed in accordance with EPA guidance and recommendations of the scientific community. Dibutyl phthalate (DBP) was selected for the case study exercise. The scoping phase of the dibutyl phthalate case study was conducted by considering the available DBP genomic data, taken together with the entire data set, for whether they could inform various risk assessment aspects, such as toxicodynamics, toxicokinetics, and dose-response. A description of weighing the available dibutyl phthalate data set for utility in risk assessment provides an example for considering genomic data for future chemical assessments. As a result of conducting the scoping process, two questions--Do the DBP toxicogenomic data inform 1) the mechanisms or modes of action?, and 2) the interspecies differences in toxicodynamics?--were selected to focus the case study exercise. Principles of the general approach include considering the genomics data in conjunction with all other data to determine their ability to inform the various qualitative and/or quantitative aspects of risk assessment, and evaluating the relationship between the available genomic and toxicity outcome data with respect to study comparability and phenotypic anchoring. Based on experience from the DBP case study, recommendations and a general approach for integrating genomic data in chemical assessment were developed to advance the broader effort to utilize 21st century data in risk assessment.

Evaluation of the dibutyl phthalate toxicogenomic data set

An evaluation of the toxicogenomic data set for dibutyl phthalate (DBP) and male reproductive developmental effects was performed as part of a larger case study to test an approach for incorporating genomic data in risk assessment. The DBP toxicogenomic data set is composed of nine in vivo studies from the published literature that exposed rats to DBP during gestation and evaluated gene expression changes in testes or Wolffian ducts of male fetuses. The exercise focused on qualitative evaluation, based on a lack of available dose-response data, of the DBP toxicogenomic data set to postulate modes and mechanisms of action for the male reproductive developmental outcomes, which occur in the lower dose range. A weight-of-evidence evaluation was performed on the eight DBP toxicogenomic studies of the rat testis at the gene and pathway levels. The results showed relatively strong evidence of DBP-induced downregulation of genes in the steroidogenesis pathway and lipid/sterol/cholesterol transport pathway as well as effects on immediate early gene/growth/differentiation, transcription, peroxisome proliferator-activated receptor signaling and apoptosis pathways in the testis. Since two established modes of action (MOAs), reduced fetal testicular testosterone production and Insl3 gene expression, explain some but not all of the testis effects observed in rats after in utero DBP exposure, other MOAs are likely to be operative. A reanalysis of one DBP microarray study identified additional pathways within cell signaling, metabolism, hormone, disease, and cell adhesion biological processes. These putative new pathways may be associated with DBP effects on the testes that are currently unexplained. This case study on DBP identified data gaps and research needs for the use of toxicogenomic data in risk assessment. Furthermore, this study demonstrated an approach for evaluating toxicogenomic data in human health risk assessment that could be applied to future chemicals.

Use of genomic data in risk assessment case study: I. Evaluation of the dibutyl phthalate male reproductive development toxicity data set

A case study was conducted, using dibutyl phthalate (DBP), to explore an approach to using toxicogenomic data in risk assessment. The toxicity and toxicogenomic data sets relative to DBP-related male reproductive developmental outcomes were considered conjointly to derive information about mode and mechanism of action. In this manuscript, we describe the case study evaluation of the toxicological database for DBP, focusing on identifying the full spectrum of male reproductive developmental effects. The data were assessed to 1) evaluate low dose and low incidence findings and 2) identify male reproductive toxicity endpoints without well-established modes of action (MOAs). These efforts led to the characterization of data gaps and research needs for the toxicity and toxicogenomic studies in a risk assessment context. Further, the identification of endpoints with unexplained MOAs in the toxicity data set was useful in the subsequent evaluation of the mechanistic information that the toxicogenomic data set evaluation could provide. The extensive analysis of the toxicology data set within the MOA context provided a resource of information for DBP in attempts to hypothesize MOAs (for endpoints without a well-established MOA) and to phenotypically anchor toxicogenomic and other mechanistic data both to toxicity endpoints and to available toxicogenomic data. This case study serves as an example of the steps that can be taken to develop a toxicological data source for a risk assessment, both in general and especially for risk assessments that include toxicogenomic data.