Elements and development processes for test methods in toxicology and human health-relevant life science research.

Many laboratory procedures generate data on properties of chemicals, but they cannot be equated with toxicological "test methods". This apparent discrepancy is not limited to in vitro testing, using animal-free new approach methods (NAM), but also applies to animal-based testing approaches. Here, we give a brief overview of the differences between data generation and the setup or use of a complete test method. While there is excellent literature available on this topic for specialists (GIVIMP guidance; ToxTemp overview), a brief overview and easily-accessible entry point may be useful for a broader community. We provide a single figure to summarize all test method elements and processes required in the development (setup and adaptation) of a test method. The exposure scheme, the endpoint, and the test system are briefly outlined as fundamental elements of any test method. A rationale is provided, why they are not sufficient. We then explain the importance and role of purpose definition (including some information on what is modelled) and the prediction model, aka data interpretation procedure, which depends on the purpose definition, as further essential elements. This connection exemplifies that all fundamental elements are interdependent, and none can be omitted. Finally, discussion is provided on validation as a measure to provide confidence in the reliability, performance, and relevance of a test method. In this sense, validation may be considered a sixth fundamental element for practical use of test methods.

Many laboratory procedures generate data on chemicals, but they cannot be considered complete toxicological "test methods". Here, we give a brief explanation of the fundamental elements of a toxicological test method. We provide an illustration that gives a complete overview of the devel­opment of a test method for non-specialists. We introduce the six fundamental elements, i.e., the exposure scheme, the test endpoint, the test system, the purpose definition and the prediction model and describe how they work together. Finally, we discuss the concept of validation. An understanding of these concepts is important for good-quality scientific research and especially for the development and acceptance of alternatives to animal experiments.

Identification of the bacterial metabolite aerugine as potential trigger of human dopaminergic neurodegeneration.

The causes of nigrostriatal cell death in idiopathic Parkinson's disease are unknown, but exposure to toxic chemicals may play some role. We followed up here on suggestions that bacterial secondary metabolites might be selectively cytotoxic to dopaminergic neurons. Extracts from Streptomyces venezuelae were found to kill human dopaminergic neurons (LUHMES cells). Utilizing this model system as a bioassay, we identified a bacterial metabolite known as aerugine (C10H11NO2S; 2-[4-(hydroxymethyl)-4,5-dihydro-1,3-thiazol-2-yl]phenol) and confirmed this finding by chemical re-synthesis. This 2-hydroxyphenyl-thiazoline compound was previously shown to be a product of a wide-spread biosynthetic cluster also found in the human microbiome and in several pathogens. Aerugine triggered half-maximal dopaminergic neurotoxicity at 3-4 µM. It was less toxic for other neurons (10-20 µM), and non-toxic (at <100 µM) for common human cell lines. Neurotoxicity was completely prevented by several iron chelators, by distinct anti-oxidants and by a caspase inhibitor. In the Caenorhabditis elegans model organism, general survival was not affected by aerugine concentrations up to 100 µM. When transgenic worms, expressing green fluorescent protein only in their dopamine neurons, were exposed to aerugine, specific neurodegeneration was observed. The toxicant also exerted functional dopaminergic toxicity in nematodes as determined by the "basal slowing response" assay. Thus, our research has unveiled a bacterial metabolite with a remarkably selective toxicity toward human dopaminergic neurons in vitro and for the dopaminergic nervous system of Caenorhabditis elegans in vivo. These findings suggest that microbe-derived environmental chemicals should be further investigated for their role in the pathogenesis of Parkinson's disease.

Acceptance criteria for new approach methods in toxicology and human health-relevant life science research - part I.

Every test procedure, scientific and non-scientific, has inherent uncertainties, even when performed according to a standard operating procedure (SOP). In addition, it is prone to errors, defects, and mistakes introduced by operators, laboratory equipment, or materials used. Adherence to an SOP and comprehensive validation of the test method cannot guarantee that each test run produces data within the acceptable range of variability and with the precision and accuracy determined during the method validation. We illustrate here (part I) why controlling the validity of each test run is an important element of experimental design. The definition and application of acceptance criteria (AC) for the validity of test runs is important for the setup and use of test methods, particularly for the use of new approach methods (NAM) in toxicity testing. AC can be used for decision rules on how to handle data, e.g., to accept the data for further use (AC fulfilled) or to reject the data (AC not fulfilled). The adherence to AC has important requirements and consequences that may seem surprising at first sight: (i) AC depend on a test method's objectives, e.g., on the types/concentrations of chemicals tested, the regulatory context, the desired throughput; (ii) AC are applied and documented at each test run, while validation of a method (including the definition of AC) is only performed once; (iii) if AC are altered, then the set of data produced by a method can change. AC, if missing, are the blind spot of quality assurance: Test results may not be reliable and comparable. The establishment and uses of AC will be further detailed in part II of this series.

