Integrated Omic Analyses Identify Pathways and Transcriptomic Regulators Associated With Chemical Alterations of In Vitro Neural Network Formation.

Development of in vitro new approach methodologies has been driven by the need for developmental neurotoxicity (DNT) hazard data on thousands of chemicals. The network formation assay characterizes DNT hazard based on changes in network formation but provides no mechanistic information. This study investigated nervous system signaling pathways and upstream physiological regulators underlying chemically induced neural network dysfunction. Rat primary cortical neural networks grown on microelectrode arrays were exposed for 12 days in vitro to cytosine arabinoside, 5-fluorouracil, domoic acid, cypermethrin, deltamethrin, or haloperidol as these exposures altered network formation in previous studies. RNA-seq from cells and gas chromatography/mass spectrometry analysis of media extracts collected on days in vitro 12 provided gene expression and metabolomic identification, respectively. The integration of differentially expressed genes and metabolites for each neurotoxicant was analyzed using ingenuity pathway analysis. All 6 compounds altered gene expression that linked to developmental disorders and neurological diseases. Other enriched canonical pathways overlapped among compounds of the same class; eg, genes and metabolites altered by both cytosine arabinoside and 5-fluorouracil exposures are enriched in axonal guidance pathways. Integrated analysis of upstream regulators was heterogeneous across compounds, but identified several transcriptomic regulators including CREB1, SOX2, NOTCH1, and PRODH. These results demonstrate that changes in network formation are accompanied by transcriptomic and metabolomic changes and that different classes of compounds produce differing responses. This approach can enhance information obtained from new approach methodologies and contribute to the identification and development of adverse outcome pathways associated with DNT.

Applications of Functional Genomics for Drug Discovery.

Many diseases, such as diabetes, autoimmune diseases, cancer, and neurological disorders, are caused by a dysregulation of a complex interplay of genes. Genome-wide association studies have identified thousands of disease-linked polymorphisms in the human population. However, detailing the causative gene expression or functional changes underlying those associations has been elusive in many cases. Functional genomics is an emerging field of research that aims to deconvolute the link between genotype and phenotype by making use of large -omic data sets and next-generation gene and epigenome editing tools to perturb genes of interest. Here we review how functional genomic tools can be used to better understand the biological interplay between genes, improve disease modeling, and identify novel drug targets. Incorporation of functional genomic capabilities into conventional drug development pipelines is predicted to expedite the development of first-in-class therapeutics.

Defining toxicological tipping points in neuronal network development.

Measuring electrical activity of neural networks by microelectrode array (MEA) has recently shown promise for screening level assessments of chemical toxicity on network development and function. Important aspects of interneuronal communication can be quantified from a single MEA recording, including individual firing rates, coordinated bursting, and measures of network synchrony, providing rich datasets to evaluate chemical effects. Further, multiple recordings can be made from the same network, including during the formation of these networks in vitro. The ability to perform multiple recording sessions over the in vitro development of network activity may provide further insight into developmental effects of neurotoxicants. In the current study, a recently described MEA-based screen of 86 compounds in primary rat cortical cultures over 12 days in vitro was revisited to establish a framework that integrates all available primary measures of electrical activity from MEA recordings into a composite metric for deviation from normal activity (total scalar perturbation). Examining scalar perturbations over time and increasing concentration of compound allowed for definition of critical concentrations or "tipping points" at which the neural networks switched from recovery to non-recovery trajectories for 42 compounds. These tipping point concentrations occurred at predominantly lower concentrations than those causing overt cell viability loss or disrupting individual network parameters, suggesting tipping points may be a more sensitive measure of network functional loss. Comparing tipping points for six compounds with plasma concentrations known to cause developmental neurotoxicity in vivo demonstrated strong concordance and suggests there is potential for using tipping points for chemical prioritization.

From the Cover: Developmental Neurotoxicants Disrupt Activity in Cortical Networks on Microelectrode Arrays: Results of Screening 86 Compounds During Neural Network Formation.

