Background Human embryonic stem cells (hESC) have been an invaluable research tool to study motor neuron development and disorders. However, transcriptional regulation of multiple temporal stages from ESCs to spinal motor neurons (MNs) has not yet been fully elucidated. Thus, the goals of this study were to profile the time-course expression patterns of lncRNAs during MN differentiation of ESCs and to clarify the potential mechanisms of the lncRNAs that are related to MN differentiation. Methods We utilized our previous protocol which can harvest motor neuron in more than 90% purity from hESCs. Then, differentially expressed lncRNAs (DElncRNAs) and mRNAs (DEmRNAs) during MN differentiation were identified through RNA sequencing. Bioinformatic analyses were performed to assess potential biological functions of genes. We also performed qRT-PCR to validate the DElncRNAs and DEmRNAs. Results A total of 441 lncRNAs and 1,068 mRNAs at day 6, 443 and 1,175 at day 12, and 338 lncRNAs and 68 mRNAs at day 18 were differentially expressed compared with day 0. Bioinformatic analyses identified that several key regulatory genes including POU5F1, TDGF1, SOX17, LEFTY2 and ZSCAN10, which involved in the regulation of embryonic development. We also predicted 283 target genes of DElncRNAs, in which 6 mRNAs were differentially expressed. Significant fold changes in lncRNAs (NCAM1-AS) and mRNAs (HOXA3) were confirmed by qRT-PCR. Then, through predicted overlapped miRNA verification, we constructed a lncRNA NCAM1-AS-miRNA-HOXA3 network.
