Toxicogenomics and systems toxicology: aims and prospects

Inhibitory effects of black phosphorus nanosheets on tumor cell proliferation through downregulation of ADIPOQ and downstream signaling pathways

Exposure to environmental pollutants, including nanomaterials, has a significant impact on tumor progression. The increased demand for black phosphorus nanosheets (BPNSs), driven by their exceptional properties, raises concerns about potential environmental contamination. Assessing their toxicity on tumor growth is essential. Herein, we employed a range of biological techniques, including cytotoxicity measurement, bioinformatics tools, proteomics, target gene overexpression, Western blot analysis, and apoptosis detection, to investigate the toxicity of BPNSs across A549, HepG-2, MCF-7, and Caco-2 cell lines. Our results demonstrated that BPNSs downregulated the expression of ADIPOQ and its associated downstream pathways, such as AMP-activated protein kinase (AMPK), nuclear factor erythroid 2-related factor 2 (Nrf2), and other unidentified pathways. These downregulated pathways ultimately led to mitochondria-dependent apoptosis. Notably, the specific downstream pathways involved varied depending on the type of tumors. These insightful findings not only confirm the consistent inhibitory effects of BPNSs across different tumor cells, but also elucidate the cytotoxicity mechanisms of BPNSs in different tumors, providing valuable information for their safe application and health risk assessment.

Curcumin alleviates Aflatoxin B1-triggered chicken liver necroptosis by targeting the LOC769044/miR-1679/STAT1 axis

Aflatoxin B1 (AFB1) is an unavoidable environmental toxin. The accumulation of AFB1 and its metabolites in the liver poses a threat to both human and animal health. Curcumin exhibits anti-oxidative, anti-tumor, and anti-inflammatory properties. There is no report on the mechanism regarding how curcumin relived liver necroptosis in chickens induced by AFB1 based on the regulatory network of ceRNA. To explore this, we performed transmission electron microscopy and sequenced lncRNA and mRNA in chicken livers treated with AFB1 and/or curcumin for 28 d in vivo. We observed substantial alterations in the lncRNA and mRNA expression profiles within the chicken liver, indicating that curcumin can mitigate AFB1-induced necroptosis both in vivo and in vitro. Further analysis, including the establishment of an lncRNA-miRNA-mRNA network and the utilization of a dual luciferase reporter assay, revealed that LOC769044 acts as a competing endogenous RNA (ceRNA) for miR-1679. In addition, STAT1 was identified as a direct target of miR-1679. Modulating miR-1679 levels through overexpression, and silencing LOC769044 and STAT1, effectively reversed the necroptotic effects induced by AFB1, a reversal that was also observed with curcumin supplementation. In conclusion, our data demonstrate that curcumin alleviates AFB1-induced liver necroptosis through the LOC769044/miR-1679/STAT1 signaling axis. This study suggests that LOC769044 may serve as a novel therapeutic target for managing AFB1-mediated liver toxicity.

The application of natural language processing for the extraction of mechanistic information in toxicology

