LLM3D: a log-linear modeling-based method to predict functional gene regulatory interactions from genome-wide expression data

Aryl hydrocarbon receptor restricts axon regeneration of DRG neurons in response to injury

Injured neurons sense environmental cues to balance neural protection and axon regeneration, but the mechanisms are unclear. Here, we unveil aryl hydrocarbon receptor (AhR), a ligand-activated bHLH-PAS transcription factor, as a molecular sensor and key regulator of acute stress response at the expense of axon regeneration. We demonstrate responsiveness of DRG sensory neurons to AhR signaling, which functions to inhibit axon regeneration. Conditional Ahr deletion in neurons accelerates axon regeneration after sciatic nerve injury. Ahr deletion partially mimics the conditioning lesion in priming DRG to initiate axonogenesis gene programs; upon peripheral axotomy, Ahr ablation suppresses inflammation and stress signaling while augmenting pro-growth pathways. Moreover, comparative transcriptomics revealed signaling interactions between AhR and HIF-1α, two structurally related bHLH-PAS α units that share the dimerization partner Arnt/HIF-1β. Functional assays showed that the growth advantage of AhR-deficient DRG neurons requires HIF-1α; but in the absence of Arnt, DRG neurons can still mount a regenerative response. We further unveil a link between bHLH-PAS transcription factors and DNA hydroxymethylation in response to peripheral axotomy, while RNA-seq of DRG neurons and neuronal single cell RNA-seq analysis revealed a link of AhR regulon to RNA regulation and integrated stress response (ISR). Altogether, AhR activation favors stress coping and inflammation at the expense of axon regeneration; targeting AhR has the potential to enhance nerve repair.

The Jun-dependent axon regeneration gene program: Jun promotes regeneration over plasticity

The regeneration-associated gene (RAG) expression program is activated in injured peripheral neurons after axotomy and enables long-distance axon re-growth. Over 1000 genes are regulated, and many transcription factors are upregulated or activated as part of this response. However, a detailed picture of how RAG expression is regulated is lacking. In particular, the transcriptional targets and specific functions of the various transcription factors are unclear. Jun was the first-regeneration-associated transcription factor identified and the first shown to be functionally important. Here we fully define the role of Jun in the RAG expression program in regenerating facial motor neurons. At 1, 4 and 14 days after axotomy, Jun upregulates 11, 23 and 44% of the RAG program, respectively. Jun functions relevant to regeneration include cytoskeleton production, metabolic functions and cell activation, and the downregulation of neurotransmission machinery. In silico analysis of promoter regions of Jun targets identifies stronger over-representation of AP1-like sites than CRE-like sites, although CRE sites were also over-represented in regions flanking AP1 sites. Strikingly, in motor neurons lacking Jun, an alternative SRF-dependent gene expression program is initiated after axotomy. The promoters of these newly expressed genes exhibit over-representation of CRE sites in regions near to SRF target sites. This alternative gene expression program includes plasticity-associated transcription factors and leads to an aberrant early increase in synapse density on motor neurons. Jun thus has the important function in the early phase after axotomy of pushing the injured neuron away from a plasticity response and towards a regenerative phenotype.

Age-related loss of axonal regeneration is reflected by the level of local translation

Regeneration capacity is reduced as CNS axons mature. Using laser-mediated axotomy, proteomics and puromycin-based tagging of newly-synthesized proteins in a human embryonic stem cell-derived neuron culture system that allows isolation of axons from cell bodies, we show here that efficient regeneration in younger axons (d45 in culture) is associated with local axonal protein synthesis (local translation). Enhanced regeneration, promoted by co-culture with human glial precursor cells, is associated with increased axonal synthesis of proteins, including those constituting the translation machinery itself. Reduced regeneration, as occurs with the maturation of these axons by d65 in culture, correlates with reduced levels of axonal proteins involved in translation and an inability to respond by increased translation of regeneration promoting axonal mRNAs released from stress granules. Together, our results provide evidence that, as in development and in the PNS, local translation contributes to CNS axon regeneration.

Myelin basic protein enhances axonal regeneration from neural progenitor cells

Introduction

Stem cell therapy using neural progenitor cells (NPCs) shows promise in mitigating the debilitating effects of spinal cord injury (SCI). Notably, myelin stimulates axonal regeneration from mammalian NPCs. This led us to hypothesize that myelin-associated proteins may contribute to axonal regeneration from NPCs.

Methods

