ZSCAN4-binding motif-TGCACAC is conserved and enriched in CA/TG microsatellites in both mouse and human genomes.

The Zinc finger and SCAN domain containing 4 (ZSCAN4) protein, expressed transiently in pluripotent stem cells, gametes, and early embryos, extends telomeres, enhances genome stability, and improves karyotypes in mouse embryonic stem (mES) cells. To gain insights into the mechanism of ZSCAN4 function, we identified genome-wide binding sites of endogenous ZSCAN4 protein using ChIP-seq technology in mouse and human ES cells, where the expression of endogenous ZSCAN4 was induced by treating cells with retinoic acids or by overexpressing DUX4. We revealed that both mouse and human ZSCAN4 bind to the TGCACAC motif located in CA/TG microsatellite repeats, which are known to form unstable left-handed duplexes called Z-DNA that can induce double-strand DNA breaks and mutations. These ZSCAN4 binding sites are mostly located in intergenic and intronic regions of the genomes. By generating ZSCAN4 knockout in human ES cells, we showed that ZSCAN4 does not seem to be involved in transcriptional regulation. We also found that ectopic expression of mouse ZSCAN4 enhances the suppression of chromatin at ZSCAN4-binding sites. These results together suggest that some of the ZSCAN4 functions are mediated by binding to the error-prone regions in mouse and human genomes.

A Phase I/II Clinical Trial of Intradermal, Controllable Self-Replicating Ribonucleic Acid Vaccine EXG-5003 against SARS-CoV-2.

mRNA vaccines against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have played a key role in reducing morbidity and mortality from coronavirus disease 2019 (COVID-19). We conducted a double-blind, placebo-controlled phase I/II trial to evaluate the safety, tolerability, and immunogenicity of EXG-5003, a two-dose, controllable self-replicating RNA vaccine against SARS-CoV-2. EXG-5003 encodes the receptor binding domain (RBD) of SARS-CoV-2 and was administered intradermally without lipid nanoparticles (LNPs). The participants were followed for 12 months. Forty healthy participants were enrolled in Cohort 1 (5 µg per dose, n = 16; placebo, n = 4) and Cohort 2 (25 µg per dose, n = 16; placebo, n = 4). No safety concerns were observed with EXG-5003 administration. SARS-CoV-2 RBD antibody titers and neutralizing antibody titers were not elevated in either cohort. Elicitation of antigen-specific cellular immunity was observed in the EXG-5003 recipients in Cohort 2. At the 12-month follow-up, participants who had received an approved mRNA vaccine (BNT162b2 or mRNA-1273) >1 month after receiving the second dose of EXG-5003 showed higher cellular responses compared with equivalently vaccinated participants in the placebo group. The findings suggest a priming effect of EXG-5003 on the long-term cellular immunity of approved SARS-CoV-2 mRNA vaccines.

SARS-CoV-2 infection establishes a stable and age-independent CD8+ T cell response against a dominant nucleocapsid epitope using restricted T cell receptors.

The resolution of SARS-CoV-2 replication hinges on cell-mediated immunity, wherein CD8+ T cells play a vital role. Nonetheless, the characterization of the specificity and TCR composition of CD8+ T cells targeting non-spike protein of SARS-CoV-2 before and after infection remains incomplete. Here, we analyzed CD8+ T cells recognizing six epitopes from the SARS-CoV-2 nucleocapsid (N) protein and found that SARS-CoV-2 infection slightly increased the frequencies of N-recognizing CD8+ T cells but significantly enhanced activation-induced proliferation compared to that of the uninfected donors. The frequencies of N-specific CD8+ T cells and their proliferative response to stimulation did not decrease over one year. We identified the N222-230 peptide (LLLDRLNQL, referred to as LLL thereafter) as a dominant epitope that elicited the greatest proliferative response from both convalescent and uninfected donors. Single-cell sequencing of T cell receptors (TCR) from LLL-specific CD8+ T cells revealed highly restricted Vα gene usage (TRAV12-2) with limited CDR3α motifs, supported by structural characterization of the TCR-LLL-HLA-A2 complex. Lastly, transcriptome analysis of LLL-specific CD8+ T cells from donors who had expansion (expanders) or no expansion (non-expanders) after in vitro stimulation identified increased chromatin modification and innate immune functions of CD8+ T cells in non-expanders. These results suggests that SARS-CoV-2 infection induces LLL-specific CD8+ T cell responses with a restricted TCR repertoire.

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.

