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

Integrating Multi-Omics with environmental data for precision health: A novel analytic framework and case study on prenatal mercury induced childhood fatty liver disease

BACKGROUND

Precision Health aims to revolutionize disease prevention by leveraging information across multiple omic datasets (multi-omics). However, existing methods generally do not consider personalized environmental risk factors (e.g., environmental pollutants).

OBJECTIVE

To develop and apply a precision health framework which combines multiomic integration (including early, intermediate, and late integration, representing sequential stages at which omics layers are combined for modeling) with mediation approaches (including high-dimensional mediation to identify biomarkers, mediation with latent factors to identify pathways, and integrated/quasi-mediation to identify high-risk subpopulations) to identify novel biomarkers of prenatal mercury induced metabolic dysfunction-associated fatty liver disease (MAFLD), elucidate molecular pathways linking prenatal mercury with MAFLD in children, and identify high-risk children based on integrated exposure and multiomics data.

METHODS

This prospective cohort study used data from 420 mother-child pairs from the Human Early Life Exposome (HELIX) project. Mercury concentrations were determined in maternal or cord blood from pregnancy. Cytokeratin 18 (CK-18; a MAFLD biomarker) and five omics layers (DNA Methylation, gene transcription, microRNA, proteins, and metabolites) were measured in blood in childhood (age 6-10 years).

RESULTS

Each standard deviation increase in prenatal mercury was associated with a 0.11 [95% confidence interval: 0.02-0.21] standard deviation increase in CK-18. High dimensional mediation analysis identified 10 biomarkers linking prenatal mercury and CK-18, including six CpG sites and four transcripts. Mediation with latent factors identified molecular pathways linking mercury and MAFLD, including altered cytokine signaling and hepatic stellate cell activation. Integrated/quasi-mediation identified high risk subgroups of children based on unique combinations of exposure levels, omics profiles (driven by epigenetic markers), and MAFLD.

CONCLUSIONS

Prenatal mercury exposure is associated with elevated liver enzymes in childhood, likely through alterations in DNA methylation and gene expression. Our analytic framework can be applied across many different fields and serve as a resource to help guide future precision health investigations.

Changes in transcriptome regulations of a marine rotifer Brachionus plicatilis under methylmercury stress

Mercury (Hg), a heavy metal pollutant worldwide, can be transformed into methylmercury (MeHg) by various aquatic microorganisms in water, thus accumulating along the aquatic food chain and posing a particular challenge to human health. Zooplankton plays a crucial role in aquatic ecosystems and serves as a major component of the food chain. To evaluate the effects of MeHg on the rotifer Brachionus plicatilis and reveal the underlying mechanism of these effects, we exposed B. plicatilis to MeHg by either direct immersion or by feeding with MeHg-poisoned Chlorella pyrenoidesa, respectively, and conducted a transcriptomic analysis. The results showed that B. plicatilis directly exposed to MeHg by immersion showed significant enrichment of the glutathione metabolism pathway for detoxification of MeHg. In addition, the exposure to MeHg by feeding induced a significant enrichment of lysosome and notch signaling pathways of rotifers, supporting the hypothesis that MeHg can induce autophagy dysfunction in cells and disturb the nervous system of rotifers. In two different routes of MeHg exposure, the pathway of cytochrome P450 in rotifers showed significant enrichment for resisting MeHg toxicity. Our results suggest further studies on the potential mechanism and biological responses of MeHg toxicity in other links of the aquatic food chain.

Identification of Potential Biomarkers for Group I Pulmonary Hypertension Based on Machine Learning and Bioinformatics Analysis