To study the ways in which compounds can induce adverse effects, toxicologists have been constructing Adverse Outcome Pathways (AOPs). An AOP can be considered as a pragmatic tool to capture and visualize mechanisms underlying different types of toxicity inflicted by any kind of stressor, and describes the interactions between key entities that lead to the adverse outcome on multiple biological levels of organization. The construction or optimization of an AOP is a labor intensive process, which currently depends on the manual search, collection, reviewing and synthesis of available scientific literature. This process could however be largely facilitated using Natural Language Processing (NLP) to extract information contained in scientific literature in a systematic, objective, and rapid manner that would lead to greater accuracy and reproducibility. This would support researchers to invest their expertise in the substantive assessment of the AOPs by replacing the time spent on evidence gathering by a critical review of the data extracted by NLP. As case examples, we selected two frequent adversities observed in the liver: namely, cholestasis and steatosis denoting accumulation of bile and lipid, respectively. We used deep learning language models to recognize entities of interest in text and establish causal relationships between them. We demonstrate how an NLP pipeline combining Named Entity Recognition and a simple rules-based relationship extraction model helps screen compounds related to liver adversities in the literature, but also extract mechanistic information for how such adversities develop, from the molecular to the organismal level. Finally, we provide some perspectives opened by the recent progress in Large Language Models and how these could be used in the future. We propose this work brings two main contributions: 1) a proof-of-concept that NLP can support the extraction of information from text for modern toxicology and 2) a template open-source model for recognition of toxicological entities and extraction of their relationships. All resources are openly accessible via GitHub (https://github.com/ontox-project/en-tox).

The alternation of halobenzoquinone disinfection byproduct on toxicogenomics of DNA damage and repair in uroepithelial cells

Halobenzoquinones (HBQs) were recently discovered as an emerging class of drinking water disinfection byproducts with carcinogenic concern. However, the molecular mechanism underlying HBQs-induced DNA damage is not clear. In this study, we integrated in vitro genotoxicity, computational toxicology, and the quantitative toxicogenomic analysis of HBQs on DNA damage/repair pathways in human bladder epithelial cells SV-HUC-1. The results showed that HBQs could induce cytotoxicity with the descending order as 2,6-DIBQ > 2,6-DCBQ ≈ 2,6-DBBQ. Also, HBQs can increase DNA damage in SV-HUC-1 cells and thus generate genotoxicity. However, there is no significant difference in genotoxicity among the three HBQs. The results of molecular docking and molecular dynamics simulation further confirmed that HBQs had high binding fractions and stability to DNA. Toxicogenomic analysis indicated that HBQs interfered with DNA repair pathways, mainly affecting base excision repair, nucleotide excision repair and homologous recombination repair. These results have provided new insights into the underlying molecular mechanisms of HBQs-induced DNA damage, and contributed to the understanding of the relationship between exposure to DBPs and risks of developing bladder cancer.

High-Throughput Transcriptomics Differentiates Toxic versus Non-Toxic Chemical Exposures Using a Rat Liver Model

To address the challenge of limited throughput with traditional toxicity testing, a newly developed high-throughput transcriptomics (HTT) platform, together with a 5-day in vivo rat model, offers an alternative approach to estimate chemical exposures and provide reasonable estimates of toxicological endpoints. This study contains an HTT analysis of 18 environmental chemicals with known liver toxicity. They were evaluated using male Sprague Dawley rats exposed to various concentrations daily for five consecutive days via oral gavage, with data collected on the sixth day. Here, we further explored the 5-day rat model to identify potential gene signatures that can differentiate between toxic and non-toxic liver responses and provide us with a potential histopathological endpoint of chemical exposure. We identified a distinct gene expression pattern that differentiated non-hepatotoxic compounds from hepatotoxic compounds in a dose-dependent manner, and an analysis of the significantly altered common genes indicated that toxic chemicals predominantly upregulated most of the genes and several pathways in amino acid and lipid metabolism. Finally, our liver injury module analysis revealed that several liver-toxic compounds showed similarities in the key injury phenotypes of cellular inflammation and proliferation, indicating potential molecular initiating processes that may lead to a specific end-stage liver disease.

Advancing Computational Toxicology by Interpretable Machine Learning

Chemical toxicity evaluations for drugs, consumer products, and environmental chemicals have a critical impact on human health. Traditional animal models to evaluate chemical toxicity are expensive, time-consuming, and often fail to detect toxicants in humans. Computational toxicology is a promising alternative approach that utilizes machine learning (ML) and deep learning (DL) techniques to predict the toxicity potentials of chemicals. Although the applications of ML- and DL-based computational models in chemical toxicity predictions are attractive, many toxicity models are "black boxes" in nature and difficult to interpret by toxicologists, which hampers the chemical risk assessments using these models. The recent progress of interpretable ML (IML) in the computer science field meets this urgent need to unveil the underlying toxicity mechanisms and elucidate the domain knowledge of toxicity models. In this review, we focused on the applications of IML in computational toxicology, including toxicity feature data, model interpretation methods, use of knowledge base frameworks in IML development, and recent applications. The challenges and future directions of IML modeling in toxicology are also discussed. We hope this review can encourage efforts in developing interpretable models with new IML algorithms that can assist new chemical assessments by illustrating toxicity mechanisms in humans.

Application of CRISPR/Cas9-based genome editing in ecotoxicology

Chemicals are widely used and released into the environment, and their degradation, accumulation, migration, and transformation processes in the environment can pose a threat to the ecosystem. The advancement in analytical methods with high-throughput screening of biomolecules has revolutionized the way toxicologists used to explore the effects of chemicals on organisms. CRISPR/Cas is a newly developed tool, widely used in the exploration of basic science and biologically engineered products given its high efficiency and low cost. For example, it can edit target genes efficiently, and save loss of the crop yield caused by environmental pollution as well as gain a better understanding of the toxicity mechanisms from various chemicals. This review briefly introduces the development history of CRISPR/Cas and summarizes the current application of CRISPR/Cas in ecotoxicology, including its application on improving crop yield and drug resistance towards agricultural pollution, antibiotic pollution and other threats. The benefits by applying the CRISPR/Cas9 system in conventional toxicity mechanism studies are fully demonstrated here together with its foreseeable expansions in other area of ecotoxicology. Finally, the prospects and disadvantages of CRISPR/Cas system in the field of ecotoxicology are also discussed.

Transcriptomic analysis reveals responses to a polluted sediment in the Mediterranean copepod Acartia clausi

Marine sediments are regarded as sinks for several classes of contaminants. Characterization and effects of sediments on marine biota now require a multidisciplinary approach, which includes chemical and ecotoxicological analyses and molecular biomarkers. Here, a gene expression study was performed to measure the response of adult females of the Mediterranean copepod Acartia clausi to elutriates of polluted sediments (containing high concentrations of polycyclic aromatic hydrocarbons, PAHs, and heavy metals) from an industrial area in the Southern Tyrrhenian Sea (Bagnoli-Coroglio). Functional annotation of the A. clausi transcriptome generated as reference here, showed a good quality of the assembly and great homology with other copepod and crustacean sequences in public databases. This is one of the few available transcriptomic resources for this widespread copepod species of great ecological relevance in temperate coastal areas. Differential expression analysis between females exposed to the elutriate and those in control seawater identified 1000 differentially expressed genes, of which 743 up- and 257 down-regulated. Within the up-regulated genes, the most represented functions were related to proteolysis (lysosomal protease, peptidase, cathepsin), response to stress and detoxification (heat-shock protein, superoxide dismutase, glutathione-S-transferase, cytochrome P450), and cytoskeleton structure (α- and β-tubulin). Down-regulated genes were mostly involved with ribosome structure (ribosomal proteins) and DNA binding (histone proteins, transcription factors). Overall, these results suggest that processes such as transcription, translation, protein degradation, metabolism of biomolecules, reproduction, and xenobiotic detoxification were altered in the copepod in response to polluted elutriates. In conclusion, our results contribute to gaining information on the transcriptomic responses of copepods to polluted sediments. They will also prompt the selection of genes of interest to be used as biomarkers of exposure to PAHs and heavy metals in molecular toxicology studies on copepods, and in general, in comparative functional genomic studies on marine zooplankton.

Systems toxicology of complex wood combustion aerosol reveals gaseous carbonyl compounds as critical constituents

Epidemiological studies identified air pollution as one of the prime causes for human morbidity and mortality, due to harmful effects mainly on the cardiovascular and respiratory systems. Damage to the lung leads to several severe diseases such as fibrosis, chronic obstructive pulmonary disease and cancer. Noxious environmental aerosols are comprised of a gas and particulate phase representing highly complex chemical mixtures composed of myriads of compounds. Although some critical pollutants, foremost particulate matter (PM), could be linked to adverse health effects, a comprehensive understanding of relevant biological mechanisms and detrimental aerosol constituents is still lacking. Here, we employed a systems toxicology approach focusing on wood combustion, an important source for air pollution, and demonstrate a key role of the gas phase, specifically carbonyls, in driving adverse effects. Transcriptional profiling and biochemical analysis of human lung cells exposed at the air-liquid-interface determined DNA damage and stress response, as well as perturbation of cellular metabolism, as major key events. Connectivity mapping revealed a high similarity of gene expression signatures induced by wood smoke and agents prompting DNA-protein crosslinks (DPCs). Indeed, various gaseous aldehydes were detected in wood smoke, which promote DPCs, initiate similar genomic responses and are responsible for DNA damage provoked by wood smoke. Hence, systems toxicology enables the discovery of critical constituents of complex mixtures i.e. aerosols and highlights the role of carbonyls on top of particulate matter as an important health hazard.

Expanding adverse outcome pathways towards one health models for nanosafety

The ever-growing production of nano-enabled products has generated the need for dedicated risk assessment strategies that ensure safety for humans and the environment. Transdisciplinary approaches are needed to support the development of new technologies while respecting environmental limits, as also highlighted by the EU Green Deal Chemicals Strategy for Sustainability and its safe and sustainable by design (SSbD) framework. The One Health concept offers a holistic multiscale approach for the assessment of nanosafety. However, toxicology is not yet capable of explaining the interaction between chemicals and biological systems at the multiscale level and in the context of the One Health framework. Furthermore, there is a disconnect between chemical safety assessment, epidemiology, and other fields of biology that, if unified, would enable the adoption of the One Health model. The development of mechanistic toxicology and the generation of omics data has provided important biological knowledge of the response of individual biological systems to nanomaterials (NMs). On the other hand, epigenetic data have the potential to inform on interspecies mechanisms of adaptation. These data types, however, need to be linked to concepts that support their intuitive interpretation. Adverse Outcome Pathways (AOPs) represent an evolving framework to anchor existing knowledge to chemical risk assessment. In this perspective, we discuss the possibility of integrating multi-level toxicogenomics data, including toxicoepigenetic insights, into the AOP framework. We anticipate that this new direction of toxicogenomics can support the development of One Health models applicable to groups of chemicals and to multiple species in the tree of life.

Toxicogenomic analysis of physiologically important metals: An integrated in silico approach

Biologically important metals regulate cellular homeostasis in living systems. Anthropogenic exposure to these metals can cause adverse effects, including an increased incidence of diseases in humans such as cancer, lung, and cardiovascular defects. However, the effects of metals and the common genes/signaling pathways involved in metal toxicity have not been elucidated. Hence, the present study used toxicogenomic data mining with the comparative toxicogenomics database to explore the impact of these metals. The metals were categorized into transition, alkali, and alkali earth. The common genes were identified and subjected to functional enrichment analysis. Further, gene-gene and protein-protein interactions were assessed. Also, the top ten transcription factors and miRNAs that regulate the genes were identified. The phenotypes and diseases that have increased incidence upon alterations of these genes were detected. Overall, we were able to identify IL1B and SOD2 as the common genes, along with the AGE-RAGE signaling pathway in diabetic complications as the common pathway altered. Enriched genes and pathways specific to each metal category were also found. Further, we identified heart failure as the major diseases that have increased the incidence of these metals' exposure. In conclusion, exposure to essential metals might cause adverse effects via inflammation and oxidative stress.

Comprehensive investigation of hepatotoxicity of the mixture containing phthalates and bisphenol A

Connections between the mixture containing bis(2- ethylhexyl) phthalate (DEHP), dibutyl phthalate (DBP) and bisphenol A (BPA) and liver injury were explored through in silico investigation and 2 in vivo models. Comparative Toxicogenomics Database (CTD), ShinyGO, ToppCluster and Cytoscape were used for bioinformatic analysis. In vivo subacute study was performed on rats - five groups (n = 6): (1) Control: corn oil, (2) DEHP: 50 mg/kg b.w./day, (3) DBP: 50 mg/kg b.w./day, (4) BPA: 25 mg/kg b.w./day, (5) MIX: DEHP + DBP + BPA. Zebrafish embryos were exposed to the investigated substances in different doses, singularly and combined (binary and ternary mixtures). Liver injury was linked to 75 DEHP, DBP, and BPA genes, mostly connected to inflammation/oxidative stress. In rats, significant alterations in redox status/bioelements and pathohistology were most notable or exclusively present in MIX (probable additive effects). BPA decreased liver area (LA) index in dose-dependent manner. DEHP (< 2 µg/mL) and DBP (≤ 5 µg/mL) reduced LA values, while their higher doses increased LA index. The effect of DBP in binary mixtures led to a lethal outcome at the two highest concentrations, while the hepatotoxicity of DEHP/DBP/BPA mixture was dictated by BPA (confirmed by the benchmark dose analysis).

Astrocyte-Like Cells Transcriptome Changes After Exposure to a Low and Non-cytotoxic MeHg Concentration

The central nervous system is the main target of MeHg toxicity and glial cells are the first line of defense; however, their true role remains unclear. This study aimed to identify the global map of human glial-like (U87) cells transcriptome after exposure to a non-toxic and non-lethal MeHg concentration and to investigate the related molecular changes. U87 cells were exposed upon 0.1, 0.5, and 1 µM MeHg for 4 and 24 h. Although no changes were observed in the percentage of viable cells, the metabolic viability was significantly decreased after exposure to 1 µM MeHg for 24 h; thus, the non-toxic concentration of 0.1 µM MeHg was chosen to perform microarray analysis. Significant changes in U87 cells transcriptome were observed only after 24 h. The expression of 392 genes was down regulated while 431 genes were up-regulated. Gene ontology showed alterations in biological processes (75%), cellular components (21%), and molecular functions (4%). The main pathways showed by KEGG and Reactome were cell cycle regulation and Rho GTPase signaling. The complex mechanism of U87 cells response against MeHg exposure indicates that even a low and non-toxic concentration is able to alter the gene expression profile.

How the exposure to environmentally relevant pesticide formulations affects the expression of stress response genes and its relation to oxidative damage and genotoxicity in Caiman latirostris

This study aimed to analyze the molecular stress responses thought the expression levels of catalase (cat), superoxide dismutase (sod) and heat shock protein 70 (hsp70) genes, and how these relate with cellular stress response considering oxidative damage to lipids, DNA and genotoxicity in blood of Caiman latirostris hatchlings exposed to pesticide formulations under ex situ conditions. Treatments were: negative control (NC-tap water), glyphosate 2% (GLY), cypermethrin 0.12% (CYP), chlorpyrifos 0.8% (CPF), and their ternary mixture (Mx3). The concentrations and schedule of application were those recommended in soybean crops. Soil and water showed pesticides residues in all exposed groups. Results showed a statistically significant increase in the micronucleus frequency and DNA damage, with an important oxidation in all exposed groups. The expression level of cat gene was significantly higher in CYP while the expression of hsp70 was significantly lower in GLY, CYP and Mx3, compared to NC. Pesticides tested showed alterations in expression levels, growth parameters, DNA damage and base oxidation under realistic exposure conditions, and can threaten, in the long term, the health status of wild populations.

Toxicogenomics Data for Chemical Safety Assessment and Development of New Approach Methodologies: An Adverse Outcome Pathway-Based Approach

Mechanistic toxicology provides a powerful approach to inform on the safety of chemicals and the development of safe-by-design compounds. Although toxicogenomics supports mechanistic evaluation of chemical exposures, its implementation into the regulatory framework is hindered by uncertainties in the analysis and interpretation of such data. The use of mechanistic evidence through the adverse outcome pathway (AOP) concept is promoted for the development of new approach methodologies (NAMs) that can reduce animal experimentation. However, to unleash the full potential of AOPs and build confidence into toxicogenomics, robust associations between AOPs and patterns of molecular alteration need to be established. Systematic curation of molecular events to AOPs will create the much-needed link between toxicogenomics and systemic mechanisms depicted by the AOPs. This, in turn, will introduce novel ways of benefitting from the AOPs, including predictive models and targeted assays, while also reducing the need for multiple testing strategies. Hence, a multi-step strategy to annotate AOPs is developed, and the resulting associations are applied to successfully highlight relevant adverse outcomes for chemical exposures with strong in vitro and in vivo convergence, supporting chemical grouping and other data-driven approaches. Finally, a panel of AOP-derived in vitro biomarkers for pulmonary fibrosis (PF) is identified and experimentally validated.

Dietary Selenomethionine Reduce Mercury Tissue Levels and Modulate Methylmercury Induced Proteomic and Transcriptomic Alterations in Hippocampi of Adolescent BALB/c Mice

Methylmercury (MeHg) is a well-known environmental contaminant, particularly harmful to the developing brain. The main human dietary exposure to MeHg occurs through seafood consumption. However, seafood also contains several nutrients, including selenium, which has been shown to interact with MeHg and potentially ameliorate its toxicity. The aim of this study was to investigate the combined effects of selenium (as selenomethionine; SeMet) and MeHg on mercury accumulation in tissues and the effects concomitant dietary exposure of these compounds exert on the hippocampal proteome and transcriptome in mice. Adolescent male BALB/c mice were exposed to SeMet and two different doses of MeHg through their diet for 11 weeks. Organs, including the brain, were sampled for mercury analyses. Hippocampi were collected and analyzed using proteomics and transcriptomics followed by multi-omics bioinformatics data analysis. The dietary presence of SeMet reduced the amount of mercury in several organs, including the brain. Proteomic and RNA-seq analyses showed that both protein and RNA expression patterns were inversely regulated in mice receiving SeMet together with MeHg compared to MeHg alone. Several pathways, proteins and RNA transcripts involved in conditions such as immune responses and inflammation, oxidative stress, cell plasticity and Alzheimer’s disease were affected inversely by SeMet and MeHg, indicating that SeMet can ameliorate several toxic effects of MeHg in mice.

Transcriptomic and Histological Analysis of the Greentail Prawn (Metapenaeus bennettae) Following Light Crude Oil Exposure

Oil spills pose a significant threat to marine biodiversity. Crude oil can partition into sediments where it may be persistent, placing benthic species such as decapods at particular risk of exposure. Transcriptomic and histological tools are often used to investigate the effects of hydrocarbon exposure on marine organisms following oil spill events, allowing for the identification of metabolic pathways impacted by oil exposure. However, there is limited information available for decapod crustaceans, many of which carry significant economic value. In the present study, we assess the sublethal impacts of crude oil exposure in the commercially important Australian greentail prawn (Metapenaeus bennettae) using transcriptomic and histological analyses. Prawns exposed to light, unweathered crude oil "spiked" sediments for 90 h were transferred to clean sediments for a further 72 h to assess recovery. Chemical analyses indicated that polycyclic aromatic hydrocarbons increased by approximately 65% and 91% in prawn muscle following 24 and 90 h of exposure, respectively, and significantly decreased during 24- and 72-h recovery periods. Transcriptomic responses followed an exposure and recovery pattern with innate immunity and nutrient metabolism transcripts significantly lowered in abundance after 24 h of exposure and were higher in abundance after 72 h of recovery. In addition, transcription/translation, cellular responses, and DNA repair pathways were significantly impacted after 24 h of exposure and recovered after 72 h of recovery. However, histological alterations such as tubule atrophy indicated an increase in severity after 24 and 72 h of recovery. The present study provides new insights into the sublethal impacts of crude oil exposure in greentail prawns and identifies molecular pathways altered by exposure. We expect these findings to inform future management associated with oil extraction activity and spills. Environ Toxicol Chem 2022;41:2162-2180. © 2022 John Wiley & Sons Ltd. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Joint impact of key air pollutants on COVID-19 severity: prediction based on toxicogenomic data analysis

Considering that some researchers point to a possible influence of air pollution on COVID-19 transmission, severity, and death rate, the aim of our in silico study was to determine the relationship between the key air pollutants [sulphur dioxide (SO), carbon monoxide (CO), ₂particulate matter (PM_x), nitrogen dioxide (NO₂), and ozone (O₃)] and COVID-19 complications using the publicly available toxicogenomic analytical and prediction tools: (i) Comparative Toxicogenomic Database (CTD) to identify genes common to air pollutants and COVID-19 complications; (ii) GeneMANIA to construct a network of these common and related genes; (iii) ToppGene Suite to extract the most important biological processes and molecular pathways; and (iv) DisGeNET to search for the top gene-disease pairs. SO₂, CO, PM_x, NO₂, and O₃ interacted with 6, 6, 18, 9, and 12 COVID-19-related genes, respectively. Four of these are common for all pollutants (IL10, IL6, IL1B, and TNF) and participate in most (77.64 %) physical interactions. Further analysis pointed to cytokine binding and cytokine-mediated signalling pathway as the most important molecular function and biological process, respectively. Other molecular functions and biological processes are mostly related to cytokine activity and inflammation, which might be connected to the cytokine storm and resulting COVID-19 complications. The final step singled out the link between the CEBPA gene and acute myelocytic leukaemia and between TNFRSF1A and TNF receptor-associated periodic fever syndrome. This indicates possible complications in COVID-19 patients suffering from these diseases, especially those living in urban areas with poor air quality.

First transcriptome of the copepod Gladioferens pectinatus subjected to chronic contaminant exposures

Contaminants are often at low concentrations in ecosystems and their effects on exposed organisms can occur over long periods of time and across multiple generations. Alterations to subcellular mechanistic pathways in response to exposure to contaminants can provide insights into mechanisms of toxicity that methods measuring higher levels of biological may miss. Analysis of the whole transcriptome can identify novel mechanisms of action leading to impacts in exposed biota. The aim of this study was to characterise how exposures to copper, benzophenone and diclofenac across multiple generations altered molecular expression pathways in the marine copepod Gladioferens pectinatus. Results of the study demonstrated differential gene expression was observed in cultures exposure to diclofenac (569), copper (449) and benzophenone (59). Pathways linked to stress, growth, cellular and metabolic processes were altered by exposure to all three contaminants with genes associated with oxidative stress and xenobiotic regulation also impacted. Protein kinase functioning, cytochrome P450, transcription, skeletal muscle contraction/relaxation, mitochondrial phosphate translocator, protein synthesis and mitochondrial methylation were all differentially expressed with all three chemicals. The results of the study also suggested that using dimethyl sulfoxide as a dispersant influenced the transcriptome and future research may want to investigate it's use in molecular studies. Data generated in this study provides a first look at transcriptomic response of G. pectinatus exposed to contaminants across multiple generations, future research is needed to validate the identified biomarkers and link these results to apical responses such as population growth to demonstrate the predictive capacity of molecular tools.

Single-Larva RNA Sequencing Identifies Markers of Copper Toxicity and Exposure in Early Mytilus californianus Larvae

One of the challenges facing efforts to generate molecular biomarkers for toxins is distinguishing between markers that are indicative of exposure and markers that provide evidence of the effects of toxicity. Phenotypic anchoring provides an approach to help segregate markers into these categories based on some phenotypic index of toxicity. Here we leveraged the mussel embryo-larval toxicity assay in which toxicity is estimated by the fraction of larvae that exhibit an abnormal morphology, to isolate subsets of larvae that were abnormal and thus showed evidence of copper-toxicity, versus others that while exposed to copper exhibited normal morphology. Mussel larvae reared under control conditions or in the presence of increasing levels of copper (3-15 μg/L Cu²⁺) were physically sorted according to whether their morphology was normal or abnormal, and then profiled using RNAseq. Supervised differential expression analysis identified sets of genes whose differential expression was specific to the pools of abnormal larvae versus normal larvae, providing putative markers of copper toxicity versus exposure. Markers of copper exposure and copper-induced abnormality were involved in many of the same pathways, including development, shell formation, cell adhesion, and oxidative stress, yet unique markers were detected in each gene set. Markers of effect appeared to be more resolving between phenotypes at the lower copper concentration, while markers of exposure were informative at both copper concentrations.

An integrative approach to define chemical exposure threshold limits for endangered sea turtles

Environmental contaminants pose serious health threats to marine megafauna species, yet methods defining exposure threshold limits are lacking. Here, a three-pillar chemical risk assessment framework is presented based on (1) species- and chemical-specific lifetime bioaccumulation modelling, (2) non-destructive in vitro and in vivo toxicity threshold assessment, and (3) chemical risk quantification. We used the effects of cadmium (Cd) in green sea turtles (Chelonia mydas) as a proof of concept to evaluate the quantitative mechanistic modelling approach. A physiologically-based kinetic (PBK) model simulated Cd tissue concentrations (liver, kidney, muscle, fat, brain, scute, and 'rest of the body') in C.mydas. The validated PBK model then translated species-specific in vitro results to in vivo effects. The results showed that the resilience of C.mydas towards Cd kidney toxicity is age-dependent and differs with changing physiology and feeding ecology. Using the model in reverse mode, a steady-state exposure threshold of 0.1 µg/g dry weight Cd in forage was derived and compared to real-world exposure scenarios. Three out of the four globally distinct C.mydas populations assessed are exposed to Cd levels above this threshold limit. This approach can be adapted to other marine species and chemicals to prioritize measures for managing potentially harmful chemical exposures.

Trends in the Application of "Omics" to Ecotoxicology and Stress Ecology

Our ability to predict and assess how environmental changes such as pollution and climate change affect components of the Earth's biome is of paramount importance. This need positioned the fields of ecotoxicology and stress ecology at the center of environmental monitoring efforts. Advances in these interdisciplinary fields depend not only on conceptual leaps but also on technological advances and data integration. High-throughput "omics" technologies enabled the measurement of molecular changes at virtually all levels of an organism's biological organization and thus continue to influence how the impacts of stressors are understood. This bibliometric review describes literature trends (2000-2020) that indicate that more different stressors than species are studied each year but that only a few stressors have been studied in more than two phyla. At the same time, the molecular responses of a diverse set of non-model species have been investigated, but cross-species comparisons are still rare. While transcriptomics studies dominated until 2016, a shift towards proteomics and multiomics studies is apparent. There is now a wealth of data at functional omics levels from many phylogenetically diverse species. This review, therefore, addresses the question of how to integrate omics information across species.

Ecotoxicogenomic analysis of zebrafish embryos exposed to triclosan and mixture triclosan and methyl triclosan using suppression subtractive hybridization and next-generation sequencing

Triclosan (TCS) and methyl-triclosan (MTCS), an environmental transformation product of biocide of TCS, have been detected in water, sediment, fish, and invertebrates. In this study, the key pathway perturbation in zebrafish (Danio rerio) embryos exposed to TCS (300 μg/L) and TCS/MTCS mixture (300 μg/L TCS + 30 μg/L MTCS) was assessed by integrating the metabolomic and transcriptomic dysregulation. The differential expressed genes (DEGs) were obtained from the subtracted cDNA libraries by using the suppression subtractive hybridization and next-generation sequencing approach. The dysregulation of twenty-eight GO terms and four KEGG pathways, including oxidative phosphorylation and cardiac muscle contraction, were shown in the TCS treatment group, indicating that TCS could disrupt the mitochondrial inner membrane function by downshifting the electrochemical gradient. Meanwhile, the addition of MTCS in the exposure would cause fourteen additional significant KEGG pathway changes, demonstrating the different effects between two exposure. A pathway-based analysis using the identified DEGs and the altered metabolites in zebrafish embryos treated with TCS and TCS/MTCS mixture, collectively, has been applied. This study demonstrated that the integration of SSH-NGS and metabolomics could reveal toxic effects and potential diseases associated with the exposures of TCS and MTCS in aquatic environments.

Determining the Biological Mechanisms of Action for Environmental Exposures: Applying CRISPR/Cas9 to Toxicological Assessments

Toxicology is a constantly evolving field, especially in the area of developing alternatives to animal testing. Toxicological research must evolve and utilize adaptive technologies in an effort to improve public, environmental, and occupational health. The most commonly cited mechanisms of toxic action after exposure to a chemical or particle test substance is oxidative stress. However, because oxidative stress involves a plethora of genes and proteins, the exact mechanism(s) are not commonly defined. Exact mechanisms of toxicity can be revealed using an emerging laboratory technique referred to as CRISPR (clustered regularly interspaced short palindromic repeats). This article reviews the most common CRISPR techniques utilized today and how each may be applied in Toxicological Sciences. Specifically, the CRISPR/CRISPR-associated protein complex is used for single gene knock-outs, whereas CRISPR interference/activation is used for silencing or activating (respectively) ribonucleic acid. Finally, CRISPR libraries are used for knocking-out entire gene pathways. This review highlights the application of CRISPR in toxicology to elucidate the exact mechanism through which toxicants perturb normal cellular functions.

Toxicogenomic Profiling of 28 Nanomaterials in Mouse Airways

Toxicogenomics opens novel opportunities for hazard assessment by utilizing computational methods to map molecular events and biological processes. In this study, the transcriptomic and immunopathological changes associated with airway exposure to a total of 28 engineered nanomaterials (ENM) are investigated. The ENM are selected to have different core (Ag, Au, TiO2, CuO, nanodiamond, and multiwalled carbon nanotubes) and surface chemistries (COOH, NH2, or polyethylene glycosylation (PEG)). Additionally, ENM with variations in either size (Au) or shape (TiO2) are included. Mice are exposed to 10 µg of ENM by oropharyngeal aspiration for 4 consecutive days, followed by extensive histological/cytological analyses and transcriptomic characterization of lung tissue. The results demonstrate that transcriptomic alterations are correlated with the inflammatory cell infiltrate in the lungs. Surface modification has varying effects on the airways with amination rendering the strongest inflammatory response, while PEGylation suppresses toxicity. However, toxicological responses are also dependent on ENM core chemistry. In addition to ENM-specific transcriptional changes, a subset of 50 shared differentially expressed genes is also highlighted that cluster these ENM according to their toxicity. This study provides the largest in vivo data set currently available and as such provides valuable information to be utilized in developing predictive models for ENM toxicity.

Proteomics in systems toxicology

Proteins are the ultimate product of gene expression. As they hinge between gene transcription and phenotype, they offer a more realistic perspective of toxicopathic effects, responses and even susceptibility to insult than targeting genes and mRNAs while dodging some inter-individual variability that hinders measuring downstream endpoints like metabolites or enzyme activity. Toxicologists have long focused on proteins as biomarkers but the advent of proteomics shifted risk assessment from narrow single-endpoint analyses to whole-proteome screening, enabling deriving protein-centric adverse outcome pathways (AOPs), which are pivotal for the derivation of Systems Biology informally named Systems Toxicology. Especially if coupled pathology, the identification of molecular initiating events (MIEs) and AOPs allow predictive modeling of toxicological pathways, which now stands as the frontier for the next generation of toxicologists. Advances in mass spectrometry, bioinformatics, protein databases and top-down proteomics create new opportunities for mechanistic and effects-oriented research in all fields, from ecotoxicology to pharmacotoxicology.

Integrated RNA Sequencing Reveals Epigenetic Impacts of Diesel Particulate Matter Exposure in Human Cerebral Organoids

Autism spectrum disorder (ASD) manifests early in childhood. While genetic variants increase risk for ASD, a growing body of literature has established that in utero chemical exposures also contribute to ASD risk. These chemicals include air-based pollutants like diesel particulate matter (DPM). A combination of single-cell and direct transcriptomics of DPM-exposed human-induced pluripotent stem cell-derived cerebral organoids revealed toxicogenomic effects of DPM exposure during fetal brain development. Direct transcriptomics, sequencing RNA bases via Nanopore, revealed that cerebral organoids contain extensive RNA modifications, with DPM-altering cytosine methylation in oxidative mitochondrial transcripts expressed in outer radial glia cells. Single-cell transcriptomics further confirmed an oxidative phosphorylation change in cell groups such as outer radial glia upon DPM exposure. This approach highlights how DPM exposure perturbs normal mitochondrial function and cellular respiration during early brain development, which may contribute to developmental disorders like ASD by altering neurodevelopment.

Combining in vivo pathohistological and redox status analysis with in silico toxicogenomic study to explore the phthalates and bisphenol A mixture-induced testicular toxicity

The aim of this study was to: (i) determine and compare the capacity of bis (2 -ethylhexyl) phthalate (DEHP), dibutyl phthalate (DBP), bisphenol A (BPA), and their mixture to produce testicular toxicity after the subacute exposure; (ii) explore the mechanisms behind the observed changes using in silico toxicogenomic approach. Male rats were randomly split into groups (n = 6): (1) Control (corn oil); (2) DEHP (50 mg/kg b.w./day); (3) DBP (50 mg/kg b.w./day); (4) BPA (25 mg/kg b.w./day); and (5) MIX (50 mg/kg b.w./day DEHP + 50 mg/kg b.w/day DBP + 25 mg/kg b.w./day BPA). Animals were sacrificed after 28 days of oral exposure, testes were extracted and prepared for histological assessments under the light microscope (haematoxylin and eosin staining) and redox status analysis. The Comparative Toxicogenomics Database (CTD; http://CTD.mdibl.org), Cytoscape software (https://cytoscape.org) and ToppGene Suite (https://toppgene.cchmc.org) were used for data-mining. Present pathohistological study has demonstrated more pronounced testicular toxicity of the MIX group (desquamated germinal epithelium cells, enlarged cells with hyperchromatic nuclei, multinucleated cell forms and intracytoplasmic vacuoles) in comparison with the single substances, while effects on redox status parameters were either more prominent, or present only in the MIX group. In silico investigation revealed 20 genes linked to male reproductive disorders, affected by all three investigated substances. Effects on metabolism, AhR pathway, apoptosis and oxidative stress could be singled out as the most probable mechanisms involved in the subacute DEHP, DBP and BPA mixture testicular toxicity, while the effect on oxidative stress parameters was confirmed by in vivo experiment.

Mouse Genetic Reference Populations: Cellular Platforms for Integrative Systems Genetics

Interrogation of disease-relevant cellular and molecular traits exhibited by genetically diverse cell populations enables in vitro systems genetics approaches for uncovering the basic properties of cellular function and identity. Primary cells, stem cells, and organoids derived from genetically diverse mouse strains, such as Collaborative Cross and Diversity Outbred populations, offer the opportunity for parallel in vitro/in vivo screening. These panels provide genetic resolution for variant discovery and functional characterization, as well as disease modeling and in vivo validation capabilities. Here we review mouse cellular systems genetics approaches for characterizing the influence of genetic variation on signaling networks and phenotypic diversity, and we discuss approaches for data integration and cross-species validation.

Neurotoxicity in Marine Invertebrates: An Update

Invertebrates represent about 95% of existing species, and most of them belong to aquatic ecosystems. Marine invertebrates are found at intermediate levels of the food chain and, therefore, they play a central role in the biodiversity of ecosystems. Furthermore, these organisms have a short life cycle, easy laboratory manipulation, and high sensitivity to marine pollution and, therefore, they are considered to be optimal bioindicators for assessing detrimental chemical agents that are related to the marine environment and with potential toxicity to human health, including neurotoxicity. In general, albeit simple, the nervous system of marine invertebrates is composed of neuronal and glial cells, and it exhibits biochemical and functional similarities with the vertebrate nervous system, including humans. In recent decades, new genetic and transcriptomic technologies have made the identification of many neural genes and transcription factors homologous to those in humans possible. Neuroinflammation, oxidative stress, and altered levels of neurotransmitters are some of the aspects of neurotoxic effects that can also occur in marine invertebrate organisms. The purpose of this review is to provide an overview of major marine pollutants, such as heavy metals, pesticides, and micro and nano-plastics, with a focus on their neurotoxic effects in marine invertebrate organisms. This review could be a stimulus to bio-research towards the use of invertebrate model systems other than traditional, ethically questionable, time-consuming, and highly expensive mammalian models.

Safety assessment of drug combinations used in COVID-19 treatment: in silico toxicogenomic data-mining approach

Improvement of COVID-19 clinical condition was seen in studies where combination of antiretroviral drugs, lopinavir and ritonavir, as well as immunomodulant antimalaric, chloroquine/hydroxychloroquine together with the macrolide-type antibiotic, azithromycin, was used for patient's treatment. Although these drugs are "old", their pharmacological and toxicological profile in SARS-CoV-2 - infected patients are still unknown. Thus, by using in silico toxicogenomic data-mining approach, we aimed to assess both risks and benefits of the COVID-19 treatment with the most promising candidate drugs combinations: lopinavir/ritonavir and chloroquine/hydroxychloroquine + azithromycin. The Comparative Toxicogenomics Database (CTD; http://CTD.mdibl.org), Cytoscape software (https://cytoscape.org) and ToppGene Suite portal (https://toppgene.cchmc.org) served as a foundation in our research. Our results have demonstrated that lopinavir/ritonavir increased the expression of the genes involved in immune response and lipid metabolism (IL6, ICAM1, CCL2, TNF, APOA1, etc.). Chloroquine/hydroxychloroquine + azithromycin interacted with 6 genes (CCL2, CTSB, CXCL8, IL1B, IL6 and TNF), whereas chloroquine and azithromycin affected two additional genes (BCL2L1 and CYP3A4), which might be a reason behind a greater number of consequential diseases. In contrast to lopinavir/ritonavir, chloroquine/hydroxychloroquine + azithromycin downregulated the expression of TNF and IL6. As expected, inflammation, cardiotoxicity, and dyslipidaemias were revealed as the main risks of lopinavir/ritonavir treatment, while chloroquine/hydroxychloroquine + azithromycin therapy was additionally linked to gastrointestinal and skin diseases. According to our results, these drug combinations should be administrated with caution to patients suffering from cardiovascular problems, autoimmune diseases, or acquired and hereditary lipid disorders.

Transcriptomics in Toxicogenomics, Part II: Preprocessing and Differential Expression Analysis for High Quality Data

Preprocessing of transcriptomics data plays a pivotal role in the development of toxicogenomics-driven tools for chemical toxicity assessment. The generation and exploitation of large volumes of molecular profiles, following an appropriate experimental design, allows the employment of toxicogenomics (TGx) approaches for a thorough characterisation of the mechanism of action (MOA) of different compounds. To date, a plethora of data preprocessing methodologies have been suggested. However, in most cases, building the optimal analytical workflow is not straightforward. A careful selection of the right tools must be carried out, since it will affect the downstream analyses and modelling approaches. Transcriptomics data preprocessing spans across multiple steps such as quality check, filtering, normalization, batch effect detection and correction. Currently, there is a lack of standard guidelines for data preprocessing in the TGx field. Defining the optimal tools and procedures to be employed in the transcriptomics data preprocessing will lead to the generation of homogeneous and unbiased data, allowing the development of more reliable, robust and accurate predictive models. In this review, we outline methods for the preprocessing of three main transcriptomic technologies including microarray, bulk RNA-Sequencing (RNA-Seq), and single cell RNA-Sequencing (scRNA-Seq). Moreover, we discuss the most common methods for the identification of differentially expressed genes and to perform a functional enrichment analysis. This review is the second part of a three-article series on Transcriptomics in Toxicogenomics.

The Power of Resolution: Contextualized Understanding of Biological Responses to Liver Injury Chemicals Using High-throughput Transcriptomics and Benchmark Concentration Modeling

Prediction of human response to chemical exposures is a major challenge in both pharmaceutical and toxicological research. Transcriptomics has been a powerful tool to explore chemical-biological interactions, however, limited throughput, high-costs, and complexity of transcriptomic interpretations have yielded numerous studies lacking sufficient experimental context for predictive application. To address these challenges, we have utilized a novel high-throughput transcriptomics (HTT) platform, TempO-Seq, to apply the interpretive power of concentration-response modeling with exposures to 24 reference compounds in both differentiated and non-differentiated human HepaRG cell cultures. Our goals were to (1) explore transcriptomic characteristics distinguishing liver injury compounds, (2) assess impacts of differentiation state of HepaRG cells on baseline and compound-induced responses (eg, metabolically-activated), and (3) identify and resolve reference biological-response pathways through benchmark concentration (BMC) modeling. Study data revealed the predictive utility of this approach to identify human liver injury compounds by their respective BMCs in relation to human internal exposure plasma concentrations, and effectively distinguished drug analogs with varied associations of human liver injury (eg, withdrawn therapeutics trovafloxacin and troglitazone). Impacts of cellular differentiation state (proliferated vs differentiated) were revealed on baseline drug metabolizing enzyme expression, hepatic receptor signaling, and responsiveness to metabolically-activated toxicants (eg, cyclophosphamide, benzo(a)pyrene, and aflatoxin B1). Finally, concentration-response modeling enabled efficient identification and resolution of plausibly-relevant biological-response pathways through their respective pathway-level BMCs. Taken together, these findings revealed HTT paired with differentiated in vitro liver models as an effective tool to model, explore, and interpret toxicological and pharmacological interactions.

Transcriptomics in Toxicogenomics, Part III: Data Modelling for Risk Assessment

Transcriptomics data are relevant to address a number of challenges in Toxicogenomics (TGx). After careful planning of exposure conditions and data preprocessing, the TGx data can be used in predictive toxicology, where more advanced modelling techniques are applied. The large volume of molecular profiles produced by omics-based technologies allows the development and application of artificial intelligence (AI) methods in TGx. Indeed, the publicly available omics datasets are constantly increasing together with a plethora of different methods that are made available to facilitate their analysis, interpretation and the generation of accurate and stable predictive models. In this review, we present the state-of-the-art of data modelling applied to transcriptomics data in TGx. We show how the benchmark dose (BMD) analysis can be applied to TGx data. We review read across and adverse outcome pathways (AOP) modelling methodologies. We discuss how network-based approaches can be successfully employed to clarify the mechanism of action (MOA) or specific biomarkers of exposure. We also describe the main AI methodologies applied to TGx data to create predictive classification and regression models and we address current challenges. Finally, we present a short description of deep learning (DL) and data integration methodologies applied in these contexts. Modelling of TGx data represents a valuable tool for more accurate chemical safety assessment. This review is the third part of a three-article series on Transcriptomics in Toxicogenomics.

Factors Affecting the Perception of New Approach Methodologies (NAMs) in the Ecotoxicology Community

Given current legislative mandates to assess the safety of thousands of chemicals and the slow pace at which conventional testing proceeds, there is a need to accelerate chemical risk assessment. Governments and businesses are increasingly interested in new approach methodologies (NAMs) that promise to reduce costs and delays. We explore 5 sociological factors within the ecotoxicology community that can influence the perception of NAMs: 1) professional profile (educational cohort, employer), 2) internal science communication within professional forums, 3) concern for "error cost," 4) collaboration across stakeholders, and 5) fundamental beliefs regarding toxicology. We conducted an online survey (n = 171; 2018) asking participants about their experiences and perspectives at events of the Society of Environmental Toxicology and Chemistry (SETAC) to assess 1) how NAMs are discussed compared to conventional testing and 2) how respondents perceive their viability. We developed ordered logistic regression (OLR) models to understand the influence of exploratory variables (cohort, core views on toxicology, frequency of collaboration) on respondents' evaluation of the viability of different NAMs. Our results showed that 1) NAMs were more likely than conventional methods to be challenged in forum discussions, which may be fueled by concerns for error costs in regulatory decision making; 2) perceptions of the viability of NAMs tended to follow a "pattern of familiarity," whereby respondents that were more knowledgeable about a test method tended to find it more viable; 3) respondents who agreed with the Paracelsus maxim had a greater likelihood of finding conventional testing viable; and 4) the more a respondent reported collaborating with industry on alternative testing strategies, the more likely she or he was to report that NAMs were less viable. These results suggest that there are professional and organizational barriers to greater acceptance of NAMs that can be addressed through a social learning process within the professional community. Integr Environ Assess Manag 2020;16:269-281. © 2020 SETAC.