We conducted an R-based bioinformatics analysis to identify key gene(s) that may participate in myelin-associated axonal regeneration from murine NPCs, which identified the serine protease myelin basic protein (Mbp). We employed E12 murine NPCs, E14 rat NPCs, and human iPSC-derived Day 1 NPCs (D1 hNPCs) with or without CRISPR/Cas9-mediated Mbp knockout in combination with rescue L1-70 overexpression, constitutively-active VP16-PPARγ2, or the PPARγ agonist ciglitazone. A murine dorsal column crush model of SCI utilizing porous collagen-based scaffolding (PCS)-seeded murine NPCs with or without stable Mbp overexpression was used to assess locomotive recovery and axonal regeneration in vivo.

Results

Myelin promotes axonal outgrowth from NPCs in an Mbp-dependent manner and that Mbp’s stimulatory effects on NPC neurite outgrowth are mediated by Mbp’s production of L1-70. Furthermore, we determined that Mbp/L1-70’s stimulatory effects on NPC neurite outgrowth are mediated by PPARγ-based repression of neuron differentiation-associated gene expression and PPARγ-based Erk1/2 activation. In vivo, PCS-seeded murine NPCs stably overexpressing Mbp significantly enhanced locomotive recovery and axonal regeneration in post-SCI mice.

Conclusions

We discovered that Mbp supports axonal regeneration from mammalian NPCs through the novel Mbp/L1cam/Pparγ signaling pathway. This study suggests that bioengineered, NPC-based interventions can promote axonal regeneration and functional recovery post-SCI.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13578-021-00584-7.

Recently Evolved Enhancers Emerge with High Interindividual Variability and Less Frequently Associate with Disease

Mutations in non-coding regulatory DNA such as enhancers underlie a wide variety of diseases including developmental disorders and cancer. As enhancers rapidly evolve, understanding their function and configuration in non-human disease models can have important clinical applications. Here, we analyze enhancer configurations in tissues isolated from the common marmoset, a widely used primate model for human disease. Integrating these data with human and mouse data, we find that enhancers containing trait-associated variants are preferentially conserved. In contrast, most human-specific enhancers are highly variable between individuals, with a subset failing to contact promoters. These are located further away from genes and more often reside in inactive B-compartments. Our data show that enhancers typically emerge as instable elements with minimal biological impact prior to their integration in a transcriptional program. Furthermore, our data provide insight into which trait variations in enhancers can be faithfully modeled using the common marmoset.

Hominin-specific regulatory elements selectively emerged in oligodendrocytes and are disrupted in autism patients

Speciation is associated with substantial rewiring of the regulatory circuitry underlying the expression of genes. Determining which changes are relevant and underlie the emergence of the human brain or its unique susceptibility to neural disease has been challenging. Here we annotate changes to gene regulatory elements (GREs) at cell type resolution in the brains of multiple primate species spanning most of primate evolution. We identify a unique set of regulatory elements that emerged in hominins prior to the separation of humans and chimpanzees. We demonstrate that these hominin gains perferentially affect oligodendrocyte function postnatally and are preferentially affected in the brains of autism patients. This preference is also observed for human-specific GREs suggesting this system is under continued selective pressure. Our data provide a roadmap of regulatory rewiring across primate evolution providing insight into the genomic changes that underlie the emergence of the brain and its susceptibility to neural disease.

The understanding of the changes regulating gene expression relevant for the emergence of the human brain and its susceptibility to disease is limited. Here, the authors identified a set of regulatory elements that evolved in hominins affecting oligodendrocyte function, and link these to autism.

Epigenomic signatures underpin the axonal regenerative ability of dorsal root ganglia sensory neurons

Axonal injury results in regenerative success or failure, depending on whether the axon lies in the peripheral or the CNS, respectively. The present study addresses whether epigenetic signatures in dorsal root ganglia discriminate between regenerative and non-regenerative axonal injury. Chromatin immunoprecipitation for the histone 3 (H3) post-translational modifications H3K9ac, H3K27ac and H3K27me3; an assay for transposase-accessible chromatin; and RNA sequencing were performed in dorsal root ganglia after sciatic nerve or dorsal column axotomy. Distinct histone acetylation and chromatin accessibility signatures correlated with gene expression after peripheral, but not central, axonal injury. DNA-footprinting analyses revealed new transcriptional regulators associated with regenerative ability. Machine-learning algorithms inferred the direction of most of the gene expression changes. Neuronal conditional deletion of the chromatin remodeler CCCTC-binding factor impaired nerve regeneration, implicating chromatin organization in the regenerative competence. Altogether, the present study offers the first epigenomic map providing insight into the transcriptional response to injury and the differential regenerative ability of sensory neurons.