A number of processes and pathways have been reported in the development of Group I pulmonary hypertension (Group I PAH); however, novel biomarkers need to be identified for a better diagnosis and management. We employed a robust rank aggregation (RRA) algorithm to shortlist the key differentially expressed genes (DEGs) between Group I PAH patients and controls. An optimal diagnostic model was obtained by comparing seven machine learning algorithms and was verified in an independent dataset. The functional roles of key DEGs and biomarkers were analyzed using various in silico methods. Finally, the biomarkers and a set of key candidates were experimentally validated using patient samples and a cell line model. A total of 48 key DEGs with preferable diagnostic value were identified. A gradient boosting decision tree algorithm was utilized to build a diagnostic model with three biomarkers, PBRM1, CA1, and TXLNG. An immune-cell infiltration analysis revealed significant differences in the relative abundances of seven immune cells between controls and PAH patients and a correlation with the biomarkers. Experimental validation confirmed the upregulation of the three biomarkers in Group I PAH patients. In conclusion, machine learning and a bioinformatics analysis along with experimental techniques identified PBRM1, CA1, and TXLNG as potential biomarkers for Group I PAH.

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.

High Accuracy Classification of Developmental Toxicants by In Vitro Tests of Human Neuroepithelial and Cardiomyoblast Differentiation

Human-relevant tests to predict developmental toxicity are urgently needed. A currently intensively studied approach makes use of differentiating human stem cells to measure chemically-induced deviations of the normal developmental program, as in a recent study based on cardiac differentiation (UKK2). Here, we (i) tested the performance of an assay modeling neuroepithelial differentiation (UKN1), and (ii) explored the benefit of combining assays (UKN1 and UKK2) that model different germ layers. Substance-induced cytotoxicity and genome-wide expression profiles of 23 teratogens and 16 non-teratogens at human-relevant concentrations were generated and used for statistical classification, resulting in accuracies of the UKN1 assay of 87-90%. A comparison to the UKK2 assay (accuracies of 90-92%) showed, in general, a high congruence in compound classification that may be explained by the fact that there was a high overlap of signaling pathways. Finally, the combination of both assays improved the prediction compared to each test alone, and reached accuracies of 92-95%. Although some compounds were misclassified by the individual tests, we conclude that UKN1 and UKK2 can be used for a reliable detection of teratogens in vitro, and that a combined analysis of tests that differentiate hiPSCs into different germ layers and cell types can even further improve the prediction of developmental toxicants.

Methylmercury-induced DNA methylation—From epidemiological observations to experimental evidence

Methylmercury (MeHg) is a developmental neurotoxicant, and one potential mechanism of MeHg toxicity is epigenetic dysregulation. In a recent meta-analysis of epigenome-wide association studies (EWAS), associations between prenatal MeHg exposure and DNA methylation at several genomic sites were identified in blood from newborns and children. While EWASs reveal human-relevant associations, experimental studies are required to validate the relationship between exposure and DNA methylation changes, and to assess if such changes have implications for gene expression. Herein, we studied DNA methylation and gene expression of five of the top genes identified in the EWAS meta-analysis, MED31, MRPL19, GGH, GRK1, and LYSMD3, upon MeHg exposure in human SH-SY5Y cells exposed to 8 or 40 nM of MeHg during differentiation, using bisulfite-pyrosequencing and qPCR, respectively. The concentrations were selected to cover the range of MeHg concentrations in cord blood (2–8.5 μg/L) observed in the cohorts included in the EWAS. Exposure to MeHg increased DNA methylation at MED31, a transcriptional regulator essential for fetal development. The results were in concordance with the epidemiological findings where more MED31 methylation was associated with higher concentrations of MeHg. Additionally, we found a non-significant decrease in DNA methylation at GGH, which corresponds to the direction of change observed in the EWAS, and a significant correlation of GGH methylation with its expression. In conclusion, this study corroborates some of the EWAS findings and puts forward candidate genes involved in MeHg’s effects on the developing brain, thus highlighting the value of experimental validation of epidemiological association studies.

Classification of Developmental Toxicants in a Human iPSC Transcriptomics-Based Test