Narrowing the Gap Between In Vitro and In Vivo Genetic Profiles by Deconvoluting Toxicogenomic Data In Silico

Toxicogenomics (TGx) is a powerful method to evaluate toxicity and is widely used in both in vivo and in vitro assays. For in vivo TGx, reduction, refinement, and replacement represent the unremitting pursuit of live-animal tests, but in vitro assays, as alternatives, usually demonstrate poor correlation with real in vivo assays. In living subjects, in addition to drug effects, inner-environmental reactions also affect genetic variation, and these two factors are further jointly reflected in gene abundance. Thus, finding a strategy to factorize inner-environmental factor from in vivo assays based on gene expression levels and to further utilize in vitro data to better simulate in vivo data is needed. We proposed a strategy based on post-modified non-negative matrix factorization, which can estimate the gene expression profiles and contents of major factors in samples. The applicability of the strategy was first verified, and the strategy was then utilized to simulate in vivo data by correcting in vitro data. The similarities between real in vivo data and simulated data (single-dose 0.72, repeat-doses 0.75) were higher than those observed when directly comparing real in vivo data with in vitro data (single-dose 0.56, repeat-doses 0.70). Moreover, by keeping environment-related factor, a simulation can always be generated by using in vitro data to provide potential substitutions for in vivo TGx and to reduce the launch of live-animal tests.