Establishment of a human cell-based in vitro battery to assess developmental neurotoxicity hazard of chemicals.

Developmental neurotoxicity (DNT) is a major safety concern for all chemicals of the human exposome. However, DNT data from animal studies are available for only a small percentage of manufactured compounds. Test methods with a higher throughput than current regulatory guideline methods, and with improved human relevance are urgently needed. We therefore explored the feasibility of DNT hazard assessment based on new approach methods (NAMs). An in vitro battery (IVB) was assembled from ten individual NAMs that had been developed during the past years to probe effects of chemicals on various fundamental neurodevelopmental processes. All assays used human neural cells at different developmental stages. This allowed us to assess disturbances of: (i) proliferation of neural progenitor cells (NPC); (ii) migration of neural crest cells, radial glia cells, neurons and oligodendrocytes; (iii) differentiation of NPC into neurons and oligodendrocytes; and (iv) neurite outgrowth of peripheral and central neurons. In parallel, cytotoxicity measures were obtained. The feasibility of concentration-dependent screening and of a reliable biostatistical processing of the complex multi-dimensional data was explored with a set of 120 test compounds, containing subsets of pre-defined positive and negative DNT compounds. The battery provided alerts (hit or borderline) for 24 of 28 known toxicants (82% sensitivity), and for none of the 17 negative controls. Based on the results from this screen project, strategies were developed on how IVB data may be used in the context of risk assessment scenarios employing integrated approaches for testing and assessment (IATA).

Transcriptomic-based evaluation of trichloroethylene glutathione and cysteine conjugates demonstrate phenotype-dependent stress responses in a panel of human in vitro models.

Environmental or occupational exposure of humans to trichloroethylene (TCE) has been associated with different extrahepatic toxic effects, including nephrotoxicity and neurotoxicity. Bioactivation of TCE via the glutathione (GSH) conjugation pathway has been proposed as underlying mechanism, although only few mechanistic studies have used cell models of human origin. In this study, six human derived cell models were evaluated as in vitro models representing potential target tissues of TCE-conjugates: RPTEC/TERT1 (kidney), HepaRG (liver), HUVEC/TERT2 (vascular endothelial), LUHMES (neuronal, dopaminergic), human induced pluripotent stem cells (hiPSC) derived peripheral neurons (UKN5) and hiPSC-derived differentiated brain cortical cultures containing all subtypes of neurons and astrocytes (BCC42). A high throughput transcriptomic screening, utilizing mRNA templated oligo-sequencing (TempO-Seq), was used to study transcriptomic effects after exposure to TCE-conjugates. Cells were exposed to a wide range of concentrations of S-(1,2-trans-dichlorovinyl)glutathione (1,2-DCVG), S-(1,2-trans-dichlorovinyl)-L-cysteine (1,2-DCVC), S-(2,2-dichlorovinyl)glutathione (2,2-DCVG), and S-(2,2-dichlorovinyl)-L-cysteine (2,2-DCVC). 1,2-DCVC caused stress responses belonging to the Nrf2 pathway and Unfolded protein response in all the tested models but to different extents. The renal model was the most sensitive model to both 1,2-DCVC and 1,2-DCVG, with an early Nrf2-response at 3 µM and hundreds of differentially expressed genes at higher concentrations. Exposure to 2,2-DCVG and 2,2-DCVC also resulted in the upregulation of Nrf2 pathway genes in RPTEC/TERT1 although at higher concentrations. Of the three neuronal models, both the LUHMES and BCC42 showed significant Nrf2-responses and at higher concentration UPR-responses, supporting recent hypotheses that 1,2-DCVC may be involved in neurotoxic effects of TCE. The cell models with the highest expression of γ-glutamyltransferase (GGT) enzymes, showed cellular responses to both 1,2-DCVG and 1,2-DCVC. Little to no effects were found in the neuronal models from 1,2-DCVG exposure due to their low GGT-expression. This study expands our knowledge on tissue specificity of TCE S-conjugates and emphasizes the value of human cell models together with transcriptomics for such mechanistic studies.

