Reconstructed cell fate-regulatory programs in stem cells reveal hierarchies and key factors of neurogenesis.

Cell lineages, which shape the body architecture and specify cell functions, derive from the integration of a plethora of cell intrinsic and extrinsic signals. These signals trigger a multiplicity of decisions at several levels to modulate the activity of dynamic gene regulatory networks (GRNs), which ensure both general and cell-specific functions within a given lineage, thereby establishing cell fates. Significant knowledge about these events and the involved key drivers comes from homogeneous cell differentiation models. Even a single chemical trigger, such as the morphogen all-trans retinoic acid (RA), can induce the complex network of gene-regulatory decisions that matures a stem/precursor cell to a particular step within a given lineage. Here we have dissected the GRNs involved in the RA-induced neuronal or endodermal cell fate specification by integrating dynamic RXRA binding, chromatin accessibility, epigenetic promoter epigenetic status, and the transcriptional activity inferred from RNA polymerase II mapping and transcription profiling. Our data reveal how RA induces a network of transcription factors (TFs), which direct the temporal organization of cognate GRNs, thereby driving neuronal/endodermal cell fate specification. Modeling signal transduction propagation using the reconstructed GRNs indicated critical TFs for neuronal cell fate specification, which were confirmed by CRISPR/Cas9-mediated genome editing. Overall, this study demonstrates that a systems view of cell fate specification combined with computational signal transduction models provides the necessary insight in cellular plasticity for cell fate engineering. The present integrated approach can be used to monitor the in vitro capacity of (engineered) cells/tissues to establish cell lineages for regenerative medicine.

The inactive X chromosome is epigenetically unstable and transcriptionally labile in breast cancer.

Disappearance of the Barr body is considered a hallmark of cancer, although whether this corresponds to genetic loss or to epigenetic instability and transcriptional reactivation is unclear. Here we show that breast tumors and cell lines frequently display major epigenetic instability of the inactive X chromosome, with highly abnormal 3D nuclear organization and global perturbations of heterochromatin, including gain of euchromatic marks and aberrant distributions of repressive marks such as H3K27me3 and promoter DNA methylation. Genome-wide profiling of chromatin and transcription reveal modified epigenomic landscapes in cancer cells and a significant degree of aberrant gene activity from the inactive X chromosome, including several genes involved in cancer promotion. We demonstrate that many of these genes are aberrantly reactivated in primary breast tumors, and we further demonstrate that epigenetic instability of the inactive X can lead to perturbed dosage of X-linked factors. Taken together, our study provides the first integrated analysis of the inactive X chromosome in the context of breast cancer and establishes that epigenetic erosion of the inactive X can lead to the disappearance of the Barr body in breast cancer cells. This work offers new insights and opens up the possibility of exploiting the inactive X chromosome as an epigenetic biomarker at the molecular and cytological levels in cancer.

Development of biotin-retinoid conjugates as chemical probes for analysis of retinoid function.

Herein, we report the rational design, synthesis and biological evaluation of conjugates consisting of the synthetic retinoid Am580 and biotin connected via a linker moiety. We found that the linking substructure between the retinoid part and the biotin part is critical for retaining the biological activity. Conjugate 4 with a shorter linker showed similar potency to endogenous retinoid ATRA (1) and the parent compound Am580 (2) for neural differentiation of mouse embryotic carcinoma P19 cells, and showed the same pattern of induction of gene expression. It is expected to be useful as a probe for investigations of retinoid function. The design rationale and structure-activity relationship of the linker moiety are expected to be helpful for developing biotin conjugates of other nuclear receptor ligands.

Epimetheus - a multi-profile normalizer for epigenomic sequencing data.

BACKGROUND: Exponentially increasing numbers of NGS-based epigenomic datasets in public repositories like GEO constitute an enormous source of information that is invaluable for integrative and comparative studies of gene regulatory mechanisms. One of today's challenges for such studies is to identify functionally informative local and global patterns of chromatin states in order to describe the regulatory impact of the epigenome in normal cell physiology and in case of pathological aberrations. Critically, the most preferred Chromatin ImmunoPrecipitation-Sequencing (ChIP-Seq) is inherently prone to significant variability between assays, which poses significant challenge on comparative studies. One challenge concerns data normalization to adjust sequencing depth variation. RESULTS: Currently existing tools either apply linear scaling corrections and/or are restricted to specific genomic regions, which can be prone to biases. To overcome these restrictions without any external biases, we developed Epimetheus, a genome-wide quantile-based multi-profile normalization tool for histone modification data and related datasets. CONCLUSIONS: Epimetheus has been successfully used to normalize epigenomics data in previous studies on X inactivation in breast cancer and in integrative studies of neuronal cell fate acquisition and tumorigenic transformation; Epimetheus is freely available to the scientific community.