Assessment of Drugs Toxicity and Associated Biomarker Genes Using Hierarchical Clustering

Background and objectives: Assessment of drugs toxicity and associated biomarker genes is one of the most important tasks in the pre-clinical phase of drug development pipeline as well as in toxicogenomic studies. There are few statistical methods for the assessment of doses of drugs (DDs) toxicity and their associated biomarker genes. However, these methods consume more time for computation of the model parameters using the EM (expectation-maximization) based iterative approaches. To overcome this problem, in this paper, an attempt is made to propose an alternative approach based on hierarchical clustering (HC) for the same purpose. Methods and materials: There are several types of HC approaches whose performance depends on different similarity/distance measures. Therefore, we explored suitable combinations of distance measures and HC methods based on Japanese Toxicogenomics Project (TGP) datasets for better clustering/co-clustering between DDs and genes as well as to detect toxic DDs and their associated biomarker genes. Results: We observed that Word’s HC method with each of Euclidean, Manhattan, and Minkowski distance measures produces better clustering/co-clustering results. For an example, in the case of the glutathione metabolism pathway (GMP) dataset LOC100359539/Rrm2, Gpx6, RGD1562107, Gstm4, Gstm3, G6pd, Gsta5, Gclc, Mgst2, Gsr, Gpx2, Gclm, Gstp1, LOC100912604/Srm, Gstm4, Odc1, Gsr, Gss are the biomarker genes and Acetaminophen_Middle, Acetaminophen_High, Methapyrilene_High, Nitrofurazone_High, Nitrofurazone_Middle, Isoniazid_Middle, Isoniazid_High are their regulatory (associated) DDs explored by our proposed co-clustering algorithm based on the distance and HC method combination Euclidean: Word. Similarly, for the peroxisome proliferator-activated receptor signaling pathway (PPAR-SP) dataset Cpt1a, Cyp8b1, Cyp4a3, Ehhadh, Plin5, Plin2, Fabp3, Me1, Fabp5, LOC100910385, Cpt2, Acaa1a, Cyp4a1, LOC100365047, Cpt1a, LOC100365047, Angptl4, Aqp7, Cpt1c, Cpt1b, Me1 are the biomarker genes and Aspirin_Low, Aspirin_Middle, Aspirin_High, Benzbromarone_Middle, Benzbromarone_High, Clofibrate_Middle, Clofibrate_High, WY14643_Low, WY14643_High, WY14643_Middle, Gemfibrozil_Middle, Gemfibrozil_High are their regulatory DDs. Conclusions: Overall, the methods proposed in this article, co-cluster the genes and DDs as well as detect biomarker genes and their regulatory DDs simultaneously consuming less time compared to other mentioned methods. The results produced by the proposed methods have been validated by the available literature and functional annotation.