Less than 1% of environmental chemicals have been evaluated for developmental neurotoxicity (DNT). Current guideline DNT studies are resource intensive and not amenable to screening large numbers of compounds for hazard. As part of evaluating a battery of more rapid and scalable in vitro assays for DNT hazard, 86 compounds were screened for their ability to alter function during cortical network development. Developing rat cortical networks were treated with a concentration series (usually 0.03-30 µM) of 86 compounds, 60 of which have known in vivo DNT effects ("DNT Reference Set"). Spontaneous network activity was monitored by microelectrode array recordings over 12 days in vitro, and 17 measures of network activity and synchrony were quantified. Following recordings on days in vitro 12, in-well cell assessment of metabolic activity (Alamar blue) and total cellular content (lactase dehydrogenase) were conducted. Of the 86 compounds tested, 64 perturbed cortical network function in a concentration-dependent manner; 49 of the 60 DNT Reference Set compounds (81.7%) altered network formation. Compounds were ranked by potency (network effect EC50) and selectivity (separation of network and cell viability EC50) for hazard prioritization. Machine learning indicates a combination of an overall network activity metric with a measure of network coordination is key in distinguishing network-disruptive from benign treatments. These data demonstrate that this microelectrode array-based assay for developing cortical network function is amenable to medium-throughput evaluation of environmental substances for DNT hazard and further prioritization. For comprehensive identification of compounds of concern, this assay will be a useful component of a battery of assays targeting independent neurodevelopmental processes.

Editor's Highlight: Evaluation of a Microelectrode Array-Based Assay for Neural Network Ontogeny Using Training Set Chemicals.

Thousands of compounds in the environment have not been characterized for developmental neurotoxicity (DNT) hazard. To address this issue, methods to screen compounds rapidly for DNT hazard evaluation are necessary and are being developed for key neurodevelopmental processes. In order to develop an assay for network formation, this study evaluated effects of a training set of chemicals on network ontogeny by measuring spontaneous electrical activity in neural networks grown on microelectrode arrays (MEAs). Rat (0-24 h old) primary cortical cells were plated in 48 well-MEA plates and exposed to 6 compounds: acetaminophen, bisindolylmaleimide-1 (Bis-1), domoic acid, mevastatin, sodium orthovanadate, and loperamide for a period of 12 days. Spontaneous network activity was recorded on days 2, 5, 7, 9, and 12 and viability was assessed using the Cell Titer Blue assay on day 12. Network activity (e.g. mean firing rate [MFR], burst rate [BR], etc), increased between days 5 and 12. Random Forest analysis indicated that across all compounds and times, temporal correlation of firing patterns (r), MFR, BR, number of active electrodes and % of spikes in a burst were the most influential parameters in separating control from treated wells. All compounds except acetaminophen (≤ 30 µM) caused concentration-related effects on one or more of these parameters. Domoic acid and sodium orthovanadate altered several of these parameters in the absence of cytotoxicity. Although cytotoxicity was observed with Bis1, mevastatin, and loperamide, some parameters were affected by these compounds at concentrations below those resulting in cytotoxicity. These results demonstrate that this assay may be suitable for screening of compounds for DNT hazard identification.

HDAC inhibitors cause site-specific chromatin remodeling at PU.1-bound enhancers in K562 cells.

BACKGROUND: Small molecule inhibitors of histone deacetylases (HDACi) hold promise as anticancer agents for particular malignancies. However, clinical use is often confounded by toxicity, perhaps due to indiscriminate hyperacetylation of cellular proteins. Therefore, elucidating the mechanisms by which HDACi trigger differentiation, cell cycle arrest, or apoptosis of cancer cells could inform development of more targeted therapies. We used the myelogenous leukemia line K562 as a model of HDACi-induced differentiation to investigate chromatin accessibility (DNase-seq) and expression (RNA-seq) changes associated with this process. RESULTS: We identified several thousand specific regulatory elements [~10 % of total DNase I-hypersensitive (DHS) sites] that become significantly more or less accessible with sodium butyrate or suberanilohydroxamic acid treatment. Most of the differential DHS sites display hallmarks of enhancers, including being enriched for non-promoter regions, associating with nearby gene expression changes, and increasing luciferase reporter expression in K562 cells. Differential DHS sites were enriched for key hematopoietic lineage transcription factor motifs, including SPI1 (PU.1), a known pioneer factor. We found PU.1 increases binding at opened DHS sites with HDACi treatment by ChIP-seq, but PU.1 knockdown by shRNA fails to block the chromatin accessibility and expression changes. A machine-learning approach indicates H3K27me3 initially marks PU.1-bound sites that open with HDACi treatment, suggesting these sites are epigenetically poised. CONCLUSIONS: We find HDACi treatment of K562 cells results in site-specific chromatin remodeling at epigenetically poised regulatory elements. PU.1 shows evidence of a pioneer role in this process by marking poised enhancers but is not required for transcriptional activation.

