Lineage-specific regulation of epigenetic modifier genes in human liver and brain

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.

hMSCs-derived exosome circCDK13 inhibits liver fibrosis by regulating the expression of MFGE8 through miR-17-5p/KAT2B

OBJECTIVE

To investigate the effects of human bone marrow mesenchymal stem cells (hMSCs)-derived exosome circCDK13 on liver fibrosis and its mechanism.

METHODS

Exosomes derived from hMSCs were extracted and identified by flow cytometry and osteogenic and adipogenic induction, and the expressions of marker proteins on the surface of exosomes were detected by western blot. Cell proliferation was measured by CCK8 assay, the expression of active markers of HSCs by immunofluorescence, and the expressions of fibrosis-related factors by western blot. A mouse model of liver fibrosis was established by intraperitoneal injection of thioacetamide (TAA). Fibrosis was detected by HE staining, Masson staining, and Sirius red staining. Western blot was utilized to test the expressions of PI3K/AKT and NF-κB pathway related proteins, dual-luciferase reporter assay and RIP assay to validate the binding between circCDK13 and miR-17-5p as well as between miR-17-5p and KAT2B, and ChIP to validate the effect of KAT2B on H3 acetylation and MFGE8 transcription.

RESULTS

hMSCs-derived exosomes inhibited liver fibrosis mainly through circCDK13. Dual-luciferase reporter assay and RIP assay demonstrated the binding between circCDK13 and miR-17-5p as well as between miR-17-5p and KAT2B. Further experimental results indicated that circCDK13 mediated liver fibrosis by regulating the miR-17-5p/KAT2B axis, and KAT2B promoted MFGE8 transcription by H3 acetylation. Exo-circCDK13 inhibited PI3K/AKT and NF-κB signaling pathways activation through regulating the miR-17-5p/KAT2B axis.

CONCLUSION

hMSCs-derived exosome circCDK13 inhibited liver fibrosis by regulating the expression of MFGE8 through miR-17-5p/KAT2B axis.

Building a schizophrenia genetic network: transcription factor 4 regulates genes involved in neuronal development and schizophrenia risk

The transcription factor 4 (TCF4) locus is a robust association finding with schizophrenia (SCZ), but little is known about the genes regulated by the encoded transcription factor. Therefore, we conducted chromatin immunoprecipitation sequencing (ChIP-seq) of TCF4 in neural-derived (SH-SY5Y) cells to identify genome-wide TCF4 binding sites, followed by data integration with SCZ association findings. We identified 11 322 TCF4 binding sites overlapping in two ChIP-seq experiments. These sites are significantly enriched for the TCF4 Ebox binding motif (>85% having ≥1 Ebox) and implicate a gene set enriched for genes downregulated in TCF4 small-interfering RNA (siRNA) knockdown experiments, indicating the validity of our findings. The TCF4 gene set was also enriched among (1) gene ontology categories such as axon/neuronal development, (2) genes preferentially expressed in brain, in particular pyramidal neurons of the somatosensory cortex and (3) genes downregulated in postmortem brain tissue from SCZ patients (odds ratio, OR = 2.8, permutation P < 4x10-5). Considering genomic alignments, TCF4 binding sites significantly overlapped those for neural DNA-binding proteins such as FOXP2 and the SCZ-associated EP300. TCF4 binding sites were modestly enriched among SCZ risk loci from the Psychiatric Genomic Consortium (OR = 1.56, P = 0.03). In total, 130 TCF4 binding sites occurred in 39 of the 108 regions published in 2014. Thirteen genes within the 108 loci had both a TCF4 binding site ±10kb and were differentially expressed in siRNA knockdown experiments of TCF4, suggesting direct TCF4 regulation. These findings confirm TCF4 as an important regulator of neural genes and point toward functional interactions with potential relevance for SCZ.

Epigenetic Modifications of the Liver Tumor Cell Line HepG2 Increase Their Drug Metabolic Capacity

