A Multi-step Transcriptional and Chromatin State Cascade Underlies Motor Neuron Programming from Embryonic Stem Cells

Optimization of a Cas12a-Driven Synthetic Gene Regulatory Network System

Gene regulatory networks, which control gene expression patterns in development and in response to stimuli, use regulatory logic modules to coordinate inputs and outputs. One example of a regulatory logic module is the gene regulatory cascade (GRC), where a series of transcription factor genes turn on in order. Synthetic biologists have derived artificial systems that encode regulatory rules, including GRCs. Furthermore, the development of single-cell approaches has enabled the discovery of gene regulatory modules in a variety of experimental settings. However, the tools available for validating these observations remain limited. Based on a synthetic GRC using DNA cutting-defective Cas9 (dCas9), we designed and implemented an alternative synthetic GRC utilizing DNA cutting-defective Cas12a (dCas12a). Comparing the ability of these two systems to express a fluorescent reporter, the dCas9 system was initially more active, while the dCas12a system was more streamlined. Investigating the influence of individual components of the systems identified nuclear localization as a major driver of differences in activity. Improving nuclear localization for the dCas12a system resulted in 1.5-fold more reporter-positive cells and a 15-fold increase in reporter intensity relative to the dCas9 system. We call this optimized system the "Synthetic Gene Regulatory Network" (SGRN, pronounced "sojourn").

Modeling ALS with Patient-Derived iPSCs: Recent Advances and Future Potentials

Amyotrophic lateral sclerosis (ALS) is a terminal complex neurodegenerative disease, with 10-15% of cases being familial and the majority being sporadic with no known cause. There are no animal models for the 85-90% of sporadic ALS cases. More creative, sophisticated models of ALS disease are required to unravel the mysteries of this complicated disease. While ALS patients urgently require new medications and treatments, suitable preclinical in vitro models for drug screening are lacking. Therefore, human-derived induced pluripotent stem cell (hiPSC) technology offers the opportunity to model diverse and unreachable cell types in a culture dish. In this review, we focus on recent hiPSC-derived ALS neuronal and non-neuronal models to examine the research progress of current ALS 2D monocultures, co-cultures, and more complex 3D-model organoids. Despite the challenges inherent to hiPSC-based models, their application to preclinical drug studies is enormous.

Recruitment of homodimeric proneural factors by conserved CAT-CAT E-boxes drives major epigenetic reconfiguration in cortical neurogenesis

Proneural factors of the basic helix-loop-helix family coordinate neurogenesis and neurodifferentiation. Among them, NEUROG2 and NEUROD2 subsequently act to specify neurons of the glutamatergic lineage. Disruption of these factors, their target genes and binding DNA motifs has been linked to various neuropsychiatric disorders. Proneural factors bind to specific DNA motifs called E-boxes (hexanucleotides of the form CANNTG, composed of two CAN half sites on opposed strands). While corticogenesis heavily relies on E-box activity, the collaboration of proneural factors on different E-box types and their chromatin remodeling mechanisms remain largely unknown. Here, we conducted a comprehensive analysis using chromatin immunoprecipitation followed by sequencing (ChIP-seq) data for NEUROG2 and NEUROD2, along with time-matched single-cell RNA-seq, ATAC-seq and DNA methylation data from the developing mouse cortex. Our findings show that these factors are highly enriched in transiently active genomic regions during intermediate stages of neuronal differentiation. Although they primarily bind CAG-containing E-boxes, their binding in dynamic regions is notably enriched in CAT-CAT E-boxes (i.e. CATATG, denoted as 5'3' half sites for dimers), which undergo significant DNA demethylation and exhibit the highest levels of evolutionary constraint. Aided by HT-SELEX data reanalysis, structural modeling and DNA footprinting, we propose that these proneural factors exert maximal chromatin remodeling influence during intermediate stages of neurogenesis by binding as homodimers to CAT-CAT motifs. This study provides an in-depth integrative analysis of the dynamic regulation of E-boxes during neuronal development, enhancing our understanding of the mechanisms underlying the binding specificity of critical proneural factors.

Multimodal Hox5 activity generates motor neuron diversity

Motor neurons (MNs) are the final output of circuits driving fundamental behaviors, such as respiration and locomotion. Hox proteins are essential in generating the MN diversity required for accomplishing these functions, but the transcriptional mechanisms that enable Hox paralogs to assign distinct MN subtype identities despite their promiscuous DNA binding motif are not well understood. Here we show that Hoxa5 modifies chromatin accessibility in all mouse spinal cervical MN subtypes and engages TALE co-factors to directly bind and regulate subtype-specific genes. We identify a paralog-specific interaction of Hoxa5 with the phrenic MN-specific transcription factor Scip and show that heterologous expression of Hoxa5 and Scip is sufficient to suppress limb-innervating MN identity. We also demonstrate that phrenic MN identity is stable after Hoxa5 downregulation and identify Klf proteins as potential regulators of phrenic MN maintenance. Our data identify multiple modes of Hoxa5 action that converge to induce and maintain MN identity.

Transcriptional Dynamics and Key Regulators of Adipogenesis in Mouse Embryonic Stem Cells: Insights from Robust Rank Aggregation Analysis

Embryonic stem cells are crucial for studying developmental biology due to their self-renewal and pluripotency capabilities. This research investigates the differentiation of mouse ESCs into adipocytes, offering insights into obesity and metabolic disorders. Using a monolayer differentiation approach over 30 days, lipid accumulation and adipogenic markers, such as Cebpb, Pparg, and Fabp4, confirmed successful differentiation. RNA sequencing revealed extensive transcriptional changes, with over 15,000 differentially expressed genes linked to transcription regulation, cell cycle, and DNA repair. This study utilized Robust Rank Aggregation to identify critical regulatory genes like PPARG, CEBPA, and EP300. Network analysis further highlighted Atf5, Ccnd1, and Nr4a1 as potential key players in adipogenesis and its mature state, validated through RT-PCR. While key adipogenic factors showed plateaued expression levels, suggesting early differentiation events, this study underscores the value of ESCs in modeling adipogenesis. These findings contribute to our understanding of adipocyte differentiation and have significant implications for therapeutic strategies targeting metabolic diseases.

Lhx4 surpasses its paralog Lhx3 in promoting the differentiation of spinal V2a interneurons

Paralog factors are considered to ensure the robustness of biological processes by providing redundant activity in cells where they are co-expressed. However, the specific contribution of each factor is frequently underestimated. In the developing spinal cord, multiple families of transcription factors successively contribute to differentiate an initially homogenous population of neural progenitors into a myriad of neuronal subsets with distinct molecular, morphological, and functional characteristics. The LIM-homeodomain transcription factors Lhx3, Lhx4, Isl1 and Isl2 promote the segregation and differentiation of spinal motor neurons and V2 interneurons. Based on their high sequence identity and their similar distribution, the Lhx3 and Lhx4 paralogs are considered to contribute similarly to these processes. However, the specific contribution of Lhx4 has never been studied. Here, we provide evidence that Lhx3 and Lhx4 are present in the same cell populations during spinal cord development. Similarly to Lhx3, Lhx4 can form multiproteic complexes with Isl1 or Isl2 and the nuclear LIM interactor NLI. Lhx4 can stimulate a V2-specific enhancer more efficiently than Lhx3 and surpasses Lhx3 in promoting the differentiation of V2a interneurons in chicken embryo electroporation experiments. Finally, Lhx4 inactivation in mice results in alterations of differentiation of the V2a subpopulation, but not of motor neuron production, suggesting that Lhx4 plays unique roles in V2a differentiation that are not compensated by the presence of Lhx3. Thus, Lhx4 could be the major LIM-HD factor involved in V2a interneuron differentiation during spinal cord development and should be considered for in vitro differentiation of spinal neuronal populations.

