Immunoresponsive gene 1 (IRG1) is one of the highest induced genes in macrophages under pro-inflammatory conditions. Its function has been recently described: it codes for immune-responsive gene 1 protein/cis-aconitic acid decarboxylase (IRG1/CAD), an enzyme catalysing the production of itaconic acid from cis-aconitic acid, a tricarboxylic acid (TCA) cycle intermediate. Itaconic acid possesses specific antimicrobial properties inhibiting isocitrate lyase, the first enzyme of the glyoxylate shunt, an anaplerotic pathway that bypasses the TCA cycle and enables bacteria to survive on limited carbon conditions. To elucidate the mechanisms underlying itaconic acid production through IRG1 induction in macrophages, we examined the transcriptional regulation of IRG1. To this end, we studied IRG1 expression in human immune cells under different inflammatory stimuli, such as TNFα and IFNγ, in addition to lipopolysaccharides. Under these conditions, as previously shown in mouse macrophages, IRG1/CAD accumulates in mitochondria. Furthermore, using literature information and transcription factor prediction models, we re-constructed raw gene regulatory networks (GRNs) for IRG1 in mouse and human macrophages. We further implemented a contextualization algorithm that relies on genome-wide gene expression data to infer putative cell type-specific gene regulatory interactions in mouse and human macrophages, which allowed us to predict potential transcriptional regulators of IRG1. Among the computationally identified regulators, siRNA-mediated gene silencing of interferon regulatory factor 1 (IRF1) in macrophages significantly decreased the expression of IRG1/CAD at the gene and protein level, which correlated with a reduced production of itaconic acid. Using a synergistic approach of both computational and experimental methods, we here shed more light on the transcriptional machinery of IRG1 expression and could pave the way to therapeutic approaches targeting itaconic acid levels.

Neural stem cells (NSCs) in the adult mouse hippocampus occur in a specific neurogenic niche, where a multitude of extracellular signaling molecules converges to regulate NSC proliferation as well as fate and functional integration. However, the underlying mechanisms how NSCs react to extrinsic signals and convert them to intracellular responses still remains elusive. NSCs contain a functional endocannabinoid system, including the cannabinoid type-1 receptor (CB1). To decipher whether CB1 regulates adult neurogenesis directly or indirectly in vivo, we performed NSC-specific conditional inactivation of CB1 by using triple-transgenic mice. Here, we show that lack of CB1 in NSCs is sufficient to decrease proliferation of the stem cell pool, which consequently leads to a reduction in the number of newborn neurons. Furthermore, neuronal differentiation was compromised at the level of dendritic maturation pointing towards a postsynaptic role of CB1 in vivo. Deteriorated neurogenesis in NSC-specific CB1 knock-outs additionally resulted in reduced long-term potentiation in the hippocampal formation. The observed cellular and physiological alterations led to decreased short-term spatial memory and increased depression-like behavior. These results demonstrate that CB1 expressed in NSCs and their progeny controls neurogenesis in adult mice to regulate the NSC stem cell pool, dendritic morphology, activity-dependent plasticity, and behavior.

The function of somatic stem cells declines with age. Understanding the molecular underpinnings of this decline is key to counteract age-related disease. Here, we report a dramatic drop in the neural stem cells (NSCs) number in the aging murine brain. We find that this smaller stem cell reservoir is protected from full depletion by an increase in quiescence that makes old NSCs more resistant to regenerate the injured brain. Once activated, however, young and old NSCs show similar proliferation and differentiation capacity. Single-cell transcriptomics of NSCs indicate that aging changes NSCs minimally. In the aging brain, niche-derived inflammatory signals and the Wnt antagonist sFRP5 induce quiescence. Indeed, intervention to neutralize them increases activation of old NSCs during homeostasis and following injury. Our study identifies quiescence as a key feature of old NSCs imposed by the niche and uncovers ways to activate NSCs to repair the aging brain.

