Toxicogenomics: the challenges and opportunities to identify biomarkers, signatures and thresholds to support mode-of-action

Genome-wide expression screening in the cardiac embryonic stem cell test shows additional differentiation routes that are regulated by morpholines and piperidines

•

The cardiac embryonic stem cell test showed additional differentiation routes.

•

Morpholines and piperidines regulated the alternative differentiation routes.

•

The gene expression levels help in understanding the applicability domain.

The cardiac embryonic stem cell test (ESTc) is a well-studied non-animal alternative test method based on cardiac cell differentiation inhibition as a measure for developmental toxicity of tested chemicals. In the ESTc, a heterogenic cell population is generated besides cardiomyocytes. Using the full biological domain of ESTc may improve the sensitivity of the test system, possibly broadening the range of chemicals for which developmental effects can be detected in the test. In order to improve our knowledge of the biological and chemical applicability domains of the ESTc, we applied a hypothesis-generating data-driven approach on control samples as follows. A genome-wide expression screening was performed, using Next Generation Sequencing (NGS), to map the range of developmental pathways in the ESTc and to search for a predictive embryotoxicity biomarker profile, instead of the conventional read-out of beating cardiomyocytes. The detected developmental pathways included circulatory system development, skeletal system development, heart development, muscle and organ tissue development, and nervous system and cell development. Two pesticidal chemical classes, the morpholines and piperidines, were assessed for perturbation of differentiation in the ESTc using NGS. In addition to the anticipated impact on cardiomyocyte differentiation, the other developmental pathways were also regulated, in a concentration–response fashion. Despite the structural differences between the morpholine and piperidine pairs, their gene expression effect patterns were largely comparable. In addition, some chemical-specific gene regulation was also observed, which may help with future mechanistic understanding of specific effects with individual test compounds. These similar and unique regulations of gene expression profiles by the test compounds, adds to our knowledge of the chemical applicability domain, specificity and sensitivity of the ESTc. Knowledge of both the biological and chemical applicability domain contributes to the optimal placement of ESTc in test batteries and in Integrated Approaches to Testing and Assessment (IATA).

Functional differences and similarities in activated peripheral blood mononuclear cells by lipopolysaccharide or phytohemagglutinin stimulation between human and cynomolgus monkeys

Background

The monkey is a primary species used in toxicological research. However, the failures of preclinical studies to predict a life-threatening “cytokine storm”, which, for instance, rapidly occurred in six healthy volunteers with the CD28 superagonist monoclonal antibody (mAb) TGN1412 in the first-in-human phase I clinical trial, have emphasized a need to clarify the differences between human and monkey immune systems.

Methods

In the present study, we analyzed and compared the lymphocyte proliferation, cytokine secretion, and gene expression profiles after phytohemagglutinin (PHA) and lipopolysaccharide (LPS) stimulation of peripheral blood mononuclear cells (PBMCs) from three healthy humans and cynomolgus monkeys (Macaca fascicularis).

Results

The results derived from comparison with the corresponding control groups showed that PHA in humans induced a stronger proliferation and wider range of cytokine secretion, along with a greater number of differently expressed genes (DEGs), than when PHA was applied in cynomolgus monkeys. The significant upregulation of genes involved in the mitotic cell cycle, including cyclin B2, TOP2A, TYMS, and CEP55, was observed in human PBMCs with PHA stimulation, while only infrequent or slight upregulation occurred in cynomolgus monkey PBMCs, which may be one of the reasons for a stronger response to PHA in humans. In contrast to PHA, LPS in both species induced a similar proliferation ratio, cytokine profile, and DEG count, suggesting that human and cynomolgus monkeys have a similar response intensity for innate immune responses. Furthermore, 38 and 20 overlapped genes under PHA and LPS stimulation, respectively, were found in both species. These overlapped DEGs were associated with the same biological functions, including DNA replication, mitosis, immune response, chemotaxis, and inflammatory response. Thus, these results might reflect the highly conserved signatures of immune responses to PHA/LPS stimulation across the primates. Moreover, there were some differences in antigen processing and presentation, and the interferon gamma (INF-γ)–mediated signaling pathway in these species detected by gene expression profile study.

Conclusions

In conclusion, this is the first study to compare data on the responses of PBMCs to PHA and LPS in humans versus cynomolgus monkeys, and these findings may provide crucial insights into translating non-human primate (NHP) studies into human trials.

Quantification of cerivastatin toxicity supports organismal performance assays as an effective tool during pharmaceutical safety assessment