Uncovering the dynamics and consequences of RNA isoform changes during neuronal differentiation

Static gene expression programs have been extensively characterized in stem cells and mature human cells. However, the dynamics of RNA isoform changes upon cell-state-transitions during cell differentiation, the determinants and functional consequences have largely remained unclear. Here, we established an improved model for human neurogenesis in vitro that is amenable for systems-wide analyses of gene expression. Our multi-omics analysis reveals that the pronounced alterations in cell morphology correlate strongly with widespread changes in RNA isoform expression. Our approach identifies thousands of new RNA isoforms that are expressed at distinct differentiation stages. RNA isoforms mainly arise from exon skipping and the alternative usage of transcription start and polyadenylation sites during human neurogenesis. The transcript isoform changes can remodel the identity and functions of protein isoforms. Finally, our study identifies a set of RNA binding proteins as a potential determinant of differentiation stage-specific global isoform changes. This work supports the view of regulated isoform changes that underlie state-transitions during neurogenesis.

The multi-lineage transcription factor ISL1 controls cardiomyocyte cell fate through interaction with NKX2.5

Congenital heart disease often arises from perturbations of transcription factors (TFs) that guide cardiac development. ISLET1 (ISL1) is a TF that influences early cardiac cell fate, as well as differentiation of other cell types including motor neuron progenitors (MNPs) and pancreatic islet cells. While lineage specificity of ISL1 function is likely achieved through combinatorial interactions, its essential cardiac interacting partners are unknown. By assaying ISL1 genomic occupancy in human induced pluripotent stem cell-derived cardiac progenitors (CPs) or MNPs and leveraging the deep learning approach BPNet, we identified motifs of other TFs that predicted ISL1 occupancy in each lineage, with NKX2.5 and GATA motifs being most closely associated to ISL1 in CPs. Experimentally, nearly two-thirds of ISL1-bound loci were co-occupied by NKX2.5 and/or GATA4. Removal of NKX2.5 from CPs led to widespread ISL1 redistribution, and overexpression of NKX2.5 in MNPs led to ISL1 occupancy of CP-specific loci. These results reveal how ISL1 guides lineage choices through a combinatorial code that dictates genomic occupancy and transcription.

CRISPR/Cas9 Directed Reprogramming of iPSC for Accelerated Motor Neuron Differentiation Leads to Dysregulation of Neuronal Fate Patterning and Function

Neurodegeneration causes a significant disease burden and there are few therapeutic interventions available for reversing or slowing the disease progression. Induced pluripotent stem cells (iPSCs) hold significant potential since they are sourced from adult tissue and have the capacity to be differentiated into numerous cell lineages, including motor neurons. This differentiation process traditionally relies on cell lineage patterning factors to be supplied in the differentiation media. Genetic engineering of iPSC with the introduction of recombinant master regulators of motor neuron (MN) differentiation has the potential to shorten and streamline cell developmental programs. We have established stable iPSC cell lines with transient induction of exogenous LHX3 and ISL1 from the Tet-activator regulatory region and have demonstrated that induction of the transgenes is not sufficient for the development of mature MNs in the absence of neuron patterning factors. Comparative global transcriptome analysis of MN development from native and Lhx-ISL1 modified iPSC cultures demonstrated that the genetic manipulation helped to streamline the neuronal patterning process. However, leaky gene expression of the exogenous MN master regulators in iPSC resulted in the premature activation of genetic pathways characteristic of the mature MN function. Dysregulation of metabolic and regulatory pathways within the developmental process affected the MN electrophysiological responses.

Joint sequence & chromatin neural networks characterize the differential abilities of Forkhead transcription factors to engage inaccessible chromatin

The DNA-binding activities of transcription factors (TFs) are influenced by both intrinsic sequence preferences and extrinsic interactions with cell-specific chromatin landscapes and other regulatory proteins. Disentangling the roles of these binding determinants remains challenging. For example, the FoxA subfamily of Forkhead domain (Fox) TFs are known pioneer factors that can bind to relatively inaccessible sites during development. Yet FoxA TF binding also varies across cell types, pointing to a combination of intrinsic and extrinsic forces guiding their binding. While other Forkhead domain TFs are often assumed to have pioneering abilities, how sequence and chromatin features influence the binding of related Fox TFs has not been systematically characterized. Here, we present a principled approach to compare the relative contributions of intrinsic DNA sequence preference and cell-specific chromatin environments to a TF's DNA-binding activities. We apply our approach to investigate how a selection of Fox TFs (FoxA1, FoxC1, FoxG1, FoxL2, and FoxP3) vary in their binding specificity. We over-express the selected Fox TFs in mouse embryonic stem cells, which offer a platform to contrast each TF's binding activity within the same preexisting chromatin background. By applying a convolutional neural network to interpret the Fox TF binding patterns, we evaluate how sequence and preexisting chromatin features jointly contribute to induced TF binding. We demonstrate that Fox TFs bind different DNA targets, and drive differential gene expression patterns, even when induced in identical chromatin settings. Despite the association between Forkhead domains and pioneering activities, the selected Fox TFs display a wide range of affinities for preexiting chromatin states. Using sequence and chromatin feature attribution techniques to interpret the neural network predictions, we show that differential sequence preferences combined with differential abilities to engage relatively inaccessible chromatin together explain Fox TF binding patterns at individual sites and genome-wide.

CUT homeobox genes: transcriptional regulation of neuronal specification and beyond

CUT homeobox genes represent a captivating gene class fulfilling critical functions in the development and maintenance of multiple cell types across a wide range of organisms. They belong to the larger group of homeobox genes, which encode transcription factors responsible for regulating gene expression patterns during development. CUT homeobox genes exhibit two distinct and conserved DNA binding domains, a homeodomain accompanied by one or more CUT domains. Numerous studies have shown the involvement of CUT homeobox genes in diverse developmental processes such as body axis formation, organogenesis, tissue patterning and neuronal specification. They govern these processes by exerting control over gene expression through their transcriptional regulatory activities, which they accomplish by a combination of classic and unconventional interactions with the DNA. Intriguingly, apart from their roles as transcriptional regulators, they also serve as accessory factors in DNA repair pathways through protein-protein interactions. They are highly conserved across species, highlighting their fundamental importance in developmental biology. Remarkably, evolutionary analysis has revealed that CUT homeobox genes have experienced an extraordinary degree of rearrangements and diversification compared to other classes of homeobox genes, including the emergence of a novel gene family in vertebrates. Investigating the functions and regulatory networks of CUT homeobox genes provides significant understanding into the molecular mechanisms underlying embryonic development and tissue homeostasis. Furthermore, aberrant expression or mutations in CUT homeobox genes have been associated with various human diseases, highlighting their relevance beyond developmental processes. This review will overview the well known roles of CUT homeobox genes in nervous system development, as well as their functions in other tissues across phylogeny.