The Dorsal Column Lesion Model of Spinal Cord Injury and Its Use in Deciphering the Neuron-Intrinsic Injury Response

The neuron‐intrinsic response to axonal injury differs markedly between neurons of the peripheral and central nervous system. Following a peripheral lesion, a robust axonal growth program is initiated, whereas neurons of the central nervous system do not mount an effective regenerative response. Increasing the neuron‐intrinsic regenerative response would therefore be one way to promote axonal regeneration in the injured central nervous system. The large‐diameter sensory neurons located in the dorsal root ganglia are pseudo‐unipolar neurons that project one axon branch into the spinal cord, and, via the dorsal column to the brain stem, and a peripheral process to the muscles and skin. Dorsal root ganglion neurons are ideally suited to study the neuron‐intrinsic injury response because they exhibit a successful growth response following peripheral axotomy, while they fail to do so after a lesion of the central branch in the dorsal column. The dorsal column injury model allows the neuron‐intrinsic regeneration response to be studied in the context of a spinal cord injury. Here we will discuss the advantages and disadvantages of this model. We describe the surgical methods used to implement a lesion of the ascending fibers, the anatomy of the sensory afferent pathways and anatomical, electrophysiological, and behavioral techniques to quantify regeneration and functional recovery. Subsequently we review the results of experimental interventions in the dorsal column lesion model, with an emphasis on the molecular mechanisms that govern the neuron‐intrinsic injury response and manipulations of these after central axotomy. Finally, we highlight a number of recent advances that will have an impact on the design of future studies in this spinal cord injury model, including the continued development of adeno‐associated viral vectors likely to improve the genetic manipulation of dorsal root ganglion neurons and the use of tissue clearing techniques enabling 3D reconstruction of regenerating axon tracts. © 2018 The Authors. Developmental Neurobiology Published by Wiley Periodicals, Inc. Develop Neurobiol 00: 000–000, 2018

Combined chondroitinase and KLF7 expression reduce net retraction of sensory and CST axons from sites of spinal injury

Axon regeneration in the central nervous system is limited both by inhibitory extracellular cues and by an intrinsically low capacity for axon growth in some CNS populations. Chondroitin sulfate proteoglycans (CSPGs) are well-studied inhibitors of axon growth in the CNS, and degradation of CSPGs by chondroitinase has been shown to improve the extension of injured axons. Alternatively, axon growth can be improved by targeting the neuron-intrinsic growth capacity through forced expression of regeneration-associated transcription factors. For example, a transcriptionally active chimera of Krüppel-like Factor 7 (KLF7) and a VP16 domain improves axon growth when expressed in corticospinal tract neurons. Here we tested the hypothesis that combined expression of chondroitinase and VP16-KLF7 would lead to further improvements in axon growth after spinal injury. Chondroitinase was expressed by viral transduction of cells in the spinal cord, while VP16-KLF7 was virally expressed in sensory neurons of the dorsal root ganglia or corticospinal tract (CST) neurons. After transection of the dorsal columns, both chondroitinase and VP16-KLF7 increased the proximity of severed sensory axons to the injury site. Similarly, after complete crush injuries, VP16-KLF7 expression increased the approach of CST axons to the injury site. In neither paradigm however, did single or combined treatment with chondroitinase or VP16-KLF7 enable regenerative growth distal to the injury. These results substantiate a role for CSPG inhibition and low KLF7 activity in determining the net retraction of axons from sites of spinal injury, while suggesting that additional factors act to limit a full regenerative response.

Epigenetic profiling reveals a developmental decrease in promoter accessibility during cortical maturation in vivo

Axon regeneration in adult central nervous system (CNS) is limited in part by a developmental decline in the ability of injured neurons to re-express needed regeneration associated genes (RAGs). Adult CNS neurons may lack appropriate pro-regenerative transcription factors, or may display chromatin structure that restricts transcriptional access to RAGs. Here we performed epigenetic profiling around the promoter regions of key RAGs, and found progressive restriction across a time course of cortical maturation. These data identify a potential intrinsic constraint to axon growth in adult CNS neurons. Neurite outgrowth from cultured postnatal cortical neurons, however, proved insensitive to treatments that improve axon growth in other cell types, including combinatorial overexpression of AP1 factors, overexpression of histone acetyltransferases, and pharmacological inhibitors of histone deacetylases. This insensitivity could be due to intermediate chromatin closure at the time of culture, and highlights important differences in cell culture models used to test potential pro-regenerative interventions.