Generation of Human Nociceptor-Enriched Sensory Neurons for the Study of Pain-Related Dysfunctions.

In vitro models of the peripheral nervous system would benefit from further refinements to better support studies on neuropathies. In particular, the assessment of pain-related signals is still difficult in human cell cultures. Here, we harnessed induced pluripotent stem cells (iPSCs) to generate peripheral sensory neurons enriched in nociceptors. The objective was to generate a culture system with signaling endpoints suitable for pharmacological and toxicological studies. Neurons generated by conventional differentiation protocols expressed moderate levels of P2X3 purinergic receptors and only low levels of TRPV1 capsaicin receptors, when maturation time was kept to the upper practically useful limit of 6 weeks. As alternative approach, we generated cells with an inducible NGN1 transgene. Ectopic expression of this transcription factor during a defined time window of differentiation resulted in highly enriched nociceptor cultures, as determined by functional (P2X3 and TRPV1 receptors) and immunocytochemical phenotyping, complemented by extensive transcriptome profiling. Single cell recordings of Ca2+-indicator fluorescence from >9000 cells were used to establish the "fraction of reactive cells" in a stimulated population as experimental endpoint, that appeared robust, transparent and quantifiable. To provide an example of application to biomedical studies, functional consequences of prolonged exposure to the chemotherapeutic drug oxaliplatin were examined at non-cytotoxic concentrations. We found (i) neuronal (allodynia-like) hypersensitivity to otherwise non-activating mechanical stimulation that could be blocked by modulators of voltage-gated sodium channels; (ii) hyper-responsiveness to TRPV1 receptor stimulation. These findings and several other measured functional alterations indicate that the model is suitable for pharmacological and toxicological studies related to peripheral neuropathies.

Specific Attenuation of Purinergic Signaling during Bortezomib-Induced Peripheral Neuropathy In Vitro.

Human peripheral neuropathies are poorly understood, and the availability of experimental models limits further research. The PeriTox test uses immature dorsal root ganglia (DRG)-like neurons, derived from induced pluripotent stem cells (iPSC), to assess cell death and neurite damage. Here, we explored the suitability of matured peripheral neuron cultures for the detection of sub-cytotoxic endpoints, such as altered responses of pain-related P2X receptors. A two-step differentiation protocol, involving the transient expression of ectopic neurogenin-1 (NGN1) allowed for the generation of homogeneous cultures of sensory neurons. After >38 days of differentiation, they showed a robust response (Ca2+-signaling) to the P2X3 ligand α,ß-methylene ATP. The clinical proteasome inhibitor bortezomib abolished the P2X3 signal at ≥5 nM, while 50−200 nM was required in the PeriTox test to identify neurite damage and cell death. A 24 h treatment with low nM concentrations of bortezomib led to moderate increases in resting cell intracellular Ca2+ concentration but signaling through transient receptor potential V1 (TRPV1) receptors or depolarization-triggered Ca2+ influx remained unaffected. We interpreted the specific attenuation of purinergic signaling as a functional cell stress response. A reorganization of tubulin to form dense structures around the cell somata confirmed a mild, non-cytotoxic stress triggered by low concentrations of bortezomib. The proteasome inhibitors carfilzomib, delanzomib, epoxomicin, and MG-132 showed similar stress responses. Thus, the model presented here may be used for the profiling of new proteasome inhibitors in regard to their side effect (neuropathy) potential, or for pharmacological studies on the attenuation of their neurotoxicity. P2X3 signaling proved useful as endpoint to assess potential neurotoxicants in peripheral neurons.

Neurodevelopmental toxicity assessment of flame retardants using a human DNT in vitro testing battery.