Ascl1 and Ngn2 convert mouse embryonic stem cells to neurons via functionally distinct paths

Ascl1 and Ngn2, closely related proneural transcription factors, are able to convert mouse embryonic stem cells into induced neurons. Despite their similarities, these factors elicit only partially overlapping transcriptional programs, and it remains unknown whether cells are converted via distinct mechanisms. Here we show that Ascl1 and Ngn2 induce mutually exclusive side populations by binding and activating distinct lineage drivers. Furthermore, Ascl1 rapidly dismantles the pluripotency network and installs neuronal and trophoblast cell fates, while Ngn2 generates a neural stem cell-like intermediate supported by incomplete shutdown of the pluripotency network. Using CRISPR-Cas9 knockout screening, we find that Ascl1 relies more on factors regulating pluripotency and the cell cycle, such as Tcf7l1. In the absence of Tcf7l1, Ascl1 still represses core pluripotency genes but fails to exit the cell cycle. However, overexpression of Cdkn1c induces cell cycle exit and restores the generation of neurons. These findings highlight that cell type conversion can occur through two distinct mechanistic paths, even when induced by closely related transcription factors.

The Troyer syndrome protein spartin mediates selective autophagy of lipid droplets

Lipid droplets (LDs) are crucial organelles for energy storage and lipid homeostasis. Autophagy of LDs is an important pathway for their catabolism, but the molecular mechanisms mediating LD degradation by selective autophagy (lipophagy) are unknown. Here we identify spartin as a receptor localizing to LDs and interacting with core autophagy machinery, and we show that spartin is required to deliver LDs to lysosomes for triglyceride mobilization. Mutations in SPART (encoding spartin) lead to Troyer syndrome, a form of complex hereditary spastic paraplegia1. Interfering with spartin function in cultured human neurons or murine brain neurons leads to LD and triglyceride accumulation. Our identification of spartin as a lipophagy receptor, thus, suggests that impaired LD turnover contributes to Troyer syndrome development.

Metabolic regulation of CTCF expression and chromatin association dictates starvation response in mice and flies

Coordinated temporal control of gene expression is essential for physiological homeostasis, especially during metabolic transitions. However, the interplay between chromatin architectural proteins and metabolism in regulating transcription is less understood. Here, we demonstrate a conserved bidirectional interplay between CTCF (CCCTC-binding factor) expression/function and metabolic inputs during feed-fast cycles. Our results indicate that its loci-specific functional diversity is associated with physiological plasticity in mouse hepatocytes. CTCF differential expression and long non-coding RNA-Jpx mediated changes in chromatin occupancy, unraveled its paradoxical yet tuneable functions, which are governed by metabolic inputs. We illustrate the key role of CTCF in controlling temporal cascade of transcriptional response, with effects on hepatic mitochondrial energetics and lipidome. Underscoring the evolutionary conservation of CTCF-dependent metabolic homeostasis, CTCF knockdown in flies abrogated starvation resistance. In summary, we demonstrate the interplay between CTCF and metabolic inputs that highlights the coupled plasticity of physiological responses and chromatin function.

The NuRD complex cooperates with SALL4 to orchestrate reprogramming

Cell fate decision involves rewiring of the genome, but remains poorly understood at the chromatin level. Here, we report that chromatin remodeling complex NuRD participates in closing open chromatin in the early phase of somatic reprogramming. Sall4, Jdp2, Glis1 and Esrrb can reprogram MEFs to iPSCs efficiently, but only Sall4 is indispensable capable of recruiting endogenous components of NuRD. Yet knocking down NuRD components only reduces reprogramming modestly, in contrast to disrupting the known Sall4-NuRD interaction by mutating or deleting the NuRD interacting motif at its N-terminus that renders Sall4 inept to reprogram. Remarkably, these defects can be partially rescured by grafting NuRD interacting motif onto Jdp2. Further analysis of chromatin accessibility dynamics demonstrates that the Sall4-NuRD axis plays a critical role in closing the open chromatin in the early phase of reprogramming. Among the chromatin loci closed by Sall4-NuRD encode genes resistant to reprogramming. These results identify a previously unrecognized role of NuRD in reprogramming, and may further illuminate chromatin closing as a critical step in cell fate control.

SGRN: A Cas12a-driven Synthetic Gene Regulatory Network System

Gene regulatory networks, which control gene expression patterns in development and in response to stimuli, use regulatory logic modules to coordinate inputs and outputs. One example of a regulatory logic module is the gene regulatory cascade (GRC), where a series of transcription factor genes turn on in order. Synthetic biologists have derived artificial systems that encode regulatory rules, including GRCs. Furthermore, the development of single-cell approaches has enabled the discovery of gene regulatory modules in a variety of experimental settings. However, the tools available for validating these observations remain limited. Based on a synthetic GRC using DNA cutting-defective Cas9 (dCas9), we designed and implemented an alternative synthetic GRC utilizing DNA cutting-defective Cas12a (dCas12a). Comparing the ability of these two systems to express a fluorescent reporter, the dCas9 system was initially more active, while the dCas12a system was more streamlined. Investigating the influence of individual components of the systems identified nuclear localization as a major driver of differences in activity. Improving nuclear localization for the dCas12a system resulted in 1.5-fold more reporter-positive cells and a 15-fold increase in reporter intensity relative to the dCas9 system. We call this optimized system the "Synthetic Gene Regulatory Network" (SGRN, pronounced "sojourn").

Regulation of Cell Plasticity by Bromodomain and Extraterminal Domain (BET) Proteins: A New Perspective in Glioblastoma Therapy

BET proteins are a family of multifunctional epigenetic readers, mainly involved in transcriptional regulation through chromatin modelling. Transcriptome handling ability of BET proteins suggests a key role in the modulation of cell plasticity, both in fate decision and in lineage commitment during embryonic development and in pathogenic conditions, including cancerogenesis. Glioblastoma is the most aggressive form of glioma, characterized by a very poor prognosis despite the application of a multimodal therapy. Recently, new insights are emerging about the glioblastoma cellular origin, leading to the hypothesis that several putative mechanisms occur during gliomagenesis. Interestingly, epigenome dysregulation associated with loss of cellular identity and functions are emerging as crucial features of glioblastoma pathogenesis. Therefore, the emerging roles of BET protein in glioblastoma onco-biology and the compelling demand for more effective therapeutic strategies suggest that BET family members could be promising targets for translational breakthroughs in glioblastoma treatment. Primarily, “Reprogramming Therapy”, which is aimed at reverting the malignant phenotype, is now considered a promising strategy for GBM therapy.

Molecular Dissection of Somatic Skeletal Disease in Neurofibromatosis Type 1