The balance between stem cell maintenance and differentiation has been proposed to depend on antagonizing ubiquitination and deubiquitination reactions of key stem cell transcription factors (SCTFs) mediated by pairs of E3 ubiquitin ligases and deubiquitinating enzymes. Accordingly, increased ubiquitination results in proteasomal degradation of the SCTF, thereby inducing cellular differentiation, whereas increased deubiquitination stabilizes the SCTF, leading to maintenance of the stem cell fate. In neural stem cells, one of the key SCTFs is c-Myc. Previously, it has been shown that c-Myc is ubiquitinated by the E3 ligase TRIM32, thereby targeting c-Myc for proteasomal degradation and inducing neuronal differentiation. Accordingly, TRIM32 becomes upregulated during adult neurogenesis. This upregulation is accompanied by subcellular translocation of TRIM32 from the cytoplasm of neuroblasts to the nucleus of neurons. However, we observed that a subpopulation of proliferative type C cells already contains nuclear TRIM32. As these cells do not undergo neuronal differentiation, despite containing TRIM32 in the nucleus, where it can ubiquitinate c-Myc, we hypothesize that antagonizing factors, specifically deubiquitinating enzymes, are present in these particular cells. Here we show that TRIM32 associates with the deubiquitination enzyme USP7, which previously has been implicated in neural stem cell maintenance. USP7 and TRIM32 were found to exhibit a dynamic and partially overlapping expression pattern during neuronal differentiation both in vitro and in vivo. Most importantly, we are able to demonstrate that USP7 deubiquitinates and thereby stabilizes c-Myc and that this function is required to maintain neural stem cell fate. Accordingly, we propose the balanced ubiquitination and deubiquitination of c-Myc by TRIM32 and USP7 as a novel mechanism for stem cell fate determination.

Cellular differentiation is a complex process where a less specialized cell evolves into a more specialized cell. Despite the increasing research effort, identification of cell-fate determinants (transcription factors (TFs) determining cell fates during differentiation) still remains a challenge, especially when closely related cell types from a common progenitor are considered. Here, we develop SeesawPred, a web application that, based on a gene regulatory network (GRN) model of cell differentiation, can computationally predict cell-fate determinants from transcriptomics data. Unlike previous approaches, it allows the user to upload gene expression data and does not rely on pre-compiled reference data sets, enabling its application to novel differentiation systems. SeesawPred correctly predicted known cell-fate determinants on various cell differentiation examples in both mouse and human, and also performed better compared to state-of-the-art methods. The application is freely available for academic, non-profit use at http://seesaw.lcsb.uni.lu.

Human FOXP3(+)CD25(+)CD4(+) regulatory T cells (Tregs) are essential to the maintenance of immune homeostasis. Several genes are known to be important for murine Tregs, but for human Tregs the genes and underlying molecular networks controlling the suppressor function still largely remain unclear. Here, we describe a strategy to identify the key genes directly from an undirected correlation network which we reconstruct from a very high time-resolution (HTR) transcriptome during the activation of human Tregs/CD4(+) T-effector cells. We show that a predicted top-ranked new key gene PLAU (the plasminogen activator urokinase) is important for the suppressor function of both human and murine Tregs. Further analysis unveils that PLAU is particularly important for memory Tregs and that PLAU mediates Treg suppressor function via STAT5 and ERK signaling pathways. Our study demonstrates the potential for identifying novel key genes for complex dynamic biological processes using a network strategy based on HTR data, and reveals a critical role for PLAU in Treg suppressor function.

Allosteric communications are vital for cellular signaling. Here we explore a relationship between protein architectural organization and shortcuts in signaling pathways.

A plethora of computational approaches have been proposed for reconstructing gene regulatory networks (GRNs) from gene expression data. However, gene regulatory processes are often too complex to predict from the transcriptome alone. Here, we present a computational method, Moni, that systematically integrates epigenetics, transcriptomics, and protein-protein interactions to reconstruct GRNs among core transcription factors and their co-factors governing cell identity. We applied Moni to 57 datasets of human cell types and lines and demonstrate that it can accurately infer GRNs, thereby outperforming state-of-the-art methods.