Functional and long-lived melanocytes from human pluripotent stem cells with transient ectopic expression of JMJD3.

BACKGROUND: Melanocytes are an essential part of the epidermis, and their regeneration has received much attention because propagation of human adult melanocytes in vitro is too slow for clinical use. Differentiation from human pluripotent stem cells to melanocytes has been reported, but the protocols to produce them require multiple and complex differentiation steps. METHOD: We differentiated human embryonic stem cells (hESCs) that transiently express JMJD3 to pigmented cells. We investigated whether the pigmented cells have melanocytic characteristics and functions by qRT-PCR, immunocytochemical analysis and flow cytometry. We also investigated their biocompatibility by injecting the cells into immunodeficient mice for clinical use. RESULT: We successfully differentiated and established a pure culture of melanocytes. The melanocytes maintained their growth rate for a long time, approximately 200 days, and were functional. They exhibited melanogenesis and transfer of melanin to peripheral keratinocytes. Moreover, melanocytes simulated the developmental processes from melanoblasts to melanocytes. The melanocytes had high engraftability and biocompatibility in the immunodeficient mice. CONCLUSION: The robust generation of functional and long-lived melanocytes are key to developing clinical applications for the treatment of pigmentary skin disorders.

Controllable self-replicating RNA vaccine delivered intradermally elicits predominantly cellular immunity.

Intradermal delivery of self-replicating RNA (srRNA) is a promising vaccine platform. We have developed an srRNA that functions optimally at around 33°C (skin temperature) and is inactivated at or above 37°C (core body temperature) as a safety switch. This temperature-controllable srRNA (c-srRNA), when tested as an intradermal vaccine against SARS-CoV-2, functions when injected naked without lipid nanoparticles. Unlike most currently available vaccines, c-srRNA vaccines predominantly elicit cellular immunity with little or no antibody production. Interestingly, c-srRNA-vaccinated mice produced antigen-specific antibodies upon subsequent stimulation with antigen protein. Antigen-specific antibodies were also produced when B cell stimulation using antigen protein was followed by c-srRNA booster vaccination. We have thus designed a pan-coronavirus booster vaccine that incorporates both spike-receptor-binding domains as viral surface proteins and evolutionarily conserved nucleoproteins as viral internal proteins, from both severe acute respiratory syndrome coronavirus 2 and Middle East respiratory syndrome coronavirus. c-srRNA may provide a route to activate cellular immunity against a wide variety of pathogens.

Controllable self-replicating RNA vaccine delivered intradermally elicits predominantly cellular immunity.

Intradermal delivery of self-replicating RNA (srRNA) is a promising vaccine platform. Considering that human skin temperature is around 33°C, lower than core body temperature of 37°C, we have developed an srRNA that functions optimally at skin temperature and is inactivated at or above 37°C as a safety switch. This temperature- c ontrollable srRNA (c-srRNA), when tested as an intradermal vaccine against SARS-CoV-2, functions when injected naked without lipid nanoparticles. Unlike most currently available vaccines, c-srRNA vaccines predominantly elicit cellular immunity with little or no antibody production. Interestingly, c-srRNA-vaccinated mice produced antigen-specific antibodies upon subsequent stimulation with antigen protein. Antigen-specific antibodies were also produced when B-cell stimulation using antigen protein was followed by c-srRNA booster vaccination. Using c-srRNA, we have designed a pan-coronavirus booster vaccine that incorporates both spike receptor binding domains as viral surface proteins and evolutionarily conserved nucleoproteins as viral non-surface proteins, from both SARS-CoV-2 and MERS-CoV. It can thereby potentially immunize against SARS-CoV-2, SARS-CoV, MERS-CoV, and their variants. c-srRNA may provide a route to activate cellular immunity against a wide variety of pathogens.

Purification of cardiomyocytes and neurons derived from human pluripotent stem cells by inhibition of de novo fatty acid synthesis.