Neurofibromatosis type 1 (NF1) is a tumor predisposition syndrome caused by heterozygous NF1 gene mutations. Patients with NF1 present with pleiotropic somatic secondary manifestations, including development of bone pseudarthrosis after fracture. Somatic NF1 gene mutations were reproducibly identified in patient-derived pseudarthrosis specimens, suggesting a local mosaic cell population including somatic pathologic cells. The somatic cellular pathogenesis of NF1 pseudarthroses remains unclear, though defects in osteogenesis have been posited. Here, we applied time-series single-cell RNA-sequencing (scRNA-seq) to patient-matched control and pseudarthrosis-derived primary bone stromal cells (BSCs). We show that osteogenic specification to an osteoblast progenitor cell population was evident for control bone-derived cells and haploinsufficient pseudarthrosis-derived cells. Similar results were observed for somatic patient fracture-derived NF1^-/- cells; however, expression of genetic pathways associated with skeletal mineralization were significantly reduced in NF1^-/- cells compared with fracture-derived NF1^+/- cells. In mice, we show that Nf1 expressed in bone marrow osteoprogenitors is required for the maintenance of the adult skeleton. Results from our study implicate impaired Clec11a-Itga11-Wnt signaling in the pathogenesis of NF1-associated skeletal disease. © 2022 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Efficient generation of lower induced motor neurons by coupling Ngn2 expression with developmental cues

Human pluripotent stem cells (hPSCs) are a powerful tool for disease modeling of hard-to-access tissues (such as the brain). Current protocols either direct neuronal differentiation with small molecules or use transcription-factor-mediated programming. In this study, we couple overexpression of transcription factor Neurogenin2 (Ngn2) with small molecule patterning to differentiate hPSCs into lower induced motor neurons (liMoNes/liMNs). This approach induces canonical MN markers including MN-specific Hb9/MNX1 in more than 95% of cells. liMNs resemble bona fide hPSC-derived MN, exhibit spontaneous electrical activity, express synaptic markers, and can contact muscle cells in vitro. Pooled, multiplexed single-cell RNA sequencing on 50 hPSC lines reveals reproducible populations of distinct subtypes of cervical and brachial MNs that resemble their in vivo, embryonic counterparts. Combining small molecule patterning with Ngn2 overexpression facilitates high-yield, reproducible production of disease-relevant MN subtypes, which is fundamental in propelling our knowledge of MN biology and its disruption in disease.

Single-cell transcriptomic analysis reveals diversity within mammalian spinal motor neurons

Spinal motor neurons (MNs) integrate sensory stimuli and brain commands to generate movements. In vertebrates, the molecular identities of the cardinal MN types such as those innervating limb versus trunk muscles are well elucidated. Yet the identities of finer subtypes within these cell populations that innervate individual muscle groups remain enigmatic. Here we investigate heterogeneity in mouse MNs using single-cell transcriptomics. Among limb-innervating MNs, we reveal a diverse neuropeptide code for delineating putative motor pool identities. Additionally, we uncover that axial MNs are subdivided into three molecularly distinct subtypes, defined by mediolaterally-biased Satb2, Nr2f2 or Bcl11b expression patterns with different axon guidance signatures. These three subtypes are present in chicken and human embryos, suggesting a conserved axial MN expression pattern across higher vertebrates. Overall, our study provides a molecular resource of spinal MN types and paves the way towards deciphering how neuronal subtypes evolved to accommodate vertebrate motor behaviors.

Profiling Dynamic Changes in DNA Accessibility During Axon Regeneration After Optic Nerve Crush in Adult Zebrafish

A time-course series utilizing assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) can be used to detect changes in accessibility of DNA regulatory elements such as promoters and enhancers over the course of regeneration. This chapter describes methods for preparing ATAC-seq libraries from isolated zebrafish retinal ganglion cells (RGCs) following optic nerve crush at selected post-injury time points. These methods have been used for identifying dynamic changes in DNA accessibility that govern successful optic nerve regeneration in zebrafish. This method may be adapted to identify changes in DNA accessibility that accompany other types of insults to RGCs or to identify changes that occur over the course of development.

Multiomic Analysis of Neurons with Divergent Projection Patterns Identifies Novel Regulators of Axon Pathfinding

Axon pathfinding is a key step in neural circuits formation. However, the transcriptional mechanisms regulating its progression remain poorly understood. The binary decision of crossing or avoiding the midline taken by some neuronal axons during development represents a robust model to investigate the mechanisms that control the selection of axonal trajectories. Here, to identify novel regulators of axon guidance, this work compares the transcriptome and chromatin occupancy profiles of two neuronal subpopulations, ipsilateral (iRGC) and contralateral retinal ganglion cells (cRGC), with similar functions but divergent axon trajectories. These analyses retrieved a number of genes encoding for proteins not previously implicated in axon pathfinding. In vivo functional experiments confirm the implication of some of these candidates in axonal navigation. Among the candidate genes, γ-synuclein is identified as essential for inducing midline crossing. Footprint and luciferase assays demonstrate that this small-sized protein is regulated by the transcription factor (TF) Pou4f1 in cRGCs. It is also shown that Lhx2/9 are specifically expressed in iRGCs and control a program that partially overlaps with that regulated by Zic2, previously described as essential for iRGC specification. Overall, the analyses identify dozens of new molecules potentially involved in axon guidance and reveal the regulatory logic behind the selection of axonal trajectories.

Transcriptional dynamics of murine motor neuron maturation in vivo and in vitro

Neurons born in the embryo can undergo a protracted period of maturation lasting well into postnatal life. How gene expression changes are regulated during maturation and whether they can be recapitulated in cultured neurons remains poorly understood. Here, we show that mouse motor neurons exhibit pervasive changes in gene expression and accessibility of associated regulatory regions from embryonic till juvenile age. While motifs of selector transcription factors, ISL1 and LHX3, are enriched in nascent regulatory regions, motifs of NFI factors, activity-dependent factors, and hormone receptors become more prominent in maturation-dependent enhancers. Notably, stem cell-derived motor neurons recapitulate ~40% of the maturation expression program in vitro, with neural activity playing only a modest role as a late-stage modulator. Thus, the genetic maturation program consists of a core hardwired subprogram that is correctly executed in vitro and an extrinsically-controlled subprogram that is dependent on the in vivo context of the maturing organism.

DiffChIPL: a differential peak analysis method for high-throughput sequencing data with biological replicates based on limma

MOTIVATION

ChIP-seq detects protein-DNA interactions within chromatin, such as that of chromatin structural components and transcription machinery. ChIP-seq profiles are often noisy and variable across replicates, posing a challenge to the development of effective algorithms to accurately detect differential peaks. Methods have recently been designed for this purpose but sometimes yield conflicting results that are inconsistent with the underlying biology. Most existing algorithms perform well on limited datasets. To improve differential analysis of ChIP-seq, we present a novel Differential analysis method for ChIP-seq based on Limma (DiffChIPL).

RESULTS

DiffChIPL is adaptive to asymmetrical or symmetrical data and can accurately report global differences. We used simulated and real datasets for transcription factors (TFs) and histone modification marks to validate and benchmark our algorithm. DiffChIPL shows superior performance in sensitivity and false positive rate in different simulations and control datasets. DiffChIPL also performs well on real ChIP-seq, CUT&RUN, CUT&Tag and ATAC-seq datasets. DiffChIPL is an accurate and robust method, exhibiting better performance in differential analysis for a variety of applications including TF binding, histone modifications and chromatin accessibility.