Despite the progress made in developmental toxicology, there is a great need for in vitro tests that identify developmental toxicants in relation to human oral doses and blood concentrations. In the present study, we established the hiPSC-based UKK2 in vitro test and analyzed genome-wide expression profiles of 23 known teratogens and 16 non-teratogens. Compounds were analyzed at the maximal plasma concentration (Cmax) and at 20-fold Cmax for a 24 h incubation period in three independent experiments. Based on the 1000 probe sets with the highest variance and including information on cytotoxicity, penalized logistic regression with leave-one-out cross-validation was used to classify the compounds as test-positive or test-negative, reaching an area under the curve (AUC), accuracy, sensitivity, and specificity of 0.96, 0.92, 0.96, and 0.88, respectively. Omitting the cytotoxicity information reduced the test performance to an AUC of 0.94, an accuracy of 0.79, and a sensitivity of 0.74. A second method, which used the number of significantly deregulated probe sets to classify the compounds, resulted in a specificity of 1; however, the AUC (0.90), accuracy (0.90), and sensitivity (0.83) were inferior compared to those of the logistic regression-based procedure. Finally, no increased performance was achieved when the high test concentrations (20-fold Cmax) were used, in comparison to testing within the realistic clinical range (1-fold Cmax). In conclusion, although further optimization is required, for example, by including additional readouts and cell systems that model different developmental processes, the UKK2-test in its present form can support the early discovery-phase detection of human developmental toxicants.

Optimization of the TeraTox assay for preclinical teratogenicity assessment

Current animal-free methods to assess teratogenicity of drugs under development still deliver high numbers of false negatives. To improve the sensitivity of human teratogenicity prediction, we characterized the TeraTox test, a newly developed multilineage differentiation assay using 3D human-induced pluripotent stem cells. TeraTox produces primary output concentration-dependent cytotoxicity and altered gene expression induced by each test compound. These data are fed into an interpretable machine-learning model to perform prediction, which relates to the concentration-dependent human teratogenicity potential of drug candidates. We applied TeraTox to profile 33 approved pharmaceuticals and 12 proprietary drug candidates with known in vivo data. Comparing TeraTox predictions with known human or animal toxicity, we report an accuracy of 69% (specificity: 53%, sensitivity: 79%). TeraTox performed better than 2 quantitative structure-activity relationship models and had a higher sensitivity than the murine embryonic stem cell test (accuracy: 58%, specificity: 76%, and sensitivity: 46%) run in the same laboratory. The overall prediction accuracy could be further improved by combining TeraTox and mouse embryonic stem cell test results. Furthermore, patterns of altered gene expression revealed by TeraTox may help grouping toxicologically similar compounds and possibly deducing common modes of action. The TeraTox assay and the dataset described here therefore represent a new tool and a valuable resource for drug teratogenicity assessment.

DNA methylation changes associated with prenatal mercury exposure: A meta-analysis of prospective cohort studies from PACE consortium

Mercury (Hg) is a ubiquitous heavy metal that originates from both natural and anthropogenic sources and is transformed in the environment to its most toxicant form, methylmercury (MeHg). Recent studies suggest that MeHg exposure can alter epigenetic modifications during embryogenesis. In this study, we examined associations between prenatal MeHg exposure and levels of cord blood DNA methylation (DNAm) by meta-analysis in up to seven independent studies (n = 1462) as well as persistence of those relationships in blood from 7 to 8 year-old children (n = 794). In cord blood, we found limited evidence of differential DNAm at cg24184221 in MED31 (β = 2.28 × 10-4, p-value = 5.87 × 10-5) in relation to prenatal MeHg exposure. In child blood, we identified differential DNAm at cg15288800 (β = 0.004, p-value = 4.97 × 10-5), also located in MED31. This repeated link to MED31, a gene involved in lipid metabolism and RNA Polymerase II transcription function, may suggest a DNAm perturbation related to MeHg exposure that persists into early childhood. Further, we found evidence for association between prenatal MeHg exposure and child blood DNAm levels at two additional CpGs: cg12204245 (β = 0.002, p-value = 4.81 × 10-7) in GRK1 and cg02212000 (β = -0.001, p-value = 8.13 × 10-7) in GGH. Prenatal MeHg exposure was associated with DNAm modifications that may influence health outcomes, such as cognitive or anthropometric development, in different populations.