LOGIQA: a database dedicated to long-range genome interactions quality assessment.

BACKGROUND: Proximity ligation-mediated methods are essential to study the impact of three-dimensional chromatin organization on gene programming. Albeit significant progress has been made in the development of computational tools that assess long-range chromatin interactions, next to nothing is known about the quality of the generated datasets. METHOD: We have developed LOGIQA ( www.ngs-qc.org/logiqa ), a database hosting quality scores for long-range genome interaction assays, accessible through a user-friendly web-based environment. RESULTS: Currently, LOGIQA harbors QC scores for >900 datasets, which provides a global view of their relative quality and reveals the impact of genome size, coverage and other technical aspects. LOGIQA provides a user-friendly dataset query panel and a genome viewer to assess local genome-interaction maps at different resolution and quality-assessment conditions. CONCLUSIONS: LOGIQA is the first database hosting quality scores dedicated to long-range chromatin interaction assays, which in addition provides a platform for visualizing genome interactions made available by the scientific community.

Genome-wide studies of nuclear receptors in cell fate decisions.

Nuclear receptors (NRs) are important mediators of the information encoded in the chemical structure of its corresponding ligand, as they interpret such information in the context of the cell identity and physiological status and convert it into sequential transcription regulatory events. At the cell level this can result in temporally coordinated processes such as cell fate transitions, which comprise the regulation of a plethora of gene programs including among others regulation of cell proliferation, metabolism and specific functionalities that are acquired by the differentiated cell. While both the early steps of nuclear receptor function and their impact on animal/organ physiology is rather well understood, little is known about the dynamic gene networks that ultimately cause a particular (cell) physiological phenomenon induced by the cognate NR ligand/hormone. Thanks to advances in massive parallel sequencing and bioinformatics analyses of genome-wide data sets, time has come for the development of NR systems biology. Indeed it is now possible to integrate global transcription factor binding, epigenetic chromatin histone and DNA modification patterns with transcriptomes and 3-dimensional chromatin structures, extract decision points in temporal studies and decipher the temporal control of gene networks that are the ultimate genetic readouts of NR ligand-induced physiological phenomena. In this review we will summarize the chronology of the development of increasingly larger data sets for NR action, with a particular focus on studies performed with the RAR/RXR nuclear receptor family, and discuss the present attempts to integrate a multitude of genome-wide data sets in the ultimate context of the temporal 3-dimensional chromatin structure.

A quality control system for profiles obtained by ChIP sequencing.

The absence of a quality control (QC) system is a major weakness for the comparative analysis of genome-wide profiles generated by next-generation sequencing (NGS). This concerns particularly genome binding/occupancy profiling assays like chromatin immunoprecipitation (ChIP-seq) but also related enrichment-based studies like methylated DNA immunoprecipitation/methylated DNA binding domain sequencing, global run on sequencing or RNA-seq. Importantly, QC assessment may significantly improve multidimensional comparisons that have great promise for extracting information from combinatorial analyses of the global profiles established for chromatin modifications, the bindings of epigenetic and chromatin-modifying enzymes/machineries, RNA polymerases and transcription factors and total, nascent or ribosome-bound RNAs. Here we present an approach that associates global and local QC indicators to ChIP-seq data sets as well as to a variety of enrichment-based studies by NGS. This QC system was used to certify >5600 publicly available data sets, hosted in a database for data mining and comparative QC analyses.

Integrative genomics to dissect retinoid functions.