Here we describe a protocol to obtain highly pure cardiomyocytes and neurons from human induced pluripotent stem cells (hiPSCs) via metabolic selection processes. Compared to conventional purification protocols, this approach is easier to perform and scale up and more cost-efficient. The protocol can be applied to hiPSCs and human embryonic stem cells. For complete details on the use and execution of this protocol, please refer to Tohyama et al. (2016) and Tanosaki et al. (2020).

Synthetic mRNA-based differentiation method enables early detection of Parkinson's phenotypes in neurons derived from Gaucher disease-induced pluripotent stem cells.

Gaucher disease, the most prevalent metabolic storage disorder, is caused by mutations in the glucocerebrosidase gene GBA1, which lead to the accumulation of glucosylceramide (GlcCer) in affected cells. Gaucher disease type 1 (GD1), although defined as a nonneuronopathic subtype, is accompanied by an increased risk of Parkinson's disease. To gain insights into the association of progressive accumulation of GlcCer and the Parkinson's disease phenotypes, we generated dopaminergic (DA) neurons from induced pluripotent stem cells (iPSCs) derived from a GD1 patient and a healthy donor control, and measured GlcCer accumulation by liquid chromatography-mass spectrometry. We tested two DA neuron differentiation methods: a well-established method that mimics a step-wise developmental process from iPSCs to neural progenitor cells, and to DA neurons; and a synthetic mRNA-based method that overexpresses a transcription factor in iPSCs. GD1-specific accumulation of GlcCer was detected after 60 days of differentiation by the former method, whereas it was detected after only 10 days by the latter method. With this synthetic mRNA-based rapid differentiation method, we found that the metabolic defect in GD1 patient cells can be rescued by the overexpression of wild-type GBA1 or the treatment with an inhibitor for GlcCer synthesis. Furthermore, we detected the increased phosphorylation of α-synuclein, a biomarker for Parkinson's disease, in DA neurons derived from a GD1 patient, which was significantly decreased by the overexpression of wild-type GBA1. These results suggest that synthetic mRNA-based method accelerates the analyses of the pathological mechanisms of Parkinson's disease in GD1 patients and possibly facilitates drug discovery processes.

Fatty Acid Synthesis Is Indispensable for Survival of Human Pluripotent Stem Cells.

The role of lipid metabolism in human pluripotent stem cells (hPSCs) is poorly understood. We have used large-scale targeted proteomics to demonstrate that undifferentiated hPSCs express different fatty acid (FA) biosynthesis-related enzymes, including ATP citrate lyase and FA synthase (FASN), than those expressed in hPSC-derived cardiomyocytes (hPSC-CMs). Detailed lipid profiling revealed that inhibition of FASN resulted in significant reduction of sphingolipids and phosphatidylcholine (PC); moreover, we found that PC was the key metabolite for cell survival in hPSCs. Inhibition of FASN induced cell death in undifferentiated hPSCs via mitochondria-mediated apoptosis; however, it did not affect cell survival in hPSC-CMs, neurons, or hepatocytes as there was no significant reduction of PC. Furthermore, we did not observe tumor formation following transplantation of FASN inhibitor-treated cells. Our findings demonstrate the importance of de novo FA synthesis in the survival of undifferentiated hPSCs and suggest applications for FASN inhibition in regenerative medicine.

Generation and Profiling of 2,135 Human ESC Lines for the Systematic Analyses of Cell States Perturbed by Inducing Single Transcription Factors.

Transcription factors (TFs) play a pivotal role in determining cell states, yet our understanding of the causative relationship between TFs and cell states is limited. Here, we systematically examine the state changes of human pluripotent embryonic stem cells (hESCs) by the large-scale manipulation of single TFs. We establish 2,135 hESC lines, representing three clones each of 714 doxycycline (Dox)-inducible genes including 481 TFs, and obtain 26,998 microscopic cell images and 2,174 transcriptome datasets-RNA sequencing (RNA-seq) or microarrays-48 h after the presence or absence of Dox. Interestingly, the expression of essentially all the genes, including genes located in heterochromatin regions, are perturbed by these TFs. TFs are also characterized by their ability to induce differentiation of hESCs into specific cell lineages. These analyses help to provide a way of classifying TFs and identifying specific sets of TFs for directing hESC differentiation into desired cell types.