OUP accepted manuscript

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

RNA modifications as emerging therapeutic targets

The field of epitranscriptome, posttranscriptional modifications to RNAs, is still growing up and has presented substantial evidences for the role of RNA modifications in diseases. In terms of new drug development, RNA modifications have a great promise for therapy. For example, more than 170 type of modifications exist in various types of RNAs. Regulatory genes and their roles in critical biological process have been identified and they are associated with several diseases. Current data, for example, identification of inhibitors targeting RNA modifications regulatory genes, strongly support the idea that RNA modifications have potential as emerging therapeutic targets. Therefore, in this review, RNA modifications and regulatory genes were comprehensively documented in terms of drug development by summarizing the findings from previous studies. It was discussed how RNA modifications or regulatory genes can be targeted by altering molecular mechanisms. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA Processing > RNA Editing and Modification.

Non-Lethal Concentration of MeHg Causes Marked Responses in the DNA Repair, Integrity, and Replication Pathways in the Exposed Human Salivary Gland Cell Line

In Brazilian northern Amazon, communities are potentially exposed and vulnerable to methylmercury (MeHg) toxicity through the vast ingestion of fish. In vivo and in vitro studies demonstrated that the salivary glands as a susceptible organ to this potent environmental pollutant, reporting alterations on physiological, biochemical, and proteomic parameters. However, the alterations caused by MeHg on the gene expression of the exposed human salivary gland cells are still unknown. Therefore, the goal was to perform the transcriptome profile of the human salivary gland cell line after exposure to MeHg, using the microarray technique and posterior bioinformatics analysis. The cell exposure was performed using 2.5 µM MeHg. A previously published study demonstrated that this concentration belongs to a range of concentrations that caused biochemical and metabolic alterations in this linage. As a result, the MeHg exposure did not cause lethality in the human salivary gland cells line but was able to alter the expression of 155 genes. Downregulated genes (15) are entirety relating to the cell metabolism impairment, and according to KEGG analysis, they belong to the glycosphingolipid (GSL) biosynthesis pathway. On the other hand, most of the 140 upregulated genes were related to cell-cycle progression, DNA repair, and replication pathway, or cellular defenses through the GSH basal metabolism. These genomic changes revealed the effort to the cell to maintain physiological and genomic stability to avoid cell death, being in accordance with the nonlethality in the toxicity test. Last, the results support in-depth studies on nonlethal MeHg concentrations for biomarkers identification that interpret transcriptomics data in toxicological tests serving as an early alert of physiological changes in vitro biological models.

Bioaccumulation of mercury and transcriptional responses in tusk (Brosme brosme), a deep-water fish from a Norwegian fjord

High concentrations of mercury (Hg) have been documented in deep-water fish species from some Norwegian fjords. In this study, tusk (Brosme brosme) was sampled from four locations in the innermost parts of Sognefjorden in Western Norway. Total Hg and methylmercury (MeHg) levels were measured in liver tissue. To search for potential sublethal effects of Hg, we characterized the hepatic transcriptome in tusk with high and low levels of Hg bioaccumulation using global transcriptomics analysis (RNA-seq). The results showed that there was a significant correlation between fish weight and accumulated concentrations of MeHg but not total Hg. MeHg accounted for 30-40% of total Hg in liver of most of the fish, although at concentrations above 2-3 mg Hg/kg wet weight the percentage of MeHg dropped considerably. Transcriptome analysis resulted in hundreds of differentially expressed genes in the liver of tusk with high Hg levels. Functional enrichment analysis suggested that the top affected pathways are associated with protein folding, adipogenesis, notch signaling, and lipid metabolism (beta-oxidation and phospholipids). Based on transcriptional responses pointing to well-known effects of Hg compounds in fish, the study suggests that tusk in Sognefjorden could be negatively impacted by Hg bioaccumulation.