What makes a RAG regeneration associated?

Regenerative failure remains a significant barrier for functional recovery after central nervous system (CNS) injury. As such, understanding the physiological processes that regulate axon regeneration is a central focus of regenerative medicine. Studying the gene transcription responses to axon injury of regeneration competent neurons, such as those of the peripheral nervous system (PNS), has provided insight into the genes associated with regeneration. Though several individual “regeneration-associated genes” (RAGs) have been identified from these studies, the response to injury likely regulates the expression of functionally coordinated and complementary gene groups. For instance, successful regeneration would require the induction of genes that drive the intrinsic growth capacity of neurons, while simultaneously downregulating the genes that convey environmental inhibitory cues. Thus, this view emphasizes the transcriptional regulation of gene “programs” that contribute to the overall goal of axonal regeneration. Here, we review the known RAGs, focusing on how their transcriptional regulation can reveal the underlying gene programs that drive a regenerative phenotype. Finally, we will discuss paradigms under which we can determine whether these genes are injury-associated, or indeed necessary for regeneration.

Understanding the neural repair-promoting properties of olfactory ensheathing cells

Olfactory ensheathing glial cells (OECs) are a specialized type of glia that form a continuously aligned cellular pathway that actively supports unprecedented regeneration of primary olfactory axons from the periphery into the central nervous system. Implantation of OECs stimulates neural repair in experimental models of spinal cord, brain and peripheral nerve injury and delays disease progression in animal models for neurodegenerative diseases like amyotrophic lateral sclerosis. OECs implanted in the injured spinal cord display a plethora of pro-regenerative effects; they promote axonal regeneration, reorganize the glial scar, remyelinate axons, stimulate blood vessel formation, have phagocytic properties and modulate the immune response. Recently genome wide transcriptional profiling and proteomics analysis combined with classical or larger scale "medium-throughput" bioassays have provided novel insights into the molecular mechanism that endow OECs with their pro-regenerative properties. Here we review these studies and show that the gaps that existed in our understanding of the molecular basis of the reparative properties of OECs are narrowing. OECs express functionally connected sets of genes that can be linked to at least 10 distinct processes directly relevant to neural repair. The data indicate that OECs exhibit a range of synergistic cellular activities, including active and passive stimulation of axon regeneration (by secretion of growth factors, axon guidance molecules and basement membrane components) and critical aspects of tissue repair (by structural remodeling and support, modulation of the immune system, enhancement of neurotrophic and antigenic stimuli and by metabolizing toxic macromolecules). Future experimentation will have to further explore the newly acquired knowledge to enhance the therapeutic potential of OECs.

Comparison of targeted maximum likelihood and shrinkage estimators of parameters in gene networks

Gene regulatory networks, in which edges between nodes describe interactions between transcription factors (TFs) and their target genes, model regulatory interactions that determine the cell-type and condition-specific expression of genes. Regression methods can be used to identify TF-target gene interactions from gene expression and DNA sequence data. The response variable, i.e. observed gene expression, is modeled as a function of many predictor variables simultaneously. In practice, it is generally not possible to select a single model that clearly achieves the best fit to the observed experimental data and the selected models typically contain overlapping sets of predictor variables. Moreover, parameters that represent the marginal effect of the individual predictors are not always present. In this paper, we use the statistical framework of estimation of variable importance to define variable importance as a parameter of interest and study two different estimators of this parameter in the context of gene regulatory networks. On yeast data we show that the resulting parameter has a biologically appealing interpretation. We apply the proposed methodology on mammalian gene expression data to gain insight into the temporal activity of TFs that underly gene expression changes in F11 cells in response to Forskolin stimulation.

Achieving CNS axon regeneration by manipulating convergent neuro-immune signaling