Although human liver tumor cells have reduced metabolic functions as compared to primary human hepatocytes (PHH) they are widely used for pre-screening tests of drug metabolism and toxicity. The aim of the present study was to modify liver cancer cell lines in order to improve their drug-metabolizing activities towards PHH. It is well-known that epigenetics is strongly modified in tumor cells and that epigenetic regulators influence the expression and function of Cytochrome P450 (CYP) enzymes through altering crucial transcription factors responsible for drug-metabolizing enzymes. Therefore, we screened the epigenetic status of four different liver cancer cell lines (Huh7, HLE, HepG2 and AKN-1) which were reported to have metabolizing drug activities. Our results showed that HepG2 cells demonstrated the highest similarity compared to PHH. Thus, we modified the epigenetic status of HepG2 cells towards 'normal' liver cells by 5-Azacytidine (5-AZA) and Vitamin C exposure. Then, mRNA expression of Epithelial-mesenchymal transition (EMT) marker SNAIL and CYP enzymes were measured by PCR and determinate specific drug metabolites, associated with CYP enzymes by LC/MS. Our results demonstrated an epigenetic shift in HepG2 cells towards PHH after exposure to 5-AZA and Vitamin C which resulted in a higher expression and activity of specific drug metabolizing CYP enzymes. Finally, we observed that 5-AZA and Vitamin C led to an increased expression of Hepatocyte nuclear factor 4α (HNF4α) and E-Cadherin and a significant down regulation of Snail1 (SNAIL), the key transcriptional repressor of E-Cadherin. Our study shows, that certain phase I genes and their enzyme activities are increased by epigenetic modification in HepG2 cells with a concomitant reduction of EMT marker gene SNAIL. The enhancing of liver specific functions in hepatoma cells using epigenetic modifiers opens new opportunities for the usage of cell lines as a potential liver in vitro model for drug testing and development.

EZH1 is an antipsychotic-sensitive epigenetic modulator of social and motivational behavior that is dysregulated in schizophrenia

BACKGROUND

With the capacity to modulate gene networks in an environmentally-sensitive manner, the role of epigenetic systems in mental disorders has come under intense investigation. Dysregulation of epigenetic effectors, including microRNAs and histone-modifying enzymes, may better explain the role of environmental risk factors and the observed heritability rate that cannot be fully attributed to known genetic risk alleles. Here, we aimed to identify novel epigenetic targets of the schizophrenia-associated microRNA 132 (miR-132).

METHODS

Histone modifications were quantified by immunodetection in response to viral-mediated overexpression of miR-132 while a luminescent reporter system was used to validate targets of miR-132 in vitro. Genome-wide profiling, quantitative PCR and NanoSting were used to quantify gene expression in post-mortem human brains, neuronal cultures and prefrontal cortex (PFC) of mice chronically exposed to antipsychotics. Following viral-mediated depletion of Enhancer of Zeste 1 (EZH1) in the murine PFC, behaviors including sociability and motivation were assessed using a 3-chambered apparatus and forced-swim test, respectively.

RESULTS

Overexpression of miR-132 decreased global histone 3 lysine 27 tri-methylation (H3K27me3), a repressive epigenetic mark. Moreover, the polycomb-associated H3K27 methyltransferase, EZH1, is regulated by miR-132 and upregulated in the PFC of schizophrenics. Unlike its homolog EZH2, expression of EZH1 in the murine PFC decreased following chronic exposure to antipsychotics. Viral-mediated depletion of EZH1 in the mouse PFC attenuated sociability, enhanced motivational behaviors, and affected gene expression pathways related to neurotransmission and behavioral phenotypes.

CONCLUSIONS

EZH1 is dysregulated in schizophrenia, sensitive to antipsychotic medications, and a brain-enriched miR-132 target that controls neurobehavioral phenotypes.

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).

Reverse-transcription quantitative PCR directly from cells without RNA extraction and without isothermal reverse-transcription: a 'zero-step' RT-qPCR protocol

We describe an ultra-rapid and sensitive method to quantify gene expression levels in cultured cells. The procedure is based on reverse-transcription quantitative PCR (RT-qPCR) directly from cells, without RNA extraction and without an isothermal reverse-transcription step. Human neurons (Lund human mesencephalic cells) were lysed at different stages of differentiation, and the lysates were used directly as template for the combined RT-qPCR reaction. We detected a down-regulation of a proliferation marker and an up-regulation of neuronal dopaminergic genes expression. We were able to detect the reference gene target from as few as a single cell, demonstrating the application of the method for efficient amplification from small cell numbers. The data were fully in line with those obtained by the standard two-step RT-qPCR from the extracted total RNA. Our ‘zero-step’ RT-qPCR method proved to be simple and reliable with a total time from cell lysis to the end of the qPCR as short as 1.5 h. It is therefore particularly suitable for RT-qPCRs where large numbers of samples must be handled, or where data are required within short time.