AVAILABILITY AND IMPLEMENTATION

https://github.com/yancychy/DiffChIPL.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Foxa2 and Pet1 Direct and Indirect Synergy Drive Serotonergic Neuronal Differentiation

Neuronal programming by forced expression of transcription factors (TFs) holds promise for clinical applications of regenerative medicine. However, the mechanisms by which TFs coordinate their activities on the genome and control distinct neuronal fates remain obscure. Using direct neuronal programming of embryonic stem cells, we dissected the contribution of a series of TFs to specific neuronal regulatory programs. We deconstructed the Ascl1-Lmx1b-Foxa2-Pet1 TF combination that has been shown to generate serotonergic neurons and found that stepwise addition of TFs to Ascl1 canalizes the neuronal fate into a diffuse monoaminergic fate. The addition of pioneer factor Foxa2 represses Phox2b to induce serotonergic fate, similar to in vivo regulatory networks. Foxa2 and Pet1 appear to act synergistically to upregulate serotonergic fate. Foxa2 and Pet1 co-bind to a small fraction of genomic regions but mostly bind to different regulatory sites. In contrast to the combinatorial binding activities of other programming TFs, Pet1 does not strictly follow the Foxa2 pioneer. These findings highlight the challenges in formulating generalizable rules for describing the behavior of TF combinations that program distinct neuronal subtypes.

Capybara: A computational tool to measure cell identity and fate transitions

Measuring cell identity in development, disease, and reprogramming is challenging as cell types and states are in continual transition. Here, we present Capybara, a computational tool to classify discrete cell identity and intermediate "hybrid" cell states, supporting a metric to quantify cell fate transition dynamics. We validate hybrid cells using experimental lineage tracing data to demonstrate the multi-lineage potential of these intermediate cell states. We apply Capybara to diagnose shortcomings in several cell engineering protocols, identifying hybrid states in cardiac reprogramming and off-target identities in motor neuron programming, which we alleviate by adding exogenous signaling factors. Further, we establish a putative in vivo correlate for induced endoderm progenitors. Together, these results showcase the utility of Capybara to dissect cell identity and fate transitions, prioritizing interventions to enhance the efficiency and fidelity of stem cell engineering.

The BTB transcription factors ZBTB11 and ZFP131 maintain pluripotency by repressing pro-differentiation genes

In pluripotent cells, a delicate activation-repression balance maintains pro-differentiation genes ready for rapid activation. The identity of transcription factors (TFs) that specifically repress pro-differentiation genes remains obscure. By targeting ∼1,700 TFs with CRISPR loss-of-function screen, we found that ZBTB11 and ZFP131 are required for embryonic stem cell (ESC) pluripotency. ESCs without ZBTB11 or ZFP131 lose colony morphology, reduce proliferation rate, and upregulate transcription of genes associated with three germ layers. ZBTB11 and ZFP131 bind proximally to pro-differentiation genes. ZBTB11 or ZFP131 loss leads to an increase in H3K4me3, negative elongation factor (NELF) complex release, and concomitant transcription at associated genes. Together, our results suggest that ZBTB11 and ZFP131 maintain pluripotency by preventing premature expression of pro-differentiation genes and present a generalizable framework to maintain cellular potency.

The histone demethylase Kdm6b regulates subtype diversification of mouse spinal motor neurons during development

How a single neuronal population diversifies into subtypes with distinct synaptic targets is a fundamental topic in neuroscience whose underlying mechanisms are unclear. Here, we show that the histone H3-lysine 27 demethylase Kdm6b regulates the diversification of motor neurons to distinct subtypes innervating different muscle targets during spinal cord development. In mouse embryonic motor neurons, Kdm6b promotes the medial motor column (MMC) and hypaxial motor column (HMC) fates while inhibiting the lateral motor column (LMC) and preganglionic motor column (PGC) identities. Our single-cell RNA-sequencing analyses reveal the heterogeneity of PGC, LMC, and MMC motor neurons. Further, our single-cell RNA-sequencing data, combined with mouse model studies, demonstrates that Kdm6b acquires cell fate specificity together with the transcription factor complex Isl1-Lhx3. Our study provides mechanistic insight into the gene regulatory network regulating neuronal cell-type diversification and defines a regulatory role of Kdm6b in the generation of motor neuron subtypes in the mouse spinal cord.

Establishing the Molecular and Functional Diversity of Spinal Motoneurons

Spinal motoneurons are a remarkably diverse class of neurons responsible for facilitating a broad range of motor behaviors and autonomic functions. Studies of motoneuron differentiation have provided fundamental insights into the developmental mechanisms of neuronal diversification, and have illuminated principles of neural fate specification that operate throughout the central nervous system. Because of their relative anatomical simplicity and accessibility, motoneurons have provided a tractable model system to address multiple facets of neural development, including early patterning, neuronal migration, axon guidance, and synaptic specificity. Beyond their roles in providing direct communication between central circuits and muscle, recent studies have revealed that motoneuron subtype-specific programs also play important roles in determining the central connectivity and function of motor circuits. Cross-species comparative analyses have provided novel insights into how evolutionary changes in subtype specification programs may have contributed to adaptive changes in locomotor behaviors. This chapter focusses on the gene regulatory networks governing spinal motoneuron specification, and how studies of spinal motoneurons have informed our understanding of the basic mechanisms of neuronal specification and spinal circuit assembly.

Making neurons, made easy: The use of Neurogenin-2 in neuronal differentiation

Directed neuronal differentiation of human pluripotent stem cells (hPSCs), neural progenitors, or fibroblasts using transcription factors has allowed for the rapid and highly reproducible differentiation of mature and functional neurons. Exogenous expression of the transcription factor Neurogenin-2 (NGN2) has been widely used to generate different populations of neurons, which have been used in neurodevelopment studies, disease modeling, drug screening, and neuronal replacement therapies. Could NGN2 be a “one-glove-fits-all” approach for neuronal differentiations? This review summarizes the cellular roles of NGN2 and describes the applications and limitations of using NGN2 for the rapid and directed differentiation of neurons.

Transcriptional regulation of neuronal identity

Neuronal diversity is an intrinsic feature of the nervous system. Transcription factors (TFs) are key regulators in the establishment of different neuronal identities; how are the actions of different TFs coordinated to orchestrate this diversity? Are there common features shared among the different neuron types of an organism or even among different animal groups? In this review, we provide a brief overview on common traits emerging on the transcriptional regulation of neuron type diversification with a special focus on the comparison between mouse and Caenorhabditis elegans model systems. In the first part, we describe general concepts on neuronal identity and transcriptional regulation of gene expression. In the second part of the review, TFs are classified in different categories according to their key roles at specific steps along the protracted process of neuronal specification and differentiation. The same TF categories can be identified both in mammals and nematodes. Importantly, TFs are very pleiotropic: Depending on the neuron type or the time in development, the same TF can fulfil functions belonging to different categories. Finally, we describe the key role of transcriptional repression at all steps controlling neuronal diversity and propose that acquisition of neuronal identities could be considered a metastable process.

A shared transcriptional code orchestrates temporal patterning of the central nervous system