Retinoids and rexinoids, as all other ligands of the nuclear receptor (NR) family, act as ligand-regulated trans-acting transcription factors that bind to cis-acting DNA regulatory elements in the promoter regions of target genes (for reviews see [12, 22, 23, 26, 36]). Ligand binding modulates the communication functions of the receptor with the intracellular environment, which essentially entails receptor-protein and receptor-DNA or receptor-chromatin interactions. In this communication network, the receptor simultaneously serves as both intracellular sensor and regulator of cell/organ functions. Receptors are "intelligent" mediators of the information encoded in the chemical structure of a nuclear receptor ligand, as they interpret this information in the context of cellular identity and cell-physiological status and convert it into a dynamic chain of receptor-protein and receptor-DNA interactions. To process input and output information, they are composed of a modular structure with several domains that have evolved to exert particular molecular recognition functions. As detailed in other chapters in this volume, the main functional domains are the DNA-binding (DBD) and ligand-binding (LBD) [5-7, 38, 56, 71]. The LBD serves as a dual input-output information processor. Inputs, such as ligand binding or receptor phosphorylations, induce allosteric changes in receptor surfaces that serve as docking sites for outputs, such as subunits of transcription and epigenetic machineries or enzyme complexes. The complexity of input and output signals and their interdependencies is far from being understood.

Single-tube linear DNA amplification (LinDA) for robust ChIP-seq.

Genome-wide profiling of transcription factors based on massive parallel sequencing of immunoprecipitated chromatin (ChIP-seq) requires nanogram amounts of DNA. Here we describe a high-fidelity, single-tube linear DNA amplification method (LinDA) for ChIP-seq and reChIP-seq with picogram DNA amounts obtained from a few thousand cells. This amplification technology will facilitate global analyses of transcription-factor binding and chromatin with very small cell populations, such as stem or cancer-initiating cells.

C286, an orally available retinoic acid receptor ß agonist drug, regulates multiple pathways to achieve spinal cord injury repair.

Retinoic acid receptor ß2 (RARß2) is an emerging therapeutic target for spinal cord injuries (SCIs) with a unique multimodal regenerative effect. We have developed a first-in-class RARß agonist drug, C286, that modulates neuron-glial pathways to induce functional recovery in a rodent model of sensory root avulsion. Here, using genome-wide and pathway enrichment analysis of avulsed rats' spinal cords, we show that C286 also influences the extracellular milieu (ECM). Protein expression studies showed that C286 upregulates tenascin-C, integrin-α9, and osteopontin in the injured cord. Similarly, C286 remodulates these ECM molecules, hampers inflammation and prevents tissue loss in a rodent model of spinal cord contusion C286. We further demonstrate C286's efficacy in human iPSC-derived neurons, with treatment resulting in a significant increase in neurite outgrowth. Additionally, we identify a putative efficacy biomarker, S100B, which plasma levels correlated with axonal regeneration in nerve-injured rats. We also found that other clinically available retinoids, that are not RARß specific agonists, did not lead to functional recovery in avulsed rats, demonstrating the requirement for RARß specific pathways in regeneration. In a Phase 1 trial, the single ascending dose (SAD) cohorts showed increases in expression of RARß2 in white blood cells correlative to increased doses and at the highest dose administered, the pharmacokinetics were similar to the rat proof of concept (POC) studies. Collectively, our data suggests that C286 signalling in neurite/axonal outgrowth is conserved between species and across nerve injuries. This warrants further clinical testing of C286 to ascertain POC in a broad spectrum of neurodegenerative conditions.

Characterising ChIP-seq binding patterns by model-based peak shape deconvolution.

BACKGROUND: Chromatin immunoprecipitation combined with massive parallel sequencing (ChIP-seq) is widely used to study protein-chromatin interactions or chromatin modifications at genome-wide level. Sequence reads that accumulate locally at the genome (peaks) reveal loci of selectively modified chromatin or specific sites of chromatin-binding factors. Computational approaches (peak callers) have been developed to identify the global pattern of these sites, most of which assess the deviation from background by applying distribution statistics. RESULTS: We have implemented MeDiChISeq, a regression-based approach, which--by following a learning process--defines a representative binding pattern from the investigated ChIP-seq dataset. Using this model MeDiChISeq identifies significant genome-wide patterns of chromatin-bound factors or chromatin modification. MeDiChISeq has been validated for various publicly available ChIP-seq datasets and extensively compared with other peak callers. CONCLUSIONS: MeDiChI-Seq has a high resolution when identifying binding events, a high degree of peak-assessment reproducibility in biological replicates, a low level of false calls and a high true discovery rate when evaluated in the context of gold-standard benchmark datasets. Importantly, this approach can be applied not only to 'sharp' binding patterns--like those retrieved for transcription factors (TFs)--but also to the broad binding patterns seen for several histone modifications. Notably, we show that at high sequencing depths, MeDiChISeq outperforms other algorithms due to its powerful peak shape recognition capacity which facilitates discerning significant binding events from spurious background enrichment patterns that are enhanced with increased sequencing depths.