Due to their neurodevelopmental toxicity, flame retardants (FRs) like polybrominated diphenyl ethers are banned from the market and replaced by alternative FRs, like organophosphorus FRs, that have mostly unknown toxicological profiles. To study their neurodevelopmental toxicity, we evaluated the hazard of several FRs including phased-out polybrominated FRs and organophosphorus FRs: 2,2',4,4'-tetrabromodiphenylether (BDE-47), 2,2',4,4',5-pentabromodiphenylether (BDE-99), tetrabromobisphenol A, triphenyl phosphate, tris(2-butoxyethyl) phosphate and its metabolite bis-(2-butoxyethyl) phosphate, isodecyl diphenyl phosphate, triphenyl isopropylated phosphate, tricresyl phosphate, tris(1,3-dichloro-2-propyl) phosphate, tert-butylphenyl diphenyl phosphate, 2-ethylhexyl diphenyl phosphate, tris(1-chloroisopropyl) phosphate, and tris(2-chloroethyl) phosphate. Therefore, we used a human cell-based developmental neurotoxicity (DNT) in vitro battery covering a large variety of neurodevelopmental endpoints. Potency according to the respective most sensitive benchmark concentration (BMC) across the battery ranked from <1 µM (5 FRs), 1<10 µM (7 FRs) to the >10 µM range (3 FRs). Evaluation of the data with the ToxPi tool revealed a distinct ranking (a) than with the BMC and (b) compared to the ToxCast data, suggesting that DNT hazard of these FRs is not well predicted by ToxCast assays. Extrapolating the DNT in vitro battery BMCs to human FR exposure via breast milk suggests low risk for individual compounds. However, it raises a potential concern for real-life mixture exposure, especially when different compounds converge through diverse modes-of-action on common endpoints, like oligodendrocyte differentiation in this study. This case study using FRs suggests that human cell-based DNT in vitro battery is a promising approach for neurodevelopmental hazard assessment and compound prioritization in risk assessment.

The EU-ToxRisk method documentation, data processing and chemical testing pipeline for the regulatory use of new approach methods.

Hazard assessment, based on new approach methods (NAM), requires the use of batteries of assays, where individual tests may be contributed by different laboratories. A unified strategy for such collaborative testing is presented. It details all procedures required to allow test information to be usable for integrated hazard assessment, strategic project decisions and/or for regulatory purposes. The EU-ToxRisk project developed a strategy to provide regulatorily valid data, and exemplified this using a panel of > 20 assays (with > 50 individual endpoints), each exposed to 19 well-known test compounds (e.g. rotenone, colchicine, mercury, paracetamol, rifampicine, paraquat, taxol). Examples of strategy implementation are provided for all aspects required to ensure data validity: (i) documentation of test methods in a publicly accessible database; (ii) deposition of standard operating procedures (SOP) at the European Union DB-ALM repository; (iii) test readiness scoring accoding to defined criteria; (iv) disclosure of the pipeline for data processing; (v) link of uncertainty measures and metadata to the data; (vi) definition of test chemicals, their handling and their behavior in test media; (vii) specification of the test purpose and overall evaluation plans. Moreover, data generation was exemplified by providing results from 25 reporter assays. A complete evaluation of the entire test battery will be described elsewhere. A major learning from the retrospective analysis of this large testing project was the need for thorough definitions of the above strategy aspects, ideally in form of a study pre-registration, to allow adequate interpretation of the data and to ensure overall scientific/toxicological validity.

CaFFEE: A program for evaluating time courses of Ca2+ dependent signal changes of complex cells loaded with fluorescent indicator dyes.

Quantification of changes in intracellular free Ca2+ concentrations [Ca2+]i is fundamental to the understanding of the physiology of single cells in response to both environmental and endogenous stimuli. Here we present easy to use freeware that allows especially the evaluation of [Ca2+]i signals in complex and mixed cultures. The program CaFFEE (Calcium Fluorescent Flash Evaluating Engine) enables the user to evaluate the response of hundreds of cells to treat-ments that influence [Ca2+]i. CaFFEE processes large quantities of image data, automatically identifies individual cells in mixed, heterogeneous populations, and evaluates their fluorescence signal. All data are exported in spreadsheet format, and data on thousands of cells can be batch-processed. Moreover, the program optimizes the visual representation of time-lapse image data for user-guided data exploration (setting of parameters for semi-automated data processing). The freeware allows the standardized and transparent processing of imaging data independent of the platform used to gen-erate the data.

Lapachol acetylglycosylation enhances its cytotoxic and pro-apoptotic activities in HL60 cells.