The molecular mechanisms that produce the full array of neuronal subtypes in the vertebrate nervous system are incompletely understood. Here, we provide evidence of a global temporal patterning program comprising sets of transcription factors that stratifies neurons based on the developmental time at which they are generated. This transcriptional code acts throughout the central nervous system, in parallel to spatial patterning, thereby increasing the diversity of neurons generated along the neuraxis. We further demonstrate that this temporal program operates in stem cell-derived neurons and is under the control of the TGFβ signaling pathway. Targeted perturbation of components of the temporal program, Nfia and Nfib, reveals their functional requirement for the generation of late-born neuronal subtypes. Together, our results provide evidence for the existence of a previously unappreciated global temporal transcriptional program of neuronal subtype identity and suggest that the integration of spatial and temporal patterning mechanisms diversifies and organizes neuronal subtypes in the vertebrate nervous system.

Discovering differential genome sequence activity with interpretable and efficient deep learning

Discovering sequence features that differentially direct cells to alternate fates is key to understanding both cellular development and the consequences of disease related mutations. We introduce Expected Pattern Effect and Differential Expected Pattern Effect, two black-box methods that can interpret genome regulatory sequences for cell type-specific or condition specific patterns. We show that these methods identify relevant transcription factor motifs and spacings that are predictive of cell state-specific chromatin accessibility. Finally, we integrate these methods into framework that is readily accessible to non-experts and available for download as a binary or installed via PyPI or bioconda at https://cgs.csail.mit.edu/deepaccess-package/.

Direct Neuronal Reprogramming: Bridging the Gap Between Basic Science and Clinical Application

Direct neuronal reprogramming is an innovative new technology that involves the conversion of somatic cells to induced neurons (iNs) without passing through a pluripotent state. The capacity to make new neurons in the brain, which previously was not achievable, has created great excitement in the field as it has opened the door for the potential treatment of incurable neurodegenerative diseases and brain injuries such as stroke. These neurological disorders are associated with frank neuronal loss, and as new neurons are not made in most of the adult brain, treatment options are limited. Developmental biologists have paved the way for the field of direct neuronal reprogramming by identifying both intrinsic cues, primarily transcription factors (TFs) and miRNAs, and extrinsic cues, including growth factors and other signaling molecules, that induce neurogenesis and specify neuronal subtype identities in the embryonic brain. The striking observation that postmitotic, terminally differentiated somatic cells can be converted to iNs by mis-expression of TFs or miRNAs involved in neural lineage development, and/or by exposure to growth factors or small molecule cocktails that recapitulate the signaling environment of the developing brain, has opened the door to the rapid expansion of new neuronal reprogramming methodologies. Furthermore, the more recent applications of neuronal lineage conversion strategies that target resident glial cells in situ has expanded the clinical potential of direct neuronal reprogramming techniques. Herein, we present an overview of the history, accomplishments, and therapeutic potential of direct neuronal reprogramming as revealed over the last two decades.

Urine-derived induced pluripotent/neural stem cells for modeling neurological diseases

Neurological diseases are mainly modeled using rodents through gene editing, surgery or injury approaches. However, differences between humans and rodents in terms of genetics, neural development, and physiology pose limitations on studying disease pathogenesis in rodent models for neuroscience research. In the past decade, the generation of induced pluripotent stem cells (iPSCs) and induced neural stem cells (iNSCs) by reprogramming somatic cells offers a powerful alternative for modeling neurological diseases and for testing regenerative medicines. Among the different somatic cell types, urine-derived stem cells (USCs) are an ideal cell source for iPSC and iNSC reprogramming, as USCs are highly proliferative, multipotent, epithelial in nature, and easier to reprogram than skin fibroblasts. In addition, the use of USCs represents a simple, low-cost and non-invasive procedure for generating iPSCs/iNSCs. This review describes the cellular and molecular properties of USCs, their differentiation potency, different reprogramming methods for the generation of iPSCs/iNSCs, and their potential applications in modeling neurological diseases.

Unraveling Root Development Through Single-Cell Omics and Reconstruction of Gene Regulatory Networks

Over the last decades, research on postembryonic root development has been facilitated by "omics" technologies. Among these technologies, microarrays first, and RNA sequencing (RNA-seq) later, have provided transcriptional information on the underlying molecular processes establishing the basis of System Biology studies in roots. Cell fate specification and development have been widely studied in the primary root, which involved the identification of many cell type transcriptomes and the reconstruction of gene regulatory networks (GRN). The study of lateral root (LR) development has not been an exception. However, the molecular mechanisms regulating cell fate specification during LR formation remain largely unexplored. Recently, single-cell RNA-seq (scRNA-seq) studies have addressed the specification of tissues from stem cells in the primary root. scRNA-seq studies are anticipated to be a useful approach to decipher cell fate specification and patterning during LR formation. In this review, we address the different scRNA-seq strategies used both in plants and animals and how we could take advantage of scRNA-seq to unravel new regulatory mechanisms and reconstruct GRN. In addition, we discuss how to integrate scRNA-seq results with previous RNA-seq datasets and GRN. We also address relevant findings obtained through single-cell based studies and how LR developmental studies could be facilitated by scRNA-seq approaches and subsequent GRN inference. The use of single-cell approaches to investigate LR formation could help to decipher fundamental biological mechanisms such as cell memory, synchronization, polarization, or pluripotency.

Onecut Regulates Core Components of the Molecular Machinery for Neurotransmission in Photoreceptor Differentiation

Photoreceptor cells (PRC) are neurons highly specialized for sensing light stimuli and have considerably diversified during evolution. The genetic mechanisms that underlie photoreceptor differentiation and accompanied the progressive increase in complexity and diversification of this sensory cell type are a matter of great interest in the field. A role of the homeodomain transcription factor Onecut (Oc) in photoreceptor cell formation is proposed throughout multicellular organisms. However, knowledge of the identity of the Oc downstream-acting factors that mediate specific tasks in the differentiation of the PRC remains limited. Here, we used transgenic perturbation of the Ciona robusta Oc protein to show its requirement for ciliary PRC differentiation. Then, transcriptome profiling between the trans-activation and trans-repression Oc phenotypes identified differentially expressed genes that are enriched in exocytosis, calcium homeostasis, and neurotransmission. Finally, comparison of RNA-Seq datasets in Ciona and mouse identifies a set of Oc downstream genes conserved between tunicates and vertebrates. The transcription factor Oc emerges as a key regulator of neurotransmission in retinal cell types.

Transcriptional Programming of Human Mechanosensory Neuron Subtypes from Pluripotent Stem Cells

Efficient and homogeneous in vitro generation of peripheral sensory neurons may provide a framework for novel drug screening platforms and disease models of touch and pain. We discover that, by ovesssrexpressing NGN2 and BRN3A, human pluripotent stem cells can be transcriptionally programmed to differentiate into a surprisingly uniform culture of cold- and mechano-sensing neurons. Although such a neuronal subtype is not found in mice, we identify molecular evidence for its existence in human sensory ganglia. Combining NGN2 and BRN3A programming with neural crest patterning, we produce two additional populations of sensory neurons, including a specialized touch receptor neuron subtype. Finally, we apply this system to model a rare inherited sensory disorder of touch and proprioception caused by inactivating mutations in PIEZO2. Together, these findings establish an approach to specify distinct sensory neuron subtypes in vitro, underscoring the utility of stem cell technology to capture human-specific features of physiology and disease.