Single-cell RNA-sequencing is increasingly employed to characterize disease or ageing cell subpopulation phenotypes. Despite exponential increase in data generation, systematic identification of key regulatory factors for controlling cellular phenotype to enable cell rejuvenation in disease or ageing remains a challenge. Here, we present SigHotSpotter, a computational tool to predict hotspots of signaling pathways responsible for the stable maintenance of cell subpopulation phenotypes, by integrating signaling and transcriptional networks. Targeted perturbation of these signaling hotspots can enable precise control of cell subpopulation phenotypes. SigHotSpotter correctly predicts the signaling hotspots with known experimental validations in different cellular systems. The tool is simple, user-friendly and is available as web-server or as stand-alone software. We believe SigHotSpotter will serve as a general purpose tool for the systematic prediction of signaling hotspots based on single-cell RNA-seq data, and potentiate novel cell rejuvenation strategies in the context of disease and ageing.

The functional specialization of cell types arises during development and is shaped by cell-cell communication networks determining a distribution of functional cell states that are collectively important for tissue functioning. However, the identification of these tissue-specific functional cell states remains challenging. Although a plethora of computational approaches have been successful in detecting cell types and subtypes, they fail in resolving tissue-specific functional cell states. To address this issue, we present FunRes, a computational method designed for the identification of functional cell states. FunRes relies on scRNA-seq data of a tissue to initially reconstruct the functional cell-cell communication network, which is leveraged for partitioning each cell type into functional cell states. We applied FunRes to 177 cell types in 10 different tissues and demonstrated that the detected states correspond to known functional cell states of various cell types, which cannot be recapitulated by existing computational tools. Finally, we characterize emerging and vanishing functional cell states in aging and disease, and demonstrate their involvement in key tissue functions. Thus, we believe that FunRes will be of great utility in the characterization of the functional landscape of cell types and the identification of dysfunctional cell states in aging and disease.

Asymmetric neuronal expansion is thought to drive evolutionary transitions between lissencephalic and gyrencephalic cerebral cortices. We report that Neurog2 and Ascl1 proneural genes together sustain neurogenic continuity and lissencephaly in rodent cortices. Using transgenic reporter mice and human cerebral organoids, we found that Neurog2 and Ascl1 expression defines a continuum of four lineage-biased neural progenitor cell (NPC) pools. Double+ NPCs, at the hierarchical apex, are least lineage restricted due to Neurog2-Ascl1 cross-repression and display unique features of multipotency (more open chromatin, complex gene regulatory network, G2 pausing). Strikingly, selectively eliminating double+ NPCs by crossing Neurog2-Ascl1 split-Cre mice with diphtheria toxin-dependent "deleter" strains locally disrupts Notch signaling, perturbs neurogenic symmetry, and triggers cortical folding. In support of our discovery that double+ NPCs are Notch-ligand-expressing "niche" cells that control neurogenic periodicity and cortical folding, NEUROG2, ASCL1, and HES1 transcript distribution is modular (adjacent high/low zones) in gyrencephalic macaque cortices, prefiguring future folds.

Tissue regeneration requires coordination between resident stem cells and local niche cells1,2. Here we identify that senescent cells are integral components of the skeletal muscle regenerative niche that repress regeneration at all stages of life. The technical limitation of senescent-cell scarcity3 was overcome by combining single-cell transcriptomics and a senescent-cell enrichment sorting protocol. We identified and isolated different senescent cell types from damaged muscles of young and old mice. Deeper transcriptome, chromatin and pathway analyses revealed conservation of cell identity traits as well as two universal senescence hallmarks (inflammation and fibrosis) across cell type, regeneration time and ageing. Senescent cells create an aged-like inflamed niche that mirrors inflammation associated with ageing (inflammageing4) and arrests stem cell proliferation and regeneration. Reducing the burden of senescent cells, or reducing their inflammatory secretome through CD36 neutralization, accelerates regeneration in young and old mice. By contrast, transplantation of senescent cells delays regeneration. Our results provide a technique for isolating in vivo senescent cells, define a senescence blueprint for muscle, and uncover unproductive functional interactions between senescent cells and stem cells in regenerative niches that can be overcome. As senescent cells also accumulate in human muscles, our findings open potential paths for improving muscle repair throughout life.