Retinoids: Mechanisms of Action in Neuronal Cell Fate Acquisition.

Neuronal differentiation has been shown to be directed by retinoid action during embryo development and has been exploited in various in vitro cell differentiation systems. In this review, we summarize the role of retinoids through the activation of their specific retinoic acid nuclear receptors during embryo development and also in a variety of in vitro strategies for neuronal differentiation, including recent efforts in driving cell specialization towards a range of neuronal subtypes and glial cells. Finally, we highlight the role of retinoic acid in recent protocols recapitulating nervous tissue complexity (cerebral organoids). Overall, we expect that this effort might pave the way for exploring the usage of specific synthetic retinoids for directing complex nervous tissue differentiation.

Three-dimensional molecular cartography of human cerebral organoids revealed by double-barcoded spatial transcriptomics.

Spatially resolved transcriptomics is revolutionizing our understanding of complex tissues, but scaling these approaches to multiple tissue sections and three-dimensional tissue reconstruction remains challenging and cost prohibitive. In this work, we present a low-cost strategy for manufacturing molecularly double-barcoded DNA arrays, enabling large-scale spatially resolved transcriptomics studies. We applied this technique to spatially resolve gene expression in several human brain organoids, including the reconstruction of a three-dimensional view from multiple consecutive sections, revealing gene expression heterogeneity throughout the tissue.

Synergistic activation of RARß and RARγ nuclear receptors restores cell specialization during stem cell differentiation by hijacking RARα-controlled programs.

How cells respond to different external cues to develop along defined cell lineages to form complex tissues is a major question in systems biology. Here, we investigated the potential of retinoic acid receptor (RAR)-selective synthetic agonists to activate the gene regulatory programs driving cell specialization during nervous tissue formation from embryonic carcinoma (P19) and mouse embryonic (E14) stem cells. Specifically, we found that the synergistic activation of the RARß and RARγ by selective ligands (BMS641 or BMS961) induces cell maturation to specialized neuronal subtypes, and to astrocytes and oligodendrocyte precursors. Using RAR isotype knockout lines exposed to RAR-specific agonists, interrogated by global transcriptome landscaping and in silico modeling of transcription regulatory signal propagation, revealed major RARα-driven gene programs essential for optimal neuronal cell specialization and hijacked by the synergistic activation of the RARß and RARγ receptors. Overall, this study provides a systems biology view of the gene programs accounting for the previously observed redundancy between RARs, paving the way toward their potential use for directing cell specialization during nervous tissue formation.

Real-time SARS-CoV-2 diagnostic and variants tracking over multiple candidates using nanopore DNA sequencing.

Since December 2019, a novel coronavirus responsible for a severe acute respiratory syndrome (SARS-CoV-2) is accountable for a major pandemic situation. The emergence of the B.1.1.7 strain, as a highly transmissible variant has accelerated the world-wide interest in tracking SARS-CoV-2 variants' occurrence. Similarly, other extremely infectious variants, were described and further others are expected to be discovered due to the long period of time on which the pandemic situation is lasting. All described SARS-CoV-2 variants present several mutations within the gene encoding the Spike protein, involved in host receptor recognition and entry into the cell. Hence, instead of sequencing the whole viral genome for variants' tracking, herein we propose to focus on the SPIKE region to increase the number of candidate samples to screen at once; an essential aspect to accelerate diagnostics, but also variants' emergence/progression surveillance. This proof of concept study accomplishes both at once, population-scale diagnostics and variants' tracking. This strategy relies on (1) the use of the portable MinION DNA sequencer; (2) a DNA barcoding and a SPIKE gene-centered variant's tracking, increasing the number of candidates per assay; and (3) a real-time diagnostics and variant's tracking monitoring thanks to our software RETIVAD. This strategy represents an optimal solution for addressing the current needs on SARS-CoV-2 progression surveillance, notably due to its affordable implementation, allowing its implantation even in remote places over the world.