Regulation of chromatin accessibility and Zic binding at enhancers in the developing cerebellum.

To identify chromatin mechanisms of neuronal differentiation, we characterized chromatin accessibility and gene expression in cerebellar granule neurons (CGNs) of the developing mouse. We used DNase-seq to map accessibility of cis-regulatory elements and RNA-seq to profile transcript abundance across postnatal stages of neuronal differentiation in vivo and in culture. We observed thousands of chromatin accessibility changes as CGNs differentiated, and verified, using H3K27ac ChIP-seq, reporter gene assays and CRISPR-mediated activation, that many of these regions function as neuronal enhancers. Motif discovery in differentially accessible chromatin regions suggested a previously unknown role for the Zic family of transcription factors in CGN maturation. We confirmed the association of Zic with these elements by ChIP-seq and found, using knockdown, that Zic1 and Zic2 are required for coordinating mature neuronal gene expression patterns. Together, our data reveal chromatin dynamics at thousands of gene regulatory elements that facilitate the gene expression patterns necessary for neuronal differentiation and function.

Genome accessibility is widely preserved and locally modulated during mitosis.

Mitosis entails global alterations to chromosome structure and nuclear architecture, concomitant with transient silencing of transcription. How cells transmit transcriptional states through mitosis remains incompletely understood. While many nuclear factors dissociate from mitotic chromosomes, the observation that certain nuclear factors and chromatin features remain associated with individual loci during mitosis originated the hypothesis that such mitotically retained molecular signatures could provide transcriptional memory through mitosis. To understand the role of chromatin structure in mitotic memory, we performed the first genome-wide comparison of DNase I sensitivity of chromatin in mitosis and interphase, using a murine erythroblast model. Despite chromosome condensation during mitosis visible by microscopy, the landscape of chromatin accessibility at the macromolecular level is largely unaltered. However, mitotic chromatin accessibility is locally dynamic, with individual loci maintaining none, some, or all of their interphase accessibility. Mitotic reduction in accessibility occurs primarily within narrow, highly DNase hypersensitive sites that frequently coincide with transcription factor binding sites, whereas broader domains of moderate accessibility tend to be more stable. In mitosis, proximal promoters generally maintain their accessibility more strongly, whereas distal regulatory elements tend to lose accessibility. Large domains of DNA hypomethylation mark a subset of promoters that retain accessibility during mitosis and across many cell types in interphase. Erythroid transcription factor GATA1 exerts site-specific changes in interphase accessibility that are most pronounced at distal regulatory elements, but has little influence on mitotic accessibility. We conclude that features of open chromatin are remarkably stable through mitosis, but are modulated at the level of individual genes and regulatory elements.

Explicit DNase sequence bias modeling enables high-resolution transcription factor footprint detection.

DNaseI footprinting is an established assay for identifying transcription factor (TF)-DNA interactions with single base pair resolution. High-throughput DNase-seq assays have recently been used to detect in vivo DNase footprints across the genome. Multiple computational approaches have been developed to identify DNase-seq footprints as predictors of TF binding. However, recent studies have pointed to a substantial cleavage bias of DNase and its negative impact on predictive performance of footprinting. To assess the potential for using DNase-seq to identify individual binding sites, we performed DNase-seq on deproteinized genomic DNA and determined sequence cleavage bias. This allowed us to build bias corrected and TF-specific footprint models. The predictive performance of these models demonstrated that predicted footprints corresponded to high-confidence TF-DNA interactions. DNase-seq footprints were absent under a fraction of ChIP-seq peaks, which we show to be indicative of weaker binding, indirect TF-DNA interactions or possible ChIP artifacts. The modeling approach was also able to detect variation in the consensus motifs that TFs bind to. Finally, cell type specific footprints were detected within DNase hypersensitive sites that are present in multiple cell types, further supporting that footprints can identify changes in TF binding that are not detectable using other strategies.