Lapachol is a plant-derived naphthoquinone that kills several types of cancer cells. Derivatives of this molecule may therefore prove to be useful chemotherapeutic agents. In this study, we explored whether glycosylation increases the cytotoxic potency of lapachol towards HL-60 human leukemia cells. Two beta-glycosides were synthesized and characterized: LA4A (lapachol-ß-glucoside) and LA4C (lapachol-N-acetylglucosamine-ß-glucoside). The sugar moieties of both novel molecules were per-acetylated to facilitate cellular uptake. The IC50 values (in µM) for LA4A (5.7) and LA4C (5.3) were lower than those for lapachol (25). LA4A and LA4C triggered typical signs of apoptosis, such as the exposure of phosphatidylserine on the outside of cells, chromatin condensation, DNA fragmentation and a decrease of the mitochondrial transmembrane potential (ΔΨm) prior to cell lysis. Moreover, DNA fragmentation triggered by the lapachol-glycosides was reduced by pre-treatment with the caspase inhibitor, z-VAD-fmk. While LA4A and LA4C activated caspases-3, -8 and -9, lapachol failed to activate these apoptotic proteases, even when used at high concentrations. Finally, the toxicity of lapachol and its derivatives was also tested on non-tumor cells. We used human peripheral neurons (PeriTox test) to evaluate the side effect potential of these compounds. LA4C was clearly less toxic than LA4A. We conclude that LA4C had the most favorable profile as drug candidate (high tumor cell toxicity, reduced neurotoxicity). In general, this study shows that the cytotoxicity of lapachol towards HL-60 can be enhanced by glycosylation, and that the therapeutic ratio may be modified by the type of sugar added.

The influence of structural gradients in large pore organosilica materials on the capabilities for hosting cellular communities.

Cells exist in the so-called extracellular matrix (ECM) in their native state, and numerous future applications require reliable and potent ECM-mimics. A perspective, which goes beyond ECM emulation, is the design of a host-material with features which are not accessible in the biological portfolio. Such a feature would, for instance, be the creation of a structural or chemical gradient, and to explore how this special property influences the biological processes. First, we wanted to test if macroporous organosilica materials with appropriate surface modification can act as a host for the implementation of human cells like HeLa or LUHMES. It was possible to use a commercially available polymeric foam as a scaffold and coat it with a thiophenol-containing organosilica layer, followed by biofunctionalization with biotin using click chemistry and the subsequent coupling of streptavidin-fibronectin to it. More importantly, deformation of the scaffold allowed the generation of a permanent structural gradient. In this work, we show that the structural gradient has a tremendous influence on the capability of the described material for the accommodation of living cells. The introduction of a bi-directional gradient enabled the establishment of a cellular community comprising different cell types in spatially distinct regions of the material. An interesting perspective is to study communication between cell types or to create cellular communities, which can never exist in a natural environment.

Time and space-resolved quantification of plasma membrane sialylation for measurements of cell function and neurotoxicity.

While there are many methods to quantify the synthesis, localization, and pool sizes of proteins and DNA during physiological responses and toxicological stress, only few approaches allow following the fate of carbohydrates. One of them is metabolic glycoengineering (MGE), which makes use of chemically modified sugars (CMS) that enter the cellular biosynthesis pathways leading to glycoproteins and glycolipids. The CMS can subsequently be coupled (via bio-orthogonal chemical reactions) to tags that are quantifiable by microscopic imaging. We asked here, whether MGE can be used in a quantitative and time-resolved way to study neuronal glycoprotein synthesis and its impairment. We focused on the detection of sialic acid (Sia), by feeding human neurons the biosynthetic precursor N-acetyl-mannosamine, modified by an azide tag. Using this system, we identified non-toxic conditions that allowed live cell labeling with high spatial and temporal resolution, as well as the quantification of cell surface Sia. Using combinations of immunostaining, chromatography, and western blotting, we quantified the percentage of cellular label incorporation and effects on glycoproteins such as polysialylated neural cell adhesion molecule. A specific imaging algorithm was used to quantify Sia incorporation into neuronal projections, as potential measure of complex cell function in toxicological studies. When various toxicants were studied, we identified a subgroup (mitochondrial respiration inhibitors) that affected neurite glycan levels several hours before any other viability parameter was affected. The MGE-based neurotoxicity assay, thus allowed the identification of subtle impairments of neurochemical function with very high sensitivity.

Development of a neural rosette formation assay (RoFA) to identify neurodevelopmental toxicants and to characterize their transcriptome disturbances.