Variation in transcriptional responses to copper exposure across Daphnia pulex lineages

Copper pollution is pervasive in aquatic habitats and is particularly harmful to invertebrates sensitive to environmental changes such as Daphnia pulex. Mechanisms of toxicity and tolerance to copper are not well understood. We used RNA-sequencing to investigate these mechanisms in three genetically distinct D. pulex clonal lineages with different histories of copper exposure. Upregulated genes after copper exposure were enriched with Gene Ontology (GO) categories involved in digestion, molting and growth, whereas downregulated genes after copper exposure were enriched in the metal-regulatory system, immune response and epigenetic modifications. The three D. pulex clones in our study show largely similar transcriptional patterns in response to copper, with only a total of twenty genes differentially expressed in a single clonal lineages. We also detected lower relative expression of some genes known to be important for copper tolerance, metallothionein and glutathione-S-transferase, in a sensitive lineage sampled from an uncontaminated habitat. Daphnia-specific genes (without orthologs outside the genus) and Daphnia-specific duplications (genes duplicated in the Daphnia lineage) were overrepresented in differentially expressed genes, highlighting an important role for newly emerged genes in tolerating environmental stressors. The results indicate that the D. pulex lineages tested in this study generally respond to copper stress using the same major pathways, but that the more resistant clone with previous copper exposure might be better able to regulate key genes. This finding highlights the important nuances in gene expression among clones, shaped by historical exposure and influencing copper tolerance.