Nickolls et al. develop a method, using human stem cells, to generate specific types of sensory neurons that detect cold temperature and mechanical force. This approach uncovers a class of neuron found in humans, but not mice, and enables the modeling of a rare sensory disorder of touch and proprioception.

A cis-regulatory change underlying the motor neuron-specific loss of Ebf expression in immotile tunicate larvae

Many species in the tunicate family Molgulidae have independently lost their swimming larval form and instead develop as tailless, immotile larvae. These larvae do not develop structures that are essential for swimming such as the notochord, otolith, and tail muscles. However, little is known about neural development in these nonswimming larvae. Here, we studied the patterning of the Motor Ganglion (MG) of Molgula occulta, a nonswimming species. We found that spatial patterns of MG neuron regulators in this species are conserved, compared with species with swimming larvae, suggesting that the gene networks regulating their expression are intact despite the loss of swimming. However, expression of the key motor neuron regulatory gene Ebf (Collier/Olf/EBF) was reduced in the developing MG of M. occulta when compared with molgulid species with swimming larvae. This was corroborated by measuring allele-specific expression of Ebf in hybrid embryos from crosses of M. occulta with the swimming species M. oculata. Heterologous reporter construct assays in the model tunicate species Ciona robusta revealed a specific cis-regulatory sequence change that reduces expression of Ebf in the MG, but not in other cells. Taken together, these data suggest that MG neurons are still specified in M. occulta larvae, but their differentiation might be impaired due to reduction of Ebf expression levels.

Inferring time series chromatin states for promoter-enhancer pairs based on Hi-C data

BACKGROUND

Co-localized combinations of histone modifications ("chromatin states") have been shown to correlate with promoter and enhancer activity. Changes in chromatin states over multiple time points ("chromatin state trajectories") have previously been analyzed at promoter and enhancers separately. With the advent of time series Hi-C data it is now possible to connect promoters and enhancers and to analyze chromatin state trajectories at promoter-enhancer pairs.

RESULTS

We present TimelessFlex, a framework for investigating chromatin state trajectories at promoters and enhancers and at promoter-enhancer pairs based on Hi-C information. TimelessFlex extends our previous approach Timeless, a Bayesian network for clustering multiple histone modification data sets at promoter and enhancer feature regions. We utilize time series ATAC-seq data measuring open chromatin to define promoters and enhancer candidates. We developed an expectation-maximization algorithm to assign promoters and enhancers to each other based on Hi-C interactions and jointly cluster their feature regions into paired chromatin state trajectories. We find jointly clustered promoter-enhancer pairs showing the same activation patterns on both sides but with a stronger trend at the enhancer side. While the promoter side remains accessible across the time series, the enhancer side becomes dynamically more open towards the gene activation time point. Promoter cluster patterns show strong correlations with gene expression signals, whereas Hi-C signals get only slightly stronger towards activation. The code of the framework is available at https://github.com/henriettemiko/TimelessFlex .

CONCLUSIONS

TimelessFlex clusters time series histone modifications at promoter-enhancer pairs based on Hi-C and it can identify distinct chromatin states at promoter and enhancer feature regions and their changes over time.

An interpretable bimodal neural network characterizes the sequence and preexisting chromatin predictors of induced transcription factor binding

BACKGROUND

Transcription factor (TF) binding specificity is determined via a complex interplay between the transcription factor's DNA binding preference and cell type-specific chromatin environments. The chromatin features that correlate with transcription factor binding in a given cell type have been well characterized. For instance, the binding sites for a majority of transcription factors display concurrent chromatin accessibility. However, concurrent chromatin features reflect the binding activities of the transcription factor itself and thus provide limited insight into how genome-wide TF-DNA binding patterns became established in the first place. To understand the determinants of transcription factor binding specificity, we therefore need to examine how newly activated transcription factors interact with sequence and preexisting chromatin landscapes.

RESULTS

Here, we investigate the sequence and preexisting chromatin predictors of TF-DNA binding by examining the genome-wide occupancy of transcription factors that have been induced in well-characterized chromatin environments. We develop Bichrom, a bimodal neural network that jointly models sequence and preexisting chromatin data to interpret the genome-wide binding patterns of induced transcription factors. We find that the preexisting chromatin landscape is a differential global predictor of TF-DNA binding; incorporating preexisting chromatin features improves our ability to explain the binding specificity of some transcription factors substantially, but not others. Furthermore, by analyzing site-level predictors, we show that transcription factor binding in previously inaccessible chromatin tends to correspond to the presence of more favorable cognate DNA sequences.

CONCLUSIONS

Bichrom thus provides a framework for modeling, interpreting, and visualizing the joint sequence and chromatin landscapes that determine TF-DNA binding dynamics.

Neuronal differentiation strategies: insights from single-cell sequencing and machine learning

Neuronal replacement therapies rely on the in vitro differentiation of specific cell types from embryonic or induced pluripotent stem cells, or on the direct reprogramming of differentiated adult cells via the expression of transcription factors or signaling molecules. The factors used to induce differentiation or reprogramming are often identified by informed guesses based on differential gene expression or known roles for these factors during development. Moreover, differentiation protocols usually result in partly differentiated cells or the production of a mix of cell types. In this Hypothesis article, we suggest that, to overcome these inefficiencies and improve neuronal differentiation protocols, we need to take into account the developmental history of the desired cell types. Specifically, we present a strategy that uses single-cell sequencing techniques combined with machine learning as a principled method to select a sequence of programming factors that are important not only in adult neurons but also during differentiation.

Synthetic Developmental Biology: Understanding Through Reconstitution

Reconstitution is an experimental strategy that seeks to recapitulate biological events outside their natural contexts using a reduced set of components. Classically, biochemical reconstitution has been extensively applied to identify the minimal set of molecules sufficient for recreating the basic chemistry of life. By analogy, reconstitution approaches to developmental biology recapitulate aspects of developmental events outside an embryo, with the goal of revealing the basic genetic circuits or physical cues sufficient for recreating developmental decisions. The rapidly growing repertoire of genetic, molecular, microscopic, and bioengineering tools is expanding the complexity and precision of reconstitution experiments. We review the emerging field of synthetic developmental biology, with a focus on the ways in which reconstitution strategies and new biological tools have enhanced our modern understanding of fundamental questions in developmental biology.

Vsx1 and Chx10 paralogs sequentially secure V2 interneuron identity during spinal cord development