Transcriptional Analyses of Acute Exposure to Methylmercury on Erythrocytes of Loggerhead Sea Turtle

To understand changes in enzyme activity and gene expression as biomarkers of exposure to methylmercury, we exposed loggerhead turtle erythrocytes (RBCs) to concentrations of 0, 1, and 5 mg L-1 of MeHg and de novo transcriptome were assembled using RNA-seq. The analysis of differentially expressed genes (DEGs) indicated that 79 unique genes were dysregulated (39 upregulated and 44 downregulated genes). The results showed that MeHg altered gene expression patterns as a response to the cellular stress produced, reflected in cell cycle regulation, lysosomal activity, autophagy, calcium regulation, mitochondrial regulation, apoptosis, and regulation of transcription and translation. The analysis of DEGs showed a low response of the antioxidant machinery to MeHg, evidenced by the fact that genes of early response to oxidative stress were not dysregulated. The RBCs maintained a constitutive expression of proteins that represented a good part of the defense against reactive oxygen species (ROS) induced by MeHg.

Kinetic modeling of stem cell transcriptome dynamics to identify regulatory modules of normal and disturbed neuroectodermal differentiation

Thousands of transcriptome data sets are available, but approaches for their use in dynamic cell response modelling are few, especially for processes affected simultaneously by two orthogonal influencing variables. We approached this problem for neuroepithelial development of human pluripotent stem cells (differentiation variable), in the presence or absence of valproic acid (signaling variable). Using few basic assumptions (sequential differentiation states of cells; discrete on/off states for individual genes in these states), and time-resolved transcriptome data, a comprehensive model of spontaneous and perturbed gene expression dynamics was developed. The model made reliable predictions (average correlation of 0.85 between predicted and subsequently tested expression values). Even regulations predicted to be non-monotonic were successfully validated by PCR in new sets of experiments. Transient patterns of gene regulation were identified from model predictions. They pointed towards activation of Wnt signaling as a candidate pathway leading to a redirection of differentiation away from neuroepithelial cells towards neural crest. Intervention experiments, using a Wnt/beta-catenin antagonist, led to a phenotypic rescue of this disturbed differentiation. Thus, our broadly applicable model allows the analysis of transcriptome changes in complex time/perturbation matrices.

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

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

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

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

Accuracy, discriminative properties and reliability of a human ESC-based in vitro toxicity assay to distinguish teratogens responsible for neural tube defects

The poor correlation of developmental toxicity studies in animals with human outcome data has emphasized the need for complementary assays based on human cells and tissues. As neural tube defects represent an important proportion of congenital malformations, we evaluated here the accuracy of a human embryonic stem cell (hESC)-based assay to predict chemically induced disruption of neural tube formation. As teratogenic compounds, we used cyclopamine (CPA), valproic acid (VPA), ochratoxin A (OTA) and mycophenolic acid (MMF), all suspected or known inducers of human neural tube defects, as well as theophylline and saccharin as negative control compounds. We analyzed their effects on the ability of hES cells to give rise to neural precursors (expressing specific marker Nestin), to form neural tube-like structures (rosettes), and to express specific markers (Sox1, Otx2, Lix1, EvI1, Rspo3) during rosette formation. The results showed that various effects of the selected compounds on early neural development could be specifically revealed in vitro through related alterations of neurogenic differentiation of hESC. Furthermore, it was possible to discriminate toxicants acting at different time points during embryonic development and, therefore, responsible for distinct adverse effects on neural tube formation. By comparing four different hESC lines, we observed a significant (up to fivefold) variability of the line-dependent response to toxicants. We highlight at least two sources of variability: one related to the heterogeneity of hESC lines in culture (stemness/commitment profiles); the second to possible genetically determined differences in individual sensitivity to teratogens.