384-Channel electrochemical sensor array chips based on hybridization-triggered switching for simultaneous oligonucleotide detection

We investigated the feasibility of simultaneous detection of multiple environmentally- and biomedically-relevant RNA biomarker target sequences on a single newly fabricated 384-ch sensor array chip aiming at practical application. The individual sensor is composed of a photolithographically-fabricated Au/Cr-based electrode modified with peptide nucleic acid (PNA) probes. The sensor array chips showed sequence-specific responses upon hybridization of the probes with target sequences complementary to the probes in contrast to mismatch versions. The target oligonucleotides have 15-22 mer sequences from messenger RNAs for estrogen-responsive genes and microRNAs for lung cancer biomarkers. The dependence on target concentrations of sensor responses was observed by using a single chip on which experiments for detection of several target concentrations proceeded simultaneously, with the detection limit of 7.33 × 10^-8 M. As more realistic samples, oligonucleotide samples amplified by PCR from a synthesized template sequence were applied to the chip. They showed sequence-specific responses, revealing the potential for fabricated sensor array chips to be utilized to analyze PCR samples. Unlike complicated and expensive chips that require nanofabrication, our sensor array chips based on glass coated with gold thin films are simple and can be fabricated from inexpensive and readily available materials.

Expression of immune, antioxidant and stress related genes in different organs of common carp exposed to indoxacarb