After central nervous system (CNS) trauma, axons have a low capacity for regeneration. Regeneration failure is associated with a muted regenerative response of the neuron itself, combined with a growth-inhibitory and cytotoxic post-injury environment. After spinal cord injury (SCI), resident and infiltrating immune cells (especially microglia/macrophages) contribute significantly to the growth-refractory milieu near the lesion. By targeting both the regenerative potential of the axon and the cytotoxic phenotype of microglia/macrophages, we may be able to improve CNS repair after SCI. In this review, we discuss molecules shown to impact CNS repair by affecting both immune cells and neurons. Specifically, we provide examples of pattern recognition receptors, integrins, cytokines/chemokines, nuclear receptors and galectins that could improve CNS repair. In many cases, signaling by these molecules is complex and may have contradictory effects on recovery depending on the cell types involved or the model studied. Despite this caveat, deciphering convergent signaling pathways on immune cells (which affect axon growth indirectly) and neurons (direct effects on axon growth) could improve repair and recovery after SCI. Future studies must continue to consider how regenerative therapies targeting neurons impact other cells in the pathological CNS. By identifying molecules that simultaneously improve axon regenerative capacity and drive the protective, growth-promoting phenotype of immune cells, we may discover SCI therapies that act synergistically to improve CNS repair and functional recovery.

Molecular target discovery for neural repair in the functional genomics era

A comprehensive understanding of the molecular pathways activated by traumatic neural injury is of major importance for the development of treatments for spinal cord injury (SCI). High-throughput gene expression profiling is a powerful approach to reveal genome-wide changes in gene expression during a specific biological process. Microarray analysis of injured nerves or neurons would ideally generate new hypotheses concerning the progression or deregulation of injury- and repair-related biological processes, such as neural scar formation and axon regeneration. These hypotheses should subsequently be tested experimentally and would eventually provide the molecular substrates for the development of novel therapeutics. Over the last decade, this approach has elucidated numerous extrinsic (mostly neural scar-associated) as well as neuron-intrinsic genes that are regulated following an injury. To date, the main challenge is to translate the observed injury-induced gene expression changes into a mechanistic framework to understand their functional implications. To achieve this, research on neural repair will have to adopt the conceptual advances and analytical tools provided by the functional genomics and systems biology revolution. Based on progress made in bioinformatics, high-throughput and high-content functional cellular screening, and in vivo gene transfer technology, we propose a multistep "roadmap" that provides an integrated strategy for molecular target discovery for repair of the injured spinal cord.

A gene network perspective on axonal regeneration

The regenerative capacity of injured neurons in the central nervous system is limited due to the absence of a robust neuron-intrinsic injury-induced gene response that supports axon regeneration. In peripheral neurons axotomy induces a large cohort of regeneration-associated genes (RAGs). The forced expression of some of these RAGs in injured neurons has some beneficial effect on axon regeneration, but the reported effects are rather small. Transcription factors (TFs) provide a promising class of RAGs. TFs are hubs in the regeneration-associated gene network, and potentially control the coordinate expression of many RAGs simultaneously. Here we discuss the use of combined experimental and computational methods to identify novel regeneration-associated TFs with a key role in initiating and maintaining the RAG-response in injured neurons. We propose that a relatively small number of hub TFs with multiple functional connections in the RAG network might provide attractive new targets for gene-based and/or pharmacological approaches to promote axon regeneration in the central nervous system.

Identification of context-specific gene regulatory networks with GEMULA--gene expression modeling using LAsso

MOTIVATION

Gene regulatory networks, in which edges between nodes describe interactions between transcriptional regulators and their target genes, determine the coordinated spatiotemporal expression of genes. Especially in higher organisms, context-specific combinatorial regulation by transcription factors (TFs) is believed to determine cellular states and fates. TF-target gene interactions can be studied using high-throughput techniques such as ChIP-chip or ChIP-Seq. These experiments are time and cost intensive, and further limited by, for instance, availability of high affinity TF antibodies. Hence, there is a practical need for methods that can predict TF-TF and TF-target gene interactions in silico, i.e. from gene expression and DNA sequence data alone. We propose GEMULA, a novel approach based on linear models to predict TF-gene expression associations and TF-TF interactions from experimental data. GEMULA is based on linear models, fast and considers a wide range of biologically plausible models that describe gene expression data as a function of predicted TF binding to gene promoters.

RESULTS

We show that models inferred with GEMULA are able to explain roughly 70% of the observed variation in gene expression in the yeast heat shock response. The functional relevance of the inferred TF-TF interactions in these models are validated by different sources of independent experimental evidence. We also have applied GEMULA to an in vitro model of neuronal outgrowth. Our findings confirm existing knowledge on gene regulatory interactions underlying neuronal outgrowth, but importantly also generate new insights into the temporal dynamics of this gene regulatory network that can now be addressed experimentally.

AVAILABILITY

The GEMULA R-package is available from http://www.few.vu.nl/~degunst/gemula_1.0.tar.gz.