The histone lysine demethylase Kdm6b is required for activity-dependent preconditioning of hippocampal neuronal survival.

Enzymes that regulate histone lysine methylation play important roles in neuronal differentiation, but little is known about their contributions to activity-regulated gene transcription in differentiated neurons. We characterized activity-regulated expression of lysine demethylases and lysine methyltransferases in the hippocampus of adult male mice following pilocarpine-induced seizure. Pilocarpine drove a 20-fold increase in mRNA encoding the histone H3 lysine 27-specific demethylase Kdm6b selectively in granule neurons of the dentate gyrus, and this induction was recapitulated in cultured hippocampal neurons by bicuculline and 4-aminopyridine (Bic + 4AP) stimulation of synaptic activity. Because activity-regulated gene expression is highly correlated with neuronal survival, we tested the requirement for Kdm6b expression in Bic + 4AP induced preconditioning of neuronal survival. Prior exposure to Bic + 4AP promoted neuronal survival in control neurons upon growth factor withdrawal; however, this effect was ablated when we knocked down Kdm6b expression. Loss of Kdm6b did not disrupt activity-induced expression of most genes, including that of a gene set previously established to promote neuronal survival in this assay. However, using bioinformatic analysis of RNA sequencing data, we discovered that Kdm6b knockdown neurons showed impaired inducibility of a discrete set of genes annotated for their function in inflammation. These data reveal a novel function for Kdm6b in activity-regulated neuronal survival, and they suggest that activity- and Kdm6b-dependent regulation of inflammatory gene pathways may serve as an adaptive pro-survival response to increased neuronal activity.

Microbiota modulate transcription in the intestinal epithelium without remodeling the accessible chromatin landscape.

Microbiota regulate intestinal physiology by modifying host gene expression along the length of the intestine, but the underlying regulatory mechanisms remain unresolved. Transcriptional specificity occurs through interactions between transcription factors (TFs) and cis-regulatory regions (CRRs) characterized by nucleosome-depleted accessible chromatin. We profiled transcriptome and accessible chromatin landscapes in intestinal epithelial cells (IECs) from mice reared in the presence or absence of microbiota. We show that regional differences in gene transcription along the intestinal tract were accompanied by major alterations in chromatin accessibility. Surprisingly, we discovered that microbiota modify host gene transcription in IECs without significantly impacting the accessible chromatin landscape. Instead, microbiota regulation of host gene transcription might be achieved by differential expression of specific TFs and enrichment of their binding sites in nucleosome-depleted CRRs near target genes. Our results suggest that the chromatin landscape in IECs is preprogrammed by the host in a region-specific manner to permit responses to microbiota through binding of open CRRs by specific TFs.

The viral and cellular microRNA targetome in lymphoblastoid cell lines.

Epstein-Barr virus (EBV) is a ubiquitous human herpesvirus linked to a number of B cell cancers and lymphoproliferative disorders. During latent infection, EBV expresses 25 viral pre-microRNAs (miRNAs) and induces the expression of specific host miRNAs, such as miR-155 and miR-21, which potentially play a role in viral oncogenesis. To date, only a limited number of EBV miRNA targets have been identified; thus, the role of EBV miRNAs in viral pathogenesis and/or lymphomagenesis is not well defined. Here, we used photoactivatable ribonucleoside-enhanced crosslinking and immunoprecipitation (PAR-CLIP) combined with deep sequencing and computational analysis to comprehensively examine the viral and cellular miRNA targetome in EBV strain B95-8-infected lymphoblastoid cell lines (LCLs). We identified 7,827 miRNA-interaction sites in 3,492 cellular 3'UTRs. 531 of these sites contained seed matches to viral miRNAs. 24 PAR-CLIP-identified miRNA:3'UTR interactions were confirmed by reporter assays. Our results reveal that EBV miRNAs predominantly target cellular transcripts during latent infection, thereby manipulating the host environment. Furthermore, targets of EBV miRNAs are involved in multiple cellular processes that are directly relevant to viral infection, including innate immunity, cell survival, and cell proliferation. Finally, we present evidence that myc-regulated host miRNAs from the miR-17/92 cluster can regulate latent viral gene expression. This comprehensive survey of the miRNA targetome in EBV-infected B cells represents a key step towards defining the functions of EBV-encoded miRNAs, and potentially, identifying novel therapeutic targets for EBV-associated malignancies.