The aim of the present study was to investigate the effects of chronic exposure of common carp (Cyprinus carpio) to indoxacarb on immune, antioxidant and stress gene expression. After 21 days exposure to 0, 0.75, 1.5 and 3 ppm indoxacarb, expression of IL-1β, IL-8, IL-10, TNF-α, IFN-γ, SOD, CAT, HSP70, IGF-I and IGF-II were assessed in liver, kidney and gills. In general, exposure to low concentration of indoxacarb increased inflammatory cytokine gene expression (IL-1β, IL-8, IL-10, TNF-α and IFN-γ) and inhibits inflammatory cytokines' expression at higher concentrations. The assessment of antioxidant gene expression (SOD and CAT) in different organs indicate that they were increased by low concentrations of indoxacarb to deal with primary oxidative situation. However, higher concentrations of indoxacarb caused reduction in oxidative gene expression. IGF genes expression in liver significantly increased at a concentration of 0.75 ppm treatment, then it decreased at 1.5 ppm indoxacarb and increased again by increasing in the indoxacarb concentration to 3 ppm. The expression of HSP70 in kidney showed a significant elevation in 0.75 and 1.5 ppm treatments compared with 3 ppm treatment and the control group. The expression of this gene in liver was significantly increased in 1.5 and 3 ppm treatments. The same pattern of expression was also observed in gill. Overall, indoxacarb exposure affects common carp health at transcription levels. Changes in the genes expression generally suggest that indoxacarb exposure led to interference in inflammation, oxidative stress and tissue damage.

GWAS and Beyond: Using Omics Approaches to Interpret SNP Associations

PURPOSE OF REVIEW

Neurodegenerative diseases, neuropsychiatric disorders, and related traits have highly complex etiologies but are also highly heritable and identifying the causal genes and biological pathways underlying these traits may advance the development of treatments and preventive strategies. While many genome-wide association studies (GWAS) have successfully identified variants contributing to polygenic neurodegenerative and neuropsychiatric phenotypes including Alzheimer's disease (AD), schizophrenia (SCZ), and bipolar disorder (BPD) amongst others, interpreting the biological roles of significantly-associated variants in the genetic architecture of these traits remains a significant challenge. Here we review several 'omics' approaches which attempt to bridge the gap from associated genetic variants to phenotype by helping define the functional roles of GWAS loci in the development of neuropsychiatric disorders and traits.

RECENT FINDINGS

Several common 'omics' approaches have been applied to examine neuropsychiatric traits, such as nearest-gene mapping, trans-ethnic fine mapping, annotation enrichment analysis, transcriptomic analysis, and pathway analysis, and each of these approaches has strengths and limitations in providing insight into biological mechanisms. One popular emerging method is the examination of tissue-specific genetically-regulated gene expression (GReX), which aggregates the genetic variants' effects at the gene-level. Furthermore, proteomic, metabolomic, and microbiomic studies and phenome-wide association studies will further enhance our understanding of neuropsychiatric traits.

SUMMARY