MEIOSIN Directs the Switch from Mitosis to Meiosis in Mammalian Germ Cells.

The mechanisms regulating meiotic initiation in mammals are enigmatic. It is known that retinoic acid (RA) signaling plays a pivotal role during meiotic initiation. STRA8, which is expressed in response to RA, is thought to be a key factor promoting meiotic initiation. However, the specific role of STRA8 in meiotic initiation has remained elusive. Here, we identified MEIOSIN as a germ-cell-specific factor that associates with STRA8. MEIOSIN, like STRA8, is expressed in response to RA and plays an essential role in meiotic initiation in both males and females. Functional analyses revealed that MEIOSIN acts as a transcription factor together with STRA8, and that both factors are critical for driving meiotic gene activation. Furthermore, temporally restricted expression of MEIOSIN leads to meiotic entry decision during spermatogenesis. The present study demonstrates that MEIOSIN, in collaboration with STRA8, plays a central role in regulating the mitosis to meiosis germ cell fate decision in mammals.

Induced Pluripotent Stem Cells Reprogrammed with Three Inhibitors Show Accelerated Differentiation Potentials with High Levels of 2-Cell Stage Marker Expression.

Although pluripotent stem cells can generate various types of differentiated cells, it is unclear why lineage-committed stem/progenitor cells derived from pluripotent stem cells are decelerated and why the differentiation-resistant propensity of embryonic stem cell (ESC)/induced pluripotent stem cell (iPSC)-derived cells is predominant compared with the in vivo equivalents derived from embryonic/adult tissues. In this study, we demonstrated that iPSCs reprogrammed and maintained with three chemical inhibitors of the fibroblast growth factor 4-mitogen-activated protein kinase cascade and GSK3ß (3i) could be differentiated into all three germ layers more efficiently than the iPSCs reprogrammed without the 3i chemicals, even though they were maintained with 3i chemicals once they were reprogrammed. Although the iPSCs reprogrammed with 3i had increased numbers of Zscan4-positive cells, the Zscan4-positive cells among iPSCs that were reprogrammed without 3i did not have an accelerated differentiation ability. These observations suggest that 3i exposure during the reprogramming period determines the accelerated differentiation/maturation potentials of iPSCs that are stably maintained at the distinct state.

Induction of human pluripotent stem cells into kidney tissues by synthetic mRNAs encoding transcription factors.

The derivation of kidney tissues from human pluripotent stem cells (hPSCs) and its application for replacement therapy in end-stage renal disease have been widely discussed. Here we report that consecutive transfections of two sets of synthetic mRNAs encoding transcription factors can induce rapid and efficient differentiation of hPSCs into kidney tissues, termed induced nephron-like organoids (iNephLOs). The first set - FIGLA, PITX2, ASCL1 and TFAP2C, differentiated hPSCs into SIX2+SALL1+ nephron progenitor cells with 92% efficiency within 2 days. Subsequently, the second set - HNF1A, GATA3, GATA1 and EMX2, differentiated these cells into PAX8+LHX1+ pretubular aggregates in another 2 days. Further culture in both 2-dimensional and 3-dimensional conditions produced iNephLOs containing cells characterized as podocytes, proximal tubules, and distal tubules in an additional 10 days. Global gene expression profiles showed similarities between iNephLOs and the human adult kidney, suggesting possible uses of iNephLOs as in vitro models for kidneys.

Establishment of a rapid and footprint-free protocol for differentiation of human embryonic stem cells into pancreatic endocrine cells with synthetic mRNAs encoding transcription factors.