Association between methylmercury environmental exposure and neurological disorders: A systematic review

The mercury-related central nervous system disorders have been extensively studied on animal models and human beings. However, clinical evidences of which neurological changes are in fact associated with mercury exposure remains controversial. This systematic review (Prospero registration under the number CRD42016041760) aimed to elucidate the association of methylmercury (MeHg) exposure with neurological alteration in populations living in MeHg-endemic risk area. A systematic search was performed according to the Preferred Reporting Items for Systematic Review and Meta-Analysis criteria using available databases PubMed, LILACS, Scopus, Web of Science, The Cochrane Library, OpenGrey and Google Scholar. A search of the following terms: "methylmercury compounds", "organomercury compounds", "neurologic manifestations", "memory disorders", "neurobehavioral manifestations" and "communication disorders" were performed in a systematic way. Studies focusing on MeHg exposure and subsequent neurological alteration on humans (>13 years) were included. Evaluation of methodological quality and risk of bias as well as the level of evidence was performed. Our results have identified 470 studies and six articles were eligible for systematic review inclusion criteria. The studies suggested alterations related to the psychosensory, motor and coordination system, as well as motor speech, hearing, visual impairment, mood alterations and loss of intelligent quotient. Of all the six studies, two presented a high risk of bias, with methodological problems related to the confounding factors and all studies presented evidence level ranged from very low to low. In this way our results revealed that a definitive demonstration of an association of MeHg and neurological alterations in human beings is still a pending subject. Future studies in this topic should take into consideration more confident and reliable methods to answer this question.

Memory effect of arsenic-induced cellular response and its influences on toxicity of titanium dioxide nanoparticle

Toxicity of arsenic (As) has been widely characterized. However, few studies focus on whether cell responses induced by As at nontoxic concentration could be inherited and further change cell tolerance to another pollutant. In this study, human A549 and HeLa cells were exposed to As at nontoxic concentrations for 10 or 15 passages, then the cells were recovered in the cell medium without As. At 25th passage, residual As in both type of cells was completely removed through the recovery process. And no abnormity in cell viability was identified in both type of cells between control and As-treated groups. Above results indicated that As exposure-recovery treatment had limited influence on phenotype of the cells. However, gene expression profiles determined by high-throughput sequencing showed that As exposure-recovery treatment induced similar expression modification of genes related to inflammation, oxidative stress and epigenetic modulation in the A549 and HeLa cells after recovery of 10 or 15 passages, indicating that As-induced cellular responses have been partially memorized at transcriptional level. The memory effect might play key roles in increased tolerance of the A549 and HeLa cells to adverse effects (cell viability, intracellular reactive oxygen species (ROS) generation and plasma membrane damage) induced by titanium dioxide nanoparticles (as representative pollutant). This study shed new lights on toxic effects induced by As at nontoxic concentration, which is useful for risk assessment of combined effects of As and other pollutants.

Look-alikes may not act alike: Gene expression regulation and cell-type-specific responses of three valproic acid analogues in the neural embryonic stem cell test (ESTn)

In vitro assays to assess developmental neurotoxicity of chemicals are highly desirable. The murine neural embryonic stem cell test (ESTn) can mimic parts of early differentiation of embryonic brain and may therefore be useful for this purpose. The aim of this study was to investigate whether this test is able to rank the toxic potencies of three valproic acid analogues and to study their mode of action by investigating their individual effects on four cell types: stem cells, neurons, astrocytes and neural crest cells. Using immunocytochemical read-outs and qPCR for cell type-specific genes, the effects of valproic acid (VPA), 2-ethylhexanoic acid (EHA) and 2-ethyl-4-methylpentanoic (EMPA) were assessed. VPA and EHA but not EMPA downregulated cell type-specific differentiation makers and upregulated stem cell related markers (Fut4, Cdh1) at different time points during differentiation. Expression of Gfap, a marker for astrocytes, was dramatically downregulated by VPA and EHA, but not by EMPA. This finding was verified using immunostainings. Based on the number and extent of genes regulated by the three compounds, relative potencies were determined as VPA > EHA > EMPA, which is consistent with in vivo developmental toxicity potency ranking of these compounds. Thus, ESTn using a combination of morphology, gene and protein expression readouts, may provide a medium-throughput system for monitoring the effects of compounds on differentiation of cell types in early brain development.