The different actions of abscisic acid (ABA) in the aboveground and belowground parts of plants suggest the existence of a distinct perception mechanism between these organs. Although characterization of the soluble ABA receptors PYR1/PYL/RCAR as well as core signaling components has greatly advanced our understanding of ABA perception, signal transduction, and responses, the environment-dependent organ-specific sensitivity of plants to ABA is less well understood.

With defined culture protocol, human embryonic stem cells (hESCs) are able to generate cardiomyocytes in vitro, therefore providing a great model for human heart development, and holding great potential for cardiac disease therapies. In this study, we successfully generated a highly pure population of human cardiomyocytes (hCMs) (>95% cTnT(+)) from hESC line, which enabled us to identify and characterize an hCM-specific signature, at both the gene expression and DNA methylation levels. Gene functional association network and gene-disease network analyses of these hCM-enriched genes provide new insights into the mechanisms of hCM transcriptional regulation, and stand as an informative and rich resource for investigating cardiac gene functions and disease mechanisms. Moreover, we show that cardiac-structural genes and cardiac-transcription factors have distinct epigenetic mechanisms to regulate their gene expression, providing a better understanding of how the epigenetic machinery coordinates to regulate gene expression in different cell types.

MicroRNAs (miRNA) are a class of small, highly conserved noncoding RNA molecules, which contain 18-28 nucleotides and are involved in the regulation of gene expression. It has been proved that microRNAs play a very important role in several key cellular processes, such as cell differentiation, cell cycle progression, and apoptosis, as well as in autoimmune disease. One recently identified miRNA, miR-708-5p, has been demonstrated to have profound roles in suppressing oncogenesis in different types of tumors. However, the role of miR-708-5p in rheumatoid arthritis (RA) remains to be fully elucidated. Therefore, in this study, we are aiming to identify the role of miR-708-5p in RA.

The juvenile form of neuronal ceroid Lipofuscinosis (JNCL) is the most common form within this group of rare lysosomal storage disorders, causing pediatric neurodegeneration. The genetic disorder, which is caused by recessive mutations affecting the CLN3 gene, features progressive vision loss, cognitive and motor decline and other psychiatric conditions, seizure episodes, leading to premature death. Animal models have traditionally aid the understanding of the disease mechanisms and pathology and are very relevant for biomarker research and therapeutic testing. Nevertheless, there is a need for establishing reliable and predictive human cellular models to study the disease. Since patient material, particularly from children, is scarce and difficult to obtain, we generated an engineered a CLN3-mutant isogenic human induced pluripotent stem cell (hiPSC) line carrying the c.1054C → T pathologic variant, using state of the art CRISPR/Cas9 technology. To prove the suitability of the isogenic pair to model JNCL, we screened for disease-specific phenotypes in non-neuronal two-dimensional cell culture models as well as in cerebral brain organoids. Our data demonstrates that the sole introduction of the pathogenic variant gives rise to classical hallmarks of JNCL in vitro. Additionally, we discovered an alteration of the splicing caused by this particular mutation. Next, we derived cerebral organoids and used them as a neurodevelopmental model to study the particular effects of the CLN3Q352X mutation during brain formation in the disease context. About half of the mutation -carrying cerebral organoids completely failed to develop normally. The other half, which escaped this severe defect were used for the analysis of more subtle alterations. In these escapers, whole-transcriptome analysis demonstrated early disease signatures, affecting pathways related to development, corticogenesis and synapses. Complementary metabolomics analysis confirmed decreased levels of cerebral tissue metabolites, some particularly relevant for synapse formation and neurotransmission, such as gamma-amino butyric acid (GABA). Our data suggests that a mutation in CLN3 severely affects brain development. Furthermore, before disease onset, disease -associated neurodevelopmental changes, particular concerning synapse formation and function, occur.