Control of activating transcription factor 4 (ATF4) persistence by multisite phosphorylation impacts cell cycle progression and neurogenesis.

Organogenesis is a highly integrated process with a fundamental requirement for precise cell cycle control. Mechanistically, the cell cycle is composed of transitions and thresholds that are controlled by coordinated post-translational modifications. In this study, we describe a novel mechanism controlling the persistence of the transcription factor ATF4 by multisite phosphorylation. Proline-directed phosphorylation acted additively to regulate multiple aspects of ATF4 degradation. Stabilized ATF4 mutants exhibit decreased ß-TrCP degron phosphorylation, ß-TrCP interaction, and ubiquitination, as well as elicit early G(1) arrest. Expression of stabilized ATF4 also had significant consequences in the developing neocortex. Mutant ATF4 expressing cells exhibited positioning and differentiation defects that were attributed to early G(1) arrest, suggesting that neurogenesis is sensitive to ATF4 dosage. We propose that precise regulation of the ATF4 dosage impacts cell cycle control and impinges on neurogenesis.

Cdk5rap2 interacts with pericentrin to maintain the neural progenitor pool in the developing neocortex.

Primary autosomal-recessive microcephaly (MCPH) and Majewski osteodysplastic primordial dwarfism type II (MOPDII) are both genetic diseases that result in decreased brain size at birth. MCPH is thought to arise from alterations in the size of the neural progenitor pool, but the cause of this defect has not been thoroughly explored. We find that one of the genes associated with MCPH, Cdk5rap2, is highly expressed in the neural progenitor pool and that its loss results in a depletion of apical progenitors and increased cell-cycle exit leading to premature neuronal differentiation. We link Cdk5rap2 function to the pericentriolar material protein pericentrin, loss of function of which is associated with MOPDII. Depletion of pericentrin in neural progenitors phenocopies effects of Cdk5rap2 knockdown and results in decreased recruitment of Cdk5rap2 to the centrosome. Our findings uncover a common mechanism, involving aberrations in the neurogenesis program, that may underlie the development of microcephaly in multiple diseases.

Hook3 interacts with PCM1 to regulate pericentriolar material assembly and the timing of neurogenesis.

Centrosome functions are important in multiple brain developmental processes. Proper functioning of the centrosome relies on assembly of protein components into the pericentriolar material. This dynamic assembly is mediated by the trafficking of pericentriolar satellites, which are comprised of centrosomal proteins. Here we demonstrate that trafficking of pericentriolar satellites requires the interaction between Hook3 and Pericentriolar Material 1 (PCM1). Hook3, previously shown to link the centrosome and the nucleus in C. elegans, is recruited to pericentriolar satellites through interaction with PCM1, a protein associated with schizophrenia. Disruption of the Hook3-PCM1 interaction in vivo impairs interkinetic nuclear migration, a featured behavior of embryonic neural progenitors. This in turn leads to overproduction of neurons and premature depletion of the neural progenitor pool in the developing neocortex. These results underscore the importance of centrosomal assembly in neurogenesis and provide potential insights into the etiology of brain developmental diseases related to the centrosome dysfunction.

Alternative functions of core cell cycle regulators in neuronal migration, neuronal maturation, and synaptic plasticity.