Recommendation on test readiness criteria for new approach methods in toxicology: Exemplified for developmental neurotoxicity

Multiple non-animal-based test methods have never been formally validated. In order to use such new approach methods (NAMs) in a regulatory context, criteria to define their readiness are necessary. The field of developmental neurotoxicity (DNT) testing is used to exemplify the application of readiness criteria. The costs and number of untested chemicals are overwhelming for in vivo DNT testing. Thus, there is a need for inexpensive, high-throughput NAMs, to obtain initial information on potential hazards, and to allow prioritization for further testing. A background on the regulatory and scientific status of DNT testing is provided showing different types of test readiness levels, depending on the intended use of data from NAMs. Readiness criteria, compiled during a stakeholder workshop, uniting scientists from academia, industry and regulatory authorities are presented. An important step beyond the listing of criteria, was the suggestion for a preliminary scoring scheme. On this basis a (semi)-quantitative analysis process was assembled on test readiness of 17 NAMs with respect to various uses (e.g. prioritization/screening, risk assessment). The scoring results suggest that several assays are currently at high readiness levels. Therefore, suggestions are made on how DNT NAMs may be assembled into an integrated approach to testing and assessment (IATA). In parallel, the testing state in these assays was compiled for more than 1000 compounds. Finally, a vision is presented on how further NAM development may be guided by knowledge of signaling pathways necessary for brain development, DNT pathophysiology, and relevant adverse outcome pathways (AOP).

Stage-specific metabolic features of differentiating neurons: Implications for toxicant sensitivity

Developmental neurotoxicity (DNT) may be induced when chemicals disturb a key neurodevelopmental process, and many tests focus on this type of toxicity. Alternatively, DNT may occur when chemicals are cytotoxic only during a specific neurodevelopmental stage. The toxicant sensitivity is affected by the expression of toxicant targets and by resilience factors. Although cellular metabolism plays an important role, little is known how it changes during human neurogenesis, and how potential alterations affect toxicant sensitivity of mature vs. immature neurons. We used immature (d0) and mature (d6) LUHMES cells (dopaminergic human neurons) to provide initial answers to these questions. Transcriptome profiling and characterization of energy metabolism suggested a switch from predominantly glycolytic energy generation to a more pronounced contribution of the tricarboxylic acid cycle (TCA) during neuronal maturation. Therefore, we used pulsed stable isotope-resolved metabolomics (pSIRM) to determine intracellular metabolite pool sizes (concentrations), and isotopically non-stationary ¹³C-metabolic flux analysis (INST ¹³C-MFA) to calculate metabolic fluxes. We found that d0 cells mainly use glutamine to fuel the TCA. Furthermore, they rely on extracellular pyruvate to allow continuous growth. This metabolic situation does not allow for mitochondrial or glycolytic spare capacity, i.e. the ability to adapt energy generation to altered needs. Accordingly, neuronal precursor cells displayed a higher sensitivity to several mitochondrial toxicants than mature neurons differentiated from them. In summary, this study shows that precursor cells lose their glutamine dependency during differentiation while they gain flexibility of energy generation and thereby increase their resistance to low concentrations of mitochondrial toxicants.

New animal-free concepts and test methods for developmental toxicity and peripheral neurotoxicity