Impairment of human neural crest cell migration by prolonged exposure to interferon-beta

Human cell-based toxicological assays have been used successfully to detect known toxicants, and to distinguish them from negative controls. However, there is at present little experience on how to deal with hits from screens of compounds with yet unknown hazard. As a case study to this issue, we characterized human interferon-beta (IFNβ) as potential developmental toxicant affecting neural crest cells (NCC). The protein was identified as a hit during a screen of clinically used drugs in the ‘migration inhibition of neural crest’ (MINC) assay. Concentration–response studies in the MINC combined with immunocytochemistry and mRNA quantification of cellular markers showed that IFNβ inhibited NCC migration at concentrations as low as 20 pM. The effective concentrations found here correspond to levels found in human plasma, and they were neither cytostatic nor cytotoxic nor did they did they affect the differentiation state and overall phenotype of NCC. Data from two other migration assays confirmed that picomolar concentration of IFNβ reduced the motility of NCC, while other interferons were less potent. The activation of JAK kinase by IFNβ, as suggested by bioinformatics analysis of the transcriptome changes, was confirmed by biochemical methods. The degree and duration of pathway activation correlated with the extent of migration inhibition, and pharmacological block of this signaling pathway before, or up to 6 h after exposure to the cytokine prevented the effects of IFNβ on migration. Thus, the reduction of vital functions of human NCC is a hitherto unknown potential hazard of endogenous or pharmacologically applied interferons.

Definition of transcriptome-based indices for quantitative characterization of chemically disturbed stem cell development: introduction of the STOP-Toxukn and STOP-Toxukk tests

Stem cell-based in vitro test systems can recapitulate specific phases of human development. In the UKK test system, human pluripotent stem cells (hPSCs) randomly differentiate into cells of the three germ layers and their derivatives. In the UKN1 test system, hPSCs differentiate into early neural precursor cells. During the normal differentiation period (14 days) of the UKK system, 570 genes [849 probe sets (PSs)] were regulated >fivefold; in the UKN1 system (6 days), 879 genes (1238 PSs) were regulated. We refer to these genes as 'developmental genes'. In the present study, we used genome-wide expression data of 12 test substances in the UKK and UKN1 test systems to understand the basic principles of how chemicals interfere with the spontaneous transcriptional development in both test systems. The set of test compounds included six histone deacetylase inhibitors (HDACis), six mercury-containing compounds ('mercurials') and thalidomide. All compounds were tested at the maximum non-cytotoxic concentration, while valproic acid and thalidomide were additionally tested over a wide range of concentrations. In total, 242 genes (252 PSs) in the UKK test system and 793 genes (1092 PSs) in the UKN1 test system were deregulated by the 12 test compounds. We identified sets of 'diagnostic genes' appropriate for the identification of the influence of HDACis or mercurials. Test compounds that interfered with the expression of developmental genes usually antagonized their spontaneous development, meaning that up-regulated developmental genes were suppressed and developmental genes whose expression normally decreases were induced. The fraction of compromised developmental genes varied widely between the test compounds, and it reached up to 60 %. To quantitatively describe disturbed development on a genome-wide basis, we recommend a concept of two indices, 'developmental potency' (D _p) and 'developmental index' (D _i), whereby D _p is the fraction of all developmental genes that are up- or down-regulated by a test compound, and D _i is the ratio of overrepresentation of developmental genes among all genes deregulated by a test compound. The use of D _i makes hazard identification more sensitive because some compounds compromise the expression of only a relatively small number of genes but have a high propensity to deregulate developmental genes specifically, resulting in a low D _p but a high D _i. In conclusion, the concept based on the indices D _p and D _i offers the possibility to quantitatively express the propensity of test compounds to interfere with normal development.

A transcriptome-based classifier to identify developmental toxicants by stem cell testing: design, validation and optimization for histone deacetylase inhibitors