Recent studies have demonstrated that boundaries separating a cycling cell from a postmitotic neuron are not as concrete as expected. Novel and unique physiological functions in neurons have been ascribed for proteins fundamentally required for cell cycle progression and control. These "core" cell cycle regulators serve diverse postmitotic functions that span various developmental stages of a neuron, including neuronal migration, axonal elongation, axon pruning, dendrite morphogenesis, and synaptic maturation and plasticity. In this review, we detail the nonproliferative postmitotic roles that these cell cycle proteins have recently been reported to play, the significance of their expression in neurons, mechanistic insight when available, and future prospects.

Disrupted in schizophrenia 1 regulates neuronal progenitor proliferation via modulation of GSK3beta/beta-catenin signaling.

The Disrupted in Schizophrenia 1 (DISC1) gene is disrupted by a balanced chromosomal translocation (1; 11) (q42; q14.3) in a Scottish family with a high incidence of major depression, schizophrenia, and bipolar disorder. Subsequent studies provided indications that DISC1 plays a role in brain development. Here, we demonstrate that suppression of DISC1 expression reduces neural progenitor proliferation, leading to premature cell cycle exit and differentiation. Several lines of evidence suggest that DISC1 mediates this function by regulating GSK3beta. First, DISC1 inhibits GSK3beta activity through direct physical interaction, which reduces beta-catenin phosphorylation and stabilizes beta-catenin. Importantly, expression of stabilized beta-catenin overrides the impairment of progenitor proliferation caused by DISC1 loss of function. Furthermore, GSK3 inhibitors normalize progenitor proliferation and behavioral defects caused by DISC1 loss of function. Together, these results implicate DISC1 in GSK3beta/beta-catenin signaling pathways and provide a framework for understanding how alterations in this pathway may contribute to the etiology of psychiatric disorders.

Deregulation of HDAC1 by p25/Cdk5 in neurotoxicity.

Aberrant cell-cycle activity and DNA damage are emerging as important pathological components in various neurodegenerative conditions. However, their underlying mechanisms are poorly understood. Here, we show that deregulation of histone deacetylase 1 (HDAC1) activity by p25/Cdk5 induces aberrant cell-cycle activity and double-strand DNA breaks leading to neurotoxicity. In a transgenic model for neurodegeneration, p25/Cdk5 activity elicited cell-cycle activity and double-strand DNA breaks that preceded neuronal death. Inhibition of HDAC1 activity by p25/Cdk5 was identified as an underlying mechanism for these events, and HDAC1 gain of function provided potent protection against DNA damage and neurotoxicity in cultured neurons and an in vivo model for ischemia. Our findings outline a pathological signaling pathway illustrating the importance of maintaining HDAC1 activity in the adult neuron. This pathway constitutes a molecular link between aberrant cell-cycle activity and DNA damage and is a potential target for therapeutics against diseases and conditions involving neuronal death.

Control of PERK eIF2alpha kinase activity by the endoplasmic reticulum stress-induced molecular chaperone P58IPK.

P58(IPK) is an Hsp40 family member known to inhibit the interferon (IFN)-induced, double-stranded RNA-activated, eukaryotic initiation factor 2alpha (eIF2alpha) protein kinase R (PKR) by binding to its kinase domain. We find that the stress of unfolded proteins in the endoplasmic reticulum (ER) activates P58(IPK) gene transcription through an ER stress-response element in its promoter region. P58(IPK) interacts with and inhibits the PKR-like ER-localized eIF2alpha kinase PERK, which is normally activated during the ER-stress response to protect cells from ER stress by attenuating protein synthesis and reducing ER client protein load. Levels of phosphorylated eIF2alpha were lower in ER-stressed P58(IPK)-overexpressing cells and were enhanced in P58(IPK) mutant cells. In the ER-stress response, PKR-like ER kinase (PERK)-mediated translational repression is transient and is followed by translational recovery and enhanced expression of genes that increase the capacity of the ER to process client proteins. The absence of P58(IPK) resulted in increased expression levels of two ER stress-inducible genes, BiP and Chop, consistent with the enhanced eIF2alpha phosphorylation in the P58(IPK) deletion cells. Our studies suggest that P58(IPK) induction during the ER-stress response represses PERK activity and plays a functional role in the expression of downstream markers of PERK activity in the later phase of the ER-stress response.