GWAS has been applied to neuropsychiatric traits for a decade, but our understanding about the biological function of identified variants remains limited. Today, technological advancements have created analytical approaches for integrating transcriptomics, metabolomics, proteomics, pharmacology and toxicology as tools for understanding the functional roles of genetics variants. These data, as well as the broader clinical information provided by electronic health records, can provide additional insight and complement genomic analyses.

Differences among Unique Nanoparticle Protein Corona Constructs: A Case Study Using Data Analytics and Multi-Variant Visualization to Describe Physicochemical Characteristics

Gold nanoparticles (AuNPs) used in pharmaceutical treatments have been shown to effectively deliver a payload, such as an active pharmaceutical ingredient or image contrast agent, to targeted tissues in need of therapy or diagnostics while minimizing exposure, availability, and accumulation to surrounding biological compartments. Data sets collected in this field of study include some toxico- and pharmacodynamic properties (e.g., distribution and metabolism) but many studies lack information about adsorption of biological molecules or absorption into cells. When nanoparticles are suspended in blood serum, a protein corona cloud forms around its surface. The extent of the applications and implications of this formed cloud are unknown. Some researchers have speculated that the successful use of nanoparticles in pharmaceutical treatments relies on a comprehensive understanding of the protein corona composition. The work presented in this paper uses a suite of data analytics and multi-variant visualization techniques to elucidate particle-to-protein interactions at the molecular level. Through mass spectrometry analyses, corona proteins were identified through large and complex datasets. With such high-output analyses, complex datasets pose a challenge when visualizing and communicating nanoparticle-protein interactions. Thus, the creation of a streamlined visualization method is necessary. A series of user-friendly data informatics techniques were used to demonstrate the data flow of protein corona characteristics. Multi-variant heat maps, pie charts, tables, and three-dimensional regression analyses were used to improve results interpretation, facilitate an iterative data transfer process, and emphasize features of the nanoparticle-protein corona system that might be controllable. Data informatics successfully highlights the differences between protein corona compositions and how they relate to nanoparticle surface charge.

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.

Gene expression profiling toward the next generation safety control of influenza vaccines and adjuvants in Japan

Influenza becomes epidemic worldwide every year, and many individuals receive vaccination annually. Quality control relating to safety and potency of influenza vaccines is important to maintain public confidence. The safety of influenza vaccines has been assessed by clinical trials, and animal safety tests are performed to monitor the consistent quality between vaccines used for clinical trials and marketing; the biological responses in vaccinated animals are evaluated, including changes in body weight and white blood cell count. Animal safety tests have been contributing to the quality relating to the safety of influenza vaccines for decades, but improvements are needed. Although precise mechanisms involving biological changes in animal safety tests have not been fully elucidated, the application of cDNA microarray technology make it possible to reliably identify genes related to biological responses in vaccinated animals. From analysis of the expression profile of >10,000 genes of lung in animals treated with an inactivated whole virion influenza vaccine, we identified 17 marker genes whose expression patterns correlated well to changes in body weight and leukocyte count in vaccinated animals. In influenza HA vaccine-treated animals exhibiting subtle changes in biological responses, a robust expression pattern of marker genes was found. Furthermore, these marker genes could also be used in the evaluation of adjuvanted influenza vaccines. The expression profile of marker genes is expected to be an alternative indicator for safety control of various influenza vaccines conferring high sensitivity and short turnaround time. Thus, gene expression profiling may be a powerful tool for safety control of vaccines in the future.

Transcriptomic, lipid, and histological profiles suggest changes in health in fish from a pesticide hot spot

Barramundi (Lates calcarifer) were collected at the beginning (1st sampling) and end (2nd sampling) of the wet season from Sandy Creek, an agriculturally impacted catchment in the Mackay Whitsundays region of the Great Barrier Reef catchment area, and from Repulse Creek, located approximately 100 km north in Conway National Park, to assess the impacts of pesticide exposure. Gill and liver histology, lipid class composition in muscle, and the hepatic transcriptome were examined. The first sample of Repulse Creek fish showed little tissue damage and low transcript levels of xenobiotic metabolism enzymes. Sandy Creek fish showed altered transcriptomic patterns, including those that regulate lipid metabolism, xenobiotic metabolism, and immune response; gross histological alterations including lipidosis; and differences in some lipid classes. The second sampling of Repulse Creek fish showed similar alterations in hepatic transcriptome and tissue structure as fish from Sandy Creek. These changes may indicate a decrease in health of pesticide exposed fish.

XCMS-MRM and METLIN-MRM: a cloud library and public resource for targeted analysis of small molecules

We report XCMS-MRM and METLIN-MRM ( http://xcmsonline-mrm.scripps.edu/ and http://metlin.scripps.edu/ ), a cloud-based data-analysis platform and a public multiple-reaction monitoring (MRM) transition repository for small-molecule quantitative tandem mass spectrometry. This platform provides MRM transitions for more than 15,500 molecules and facilitates data sharing across different instruments and laboratories.

Halobenzoquinone-Induced Alteration of Gene Expression Associated with Oxidative Stress Signaling Pathways

Halobenzoquinones (HBQs) are emerging disinfection byproducts (DBPs) that effectively induce reactive oxygen species and oxidative damage in vitro. However, the impacts of HBQs on oxidative-stress-related gene expression have not been investigated. In this study, we examined alterations in the expression of 44 genes related to oxidative-stress-induced signaling pathways in human uroepithelial cells (SV-HUC-1) upon exposure to six HBQs. The results show the structure-dependent effects of HBQs on the studied gene expression. After 2 h of exposure, the expression levels of 9 to 28 genes were altered, while after 8 h of exposure, the expression levels of 29 to 31 genes were altered. Four genes ( HMOX1, NQO1, PTGS2, and TXNRD1) were significantly upregulated by all six HBQs at both exposure time points. Ingenuity pathway analysis revealed that the Nrf2 pathway was significantly responsive to HBQ exposure. Other canonical pathways responsive to HBQ exposure included GSH redox reductions, superoxide radical degradation, and xenobiotic metabolism signaling. This study has demonstrated that HBQs significantly alter the gene expression of oxidative-stress-related signaling pathways and contributes to the understanding of HBQ-DBP-associated toxicity.

Transcriptome analysis to assess the cholestatic hepatotoxicity induced by Polygoni Multiflori Radix: Up-regulation of key enzymes of cholesterol and bile acid biosynthesis

Polygoni Multiflori Radix (PMR) has been commonly used as a tonic in China for centuries. However, PMR-associated hepatotoxicity is becoming a safety issue. Cholestasis often occurs in PMR-induced hepatotoxicity in clinical medicine, but the exact mechanism is not completely understood. An RNA-Seq method was employed, in the present study, to explore the molecular mechanism of cholestatic liver injury induced by PMR, characterized by the hepatic transcriptional response in rats exposed to 1 and 20 g/kg PMR for 90 days. Pathological changes seen in rat livers exposed to PMR included increased bile ducts in portal areas and biliary epithelial cell hyperplasia, which were accompanied by the elevation of serum biochemistries. Dose-dependent increases in the expression of 14 transcripts encoding enzymes involved in the cholesterol biosynthetic pathway were identified. Furthermore, cholesterol 7-alpha hydroxylase (Cyp7a1), a rate-limiting enzyme in the synthesis of bile acids (BAs) from cholesterol, was found to be upregulated by PMR treatment. Protein analysis by western blot suggested that expression of 3-hydroxy-3-methylglutaryl CoA reductase (Hmgcr) and Cyp7a1 were increased in a dose-dependent manner. Collectively, the present study demonstrates that PMR upregulates key enzymes for biosynthesis of cholesterol and BA, which poses the risk of cholestatic liver injury.

SIGNIFICANCE

To the best of our knowledge, this is the first transcriptome analysis to highlight the main molecular changes occurring in rats chronic exposed to PMR. We have identified 39 specific differentially expressed genes (DEGs) that were present in various comparisons. A total of 14 of these altered gene transcripts were associated with cholesterol biosynthesis. Another factor of great importance in our opinion seemed to be the enhancement of bile acid (BA) biosynthesis, which were closely linked to cholesterol biosynthesis or metabolism. Our findings suggested that the disturbance on balance of BA formation and elimination might lead to a BA overload in hepatocytes, thereby resulting in liver injury.