Test systems to identify developmental toxicants are urgently needed. A combination of human stem cell technology and transcriptome analysis was to provide a proof of concept that toxicants with a related mode of action can be identified and grouped for read-across. We chose a test system of developmental toxicity, related to the generation of neuroectoderm from pluripotent stem cells (UKN1), and exposed cells for 6 days to the histone deacetylase inhibitors (HDACi) valproic acid, trichostatin A, vorinostat, belinostat, panobinostat and entinostat. To provide insight into their toxic action, we identified HDACi consensus genes, assigned them to superordinate biological processes and mapped them to a human transcription factor network constructed from hundreds of transcriptome data sets. We also tested a heterogeneous group of ‘mercurials’ (methylmercury, thimerosal, mercury(II)chloride, mercury(II)bromide, 4-chloromercuribenzoic acid, phenylmercuric acid). Microarray data were compared at the highest non-cytotoxic concentration for all 12 toxicants. A support vector machine (SVM)-based classifier predicted all HDACi correctly. For validation, the classifier was applied to legacy data sets of HDACi, and for each exposure situation, the SVM predictions correlated with the developmental toxicity. Finally, optimization of the classifier based on 100 probe sets showed that eight genes (F2RL2, TFAP2B, EDNRA, FOXD3, SIX3, MT1E, ETS1 and LHX2) are sufficient to separate HDACi from mercurials. Our data demonstrate how human stem cells and transcriptome analysis can be combined for mechanistic grouping and prediction of toxicants. Extension of this concept to mechanisms beyond HDACi would allow prediction of human developmental toxicity hazard of unknown compounds with the UKN1 test system.

Differences in the Epigenetic Regulation of Cytochrome P450 Genes between Human Embryonic Stem Cell-Derived Hepatocytes and Primary Hepatocytes

Human pluripotent stem cell-derived hepatocytes have the potential to provide in vitro model systems for drug discovery and hepatotoxicity testing. However, these cells are currently unsuitable for drug toxicity and efficacy testing because of their limited expression of genes encoding drug-metabolizing enzymes, especially cytochrome P450 (CYP) enzymes. Transcript levels of major CYP genes were much lower in human embryonic stem cell-derived hepatocytes (hESC-Hep) than in human primary hepatocytes (hPH). To verify the mechanism underlying this reduced expression of CYP genes, including CYP1A1, CYP1A2, CYP1B1, CYP2D6, and CYP2E1, we investigated their epigenetic regulation in terms of DNA methylation and histone modifications in hESC-Hep and hPH. CpG islands of CYP genes were hypermethylated in hESC-Hep, whereas they had an open chromatin structure, as represented by hypomethylation of CpG sites and permissive histone modifications, in hPH. Inhibition of DNA methyltransferases (DNMTs) during hepatic maturation induced demethylation of the CpG sites of CYP1A1 and CYP1A2, leading to the up-regulation of their transcription. Combinatorial inhibition of DNMTs and histone deacetylases (HDACs) increased the transcript levels of CYP1A1, CYP1A2, CYP1B1, and CYP2D6. Our findings suggest that limited expression of CYP genes in hESC-Hep is modulated by epigenetic regulatory factors such as DNMTs and HDACs.

Evolution of an epigenetic gene ensemble within the genus Anopheles

Epigenetic control of gene expression has important implications for the regulation of developmental processes, for mediating homeostasis and responses to the external environment, and for transgenerational inheritance of gene expression patterns. Genes that mediate epigenetic control have been well-characterized in Drosophila melanogaster, and we have identified and analyzed an orthologous gene ensemble in Anopheles gambiae that comprises 169 orthologs related to a 215-member epigenetic gene ensemble in D. melanogaster. We find that this ensemble is highly conserved among anopheline mosquitoes, as we identify only seven gene family expansion/contraction events within the ensemble among 12 mosquito species we have studied within the genus Anopheles. Comparative analyses of the epigenetic gene expression across the genera Drosophila and Anopheles reveal distinct tissue-associated expression patterns in the two genera, but similar temporal expression patterns. The A. gambiae complex and D. melanogaster subgroup epigenetic gene ensembles exhibit similar evolutionary rates, as assessed by their respective dN/dS values. These differences in tissue-associated expression patterns, in contrast to similarities in evolutionary rates and temporal expression patterns, may imply that some members of the epigenetic gene ensemble have been redeployed within one or both genera, in comparison to the most recent common ancestor of these two clades. Members of this epigenetic gene ensemble may constitute another set of potential targets for vector control and enable further reductions in the burden of human malaria, by analogy to recent success in development of small molecule antagonists for mammalian epigenetic machinery.