Inferring biologically relevant molecular tissue substructures by agglomerative clustering of digitized spatial transcriptomes with multilayer.

Spatially resolved transcriptomics (SrT) can investigate organ or tissue architecture from the angle of gene programs that define their molecular complexity. However, computational methods to analyze SrT data underexploit their spatial signature. Inspired by contextual pixel classification strategies applied to image analysis, we developed MULTILAYER to stratify maps into functionally relevant molecular substructures. MULTILAYER applies agglomerative clustering within contiguous locally defined transcriptomes (gene expression elements or "gexels") combined with community detection methods for graphical partitioning. MULTILAYER resolves molecular tissue substructures within a variety of SrT data with superior performance to commonly used dimensionality reduction strategies and still detects differentially expressed genes on par with existing methods. MULTILAYER can process high-resolution as well as multiple SrT data in a comparative mode, anticipating future needs in the field. MULTILAYER provides a digital image perspective for SrT analysis and opens the door to contextual gexel classification strategies for developing self-supervised molecular diagnosis solutions. A record of this paper's transparent peer review process is included in the supplemental information.

A Core Transcription Regulatory Circuitry Defining Microglia Cell Identity Inferred from the Reanalysis of Multiple Human Microglia Differentiation Protocols.

Microglia, the immune cells in the brain involved in both homeostasis and injury/infection control, play a predominant role in neurodegenerative diseases. In vivo studies on microglia are limited due to the requirement of surgical intervention, which can lead to the destruction of the tissues. Over the last few years, multiple protocols-presenting a variety of strategies-have described microglia differentiation issued from human pluripotent stem cells. Herein, we have reanalyzed the transcriptomes released on six different microglia differentiation protocols and revealed a consensus core of master transcription regulatory circuitry defining microglia identity. Furthermore, we have discussed the major divergencies among the studied protocols and have provided suggestions to further enhance microglia differentiation assays.

Protocol for using MULTILAYER to reveal molecular tissue substructures from digitized spatial transcriptomes.

Spatially resolved transcriptomics (SrT) allow researchers to explore organ/tissue architecture from the angle of the gene programs involved in their molecular complexity. Here, we describe the use of MULTILAYER to reveal molecular tissue substructures from the analysis of localized transcriptomes (defined as gexels). MULTILAYER can process low- and high-resolution SrT data but also perform comparative analyses within multiple SrT readouts. For complete details on the use and execution of this protocol, please refer to Moehlin et al., 2021.

A comprehensive resource for retrieving, visualizing, and integrating functional genomics data.

The enormous amount of freely accessible functional genomics data is an invaluable resource for interrogating the biological function of multiple DNA-interacting players and chromatin modifications by large-scale comparative analyses. However, in practice, interrogating large collections of public data requires major efforts for (i) reprocessing available raw reads, (ii) incorporating quality assessments to exclude artefactual and low-quality data, and (iii) processing data by using high-performance computation. Here, we present qcGenomics, a user-friendly online resource for ultrafast retrieval, visualization, and comparative analysis of tens of thousands of genomics datasets to gain new functional insight from global or focused multidimensional data integration.

RARß Agonist Drug (C286) Demonstrates Efficacy in a Pre-clinical Neuropathic Pain Model Restoring Multiple Pathways via DNA Repair Mechanisms.

Neuropathic pain (NP) is associated with profound gene expression alterations within the nociceptive system. DNA mechanisms, such as epigenetic remodeling and repair pathways have been implicated in NP. Here we have used a rat model of peripheral nerve injury to study the effect of a recently developed RARß agonist, C286, currently under clinical research, in NP. A 4-week treatment initiated 2 days after the injury normalized pain sensation. Genome-wide and pathway enrichment analysis showed that multiple mechanisms persistently altered in the spinal cord were restored to preinjury levels by the agonist. Concomitant upregulation of DNA repair proteins, ATM and BRCA1, the latter being required for C286-mediated pain modulation, suggests that early DNA repair may be important to prevent phenotypic epigenetic imprints in NP. Thus, C286 is a promising drug candidate for neuropathic pain and DNA repair mechanisms may be useful therapeutic targets to explore.