The first in vitro tests for developmental toxicity made use of rodent cells. Newer teratology tests, e.g. developed during the ESNATS project, use human cells and measure mechanistic endpoints (such as transcriptome changes). However, the toxicological implications of mechanistic parameters are hard to judge, without functional/morphological endpoints. To address this issue, we developed a new version of the human stem cell-based test STOP-tox(UKN). For this purpose, the capacity of the cells to self-organize to neural rosettes was assessed as functional endpoint: pluripotent stem cells were allowed to differentiate into neuroepithelial cells for 6 days in the presence or absence of toxicants. Then, both transcriptome changes were measured (standard STOP-tox(UKN)) and cells were allowed to form rosettes. After optimization of staining methods, an imaging algorithm for rosette quantification was implemented and used for an automated rosette formation assay (RoFA). Neural tube toxicants (like valproic acid), which are known to disturb human development at stages when rosette-forming cells are present, were used as positive controls. Established toxicants led to distinctly different tissue organization and differentiation stages. RoFA outcome and transcript changes largely correlated concerning (1) the concentration-dependence, (2) the time dependence, and (3) the set of positive hits identified amongst 24 potential toxicants. Using such comparative data, a prediction model for the RoFA was developed. The comparative analysis was also used to identify gene dysregulations that are particularly predictive for disturbed rosette formation. This 'RoFA predictor gene set' may be used for a simplified and less costly setup of the STOP-tox(UKN) assay.

A structure-activity relationship linking non-planar PCBs to functional deficits of neural crest cells: new roles for connexins.

Migration of neural crest cells (NCC) is a fundamental developmental process, and test methods to identify interfering toxicants have been developed. By examining cell function endpoints, as in the 'migration-inhibition of NCC (cMINC)' assay, a large number of toxicity mechanisms and protein targets can be covered. However, the key events that lead to the adverse effects of a given chemical or group of related compounds are hard to elucidate. To address this issue, we explored here, whether the establishment of two overlapping structure-activity relationships (SAR)-linking chemical structure on the one hand to a phenotypic test outcome, and on the other hand to a mechanistic endpoint-was useful as strategy to identify relevant toxicity mechanisms. For this purpose, we chose polychlorinated biphenyls (PCB) as a large group of related, but still toxicologically and physicochemically diverse structures. We obtained concentration-dependent data for 26 PCBs in the cMINC assay. Moreover, the test chemicals were evaluated by a new high-content imaging method for their effect on cellular re-distribution of connexin43 and for their capacity to inhibit gap junctions. Non-planar PCBs inhibited NCC migration. The potency (1-10 µM) correlated with the number of ortho-chlorine substituents; non-ortho-chloro (planar) PCBs were non-toxic. The toxicity to NCC partially correlated with gap junction inhibition, while it fully correlated (p < 0.0004) with connexin43 cellular re-distribution. Thus, our double-SAR strategy revealed a mechanistic step tightly linked to NCC toxicity of PCBs. Connexin43 patterns in NCC may be explored as a new endpoint relevant to developmental toxicity screening.

Stem Cell Transcriptome Responses and Corresponding Biomarkers That Indicate the Transition from Adaptive Responses to Cytotoxicity.

Analysis of transcriptome changes has become an established method to characterize the reaction of cells to toxicants. Such experiments are mostly performed at compound concentrations close to the cytotoxicity threshold. At present, little information is available on concentration-dependent features of transcriptome changes, in particular, at the transition from noncytotoxic concentrations to conditions that are associated with cell death. Thus, it is unclear in how far cell death confounds the results of transcriptome studies. To explore this gap of knowledge, we treated pluripotent stem cells differentiating to human neuroepithelial cells (UKN1 assay) for short periods (48 h) with increasing concentrations of valproic acid (VPA) and methyl mercury (MeHg), two compounds with vastly different modes of action. We developed various visualization tools to describe cellular responses, and the overall response was classified as "tolerance" (minor transcriptome changes), "functional adaptation" (moderate/strong transcriptome responses, but no cytotoxicity), and "degeneration". The latter two conditions were compared, using various statistical approaches. We identified (i) genes regulated at cytotoxic, but not at noncytotoxic, concentrations and (ii) KEGG pathways, gene ontology term groups, and superordinate biological processes that were only regulated at cytotoxic concentrations. The consensus markers and processes found after 48 h treatment were then overlaid with those found after prolonged (6 days) treatment. The study highlights the importance of careful concentration selection and of controlling viability for transcriptome studies. Moreover, it allowed identification of 39 candidate "biomarkers of cytotoxicity". These could serve to provide alerts that data sets of interest may have been affected by cell death in the model system studied.