BACKGROUND: Transplantation of pancreatic ß cells generated in vitro from pluripotent stem cells (hPSCs) such as embryonic stem cells (ESCs) or induced pluripotent stem cells (iPSCs) has been proposed as an alternative therapy for diabetes. Though many differentiation protocols have been developed for this purpose, lentivirus-mediated forced expression of transcription factors (TF)-PDX1 and NKX6.1-has been at the forefront for its relatively fast and straightforward approach. However, considering that such cells will be used for therapeutic purposes in the future, it is desirable to develop a procedure that does not leave any footprint on the genome, as any changes of DNAs could potentially be a source of unintended, concerning effects such as tumorigenicity. In this study, we attempted to establish a novel protocol for rapid and footprint-free hESC differentiation into a pancreatic endocrine lineage by using synthetic mRNAs (synRNAs) encoding PDX1 and NKX6.1. We also tested whether siPOU5F1, which reduces the expression of pluripotency gene POU5F1 (also known as OCT4), can enhance differentiation as reported previously for mesoderm and endoderm lineages. METHODS: synRNA-PDX1 and synRNA-NKX6.1 were synthesized in vitro and were transfected five times to hESCs with a lipofection reagent in a modified differentiation culture condition. siPOU5F1 was included only in the first transfection. Subsequently, cells were seeded onto a low attachment plate and aggregated by an orbital shaker. At day 13, the degree of differentiation was assessed by quantitative RT-PCR (qRT-PCR) and immunohistochemistry for endocrine hormones such as insulin, glucagon, and somatostatin. RESULTS: Both PDX1 and NKX6.1 expression were detected in cells co-transfected with synRNA-PDX1 and synRNA-NKX6.1 at day 3. Expression levels of insulin in the transfected cells at day 13 were 450 times and 14 times higher by qRT-PCR compared to the levels at day 0 and in cells cultured without synRNA transfection, respectively. Immunohistochemically, pancreatic endocrine hormones were not detected in cells cultured without synRNA transfection but were highly expressed in cells transfected with synRNA-PDX1, synRNA-NKX6.1, and siPOU5F1 at as early as day 13. CONCLUSIONS: In this study, we report a novel protocol for rapid and footprint-free differentiation of hESCs to endocrine cells.

Efficient differentiation of human pluripotent stem cells into skeletal muscle cells by combining RNA-based MYOD1-expression and POU5F1-silencing.

Direct generation of skeletal muscle cells from human pluripotent stem cells (hPSCs) would be beneficial for drug testing, drug discovery, and disease modelling in vitro. Here we show a rapid and robust method to induce myogenic differentiation of hPSCs by introducing mRNA encoding MYOD1 together with siRNA-mediated knockdown of POU5F1 (also known as OCT4 or OCT3/4). This integration-free approach generates functional skeletal myotubes with sarcomere-like structure and a fusion capacity in several days. The POU5F1 silencing facilitates MYOD1 recruitment to the target promoters, which results in the significant activation of myogenic genes in hPSCs. Furthermore, deep sequencing transcriptome analyses demonstrated that POU5F1-knockdown upregulates the genes associated with IGF- and FGF-signaling and extracellular matrix that may also support myogenic differentiation. This rapid and direct differentiation method may have potential applications in regenerative medicine and disease therapeutics for muscle disorders such as muscular dystrophy.

Identification of transcription factors that promote the differentiation of human pluripotent stem cells into lacrimal gland epithelium-like cells.

Dry eye disease is the most prevalent pathological condition in aging eyes. One potential therapeutic strategy is the transplantation of lacrimal glands, generated in vitro from pluripotent stem cells such as human embryonic stem cells, into patients. One of the preceding requirements is a method to differentiate human embryonic stem cells into lacrimal gland epithelium cells. As the first step for this approach, this study aims to identify a set of transcription factors whose overexpression can promote the differentiation of human embryonic stem cells into lacrimal gland epithelium-like cells. We performed microarray analyses of lacrimal glands and lacrimal glands-related organs obtained from mouse embryos and adults, and identified transcription factors enriched in lacrimal gland epithelium cells. We then transfected synthetic messenger RNAs encoding human orthologues of these transcription factors into human embryonic stem cells and examined whether the human embryonic stem cells differentiate into lacrimal gland epithelium-like cells by assessing cell morphology and marker gene expression. The microarray analysis of lacrimal glands tissues identified 16 transcription factors that were enriched in lacrimal gland epithelium cells. We focused on three of the transcription factors, because they are expressed in other glands such as salivary glands and are also known to be involved in the development of lacrimal glands. We tested the overexpression of various combinations of the three transcription factors and PAX6, which is an indispensable gene for lacrimal glands development, in human embryonic stem cells. Combining PAX6, SIX1, and FOXC1 caused significant changes in morphology, i.e., elongated cell shape and increased expression (both RNAs and proteins) of epithelial markers such as cytokeratin15, branching morphogenesis markers such as BARX2, and lacrimal glands markers such as aquaporin5 and lactoferrin. We identified a set of transcription factors enriched in lacrimal gland epithelium cells and demonstrated that the simultaneous overexpression of these transcription factors can differentiate human embryonic stem cells into lacrimal gland epithelium-like cells. This study suggests the possibility of lacrimal glands regeneration from human pluripotent stem cells.