Paralog factors are usually described as consolidating biological systems by displaying redundant functionality in the same cells. Here, we report that paralogs can also cooperate in distinct cell populations at successive stages of differentiation. In mouse embryonic spinal cord, motor neurons and V2 interneurons differentiate from adjacent progenitor domains that share identical developmental determinants. Therefore, additional strategies secure respective cell fate. In particular, Hb9 promotes motor neuron identity while inhibiting V2 differentiation, whereas Chx10 stimulates V2a differentiation while repressing motor neuron fate. However, Chx10 is not present at the onset of V2 differentiation and in other V2 populations. In the present study, we show that Vsx1, the single paralog of Chx10, which is produced earlier than Chx10 in V2 precursors, can inhibit motor neuron differentiation and promote V2 interneuron production. However, the single absence of Vsx1 does not impact on V2 fate consolidation, suggesting that lack of Vsx1 may be compensated by other factors. Nevertheless, Vsx1 cooperates with Chx10 to prevent motor neuron differentiation in early V2 precursors although these two paralog factors are not produced in the same cells. Hence, this study uncovers an original situation, namely labor division, wherein paralog genes cooperate at successive steps of neuronal development.

Identification of determinants of differential chromatin accessibility through a massively parallel genome-integrated reporter assay

A key mechanism in cellular regulation is the ability of the transcriptional machinery to physically access DNA. Transcription factors interact with DNA to alter the accessibility of chromatin, which enables changes to gene expression during development or disease or as a response to environmental stimuli. However, the regulation of DNA accessibility via the recruitment of transcription factors is difficult to study in the context of the native genome because every genomic site is distinct in multiple ways. Here we introduce the multiplexed integrated accessibility assay (MIAA), an assay that measures chromatin accessibility of synthetic oligonucleotide sequence libraries integrated into a controlled genomic context with low native accessibility. We apply MIAA to measure the effects of sequence motifs on cell type-specific accessibility between mouse embryonic stem cells and embryonic stem cell-derived definitive endoderm cells, screening 7905 distinct DNA sequences. MIAA recapitulates differential accessibility patterns of 100-nt sequences derived from natively differential genomic regions, identifying E-box motifs common to epithelial-mesenchymal transition driver transcription factors in stem cell-specific accessible regions that become repressed in endoderm. We show that a single binding motif for a key regulatory transcription factor is sufficient to open chromatin, and classify sets of stem cell-specific, endoderm-specific, and shared accessibility-modifying transcription factor motifs. We also show that overexpression of two definitive endoderm transcription factors, T and Foxa2, results in changes to accessibility in DNA sequences containing their respective DNA-binding motifs and identify preferential motif arrangements that influence accessibility.

An Amyotrophic Lateral Sclerosis-Associated Mutant of C21ORF2 Is Stabilized by NEK1-Mediated Hyperphosphorylation and the Inability to Bind FBXO3

C21ORF2 and NEK1 have been identified as amyotrophic lateral sclerosis (ALS)-associated genes. Both genes are also mutated in certain ciliopathies, suggesting that they might contribute to the same signaling pathways. Here we show that FBXO3, the substrate receptor of an SCF ubiquitin ligase complex, binds and ubiquitylates C21ORF2, thereby targeting it for proteasomal degradation. C21ORF2 stabilizes the kinase NEK1, with the result that loss of FBXO3 stabilizes not only C21ORF2 but also NEK1. Conversely, NEK1-mediated phosphorylation stabilizes C21ORF2 by attenuating its interaction with FBXO3. We found that the ALS-associated V58L mutant of C21ORF2 is more susceptible to phosphorylation by NEK1, with the result that it is not ubiquitylated by FBXO3 and therefore accumulates together with NEK1. Expression of C21ORF2(V58L) in motor neurons induced from mouse embryonic stem cells impaired neurite outgrowth. We suggest that inhibition of NEK1 activity is a potential therapeutic approach to ALS associated with C21ORF2 mutation.

Simple Method to Quantify Protein Abundances from 1000 Cells

The rise of single-cell transcriptomics has created an urgent need for similar approaches that use a minimal number of cells to quantify expression levels of proteins. We integrated and optimized multiple recent developments to establish a proteomics workflow to quantify proteins from as few as 1000 mammalian stem cells. The method uses chemical peptide labeling, does not require specific equipment other than cell lysis tools, and quantifies >2500 proteins with high reproducibility. We validated the method by comparing mouse embryonic stem cells and in vitro differentiated motor neurons. We identify differentially expressed proteins with small fold changes and a dynamic range in abundance similar to that of standard methods. Protein abundance measurements obtained with our protocol compared well to corresponding transcript abundance and to measurements using standard inputs. The protocol is also applicable to other systems, such as fluorescence-activated cell sorting (FACS)-purified cells from the tunicate Ciona. Therefore, we offer a straightforward and accurate method to acquire proteomics data from minimal input samples.

Sequence and chromatin determinants of transcription factor binding and the establishment of cell type-specific binding patterns

Transcription factors (TFs) selectively bind distinct sets of sites in different cell types. Such cell type-specific binding specificity is expected to result from interplay between the TF's intrinsic sequence preferences, cooperative interactions with other regulatory proteins, and cell type-specific chromatin landscapes. Cell type-specific TF binding events are highly correlated with patterns of chromatin accessibility and active histone modifications in the same cell type. However, since concurrent chromatin may itself be a consequence of TF binding, chromatin landscapes measured prior to TF activation provide more useful insights into how cell type-specific TF binding events became established in the first place. Here, we review the various sequence and chromatin determinants of cell type-specific TF binding specificity. We identify the current challenges and opportunities associated with computational approaches to characterizing, imputing, and predicting cell type-specific TF binding patterns. We further focus on studies that characterize TF binding in dynamic regulatory settings, and we discuss how these studies are leading to a more complex and nuanced understanding of dynamic protein-DNA binding activities. We propose that TF binding activities at individual sites can be viewed along a two-dimensional continuum of local sequence and chromatin context. Under this view, cell type-specific TF binding activities may result from either strongly favorable sequence features or strongly favorable chromatin context.

Olig2 positive cells derived from hair follicle neural crest stem cells in rats

Motor neuron disease (MND) is a kind of common clinical nervous system disease with typical characteristic of progressive motor neurons degeneration or death. Motor neuron derived from stem cells or motor neuron progenitor cells will be a good choice to be used for treatment of the disease. In this study, we used the combination of 5 small molecular including CHIR99021 (CHIR), SB431542 (SB), DMH1 (DMH), retinoic acid (RA) and Purmorphamine (Pur) to induce hair follicles neural crest stem cells (hfNCSCs) to motor neurons progenitors (MNPs). Valproic acid (VPA) was used to make MNPs proliferation. RA and Pur were used to try to induce MNPs toward motor neurons (MNs) and CpdE was tried for MNs maturation. Nestin, β-tubulin Ш (Tuj1), microtubule associated protein 2 (MAP2), Olig2, choline acetyltransferase (ChAT)and TUBB3 were examined at protein and mRNA levels by immunofluoresence cytochemistry, western blot and real time PCR at 6, 16 and 22 days. Our data showed cells changed into bipolar or multipolar shape forming the cell clusters like scattered rosettes. Nestin expression decreased significantly at 22 days. Compared to 6 days, percentage of Olig2 + MNPs was higher, (88.53 ± 6.67)%, and Olig2 expression at protein and gene level was lower at 22 days. Percentage of MAP2 positive cells increased to (90.62 ± 2.31) % and ChAT positive cells increased to (83.29 ± 6.62) % at 22 days. But no expression of ChAT was examined by western blot and real time PCR. It indicates that these 5 molecular can differentiate hfNCSCs into Olig2 positive cells with a unipotent differentiation toward motor neurons.