The complex toxicological fields of repeat dose organ toxicity (RDT) and developmental and reproductive toxicity (DART) still require new concepts and approaches to achieve a fully animal-free safety assessment of chemicals. One novel approach is the generation of relevant human cell types from pluripotent stem cells, and the use of such cells for the establishment of phenotypic test methods. Due to their broad endpoints, such tests capture multiple types of toxicants, i.e. they are a readout for the activation of many adverse outcome pathways (AOPs). The 2016 Lush Science Prize was awarded for the development of one such assay, the PeriTox test, which uses human peripheral neurons generated from stem cells. The assay endpoints measure various cell functions, and these give information on the potential neurotoxicity and developmental neurotoxicity hazard of test compounds. The PeriTox test method has a high predictivity and sensitivity for peripheral neurotoxicants, and thus addresses the inherent challenges in pesticide testing and drug development. Data from the test can be obtained quickly and at a relatively high-throughput, and thus, the assay has the potential to replace animal-based safety assessment during early product development or for screening potential environmental toxicants.

1,2-Dichloroethane Induces Reproductive Toxicity Mediated by the CREM/CREB Signaling Pathway in Male NIH Swiss Mice

1,2-Dichloroethane (1,2-DCE) is a widely used chlorinated organic toxicant but little is known about the reproductive disorders induced by its excessive exposure. To reveal 1,2-DCE-induced male reproductive toxicity and to elucidate the underlying mechanisms, we exposed male National Institutes of Health Swiss mice to 1,2-DCE by inhalation at 0, 100, 350, and 700 mg/m3 for 6 h/day, for 1 and 4 weeks. Our findings showed a significant decrease in body weight with increased testis/body weight ratio, reduced sperm concentration and induced malformation of spermatozoa, and vacuolar degeneration of germ cells in the seminiferous tubules of testes in mice exposed to 1,2-DCE. Cyclic adenosine monophosphate (cAMP)-response element binding protein (CREB) and cAMP-response element modulator (CREM) were significantly inhibited by 1,2-DCE. This is consistent with the declines in the transducer of regulated CREB activity 1 and activator of CREM in testis, which results in the decrease in lactate dehydrogenase C and testis-specific kinase 1 in the testes. Moreover, the activation of p53 and Bax with the inhibition of Bcl-2 might be the reason for the upregulation of caspase-3 in the apoptosis, as detected by TdT-mediated dUTP nick-end labeling assay in the testes induced by 1,2-DCE. Finally, elevated testosterone levels were found along with increased levels of gonadotropin-releasing hormone, cAMP, luteinizing hormone (LH), and LH receptors in the testes. These findings suggest that 1,2-DCE inhibits CREM/CREB signaling cascade and subsequently induces apoptosis associated with p53 activation and mitochondrial dysfunction. This also results in induced malformation of spermatozoa, reduced sperm concentration, and pathological impairment of the testes.

Oxidative Biochemistry Disbalance and Changes on Proteomic Profile in Salivary Glands of Rats Induced by Chronic Exposure to Methylmercury

Methylmercury (MeHg) is one of the most toxic mercury species, which can cause many systemic damages, but little is known about its effect in the salivary glands. This study aimed to analyze the mercury levels, oxidative stress, and proteomic profile in parotid, submandibular, and sublingual salivary glands of rats, after chronic MeHg intoxication. Two groups of twenty male Wistar rats (90 days of age) were used on the experiment. MeHg group was intoxicated by intragastric gavage with MeHg at a dose of 0.04 mg/kg/day for 60 days, while the control group received only oil. After the period of intoxication, the glands were collected for evaluation of total mercury levels, proteomic profile, and oxidative balance by analyzing the antioxidant capacity against peroxyl radicals (ACAP), lipid peroxidation (LPO), and nitrite levels. Our results have showed that mercury levels were significant in all three glands compared to the respective control. It also showed lower levels of ACAP, as well as higher LPO and nitrite levels. The proteomic profile presented impairments on structural components of cytoskeleton, metabolic pathways, and oxidative biochemistry. Thus, the exposure to MeHg was able to generate oxidative stress that could be associated with changes in the proteomic profile of salivary glands.