A major problem in pharmaceutical development is that adverse effects remain undetected during preclinical and clinical trials, but are later revealed after market release when prescribed to many patients. We have developed a fitness assay known as the organismal performance assay (OPA), which evaluates individual performance by utilizing outbred wild mice (Mus musculus) that are assigned to an exposed or control group, which compete against each other for resources within semi-natural enclosures. Performance measurements included reproductive success, survival, and male competitive ability. Our aim was to utilize cerivastatin (Baycol(®), Bayer), a pharmaceutical with known adverse effects, as a positive control to assess OPAs as a potential tool for evaluating the safety of compounds during preclinical trials. Mice were exposed to cerivastatin (~4.5 mg/kg/day) into early adulthood. Exposure ceased and animals were released into semi-natural enclosures. Within enclosures, cerivastatin-exposed females had 25% fewer offspring and cerivastatin-exposed males had 10% less body mass, occupied 63% fewer territories, sired 41% fewer offspring, and experienced a threefold increase in mortality when compared to controls. OPAs detected several cerivastatin-induced adverse effects indicating that fitness assays, commonly used in ecology and evolutionary biology, could be useful as an additional tool in safety testing during pharmaceutical development.

Risk assessment of Soulatrolide and Mammea (A/BA+A/BB) coumarins from Calophyllum brasiliense by a toxicogenomic and toxicological approach

Calophyllum brasiliense (Calophyllaceae) is a tropical rain forest tree distributed in Central and South America. It is an important source of tetracyclic dipyrano coumarins (Soulatrolide) and Mammea type coumarins. Soulatrolide is a potent inhibitor of HIV-1 reverse transcriptase and displays activity against Mycobacterium tuberculosis. Meanwhile, Mammea A/BA and A/BB, pure or as a mixture, are highly active against several human leukemia cell lines, Trypanosoma cruzi and Leishmania amazonensis. Nevertheless, there are few studies evaluating their safety profile. In the present work we performed toxicogenomic and toxicological analysis for both type of compounds. Soulatrolide, and the Mammea A/BA + A/BB mixture (2.1) were slightly toxic accordingly to Lorke assay classification (DL50 > 3000 mg/kg). After a short-term administration (100 mg/kg/daily, orally, 1 week) liver toxicogenomic analysis revealed 46 up and 72 downregulated genes for Mammea coumarins, and 665 up and 1077 downregulated genes for Soulatrolide. Gene enrichment analysis identified transcripts involved in drug metabolism for both compounds. In addition, network analysis through protein-protein interactions, tissue evaluation by TUNEL assay, and histological examination revealed no tissue damage on liver, kidney and spleen after treatments. Our results indicate that both type of coumarins displayed a safety profile, supporting their use in further preclinical studies to determine its therapeutic potential.

Comparison of toxicogenomics and traditional approaches to inform mode of action and points of departure in human health risk assessment of benzo[a]pyrene in drinking water

Toxicogenomics is proposed to be a useful tool in human health risk assessment. However, a systematic comparison of traditional risk assessment approaches with those applying toxicogenomics has never been done. We conducted a case study to evaluate the utility of toxicogenomics in the risk assessment of benzo[a]pyrene (BaP), a well-studied carcinogen, for drinking water exposures. Our study was intended to compare methodologies, not to evaluate drinking water safety. We compared traditional (RA1), genomics-informed (RA2) and genomics-only (RA3) approaches. RA2 and RA3 applied toxicogenomics data from human cell cultures and mice exposed to BaP to determine if these data could provide insight into BaP's mode of action (MOA) and derive tissue-specific points of departure (POD). Our global gene expression analysis supported that BaP is genotoxic in mice and allowed the development of a detailed MOA. Toxicogenomics analysis in human lymphoblastoid TK6 cells demonstrated a high degree of consistency in perturbed pathways with animal tissues. Quantitatively, the PODs for traditional and transcriptional approaches were similar (liver 1.2 vs. 1.0 mg/kg-bw/day; lungs 0.8 vs. 3.7 mg/kg-bw/day; forestomach 0.5 vs. 7.4 mg/kg-bw/day). RA3, which applied toxicogenomics in the absence of apical toxicology data, demonstrates that this approach provides useful information in data-poor situations. Overall, our study supports the use of toxicogenomics as a relatively fast and cost-effective tool for hazard identification, preliminary evaluation of potential carcinogens, and carcinogenic potency, in addition to identifying current limitations and practical questions for future work.

Highly Predictive Reprogramming of tRNA Modifications Is Linked to Selective Expression of Codon-Biased Genes