Neural differentiation of human embryonic stem cells induced by the transgene-mediated overexpression of single transcription factors.

Pluripotent human embryonic stem cells (hESCs) can differentiate into multiple cell lineages, thus, providing one of the best platforms to study molecular mechanisms during cell differentiation. Recently, we have reported rapid and efficient differentiation of hESCs into functional neurons by introducing a cocktail of synthetic mRNAs encoding five transcription factors (TFs): NEUROG1, NEUROG2, NEUROG3, NEUROD1, and NEUROD2. Here we further tested a possibility that even single transcription factors, when expressed ectopically, can differentiate hESCs into neurons. To this end, we established hESC lines in which each of these TFs can be overexpressed by the doxycycline-inducible piggyBac vector. The overexpression of any of these five TFs indeed caused a rapid and rather uniform differentiation of hESCs, which were identified as neurons based on their morphologies, qRT-PCR, and immunohistochemistry. Furthermore, calcium-imaging analyses and patch clamp recordings demonstrated that these differentiated cells are electrophysiologically functional. Interestingly, neural differentiations occurred despite the cell culture conditions that rather promote the maintenance of the undifferentiated state. These results indicate that over-expression of each of these five TFs can override the pluripotency-specific gene network and force hESCs to differentiate into neurons.

Epigenetic Manipulation Facilitates the Generation of Skeletal Muscle Cells from Pluripotent Stem Cells.

Human pluripotent stem cells (hPSCs) have the capacity to differentiate into essentially all cell types in the body. Such differentiation can be directed to specific cell types by appropriate cell culture conditions or overexpressing lineage-defining transcription factors (TFs). Especially, for the activation of myogenic program, early studies have shown the effectiveness of enforced expression of TFs associated with myogenic differentiation, such as PAX7 and MYOD1. However, the efficiency of direct differentiation was rather low, most likely due to chromatin features unique to hPSCs, which hinder the access of TFs to genes involved in muscle differentiation. Indeed, recent studies have demonstrated that ectopic expression of epigenetic-modifying factors such as a histone demethylase and an ATP-dependent remodeling factor significantly enhances myogenic differentiation from hPSCs. In this article, we review the recent progress for in vitro generation of skeletal muscles from hPSCs through forced epigenetic and transcriptional manipulation.

SCODE: an efficient regulatory network inference algorithm from single-cell RNA-Seq during differentiation.

MOTIVATION: The analysis of RNA-Seq data from individual differentiating cells enables us to reconstruct the differentiation process and the degree of differentiation (in pseudo-time) of each cell. Such analyses can reveal detailed expression dynamics and functional relationships for differentiation. To further elucidate differentiation processes, more insight into gene regulatory networks is required. The pseudo-time can be regarded as time information and, therefore, single-cell RNA-Seq data are time-course data with high time resolution. Although time-course data are useful for inferring networks, conventional inference algorithms for such data suffer from high time complexity when the number of samples and genes is large. Therefore, a novel algorithm is necessary to infer networks from single-cell RNA-Seq during differentiation. RESULTS: In this study, we developed the novel and efficient algorithm SCODE to infer regulatory networks, based on ordinary differential equations. We applied SCODE to three single-cell RNA-Seq datasets and confirmed that SCODE can reconstruct observed expression dynamics. We evaluated SCODE by comparing its inferred networks with use of a DNaseI-footprint based network. The performance of SCODE was best for two of the datasets and nearly best for the remaining dataset. We also compared the runtimes and showed that the runtimes for SCODE are significantly shorter than for alternatives. Thus, our algorithm provides a promising approach for further single-cell differentiation analyses. AVAILABILITY AND IMPLEMENTATION: The R source code of SCODE is available at https://github.com/hmatsu1226/SCODE. CONTACT: hirotaka.matsumoto@riken.jp. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.