Abstract:BackgroundDeep understanding the differentiation process of human embryonic stem cells (hESCs) is essential for developing cell-based therapeutic strategy. Substantial efforts have been made to investigate protein-coding genes, yet it remains lacking comprehensive characterization of long non-coding RNAs (lncRNAs) during this process.MethodshESCs were passaged every 5-6 days and had maintained stable karyotype even until the 50th generation. Pancreatic progenitor specification of in vitro differentiation from hESCs was performed and modified. The nuclei were stained with 4,6-Diamidino-2-phenylindole (DAPI). Droplet-based platform (10X Genomics) was applied to generate the single-cell RNA sequencing (scRNA-seq) data. The quality of the filtered read pairs was evaluated by using FastQC. Batch effects were removed using the size factor method. Dimension reduction and unsupervised clustering analyses were performed using Seurat R package. The Monocle 2 and MetaCell algorithms were used to order single cells on a pseudotime course and partition the scRNA-seq data into metacells, respectively. Co-expression network was constructed using WGCNA. Module- and hub-based methods were adopted to predict the functions of lncRNAs.ResultsA total of 77,382 cells during the differentiation process of hESCs toward pancreatic progenitors were sequenced. According to the single-cell map, the cells from different time points were authenticated to constitute a relatively homogeneous population, in which a total of 7382 lncRNAs could be detected. Through further analyzing the time course data, conserved and specific expression features of lncRNAs during hESC differentiation were revealed. Based upon pseudotime analysis, 52 pseudotime-associated lncRNAs that grouped into three distinct expression patterns were identified. We also implemented MetaCell algorithm and network-based methods to explore the functional mechanisms of these lncRNAs. Totally, 464 lncRNAs, including 49 pseudotime-associated lncRNAs were functionally annotated by either module-based or hub-based methods. Most importantly, we demonstrated that the lncRNA HOTAIRM1, which co-localized and co-expressed with several HOX genes, may play crucial role in the generation of pancreatic progenitors through regulation of exocytosis and retinoic acid receptor signaling pathway.ConclusionsOur single-cell analyses provide valuable data resources for biological researchers and novel insights into hESC differentiation processes, which will guide future endeavors to further elucidate the roles of lncRNAs.

LncRNA SOX1‐OT V1 acts as a decoy of HDAC10 to promote SOX1‐dependent hESC neuronal differentiation

Induction of Oligodendrocyte Differentiation by Olig2 and Sox10: Evidence for Reciprocal Interactions and Dosage-Dependent Mechanisms.

SOX1 Is Required for the Specification of Rostral Hindbrain Neural Progenitor Cells from Human Embryonic Stem Cells

Long Noncoding RNA-1604 Orchestrates Neural Differentiation Through the Mir-200C/zeb Axis

Comprehensive Profiling Reveals Mechanisms of SOX2-mediated Cell Fate Specification in Human ESCs and NPCs

Nucleosome Eviction along with H3K9ac Deposition Enhances Sox2 Binding During Human Neuroectodermal Commitment

Linc1548 Promotes the Transition of Epiblast Stem Cells into Neural Progenitors by Engaging OCT6 and SOX2

LincRNA1230 Inhibits the Differentiation of Mouse ES Cells Towards Neural Progenitors

Long Noncoding RNA Q Associates with Sox2 and is Involved in the Maintenance of Pluripotency in Mouse Embryonic Stem Cells

Opposing Roles for the Lncrna Haunt and Its Genomic Locus in Regulating HOXA Gene Activation During Embryonic Stem Cell Differentiation.

A Sox2 Enhancer Cluster Regulates Region-Specific Neural Fates from Mouse Embryonic Stem Cells

OCT4 and PAX6 Determine the Dual Function of SOX2 in Human ESCs As a Key Pluripotent or Neural Factor.

SOX21 Ensures Rostral Forebrain Identity by Suppression of WNT8B during Neural Regionalization of Human Embryonic Stem Cells.

<i>PAUPAR</i> and PAX6 sequentially regulate human embryonic stem cell cortical differentiation

Single-cell Analyses Reveal Distinct Expression Patterns and Roles of Long Non-Coding RNAs During Hesc Differentiation into Pancreatic Progenitors.

lncRNA Functional Networks in Oligodendrocytes Reveal Stage-Specific Myelination Control by an lncOL1/Suz12 Complex in the CNS.

SOXC Enhances NGN2‐Mediated Reprogramming of Glioblastoma Cells Into Neuron‐Like Cells by Modulating RhoA and RAC1/CDC42 Pathway Activity

LncRNA Riken-201 and Riken-203 modulates neural development by regulating the Sox6 through sequestering miRNAs.

Long non-coding RNA LHX1-DT regulates cardiomyocyte differentiation through H2A.Z-mediated LHX1 transcriptional activation

LincRNA-1614 Coordinates Sox2/PRC2-mediated Repression of Developmental Genes in Pluripotency Maintenance

Sox14 is essential for initiation of neuronal differentiation in the chick spinal cord