Engineering of biomaterials with specific biological properties has gained momentum as a means to control stem cell behavior. Here, we address the effect of bifunctionalized hydrogels comprising polylysine (PL) and a 19-mer peptide containing the laminin motif IKVAV (IKVAV) on embryonic and adult neuronal progenitor cells under different stiffness regimes. Neuronal differentiation of embryonic and adult neural progenitors was accelerated by adjusting the gel stiffness to 2 kPa and 20 kPa, respectively. While gels containing IKVAV or PL alone failed to support long-term cell adhesion, in bifunctional gels, IKVAV synergized with PL to promote differentiation and formation of focal adhesions containing β1-integrin in embryonic cortical neurons. Furthermore, in adult neural stem cell culture, bifunctionalized gels promoted neurogenesis via the expansion of neurogenic clones. These data highlight the potential of synthetic matrices to steer stem and progenitor cell behavior via defined mechano-adhesive properties.

The RNA polymerase II-associated protein 1 (RPAP1) is conserved across metazoa and required for stem cell differentiation in plants; however, very little is known about its mechanism of action or its role in mammalian cells. Here, we report that RPAP1 is essential for the expression of cell identity genes and for cell viability. Depletion of RPAP1 triggers cell de-differentiation, facilitates reprogramming toward pluripotency, and impairs differentiation. Mechanistically, we show that RPAP1 is essential for the interaction between RNA polymerase II (RNA Pol II) and Mediator, as well as for the recruitment of important regulators, such as the Mediator-specific RNA Pol II factor Gdown1 and the C-terminal domain (CTD) phosphatase RPAP2. In agreement, depletion of RPAP1 diminishes the loading of total and Ser5-phosphorylated RNA Pol II on many genes, with super-enhancer-driven genes among the most significantly downregulated. We conclude that Mediator/RPAP1/RNA Pol II is an ancient module, conserved from plants to mammals, critical for establishing and maintaining cell identity.

Prior knowledge of perturbation data can significantly assist in inferring the relationship between chemical perturbations and their specific transcriptional response. However, current databases mostly contain cancer cell lines, which are unsuitable for the aforementioned inference in non-cancer cells, such as cells related to non-cancer disease, immunology and aging. Here, we present ChemPert (https://chempert.uni.lu/), a database consisting of 82 270 transcriptional signatures in response to 2566 unique perturbagens (drugs, small molecules and protein ligands) across 167 non-cancer cell types, as well as the protein targets of 57 818 perturbagens. In addition, we develop a computational tool that leverages the non-cancer cell datasets, which enables more accurate predictions of perturbation responses and drugs in non-cancer cells compared to those based onto cancer databases. In particular, ChemPert correctly predicted drug effects for treating hepatitis and novel drugs for osteoarthritis. The ChemPert web interface is user-friendly and allows easy access of the entire datasets and the computational tool, providing valuable resources for both experimental researchers who wish to find datasets relevant to their research and computational researchers who need comprehensive non-cancer perturbation transcriptomics datasets for developing novel algorithms. Overall, ChemPert will facilitate future in silico compound screening for non-cancer cells.

Organogenesis involves integration of diverse cell types; dysregulation of cell-type-specific gene networks results in birth defects, which affect 5% of live births. Congenital heart defects are the most common malformations, and result from disruption of discrete subsets of cardiac progenitor cells1, but the transcriptional changes in individual progenitors that lead to organ-level defects remain unknown. Here we used single-cell RNA sequencing to interrogate early cardiac progenitor cells as they become specified during normal and abnormal cardiogenesis, revealing how dysregulation of specific cellular subpopulations has catastrophic consequences. A network-based computational method for single-cell RNA-sequencing analysis that predicts lineage-specifying transcription factors2,3 identified Hand2 as a specifier of outflow tract cells but not right ventricular cells, despite the failure of right ventricular formation in Hand2-null mice4. Temporal single-cell-transcriptome analysis of Hand2-null embryos revealed failure of outflow tract myocardium specification, whereas right ventricular myocardium was specified but failed to properly differentiate and migrate. Loss of Hand2 also led to dysregulation of retinoic acid signalling and disruption of anterior-posterior patterning of cardiac progenitors. This work reveals transcriptional determinants that specify fate and differentiation in individual cardiac progenitor cells, and exposes mechanisms of disrupted cardiac development at single-cell resolution, providing a framework for investigating congenital heart defects.