Cells respond to stress by controlling gene expression at several levels, with little known about the role of translation. Here, we demonstrate a coordinated translational stress response system involving stress-specific reprogramming of tRNA wobble modifications that leads to selective translation of codon-biased mRNAs representing different classes of critical response proteins. In budding yeast exposed to four oxidants and five alkylating agents, tRNA modification patterns accurately distinguished among chemically similar stressors, with 14 modified ribonucleosides forming the basis for a data-driven model that predicts toxicant chemistry with >80% sensitivity and specificity. tRNA modification subpatterns also distinguish S_N1 from S_N2 alkylating agents, with S_N2-induced increases in m³C in tRNA mechanistically linked to selective translation of threonine-rich membrane proteins from genes enriched with ACC and ACT degenerate codons for threonine. These results establish tRNA modifications as predictive biomarkers of exposure and illustrate a novel regulatory mechanism for translational control of cell stress response.

DNA microarray-based gene expression profiling of estrogenic chemicals

We summarize updated information about DNA microarray-based gene expression profiling by focusing on its application to estrogenic chemicals. First, estrogenic chemicals, including natural/industrial estrogens and phytoestrogens, and the methods for detection and evaluation of estrogenic chemicals were overviewed along with a comprehensive list of estrogenic chemicals of natural or industrial origin. Second, gene expression profiling of chemicals using a focused microarray containing estrogen-responsive genes is summarized. Third, silent estrogens, a new type of estrogenic chemicals characterized by their estrogenic gene expression profiles without growth stimulative or inhibitory effects, have been identified so far exclusively by DNA microarray assay. Lastly, the prospect of a microarray assay is discussed, including issues such as commercialization, future directions of applications and quality control methods.

Phenobarbital and propiconazole toxicogenomic profiles in mice show major similarities consistent with the key role that constitutive androstane receptor (CAR) activation plays in their mode of action

Toxicogenomics (TGx) is employed frequently to investigate underlying molecular mechanisms of the compound of interest and, thus, has become an aid to mode of action determination. However, the results and interpretation of a TGx dataset are influenced by the experimental design and methods of analysis employed. This article describes an evaluation and reanalysis, by two independent laboratories, of previously published TGx mouse liver microarray data for a triazole fungicide, propiconazole (PPZ), and the anticonvulsant drug phenobarbital (PB). Propiconazole produced an increase incidence of liver tumors in male CD-1 mice only at a dose that exceeded the maximum tolerated dose (2500 ppm). Firstly, we illustrate how experimental design differences between two in vivo studies with PPZ and PB may impact the comparisons of TGx results. Secondly, we demonstrate that different researchers using different pathway analysis tools can come to different conclusions on specific mechanistic pathways, even when using the same datasets. Finally, despite these differences the results across three different analyses also show a striking degree of similarity observed for PPZ and PB treated livers when the expression data are viewed as major signaling pathways and cell processes affected. Additional studies described here show that the postulated key event of hepatocellular proliferation was observed in CD-1 mice for both PPZ and PB, and that PPZ is also a potent activator of the mouse CAR nuclear receptor. Thus, with regard to the events which are hallmarks of CAR-induced effects that are key events in the mode of action (MOA) of mouse liver carcinogenesis with PB, PPZ-induced tumors can be viewed as being promoted by a similar PB-like CAR-dependent MOA.

Differential reconstructed gene interaction networks for deriving toxicity threshold in chemical risk assessment

Background

Pathway alterations reflected as changes in gene expression regulation and gene interaction can result from cellular exposure to toxicants. Such information is often used to elucidate toxicological modes of action. From a risk assessment perspective, alterations in biological pathways are a rich resource for setting toxicant thresholds, which may be more sensitive and mechanism-informed than traditional toxicity endpoints. Here we developed a novel differential networks (DNs) approach to connect pathway perturbation with toxicity threshold setting.

Methods

Our DNs approach consists of 6 steps: time-series gene expression data collection, identification of altered genes, gene interaction network reconstruction, differential edge inference, mapping of genes with differential edges to pathways, and establishment of causal relationships between chemical concentration and perturbed pathways. A one-sample Gaussian process model and a linear regression model were used to identify genes that exhibited significant profile changes across an entire time course and between treatments, respectively. Interaction networks of differentially expressed (DE) genes were reconstructed for different treatments using a state space model and then compared to infer differential edges/interactions. DE genes possessing differential edges were mapped to biological pathways in databases such as KEGG pathways.

Results

Using the DNs approach, we analyzed a time-series Escherichia coli live cell gene expression dataset consisting of 4 treatments (control, 10, 100, 1000 mg/L naphthenic acids, NAs) and 18 time points. Through comparison of reconstructed networks and construction of differential networks, 80 genes were identified as DE genes with a significant number of differential edges, and 22 KEGG pathways were altered in a concentration-dependent manner. Some of these pathways were perturbed to a degree as high as 70% even at the lowest exposure concentration, implying a high sensitivity of our DNs approach.

Conclusions

Findings from this proof-of-concept study suggest that our approach has a great potential in providing a novel and sensitive tool for threshold setting in chemical risk assessment. In future work, we plan to analyze more time-series datasets with a full spectrum of concentrations and sufficient replications per treatment. The pathway alteration-derived thresholds will also be compared with those derived from apical endpoints such as cell growth rate.

Analysis of the transcriptome in molecular epidemiology studies

The human transcriptome is complex, comprising multiple transcript types, mostly in the form of non-coding RNA (ncRNA). The majority of ncRNA is of the long form (lncRNA, ≥ 200 bp), which plays an important role in gene regulation through multiple mechanisms including epigenetics, chromatin modification, control of transcription factor binding, and regulation of alternative splicing. Both mRNA and ncRNA exhibit additional variability in the form of alternative splicing and RNA editing. All aspects of the human transcriptome can potentially be dysregulated by environmental exposures. Next-generation RNA sequencing (RNA-Seq) is the best available methodology to measure this although it has limitations, including experimental bias. The third phase of the MicroArray Quality Control Consortium project (MAQC-III), also called Sequencing Quality Control (SeQC), aims to address these limitations through standardization of experimental and bioinformatic methodologies. A limited number of toxicogenomic studies have been conducted to date using RNA-Seq. This review describes the complexity of the human transcriptome, the application of transcriptomics by RNA-Seq or microarray in molecular epidemiology studies, and limitations of these approaches including the type of cell or tissue analyzed, experimental variation, and confounding. By using good study designs with precise, individual exposure measurements, sufficient power and incorporation of phenotypic anchors, studies in human populations can identify biomarkers of exposure and/or early effect and elucidate mechanisms of action underlying associated diseases, even at low doses. Analysis of datasets at the pathway level can compensate for some of the limitations of RNA-Seq and, as more datasets become available, will increasingly elucidate the exposure-disease continuum.

Research ethics in the post-genomic era

New high-throughput 'omics techniques are providing exciting opportunities in clinical medicine and toxicology, especially in the development of biomarkers. In health science research there are traditional ethical considerations that are reasonably obvious, like balancing health benefits and health risks, autonomy mainly pursued by informed consent, and protecting privacy. Epidemiological studies applying new large-scale approaches (e.g., high-throughput or high-content methods and global studies that utilize biobanking of samples and produce large-scale datasets) present new challenges that call for re-evaluation of standard ethical considerations. In this context, assessment of the ethics underlying study designs, bioinformatics, and statistics applied in the generation and clinical translation of research results should also be considered. Indeed, there are ethical considerations in the research process itself, in research objectives and how research is pursued (e.g., which methodologies are selected and how they are carried out). Maintaining research integrity is critical, as demonstrated by the relatively frequent retraction of scientific papers following violations of good scientific practice. Abiding by the laws is necessary but not sufficient for good research ethics, which is and remains in the hands of the scientific community at the level of both individual scientists and organizations. Senior scientists are responsible for the transfer of research tradition to the next generation of scientists through education, mentorship, and setting an example by their own behavior, as well as by creating systems in institutions that support good research ethics.

Gene expression profiling to identify potentially relevant disease outcomes and support human health risk assessment for carbon black nanoparticle exposure

New approaches are urgently needed to evaluate potential hazards posed by exposure to nanomaterials. Gene expression profiling provides information on potential modes of action and human relevance, and tools have recently become available for pathway-based quantitative risk assessment. The objective of this study was to use toxicogenomics in the context of human health risk assessment. We explore the utility of toxicogenomics in risk assessment, using published gene expression data from C57BL/6 mice exposed to 18, 54 and 162 μg Printex 90 carbon black nanoparticles (CBNP). Analysis of CBNP-perturbed pathways, networks and transcription factors revealed concomitant changes in predicted phenotypes (e.g., pulmonary inflammation and genotoxicity), that correlated with dose and time. Benchmark doses (BMDs) for apical endpoints were comparable to minimum BMDs for relevant pathway-specific expression changes. Comparison to inflammatory lung disease models (i.e., allergic airway inflammation, bacterial infection and tissue injury and fibrosis) and human disease profiles revealed that induced gene expression changes in Printex 90 exposed mice were similar to those typical for pulmonary injury and fibrosis. Very similar fibrotic pathways were perturbed in CBNP-exposed mice and human fibrosis disease models. Our synthesis demonstrates how toxicogenomic profiles may be used in human health risk assessment of nanoparticles and constitutes an important step forward in the ultimate recognition of toxicogenomic endpoints in human health risk. As our knowledge of molecular pathways, dose-response characteristics and relevance to human disease continues to grow, we anticipate that toxicogenomics will become increasingly useful in assessing chemical toxicities and in human health risk assessment.